View Issue Details

IDProjectCategoryView StatusLast Update
0001772OpenFOAMBugpublic2015-07-02 10:25
Reportersharonyue Assigned Tohenry  
PrioritynormalSeverityminorReproducibilityN/A
Status closedResolutionno change required 
Summary0001772: d rho / d t in buoyantPimpleFoam and twoPhaseEulerFoam.dev
DescriptionIn the code of buoyantPimpleFoam and twoPhaseEulerFoam.dev. We need to calculate:

d rho/d t

This is my procedure:

rho1 = rho0 + psi*(p1 - p0) = rho0 + psi*(p_rgh1 - p_rgh0 + rho1gh + rho2gh)

But from the code, like buoyantPimpleFoam:

fvc::ddt(rho) + psi*correction(fvm::ddt(p_rgh))

Looks rho1gh + rho2gh terms were neglected. This is similar with twoPhaseEulerFoam.dev.

So does this mean that "rho1gh + rho2gh" is neglected or there is something wrong with my equations?
TagsNo tags attached.

Activities

henry

2015-07-02 09:56

manager   ~0005033

The correction(fvm::ddt(p_rgh)) term is a correction to fvc::ddt(rho) due to changes in pressure and the converged result will be the same irrespective of how this correction is formulated, i.e. the correction will converge to 0 and the results will be fvc::ddt(rho).

Currently the rho*gh terms are lagged and only the direct p_rgh term is included in the correction. Feel free to play with other forms of correction(fvm::ddt(p_rgh)) to see if you can find a form which is more convergent.

Issue History

Date Modified Username Field Change
2015-07-02 08:21 sharonyue New Issue
2015-07-02 09:56 henry Note Added: 0005033
2015-07-02 09:56 henry Status new => closed
2015-07-02 09:56 henry Assigned To => henry
2015-07-02 09:56 henry Resolution open => no change required