Foam::kineticTheoryModel::kineticTheoryModel
(
const Foam::phaseModel& phasea,
const Foam::volVectorField& Ub,
const Foam::volScalarField& alpha,
const Foam::dragModel& draga
)
:
phasea_(phasea),
Ua_(phasea.U()),
Ub_(Ub),
alpha_(alpha),
phia_(phasea.phi()),
draga_(draga),
rhoa_(phasea.rho()),
da_(phasea.d()),
nua_(phasea.nu()),
kineticTheoryProperties_
(
IOobject
(
"kineticTheoryProperties",
Ua_.time().constant(),
Ua_.mesh(),
IOobject::MUST_READ,
IOobject::NO_WRITE
)
),
kineticTheory_(kineticTheoryProperties_.lookup("kineticTheory")),
equilibrium_(kineticTheoryProperties_.lookup("equilibrium")),
viscosityModel_
(
kineticTheoryModels::viscosityModel::New
(
kineticTheoryProperties_
)
),
conductivityModel_
(
conductivityModel::New
(
kineticTheoryProperties_
)
),
radialModel_
(
radialModel::New
(
kineticTheoryProperties_
)
),
granularPressureModel_
(
granularPressureModel::New
(
kineticTheoryProperties_
)
),
frictionalStressModel_
(
frictionalStressModel::New
(
kineticTheoryProperties_
)
),
e_(kineticTheoryProperties_.lookup("e")),
alphaMax_(kineticTheoryProperties_.lookup("alphaMax")),
alphaMinFriction_(kineticTheoryProperties_.lookup("alphaMinFriction")),
Fr_(kineticTheoryProperties_.lookup("Fr")),
eta_(kineticTheoryProperties_.lookup("eta")),
p_(kineticTheoryProperties_.lookup("p")),
phi_(dimensionedScalar(kineticTheoryProperties_.lookup("phi"))*M_PI/180.0),
Theta_
(
IOobject
(
"Theta",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh()
),
mua_
(
IOobject
(
"mua",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0)
),
lambda_
(
IOobject
(
"lambda",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0)
),
pa_
(
IOobject
(
"pa",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0)
),
kappa_
(
IOobject
(
"kappa",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0)
),
gs0_
(
IOobject
(
"gs0",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 1.0)
)
{}
Foam::kineticTheoryModel::kineticTheoryModel
(
const Foam::phaseModel& phasea,
const Foam::volVectorField& Ub,
const Foam::volScalarField& alpha,
const Foam::dragModel& draga
)
:
phasea_(phasea),
Ua_(phasea.U()),
Ub_(Ub),
alpha_(alpha),
phia_(phasea.phi()),
draga_(draga),
rhoa_(phasea.rho()),
da_(phasea.d()),
nua_(phasea.nu()),
kineticTheoryProperties_
(
IOobject
(
"kineticTheoryProperties",
Ua_.time().constant(),
Ua_.mesh(),
IOobject::MUST_READ,
IOobject::NO_WRITE
)
),
kineticTheory_(kineticTheoryProperties_.lookup("kineticTheory")),
equilibrium_(kineticTheoryProperties_.lookup("equilibrium")),
viscosityModel_
(
kineticTheoryModels::viscosityModel::New
(
kineticTheoryProperties_
)
),
conductivityModel_
(
conductivityModel::New
(
kineticTheoryProperties_
)
),
radialModel_
(
radialModel::New
(
kineticTheoryProperties_
)
),
granularPressureModel_
(
granularPressureModel::New
(
kineticTheoryProperties_
)
),
frictionalStressModel_
(
frictionalStressModel::New
(
kineticTheoryProperties_
)
),
e_(kineticTheoryProperties_.lookup("e")),
alphaMax_(kineticTheoryProperties_.lookup("alphaMax")),
alphaMinFriction_(kineticTheoryProperties_.lookup("alphaMinFriction")),
Fr_(kineticTheoryProperties_.lookup("Fr")),
eta_(kineticTheoryProperties_.lookup("eta")),
p_(kineticTheoryProperties_.lookup("p")),
phi_(dimensionedScalar(kineticTheoryProperties_.lookup("phi"))*M_PI/180.0),
Theta_
(
IOobject
(
"Theta",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh()
),
mua_
(
IOobject
(
"mua",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0)
),
lambda_
(
IOobject
(
"lambda",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0)
),
pa_
(
IOobject
(
"pa",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0)
),
ppMagf_
(
IOobject
(
"ppMagf",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0)
),
kappa_
(
IOobject
(
"kappa",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -1, 0, 0), 0.0)
),
gs0_
(
IOobject
(
"gs0",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 1.0)
),
gs0Prime_
(
IOobject
(
"gs0prime",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 0.0)
),
//_AO_09/01/2014
// Shear stress ratio eta (here, we call it as upsilon)
upsilon_
(
IOobject
(
"upsilon",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(0, 0, 0, 0, 0), 0.0)
),
//_AO_09/01/2014
// Shear stress tau_ (here, we call it as upsilon)
tau_
(
IOobject
(
"tau",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedSymmTensor("zero", dimensionSet(1, -1, -2, 0, 0), symmTensor(0,0,0,0,0,0))
),
//_AO_09/01/2014
// Effective restitution coefficient (function of e and mu (friction coefficient) )
eEff_(0),
// Friction coefficient
muFric_(0),
// Dilute inertial regime 0 < alpha < alphaf (p. 2) alphaf ~ 0.49
alphaf_(0),
// Dense intertial regime
alphac_(0),
// YG 12/27/2014
alphad_(0),
// Yield stress ratio
upsilons_(0),
// Modified kinetic theory by Chialvo-Sundaresan on/off
mofidiedKineticTheoryPU_(kineticTheoryProperties_.lookup("modifiedKineticTheoryPU")),
//
//-AO, YG - Decompose particle pressure, Sundar's idea
decomposePp_(false),
//The option to use Berzi's Model in the SKT implementation
Berzi_(false),
paStar_
(
IOobject
(
"paStar",
Ua_.time().timeName(),
Ua_.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
Ua_.mesh(),
dimensionedScalar("zero", dimensionSet(1, -1, -2, 0, 0), 0.0)
)
{
if(kineticTheoryProperties_.found("Berzi"))
{
Berzi_ = true;
}
if(kineticTheoryProperties_.found("decomposePp"))
{
decomposePp_ = true;
}
}