applications/solvers/multiphase/interPhaseChangeFoam/phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2011 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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 Class
  25     Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer
  26
  27 Description
  28     SchnerrSauer cavitation model.
  29
  30     Reference:
  31     \verbatim
  32         Schnerr, G. H., And Sauer, J.,
  33         "Physical and Numerical Modeling of Unsteady Cavitation Dynamics",
  34         Proc. 4th International Conference on Multiphase Flow,
  35         New Orleans, U.S.A., 2001.
  36     \endverbatim
  37
  38 SourceFiles
  39     SchnerrSauer.C
  40
  41 \*--------------------------------------------------------------------*/
  42
  43 #ifndef SchnerrSauer_H
  44 #define SchnerrSauer_H
  45
  46 #include "phaseChangeTwoPhaseMixture.H"
  47
  48 // * * * * * * * * * * * * * * * * *  * * * * * * * * * * * * * * * * //
  49
  50 namespace Foam
  51 {
  52 namespace phaseChangeTwoPhaseMixtures
  53 {
  54
  55 /*--------------------------------------------------------------------*\
  56                               Class SchnerrSauer
  57 \*--------------------------------------------------------------------*/
  58
  59 class SchnerrSauer
  60 :
  61     public phaseChangeTwoPhaseMixture
  62 {
  63     // Private data
  64
  65         //- Bubble number density
  66         dimensionedScalar n_;
  67
  68         //- Nucleation site diameter
  69         dimensionedScalar dNuc_;
  70
  71         //- Condensation rate coefficient
  72         dimensionedScalar Cc_;
  73
  74         //- Vapourisation rate coefficient
  75         dimensionedScalar Cv_;
  76
  77         dimensionedScalar p0_;
  78
  79         //- Nucleation site volume-fraction
  80         dimensionedScalar alphaNuc() const;
  81
  82         //- Reciprocal bubble radius
  83         tmp<volScalarField>rRb(const volScalarField& limitedAlpha1) const;
  84
  85         //- Part of the condensation and vapourisation rates
  86         tmp<volScalarField> pCoeff(const volScalarField& p) const;
  87
  88
  89 public:
  90
  91     //- Runtime type information
  92     TypeName("SchnerrSauer");
  93
  94
  95     // Constructors
  96
  97         //- construct from components
  98         SchnerrSauer
  99         (
 100             const volVectorField& U,
 101             const surfaceScalarField& phi,
 102             const word& alpha1Name = "alpha1"
 103         );
 104
 105
 106     //- Destructor
 107     virtual ~SchnerrSauer()
 108     {}
 109
 110
 111     // Member Functions
 112
 113         //- Return the mass condensation and vaporisation rates as a
 114         //  coefficient to multiply (1 - alphal) for the condensation rate
 115         //  and a coefficient to multiply  alphal for the vaporisation rate
 116         virtual Pair<tmp<volScalarField> > mDotAlphal() const;
 117
 118         //- Return the mass condensation and vaporisation rates as an
 119         //  explicit term for the condensation rate and a coefficient to
 120         //  multiply (p - pSat) for the vaporisation rate
 121         virtual Pair<tmp<volScalarField> > mDotP() const;
 122
 123         //- Correct the SchnerrSauer phaseChange model
 124         virtual void correct();
 125
 126         //- Read the transportProperties dictionary and update
 127         virtual bool read();
 128 };
 129
 130
 131 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 132
 133 } // End namespace phaseChangeTwoPhaseMixtures
 134 } // End namespace Foam
 135
 136 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 137
 138 #endif
 139
 140 // ************************************************************************* //