View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001772 | OpenFOAM | Bug | public | 2015-07-02 08:21 | 2015-07-02 10:25 |
| Reporter | sharonyue | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | N/A |
| Status | closed | Resolution | no change required | ||
| Summary | 0001772: d rho / d t in buoyantPimpleFoam and twoPhaseEulerFoam.dev | ||||
| Description | In 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? | ||||
| Tags | No tags attached. | ||||
|
|
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. |
