src/thermophysicalModels/properties/liquidProperties/MB/MB.C

   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.
  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 \*---------------------------------------------------------------------------*/
  25
  26 #include "MB.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33     defineTypeNameAndDebug(MB, 0);
  34     addToRunTimeSelectionTable(liquidProperties, MB,);
  35     addToRunTimeSelectionTable(liquidProperties, MB, Istream);
  36     addToRunTimeSelectionTable(liquidProperties, MB, dictionary);
  37 }
  38
  39 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  40
  41 Foam::MB::MB()
  42 :
  43     liquidProperties
  44     (
  45         102.133,
  46         554.5,
  47         3.4734e+6,
  48         0.34,
  49         0.256,
  50         187.35,
  51         1.0102e-1,
  52         375.90,
  53         5.7373e-30,
  54         0.3807,
  55         1.7713e+4
  56     ),
  57     rho_(76.6099633, 0.257, 554.5, 0.2772),
  58     pv_(107.51, -8112.9, -12.77, 9.2919e-06, 2.0),
  59     hl_(554.5, 508307.794738233, 0.392, 0.0, 0.0, 0.0),
  60     Cp_(1135.77394182096, 2.89818178257762, 0.0, 0.0, 0.0, 0.0),
  61     h_(-5255966.14542938, 1135.77394182096, 1.44909089128881, 0.0, 0.0, 0.0),
  62     Cpg_(875.329227575808, 2849.22600922327, 1570.0, 2029.70636327142, 678.3),
  63     B_
  64     (
  65         0.00220496803188,
  66        -2.42184210783978,
  67        -401045.695318849,
  68        -2.85079259397061e+17,
  69        -3.57377145486767e+19
  70     ),
  71     mu_(-12.206, 1141.7, 0.15014, 0.0, 0.0),
  72     mug_(3.733e-07, 0.6177, 256.5, 0.0),
  73     K_(0.2298, -0.0003002, 0.0, 0.0, 0.0, 0.0),
  74     Kg_(1333.1, 0.9962, 12317000000.0, 0.0),
  75     sigma_(554.5, 0.064084, 1.2418, 0.0, 0.0, 0.0),
  76     D_(147.18, 20.1, 102.133, 28.0) // note: Same as nHeptane
  77 {}
  78
  79
  80 Foam::MB::MB
  81 (
  82     const liquidProperties& l,
  83     const NSRDSfunc5& density,
  84     const NSRDSfunc1& vapourPressure,
  85     const NSRDSfunc6& heatOfVapourisation,
  86     const NSRDSfunc0& heatCapacity,
  87     const NSRDSfunc0& enthalpy,
  88     const NSRDSfunc7& idealGasHeatCapacity,
  89     const NSRDSfunc4& secondVirialCoeff,
  90     const NSRDSfunc1& dynamicViscosity,
  91     const NSRDSfunc2& vapourDynamicViscosity,
  92     const NSRDSfunc0& thermalConductivity,
  93     const NSRDSfunc2& vapourThermalConductivity,
  94     const NSRDSfunc6& surfaceTension,
  95     const APIdiffCoefFunc& vapourDiffussivity
  96 )
  97 :
  98     liquidProperties(l),
  99     rho_(density),
 100     pv_(vapourPressure),
 101     hl_(heatOfVapourisation),
 102     Cp_(heatCapacity),
 103     h_(enthalpy),
 104     Cpg_(idealGasHeatCapacity),
 105     B_(secondVirialCoeff),
 106     mu_(dynamicViscosity),
 107     mug_(vapourDynamicViscosity),
 108     K_(thermalConductivity),
 109     Kg_(vapourThermalConductivity),
 110     sigma_(surfaceTension),
 111     D_(vapourDiffussivity)
 112 {}
 113
 114
 115 Foam::MB::MB(Istream& is)
 116 :
 117     liquidProperties(is),
 118     rho_(is),
 119     pv_(is),
 120     hl_(is),
 121     Cp_(is),
 122     h_(is),
 123     Cpg_(is),
 124     B_(is),
 125     mu_(is),
 126     mug_(is),
 127     K_(is),
 128     Kg_(is),
 129     sigma_(is),
 130     D_(is)
 131 {}
 132
 133
 134 Foam::MB::MB(const dictionary& dict)
 135 :
 136     liquidProperties(dict),
 137     rho_(dict.subDict("rho")),
 138     pv_(dict.subDict("pv")),
 139     hl_(dict.subDict("hl")),
 140     Cp_(dict.subDict("Cp")),
 141     h_(dict.subDict("h")),
 142     Cpg_(dict.subDict("Cpg")),
 143     B_(dict.subDict("B")),
 144     mu_(dict.subDict("mu")),
 145     mug_(dict.subDict("mug")),
 146     K_(dict.subDict("K")),
 147     Kg_(dict.subDict("Kg")),
 148     sigma_(dict.subDict("sigma")),
 149     D_(dict.subDict("D"))
 150 {}
 151
 152
 153 Foam::MB::MB(const MB& liq)
 154 :
 155     liquidProperties(liq),
 156     rho_(liq.rho_),
 157     pv_(liq.pv_),
 158     hl_(liq.hl_),
 159     Cp_(liq.Cp_),
 160     h_(liq.h_),
 161     Cpg_(liq.Cpg_),
 162     B_(liq.B_),
 163     mu_(liq.mu_),
 164     mug_(liq.mug_),
 165     K_(liq.K_),
 166     Kg_(liq.Kg_),
 167     sigma_(liq.sigma_),
 168     D_(liq.D_)
 169 {}
 170
 171
 172 // ************************************************************************* //