32template<
class RhoFieldType>
37 const RhoFieldType&
rho,
41 const scalar C1m1b2 = (C1_ - 1.0)/2.0;
49 const label celli =
cells[i];
59template<
class RhoFieldType>
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]
virtual void apply()=0
Apply bins.
const cellZoneMesh & cellZones() const noexcept
Return cell zone mesh.
const fvMesh & mesh_
Reference to the mesh database.
labelList cellZoneIDs_
Cell zone IDs.
scalarField Sigma_
Porosity surface area per unit volume zone field.
List< label > labelList
A List of labels.
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< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
#define forAll(list, i)
Loop across all elements in list.