19 Cs = Cs0*(1.0 +
R.component(vector::X)/
mag(
R));
42 Uinf.value()*(0.25*(3.0 +
sqr(
R[faceI].
x()/
mag(
R[faceI]))) - 1.0);
45 Uinf.value()*0.25*
R[faceI].x()*
R[faceI].y()/
sqr(
mag(
R[faceI]));
48 Uinf.value()*0.25*
R[faceI].x()*
R[faceI].z()/
sqr(
mag(
R[faceI]));
50 Us.boundaryFieldRef().evaluateCoupled<coupledFaPatch>();
52 Us -= aMesh.faceAreaNormals()*(aMesh.faceAreaNormals() &
Us);
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
const dimensionSet dimViscosity
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
Ostream & endl(Ostream &os)
Add newline and flush stream.
Info<< "Reading field Cs"<< endl;areaScalarField Cs(IOobject("Cs", runTime.timeName(), aMesh.thisDb(), IOobject::MUST_READ, IOobject::AUTO_WRITE), aMesh);dimensionedScalar Cs0("Cs0", dimMass/dimArea, 1.0);const areaVectorField &R=aMesh.areaCentres();Cs=Cs0 *(1.0+R.component(vector::X)/mag(R));dimensionedScalar Ds("Ds", dimViscosity, 1.0);areaVectorField Us(IOobject("Us", runTime.timeName(), aMesh.thisDb(), IOobject::NO_READ, IOobject::NO_WRITE), aMesh, dimensionedVector(dimVelocity, Zero));dimensioned< scalar > Uinf("Uinf", dimVelocity, 1.0)
tmp< GeometricField< Type, faePatchField, edgeMesh > > linearEdgeInterpolate(const GeometricField< Type, faPatchField, areaMesh > &vf)
GeometricField< scalar, faePatchField, edgeMesh > edgeScalarField
edgeScalarField phis(IOobject("phis", runTime.timeName(), aMesh.thisDb(), IOobject::NO_READ, IOobject::NO_WRITE), linearEdgeInterpolate(Us) &aMesh.Le())
#define R(A, B, C, D, E, F, K, M)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
messageStream Info
Information stream (stdout output on master, null elsewhere)
const dimensionSet dimMass(1, 0, 0, 0, 0, 0, 0)
GeometricField< vector, faPatchField, areaMesh > areaVectorField
const dimensionSet dimArea(sqr(dimLength))
GeometricField< scalar, faPatchField, areaMesh > areaScalarField
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity