zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

MULES: Corrected the handling of density variation for the static mesh case

Browse Source
This commit is contained in:
Henry
2011-03-12 15:49:38 +00:00
parent f404420f88
commit 448cd8b97b
1 changed files with 14 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
View File

@ -457,12 +457,20 @@ void Foam::MULES::limiter
tmp<volScalarField::DimensionedInternalField> V0 = mesh.Vsc0();
psiMaxn =
V*((rho.field()/deltaT - Sp.field())*psiMaxn - Su.field())
V
*(
(rho.field()/deltaT - Sp.field())*psiMaxn
- Su.field()
)
- (V0().field()/deltaT)*rho.oldTime().field()*psi0
+ sumPhiBD;
psiMinn =
V*(Su.field() - (rho.field()/deltaT - Sp.field())*psiMinn)
V
*(
Su.field()
- (rho.field()/deltaT - Sp.field())*psiMinn
)
+ (V0().field()/deltaT)*rho.oldTime().field()*psi0
- sumPhiBD;
}
@ -472,16 +480,17 @@ void Foam::MULES::limiter
V
*(
(rho.field()/deltaT - Sp.field())*psiMaxn
- (rho.oldTime().field()/deltaT)*psi0
- Su.field()
- (rho.oldTime().field()/deltaT)*psi0
)
+ sumPhiBD;
psiMinn =
V
*(
(rho.field()/deltaT)*psi0
- (rho.oldTime().field()/deltaT - Sp.field())*psiMinn + Su.field()
Su.field()
- (rho.field()/deltaT - Sp.field())*psiMinn
+ (rho.oldTime().field()/deltaT)*psi0
)
- sumPhiBD;
}