polynomialTransportI.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-2023 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 "specie.H"
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class Thermo, int PolySize>
32 (
33  const Thermo& t,
34  const Polynomial<PolySize>& muCoeffs,
35  const Polynomial<PolySize>& kappaCoeffs
36 )
37 :
38  Thermo(t),
39  muCoeffs_(muCoeffs),
40  kappaCoeffs_(kappaCoeffs)
41 {}
42 
43 
44 template<class Thermo, int PolySize>
46 (
47  const word& name,
48  const polynomialTransport& pt
49 )
50 :
51  Thermo(name, pt),
52  muCoeffs_(pt.muCoeffs_),
53  kappaCoeffs_(pt.kappaCoeffs_)
54 {}
55 
56 
57 template<class Thermo, int PolySize>
60 {
62  (
64  );
65 }
66 
67 
68 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
69 
70 template<class Thermo, int PolySize>
72 (
73  const scalar p,
74  const scalar T
75 ) const
76 {
77  return muCoeffs_.value(T);
78 }
79 
80 
81 template<class Thermo, int PolySize>
83 (
84  const scalar p,
85  const scalar T
86 ) const
87 {
88  return kappaCoeffs_.value(T);
89 }
90 
91 
92 // * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
93 
94 template<class Thermo, int PolySize>
96 (
98 )
99 {
100  scalar Y1 = this->Y();
101 
102  Thermo::operator+=(pt);
103 
104  if (mag(this->Y()) > small)
105  {
106  Y1 /= this->Y();
107  scalar Y2 = pt.Y()/this->Y();
108 
109  muCoeffs_ = Y1*muCoeffs_ + Y2*pt.muCoeffs_;
110  kappaCoeffs_ = Y1*kappaCoeffs_ + Y2*pt.kappaCoeffs_;
111  }
112 }
113 
114 
115 template<class Thermo, int PolySize>
117 (
118  const scalar s
119 )
120 {
121  Thermo::operator*=(s);
122 }
123 
124 
125 // * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
126 
127 template<class Thermo, int PolySize>
128 inline Foam::polynomialTransport<Thermo, PolySize> Foam::operator+
129 (
132 )
133 {
134  Thermo t
135  (
136  static_cast<const Thermo&>(pt1) + static_cast<const Thermo&>(pt2)
137  );
138 
139  if (mag(t.Y()) < small)
140  {
142  (
143  t,
144  0,
145  pt1.muCoeffs_,
146  pt1.kappaCoeffs_
147  );
148  }
149  else
150  {
151  scalar Y1 = pt1.Y()/t.Y();
152  scalar Y2 = pt2.Y()/t.Y();
153 
154  return polynomialTransport<Thermo, PolySize>
155  (
156  t,
157  Y1*pt1.muCoeffs_ + Y2*pt2.muCoeffs_,
158  Y1*pt1.kappaCoeffs_ + Y2*pt2.kappaCoeffs_
159  );
160  }
161 }
162 
163 
164 template<class Thermo, int PolySize>
165 inline Foam::polynomialTransport<Thermo, PolySize> Foam::operator*
166 (
167  const scalar s,
168  const polynomialTransport<Thermo, PolySize>& pt
169 )
170 {
171  return polynomialTransport<Thermo, PolySize>
172  (
173  s*static_cast<const Thermo&>(pt),
174  pt.muCoeffs_,
175  pt.kappaCoeffs_
176  );
177 }
178 
179 
180 // ************************************************************************* //
Polynomial templated on size (order):
Definition: Polynomial.H:84
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: autoPtr.H:51
Transport package using polynomial functions of temperature for the dynamic viscosity and thermal con...
scalar mu(const scalar p, const scalar T) const
Dynamic viscosity [kg/m/s].
polynomialTransport(const Thermo &t, const Polynomial< PolySize > &muPoly, const Polynomial< PolySize > &kappaPoly)
Construct from components.
scalar kappa(const scalar p, const scalar T) const
Thermal conductivity [W/m/K].
autoPtr< polynomialTransport > clone() const
Construct and return a clone.
A class for handling words, derived from string.
Definition: word.H:62
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.name(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
word name(const bool)
Return a word representation of a bool.
Definition: boolIO.C:39
dimensioned< scalar > mag(const dimensioned< Type > &)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
volScalarField & p
PtrList< volScalarField > & Y