CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

origin of fvm::Sp(fvc::div(phi_), epsilon_) in kepsilon Eqn?

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

Like Tree2Likes
  • 2 Post By hjasak

Reply
 
LinkBack Thread Tools Display Modes
Old   May 15, 2009, 17:55
Default origin of fvm::Sp(fvc::div(phi_), epsilon_) in kepsilon Eqn?
  #1
Bob
New Member
 
Bob De Clercq
Join Date: Apr 2009
Location: Belgium
Posts: 17
Rep Power: 9
Bob is on a distinguished road
Dear all,

I am trying to figure out what some terms are in the kepsilon turbulence equations, i.e.

- fvm::Sp(fvc::div(phi_), epsilon_)

Does someone know and can tell me where it originates from?
A similar term can be seen in the equation for kinetic energy.

Further on, I was wondering how the continuity equation is accounted for with respect to incompressible flows. I cannot find it in e.g. icoFoam.C. Is it somewhere integrated?

Many thanks in advance!

Regards,
Bob
Bob is offline   Reply With Quote

Old   May 16, 2009, 02:11
Default
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,781
Rep Power: 22
hjasak will become famous soon enough
The term you refer to is added to account for boundedness of k when the model is used with simpleFoam (and alike) and partial convergence of the pressure equation. In other cases, div(flux) will be zero, so the term disappears.

Regarding icoFoam, continuity equation is the pressure equation. Have a look at my Thesis or CFD lecture notes - the full explanation is done in about 2 pages.

Enjoy,

Hrv
DheepanPillai and MaLa like this.
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 18, 2009, 15:16
Default
  #3
Bob
New Member
 
Bob De Clercq
Join Date: Apr 2009
Location: Belgium
Posts: 17
Rep Power: 9
Bob is on a distinguished road
Dear Hrv,

Many thanks for the information. I look further into it!

Cheers,
Bob
Bob is offline   Reply With Quote

Old   May 25, 2009, 15:49
Default
  #4
Bob
New Member
 
Bob De Clercq
Join Date: Apr 2009
Location: Belgium
Posts: 17
Rep Power: 9
Bob is on a distinguished road
Dear Hrvoje,

I am still trying to figure out the term:

- fvm::Sp(fvc::div(phi_), epsilon_)

You mentioned that it is included in the equations to account for boundedness of k, epsilon, or the density (in case of compressible fluids). Do you perhaps have a reference where this is further explained so I can figure out the origin of the term?

Many thanks.

Regards,
Bob
Bob is offline   Reply With Quote

Old   May 26, 2009, 17:35
Default
  #5
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,781
Rep Power: 22
hjasak will become famous soon enough
So:

div(phi, epsilon) - in translation: div(U epsilon)

will be bounded only if div(phi) - in translation div(U)) is zero, right?

On the other hand,

U & grad epsilon

is bounded irrespective of the velocity distribution. Thus, on incomplete convergence, we prefer to solve a weak form of convection (U & grad epsilon) rather than the strong form.

Final step:

div(U epsilon) = U & grad epsilon - epsilon div(U).

Clear?

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 29, 2009, 10:07
Default
  #6
Bob
New Member
 
Bob De Clercq
Join Date: Apr 2009
Location: Belgium
Posts: 17
Rep Power: 9
Bob is on a distinguished road
Dear Hrvoje,

Thanks for the explanations, but I thought the identity of your final step would be written as something like:

div(U epsilon) = U & grad epsilon + epsilon div(U).

Is there some mistake from my side?

Cheers,
Bob
Bob is offline   Reply With Quote

Old   July 2, 2014, 02:02
Default
  #7
Member
 
Join Date: Aug 2011
Posts: 83
Rep Power: 7
idefix is on a distinguished road
Hello
I know this thread is older, but if there is someone interested:
here is an explanation for the question:
Convectional term -fvm::Sp(fvc::div(phi_), epsilon_) in RNGkEpsilon OF 1.7.1 missing
idefix is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 13 November 4, 2013 15:13
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21
A velocity dependent kepsilon boundary condition tangd OpenFOAM Pre-Processing 2 September 15, 2006 07:04
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 13, 2000 23:03
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 03:43


All times are GMT -4. The time now is 15:48.