34 template<
class RhoFieldType>
35 void Foam::MRFZone::relativeRhoFlux
37 const RhoFieldType&
rho,
50 label facei = internalFaces_[i];
51 phi[facei] -= rho[facei]*(Omega ^ (Cf[facei] - origin)) & Sf[facei];
59 label patchFacei = includedFaces_[
patchi][i];
61 phi.boundaryField()[
patchi][patchFacei] = 0.0;
70 label patchFacei = excludedFaces_[
patchi][i];
72 phi.boundaryField()[
patchi][patchFacei] -=
73 rho.boundaryField()[
patchi][patchFacei]
74 * (Omega ^ (Cf.boundaryField()[
patchi][patchFacei] - origin))
75 & Sf.boundaryField()[
patchi][patchFacei];
81 template<
class RhoFieldType>
82 void Foam::MRFZone::absoluteRhoFlux
84 const RhoFieldType& rho,
91 const vector& origin = origin_.value();
92 const vector& Omega = Omega_.value();
97 label facei = internalFaces_[i];
98 phi[facei] += rho[facei]*(Omega ^ (Cf[facei] - origin)) & Sf[facei];
102 forAll(includedFaces_, patchi)
104 forAll(includedFaces_[patchi], i)
106 label patchFacei = includedFaces_[
patchi][i];
108 phi.boundaryField()[
patchi][patchFacei] +=
109 rho.boundaryField()[
patchi][patchFacei]
110 * (Omega ^ (Cf.boundaryField()[
patchi][patchFacei] - origin))
111 & Sf.boundaryField()[
patchi][patchFacei];
116 forAll(excludedFaces_, patchi)
118 forAll(excludedFaces_[patchi], i)
120 label patchFacei = excludedFaces_[
patchi][i];
122 phi.boundaryField()[
patchi][patchFacei] +=
123 rho.boundaryField()[
patchi][patchFacei]
124 * (Omega ^ (Cf.boundaryField()[
patchi][patchFacei] - origin))
125 & Sf.boundaryField()[
patchi][patchFacei];