twoPhaseEulerFoam: added optional alphaMax to phaseModel for MULES limiter

This commit is contained in:
Henry
2014-03-07 21:54:35 +00:00
parent e43141e314
commit decd96411a
6 changed files with 21 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

@ -248,7 +248,7 @@ Foam::twoPhaseSystem::twoPhaseSystem
pair2In1_ pair2In1_
) )
); );
wallLubrication_.set wallLubrication_.set
( (
new BlendedInterfacialModel<wallLubricationModel> new BlendedInterfacialModel<wallLubricationModel>
@ -264,7 +264,7 @@ Foam::twoPhaseSystem::twoPhaseSystem
pair2In1_ pair2In1_
) )
); );
turbulentDispersion_.set turbulentDispersion_.set
( (
new BlendedInterfacialModel<turbulentDispersionModel> new BlendedInterfacialModel<turbulentDispersionModel>
@ -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

View File

@ -32,6 +32,7 @@ FoamFile
walls walls
{ {
type wall; type wall;
inGroups 1(wall);
nFaces 400; nFaces 400;
startFace 11830; startFace 11830;
} }