ENH: autoLayerDriver: better layering information message
[OpenFOAM-2.0.x.git] / src / thermophysicalModels / barotropicCompressibilityModel / Wallis / Wallis.C
blob3aa468840cf71f795f99f14cd355f38d2c7d695a
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4    \\    /   O peration     |
5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
6      \\/     M anipulation  |
7 -------------------------------------------------------------------------------
8 License
9     This file is part of OpenFOAM.
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.
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.
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/>.
24 \*---------------------------------------------------------------------------*/
26 #include "Wallis.H"
27 #include "addToRunTimeSelectionTable.H"
29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
31 namespace Foam
33     namespace compressibilityModels
34     {
35         defineTypeNameAndDebug(Wallis, 0);
36         addToRunTimeSelectionTable
37         (
38             barotropicCompressibilityModel,
39             Wallis,
40             dictionary
41         );
42     }
45 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
47 Foam::compressibilityModels::Wallis::Wallis
49     const dictionary& compressibilityProperties,
50     const volScalarField& gamma,
51     const word& psiName
54     barotropicCompressibilityModel(compressibilityProperties, gamma, psiName),
55     psiv_(compressibilityProperties_.lookup("psiv")),
56     psil_(compressibilityProperties_.lookup("psil")),
57     rhovSat_(compressibilityProperties_.lookup("rhovSat")),
58     rholSat_(compressibilityProperties_.lookup("rholSat"))
60     correct();
64 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
66 void Foam::compressibilityModels::Wallis::correct()
68     psi_ =
69         (gamma_*rhovSat_ + (scalar(1) - gamma_)*rholSat_)
70        *(gamma_*psiv_/rhovSat_ + (scalar(1) - gamma_)*psil_/rholSat_);
74 bool Foam::compressibilityModels::Wallis::read
76     const dictionary& compressibilityProperties
79     barotropicCompressibilityModel::read(compressibilityProperties);
81     compressibilityProperties_.lookup("psiv") >> psiv_;
82     compressibilityProperties_.lookup("psil") >> psil_;
83     compressibilityProperties_.lookup("rhovSat") >> rhovSat_;
84     compressibilityProperties_.lookup("rholSat") >> rholSat_;
86     return true;
90 // ************************************************************************* //