src/thermophysicalModels/liquids/IDEA/IDEA.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "IDEA.H"
  28 #include "addToRunTimeSelectionTable.H"
  29
  30 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34     defineTypeNameAndDebug(IDEA, 0);
  35     addToRunTimeSelectionTable(liquid, IDEA,);
  36     addToRunTimeSelectionTable(liquid, IDEA, Istream);
  37 }
  38
  39 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  40
  41 Foam::IDEA::IDEA()
  42 :
  43     liquid
  44     (
  45         142.26,
  46         618.074,
  47         2.11e+6,
  48         0.523,
  49         0.247,
  50         242.67,
  51         3.4929e-2,
  52         447.3,
  53         1.7012e-30,
  54         0.3478,
  55         1.57e+4
  56     ),
  57     rho_(152.012105, 3.87150382e-1, 618.073893, 4.00790044e-1),
  58     pv_
  59     (
  60         8.4817774623e+01,
  61        -8.6782398353e+03,
  62        -9.1277694857,
  63         4.6153144498e-06,
  64         2.0
  65     ),
  66     hl_
  67     (
  68         618.074,
  69         2.1671983789e+05,
  70        -4.2413153435e+00,
  71         1.1656811532e+01,
  72        -1.1656446689e+01,
  73         4.3667661492
  74     ),
  75     cp_(1.6604957e+3, -6.250871e-1, 6.1778552e-3, 0.0, 0.0, 0.0),
  76     h_(0.0, 0.0, 0.0, 0.0, 0.0, 0.0),
  77     cpg_
  78     (
  79         1.0457515243e+03,
  80         3.4410492875e+03,
  81         1.5976862298e+03,
  82         2.4697705752e+03,
  83         7.3699710536e+02
  84     ),
  85     B_
  86     (
  87         0.00337351091119935,
  88        -4.13606494008504,
  89        -534560.916470464,
  90        -1.13364022911762e+19,
  91         2.80704220402713e+21
  92     ),
  93     mu_(-6.9645853822e+01, 4.4390635942e+03, 8.4680722718e+00, 0.0, 0.0),
  94     mug_(4.2629382158e-08, 8.8144402122e-01, 9.6918097636e+01, 0.0),
  95     K_(2.03684e-01, -2.3168e-04, 0.0, 0.0, 0.0, 0.0),
  96     Kg_
  97     (
  98        -5.664925956707e+02,
  99         8.896721676320e-01,
 100        -2.849783998688e+09,
 101         6.914935658053e+05
 102     ),
 103     sigma_
 104     (
 105         618.074,
 106         8.3846525429e-03,
 107        -1.0044759047e+01,
 108         2.7261918781e+01,
 109        -2.5529134309e+01,
 110         8.6488806234
 111     ),
 112     D_(147.18, 20.1, 142.2, 28.0) // note: Same as nHeptane
 113 {}
 114
 115
 116 Foam::IDEA::IDEA
 117 (
 118     const liquid& l,
 119     const NSRDSfunc5& density,
 120     const NSRDSfunc1& vapourPressure,
 121     const NSRDSfunc6& heatOfVapourisation,
 122     const NSRDSfunc0& heatCapacity,
 123     const NSRDSfunc0& enthalpy,
 124     const NSRDSfunc7& idealGasHeatCapacity,
 125     const NSRDSfunc4& secondVirialCoeff,
 126     const NSRDSfunc1& dynamicViscosity,
 127     const NSRDSfunc2& vapourDynamicViscosity,
 128     const NSRDSfunc0& thermalConductivity,
 129     const NSRDSfunc2& vapourThermalConductivity,
 130     const NSRDSfunc6& surfaceTension,
 131     const APIdiffCoefFunc& vapourDiffussivity
 132 )
 133 :
 134     liquid(l),
 135     rho_(density),
 136     pv_(vapourPressure),
 137     hl_(heatOfVapourisation),
 138     cp_(heatCapacity),
 139     h_(enthalpy),
 140     cpg_(idealGasHeatCapacity),
 141     B_(secondVirialCoeff),
 142     mu_(dynamicViscosity),
 143     mug_(vapourDynamicViscosity),
 144     K_(thermalConductivity),
 145     Kg_(vapourThermalConductivity),
 146     sigma_(surfaceTension),
 147     D_(vapourDiffussivity)
 148 {}
 149
 150
 151 Foam::IDEA::IDEA(Istream& is)
 152 :
 153     liquid(is),
 154     rho_(is),
 155     pv_(is),
 156     hl_(is),
 157     cp_(is),
 158     h_(is),
 159     cpg_(is),
 160     B_(is),
 161     mu_(is),
 162     mug_(is),
 163     K_(is),
 164     Kg_(is),
 165     sigma_(is),
 166     D_(is)
 167 {}
 168
 169
 170 // ************************************************************************* //