32 template<
class RhoFieldType>
33 void Foam::porosityModels::powerLawLopesdaCosta::apply
37 const RhoFieldType&
rho,
41 const scalar C1m1b2 = (C1_ - 1.0)/2.0;
49 const label celli =
cells[i];
59 template<
class RhoFieldType>
60 void Foam::porosityModels::powerLawLopesdaCosta::apply
63 const RhoFieldType&
rho,
67 const scalar C1m1b2 = (C1_ - 1.0)/2.0;
69 forAll(cellZoneIDs_, zonei)
75 const label celli =
cells[i];
81 0.5*
rho[celli]*Cd_*Sigma_[i]
#define forAll(list, i)
Loop across all elements in list.
static const Identity< scalar > I
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
const fvMesh & mesh_
Reference to the mesh database.
scalarField Sigma_
Porosity surface area per unit volume zone field.
labelList cellZoneIDs_
Cell zone IDs.
Field< tensor > tensorField
Specialisation of Field<T> for tensor.
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
const cellZoneMesh & cellZones() const noexcept
Return cell zone mesh.
Field< vector > vectorField
Specialisation of Field<T> for vector.
List< label > labelList
A List of labels.
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)