| /phasta/phSolver/compressible/ |
| H A D | getthm.f | 1 subroutine getthm (pres, T, Sclr, rk, rho, argument 50 & T(npro), rk(npro), 76 pres = Rgas * rho * T 111 rho = pres / ( Rgas * T ) 135 rsrhog(:) = pres(:) / ( Rgas * T(:)) 148 ei = T * (Rgas / gamma1) 150 tmpg = T * (Rgas / gamma1) !for gas phase 151 tmpl = 3.264*1000*T !cv*T for liquid phase 163 h = T * (Rgas * gamma / gamma1) 168 alfap = one / T [all …]
|
| H A D | e3visc.f | 6 & con, rlsli, compK, T) 55 & f5(npro), f6(npro), T(npro), rk(npro) 230 compK(:, 1) = f1 * T * rlm2mu + f3 * T * rmu 231 & + f6 * T * rmu 233 compK(:, 2) = f2 * T * (rlm + rmu) 234 compK(:, 3) = f1 * T * rmu + f3 * T * rlm2mu 235 & + f6 * T * rmu 237 compK(:, 4) = f4 * T * (rlm + rmu) 238 compK(:, 5) = f5 * T * (rlm + rmu) 239 compK(:, 6) = f1 * T * rmu + f3 * T * rmu [all …]
|
| H A D | e3mtrx.f | 1 subroutine e3mtrx (rho, pres, T, argument 50 & T(npro), ei(npro), 89 A0(:,5,1) = drdp * (h + rk) - alfap * T ! e1p 239 s1 = one/(rho**2 * betaT * T) 240 cv = cp - (alfap**2 * T/rho/betaT) 243 A0DC(:,3) = rho/T 244 A0DC(:,4) = (-rho*alfap)**2 * s1 + (rho*cv/T**2) 249 fact1 = one/(rho*cv*T**2) 250 d = alfap*T/rho/betaT 253 e3bar = e2bar - cv * T [all …]
|
| H A D | e3tau.f | 284 subroutine e3tau_nd (rho, cp, rmu, T, argument 338 & T(npro), pres(npro), alphap(npro), 351 c = sqrt( (gamma * Rgas) * T ) 353 eb1 = cp*T - pt5*(u1**2+u2**2+u3**2) 370 call e3eig1 (rho, T, cp, 585 betaT = cp*T + pt5*(u1**2+u2**2+u3**2) !reuse betaT 587 Tau(:,1,1) = rho*T*(STau(:,1)+u1*STau(:,2)+ 589 Tau(:,1,2) = rho*T*(STau(:,2)+u1*STau(:,3)+ 591 Tau(:,1,3) = rho*T*(STau(:,4)+u1*STau(:,5)+ 593 Tau(:,1,4) = rho*T*(STau(:,7)+u1*STau(:,8)+ [all …]
|
| H A D | getdiff.f | 1 subroutine getDiff (T, cp, rho, ycl, argument 40 dimension T(npro), cp(npro), 156 rmu = datmat(1,2,1) * (T/datmat(2,2,1))*sqrt(T/datmat(2,2,1)) 157 & * ( datmat(2,2,1) + datmat(3,2,1) ) / (T + datmat(3,2,1)) 273 subroutine getDiffSclr (T, cp, rmu, rlm, argument 305 dimension T(npro), cp(npro), 324 rmu = datmat(1,2,1) * (T/datmat(2,2,1))*sqrt(T/datmat(2,2,1)) 325 & * ( datmat(2,2,1) + datmat(3,2,1) ) / (T + datmat(3,2,1))
|
| H A D | itrbc.f | 215 dimension T(nshg) 247 T=y(:,id) 248 call commu (T, ilwork, 1, 'out') 249 y(:,id)=T 250 T=ac(:,id) 251 call commu (T, ilwork, 1, 'out') 252 ac(:,id)=T
|
| H A D | e3qvar.f | 3 & shg, dxidx, WdetJ, T, argument 46 & T(npro), cp(npro), 64 T = zero 76 T = T + shp(:,n) * ycl(:,n,5) 95 call getthm (pres, T, Sclr,
|
| H A D | e3ivar.f | 6 & rho, pres, T, argument 73 & T(npro), ei(npro), 146 T = zero 155 T = T + shape(:,n) * ycl(:,n,5) 219 call getthm (pres, T, Sclr, 229 call getDiff (T, cp, rho, ycl, 492 & T, cp, argument 563 & T(npro), ei(npro), 592 T = zero 605 T = T + shape(:,n) * ycl( :,n,5) [all …]
|
| H A D | e3bvar.f | 3 & g3yi, WdetJb, bnorm, pres, T, argument 61 & pres(npro), T(npro), 93 T = zero 102 T = T + shpb(:,nodlcl) * yl(:,nodlcl,5) 121 call getthm (pres, T, temp, 380 & bnorm, T, rho, cp, rou, argument 421 & pres(npro), T(npro), 449 T = zero 462 T = T + shpb(:,nodlcl) * ycl(:,nodlcl,5) 474 call getthm (pres, T, Sclr,
|
| H A D | e3.f | 70 & T(npro), ei(npro), 134 & rho, pres, T, 149 call e3mtrx (rho, pres, T, 165 & rho, pres, T, 182 & con, rlsli, compK, T) 203 & cv, con, T, 230 rerrl(i,:,6)=(Tmax-Tmin)/T(i) 370 & T(npro), cp(npro), 431 & T, cp,
|
| H A D | MachControl.f90 | 104 real*8 :: u2, T local 109 T = vartsbuff(MC_iPP, 5, ivartsBuff) 114 MC_Mt = sqrt(u2/(1.4*Rgas*T))
|
| H A D | e3b.f | 74 & T(npro), ei(npro), 131 & bnorm, pres, T, 190 call getDiff (T, cp, rho, ycl, 444 & T(npro), cp(npro) 488 & bnorm, T, rho, 573 call getDiffSclr (T, cp, rmu,
|
| H A D | e3eig1.f | 1 subroutine e3eig1 (rho, T, cp, gamb, c, argument 40 dimension rho(npro), T(npro), 173 fact = (one / sqrt( two * rho * T )) / c 199 fact = (one / sqrt( cp * rho )) / T
|
| H A D | e3ql.f | 35 & T(npro), cp(npro), 77 & dxidx, WdetJ, T, 87 call getDiff (T, cp, rho, ycl,
|
| H A D | e3q.f | 36 & T(npro), cp(npro), 83 & dxidx, WdetJ, T, 92 call getDiff (T, cp, rho, ycl,
|
| H A D | e3conv.f | 3 & rho, pres, T, argument 54 & T(npro), ei(npro),
|
| H A D | e3ls.f | 3 & cv, con, T, argument 74 & PTau(npro,5,5), T(npro), 208 call e3tau_nd (rho, cp, rmu, T,
|
| /phasta/phSolver/incompressible/ |
| H A D | getfld.f | 1 subroutine getFld (T, cp, rmu, rlm, rlm2mu, argument 33 dimension T(npro), cp(npro), 48 rmu = datmat(1,2,1) * (T/datmat(2,2,1))*sqrt(T/datmat(2,2,1)) 49 & * ( datmat(2,2,1) + datmat(3,2,1) ) / (T + datmat(3,2,1))
|
| H A D | itrbc.f | 146 dimension T(nshg) 181 T=y(:,id) 182 call commu (T, ilwork, 1, 'out') 183 y(:,id)=T 184 T=ac(:,id) 185 call commu (T, ilwork, 1, 'out') 186 ac(:,id)=T
|
| H A D | les.h | 66 #define Nil(T) ((T)NULL) argument
|
| /phasta/phastaIO/ |
| H A D | phastaIO.h | 156 template<class T> 171 template<class T> 175 const T* const valueArray, 182 PhastaIO_traits<T>::type_string, 186 template<class T> 190 const T* valueArray, 199 PhastaIO_traits<T>::type_string,
|
| /phasta/phSolver/common/ |
| H A D | blowerControl.f90 | 34 real*8 :: T !uniform temperature at blower inlet component 49 & T, nu, argument 58 & T( nBlowers), nu( nblowers), local 77 blower(i)%T = T(i) 388 if(blower(iBlower)%T > 0.0) then 391 BC(imapped, 2) = blower(iBlower)%T
|
| H A D | symtri.c | 46 #define T 0.116786275726379 macro 83 static double wt12[] = {S,S,S,T,T,T,U,U,U,U,U,U};
|
| H A D | input.config | 153 Print Residual at End of Step: False # T lstres=1 F lstres=0 154 Print Error Indicators: False # F ierrcalc = 0, T ierrcalc = 1 156 Print ybar: False # F ioybar = 0, T ioybar = 1 157 Print vorticity: False # F ivort = 0, T ivort = 1 161 Print Velocity Hessian: False # F ihessian = 0, T ihessian = 1 344 Time Varying Boundary Conditions From File: False # F itvn=0 T itvn=1
|
| /phasta/M2NFixBnd/include/ |
| H A D | commonM2NFixBnd.h | 334 c mmut : pointer to table mu = mu (p,T) 335 c mrhot : pointer to table rho = rho (p,T) 336 c mxst : pointer to table xs = xs (p,T)
|