All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
hPolynomialThermo.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2019 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 "hPolynomialThermo.H"
27 #include "IOstreams.H"
28 
29 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
30 
31 template<class EquationOfState, int PolySize>
33 (
34  const dictionary& dict
35 )
36 :
37  EquationOfState(dict),
38  Hf_(dict.subDict("thermodynamics").lookup<scalar>("Hf")),
39  Sf_(dict.subDict("thermodynamics").lookup<scalar>("Sf")),
40  CpCoeffs_
41  (
42  dict.subDict("thermodynamics").lookup
43  (
44  "CpCoeffs<" + Foam::name(PolySize) + '>'
45  )
46  ),
47  hCoeffs_(),
48  sCoeffs_()
49 {
50  hCoeffs_ = CpCoeffs_.integral();
51  sCoeffs_ = CpCoeffs_.integralMinus1();
52 
53  // Offset h poly so that it is relative to the enthalpy at Tstd
54  hCoeffs_[0] += Hf_ - hCoeffs_.value(Tstd);
55 
56  // Offset s poly so that it is relative to the entropy at Tstd
57  sCoeffs_[0] += Sf_ - sCoeffs_.value(Tstd);
58 }
59 
60 
61 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
62 
63 template<class EquationOfState, int PolySize>
65 (
66  Ostream& os
67 ) const
68 {
70 
71  dictionary dict("thermodynamics");
72  dict.add("Hf", Hf_);
73  dict.add("Sf", Sf_);
74  dict.add
75  (
76  word("CpCoeffs<" + Foam::name(PolySize) + '>'),
77  CpCoeffs_
78  );
79  os << indent << dict.dictName() << dict;
80 }
81 
82 
83 // * * * * * * * * * * * * * * * Ostream Operator * * * * * * * * * * * * * //
84 
85 template<class EquationOfState, int PolySize>
86 Foam::Ostream& Foam::operator<<
87 (
88  Ostream& os,
90 )
91 {
92  pt.write(os);
93  return os;
94 }
95 
96 
97 // ************************************************************************* //
dictionary dict
virtual Ostream & write(const char)=0
Write character.
Ostream & indent(Ostream &os)
Indent stream.
Definition: Ostream.H:221
void write(Ostream &os) const
Write to Ostream.
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
const dimensionedScalar & Tstd
Standard temperature.
Thermodynamics package templated on the equation of state, using polynomial functions for cp...
bool add(entry *, bool mergeEntry=false)
Add a new entry.
Definition: dictionary.C:1056
const word dictName() const
Return the local dictionary name (final part of scoped name)
Definition: dictionary.H:123
const dictionary & subDict(const word &) const
Find and return a sub-dictionary.
Definition: dictionary.C:934
Useful combination of include files which define Sin, Sout and Serr and the use of IO streams general...
A class for handling words, derived from string.
Definition: word.H:59
void write(std::ostream &os, const bool binary, List< floatScalar > &fField)
Write floats ascii or binary.
An Ostream is an abstract base class for all output systems (streams, files, token lists...
Definition: Ostream.H:54
word name(const complex &)
Return a string representation of a complex.
Definition: complex.C:47
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Definition: dictionary.C:812