exponentialSolidTransportI.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2013 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 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
27 
28 template<class Thermo>
30 (
31  const Thermo& t,
32  const scalar kappa0,
33  const scalar n0,
34  const scalar Tref
35 )
36 :
37  Thermo(t),
38  kappa0_(kappa0),
39  n0_(n0),
40  Tref_(Tref)
41 {}
42 
43 
44 template<class Thermo>
46 (
47  const word& name,
49 )
50 :
51  Thermo(name, ct),
52  kappa0_(ct.kappa0_),
53  n0_(ct.n0_),
54  Tref_(ct.Tref_)
55 {}
56 
57 
58 template<class Thermo>
61 (
62  const dictionary& dict
63 )
64 {
66  (
68  );
69 }
70 
71 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
72 
73 template<class Thermo>
75 (
76  const scalar p, const scalar T
77 ) const
78 {
79  return (kappa0_*pow(T/Tref_, n0_));
80 }
81 
82 
83 template<class Thermo>
85 (
86  const scalar p, const scalar T
87 ) const
88 {
89  const scalar kappa(kappa0_*pow(T/Tref_, n0_));
90  return vector(kappa, kappa, kappa);
91 }
92 
93 
94 template<class Thermo>
96 mu(const scalar p, const scalar T) const
97 {
99  (
100  "Foam::scalar Foam::exponentialSolidTransport<Thermo>mu::"
101  "("
102  " const scalar p, const scalar T"
103  ") const"
104  );
105  return scalar(0);
106 }
107 
108 
109 template<class Thermo>
111 alphah(const scalar p, const scalar T) const
112 {
113  return kappa(p, T)/this->Cpv(p, T);
114 }
115 
116 // * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
117 
118 
119 template<class Thermo>
121 Foam::exponentialSolidTransport<Thermo>::operator=
122 (
124 )
125 {
126  kappa0_ = ct.kappa0_;
127  n0_ = ct.n0_;
128  Tref_ = ct.Tref_;
129  return *this;
130 }
131 
132 
133 template<class Thermo>
134 inline void Foam::exponentialSolidTransport<Thermo>::operator+=
135 (
137 )
138 {
139  scalar molr1 = this->nMoles();
140 
141  molr1 /= this->nMoles();
142  scalar molr2 = ct.nMoles()/this->nMoles();
143 
144  kappa0_ = molr1*kappa0_ + molr2*ct.kappa0_;
145  n0_ = (molr1*n0_ + molr2*ct.n0_);
146  Tref_ = (molr1*Tref_ + molr2*ct.Tref_);
147 }
148 
149 
150 template<class Thermo>
151 inline void Foam::exponentialSolidTransport<Thermo>::operator-=
152 (
154 )
155 {
156  scalar molr1 = this->nMoles();
157 
158  molr1 /= this->nMoles();
159  scalar molr2 = ct.nMoles()/this->nMoles();
160 
161  kappa0_ = (molr1*kappa0_ - molr2*ct.kappa0_);
162  n0_ = (molr1*n0_ - molr2*ct.n0_);
163  Tref_ = (molr1*Tref_ - molr2*ct.Tref_);
164 }
165 
166 
167 // * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
168 
169 template<class Thermo>
170 inline Foam::exponentialSolidTransport<Thermo> Foam::operator*
171 (
172  const scalar s,
174 )
175 {
177  (
178  s*static_cast<const Thermo&>(ct),
179  ct.kappa0_,
180  ct.n0_,
181  ct.Tref_
182  );
183 }
184 
185 // ************************************************************************* //
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
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 ))
Exponential properties for solid heat transport Templated into a given thermodynamics package...
scalar kappa(const scalar p, const scalar T) const
Thermal conductivity [W/mK].
A class for handling words, derived from string.
Definition: word.H:59
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
scalar mu(const scalar p, const scalar T) const
Dynamic viscosity [kg/ms].
dictionary dict
scalar alphah(const scalar p, const scalar T) const
Thermal diffusivity of enthalpy [kg/ms].
static autoPtr< exponentialSolidTransport > New(const dictionary &dict)
Vector< scalar > vector
A scalar version of the templated Vector.
Definition: vector.H:49
vector Kappa(const scalar p, const scalar T) const
Thermal conductivity [W/mK].
#define notImplemented(functionName)
Issue a FatalErrorIn for a function not currently implemented.
Definition: error.H:356
An auto-pointer similar to the STL auto_ptr but with automatic casting to a reference to the type and...
Definition: PtrList.H:117