diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C index 2bc5d377de989db355cb865ab3ee430762419818..3fff5e49150dbaab8632396c1d38cdfa2966c75b 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C @@ -63,7 +63,7 @@ Foam::dragModels::Ergun::~Ergun() Foam::tmp<Foam::volScalarField> Foam::dragModels::Ergun::CdRe() const { return - (4/3) + (4.0/3.0) *( 150 *max diff --git a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C index ca378985ea568e1cf05c9c5b39d701e2ae52a2c5..cb318e5d075ecb0924f2a32d4e533b2fe3240b05 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C +++ b/applications/solvers/multiphase/reactingEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C @@ -68,7 +68,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::Gibilaro::CdRe() const ); return - (4/3) + (4.0/3.0) *(17.3/alpha2 + 0.336*pair_.Re()) *max(pair_.continuous(), pair_.continuous().residualAlpha()) *pow(alpha2, -2.8); diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C index 2bc5d377de989db355cb865ab3ee430762419818..3fff5e49150dbaab8632396c1d38cdfa2966c75b 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Ergun/Ergun.C @@ -63,7 +63,7 @@ Foam::dragModels::Ergun::~Ergun() Foam::tmp<Foam::volScalarField> Foam::dragModels::Ergun::CdRe() const { return - (4/3) + (4.0/3.0) *( 150 *max diff --git a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C index ca378985ea568e1cf05c9c5b39d701e2ae52a2c5..cb318e5d075ecb0924f2a32d4e533b2fe3240b05 100644 --- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C +++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/Gibilaro/Gibilaro.C @@ -68,7 +68,7 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::Gibilaro::CdRe() const ); return - (4/3) + (4.0/3.0) *(17.3/alpha2 + 0.336*pair_.Re()) *max(pair_.continuous(), pair_.continuous().residualAlpha()) *pow(alpha2, -2.8); diff --git a/applications/utilities/preProcessing/wallFunctionTable/tabulatedWallFunction/SpaldingsLaw/SpaldingsLaw.C b/applications/utilities/preProcessing/wallFunctionTable/tabulatedWallFunction/SpaldingsLaw/SpaldingsLaw.C index 24834c227f20e1aebe55d5aca325b8855ba47f2d..27b20c59324ed9e81570b24fd30de4794632ba7c 100644 --- a/applications/utilities/preProcessing/wallFunctionTable/tabulatedWallFunction/SpaldingsLaw/SpaldingsLaw.C +++ b/applications/utilities/preProcessing/wallFunctionTable/tabulatedWallFunction/SpaldingsLaw/SpaldingsLaw.C @@ -91,7 +91,7 @@ void Foam::tabulatedWallFunctions::SpaldingsLaw::invertFunction() ( 2*E_*uPlus + exp(kUPlus)*(kUPlus + 1) - - 2/3*pow3(kUPlus) + - 2.0/3.0*pow3(kUPlus) - 1.5*sqr(kUPlus) - 2*kUPlus - 1 diff --git a/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C b/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C index 740856a5a8d15cb177c417cfcfa2180df5dc8317..6c1cfd9bab69e98fd4084d211ef1ef5f32bbf9ad 100644 --- a/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C +++ b/src/functionObjects/utilities/thermoCoupleProbes/thermoCoupleProbes.C @@ -128,13 +128,13 @@ void Foam::functionObjects::thermoCoupleProbes::derivatives scalarField Re(rhoc*Uc*d_/muc); scalarField Pr(Cpc*muc/kappac); Pr = max(ROOTVSMALL, Pr); - scalarField Nu(2.0 + (0.4*sqrt(Re) + 0.06*pow(Re, 2/3))*pow(Pr, 0.4)); + scalarField Nu(2.0 + (0.4*sqrt(Re) + 0.06*pow(Re, 2.0/3.0))*pow(Pr, 0.4)); scalarField htc(Nu*kappac/d_); const scalar sigma = physicoChemical::sigma.value(); scalar area = 4*constant::mathematical::pi*sqr(0.5*d_); - scalar volume = (4/3)*constant::mathematical::pi*pow3(0.5*d_); + scalar volume = (4.0/3.0)*constant::mathematical::pi*pow3(0.5*d_); dydx = (epsilon_*(G/4 - sigma*pow4(y))*area + htc*(Tc - y)*area) diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C index 6f9ccfe51ee29ef8cc75e5f05b0ac094b8f40e4f..08237cb39879a4e0dfdd70fbb8258ec3cec345e8 100644 --- a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C +++ b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C @@ -315,7 +315,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::wetSplashInteraction const scalar epsilon = 0.993 - theta*(1.76 - theta*(1.56 - theta*0.49)); // Update parcel velocity - U = -epsilon*(Un) + 5/7*(Ut); + U = -epsilon*(Un) + 5.0/7.0*(Ut); keepParticle = true; return;