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

Calculation of DpDt

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 21, 2007, 08:57
Default hi, after months of simply
  #1
Senior Member
 
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 8
stephan is on a distinguished road
hi,

after months of simply using compressible flow solvers (rhoTurbFOAM) i was wondering if somebody might be able to explain why DpDt=DDt(phi/rho,p) is written the way it is?
especially the "phi/rho" is pretty unclear for me?
i hope DDt(phi/rho,p) actually is really
phi/rho*DDt(p)+p*DDt(phi/rho) otherwise i am really confused.
i hope somebody might explain what this is all about or at least could give me a hint where to search in literature.
thanx in advance
regards
stephan
stephan is offline   Reply With Quote

Old   August 22, 2007, 13:49
Default hi, any help would be reall
  #2
Senior Member
 
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 8
stephan is on a distinguished road
hi,

any help would be really appreciated!
regards
stephan
stephan is offline   Reply With Quote

Old   August 22, 2007, 14:03
Default Well: DpDt = ddt(p) + U . g
  #3
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Well:

DpDt = ddt(p) + U . grad(p),

right? The way to write it in the conservative FVM is to go to the strong conservative form, because

U . grad(p)

carries a non-conservative discretisation error. Thus:

DpDt = ddt(p) + div(U p) - p div(U)

All clear so far? Now, we need the div terms. In the FVM it is the mass conservative flux that is essential for all calculations. Thus for the div(U p) you want the normal way of doing things, i.e. using the face flux instead of interpolated velocity.

But what about div(U)? For consistency, you should not just interpolate the U: instead, do flux/rho, and take a div of that. This equals div(U) but the discretisation error is consistent. To make things better, you can interpolate rho onto the face with more consistency with the density equation, but that's another story...

See it now? As an exercise, stir in the moving mesh terms now

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

Old   August 22, 2007, 15:06
Default hi, at first - thanks a lot
  #4
Senior Member
 
Stephan Gerber
Join Date: Mar 2009
Location: Germany
Posts: 118
Rep Power: 8
stephan is on a distinguished road
hi,

at first - thanks a lot...
so i guess my misunderstanding was that DDt(a,b) is something like a*DDt(b)+b*DDt(a) (this idea came from the behaivior of ddt(a,b)) instead of ddt(a)+b.grad(a).

maybe you could please say some things about "another story..." because i guess this will hit some other problems of my solvers.


regards
stephan
stephan 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
calculation of y+ value oe Main CFD Forum 1 September 23, 2008 14:41
NOx calculation Giu CFX 1 November 27, 2006 15:38
CD calculation Alireza CFX 1 May 9, 2006 08:52
udf for Cd calculation Micaela Demichela FLUENT 2 November 16, 2001 00:40
BC for fan calculation Mark Render Main CFD Forum 0 May 2, 2001 09:48


All times are GMT -4. The time now is 09:53.