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: 17	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

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	COMSOL	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: 74 Rep Power: 15	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: 17	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