const Foam::volSymmTensorField& Foam::constSolidThermo::directionalK() const
{
dimensionedSymmTensor t
(
constK_.name(),
constK_.dimensions(),
symmTensor
(
constK_.value(),
0.0,
0.0,
constK_.value(),
0.0,
constK_.value()
)
);
return tmp<volSymmTensorField>
(
new volSymmTensorField
(
IOobject
(
"K",
mesh_.time().timeName(),
mesh_,
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh_,
t
)
);
}
XPP_SE::XPP_SE
(
const word& name,
const volScalarField& alpha,
const volVectorField& U,
const surfaceScalarField& phi,
const dictionary& dict
)
:
viscoelasticLaw(name, alpha, U, phi),
tau_
(
IOobject
(
"tau" + name,
U.time().timeName(),
U.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
U.mesh()
),
I_
(
dimensionedSymmTensor
(
"I",
dimensionSet(0, 0, 0, 0, 0, 0, 0),
symmTensor
(
1, 0, 0,
1, 0,
1
)
)
),
etaS1_(dict.subDict("phase1").lookup("etaS")),
etaS2_(dict.subDict("phase2").lookup("etaS")),
etaP1_(dict.subDict("phase1").lookup("etaP")),
etaP2_(dict.subDict("phase2").lookup("etaP")),
alpha1_(dict.subDict("phase1").lookup("alpha")),
alpha2_(dict.subDict("phase2").lookup("alpha")),
lambdaOb1_(dict.subDict("phase1").lookup("lambdaOb")),
lambdaOb2_(dict.subDict("phase2").lookup("lambdaOb")),
lambdaOs1_(dict.subDict("phase1").lookup("lambdaOs")),
lambdaOs2_(dict.subDict("phase2").lookup("lambdaOs")),
q1_(dict.subDict("phase1").lookup("q")),
q2_(dict.subDict("phase2").lookup("q")),
pt_
(
dimensionedScalar
(
"zero",
lambdaOb1_.dimensions(),
scalar( SMALL )
)
)
{}
Foam::S_MDCPP::S_MDCPP
(
const word& name,
const volVectorField& U,
const surfaceScalarField& phi,
const dictionary& dict
)
:
viscoelasticLaw(name, U, phi),
tau_
(
IOobject
(
"tau" + name,
U.time().timeName(),
U.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
U.mesh()
),
I_
(
dimensionedSymmTensor
(
"I",
dimensionSet(0, 0, 0, 0, 0, 0, 0),
symmTensor
(
1, 0, 0,
1, 0,
1
)
)
),
rho_(dict.lookup("rho")),
etaS_(dict.lookup("etaS")),
etaP_(dict.lookup("etaP")),
zeta_(dict.lookup("zeta")),
lambdaOb_(dict.lookup("lambdaOb")),
lambdaOs_(dict.lookup("lambdaOs")),
q_(dict.lookup("q"))
{}
Foam::XPP_DE::XPP_DE
(
const word& name,
const volVectorField& U,
const surfaceScalarField& phi,
const dictionary& dict
)
:
viscoelasticLaw(name, U, phi),
S_
(
IOobject
(
"S" + name,
U.time().timeName(),
U.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
U.mesh()
),
Lambda_
(
IOobject
(
"Lambda" + name,
U.time().timeName(),
U.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
U.mesh()
),
tau_
(
IOobject
(
"tau" + name,
U.time().timeName(),
U.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
U.mesh(),
dimensionedSymmTensor
(
"zero",
dimensionSet(1, -1, -2, 0, 0, 0, 0),
symmTensor::zero
)
),
I_
(
dimensionedSymmTensor
(
"I",
dimensionSet(0, 0, 0, 0, 0, 0, 0),
symmTensor
(
1, 0, 0,
1, 0,
1
)
)
),
rho_(dict.lookup("rho")),
etaS_(dict.lookup("etaS")),
etaP_(dict.lookup("etaP")),
alpha_(dict.lookup("alpha")),
lambdaOb_(dict.lookup("lambdaOb")),
lambdaOs_(dict.lookup("lambdaOs")),
q_(dict.lookup("q"))
{}
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;
}
}
Foam::Leonov::Leonov
(
const word& name,
const volScalarField& alpha,
const volVectorField& U,
const surfaceScalarField& phi,
const dictionary& dict
)
:
viscoelasticLaw(name, alpha, U, phi),
sigma_
(
IOobject
(
"sigma" + name,
U.time().timeName(),
U.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
U.mesh()
),
tau_
(
IOobject
(
"tau" + name,
U.time().timeName(),
U.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
U.mesh(),
dimensionedSymmTensor
(
"zero",
dimensionSet(1, -1, -2, 0, 0, 0, 0),
symmTensor::zero
)
),
I_
(
dimensionedSymmTensor
(
"I",
dimensionSet(0, 0, 0, 0, 0, 0, 0),
symmTensor
(
1, 0, 0,
1, 0,
1
)
)
),
rho1_(dict.subDict("phase1").lookup("rho")),
rho2_(dict.subDict("phase2").lookup("rho")),
etaS1_(dict.subDict("phase1").lookup("etaS")),
etaS2_(dict.subDict("phase2").lookup("etaS")),
etaP1_(dict.subDict("phase1").lookup("etaP")),
etaP2_(dict.subDict("phase2").lookup("etaP")),
lambda1_(dict.subDict("phase1").lookup("lambda")),
lambda2_(dict.subDict("phase2").lookup("lambda"))
{}
