applications/solvers/multiphase/interPhaseChangeFoam/phaseChangeTwoPhaseMixtures/Merkle/Merkle.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "Merkle.H"
  28 #include "addToRunTimeSelectionTable.H"
  29
  30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34 namespace phaseChangeTwoPhaseMixtures
  35 {
  36     defineTypeNameAndDebug(Merkle, 0);
  37     addToRunTimeSelectionTable(phaseChangeTwoPhaseMixture, Merkle, components);
  38 }
  39 }
  40
  41 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  42
  43 Foam::phaseChangeTwoPhaseMixtures::Merkle::Merkle
  44 (
  45     const volVectorField& U,
  46     const surfaceScalarField& phi,
  47     const word& alpha1Name
  48 )
  49 :
  50     phaseChangeTwoPhaseMixture(typeName, U, phi, alpha1Name),
  51
  52     UInf_(phaseChangeTwoPhaseMixtureCoeffs_.lookup("UInf")),
  53     tInf_(phaseChangeTwoPhaseMixtureCoeffs_.lookup("tInf")),
  54     Cc_(phaseChangeTwoPhaseMixtureCoeffs_.lookup("Cc")),
  55     Cv_(phaseChangeTwoPhaseMixtureCoeffs_.lookup("Cv")),
  56
  57     p0_("0", pSat().dimensions(), 0.0),
  58
  59     mcCoeff_(Cc_/(0.5*sqr(UInf_)*tInf_)),
  60     mvCoeff_(Cv_*rho1()/(0.5*sqr(UInf_)*tInf_*rho2()))
  61 {
  62     correct();
  63 }
  64
  65
  66 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
  67
  68 Foam::Pair<Foam::tmp<Foam::volScalarField> >
  69 Foam::phaseChangeTwoPhaseMixtures::Merkle::mDotAlphal() const
  70 {
  71     const volScalarField& p = alpha1_.db().lookupObject<volScalarField>("p");
  72
  73     return Pair<tmp<volScalarField> >
  74     (
  75         mcCoeff_*max(p - pSat(), p0_),
  76         mvCoeff_*min(p - pSat(), p0_)
  77     );
  78 }
  79
  80 Foam::Pair<Foam::tmp<Foam::volScalarField> >
  81 Foam::phaseChangeTwoPhaseMixtures::Merkle::mDotP() const
  82 {
  83     const volScalarField& p = alpha1_.db().lookupObject<volScalarField>("p");
  84     volScalarField limitedAlpha1 = min(max(alpha1_, scalar(0)), scalar(1));
  85
  86     return Pair<tmp<volScalarField> >
  87     (
  88         mcCoeff_*(1.0 - limitedAlpha1)*pos(p - pSat()),
  89         (-mvCoeff_)*limitedAlpha1*neg(p - pSat())
  90     );
  91 }
  92
  93
  94 void Foam::phaseChangeTwoPhaseMixtures::Merkle::correct()
  95 {}
  96
  97
  98 bool Foam::phaseChangeTwoPhaseMixtures::Merkle::read()
  99 {
 100     if (phaseChangeTwoPhaseMixture::read())
 101     {
 102         phaseChangeTwoPhaseMixtureCoeffs_ = subDict(type() + "Coeffs");
 103
 104         phaseChangeTwoPhaseMixtureCoeffs_.lookup("UInf") >> UInf_;
 105         phaseChangeTwoPhaseMixtureCoeffs_.lookup("tInf") >> tInf_;
 106         phaseChangeTwoPhaseMixtureCoeffs_.lookup("Cc") >> Cc_;
 107         phaseChangeTwoPhaseMixtureCoeffs_.lookup("Cv") >> Cv_;
 108
 109         mcCoeff_ = Cc_/(0.5*sqr(UInf_)*tInf_);
 110         mvCoeff_ = Cv_*rho1()/(0.5*sqr(UInf_)*tInf_*rho2());
 111
 112         return true;
 113     }
 114     else
 115     {
 116         return false;
 117     }
 118 }
 119
 120
 121 // ************************************************************************* //