multiphase solvers: added support for the multiphaseFixedFluxPressure BC

applications/solvers/multiphase/compressibleTwoPhaseEulerFoam/pEqn.H

@@ -59,12 +59,13 @@
         );
     }
 
-    phi = alpha1f*phi1 + alpha2f*phi2;
+    surfaceScalarField phi0("phi0", alpha1f*phi1 + alpha2f*phi2);
+    phi = phi0;
     adjustPhi(phi, U, p);
 
     surfaceScalarField Dp
     (
-        "(rho*1|A(U))",
+        "Dp",
         mag(alpha1f*rAlphaAU1f/fvc::interpolate(rho1)
       + alpha2f*rAlphaAU2f/fvc::interpolate(rho2))
     );

applications/solvers/multiphase/twoPhaseEulerFoam/pEqn.H

@@ -50,11 +50,12 @@
     phib = (fvc::interpolate(Ub) & mesh.Sf()) + fvc::ddtPhiCorr(rUbA, Ub, phib)
          + phiDragb;
 
-    phi = alphaf*phia + betaf*phib;
+    surfaceScalarField phi0("phi0", alphaf*phia + betaf*phib);
+    phi = phi0;
 
     surfaceScalarField Dp
     (
-        "(rho*(1|A(U)))",
+        "Dp",
         alphaf*rUaAf/rhoa + betaf*rUbAf/rhob
     );
 

src/finiteVolume/Make/files

@@ -127,6 +127,7 @@ $(derivedFvPatchFields)/mappedFlowRate/mappedFlowRateFvPatchVectorField.C
 $(derivedFvPatchFields)/fan/fanFvPatchFields.C
 $(derivedFvPatchFields)/fanPressure/fanPressureFvPatchScalarField.C
 $(derivedFvPatchFields)/buoyantPressure/buoyantPressureFvPatchScalarField.C
+$(derivedFvPatchFields)/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.C
 $(derivedFvPatchFields)/fixedFluxPressure/fixedFluxPressureFvPatchScalarField.C
 $(derivedFvPatchFields)/fixedInternalValueFvPatchField/fixedInternalValueFvPatchFields.C
 $(derivedFvPatchFields)/fixedNormalSlip/fixedNormalSlipFvPatchFields.C

src/finiteVolume/fields/fvPatchFields/derived/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.C

@@ -0,0 +1,181 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*---------------------------------------------------------------------------*/
+
+#include "multiphaseFixedFluxPressureFvPatchScalarField.H"
+#include "fvPatchFieldMapper.H"
+#include "volFields.H"
+#include "surfaceFields.H"
+#include "addToRunTimeSelectionTable.H"
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::multiphaseFixedFluxPressureFvPatchScalarField::
+multiphaseFixedFluxPressureFvPatchScalarField
+(
+    const fvPatch& p,
+    const DimensionedField<scalar, volMesh>& iF
+)
+:
+    fixedGradientFvPatchScalarField(p, iF),
+    phi0Name_("phi0"),
+    phiName_("phi"),
+    rhoName_("rho")
+{}
+
+
+Foam::multiphaseFixedFluxPressureFvPatchScalarField::
+multiphaseFixedFluxPressureFvPatchScalarField
+(
+    const multiphaseFixedFluxPressureFvPatchScalarField& ptf,
+    const fvPatch& p,
+    const DimensionedField<scalar, volMesh>& iF,
+    const fvPatchFieldMapper& mapper
+)
+:
+    fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
+    phi0Name_(ptf.phi0Name_),
+    phiName_(ptf.phiName_),
+    rhoName_(ptf.rhoName_)
+{}
+
+
+Foam::multiphaseFixedFluxPressureFvPatchScalarField::
+multiphaseFixedFluxPressureFvPatchScalarField
+(
+    const fvPatch& p,
+    const DimensionedField<scalar, volMesh>& iF,
+    const dictionary& dict
+)
+:
+    fixedGradientFvPatchScalarField(p, iF),
+    phi0Name_(dict.lookupOrDefault<word>("phi0", "phi0")),
+    phiName_(dict.lookupOrDefault<word>("phi", "phi")),
+    rhoName_(dict.lookupOrDefault<word>("rho", "rho"))
+{
+    if (dict.found("gradient"))
+    {
+        gradient() = scalarField("gradient", dict, p.size());
+        fixedGradientFvPatchScalarField::updateCoeffs();
+        fixedGradientFvPatchScalarField::evaluate();
+    }
+    else
+    {
+        fvPatchField<scalar>::operator=(patchInternalField());
+        gradient() = 0.0;
+    }
+}
+
+
+Foam::multiphaseFixedFluxPressureFvPatchScalarField::
+multiphaseFixedFluxPressureFvPatchScalarField
+(
+    const multiphaseFixedFluxPressureFvPatchScalarField& wbppsf
+)
+:
+    fixedGradientFvPatchScalarField(wbppsf),
+    phi0Name_(wbppsf.phi0Name_),
+    phiName_(wbppsf.phiName_),
+    rhoName_(wbppsf.rhoName_)
+{}
+
+
+Foam::multiphaseFixedFluxPressureFvPatchScalarField::
+multiphaseFixedFluxPressureFvPatchScalarField
+(
+    const multiphaseFixedFluxPressureFvPatchScalarField& wbppsf,
+    const DimensionedField<scalar, volMesh>& iF
+)
+:
+    fixedGradientFvPatchScalarField(wbppsf, iF),
+    phi0Name_(wbppsf.phi0Name_),
+    phiName_(wbppsf.phiName_),
+    rhoName_(wbppsf.rhoName_)
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+void Foam::multiphaseFixedFluxPressureFvPatchScalarField::updateCoeffs()
+{
+    if (updated())
+    {
+        return;
+    }
+
+    const surfaceScalarField& phi0 =
+        db().lookupObject<surfaceScalarField>(phi0Name_);
+
+    const surfaceScalarField& phi =
+        db().lookupObject<surfaceScalarField>(phiName_);
+
+    fvsPatchField<scalar> phi0p =
+        patch().patchField<surfaceScalarField, scalar>(phi0);
+
+    fvsPatchField<scalar> phip =
+        patch().patchField<surfaceScalarField, scalar>(phi);
+
+    if (phi.dimensions() == dimDensity*dimVelocity*dimArea)
+    {
+        const fvPatchField<scalar>& rhop =
+            patch().lookupPatchField<volScalarField, scalar>(rhoName_);
+
+        phip /= rhop;
+    }
+
+    const fvsPatchField<scalar>& Dpp =
+        patch().lookupPatchField<surfaceScalarField, scalar>("Dp");
+
+    gradient() = (phi0p - phip)/patch().magSf()/Dpp;
+
+    fixedGradientFvPatchScalarField::updateCoeffs();
+}
+
+
+void Foam::multiphaseFixedFluxPressureFvPatchScalarField::write
+(
+    Ostream& os
+) const
+{
+    fvPatchScalarField::write(os);
+    writeEntryIfDifferent<word>(os, "phi0", "phi0", phi0Name_);
+    writeEntryIfDifferent<word>(os, "phi", "phi", phiName_);
+    writeEntryIfDifferent<word>(os, "rho", "rho", rhoName_);
+    gradient().writeEntry("gradient", os);
+    writeEntry("value", os);
+}
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    makePatchTypeField
+    (
+        fvPatchScalarField,
+        multiphaseFixedFluxPressureFvPatchScalarField
+    );
+}
+
+// ************************************************************************* //

src/finiteVolume/fields/fvPatchFields/derived/multiphaseFixedFluxPressure/multiphaseFixedFluxPressureFvPatchScalarField.H

@@ -0,0 +1,154 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::multiphaseFixedFluxPressureFvPatchScalarField
+
+Description
+    Foam::multiphaseFixedFluxPressureFvPatchScalarField
+
+SourceFiles
+    multiphaseFixedFluxPressureFvPatchScalarField.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef multiphaseFixedFluxPressureFvPatchScalarFields_H
+#define multiphaseFixedFluxPressureFvPatchScalarFields_H
+
+#include "fvPatchFields.H"
+#include "fixedGradientFvPatchFields.H"
+#include "Switch.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+              Class multiphaseFixedFluxPressureFvPatch Declaration
+\*---------------------------------------------------------------------------*/
+
+class multiphaseFixedFluxPressureFvPatchScalarField
+:
+    public fixedGradientFvPatchScalarField
+{
+    // Private data
+
+        //- Name of the predicted flux transporting the field
+        word phi0Name_;
+
+        //- Name of the flux transporting the field
+        word phiName_;
+
+        //- Name of the density field used to normalise the mass flux
+        //  if neccessary
+        word rhoName_;
+
+
+public:
+
+    //- Runtime type information
+    TypeName("multiphaseFixedFluxPressure");
+
+
+    // Constructors
+
+        //- Construct from patch and internal field
+        multiphaseFixedFluxPressureFvPatchScalarField
+        (
+            const fvPatch&,
+            const DimensionedField<scalar, volMesh>&
+        );
+
+        //- Construct from patch, internal field and dictionary
+        multiphaseFixedFluxPressureFvPatchScalarField
+        (
+            const fvPatch&,
+            const DimensionedField<scalar, volMesh>&,
+            const dictionary&
+        );
+
+        //- Construct by mapping given
+        //  multiphaseFixedFluxPressureFvPatchScalarField onto a new patch
+        multiphaseFixedFluxPressureFvPatchScalarField
+        (
+            const multiphaseFixedFluxPressureFvPatchScalarField&,
+            const fvPatch&,
+            const DimensionedField<scalar, volMesh>&,
+            const fvPatchFieldMapper&
+        );
+
+        //- Construct as copy
+        multiphaseFixedFluxPressureFvPatchScalarField
+        (
+            const multiphaseFixedFluxPressureFvPatchScalarField&
+        );
+
+        //- Construct and return a clone
+        virtual tmp<fvPatchScalarField> clone() const
+        {
+            return tmp<fvPatchScalarField>
+            (
+                new multiphaseFixedFluxPressureFvPatchScalarField(*this)
+            );
+        }
+
+        //- Construct as copy setting internal field reference
+        multiphaseFixedFluxPressureFvPatchScalarField
+        (
+            const multiphaseFixedFluxPressureFvPatchScalarField&,
+            const DimensionedField<scalar, volMesh>&
+        );
+
+        //- Construct and return a clone setting internal field reference
+        virtual tmp<fvPatchScalarField> clone
+        (
+            const DimensionedField<scalar, volMesh>& iF
+        ) const
+        {
+            return tmp<fvPatchScalarField>
+            (
+                new multiphaseFixedFluxPressureFvPatchScalarField(*this, iF)
+            );
+        }
+
+
+    // Member functions
+
+        //- Update the coefficients associated with the patch field
+        virtual void updateCoeffs();
+
+        //- Write
+        virtual void write(Ostream&) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

tutorials/multiphase/compressibleTwoPhaseEulerFoam/bubbleColumn/0/p

@@ -22,16 +22,18 @@ boundaryField
 {
     inlet
     {
-        type               zeroGradient;
+        type               multiphaseFixedFluxPressure;
+        value              $internalField;
     }
     outlet
     {
         type               fixedValue;
-        value              uniform 1e5;
+        value              $internalField;
     }
     walls
     {
-        type               zeroGradient;
+        type               multiphaseFixedFluxPressure;
+        value              $internalField;
     }
 }
 

tutorials/multiphase/twoPhaseEulerFoam/bed/0/p

@@ -22,20 +22,20 @@ boundaryField
 {
     bottom
     {
-        type            buoyantPressure;
+        type               multiphaseFixedFluxPressure;
-        value           uniform 0;
+        value              $internalField;
     }
 
     top
     {
         type               fixedValue;
-        value              uniform 0;
+        value              $internalField;
     }
 
     walls
     {
-        type            buoyantPressure;
+        type               multiphaseFixedFluxPressure;
-        value           uniform 0;
+        value              $internalField;
     }
 
     frontBack

tutorials/multiphase/twoPhaseEulerFoam/bed/system/fvSchemes

@@ -44,7 +44,7 @@ laplacianSchemes
     default         none;
     laplacian(nuEffa,Ua) Gauss linear corrected;
     laplacian(nuEffb,Ub) Gauss linear corrected;
-    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
+    laplacian(Dp,p) Gauss linear corrected;
     laplacian(alphaPpMag,alpha) Gauss linear corrected;
     laplacian(DkEff,k) Gauss linear corrected;
     laplacian(DepsilonEff,epsilon) Gauss linear corrected;

tutorials/multiphase/twoPhaseEulerFoam/bed2/0/p

@@ -22,16 +22,18 @@ boundaryField
 {
     walls
     {
-        type               zeroGradient;
+        type               multiphaseFixedFluxPressure;
+        value              $internalField;
     }
     outlet
     {
         type               fixedValue;
-        value              uniform 0;
+        value              $internalField;
     }
     inlet
     {
-        type               zeroGradient;
+        type               multiphaseFixedFluxPressure;
+        value              $internalField;
     }
     frontAndBackPlanes
     {

tutorials/multiphase/twoPhaseEulerFoam/bed2/system/fvSchemes

@@ -44,7 +44,7 @@ laplacianSchemes
     default         none;
     laplacian(nuEffa,Ua) Gauss linear corrected;
     laplacian(nuEffb,Ub) Gauss linear corrected;
-    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
+    laplacian(Dp,p) Gauss linear corrected;
     laplacian(alphaPpMag,alpha) Gauss linear corrected;
     laplacian(Galphaf,alpha) Gauss linear corrected;
     laplacian(DkEff,k) Gauss linear corrected;

tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/0/p

@@ -1901,7 +1901,8 @@ boundaryField
 {
     inlet
     {
-        type            zeroGradient;
+        type            multiphaseFixedFluxPressure;
+        value           uniform 0;
     }
 
     outlet
@@ -1912,7 +1913,7 @@ boundaryField
 
     walls
     {
-        type            buoyantPressure;
+        type            multiphaseFixedFluxPressure;
         value           uniform 0;
     }
 

tutorials/multiphase/twoPhaseEulerFoam/bubbleColumn/system/fvSchemes

@@ -44,7 +44,7 @@ laplacianSchemes
     default         none;
     laplacian(nuEffa,Ua) Gauss linear corrected;
     laplacian(nuEffb,Ub) Gauss linear corrected;
-    laplacian((rho*(1|A(U))),p) Gauss linear corrected;
+    laplacian(Dp,p) Gauss linear corrected;
     laplacian(alphaPpMag,alpha) Gauss linear corrected;
     laplacian((alphak*nuEffb),k) Gauss linear corrected;
     laplacian((alphaEps*nuEffb),epsilon) Gauss linear corrected;