polynomialSolidTransportI.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) 2013-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 "specie.H"
27 
28 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
29 
30 template<class Thermo, int PolySize>
33 (
34  const Thermo& t,
35  const Polynomial<PolySize>& kappaCoeffs
36 )
37 :
38  Thermo(t),
39  kappaCoeffs_(kappaCoeffs)
40 {}
41 
42 
43 template<class Thermo, int PolySize>
46 (
47  const word& name,
48  const polynomialSolidTransport& pt
49 )
50 :
51  Thermo(name, pt),
52  kappaCoeffs_(pt.kappaCoeffs_)
53 {}
54 
55 
56 template<class Thermo, int PolySize>
59 {
61  (
63  );
64 }
65 
66 
67 template<class Thermo, int PolySize>
70 {
72  (
74  );
75 }
76 
77 
78 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
79 
80 template<class Thermo, int PolySize>
82 (
83  const scalar p,
84  const scalar T
85 ) const
86 {
88  return scalar(0);
89 }
90 
91 
92 template<class Thermo, int PolySize>
94 (
95  const scalar p,
96  const scalar T
97 ) const
98 {
99  return kappaCoeffs_.value(T);
100 }
101 
102 
103 template<class Thermo, int PolySize>
105 (
106  const scalar p,
107  const scalar T
108 ) const
109 {
110  const scalar kappa(kappaCoeffs_.value(T));
111  return vector(kappa, kappa, kappa);
112 }
113 
114 
115 template<class Thermo, int PolySize>
117 (
118  const scalar p, const scalar T
119 ) const
120 {
121  return kappa(p, T)/this->Cp(p, T);
122 }
123 
124 
125 // * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
126 
127 template<class Thermo, int PolySize>
128 inline void Foam::polynomialSolidTransport<Thermo, PolySize>::operator+=
129 (
131 )
132 {
133  scalar Y1 = this->Y();
134 
135  Thermo::operator+=(pt);
136 
137  Y1 /= this->Y();
138  scalar Y2 = pt.Y()/this->Y();
139 
140  kappaCoeffs_ = Y1*kappaCoeffs_ + Y2*pt.kappaCoeffs_;
141 }
142 
143 
144 template<class Thermo, int PolySize>
145 inline void Foam::polynomialSolidTransport<Thermo, PolySize>::operator*=
146 (
147  const scalar s
148 )
149 {
150  Thermo::operator*=(s);
151 }
152 
153 
154 // * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
155 
156 template<class Thermo, int PolySize>
158 (
161 )
162 {
163  Thermo t
164  (
165  static_cast<const Thermo&>(pt1) + static_cast<const Thermo&>(pt2)
166  );
167 
168  scalar Y1 = pt1.Y()/t.Y();
169  scalar Y2 = pt2.Y()/t.Y();
170 
172  (
173  t,
174  Y1*pt1.kappaCoeffs_ + Y2*pt2.kappaCoeffs_
175  );
176 }
177 
178 
179 template<class Thermo, int PolySize>
181 (
182  const scalar s,
184 )
185 {
187  (
188  s*static_cast<const Thermo&>(pt),
189  pt.kappaCoeffs_
190  );
191 }
192 
193 
194 // ************************************************************************* //
dictionary dict
scalar kappa(const scalar p, const scalar T) const
Thermal conductivity [W/m/K].
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:158
Vector< scalar > vector
A scalar version of the templated Vector.
Definition: vector.H:49
Transport package using polynomial functions for solid kappa.
autoPtr< polynomialSolidTransport > clone() const
Construct and return a clone.
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))
A class for handling words, derived from string.
Definition: word.H:59
scalar mu(const scalar p, const scalar T) const
Dynamic viscosity [kg/m/s].
vector Kappa(const scalar p, const scalar T) const
Thermal conductivity [W/m/K].
scalar alphah(const scalar p, const scalar T) const
Thermal diffusivity of enthalpy [kg/m/s].
PtrList< volScalarField > & Y
scalar Cp(const scalar p, const scalar T) const
Definition: EtoHthermo.H:2
static autoPtr< polynomialSolidTransport > New(const dictionary &dict)
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:52
#define NotImplemented
Issue a FatalErrorIn for a function not currently implemented.
Definition: error.H:366