CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Calculation of DpDt (http://www.cfd-online.com/Forums/openfoam-solving/59487-calculation-dpdt.html)

stephan August 21, 2007 08:57

hi, after months of simply
 
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 August 22, 2007 13:49

hi, any help would be reall
 
hi,

any help would be really appreciated!
regards
stephan

hjasak August 22, 2007 14:03

Well: DpDt = ddt(p) + U . g
 
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 http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

Hrv

stephan August 22, 2007 15:06

hi, at first - thanks a lot
 
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


All times are GMT -4. The time now is 12:41.