30 template<
class RhoFieldType>
31 void Foam::porosityModels::powerLawLopesdaCosta::apply
35 const RhoFieldType&
rho,
39 const scalar C1m1b2 = (C1_ - 1.0)/2.0;
57 template<
class RhoFieldType>
58 void Foam::porosityModels::powerLawLopesdaCosta::apply
61 const RhoFieldType&
rho,
65 const scalar C1m1b2 = (C1_ - 1.0)/2.0;
67 forAll(cellZoneIDs_, zonei)
79 0.5*
rho[celli]*Cd_*Av_[i]
#define forAll(list, i)
Loop across all elements in list.
const meshCellZones & cellZones() const
Return cell zones.
const fvMesh & mesh_
Reference to the mesh database.
labelList cellZoneIDs_
Cell zone IDs.
scalarField Av_
Porosity surface area per unit volume zone field.
volScalarField scalarField(fieldObject, mesh)
volVectorField vectorField(fieldObject, mesh)
List< label > labelList
A List of labels.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
static const Identity< scalar > I
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Field< vector > vectorField
Specialisation of Field<T> for vector.
Field< tensor > tensorField
Specialisation of Field<T> for tensor.
dimensioned< scalar > magSqr(const dimensioned< Type > &)