icoPolynomialI.H
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-2020 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 "icoPolynomial.H"
27 
28 // * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
29 
30 template<class Specie, int PolySize>
32 (
33  const Specie& sp,
34  const Polynomial<PolySize>& rhoCoeffs
35 )
36 :
37  Specie(sp),
38  rhoCoeffs_(rhoCoeffs)
39 {}
40 
41 
42 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
43 
44 template<class Specie, int PolySize>
46 (
47  const word& name,
49 )
50 :
51  Specie(name, ip),
52  rhoCoeffs_(ip.rhoCoeffs_)
53 {}
54 
55 
56 template<class Specie, int PolySize>
59 {
61  (
63  );
64 }
65 
66 
67 template<class Specie, int PolySize>
70 {
72  (
74  );
75 }
76 
77 
78 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
79 
80 template<class Specie, int PolySize>
82 (
83  scalar p,
84  scalar T
85 ) const
86 {
87  return rhoCoeffs_.value(T);
88 }
89 
90 
91 template<class Specie, int PolySize>
93 (
94  scalar p,
95  scalar T
96 ) const
97 {
98  return p/this->rho(p, T);
99 }
100 
101 
102 template<class Specie, int PolySize>
104 (
105  scalar p,
106  scalar T
107 ) const
108 {
109  return 0;
110 }
111 
112 
113 template<class Specie, int PolySize>
115 (
116  scalar p,
117  scalar T
118 ) const
119 {
120  return 0;
121 }
122 
123 
124 template<class Specie, int PolySize>
126 (
127  scalar p,
128  scalar T
129 ) const
130 {
131  return 0;
132 }
133 
134 
135 template<class Specie, int PolySize>
137 (
138  scalar p,
139  scalar T
140 ) const
141 {
142  return 0;
143 }
144 
145 
146 template<class Specie, int PolySize>
148 (
149  scalar p,
150  scalar T
151 ) const
152 {
153  return 0;
154 }
155 
156 
157 template<class Specie, int PolySize>
159 (
160  scalar p,
161  scalar T
162 ) const
163 {
164  return 0;
165 }
166 
167 
168 template<class Specie, int PolySize>
170 (
171  scalar p,
172  scalar T
173 ) const
174 {
175  return 0;
176 }
177 
178 
179 template<class Specie, int PolySize>
181 (
182  scalar p,
183  scalar T
184 ) const
185 {
186  return 0;
187 }
188 
189 
190 // * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
191 
192 template<class Specie, int PolySize>
193 inline void Foam::icoPolynomial<Specie, PolySize>::operator+=
194 (
196 )
197 {
199 }
200 
201 
202 template<class Specie, int PolySize>
204 {
205  Specie::operator*=(s);
206 }
207 
208 
209 // * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
210 
211 template<class Specie, int PolySize>
213 (
216 )
217 {
219  return ip1;
220 }
221 
222 
223 template<class Specie, int PolySize>
225 (
226  const scalar s,
228 )
229 {
231  (
232  s*static_cast<const Specie&>(ip),
233  ip.rhoCoeffs_
234  );
235 }
236 
237 
238 template<class Specie, int PolySize>
240 (
243 )
244 {
246  return ip1;
247 }
248 
249 
250 // ************************************************************************* //
scalar H(const scalar p, const scalar T) const
Return enthalpy contribution [J/kg].
dictionary dict
scalar psi(scalar p, scalar T) const
Return compressibility [s^2/m^2].
autoPtr< icoPolynomial > clone() const
Construct and return a clone.
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
scalar Sp(const scalar p, const scalar T) const
Return entropy contribution to the integral of Cp/T [J/kg/K].
scalar E(const scalar p, const scalar T) const
Return internal energy contribution [J/kg].
Incompressible, polynomial form of equation of state, using a polynomial function for density...
Definition: icoPolynomial.H:82
scalar Z(scalar p, scalar T) const
Return compression factor [].
scalar Cp(scalar p, scalar T) const
Return Cp contribution [J/(kg K].
static autoPtr< icoPolynomial > New(const dictionary &dict)
void operator*=(const scalar)
scalar CpMCv(scalar p, scalar T) const
Return (Cp - Cv) [J/(kg K].
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
icoPolynomial(const Specie &sp, const Polynomial< PolySize > &rhoPoly)
Construct from components.
A class for handling words, derived from string.
Definition: word.H:59
scalar Cv(scalar p, scalar T) const
Return Cv contribution [J/(kg K].
scalar Sv(const scalar p, const scalar T) const
Return entropy contribution to the integral of Cv/T [J/kg/K].
Polynomial templated on size (order):
Definition: Polynomial.H:65
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:52
scalar rho(scalar p, scalar T) const
Return density [kg/m^3].
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
Definition: error.H:366