BUG: twoPhaseEulerFoam: removed multiple of the continuous phase fraction from force and heat transfer models

2025-11-28 03:28:01 +00:00 · 2014-03-17 16:34:32 +00:00
parent 354197bf33
commit 37c75a8b18
20 changed files with 34 additions and 48 deletions
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/UEqns.H
@ -23,12 +23,11 @@ volScalarField dragCoeff(fluid.dragCoeff());
          + phase1.turbulence().divDevReff(U1)
         ==
          - fvm::Sp(dragCoeff/rho1, U1)
-          - alpha1*alpha2/rho1
-           *(
+          - (
                liftForce
              + wallLubricationForce
              + turbulentDispersionForce
-            )
+            )/rho1
          - virtualMassCoeff/rho1
           *(
                fvm::ddt(U1)
@ -50,12 +49,11 @@ volScalarField dragCoeff(fluid.dragCoeff());
          + phase2.turbulence().divDevReff(U2)
         ==
          - fvm::Sp(dragCoeff/rho2, U2)
-          + alpha1*alpha2/rho2
-           *(
+          + (
                liftForce
              + wallLubricationForce
              + turbulentDispersionForce
-            )
+            )/rho2
          - virtualMassCoeff/rho2
           *(
                fvm::ddt(U2)
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/dragModel/dragModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/dragModels/dragModel/dragModel.C
@ -107,13 +107,11 @@ Foam::tmp<Foam::volScalarField> Foam::dragModel::K() const
    return
        0.75
       *CdRe()
+       *max(pair_.dispersed(), residualAlpha_)
       *swarmCorrection_->Cs()
       *pair_.continuous().rho()
       *pair_.continuous().nu()
-       /(
-            max(pair_.continuous(), residualAlpha_)
-           *sqr(pair_.dispersed().d())
-        );
+       /sqr(pair_.dispersed().d());
 }


--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C
@ -64,7 +64,12 @@ Foam::heatTransferModels::RanzMarshall::K() const
 {
    volScalarField Nu(scalar(2) + 0.6*pair_.Re()*cbrt(pair_.Pr()));

-    return 6.0*pair_.continuous().kappa()*Nu/sqr(pair_.dispersed().d());
+    return
+        6.0
+       *max(pair_.dispersed(), residualAlpha_)
+       *pair_.continuous().kappa()
+       *Nu
+       /sqr(pair_.dispersed().d());
 }


--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/heatTransferModels/heatTransferModel/heatTransferModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/heatTransferModels/heatTransferModel/heatTransferModel.C
@ -45,7 +45,8 @@ Foam::heatTransferModel::heatTransferModel
    const phasePair& pair
 )
 :
-    pair_(pair)
+    pair_(pair),
+    residualAlpha_("residualAlpha", dimless, dict.lookup("residualAlpha"))
 {}


--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/heatTransferModels/heatTransferModel/heatTransferModel.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/heatTransferModels/heatTransferModel/heatTransferModel.H
@ -58,6 +58,9 @@ protected:
        //- Phase pair
        const phasePair& pair_;

+        //- Residual phase fraction
+        const dimensionedScalar residualAlpha_;
+

 public:

--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/liftModel/liftModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/liftModels/liftModel/liftModel.C
@ -62,6 +62,7 @@ Foam::tmp<Foam::volVectorField> Foam::liftModel::F() const
 {
    return
        Cl()
+       *pair_.dispersed()
       *pair_.continuous().rho()
       *(
            pair_.Ur() ^ fvc::curl(pair_.continuous().U())
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/Gosman/Gosman.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/Gosman/Gosman.C
@ -85,6 +85,7 @@ Foam::turbulentDispersionModels::Gosman::F() const
    return
      - 0.75
       *drag.CdRe()
+       *pair_.dispersed()
       *pair_.continuous().nu()
       *pair_.continuous().turbulence().nut()
       /(
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/constantTurbulentDispersionCoefficient/constantTurbulentDispersionCoefficient.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/turbulentDispersionModels/constantTurbulentDispersionCoefficient/constantTurbulentDispersionCoefficient.C
@ -75,6 +75,7 @@ F() const
 {
    return 
      - Ctd_
+       *pair_.dispersed()
       *pair_.continuous().rho()
       *pair_.continuous().turbulence().k()
       *fvc::grad(pair_.dispersed());
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/virtualMassModels/virtualMassModel/virtualMassModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/virtualMassModels/virtualMassModel/virtualMassModel.C
@ -72,7 +72,7 @@ Foam::virtualMassModel::~virtualMassModel()

 Foam::tmp<Foam::volScalarField> Foam::virtualMassModel::K() const
 {
-    return Cvm()*pair_.continuous().rho();
+    return Cvm()*pair_.dispersed()*pair_.continuous().rho();
 }


--- a/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/wallLubricationModels/Antal/Antal.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/interfacialModels/wallLubricationModels/Antal/Antal.C
@ -79,6 +79,7 @@ Foam::tmp<Foam::volVectorField> Foam::wallLubricationModels::Antal::F() const
            dimensionedScalar("zero", dimless/dimLength, 0),
            Cw1_/pair_.dispersed().d() + Cw2_/yWall_
        )
+       *pair_.dispersed()
       *pair_.continuous().rho()
       *magSqr(Ur - (Ur & nWall)*nWall)
       *nWall;
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/BlendedInterfacialModel.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/BlendedInterfacialModel.C
@ -163,8 +163,6 @@ Foam::BlendedInterfacialModel<modelType>::K() const

    if (model_.valid() || model1In2_.valid() || model2In1_.valid())
    {
-        x() *= max(pair_.phase1()*pair_.phase2(), blending_.residualAlpha());
-
        correctFixedFluxBCs(x());
    }

--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.C
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.C
@ -40,8 +40,6 @@ Foam::blendingMethod::blendingMethod
 (
    const dictionary& dict
 )
-:
-    residualAlpha_("residualAlpha", dimless, dict.lookup("residualAlpha"))
 {}


@ -51,12 +49,4 @@ Foam::blendingMethod::~blendingMethod()
 {}


-// * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
-
-const Foam::dimensionedScalar& Foam::blendingMethod::residualAlpha() const
-{
-    return residualAlpha_;
-}
-
-
 // ************************************************************************* //
--- a/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.H
+++ b/applications/solvers/multiphase/twoPhaseEulerFoam/twoPhaseSystem/BlendedInterfacialModel/blendingMethods/blendingMethod/blendingMethod.H
@ -49,14 +49,6 @@ namespace Foam

 class blendingMethod
 {
-protected:
-
-    // Protected data
-
-        //- Residual phase fraction
-        const dimensionedScalar residualAlpha_;
-
-
 public:

    //- Runtime type information
@ -101,9 +93,6 @@ public:

    // Member Functions

-        //- Residual phase fraction
-        const dimensionedScalar& residualAlpha() const;
-
        //- Factor for first phase
        virtual tmp<volScalarField> f1
        (
--- a/tutorials/multiphase/twoPhaseEulerFoam/LES/bubbleColumn/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/LES/bubbleColumn/constant/phaseProperties
@ -41,7 +41,6 @@ blending
    default
    {
        type            linear;
-        residualAlpha   1e-6;
        maxFullyDispersedAlpha.air 0.3;
        maxPartlyDispersedAlpha.air 0.5;
        maxFullyDispersedAlpha.water 0.3;
@ -126,11 +125,13 @@ heatTransfer
    (air in water)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }

    (water in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }
 );

--- a/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/constant/phaseProperties
@ -41,7 +41,6 @@ blending
    default
    {
        type            linear;
-        residualAlpha   1e-6;
        maxFullyDispersedAlpha.air 0.3;
        maxPartlyDispersedAlpha.air 0.5;
        maxFullyDispersedAlpha.water 0.3;
@ -126,11 +125,13 @@ heatTransfer
    (air in water)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }

    (water in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }
 );

--- a/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/RAS/fluidisedBed/constant/phaseProperties
@ -42,7 +42,6 @@ blending
    default
    {
        type            none;
-        residualAlpha   1e-3;
        continuousPhase air;
    }
 }
@ -72,11 +71,6 @@ drag

 virtualMass
 (
-    (particles in air)
-    {
-        type            constantCoefficient;
-        Cvm             0;
-    }
 );

 heatTransfer
@ -84,6 +78,7 @@ heatTransfer
    (particles in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-3;
    }
 );

--- a/tutorials/multiphase/twoPhaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/laminar/bubbleColumn/constant/phaseProperties
@ -41,7 +41,6 @@ blending
    default
    {
        type            linear;
-        residualAlpha   1e-6;
        maxFullyDispersedAlpha.air 0.3;
        maxPartlyDispersedAlpha.air 0.5;
        maxFullyDispersedAlpha.water 0.3;
@ -126,11 +125,13 @@ heatTransfer
    (air in water)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }

    (water in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }
 );

--- a/tutorials/multiphase/twoPhaseEulerFoam/laminar/bubbleColumnIATE/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/laminar/bubbleColumnIATE/constant/phaseProperties
@ -62,7 +62,6 @@ blending
    default
    {
        type            linear;
-        residualAlpha   1e-6;
        maxFullyDispersedAlpha.air 0.3;
        maxPartlyDispersedAlpha.air 0.5;
        maxFullyDispersedAlpha.water 0.3;
@ -147,11 +146,13 @@ heatTransfer
    (air in water)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }

    (water in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }
 );

--- a/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/laminar/fluidisedBed/constant/phaseProperties
@ -42,7 +42,6 @@ blending
    default
    {
        type            none;
-        residualAlpha   1e-3;
        continuousPhase air;
    }
 }
@ -84,6 +83,7 @@ heatTransfer
    (particles in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-3;
    }
 );

--- a/tutorials/multiphase/twoPhaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties
+++ b/tutorials/multiphase/twoPhaseEulerFoam/laminar/mixerVessel2D/constant/phaseProperties
@ -41,7 +41,6 @@ blending
    default
    {
        type            linear;
-        residualAlpha   1e-6;
        maxFullyDispersedAlpha.air 0.3;
        maxPartlyDispersedAlpha.air 0.5;
        maxFullyDispersedAlpha.water 0.3;
@ -126,11 +125,13 @@ heatTransfer
    (air in water)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }

    (water in air)
    {
        type            RanzMarshall;
+        residualAlpha   1e-6;
    }
 );