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

multiphaseEulerFoam: Improved U correction

Browse Source
This commit is contained in:
Henry
2012-03-12 15:28:41 +00:00
parent 767b768ddc
commit 05546f46d2
1 changed files with 1 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H
View File

@ -35,7 +35,6 @@
{ {
phaseModel& phase = iter(); phaseModel& phase = iter();
phase.phi().oldTime();
mrfZones.absoluteFlux(phase.phi().oldTime()); mrfZones.absoluteFlux(phase.phi().oldTime());
mrfZones.absoluteFlux(phase.phi()); mrfZones.absoluteFlux(phase.phi());
@ -173,7 +172,6 @@
phase.phi() = phase.phi() =
phiHbyAs[phasei] phiHbyAs[phasei]
+ rAlphaAUfs[phasei]*mSfGradp/phase.rho(); + rAlphaAUfs[phasei]*mSfGradp/phase.rho();
mrfZones.relativeFlux(phase.phi());
mrfZones.relativeFlux(phase.phi().oldTime()); mrfZones.relativeFlux(phase.phi().oldTime());
phi += alphafs[phasei]*phase.phi(); phi += alphafs[phasei]*phase.phi();
@ -206,7 +204,7 @@
+ rAlphaAUfs[phasei]*mSfGradp/phase.rho() + rAlphaAUfs[phasei]*mSfGradp/phase.rho()
); );
phase.U() = fvc::reconstruct(phase.phi()); // phase.U() = fvc::reconstruct(phase.phi());
phase.U().correctBoundaryConditions(); phase.U().correctBoundaryConditions();
U += alpha*phase.U(); U += alpha*phase.U();