multiphaseEulerFoam: Improved consistency between implicit and explicit parts of the virtual mass

This commit is contained in:
Henry
2012-03-14 18:18:48 +00:00
parent dcb26898a5
commit f3fac5bbff

View File

@ -14,13 +14,18 @@ forAllIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
( (
phasei, phasei,
new fvVectorMatrix new fvVectorMatrix
(
(scalar(1) + fluid.Cvm(phase)/phase.rho())*
( (
fvm::ddt(alpha, U) fvm::ddt(alpha, U)
+ fvm::div(phase.phiAlpha(), U) + fvm::div(phase.phiAlpha(), U)
- fvm::Sp(fvc::ddt(alpha) + fvc::div(phase.phiAlpha()), U) - fvm::Sp(fvc::ddt(alpha) + fvc::div(phase.phiAlpha()), U)
+ (alpha/phase.rho())*fluid.Cvm(phase)*
(
fvm::ddt(U)
+ fvm::div(phase.phi(), U)
- fvm::Sp(fvc::div(phase.phi()), U)
) )
- fvm::laplacian(alpha*nuEff, U) - fvm::laplacian(alpha*nuEff, U)
- fvc::div - fvc::div
( (