diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H b/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H index 87a539ba75..15acdcacb7 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H +++ b/applications/solvers/multiphase/multiphaseEulerFoam/createFields.H @@ -37,7 +37,7 @@ surfaceScalarField phi IOobject::AUTO_WRITE ), mesh, - dimensionedScalar(dimArea*dimVelocity, Zero) + dimensionedScalar(dimVelocity*dimArea, Zero) ); multiphaseSystem fluid(U, phi); diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H b/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H index 4d3b927e4e..f5c7886eb2 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H +++ b/applications/solvers/multiphase/multiphaseEulerFoam/pEqn.H @@ -53,7 +53,7 @@ IOobject::AUTO_WRITE ), mesh, - dimensionedScalar(dimArea*dimVelocity, Zero) + dimensionedScalar(dimVelocity*dimArea, Zero) ); volScalarField rho("rho", fluid.rho()); diff --git a/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pU/pEqn.H b/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pU/pEqn.H index 5e488110b0..5f65eb5129 100644 --- a/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pU/pEqn.H +++ b/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pU/pEqn.H @@ -180,7 +180,7 @@ while (pimple.correct()) IOobject::AUTO_WRITE ), mesh, - dimensionedScalar("phiHbyA", dimArea*dimVelocity, 0) + dimensionedScalar(dimVelocity*dimArea, Zero) ); forAll(phases, phasei) diff --git a/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pUf/pEqn.H b/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pUf/pEqn.H index 98da1ded9b..f2e3d51067 100644 --- a/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pUf/pEqn.H +++ b/applications/solvers/multiphase/reactingMultiphaseEulerFoam/pUf/pEqn.H @@ -165,7 +165,7 @@ while (pimple.correct()) IOobject::AUTO_WRITE ), mesh, - dimensionedScalar("phiHbyA", dimArea*dimVelocity, 0) + dimensionedScalar(dimVelocity*dimArea, Zero) ); forAll(phases, phasei) diff --git a/src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C b/src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C index 4e40f829f7..d49596d030 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/advective/advectiveFvPatchField.C @@ -166,7 +166,7 @@ Foam::advectiveFvPatchField::advectionSpeed() const phiName_ ); - if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + if (phi.dimensions() == dimMass/dimTime) { const fvPatchScalarField& rhop = this->patch().template lookupPatchField diff --git a/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchFields.C b/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchFields.C index 8100ad1b24..2e7c25d284 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchFields.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/fan/fanFvPatchFields.C @@ -77,7 +77,7 @@ void Foam::fanFvPatchField::calcFanJump() } } - if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + if (phi.dimensions() == dimMass/dimTime) { Un /= patch().lookupPatchField(rhoName_); } diff --git a/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.C index 2105734dea..f339f3a854 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/fanPressure/fanPressureFvPatchScalarField.C @@ -162,11 +162,11 @@ void Foam::fanPressureFvPatchScalarField::updateCoeffs() // Average volumetric flow rate scalar volFlowRate = 0; - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { volFlowRate = dir*gSum(phip); } - else if (phi.dimensions() == dimVelocity*dimArea*dimDensity) + else if (phi.dimensions() == dimMass/dimTime) { const scalarField& rhop = patch().lookupPatchField(rhoName()); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/fluxCorrectedVelocity/fluxCorrectedVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/fluxCorrectedVelocity/fluxCorrectedVelocityFvPatchVectorField.C index f600895920..08f995dbfe 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/fluxCorrectedVelocity/fluxCorrectedVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/fluxCorrectedVelocity/fluxCorrectedVelocityFvPatchVectorField.C @@ -115,11 +115,11 @@ void Foam::fluxCorrectedVelocityFvPatchVectorField::evaluate const vectorField n(patch().nf()); const Field& magS = patch().magSf(); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { operator==(*this - n*(n & *this) + n*phip/magS); } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/mappedFlowRate/mappedFlowRateFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/mappedFlowRate/mappedFlowRateFvPatchVectorField.C index 5bdeefc980..016619b6c4 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/mappedFlowRate/mappedFlowRateFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/mappedFlowRate/mappedFlowRateFvPatchVectorField.C @@ -140,12 +140,12 @@ void Foam::mappedFlowRateFvPatchVectorField::updateCoeffs() vectorField n(patch().nf()); - if (phiName.dimensions() == dimVelocity*dimArea) + if (phiName.dimensions() == dimVolume/dimTime) { // volumetric flow-rate operator==(n*U); } - else if (phiName.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phiName.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/plenumPressure/plenumPressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/plenumPressure/plenumPressureFvPatchScalarField.C index 8497772215..770634339f 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/plenumPressure/plenumPressureFvPatchScalarField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/plenumPressure/plenumPressureFvPatchScalarField.C @@ -196,7 +196,7 @@ void Foam::plenumPressureFvPatchScalarField::updateCoeffs() // Calculate the current mass flow rate scalar massFlowRate(1.0); - if (phi.internalField().dimensions() == dimVelocity*dimArea) + if (phi.internalField().dimensions() == dimVolume/dimTime) { if (hasRho_) { @@ -209,11 +209,7 @@ void Foam::plenumPressureFvPatchScalarField::updateCoeffs() << exit(FatalError); } } - else if - ( - phi.internalField().dimensions() - == dimDensity*dimVelocity*dimArea - ) + else if (phi.internalField().dimensions() == dimMass/dimTime) { if (hasRho_) { diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletOutletVelocity/pressureDirectedInletOutletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletOutletVelocity/pressureDirectedInletOutletVelocityFvPatchVectorField.C index fddcc8a57d..e2c83fb1c1 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletOutletVelocity/pressureDirectedInletOutletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletOutletVelocity/pressureDirectedInletOutletVelocityFvPatchVectorField.C @@ -160,11 +160,11 @@ void Foam::pressureDirectedInletOutletVelocityFvPatchVectorField::updateCoeffs() tmp n = patch().nf(); tmp ndmagS = (n & inletDir_)*patch().magSf(); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { refValue() = inletDir_*phip/ndmagS; } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletVelocity/pressureDirectedInletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletVelocity/pressureDirectedInletVelocityFvPatchVectorField.C index dc7976ecb0..aa8c33a182 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletVelocity/pressureDirectedInletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressureDirectedInletVelocity/pressureDirectedInletVelocityFvPatchVectorField.C @@ -150,11 +150,11 @@ void Foam::pressureDirectedInletVelocityFvPatchVectorField::updateCoeffs() tmp n = patch().nf(); tmp ndmagS = (n & inletDir_)*patch().magSf(); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { operator==(inletDir_*phip/ndmagS); } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressureInletOutletParSlipVelocity/pressureInletOutletParSlipVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressureInletOutletParSlipVelocity/pressureInletOutletParSlipVelocityFvPatchVectorField.C index 4c89daf49b..fc21808f87 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressureInletOutletParSlipVelocity/pressureInletOutletParSlipVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressureInletOutletParSlipVelocity/pressureInletOutletParSlipVelocityFvPatchVectorField.C @@ -134,11 +134,11 @@ void Foam::pressureInletOutletParSlipVelocityFvPatchVectorField::updateCoeffs() vectorField Ut(patchInternalField()); Ut -= n()*(Ut & n()); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { refValue() = Ut + n*phip/magSf; } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressureInletVelocity/pressureInletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressureInletVelocity/pressureInletVelocityFvPatchVectorField.C index e954a2cf79..a7e1fba8fb 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressureInletVelocity/pressureInletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressureInletVelocity/pressureInletVelocityFvPatchVectorField.C @@ -118,11 +118,11 @@ void Foam::pressureInletVelocityFvPatchVectorField::updateCoeffs() tmp n = patch().nf(); const Field& magS = patch().magSf(); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { operator==(n*phip/magS); } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressureNormalInletOutletVelocity/pressureNormalInletOutletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressureNormalInletOutletVelocity/pressureNormalInletOutletVelocityFvPatchVectorField.C index 6340b222d1..c66edef9b1 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressureNormalInletOutletVelocity/pressureNormalInletOutletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressureNormalInletOutletVelocity/pressureNormalInletOutletVelocityFvPatchVectorField.C @@ -130,11 +130,11 @@ void Foam::pressureNormalInletOutletVelocityFvPatchVectorField::updateCoeffs() tmp n = patch().nf(); const Field& magS = patch().magSf(); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { refValue() = n*phip/magS; } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressurePIDControlInletVelocity/pressurePIDControlInletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressurePIDControlInletVelocity/pressurePIDControlInletVelocityFvPatchVectorField.C index 698f3dfc1b..ed7c739dfb 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressurePIDControlInletVelocity/pressurePIDControlInletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressurePIDControlInletVelocity/pressurePIDControlInletVelocityFvPatchVectorField.C @@ -238,11 +238,11 @@ void Foam::pressurePIDControlInletVelocityFvPatchVectorField::updateCoeffs() // Get the density scalar rho = 1; - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { // do nothing ... } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhoField = patch().lookupPatchField(rhoName_); @@ -254,9 +254,9 @@ void Foam::pressurePIDControlInletVelocityFvPatchVectorField::updateCoeffs() FatalErrorInFunction << "The dimensions of the field " << phiName_ << "are not recognised. The dimensions are " << phi.dimensions() - << ". The dimensions should be either " << dimVelocity*dimArea + << ". The dimensions should be either " << dimVolume/dimTime << " for an incompressible case, or " - << dimDensity*dimVelocity*dimArea << " for a compressible case." + << dimMass/dimTime << " for a compressible case." << exit(FatalError); } diff --git a/src/finiteVolume/fields/fvPatchFields/derived/pressurePermeableAlphaInletOutletVelocity/pressurePermeableAlphaInletOutletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/pressurePermeableAlphaInletOutletVelocity/pressurePermeableAlphaInletOutletVelocityFvPatchVectorField.C index 660ae711bb..8f70522e99 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/pressurePermeableAlphaInletOutletVelocity/pressurePermeableAlphaInletOutletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/pressurePermeableAlphaInletOutletVelocity/pressurePermeableAlphaInletOutletVelocityFvPatchVectorField.C @@ -138,11 +138,11 @@ updateCoeffs() const vectorField n(patch().nf()); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { refValue() = (phip/patch().magSf())*n; } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/swirlFlowRateInletVelocity/swirlFlowRateInletVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/swirlFlowRateInletVelocity/swirlFlowRateInletVelocityFvPatchVectorField.C index c0c37d2d9c..7f9092557a 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/swirlFlowRateInletVelocity/swirlFlowRateInletVelocityFvPatchVectorField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/swirlFlowRateInletVelocity/swirlFlowRateInletVelocityFvPatchVectorField.C @@ -171,12 +171,12 @@ void Foam::swirlFlowRateInletVelocityFvPatchVectorField::updateCoeffs() const surfaceScalarField& phi = db().lookupObject(phiName_); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { // volumetric flow-rate operator==(tangentialVelocity + n*avgU); } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/syringePressure/syringePressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/syringePressure/syringePressureFvPatchScalarField.C index 137350c703..82627bce14 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/syringePressure/syringePressureFvPatchScalarField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/syringePressure/syringePressureFvPatchScalarField.C @@ -206,11 +206,11 @@ void Foam::syringePressureFvPatchScalarField::updateCoeffs() const fvsPatchField& phip = patch().patchField(phi); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { ams_ = ams0_ + deltaT*sum((*this*psi_)*phip); } - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { ams_ = ams0_ + deltaT*sum(phip); } diff --git a/src/finiteVolume/fields/fvPatchFields/derived/waveSurfacePressure/waveSurfacePressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/waveSurfacePressure/waveSurfacePressureFvPatchScalarField.C index fca845562e..7e33d7d86b 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/waveSurfacePressure/waveSurfacePressureFvPatchScalarField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/waveSurfacePressure/waveSurfacePressureFvPatchScalarField.C @@ -164,7 +164,7 @@ void Foam::waveSurfacePressureFvPatchScalarField::updateCoeffs() // Change in zeta due to flux vectorField dZetap(dt*nf()*phi.boundaryField()[patchi]/patch().magSf()); - if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + if (phi.dimensions() == dimMass/dimTime) { const scalarField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/finiteVolume/fields/fvPatchFields/derived/waveTransmissive/waveTransmissiveFvPatchField.C b/src/finiteVolume/fields/fvPatchFields/derived/waveTransmissive/waveTransmissiveFvPatchField.C index 530e84ffbd..42d9cd4b21 100644 --- a/src/finiteVolume/fields/fvPatchFields/derived/waveTransmissive/waveTransmissiveFvPatchField.C +++ b/src/finiteVolume/fields/fvPatchFields/derived/waveTransmissive/waveTransmissiveFvPatchField.C @@ -121,7 +121,7 @@ Foam::waveTransmissiveFvPatchField::advectionSpeed() const this->patch().template lookupPatchField(this->phiName_); - if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + if (phi.dimensions() == dimMass/dimTime) { const fvPatchScalarField& rhop = this->patch().template diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/PhiScheme/PhiScheme.C b/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/PhiScheme/PhiScheme.C index 491230f954..da7f02d3eb 100644 --- a/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/PhiScheme/PhiScheme.C +++ b/src/finiteVolume/interpolation/surfaceInterpolation/limitedSchemes/PhiScheme/PhiScheme.C @@ -68,7 +68,7 @@ Foam::PhiScheme::limiter tmp tUflux = this->faceFlux_; - if (this->faceFlux_.dimensions() == dimDensity*dimVelocity*dimArea) + if (this->faceFlux_.dimensions() == dimMass/dimTime) { const volScalarField& rho = phi.db().objectRegistry::template lookupObject @@ -76,7 +76,7 @@ Foam::PhiScheme::limiter tUflux = this->faceFlux_/fvc::interpolate(rho); } - else if (this->faceFlux_.dimensions() != dimVelocity*dimArea) + else if (this->faceFlux_.dimensions() != dimVolume/dimTime) { FatalErrorInFunction << "dimensions of faceFlux are not correct" diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/CoBlended/CoBlended.H b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/CoBlended/CoBlended.H index 6dfe043052..aa1f10e44d 100644 --- a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/CoBlended/CoBlended.H +++ b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/CoBlended/CoBlended.H @@ -199,7 +199,7 @@ public: const fvMesh& mesh = this->mesh(); tmp tUflux = faceFlux_; - if (faceFlux_.dimensions() == dimDensity*dimVelocity*dimArea) + if (faceFlux_.dimensions() == dimMass/dimTime) { // Currently assume that the density field // corresponding to the mass-flux is named "rho" @@ -209,7 +209,7 @@ public: tUflux = faceFlux_/fvc::interpolate(rho); } - else if (faceFlux_.dimensions() != dimVelocity*dimArea) + else if (faceFlux_.dimensions() != dimVolume/dimTime) { FatalErrorInFunction << "dimensions of faceFlux are not correct" diff --git a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/cellCoBlended/cellCoBlended.H b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/cellCoBlended/cellCoBlended.H index 12b277b92f..488713445a 100644 --- a/src/finiteVolume/interpolation/surfaceInterpolation/schemes/cellCoBlended/cellCoBlended.H +++ b/src/finiteVolume/interpolation/surfaceInterpolation/schemes/cellCoBlended/cellCoBlended.H @@ -208,7 +208,7 @@ public: const fvMesh& mesh = this->mesh(); tmp tUflux = faceFlux_; - if (faceFlux_.dimensions() == dimDensity*dimVelocity*dimArea) + if (faceFlux_.dimensions() == dimMass/dimTime) { // Currently assume that the density field // corresponding to the mass-flux is named "rho" @@ -218,7 +218,7 @@ public: tUflux = faceFlux_/fvc::interpolate(rho); } - else if (faceFlux_.dimensions() != dimVelocity*dimArea) + else if (faceFlux_.dimensions() != dimVolume/dimTime) { FatalErrorInFunction << "dimensions of faceFlux are not correct" diff --git a/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C b/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C index 7671078e5e..df48d608ef 100644 --- a/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C +++ b/src/fvOptions/sources/derived/directionalPressureGradientExplicitSource/directionalPressureGradientExplicitSource.C @@ -268,7 +268,7 @@ void Foam::fv::directionalPressureGradientExplicitSource::correct { case pDarcyForchheimer: { - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { const auto& turbModel = mesh().lookupObject @@ -323,7 +323,7 @@ void Foam::fv::directionalPressureGradientExplicitSource::correct } reduce(totalphi, sumOp()); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { volFlowRate = mag(totalphi); } diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C index ae5c3c88f8..b4cf007ec5 100644 --- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C +++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/PatchFlowRateInjection/PatchFlowRateInjection.C @@ -133,7 +133,7 @@ Foam::scalar Foam::PatchFlowRateInjection::flowRate() const const scalarField& phip = phi.boundaryField()[patchId_]; scalar flowRateIn = 0.0; - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) { flowRateIn = max(0.0, -sum(phip)); } diff --git a/src/regionModels/regionCoupling/derivedFvPatchFields/filmPyrolysisVelocityCoupled/filmPyrolysisVelocityCoupledFvPatchVectorField.C b/src/regionModels/regionCoupling/derivedFvPatchFields/filmPyrolysisVelocityCoupled/filmPyrolysisVelocityCoupledFvPatchVectorField.C index 2e3e718946..210f7f942b 100644 --- a/src/regionModels/regionCoupling/derivedFvPatchFields/filmPyrolysisVelocityCoupled/filmPyrolysisVelocityCoupledFvPatchVectorField.C +++ b/src/regionModels/regionCoupling/derivedFvPatchFields/filmPyrolysisVelocityCoupled/filmPyrolysisVelocityCoupledFvPatchVectorField.C @@ -176,9 +176,9 @@ void Foam::filmPyrolysisVelocityCoupledFvPatchVectorField::updateCoeffs() const surfaceScalarField& phi = db().lookupObject(phiName_); - if (phi.dimensions() == dimVelocity*dimArea) + if (phi.dimensions() == dimVolume/dimTime) {} - else if (phi.dimensions() == dimDensity*dimVelocity*dimArea) + else if (phi.dimensions() == dimMass/dimTime) { const fvPatchField& rhop = patch().lookupPatchField(rhoName_); diff --git a/src/thermoTools/derivedFvPatchFields/outletMachNumberPressure/outletMachNumberPressureFvPatchScalarField.C b/src/thermoTools/derivedFvPatchFields/outletMachNumberPressure/outletMachNumberPressureFvPatchScalarField.C index ff96ab6487..5f69d2c0d1 100644 --- a/src/thermoTools/derivedFvPatchFields/outletMachNumberPressure/outletMachNumberPressureFvPatchScalarField.C +++ b/src/thermoTools/derivedFvPatchFields/outletMachNumberPressure/outletMachNumberPressureFvPatchScalarField.C @@ -160,7 +160,7 @@ void Foam::outletMachNumberPressureFvPatchScalarField::updateCoeffs() patch().lookupPatchField(phiName_); // Calculate the current mass flow rate - if (phi.internalField().dimensions() != dimDensity*dimVelocity*dimArea) + if (phi.internalField().dimensions() != dimMass/dimTime) { FatalErrorInFunction <<"phi is not a mass flux." << exit(FatalError);