Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 namespace Foam
00029 {
00030
00031
00032
00033
00034 template<class thermo>
00035 inline constTransport<thermo>::constTransport
00036 (
00037 const thermo& t,
00038 const scalar mu,
00039 const scalar Pr
00040 )
00041 :
00042 thermo(t),
00043 Mu(mu),
00044 rPr(1.0/Pr)
00045 {}
00046
00047
00048
00049 template<class thermo>
00050 inline constTransport<thermo>::constTransport
00051 (
00052 const word& name,
00053 const constTransport& ct
00054 )
00055 :
00056 thermo(name, ct),
00057 Mu(ct.Mu),
00058 rPr(ct.rPr)
00059 {}
00060
00061
00062
00063
00064
00065 template<class thermo>
00066 inline scalar constTransport<thermo>::mu(const scalar) const
00067 {
00068 return Mu;
00069 }
00070
00071
00072
00073 template<class thermo>
00074 inline scalar constTransport<thermo>::kappa(const scalar T) const
00075 {
00076 return this->Cp(T)*mu(T)*rPr;
00077 }
00078
00079
00080
00081 template<class thermo>
00082 inline scalar constTransport<thermo>::alpha(const scalar T) const
00083 {
00084 scalar Cp_ = this->Cp(T);
00085
00086 scalar deltaT = T - specie::Tstd;
00087 scalar CpBar =
00088 (deltaT*(this->H(T) - this->H(specie::Tstd)) + Cp_)/(sqr(deltaT) + 1);
00089
00090 return Cp_*mu(T)*rPr/CpBar;
00091 }
00092
00093
00094
00095
00096 template<class thermo>
00097 inline constTransport<thermo>& constTransport<thermo>::operator=
00098 (
00099 const constTransport<thermo>& ct
00100 )
00101 {
00102 thermo::operator=(ct);
00103
00104 Mu = ct.Mu;
00105 rPr = ct.rPr;
00106
00107 return *this;
00108 }
00109
00110
00111
00112
00113 template<class thermo>
00114 inline constTransport<thermo> operator+
00115 (
00116 const constTransport<thermo>& ct1,
00117 const constTransport<thermo>& ct2
00118 )
00119 {
00120 thermo t
00121 (
00122 static_cast<const thermo&>(ct1) + static_cast<const thermo&>(ct2)
00123 );
00124
00125 scalar molr1 = ct1.nMoles()/t.nMoles();
00126 scalar molr2 = ct2.nMoles()/t.nMoles();
00127
00128 return constTransport<thermo>
00129 (
00130 t,
00131 molr1*ct1.Mu + molr2*ct2.Mu,
00132 molr1*ct1.rPr + molr2*ct2.rPr
00133 );
00134 }
00135
00136
00137 template<class thermo>
00138 inline constTransport<thermo> operator-
00139 (
00140 const constTransport<thermo>& ct1,
00141 const constTransport<thermo>& ct2
00142 )
00143 {
00144 thermo t
00145 (
00146 static_cast<const thermo&>(ct1) - static_cast<const thermo&>(ct2)
00147 );
00148
00149 scalar molr1 = ct1.nMoles()/t.nMoles();
00150 scalar molr2 = ct2.nMoles()/t.nMoles();
00151
00152 return constTransport<thermo>
00153 (
00154 t,
00155 molr1*ct1.Mu - molr2*ct2.Mu,
00156 molr1*ct1.rPr - molr2*ct2.rPr
00157 );
00158 }
00159
00160
00161 template<class thermo>
00162 inline constTransport<thermo> operator*
00163 (
00164 const scalar s,
00165 const constTransport<thermo>& ct
00166 )
00167 {
00168 return constTransport<thermo>
00169 (
00170 s*static_cast<const thermo&>(ct),
00171 ct.Mu,
00172 ct.rPr
00173 );
00174 }
00175
00176
00177 template<class thermo>
00178 inline constTransport<thermo> operator==
00179 (
00180 const constTransport<thermo>& ct1,
00181 const constTransport<thermo>& ct2
00182 )
00183 {
00184 return ct2 - ct1;
00185 }
00186
00187
00188
00189
00190 }
00191
00192