diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C index 31a0a70328..fb547ae901 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C +++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleSlip/JohnsonJacksonParticleSlipFvPatchVectorField.C @@ -192,21 +192,6 @@ void Foam::JohnsonJacksonParticleSlipFvPatchVectorField::updateCoeffs() : alpha ); - // lookup the packed volume fraction - dimensionedScalar alphaMax - ( - "alphaMax", - dimless, - db() - .lookupObject - ( - IOobject::groupName("momentumTransport", phase.name()) - ) - .subDict("RAS") - .subDict("kineticTheoryCoeffs") - .lookup("alphaMax") - ); - // calculate the slip value fraction scalarField c ( @@ -215,7 +200,7 @@ void Foam::JohnsonJacksonParticleSlipFvPatchVectorField::updateCoeffs() *gs0 *specularityCoefficient_.value() *sqrt(3*Theta) - /max(6*(nu - nuFric)*alphaMax.value(), small) + /max(6*(nu - nuFric)*phase.alphaMax(), small) ); this->valueFraction() = c/(c + patch().deltaCoeffs()); diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C index ec376283a7..c9725d7a0b 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C +++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/derivedFvPatchFields/JohnsonJacksonParticleTheta/JohnsonJacksonParticleThetaFvPatchScalarField.C @@ -205,21 +205,6 @@ void Foam::JohnsonJacksonParticleThetaFvPatchScalarField::updateCoeffs() const scalarField Theta(patchInternalField()); - // lookup the packed volume fraction - dimensionedScalar alphaMax - ( - "alphaMax", - dimless, - db() - .lookupObject - ( - IOobject::groupName("momentumTransport", phase.name()) - ) - .subDict("RAS") - .subDict("kineticTheoryCoeffs") - .lookup("alphaMax") - ); - // calculate the reference value and the value fraction if (restitutionCoefficient_.value() != 1.0) { @@ -238,7 +223,7 @@ void Foam::JohnsonJacksonParticleThetaFvPatchScalarField::updateCoeffs() *gs0 *(scalar(1) - sqr(restitutionCoefficient_.value())) *sqrt(3*Theta) - /max(4*kappa*alphaMax.value(), small) + /max(4*kappa*phase.alphaMax(), small) ); this->valueFraction() = c/(c + patch().deltaCoeffs()); @@ -258,7 +243,7 @@ void Foam::JohnsonJacksonParticleThetaFvPatchScalarField::updateCoeffs() *gs0 *sqrt(3*Theta) *magSqr(U) - /max(6*kappa*alphaMax.value(), small); + /max(6*kappa*phase.alphaMax(), small); this->valueFraction() = 0; } diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C index 929f02b298..f0bbca4a24 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C +++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C @@ -130,7 +130,6 @@ Foam::RASModels::kineticTheoryModel::kineticTheoryModel equilibrium_(coeffDict_.lookup("equilibrium")), e_("e", dimless, coeffDict_), - alphaMax_("alphaMax", dimless, coeffDict_), alphaMinFriction_ ( "alphaMinFriction", @@ -245,7 +244,6 @@ bool Foam::RASModels::kineticTheoryModel::read() { coeffDict().lookup("equilibrium") >> equilibrium_; e_.readIfPresent(coeffDict()); - alphaMax_.readIfPresent(coeffDict()); alphaMinFriction_.readIfPresent(coeffDict()); viscosityModel_->read(); @@ -316,8 +314,18 @@ Foam::RASModels::kineticTheoryModel::pPrime() const *granularPressureModel_->granularPressureCoeffPrime ( alpha_, - radialModel_->g0(alpha_, alphaMinFriction_, alphaMax_), - radialModel_->g0prime(alpha_, alphaMinFriction_, alphaMax_), + radialModel_->g0 + ( + alpha_, + alphaMinFriction_, + phase_.alphaMax() + ), + radialModel_->g0prime + ( + alpha_, + alphaMinFriction_, + phase_.alphaMax() + ), rho, e_ ) @@ -325,7 +333,7 @@ Foam::RASModels::kineticTheoryModel::pPrime() const ( phase_, alphaMinFriction_, - alphaMax_ + phase_.alphaMax() ) ) ); @@ -413,7 +421,7 @@ void Foam::RASModels::kineticTheoryModel::correct() const volSymmTensorField D(symm(gradU)); // Calculating the radial distribution function - gs0_ = radialModel_->g0(alpha, alphaMinFriction_, alphaMax_); + gs0_ = radialModel_->g0(alpha, alphaMinFriction_, phase_.alphaMax()); if (!equilibrium_) { @@ -587,7 +595,7 @@ void Foam::RASModels::kineticTheoryModel::correct() ( phase_, alphaMinFriction_, - alphaMax_ + phase_.alphaMax() ) ); @@ -600,7 +608,7 @@ void Foam::RASModels::kineticTheoryModel::correct() ( phase_, alphaMinFriction_, - alphaMax_, + phase_.alphaMax(), pf/rho, D ), diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.H b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.H index ce7cc07822..d6932efe51 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.H +++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.H @@ -56,7 +56,6 @@ SourceFiles #include "granularPressureModel.H" #include "frictionalStressModel.H" - // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // namespace Foam @@ -74,12 +73,10 @@ class kineticTheoryModel { // Private Data - // Input Fields + const phaseModel& phase_; - const phaseModel& phase_; - - //- Name of the continuous phase - word continuousPhaseName_; + //- Name of the continuous phase + word continuousPhaseName_; // Sub-models @@ -110,9 +107,6 @@ class kineticTheoryModel //- Coefficient of restitution dimensionedScalar e_; - //- Maximum packing phase-fraction - dimensionedScalar alphaMax_; - //- Min value for which the frictional stresses are zero dimensionedScalar alphaMinFriction_; diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.C b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.C index 5123a11639..4a76e0a2dc 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.C +++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.C @@ -49,7 +49,8 @@ Foam::RASModels::phasePressureModel::phasePressureModel viscosity ), - alphaMax_(coeffDict_.lookup("alphaMax")), + phase_(refCast(viscosity)), + preAlphaExp_(coeffDict_.lookup("preAlphaExp")), expMax_(coeffDict_.lookup("expMax")), g0_ @@ -84,7 +85,6 @@ bool Foam::RASModels::phasePressureModel::read() read() ) { - coeffDict().lookup("alphaMax") >> alphaMax_; coeffDict().lookup("preAlphaExp") >> preAlphaExp_; coeffDict().lookup("expMax") >> expMax_; g0_.readIfPresent(coeffDict()); @@ -153,7 +153,7 @@ Foam::RASModels::phasePressureModel::pPrime() const g0_ *min ( - exp(preAlphaExp_*(alpha_ - alphaMax_)), + exp(preAlphaExp_*(alpha_ - phase_.alphaMax())), expMax_ ) ) @@ -185,7 +185,8 @@ Foam::RASModels::phasePressureModel::pPrimef() const g0_ *min ( - exp(preAlphaExp_*(fvc::interpolate(alpha_) - alphaMax_)), + exp(preAlphaExp_ + *(fvc::interpolate(alpha_) - phase_.alphaMax())), expMax_ ) ) diff --git a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.H b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.H index a5c1909659..b773deb104 100644 --- a/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.H +++ b/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseCompressibleMomentumTransportModels/phasePressureModel/phasePressureModel.H @@ -38,7 +38,6 @@ Description { preAlphaExp 500; expMax 1000; - alphaMax 0.62; g0 1000; } \endverbatim @@ -73,10 +72,9 @@ class phasePressureModel { // Private Data - // Kinetic Theory Model coefficients + const phaseModel& phase_; - //- Maximum packing phase-fraction - scalar alphaMax_; + // Phase pressure coefficients //- Pre-exponential factor scalar preAlphaExp_; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/momentumTransport.solids b/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/momentumTransport.solids index 1b9afa8560..50805535ac 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/momentumTransport.solids +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/momentumTransport.solids @@ -28,7 +28,6 @@ RAS equilibrium off; e 0.8; - alphaMax 0.65; alphaMinFriction 0.5; residualAlpha 1e-6; diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/phaseProperties b/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/phaseProperties index 992534d148..529d07c538 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/phaseProperties +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/LBend/constant/phaseProperties @@ -29,7 +29,8 @@ solids d 462e-6; } - residualAlpha 1e-5; + alphaMax 0.65; + residualAlpha 1e-5; } gas diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/momentumTransport.particles b/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/momentumTransport.particles index 92906d8fee..bd79d070a1 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/momentumTransport.particles +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/momentumTransport.particles @@ -28,7 +28,6 @@ RAS equilibrium off; e 0.8; - alphaMax 0.62; alphaMinFriction 0.5; residualAlpha 1e-4; @@ -52,7 +51,6 @@ RAS { preAlphaExp 500; expMax 1000; - alphaMax 0.62; g0 1000; } } diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties b/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties index defc671d89..ac958d710b 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties @@ -53,13 +53,7 @@ blending } surfaceTension -{ - air_particles - { - type constant; - sigma 0; - } -} +{} interfaceCompression {} diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/momentumTransport.particles b/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/momentumTransport.particles index 93f019571f..40fe68bb14 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/momentumTransport.particles +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/momentumTransport.particles @@ -26,7 +26,6 @@ RAS { preAlphaExp 500; expMax 1000; - alphaMax 0.62; g0 1000; } } diff --git a/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/phaseProperties b/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/phaseProperties index 0e6742fe91..785868243a 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/phaseProperties +++ b/tutorials/multiphase/multiphaseEulerFoam/RAS/pipeBend/constant/phaseProperties @@ -66,6 +66,7 @@ particles ); } + alphaMax 0.62; residualAlpha 1e-8; } @@ -176,4 +177,4 @@ turbulentDispersion } } -// ************************************************************************* // \ No newline at end of file +// ************************************************************************* // diff --git a/tutorials/multiphase/multiphaseEulerFoam/laminar/fluidisedBed/constant/momentumTransport.particles b/tutorials/multiphase/multiphaseEulerFoam/laminar/fluidisedBed/constant/momentumTransport.particles index 4556b6aef5..a779bc9176 100644 --- a/tutorials/multiphase/multiphaseEulerFoam/laminar/fluidisedBed/constant/momentumTransport.particles +++ b/tutorials/multiphase/multiphaseEulerFoam/laminar/fluidisedBed/constant/momentumTransport.particles @@ -28,7 +28,6 @@ RAS equilibrium on; e 0.8; - alphaMax 0.62; alphaMinFriction 0.5; granularViscosityModel Gidaspow; @@ -51,7 +50,6 @@ RAS { preAlphaExp 500; expMax 1000; - alphaMax 0.62; g0 1000; } }