twoPhaseEulerFoam: added optional alphaMax to phaseModel for MULES limiter

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1209
This commit is contained in:
Henry
2014-03-07 23:32:00 +00:00
committed by Andrew Heather
parent cda70b5eb7
commit dd0efd159d
5 changed files with 20 additions and 5 deletions

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -63,6 +63,7 @@ Foam::phaseModel::phaseModel
( (
phaseProperties.subDict(name_) phaseProperties.subDict(name_)
), ),
alphaMax_(phaseDict_.lookupOrDefault("alphaMax", 1.0)),
thermo_(rhoThermo::New(fluid.mesh(), name_)), thermo_(rhoThermo::New(fluid.mesh(), name_)),
U_ U_
( (

View File

@ -2,7 +2,7 @@
========= | ========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | \\ / O peration |
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation \\ / A nd | Copyright (C) 2011-2014 OpenFOAM Foundation
\\/ M anipulation | \\/ M anipulation |
------------------------------------------------------------------------------- -------------------------------------------------------------------------------
License License
@ -71,6 +71,9 @@ class phaseModel
dictionary phaseDict_; dictionary phaseDict_;
//- Optional maximum phase-fraction (e.g. packing limit)
scalar alphaMax_;
//- Thermophysical properties //- Thermophysical properties
autoPtr<rhoThermo> thermo_; autoPtr<rhoThermo> thermo_;
@ -123,6 +126,13 @@ public:
//- Return the other phase in this two-phase system //- Return the other phase in this two-phase system
const phaseModel& otherPhase() const; const phaseModel& otherPhase() const;
//- Optional maximum phase-fraction (e.g. packing limit)
// Defaults to 1
scalar alphaMax() const
{
return alphaMax_;
}
//- Return the Sauter-mean diameter //- Return the Sauter-mean diameter
tmp<volScalarField> d() const; tmp<volScalarField> d() const;

View File

@ -508,7 +508,7 @@ void Foam::twoPhaseSystem::solve()
alphaPhic1, alphaPhic1,
Sp, Sp,
Su, Su,
1, phase1_.alphaMax(),
0 0
); );

View File

@ -24,6 +24,8 @@ particles
{ {
d 3e-4; d 3e-4;
} }
alphaMax 0.62;
} }
air air

View File

@ -24,6 +24,8 @@ particles
{ {
d 3e-4; d 3e-4;
} }
alphaMax 0.62;
} }
air air