Home
last modified time | relevance | path

Searched refs:rho (Results 1 – 25 of 35) sorted by relevance

12

/phasta/phSolver/compressible/
H A De3mtrx.f1 subroutine e3mtrx (rho, pres, T, argument
49 dimension rho(npro), pres(npro),
87 drdp = rho * betaT
88 drdT = -rho * alfap
92 e3p = rho * ( h + rk)
93 e4p = drdT * (h + rk) + rho * cp
105 A0(:,2,2) = rho
112 A0(:,3,3) = rho
119 A0(:,4,4) = rho
123 A0(:,5,2) = rho * u1
[all …]
H A De3conv.f3 & rho, pres, T, argument
53 & rho(npro), pres(npro),
76 ri(:, 1) = (- u1) * rho
77 ri(:, 2) = (- u1) * rho * u1 - pres
78 ri(:, 3) = (- u1) * rho * u2
79 ri(:, 4) = (- u1) * rho * u3
80 ri(:, 5) = (- u1) * rho * (ei + rk) - u1 * pres
82 ri(:, 6) = (- u2) * rho
83 ri(:, 7) = (- u2) * rho * u1
84 ri(:, 8) = (- u2) * rho * u2 - pres
[all …]
H A De3juel.f29 dimension rk(npro), rho(npro),
59 & rk, rho, ei,
64 yl(:,inod,1) = rho
65 yl(:,inod,2) = rho * yl(:,inod,2)
66 yl(:,inod,3) = rho * yl(:,inod,3)
67 yl(:,inod,4) = rho * yl(:,inod,4)
68 yl(:,inod,5) = rho * (ei + rk)
H A Dgetthm.f1 subroutine getthm (pres, T, Sclr, rk, rho, argument
51 & rho(npro), ei(npro),
76 pres = Rgas * rho * T
111 rho = pres / ( Rgas * T )
138 rho(:) = rsrhol(:)*prop_blend(:)+rsrhog(:)
266 rho = pres / ( Rgas * T )
292 rho(:)=rsrhol(:)*prop_blend(:)
H A De3ivar.f6 & rho, pres, T, argument
72 & rho(npro), pres(npro),
121 & rk, rho, ei,
127 dui(:,1) = rho
128 dui(:,2) = rho * dui(:,2)
129 dui(:,3) = rho * dui(:,3)
130 dui(:,4) = rho * dui(:,4)
131 dui(:,5) = rho * (ei + rk)
220 & rk, rho, ei,
229 call getDiff (T, cp, rho, ycl,
[all …]
H A De3.f69 dimension rho(npro), pres(npro),
134 & rho, pres, T,
149 call e3mtrx (rho, pres, T,
165 & rho, pres, T,
190 & rho, u1, u2,
202 & rho, rmu, cp,
380 dimension rho(npro), rk(npro),
435 & rho, rmu, con,
445 & call e3sourceSclr (Sclr, rho, rmu,
459 call e3mtrxSclr (rho,
[all …]
H A De3b.f73 & rho(npro), pres(npro),
133 & rho, ei, cp,
154 rou = rho * ( un )
156 un = (rou / rho)
190 call getDiff (T, cp, rho, ycl,
443 dimension rho(npro), Sclr(npro),
488 & bnorm, T, rho,
559 rou = rho * ( un )
561 un = (rou / rho)
574 & rlm, rlm2mu, con, rho, Sclr)
[all …]
H A De3source.f2 & rho, u1, u2,
30 & u3(npro), rho(npro),
51 src(:,2) = rho(:) * datmat(1,5,1)
52 src(:,3) = rho(:) * datmat(2,5,1)
53 src(:,4) = rho(:) * datmat(3,5,1)
71 src(:,2) = rho(:) * ampl * delta
185 src(:,2) = src(:,2) + rho(:)*sforce(:,1)
186 src(:,3) = src(:,3) + rho(:)*sforce(:,2)
187 src(:,4) = src(:,4) + rho(:)*sforce(:,3)
189 & + u3*sforce(:,3))*rho(:)
[all …]
H A Dgetdiff.f1 subroutine getDiff (T, cp, rho, ycl, argument
41 & rho(npro), Sclr(npro),
247 xki = rho(e)*abs(savar)/rmu(e)
250 xmut(e) = fv1*abs(savar)*rho(e)
274 & rlm2mu, con, rho, Sclr) argument
308 & rho(npro), Sclr(npro)
H A De3tau.f1 subroutine e3tau (rho, cp, rmu, argument
39 dimension rho(npro), con(npro),
101 fact=pt5*rho*uh1/rmu
114 tau(:,2)=tau(:,2)/rho
162 tau(:,2)=rho**2*((two*dts)**2
169 tau(:,1)=pt125*fact/(rho*(gijd(:,1)+gijd(:,3)+gijd(:,6)))*taucfct
284 subroutine e3tau_nd (rho, cp, rmu, T, argument
324 dimension rho(npro), con(npro),
370 call e3eig1 (rho, T, cp,
587 Tau(:,1,1) = rho*T*(STau(:,1)+u1*STau(:,2)+
[all …]
H A De3qvar.f1 subroutine e3qvar (ycl, shp, shgl, rho, argument
40 dimension ycl(npro,nshl,ndof), rho(npro),
96 & tmp, rho, tmp,
H A De3eig1.f1 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 De3ql.f44 & rmass(npro,nshl,nshl), rho(npro)
75 & rho, xl, g1yi,
87 call getDiff (T, cp, rho, ycl,
H A De3bvar.f4 & u1, u2, u3, rho, ei, argument
63 & u3(npro), rho(npro),
122 & rk, rho, ei,
380 & bnorm, T, rho, cp, rou, argument
423 & u3(npro), rho(npro),
475 & rk, rho, ei,
480 if (iconvsclr.eq.2) rho=one
H A De3q.f40 & con(npro), rho(npro)
81 & rho, xl, g1yi,
92 call getDiff (T, cp, rho, ycl,
H A De3visc.f46 & u3(npro), rho(npro),
362 & rmu, Sclr, rho,
396 & Sclr(npro), rho(npro),
416 stifft(:,1,1)=saSigmaInv*((rmu/rho)+max(zero,Sclr))
/phasta/phSolver/incompressible/
H A De3res.f4 & rLui, rmu, rho, argument
44 & rLui(npro,nsd), rmu(npro), rho(npro),
139 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp - u1(:)*u1(:)*rho(:)
140 rGNa(:,1,2) = tmp1 - u1(:)*u2(:)*rho(:)
141 rGNa(:,1,3) = tmp3 - u1(:)*u3(:)*rho(:)
142 rGNa(:,2,1) = tmp1 - u1(:)*u2(:)*rho(:)
143 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp - u2(:)*u2(:)*rho(:)
144 rGNa(:,2,3) = tmp2 - u3(:)*u2(:)*rho(:)
145 rGNa(:,3,1) = tmp3 - u1(:)*u3(:)*rho(:)
146 rGNa(:,3,2) = tmp2 - u3(:)*u2(:)*rho(:)
[all …]
H A Dgenlmass.f179 real*8 lhsK(9,nnz_tot), gmass(nshg), rho(nshg) local
186 rho=datmat(1,1,1) ! needs to be generalized for VOF or level set
195 tmp=gmass(n)*tmp1*rho(n)
205 rho=rho*gmass*tmp1 ! reuse rho
206 res(:,1)=res(:,1)-ac(:,1)*rho(:)
207 res(:,2)=res(:,2)-ac(:,2)*rho(:)
208 res(:,3)=res(:,3)-ac(:,3)*rho(:)
221 real*8 lhsS(nnz_tot), gmass(nshg), rho(nshg) local
228 rho=datmat(1,1,1) ! needs to be generalized for VOF or level set
237 tmp=gmass(n)*tmp1*rho(n)
[all …]
H A Dgetdiff.f1 subroutine getDiff( dwl,yl, shape, xmudmi, xl, rmu, rho) argument
10 & shape(npro,nshl), rho(npro),
37 rho = datmat(1,1,1) ! single fluid model, i.e., only 1 density
64 rho = datmat(1,1,2) + (datmat(1,1,1)-datmat(1,1,2))*prop_blend
92 call AddEddyViscKE(yl, dwl, shape, rho, sigmaInv, rmu)
135 real*8 diffus(npro), rho(npro) local
138 rho(:) = datmat(1,1,1) ! single fluid model, i.e., only 1 density
151 call AddEddyViscKE(yl, dwl, shape, rho, sigmaInv, diffus)
336 subroutine AddEddyViscKE(yl, dwl, shape, rho, sigmaInv, rmu) argument
345 & rho
[all …]
H A De3.f54 & WdetJ(npro), rho(npro),
94 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho)
102 & WdetJ, rho, pres,
110 call e3stab (rho, u1, u2,
120 & rLui, rmu, rho,
130 & uBar, WdetJ, rho,
H A De3stab.f1 subroutine e3stab (rho, u1, u2, argument
30 dimension rho(npro), u1(npro),
62 rhoinv=one/rho
100 tauC =rho* pt125*fact/(gijd(:,1)+gijd(:,2)+gijd(:,3))*ff
167 tauM = rho ** 2
227 tauM = rho ** 2
329 subroutine e3uBar (rho, ui, dxidx, argument
334 real*8 rho(npro), ui(npro,nsd), local
362 tauM = rho ** 2
495 real*8 rho(npro), uMod(npro,nsd), local
H A De3ivar.f5 & WdetJ, rho, pres, argument
62 & rho(npro), pres(npro),
234 & *gradh(:,1) /rho(:)
236 & *gradh(:,2) /rho(:)
238 & *gradh(:,3) /rho(:)
246 & aci, u1, u2, u3, Temp, rho, xx,
345 & src(npro,nsd), rho(npro), local
450 rho = datmat(1,1,1)
458 & aci, u1, u2, u3, Temp, rho, x,
465 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho)
[all …]
H A De3lhs.f2 & uBar, WdetJ, rho, argument
36 & uBar(npro,3), WdetJ(npro), rho(npro),
59 tmp1 = tlW * rho
63 tsFct = lhmFct * WdetJ * rho
117 & + u3 * shg(:,b,3) )*rho
245 tauM=tauM/rho
282 & WdetJ(npro), rLS(npro), rho(npro), local
/phasta/phSolver/common/
H A Delm3keps.f122 & rho(npro) local
159 rho(:)=datmat(1,1,1)
174 nuInv = rho(e)/rmu(e)
320 mut = rho(e)*CmuKE*fmuKE* k *Tscale
322 mut_k = rho(e)*CmuKE*(
328 mut_eps = rho(e) * CmuKE *k*(
337 pk2 = rho(e)*two3rdq*k*sqrt(ss)
351 src(e,1) = pk - rho(e) * eps
353 & -rho(e)* Ce2 * f2 * eps * Tinv
355 srcRat(e,1) = -kInv*( pk - rho(e) * eps )
[all …]
H A Dbardmc.f74 & rho(nshg),
98 rho(:) = yold(:,1) / (Rgas * yold(:,5)) ! get density at nodes.
100 rho(:) = one ! unit density
174 rls(:,1) = hfres(:,6) - hfres(:,1)*hfres(:,1)/rho(:)
175 rls(:,2) = hfres(:,7) - hfres(:,2)*hfres(:,2)/rho(:)
176 rls(:,3) = hfres(:,8) - hfres(:,3)*hfres(:,3)/rho(:)
184 rls(:,4) = hfres(:,9) - hfres(:,1)*hfres(:,2)/rho(:)
185 rls(:,5) = hfres(:,10) - hfres(:,1)*hfres(:,3)/rho(:)
186 rls(:,6) = hfres(:,11) - hfres(:,2)*hfres(:,3)/rho(:)
681 & u3(npro), rho(npro),
[all …]

12