STYLE: use cbrt() instead of pow(..,1.0/3.0) in a few more places

applications/solvers/finiteArea/liquidFilmFoam/calcFrictionFactor.H

@@ -3,5 +3,5 @@
     // Friction factor is defined with standard gravity
     frictionFactor.primitiveFieldRef() =
         mag(2*9.81*sqr(manningField.primitiveField())/
-        pow(mag(h.primitiveField()) + 1e-7, 1.0/3.0));
+        cbrt(mag(h.primitiveField()) + 1e-7));
 }

applications/solvers/multiphase/compressibleMultiphaseInterFoam/multiphaseMixtureThermo/multiphaseMixtureThermo.C

@@ -107,7 +107,7 @@ Foam::multiphaseMixtureThermo::multiphaseMixtureThermo
     deltaN_
     (
         "deltaN",
-        1e-8/pow(average(mesh_.V()), 1.0/3.0)
+        1e-8/cbrt(average(mesh_.V()))
     )
 {
     rhoPhi_.setOriented();

applications/solvers/multiphase/icoReactingMultiphaseInterFoam/phasesSystem/phaseSystem/phaseSystem.C

@@ -1145,7 +1145,7 @@ Foam::tmp<Foam::surfaceVectorField> Foam::phaseSystem::nHatfv
     const dimensionedScalar deltaN
     (
         "deltaN",
-        1e-8/pow(average(mesh_.V()), 1.0/3.0)
+        1e-8/cbrt(average(mesh_.V()))
     );
 
     // Face unit interface normal

applications/solvers/multiphase/interFoam/interMixingFoam/threePhaseInterfaceProperties/threePhaseInterfaceProperties.C

@@ -166,7 +166,7 @@ Foam::threePhaseInterfaceProperties::threePhaseInterfaceProperties
     deltaN_
     (
         "deltaN",
-        1e-8/pow(average(mixture.U().mesh().V()), 1.0/3.0)
+        1e-8/cbrt(average(mixture.U().mesh().V()))
     ),
 
     nHatf_

applications/solvers/multiphase/multiphaseEulerFoam/UEqns.H

@@ -44,9 +44,9 @@ forAllIter(PtrDictionary<phaseModel>, fluid.phases(), iter)
                *max
                 (
                     mag(U()) - maxSlamVelocity,
-                    dimensionedScalar("U0", dimVelocity, 0)
+                    dimensionedScalar("U0", dimVelocity, Zero)
                 )
-               /pow(mesh.V(), 1.0/3.0),
+               /cbrt(mesh.V()),
                 U
             )
         )

applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/diameterModels/isothermalDiameter/isothermalDiameter.C

@@ -73,7 +73,7 @@ Foam::tmp<Foam::volScalarField> Foam::diameterModels::isothermal::d() const
         "p"
     );
 
-    return d0_*pow(p0_/p, 1.0/3.0);
+    return d0_*cbrt(p0_/p);
 }
 
 

applications/solvers/multiphase/multiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C

@@ -408,7 +408,7 @@ Foam::multiphaseSystem::multiphaseSystem
     deltaN_
     (
         "deltaN",
-        1e-8/pow(average(mesh_.V()), 1.0/3.0)
+        1e-8/cbrt(average(mesh_.V()))
     )
 {
     calcAlphas();

applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C

@@ -120,7 +120,7 @@ Foam::multiphaseMixture::multiphaseMixture
     deltaN_
     (
         "deltaN",
-        1e-8/pow(average(mesh_.V()), 1.0/3.0)
+        1e-8/cbrt(average(mesh_.V()))
     )
 {
     rhoPhi_.setOriented();

applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Tenneti/Tenneti.C

@@ -84,7 +84,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::Tenneti::CdRe() const
 
     volScalarField F0
     (
-        5.81*alpha1/pow3(alpha2) + 0.48*pow(alpha1, 1.0/3.0)/pow4(alpha2)
+        5.81*alpha1/pow3(alpha2) + 0.48*cbrt(alpha1)/pow4(alpha2)
     );
 
     volScalarField F1

applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/isothermalDiameter/isothermalDiameter.C

@@ -73,7 +73,7 @@ Foam::tmp<Foam::volScalarField> Foam::diameterModels::isothermal::d() const
         "p"
     );
 
-    return d0_*pow(p0_/p, 1.0/3.0);
+    return d0_*cbrt(p0_/p);
 }
 
 

applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C

@@ -503,7 +503,7 @@ Foam::multiphaseSystem::multiphaseSystem
     deltaN_
     (
         "deltaN",
-        1e-8/pow(average(mesh_.V()), 1.0/3.0)
+        1e-8/cbrt(average(mesh_.V()))
     )
 {
     forAll(phases(), phasei)

applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/diameterModels/isothermalDiameter/isothermalDiameter.C

@@ -73,7 +73,7 @@ Foam::tmp<Foam::volScalarField> Foam::diameterModels::isothermal::d() const
         "p"
     );
 
-    return d0_*pow(p0_/p, 1.0/3.0);
+    return d0_*cbrt(p0_/p);
 }
 
 

src/TurbulenceModels/turbulenceModels/LES/LESdeltas/cubeRootVolDelta/cubeRootVolDelta.C

@@ -48,7 +48,7 @@ void Foam::LESModels::cubeRootVolDelta::calcDelta()
 
     if (nD == 3)
     {
-        delta_.primitiveFieldRef() = deltaCoeff_*pow(mesh.V(), 1.0/3.0);
+        delta_.primitiveFieldRef() = deltaCoeff_ * cbrt(mesh.V());
     }
     else if (nD == 2)
     {

src/combustionModels/eddyDissipationDiffusionModel/eddyDissipationDiffusionModel.C

@@ -94,7 +94,7 @@ eddyDissipationDiffusionModel<ReactionThermo, ThermoType>::rtDiff() const
     );
 
     volScalarField& delta = tdelta.ref();
-    delta.ref() = pow(this->mesh().V(), 1.0/3.0);
+    delta.ref() = cbrt(this->mesh().V());
     delta.correctBoundaryConditions();
 
     // NOTE: Assume Prt = 1

src/dynamicMesh/polyTopoChange/polyTopoChange/edgeCollapser.C

@@ -470,8 +470,8 @@ Foam::scalarField Foam::edgeCollapser::calcTargetFaceSizes() const
         const scalar cellOwnerVol = max(0.0, V[cellOwner[intFacei]]);
         const scalar cellNeighbourVol = max(0.0, V[cellNeighbour[intFacei]]);
 
-        scalar targetFaceSizeA = Foam::pow(cellOwnerVol, 1.0/3.0);
+        scalar targetFaceSizeA = Foam::cbrt(cellOwnerVol);
-        scalar targetFaceSizeB = Foam::pow(cellNeighbourVol, 1.0/3.0);
+        scalar targetFaceSizeB = Foam::cbrt(cellNeighbourVol);
 
         targetFaceSizes[intFacei] = 0.5*(targetFaceSizeA + targetFaceSizeB);
     }
@@ -505,7 +505,7 @@ Foam::scalarField Foam::edgeCollapser::calcTargetFaceSizes() const
                 const label extFacei = patchFacei + patch.start();
                 const scalar cellOwnerVol = max(0.0, V[cellOwner[extFacei]]);
 
-                targetFaceSizes[extFacei] = Foam::pow(cellOwnerVol, 1.0/3.0);
+                targetFaceSizes[extFacei] = Foam::cbrt(cellOwnerVol);
             }
         }
     }
@@ -526,8 +526,8 @@ Foam::scalarField Foam::edgeCollapser::calcTargetFaceSizes() const
                 const scalar cellOwnerVol = max(0.0, V[cellOwner[localFacei]]);
                 const scalar cellNeighbourVol = neiCellVolumes[bFacei++];
 
-                scalar targetFaceSizeA = Foam::pow(cellOwnerVol, 1.0/3.0);
+                scalar targetFaceSizeA = Foam::cbrt(cellOwnerVol);
-                scalar targetFaceSizeB = Foam::pow(cellNeighbourVol, 1.0/3.0);
+                scalar targetFaceSizeB = Foam::cbrt(cellNeighbourVol);
 
                 targetFaceSizes[localFacei]
                     = 0.5*(targetFaceSizeA + targetFaceSizeB);

src/functionObjects/field/stabilityBlendingFactor/stabilityBlendingFactor.C

@@ -370,7 +370,7 @@ bool Foam::functionObjects::stabilityBlendingFactor::init(bool first)
         auto& Co = CoPtr.ref();
 
         Co.primitiveFieldRef() =
-            mesh_.time().deltaT()*mag(*UNamePtr)/pow(mesh_.V(), 1.0/3.0);
+            mesh_.time().deltaT()*mag(*UNamePtr)/cbrt(mesh_.V());
 
         indicator_ =
             max

src/lagrangian/intermediate/submodels/MPPIC/AveragingMethods/Moment/Moment.C

@@ -42,7 +42,7 @@ Foam::AveragingMethods::Moment<Type>::Moment
     dataY_(FieldField<Field, Type>::operator[](2)),
     dataZ_(FieldField<Field, Type>::operator[](3)),
     transform_(mesh.nCells(), Zero),
-    scale_(0.5*pow(mesh.V(), 1.0/3.0))
+    scale_(0.5*cbrt(mesh.V()))
 {
     scalar a = 1.0/24.0;
     scalar b = 0.5854101966249685;

src/lagrangian/intermediate/submodels/MPPIC/ParticleStressModels/Lun/Lun.C

@@ -89,7 +89,7 @@ Foam::ParticleStressModels::Lun::tau
         0.6
       / max
         (
-            1.0 - pow(alpha/alphaPacked_,1.0/3.0),
+            1.0 - cbrt(alpha/alphaPacked_),
             max(eps_*(1.0 - alpha), SMALL)
         )
     );

src/regionModels/surfaceFilmModels/derivedFvPatchFields/inclinedFilmNusseltHeight/inclinedFilmNusseltHeightFvPatchScalarField.C

@@ -179,7 +179,7 @@ void Foam::inclinedFilmNusseltHeightFvPatchScalarField::updateCoeffs()
 
     operator==
     (
-        pow(3*sqr(mup/rhop)/(gTan + ROOTVSMALL), 1.0/3.0)*pow(Re, 1.0/3.0)
+        cbrt(3*sqr(mup/rhop)/(gTan + ROOTVSMALL))*cbrt(Re)
     );
 
     fixedValueFvPatchScalarField::updateCoeffs();

src/regionModels/surfaceFilmModels/derivedFvPatchFields/inclinedFilmNusseltInletVelocity/inclinedFilmNusseltInletVelocityFvPatchVectorField.C

@@ -179,7 +179,7 @@ void Foam::inclinedFilmNusseltInletVelocityFvPatchVectorField::updateCoeffs()
 
     const scalarField Re(max(G, scalar(0))/mup);
 
-    operator==(n*pow(gTan*mup/(3*rhop), 1.0/3.0)*pow(Re, 2.0/3.0));
+    operator==(n*cbrt(gTan*mup/(3*rhop))*pow(Re, 2.0/3.0));
 
     fixedValueFvPatchVectorField::updateCoeffs();
 }

src/thermophysicalModels/specie/equationOfState/PengRobinsonGas/PengRobinsonGasI.H

@@ -236,15 +236,15 @@ inline Foam::scalar Foam::PengRobinsonGas<Specie>::Z
     {
         // One root is real
         const scalar D05 = sqrt(D);
-        const scalar S = pow(Rl + D05, 1.0/3.0);
+        const scalar S = cbrt(Rl + D05);
         scalar Tl = 0;
         if (D05 > Rl)
         {
-            Tl = -pow(mag(Rl - D05), 1.0/3.0);
+            Tl = -cbrt(mag(Rl - D05));
         }
         else
         {
-            Tl = pow(Rl - D05, 1.0/3.0);
+            Tl = cbrt(Rl - D05);
         }
 
         root = S + Tl - a2/3.0;

src/transportModels/interfaceProperties/interfaceProperties.C

@@ -165,7 +165,7 @@ Foam::interfaceProperties::interfaceProperties
     deltaN_
     (
         "deltaN",
-        1e-8/pow(average(alpha1.mesh().V()), 1.0/3.0)
+        1e-8/cbrt(average(alpha1.mesh().V()))
     ),
 
     alpha1_(alpha1),