48 int main(
int argc,
char *argv[])
52 "Transient solver for incompressible, laminar flow" 53 " of Newtonian fluids in liquid film formulation." 60 #include "readGravitationalAcceleration.H" 61 #include "readTransportProperties.H" 62 #include "createFaFields.H" 66 Info<<
"\nStarting time loop\n" <<
endl;
74 #include "setDeltaT.H" 93 *frictionFactor.internalField()
94 *
mag(
Us.internalField()),
108 Us.correctBoundaryConditions();
134 h.primitiveFieldRef() =
max 140 *
pos(
h0.value() -
h.primitiveField())
146 ps.correctBoundaryConditions();
150 Us.correctBoundaryConditions();
155 vsm.mapToVolume(
h,
H.boundaryFieldRef());
156 vsm.mapToVolume(
Us,
U.boundaryFieldRef());
const volSurfaceMapping vsm(aMesh)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
static tmp< GeometricField< Type, faePatchField, edgeMesh > > interpolate(const GeometricField< Type, faPatchField, areaMesh > &tvf, const edgeScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
constexpr char nl
The newline '\n' character (0x0a)
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< GeometricField< Type, faPatchField, areaMesh > > average(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Area-weighted average a edgeField creating a areaField.
Read the control parameters used by setDeltaT.
Calculates and outputs the mean and maximum Courant Numbers for the Finite Area method.
faMatrix< vector > faVectorMatrix
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
dimensionedScalar sigma("sigma", dimMass/sqr(dimTime), transportProperties)
dimensionedScalar pos(const dimensionedScalar &ds)
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
edgeScalarField phis(IOobject("phis", runTime.timeName(), aMesh.thisDb(), IOobject::NO_READ, IOobject::NO_WRITE), linearEdgeInterpolate(Us) &aMesh.Le())
tmp< GeometricField< Type, faePatchField, edgeMesh > > lnGrad(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
const dimensionedScalar h
Planck constant.
faMatrix< scalar > faScalarMatrix
dimensionedScalar rhol("rhol", dimDensity, transportProperties)
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
loopControl iters(runTime, aMesh.solutionDict(), "solution")
messageStream Info
Information stream (stdout output on master, null elsewhere)
volScalarField H(IOobject("H", runTime.timeName(), mesh.thisDb(), IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimLength, Zero))
Read the control parameters used by setDeltaT.
tmp< GeometricField< typename outerProduct< vector, Type >::type, faPatchField, areaMesh >> grad(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
GeometricField< scalar, faPatchField, areaMesh > areaScalarField