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

multiphaseInterFoam: No need for refPhase with MULES

Browse Source
This commit is contained in:
Henry
2011-09-05 19:09:32 +01:00
parent d47625903e
commit cff638a7da
2 changed files with 15 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C
View File

@ -64,7 +64,6 @@ Foam::multiphaseMixture::multiphaseMixture
:
transportModel(U, phi),
phases_(lookup("phases"), phase::iNew(U, phi)),
refPhase_(*phases_.lookup(word(lookup("refPhase")))),
mesh_(U.mesh()),
U_(U),
@ -335,7 +334,7 @@ void Foam::multiphaseMixture::correctContactAngle
) const
{
const volScalarField::GeometricBoundaryField& gbf
= refPhase_.boundaryField();
= alpha1.boundaryField();
const fvBoundaryMesh& boundary = mesh_.boundary();
@ -549,8 +548,20 @@ void Foam::multiphaseMixture::solveAlphas
MULES::limitSum(phiAlphaCorrs);
rhoPhi_ = 0.0*phi_*refPhase_.rho();
volScalarField sumAlpha("sumAlpha", 0.0*refPhase_);
rhoPhi_ = dimensionedScalar("0", dimensionSet(1, 0, -1, 0, 0), 0);
volScalarField sumAlpha
(
IOobject
(
"sumAlpha",
mesh_.time().timeName(),
mesh_
),
mesh_,
dimensionedScalar("sumAlpha", dimless, 0)
);
phasei = 0;
forAllIter(PtrDictionary<phase>, phases_, iter)

4
applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.H
View File

@ -136,10 +136,6 @@ private:
//- Dictionary of phases
PtrDictionary<phase> phases_;
//- The phase chosen as reference, the one which is derived from
// the others such thatr they sum to 1
phase& refPhase_;
const fvMesh& mesh_;
const volVectorField& U_;
const surfaceScalarField& phi_;