32 void Foam::solvers::multiphaseEuler::setRDeltaT()
40 pimpleDict.lookupOrDefault<scalar>(
"maxCo", 0.2)
43 const scalar rDeltaTSmoothingCoeff
45 pimpleDict.lookupOrDefault<scalar>(
"rDeltaTSmoothingCoeff", 0.02)
56 rDeltaT.internalFieldRef() =
60 scalar minRDeltaT =
gMin(rDeltaT.primitiveField());
61 if (pimpleDict.found(
"maxDeltaT") || minRDeltaT < rootVSmall)
63 const scalar clipRDeltaT = 1/pimpleDict.lookup<scalar>(
"maxDeltaT");
64 rDeltaT.max(clipRDeltaT);
65 minRDeltaT =
max(minRDeltaT, clipRDeltaT);
67 if (pimpleDict.found(
"minDeltaT"))
69 const scalar clipRDeltaT = 1/pimpleDict.lookup<scalar>(
"minDeltaT");
70 rDeltaT.min(clipRDeltaT);
71 minRDeltaT =
min(minRDeltaT, clipRDeltaT);
74 Info<<
"Flow time scale min/max = "
75 <<
gMin(1/rDeltaT.primitiveField()) <<
", " << 1/minRDeltaT <<
endl;
78 rDeltaT.correctBoundaryConditions();
82 Info<<
"Smoothed flow time scale min/max = "
83 <<
gMin(1/rDeltaT.primitiveField())
84 <<
", " <<
gMax(1/rDeltaT.primitiveField()) <<
endl;
#define forAll(list, i)
Loop across all elements in list.
const DimensionedField< scalar, volMesh > & V() const
Return cell volumes.
virtual const dictionary & dict() const
Return the solution dictionary.
pimpleNoLoopControl pimple
PIMPLE inner-loop controls.
const fvMesh & mesh
Region mesh.
const surfaceScalarField & phi
Reference to the mass-flux field.
tmp< volScalarField > trDeltaT
Optional LTS reciprocal time-step field.
const phaseSystem::phaseModelPartialList & movingPhases
Reference to the moving phases.
Switch faceMomentum
Cell/face momentum equation switch.
tmp< surfaceScalarField > trDeltaTf
Optional LTS reciprocal face time-step field.
Provides functions smooth spread and sweep which use the FvFaceCellWave algorithm to smooth and redis...
Surface integrate surfaceField creating a volField. Surface sum a surfaceField creating a volField.
tmp< VolField< Type > > surfaceSum(const SurfaceField< Type > &ssf)
static tmp< SurfaceField< Type > > interpolate(const VolField< Type > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
void smooth(volScalarField &field, const scalar coeff)
Ostream & endl(Ostream &os)
Add newline and flush stream.
SurfaceField< scalar > surfaceScalarField
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
VolField< scalar > volScalarField
dimensioned< scalar > mag(const dimensioned< Type > &)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
Type gMin(const FieldField< Field, Type > &f)
Type gMax(const FieldField< Field, Type > &f)