|
[Sponsors] |
Conservative formulation in OpenFOAm multiphase solvers |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 11, 2022, 12:38 |
Conservative formulation in OpenFOAm multiphase solvers
|
#1 |
New Member
Wojciech Sadowski
Join Date: Feb 2021
Posts: 1
Rep Power: 0 |
Hello everyone,
So I can't find an explanation for the way that momentum equation is coded in OF multiphase solvers. For example DPMFoam: Code:
fvVectorMatrix UcEqn ( fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc) - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) + continuousPhaseTurbulence->divDevRhoReff(Uc) == (1.0/rhoc)*cloudSU ); Code:
U1Eqn = ( fvm::ddt(alpha1, rho1, U1) + fvm::div(alphaRhoPhi1, U1) - fvm::Sp(contErr1, U1) + MRF.DDt(alpha1*rho1 + Vm, U1) + phase1.turbulence().divDevRhoReff(U1) == ... ); Here is the issue, I do recognise this as what appears from the decomposition of the conservative form and it would be logical in my opinion if this was added. My reasoning looks as follows: we have an equation in conservative form, we change it to the non-conservative, so that we solve for instead of : Instead, it is subtracted, and I don't get why. Additionally this is supposed to make the equations conservative (look at this commit by Henry, it describes the change of the UEqns.H in twoPhaseEulerFoam in OpenFOAM 2.3x). So the question is, can someone explain or provide sources that explain what this source term actually does and how does it make the equations conservative? |
|
Tags |
conservative form, momentum equation, multiphase flow |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Finding OpenFOAM solvers on two-phase cavitating flows (in cryogenic conditions) | JD_PM | OpenFOAM Running, Solving & CFD | 7 | December 21, 2022 06:54 |
OpenFOAM Training Jan-Apr 2017, Virtual, London, Houston, Berlin | cfd.direct | OpenFOAM Announcements from Other Sources | 0 | September 21, 2016 11:50 |
OpenFOAM Training, London, Chicago, Munich, Houston 2016-2017 | cfd.direct | OpenFOAM Announcements from Other Sources | 0 | September 14, 2016 03:19 |
Multiphase and FreeSurface Flows at OpenFOAM Workshop Milan 2008 | egp | OpenFOAM | 0 | March 20, 2008 06:34 |
OpenFOAM Debian packaging current status problems and TODOs | oseen | OpenFOAM Installation | 9 | August 26, 2007 13:50 |