10 const word porosityModel(
"JensenEtAl2014");
28 DarcyForchheimerCoeff +=
rho*
mag(
U)
29 *
beta*(1 +
pos(KC)*7.5/KC)*(1 - porosity)/
sqr(porosity)/d50;
47 UEqn *= scalar(1)/porosity;
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
tmp< volScalarField > tporosity
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
const dimensionSet dimless
Dimensionless.
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensionedScalar pos(const dimensionedScalar &ds)
IOdictionary porosityProperties(IOobject("porosityProperties", runTime.constant(), runTime, IOobject::READ_IF_PRESENT, IOobject::NO_WRITE))
const bool porosityEnabled(porosityProperties.getOrDefault< bool >("porosityEnabled", false))
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
const dimensionSet dimArea(sqr(dimLength))