src/thermophysicalModels/properties/liquidProperties/CH4N2O/CH4N2O.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
  27 #include "CH4N2O.H"
  28 #include "addToRunTimeSelectionTable.H"
  29
  30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34     defineTypeNameAndDebug(CH4N2O, 0);
  35     addToRunTimeSelectionTable(liquidProperties, CH4N2O,);
  36     addToRunTimeSelectionTable(liquidProperties, CH4N2O, Istream);
  37     addToRunTimeSelectionTable(liquidProperties, CH4N2O, dictionary);
  38 }
  39
  40 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  41
  42 Foam::CH4N2O::CH4N2O()
  43 :
  44     liquidProperties
  45     (
  46         60.056,
  47         705.0,
  48         9.050e+6,
  49         0.218,
  50         0.337,
  51         405.85,
  52         9.3131e+1,
  53         465.0,
  54         1.52e-29,
  55         0.3449,
  56         4.7813e+4
  57     ),
  58     rho_(1230.006936, 0.0, 0.0, 0.0, 0.0, 0.0),
  59     pv_(3015.15611544, -185497.059684, -430.223621983, 0.00017405122622, 2.0),
  60     hl_(705.0, 2534249.0, 0.5, 0.0, 0.0, 0.0),
  61     Cp_(2006.46063673904, 0.0, 0.0, 0.0, 0.0, 0.0),
  62     h_(-6154107.41641135, 2006.46063673904, 0.0, 0.0, 0.0, 0.0),
  63     Cpg_(811.875582789397, 2099.04089516451, 1627.3, 1603.63660583455, 724.41),
  64     B_
  65     (
  66        -0.000383641934194752,
  67         0.447249234048222,
  68        -469062.208605302,
  69         5.5628080458239e+18,
  70        -2.3040162514986e+21
  71     ),
  72     mu_(-51.964, 3670.6, 5.7331, -5.3495e-29, 10.0),
  73     mug_(2.6986e-06, 0.498, 1257.7, -19570.0),
  74     K_(-0.4267, 0.0056903, -8.0065e-06, 1.815e-09, 0.0, 0.0),
  75     Kg_(6.977e-05, 1.1243, 844.9, -148850.0),
  76     sigma_(705.0, 1.0, 0.0, 0.0, 0.0, 0.0), // note: set to constant
  77     D_(147.18, 20.1, 60.056, 28.0) // note: Same as nHeptane
  78 {}
  79
  80
  81 Foam::CH4N2O::CH4N2O
  82 (
  83     const liquidProperties& l,
  84     const NSRDSfunc0& density,
  85     const NSRDSfunc1& vapourPressure,
  86     const NSRDSfunc6& heatOfVapourisation,
  87     const NSRDSfunc0& heatCapacity,
  88     const NSRDSfunc0& enthalpy,
  89     const NSRDSfunc7& idealGasHeatCapacity,
  90     const NSRDSfunc4& secondVirialCoeff,
  91     const NSRDSfunc1& dynamicViscosity,
  92     const NSRDSfunc2& vapourDynamicViscosity,
  93     const NSRDSfunc0& thermalConductivity,
  94     const NSRDSfunc2& vapourThermalConductivity,
  95     const NSRDSfunc6& surfaceTension,
  96     const APIdiffCoefFunc& vapourDiffussivity
  97 )
  98 :
  99     liquidProperties(l),
 100     rho_(density),
 101     pv_(vapourPressure),
 102     hl_(heatOfVapourisation),
 103     Cp_(heatCapacity),
 104     h_(enthalpy),
 105     Cpg_(idealGasHeatCapacity),
 106     B_(secondVirialCoeff),
 107     mu_(dynamicViscosity),
 108     mug_(vapourDynamicViscosity),
 109     K_(thermalConductivity),
 110     Kg_(vapourThermalConductivity),
 111     sigma_(surfaceTension),
 112     D_(vapourDiffussivity)
 113 {}
 114
 115
 116 Foam::CH4N2O::CH4N2O(Istream& is)
 117 :
 118     liquidProperties(is),
 119     rho_(is),
 120     pv_(is),
 121     hl_(is),
 122     Cp_(is),
 123     h_(is),
 124     Cpg_(is),
 125     B_(is),
 126     mu_(is),
 127     mug_(is),
 128     K_(is),
 129     Kg_(is),
 130     sigma_(is),
 131     D_(is)
 132 {}
 133
 134
 135 Foam::CH4N2O::CH4N2O(const dictionary& dict)
 136 :
 137     liquidProperties(dict),
 138     rho_(dict.subDict("rho")),
 139     pv_(dict.subDict("pv")),
 140     hl_(dict.subDict("hl")),
 141     Cp_(dict.subDict("Cp")),
 142     h_(dict.subDict("h")),
 143     Cpg_(dict.subDict("Cpg")),
 144     B_(dict.subDict("B")),
 145     mu_(dict.subDict("mu")),
 146     mug_(dict.subDict("mug")),
 147     K_(dict.subDict("K")),
 148     Kg_(dict.subDict("Kg")),
 149     sigma_(dict.subDict("sigma")),
 150     D_(dict.subDict("D"))
 151 {}
 152
 153
 154 Foam::CH4N2O::CH4N2O(const CH4N2O& liq)
 155 :
 156     liquidProperties(liq),
 157     rho_(liq.rho_),
 158     pv_(liq.pv_),
 159     hl_(liq.hl_),
 160     Cp_(liq.Cp_),
 161     h_(liq.h_),
 162     Cpg_(liq.Cpg_),
 163     B_(liq.B_),
 164     mu_(liq.mu_),
 165     mug_(liq.mug_),
 166     K_(liq.K_),
 167     Kg_(liq.Kg_),
 168     sigma_(liq.sigma_),
 169     D_(liq.D_)
 170 {}
 171
 172
 173 // ************************************************************************* //