zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

TurbulenceModels/phaseCompressible: Up-cast the phase-system to twoPhaseSystem

Browse Source
This commit is contained in:
Henry
2015-06-07 18:53:22 +01:00
parent 57d75a691f
commit 661f5aa50c
6 changed files with 22 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
src/TurbulenceModels/phaseCompressible/LES/Niceno/NicenoKEqn.C
View File

@ -134,7 +134,8 @@ NicenoKEqn<BasicTurbulenceModel>::gasTurbulence() const
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
const transportModel& liquid = this->transport(); const transportModel& liquid = this->transport();
const twoPhaseSystem& fluid = liquid.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(liquid.fluid());
const transportModel& gas = fluid.otherPhase(liquid); const transportModel& gas = fluid.otherPhase(liquid);
gasTurbulencePtr_ = gasTurbulencePtr_ =
@ -175,7 +176,8 @@ tmp<volScalarField> NicenoKEqn<BasicTurbulenceModel>::bubbleG() const
this->gasTurbulence(); this->gasTurbulence();
const transportModel& liquid = this->transport(); const transportModel& liquid = this->transport();
const twoPhaseSystem& fluid = liquid.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(liquid.fluid());
const transportModel& gas = fluid.otherPhase(liquid); const transportModel& gas = fluid.otherPhase(liquid);
volScalarField magUr(mag(this->U_ - gasTurbulence.U())); volScalarField magUr(mag(this->U_ - gasTurbulence.U()));

3
src/TurbulenceModels/phaseCompressible/LES/SmagorinskyZhang/SmagorinskyZhang.C
View File

@ -109,7 +109,8 @@ SmagorinskyZhang<BasicTurbulenceModel>::gasTurbulence() const
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
const transportModel& liquid = this->transport(); const transportModel& liquid = this->transport();
const twoPhaseSystem& fluid = liquid.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(liquid.fluid());
const transportModel& gas = fluid.otherPhase(liquid); const transportModel& gas = fluid.otherPhase(liquid);
gasTurbulencePtr_ = gasTurbulencePtr_ =

3
src/TurbulenceModels/phaseCompressible/LES/continuousGasKEqn/continuousGasKEqn.C
View File

@ -108,7 +108,8 @@ continuousGasKEqn<BasicTurbulenceModel>::liquidTurbulence() const
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
liquidTurbulencePtr_ = liquidTurbulencePtr_ =

5
src/TurbulenceModels/phaseCompressible/RAS/LaheyKEpsilon/LaheyKEpsilon.C
View File

@ -145,7 +145,8 @@ LaheyKEpsilon<BasicTurbulenceModel>::gasTurbulence() const
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
const transportModel& liquid = this->transport(); const transportModel& liquid = this->transport();
const twoPhaseSystem& fluid = liquid.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(liquid.fluid());
const transportModel& gas = fluid.otherPhase(liquid); const transportModel& gas = fluid.otherPhase(liquid);
gasTurbulencePtr_ = gasTurbulencePtr_ =
@ -186,7 +187,7 @@ tmp<volScalarField> LaheyKEpsilon<BasicTurbulenceModel>::bubbleG() const
this->gasTurbulence(); this->gasTurbulence();
const transportModel& liquid = this->transport(); const transportModel& liquid = this->transport();
const twoPhaseSystem& fluid = liquid.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(liquid.fluid());
const transportModel& gas = fluid.otherPhase(liquid); const transportModel& gas = fluid.otherPhase(liquid);
volScalarField magUr(mag(this->U_ - gasTurbulence.U())); volScalarField magUr(mag(this->U_ - gasTurbulence.U()));

7
src/TurbulenceModels/phaseCompressible/RAS/continuousGasKEpsilon/continuousGasKEpsilon.C
View File

@ -120,7 +120,7 @@ void continuousGasKEpsilon<BasicTurbulenceModel>::correctNut()
const turbulenceModel& liquidTurbulence = this->liquidTurbulence(); const turbulenceModel& liquidTurbulence = this->liquidTurbulence();
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
volScalarField thetal(liquidTurbulence.k()/liquidTurbulence.epsilon()); volScalarField thetal(liquidTurbulence.k()/liquidTurbulence.epsilon());
@ -149,7 +149,8 @@ continuousGasKEpsilon<BasicTurbulenceModel>::liquidTurbulence() const
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
liquidTurbulencePtr_ = liquidTurbulencePtr_ =
@ -202,7 +203,7 @@ tmp<Foam::volScalarField>
continuousGasKEpsilon<BasicTurbulenceModel>::rhoEff() const continuousGasKEpsilon<BasicTurbulenceModel>::rhoEff() const
{ {
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
return tmp<volScalarField> return tmp<volScalarField>

13
src/TurbulenceModels/phaseCompressible/RAS/mixtureKEpsilon/mixtureKEpsilon.C
View File

@ -347,7 +347,8 @@ mixtureKEpsilon<BasicTurbulenceModel>::liquidTurbulence() const
const volVectorField& U = this->U_; const volVectorField& U = this->U_;
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid =
refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
liquidTurbulencePtr_ = liquidTurbulencePtr_ =
@ -375,7 +376,7 @@ tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::Ct2() const
this->liquidTurbulence(); this->liquidTurbulence();
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
const volScalarField& alphag = this->alpha_; const volScalarField& alphag = this->alpha_;
@ -399,7 +400,7 @@ template<class BasicTurbulenceModel>
tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::rholEff() const tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::rholEff() const
{ {
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
return fluid.otherPhase(gas).rho(); return fluid.otherPhase(gas).rho();
} }
@ -408,7 +409,7 @@ template<class BasicTurbulenceModel>
tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::rhogEff() const tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::rhogEff() const
{ {
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
return return
gas.rho() gas.rho()
+ fluid.virtualMass(gas).Cvm()*fluid.otherPhase(gas).rho(); + fluid.virtualMass(gas).Cvm()*fluid.otherPhase(gas).rho();
@ -484,7 +485,7 @@ tmp<volScalarField> mixtureKEpsilon<BasicTurbulenceModel>::bubbleG() const
this->liquidTurbulence(); this->liquidTurbulence();
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
const transportModel& liquid = fluid.otherPhase(gas); const transportModel& liquid = fluid.otherPhase(gas);
volScalarField magUr(mag(liquidTurbulence.U() - this->U())); volScalarField magUr(mag(liquidTurbulence.U() - this->U()));
@ -531,7 +532,7 @@ template<class BasicTurbulenceModel>
void mixtureKEpsilon<BasicTurbulenceModel>::correct() void mixtureKEpsilon<BasicTurbulenceModel>::correct()
{ {
const transportModel& gas = this->transport(); const transportModel& gas = this->transport();
const twoPhaseSystem& fluid = gas.fluid(); const twoPhaseSystem& fluid = refCast<const twoPhaseSystem>(gas.fluid());
// Only solve the mixture turbulence for the gas-phase // Only solve the mixture turbulence for the gas-phase
if (&gas != &fluid.phase1()) if (&gas != &fluid.phase1())