CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Adding additional complex source term to BuoyantBoussinesqPisoFoam solver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 6, 2012, 16:50
Default Adding additional complex source term to BuoyantBoussinesqPisoFoam solver
  #1
New Member
 
Paul
Join Date: Mar 2009
Posts: 19
Rep Power: 17
pajofego is on a distinguished road
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 F_i 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 x_1 and x_2

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
pajofego is offline   Reply With Quote

Old   April 9, 2013, 04:07
Default
  #2
Member
 
Neilson Whit
Join Date: Aug 2011
Posts: 74
Rep Power: 14
wolfindark is on a distinguished road
Hello pajofego

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

thanks
wolfindark is offline   Reply With Quote

Old   April 9, 2013, 18:00
Default
  #3
New Member
 
Paul
Join Date: Mar 2009
Posts: 19
Rep Power: 17
pajofego is on a distinguished road
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
pajofego is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 07:14.