Adding additional complex source term to BuoyantBoussinesqPisoFoam solver

pajofego · September 6, 2012, 16:50

Hello all,

I'm looking for an starting point how to add a complex source term to the momentum equation of the BuoyantBoussinesqPisoFoam solver:

As described in the openfoamwiki the governing equation for the momentum is:

$\frac{\partial \bar{u}_i}{\partial t} + \frac{\partial}{\partial x_j} \left( \bar{u}_j \bar{u}_i \right) = -\frac{\partial} {\partial{x_i}} \left( \frac{\bar{p}}{\rho_0} \right) + \frac{1}{\rho_0}\frac{\partial}{\partial x_j} \left( \tau_{ij} + \tau_{t_{ij}} \right) + \frac{\bar{\rho}}{\rho_0} g_i + \mathbf{\frac{F_i}{\rho_0}}$

The term

stays for the complex source, which is defined as follows:

$F_1 = \rho_0 2\Omega(\sin(\varphi)u_2 - \cos(\varphi)u_3x_3)$

$\varphi$ depends from the position of the coordinate

and

$F_2 = -\rho_0 2\Omega(\sin(\varphi)u_1$ and last for $F_3 = \rho_0(x_3^2-1)(\cos(\varphi)\frac{u_1}{x_3} + \beta(\bar{T}-T)g)+\rho_0\cos(\varphi)\frac{u_1}{x_3}+x_3(\bar{p}-\rho_0u^{2}_{2})$

Does anyone can give me please a hint how I should start to model such complex source terms in the openfoam library? Maybe there is a good example?

Thanks to all and best regards
pajofego

wolfindark · April 9, 2013, 04:07

Hello pajofego

I am working on similar thing. did you find a way?

thanks

pajofego · April 9, 2013, 18:00

Hi,

not really. The only way that I've found, but not tested, is to use swak4foam. http://www.cfd-online.com/Forums/ope...urce-term.html

For more complex source term handling I suppose that's necessary to get familiar with the OF code.

Best regards
pajofego

September 6, 2012, 16:50	Adding additional complex source term to BuoyantBoussinesqPisoFoam solver	#1
pajofego New Member Paul Join Date: Mar 2009 Posts: 19 Rep Power: 13	Hello all, I'm looking for an starting point how to add a complex source term to the momentum equation of the BuoyantBoussinesqPisoFoam solver: As described in the openfoamwiki the governing equation for the momentum is: $\frac{\partial \bar{u}_i}{\partial t} + \frac{\partial}{\partial x_j} \left( \bar{u}_j \bar{u}_i \right) = -\frac{\partial} {\partial{x_i}} \left( \frac{\bar{p}}{\rho_0} \right) + \frac{1}{\rho_0}\frac{\partial}{\partial x_j} \left( \tau_{ij} + \tau_{t_{ij}} \right) + \frac{\bar{\rho}}{\rho_0} g_i + \mathbf{\frac{F_i}{\rho_0}}$ The term stays for the complex source, which is defined as follows: $F_1 = \rho_0 2\Omega(\sin(\varphi)u_2 - \cos(\varphi)u_3x_3)$ $\varphi$ depends from the position of the coordinate and $F_2 = -\rho_0 2\Omega(\sin(\varphi)u_1$ and last for $F_3 = \rho_0(x_3^2-1)(\cos(\varphi)\frac{u_1}{x_3} + \beta(\bar{T}-T)g)+\rho_0\cos(\varphi)\frac{u_1}{x_3}+x_3(\bar{p}-\rho_0u^{2}_{2})$ Does anyone can give me please a hint how I should start to model such complex source terms in the openfoam library? Maybe there is a good example? Thanks to all and best regards pajofego

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Source Term Question	Blackhawks84	Fluent UDF and Scheme Programming	1	July 31, 2012 16:44
lid-driven cavity in matlab using BiCGStab	Don456	Main CFD Forum	1	January 19, 2012 16:00
Adding a volumetric source term in femlab	BS	Main CFD Forum	1	December 28, 2006 21:26
Adding a momentum source term	segersson	OpenFOAM Running, Solving & CFD	5	March 3, 2006 00:06
Adding a source term to continuity eqn(UDF)	Sebeci	FLUENT	0	April 20, 2003 04:12

April 9, 2013, 04:07		#2
wolfindark Member Neilson Whit Join Date: Aug 2011 Posts: 72 Rep Power: 11	Hello pajofego I am working on similar thing. did you find a way? thanks

April 9, 2013, 18:00		#3
pajofego New Member Paul Join Date: Mar 2009 Posts: 19 Rep Power: 13	Hi, not really. The only way that I've found, but not tested, is to use swak4foam. http://www.cfd-online.com/Forums/ope...urce-term.html For more complex source term handling I suppose that's necessary to get familiar with the OF code. Best regards pajofego