Foam::displacementLinearMotionMotionSolver::
displacementLinearMotionMotionSolver
(
const polyMesh& mesh,
const dictionary& dict
)
:
points0MotionSolver(mesh, dict, typeName),
axis_(normalised(vector(coeffDict().lookup("axis")))),
xFixed_(readScalar(coeffDict().lookup("xFixed"))),
xMoving_(readScalar(coeffDict().lookup("xMoving"))),
displacement_(Function1<scalar>::New("displacement", coeffDict()))
{}
bool dynLagrangianCsBound::read()
{
if (GenEddyVisc::read())
{
filter_.read(coeffDict());
theta_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool lowReOneEqEddy::read()
{
if (GenEddyVisc::read())
{
ck_.readIfPresent(coeffDict());
beta_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool GenEddyVisc::read()
{
if (LESModel::read())
{
ce_.readIfPresent(coeffDict());
Prt_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool SpalartAllmarasIDDES::read()
{
if (SpalartAllmaras::read())
{
fwStar_.readIfPresent(coeffDict());
cl_.readIfPresent(coeffDict());
ct_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
oneEqEddy::oneEqEddy
(
const volVectorField& U,
const surfaceScalarField& phi,
transportModel& transport
)
:
LESModel(typeName, U, phi, transport),
GenEddyVisc(U, phi, transport),
k_
(
IOobject
(
"k",
runTime_.timeName(),
mesh_,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh_
),
ck_
(
dimensioned<scalar>::lookupOrAddToDict
(
"ck",
coeffDict(),
0.094
)
)
{
printCoeffs();
}
void Foam::sixDoFRigidBodyMotionSolver::solve()
{
const Time& t = mesh().time();
if (mesh().nPoints() != points0().size())
{
FatalErrorInFunction
<< "The number of points in the mesh seems to have changed." << endl
<< "In constant/polyMesh there are " << points0().size()
<< " points; in the current mesh there are " << mesh().nPoints()
<< " points." << exit(FatalError);
}
// Store the motion state at the beginning of the time-stepbool
bool firstIter = false;
if (curTimeIndex_ != t.timeIndex())
{
motion_.newTime();
curTimeIndex_ = t.timeIndex();
firstIter = true;
}
dimensionedVector g("g", dimAcceleration, Zero);
if (mesh().foundObject<uniformDimensionedVectorField>("g"))
{
g = mesh().lookupObject<uniformDimensionedVectorField>("g");
}
else if (coeffDict().found("g"))
{
coeffDict().lookup("g") >> g;
}
locDynOneEqEddy::locDynOneEqEddy
(
const volVectorField& U,
const surfaceScalarField& phi,
transportModel& transport
)
:
LESModel(typeName, U, phi, transport),
GenEddyVisc(U, phi, transport),
k_
(
IOobject
(
"k",
runTime_.timeName(),
U_.db(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh_
),
simpleFilter_(U.mesh()),
filterPtr_(LESfilter::New(U.mesh(), coeffDict())),
filter_(filterPtr_())
{
volScalarField KK = 0.5*(filter_(magSqr(U)) - magSqr(filter_(U)));
updateSubGridScaleFields(symm(fvc::grad(U)), KK);
printCoeffs();
}
bool kEpsilon::read()
{
if (RASModel::read())
{
Cmu_.readIfPresent(coeffDict());
C1_.readIfPresent(coeffDict());
C2_.readIfPresent(coeffDict());
sigmaEps_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool LaunderSharmaKE::read()
{
if (RASModel::read())
{
Cmu_.readIfPresent(coeffDict());
C1_.readIfPresent(coeffDict());
C2_.readIfPresent(coeffDict());
sigmaEps_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool kOmega::read()
{
if (RASModel::read())
{
Cmu_.readIfPresent(coeffDict());
beta_.readIfPresent(coeffDict());
alphaK_.readIfPresent(coeffDict());
alphaOmega_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
Foam::displacementSBRStressFvMotionSolver::displacementSBRStressFvMotionSolver
(
const polyMesh& mesh,
const IOdictionary& dict
)
:
displacementMotionSolver(mesh, dict, dict.lookup("solver")),
fvMotionSolver(mesh),
cellDisplacement_
(
IOobject
(
"cellDisplacement",
mesh.time().timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
fvMesh_,
dimensionedVector
(
"cellDisplacement",
pointDisplacement().dimensions(),
Zero
),
cellMotionBoundaryTypes<vector>(pointDisplacement().boundaryField())
),
diffusivityPtr_
(
motionDiffusivity::New(fvMesh_, coeffDict().lookup("diffusivity"))
)
{}
Foam::velocityLaplacianFvMotionSolver::velocityLaplacianFvMotionSolver
(
const polyMesh& mesh,
const IOdictionary& dict
)
:
velocityMotionSolver(mesh, dict, typeName),
fvMotionSolverCore(mesh),
cellMotionU_
(
IOobject
(
"cellMotionU",
mesh.time().timeName(),
mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
fvMesh_,
dimensionedVector
(
"cellMotionU",
pointMotionU_.dimensions(),
Zero
),
cellMotionBoundaryTypes<vector>(pointMotionU_.boundaryField())
),
diffusivityPtr_
(
motionDiffusivity::New(fvMesh_, coeffDict().lookup("diffusivity"))
)
{}
bool realizableKE::read()
{
if (RASModel::read())
{
Cmu_.readIfPresent(coeffDict());
A0_.readIfPresent(coeffDict());
C2_.readIfPresent(coeffDict());
sigmak_.readIfPresent(coeffDict());
sigmaEps_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool spectEddyVisc::read()
{
if (GenEddyVisc::read())
{
cB_.readIfPresent(coeffDict());
cK1_.readIfPresent(coeffDict());
cK2_.readIfPresent(coeffDict());
cK3_.readIfPresent(coeffDict());
cK4_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool qZeta::read()
{
if (RASModel::read())
{
Cmu_.readIfPresent(coeffDict());
C1_.readIfPresent(coeffDict());
C2_.readIfPresent(coeffDict());
sigmaZeta_.readIfPresent(coeffDict());
anisotropic_.readIfPresent("anisotropic", coeffDict());
return true;
}
else
{
return false;
}
}
Foam::externalDisplacementMeshMover&
Foam::displacementMeshMoverMotionSolver::meshMover() const
{
if (!meshMoverPtr_.valid())
{
const word moverType(coeffDict().lookup("meshMover"));
meshMoverPtr_ = externalDisplacementMeshMover::New
(
moverType,
coeffDict().optionalSubDict(moverType + "Coeffs"),
localPointRegion::findDuplicateFacePairs(mesh()),
pointDisplacement_
);
}
return meshMoverPtr_();
}
rampHoldFall::rampHoldFall
(
const dictionary& relaxationDict,
const Time& runTime
)
:
relaxationModel(typeName, relaxationDict, runTime),
rampStartRelaxation_(readScalar(coeffDict().lookup("rampStartRelaxation"))),
holdRelaxation_(readScalar(coeffDict().lookup("holdRelaxation"))),
fallEndRelaxation_(readScalar(coeffDict().lookup("fallEndRelaxation"))),
rampEndFraction_(readScalar(coeffDict().lookup("rampEndFraction"))),
fallStartFraction_(readScalar(coeffDict().lookup("fallStartFraction"))),
rampGradient_((holdRelaxation_ - rampStartRelaxation_)/(rampEndFraction_)),
fallGradient_
(
(fallEndRelaxation_ - holdRelaxation_)/(1 - fallStartFraction_)
)
{}
Foam::motionDiffusivity&
Foam::displacementLaplacianFvMotionSolver::diffusivity()
{
if (!diffusivityPtr_.valid())
{
diffusivityPtr_ = motionDiffusivity::New
(
fvMesh_,
coeffDict().lookup("diffusivity")
);
}
return diffusivityPtr_();
}
bool LienCubicKE::read()
{
if (RASModel::read())
{
C1_.readIfPresent(coeffDict());
C2_.readIfPresent(coeffDict());
sigmak_.readIfPresent(coeffDict());
sigmaEps_.readIfPresent(coeffDict());
A1_.readIfPresent(coeffDict());
A2_.readIfPresent(coeffDict());
Ctau1_.readIfPresent(coeffDict());
Ctau2_.readIfPresent(coeffDict());
Ctau3_.readIfPresent(coeffDict());
alphaKsi_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool Smagorinsky::read()
{
if (GenEddyVisc::read())
{
ck_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
scaleSimilarity::scaleSimilarity
(
const volVectorField& U,
const surfaceScalarField& phi,
transportModel& transport
)
:
LESModel(typeName, U, phi, transport),
filterPtr_(LESfilter::New(U.mesh(), coeffDict())),
filter_(filterPtr_())
{
printCoeffs();
}
bool Wale::read()
{
if (GenEddyVisc::read())
{
cw_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool GenSGSStress::read()
{
if (LESModel::read())
{
ce_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool Smagorinsky2::read()
{
if (Smagorinsky::read())
{
cD2.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
bool scaleSimilarity::read()
{
if (LESModel::read())
{
filter_.read(coeffDict());
return true;
}
else
{
return false;
}
}
bool locDynOneEqEddy::read()
{
if (GenEddyVisc::read())
{
filter_.read(coeffDict());
return true;
}
else
{
return false;
}
}
bool GenSGSStress::read()
{
if (LESModel::read())
{
ce_.readIfPresent(coeffDict());
couplingFactor_.readIfPresent(coeffDict());
if (couplingFactor_.value() < 0.0 || couplingFactor_.value() > 1.0)
{
FatalErrorIn("GenSGSStress::read()")
<< "couplingFactor = " << couplingFactor_
<< " is not in range 0 - 1"
<< exit(FatalError);
}
return true;
}
else
{
return false;
}
}
bool tatunRNG::read()
{
if (RASModel::read())
{
Cmu_.readIfPresent(coeffDict());
C1_.readIfPresent(coeffDict());
C2_.readIfPresent(coeffDict());
C3_.readIfPresent(coeffDict());
Prt_.readIfPresent(coeffDict());
sigmaEps_.readIfPresent(coeffDict());
Prt_.readIfPresent(coeffDict());
eta0_.readIfPresent(coeffDict());
beta_.readIfPresent(coeffDict());
return true;
}
else
{
return false;
}
}
void Foam::velocityLaplacianFvMotionSolver::updateMesh
(
const mapPolyMesh& mpm
)
{
velocityMotionSolver::updateMesh(mpm);
// Update diffusivity. Note two stage to make sure old one is de-registered
// before creating/registering new one.
diffusivityPtr_.reset(NULL);
diffusivityPtr_ = motionDiffusivity::New
(
fvMesh_,
coeffDict().lookup("diffusivity")
);
}