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

   1 /*---------------------------------------------------------------------------*\
   2   ========Merkle=                 |
   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 Generac License as published by the
  13     Free Software Foundation; either  2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the ho it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the imarranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.he GNU General Public License
  19     for more details.
  20
  21     You should have received a copy oNU General Public License
  22     along with OpenFOAM; if not, write Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floorn, MA 02110-1301 USA
  24
  25 Class
  26     Foam::phaseChangeTwoPhaseMixtures::SchnerrSauer
  27
  28 Description
  29     SchnerrSauer cavitation model.
  30
  31     Reference:
  32     @verbatim
  33         Schnerr, G. H., And Sauer, J.,
  34         "Physical and Numerical Modeling of Unsteady Cavitation Dynamics",
  35         Proc. 4th International Conference on Multiphase Flow,
  36         New Orleans, U.S.A., 2001.
  37     @verbatim
  38
  39 SourceFiles
  40     SchnerrSauer.C
  41
  42 \*--------------------------------------------------------------------*/
  43
  44 #ifndef SchnerrSauer_H
  45 #define SchnerrSauer_H
  46
  47 #include "phaseChangeTwoPhaseMixture.H"
  48
  49 // * * * * * * * * * * * * * * * * *  * * * * * * * * * * * * * * * * //
  50
  51 namespace Foam
  52 {
  53 namespace phaseChangeTwoPhaseMixtures
  54 {
  55
  56 /*--------------------------------------------------------------------*\
  57                               Class SchnerrSauer
  58 \*--------------------------------------------------------------------*/
  59
  60 class SchnerrSauer
  61 :
  62     public phaseChangeTwoPhaseMixture
  63 {
  64     // Private data
  65
  66         //- Bubble number density
  67         dimensionedScalar n_;
  68
  69         //- Nucleation site diameter
  70         dimensionedScalar dNuc_;
  71
  72         //- Condensation rate coefficient
  73         dimensionedScalar Cc_;
  74
  75         //- Vapourisation rate coefficient
  76         dimensionedScalar Cv_;
  77
  78         dimensionedScalar p0_;
  79
  80         //- Nucleation site volume-fraction
  81         dimensionedScalar alphaNuc() const;
  82
  83         //- Reciprocal bubble radius
  84         tmp<volScalarField>rRb(const volScalarField& limitedAlpha1) const;
  85
  86         //- Part of the condensation and vapourisation rates
  87         tmp<volScalarField> pCoeff(const volScalarField& p) const;
  88
  89
  90 public:
  91
  92     //- Runtime type information
  93     TypeName("SchnerrSauer");
  94
  95
  96     // Constructors
  97
  98         //- construct from components
  99         SchnerrSauer
 100         (
 101             const volVectorField& U,
 102             const surfaceScalarField& phi,
 103             const word& alpha1Name = "alpha1"
 104         );
 105
 106
 107     // Destructor
 108
 109         virtual ~SchnerrSauer()
 110         {}
 111
 112
 113     // Member Functions
 114
 115         //- Return the mass condensation and vaporisation rates as a
 116         //  coefficient to multiply (1 - alphal) for the condensation rate
 117         //  and a coefficient to multiply  alphal for the vaporisation rate
 118         virtual Pair<tmp<volScalarField> > mDotAlphal() const;
 119
 120         //- Return the mass condensation and vaporisation rates as an
 121         //  explicit term for the condensation rate and a coefficient to
 122         //  multiply (p - pSat) for the vaporisation rate
 123         virtual Pair<tmp<volScalarField> > mDotP() const;
 124
 125         //- Correct the SchnerrSauer phaseChange model
 126         virtual void correct();
 127
 128         //- Read the transportProperties dictionary and update
 129         virtual bool read();
 130 };
 131
 132
 133 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 134
 135 } // End namespace phaseChangeTwoPhaseMixtures
 136 } // End namespace Foam
 137
 138 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 139
 140 #endif
 141
 142 // ************************************************************************* //