Home
last modified time | relevance | path

Searched refs:g1yi (Results 1 – 21 of 21) sorted by relevance

/phasta/phSolver/compressible/
H A De3dc.f1 subroutine e3DC (g1yi, g2yi, g3yi, A0, raLS, argument
34 dimension g1yi(npro,nflow), g2yi(npro,nflow),
58 A0gyi( :,1) = A0(:,1,1)*g1yi(:,1)
59 & + A0(:,1,2)*g1yi(:,2)
60 & + A0(:,1,3)*g1yi(:,3)
61 & + A0(:,1,4)*g1yi(:,4)
62 & + A0(:,1,5)*g1yi(:,5)
63 A0gyi( :,2) = A0(:,2,1)*g1yi(:,1)
64 & + A0(:,2,2)*g1yi(:,2)
65 & + A0(:,2,3)*g1yi(:,3)
[all …]
H A De3qvar.f2 & xl, g1yi, g2yi, g3yi, argument
43 & g1yi(npro,nflow), g2yi(npro,nflow),
109 g1yi = zero
120 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * ycl(:,n,1)
121 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * ycl(:,n,2)
122 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * ycl(:,n,3)
123 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * ycl(:,n,4)
124 g1yi(:,5) = g1yi(:,5) + shg(:,n,1) * ycl(:,n,5)
H A De3ql.f32 dimension g1yi(npro,nflow), g2yi(npro,nflow),
75 & rho, xl, g1yi,
95 qdi(:,1) = rlm2mu * g1yi(:,2)
98 qdi(:,2) = rmu * g1yi(:,3)
100 qdi(:,3) = rmu * g1yi(:,4)
102 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
103 & + rmu * u3 * g1yi(:,4)
106 & + con * g1yi(:,5)
111 qdi(:,5) = rmu * g1yi(:,3)
113 qdi(:,6) = rlm * g1yi(:,2)
[all …]
H A De3q.f33 dimension g1yi(npro,nflow), g2yi(npro,nflow),
81 & rho, xl, g1yi,
104 qdi(:,1) = rlm2mu * g1yi(:,2)
107 qdi(:,2) = rmu * g1yi(:,3)
109 qdi(:,3) = rmu * g1yi(:,4)
111 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
112 & + rmu * u3 * g1yi(:,4)
115 & + con * g1yi(:,5)
120 qdi(:,5) = rmu * g1yi(:,3)
122 qdi(:,6) = rlm * g1yi(:,2)
[all …]
H A De3ivar.f4 & g1yi, g2yi, g3yi, argument
68 & aci(npro,nflow), g1yi(npro,nflow),
247 g1yi = zero
260 g1yi(:,1) = g1yi(:,1) + shg(:,1,1) * yl(:,1,1)
265 g1yi(:,2) = g1yi(:,2) + shg(:,1,1) * yl(:,1,2)
270 g1yi(:,3) = g1yi(:,3) + shg(:,1,1) * yl(:,1,3)
275 g1yi(:,4) = g1yi(:,4) + shg(:,1,1) * yl(:,1,4)
280 g1yi(:,5) = g1yi(:,5) + shg(:,1,1) * yl(:,1,5)
337 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1)
338 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
[all …]
H A De3conv.f1 subroutine e3conv (g1yi, g2yi, g3yi, argument
49 dimension g1yi(npro,nflow), g2yi(npro,nflow),
101 & A1(:,1,1) * g1yi(:,1)
102 & + A1(:,1,2) * g1yi(:,2)
105 & + A1(:,1,5) * g1yi(:,5)
117 & A1(:,2,1) * g1yi(:,1)
118 & + A1(:,2,2) * g1yi(:,2)
121 & + A1(:,2,5) * g1yi(:,5)
133 & A1(:,3,1) * g1yi(:,1)
134 & + A1(:,3,2) * g1yi(:,2)
[all …]
H A Dvorticity.f100 real*8 g1yi(npro,nflow) local
125 g1yi = zero
131 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2) ! du/dx
135 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3) ! dv/dx
139 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4) ! dw/dx
158 & (g3yi(:,2)-g1yi(:,4)) !du/dz-dw/dx
162 & (g1yi(:,3)-g2yi(:,2)) !dv/dx-du/dy
168 & +0.5*(g3yi(:,2)-g1yi(:,4))**2 !du/dz-dw/dx
169 & +0.5*(g1yi(:,3)-g2yi(:,2))**2 !dv/dx-du/dy
170 & -g1yi(:,2)**2-g2yi(:,3)**2-g3yi(:,4)**2
[all …]
H A De3visc.f1 subroutine e3visc (g1yi, g2yi, g3yi, argument
42 dimension g1yi(npro,nflow), g2yi(npro,nflow),
278 rmi(:,2) = rlm2mu * g1yi(:,2)
282 rmi(:,3) = rmu * g1yi(:,3)
285 rmi(:,4) = rmu * g1yi(:,4)
288 rmi(:,5) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3)
289 & + rmu * u3 * g1yi(:,4)
292 & + con * g1yi(:,5)
303 rmi(:, 7) = rmu * g1yi(:,3)
306 rmi(:, 8) = rlm * g1yi(:,2)
[all …]
H A De3b.f63 dimension g1yi(npro,nflow), g2yi(npro,nflow),
129 & lnode, g1yi,
198 tau1n = bnorm(:,1) * (rlm2mu* g1yi(:,2) + rlm *g2yi(:,3)
200 & + bnorm(:,2) * (rmu *(g2yi(:,2) + g1yi(:,3)))
201 & + bnorm(:,3) * (rmu *(g3yi(:,2) + g1yi(:,4)))
202 tau2n = bnorm(:,1) * (rmu *(g2yi(:,2) + g1yi(:,3)))
203 & + bnorm(:,2) * (rlm * g1yi(:,2) + rlm2mu*g2yi(:,3)
206 tau3n = bnorm(:,1) * (rmu *(g3yi(:,2) + g1yi(:,4)))
208 & + bnorm(:,3) * (rlm * g1yi(:,2) + rlm *g2yi(:,3)
223 heat = -con * ( bnorm(:,1) * g1yi(:,5) +
H A De3.f63 dimension g1yi(npro,nflow), g2yi(npro,nflow),
132 & g1yi, g2yi, g3yi,
163 call e3conv (g1yi, g2yi, g3yi,
177 call e3visc (g1yi, g2yi, g3yi,
218 call e3dc (g1yi, g2yi, g3yi,
H A De3bvar.f2 & xlb, lnode, g1yi, g2yi, argument
58 & lnode(27), g1yi(npro,nflow),
256 g1yi(:,2) = dxidxb(:,1,1) * gl1yi(:,2) +
266 g1yi(:,3) = dxidxb(:,1,1) * gl1yi(:,3) +
276 g1yi(:,4) = dxidxb(:,1,1) * gl1yi(:,4) +
286 g1yi(:,5) = dxidxb(:,1,1) * gl1yi(:,5) +
/phasta/phSolver/incompressible/
H A De3res.f3 & g1yi, g2yi, g3yi, argument
43 & g1yi(npro,nflow), g2yi(npro,nflow), g3yi(npro,nflow),
102 tmp = -pres + tauC * (g1yi(:,2) + g2yi(:,3) + g3yi(:,4))
103 tmp1 = rmu * ( g2yi(:,2) + g1yi(:,3) )
105 tmp3 = rmu * ( g1yi(:,4) + g3yi(:,2) )
113 & + ubar(:,1) * g1yi(:,2)
117 & + ubar(:,1) * g1yi(:,3)
121 & + ubar(:,1) * g1yi(:,4)
125 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp
139 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp - u1(:)*u1(:)*rho(:)
[all …]
H A De3ivar.f3 & aci, g1yi, g2yi, argument
58 & aci(npro,nsd), g1yi(npro,ndof),
149 g1yi = zero
153 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1)
154 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
155 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
156 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4)
247 & g1yi, g2yi, g3yi,
343 & g1yi(npro,nflow), g2yi(npro,nflow), local
429 g1yi = zero
[all …]
H A De3qvar.f2 & xl, g1yi, argument
39 & g1yi(npro,nflow), g2yi(npro,nflow),
96 g1yi = zero
122 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
123 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
124 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4)
H A De3qGradV.f30 dimension g1yi(npro,nflow), g2yi(npro,nflow),
55 & xl, g1yi,
63 ql(:,i,1 ) = ql(:,i,1 )+ shape(:,i)*WdetJ*g1yi(:,2 ) ! du/dx
67 ql(:,i,4 ) = ql(:,i,4 )+ shape(:,i)*WdetJ*g1yi(:,3 ) ! dv/dx
71 ql(:,i,7 ) = ql(:,i,7 )+ shape(:,i)*WdetJ*g1yi(:,4 ) ! dw/dx
H A De3q.f32 dimension g1yi(npro,nflow), g2yi(npro,nflow),
74 & xl, g1yi,
91 qdi(:,1) = two * rmu * g1yi(:,2)
92 qdi(:,4) = rmu * (g1yi(:,3) + g2yi(:,2))
93 qdi(:,7) = rmu * (g1yi(:,4) + g3yi(:,2))
97 qdi(:,2) = rmu * (g1yi(:,3) + g2yi(:,2))
103 qdi(:,3) = rmu * (g1yi(:,4) + g3yi(:,2))
H A De3ql.f32 dimension g1yi(npro,ndof), g2yi(npro,ndof),
62 & xl, g1yi,
70 qdi(:,1) = two * rmu * g1yi(:,2)
71 qdi(:,4) = rmu * (g1yi(:,3) + g2yi(:,2))
72 qdi(:,7) = rmu * (g1yi(:,4) + g3yi(:,2))
76 qdi(:,2) = rmu * (g1yi(:,3) + g2yi(:,2))
82 qdi(:,3) = rmu * (g1yi(:,4) + g3yi(:,2))
H A De3.f51 dimension g1yi(npro,ndof), g2yi(npro,ndof),
100 & aci, g1yi, g2yi,
119 & g1yi, g2yi, g3yi,
H A De3bvar.f54 & lnode(27), g1yi(npro,ndof),
305 g1yi(:,2) = dxidxb(:,1,1) * gl1yi(:,2) +
315 g1yi(:,3) = dxidxb(:,1,1) * gl1yi(:,3) +
325 g1yi(:,4) = dxidxb(:,1,1) * gl1yi(:,4) +
349 tau1n = bnorm(:,1) * two * rmu * g1yi(:,2)
350 & + bnorm(:,2) * (rmu * (g2yi(:,2) + g1yi(:,3)))
351 & + bnorm(:,3) * (rmu * (g3yi(:,2) + g1yi(:,4)))
352 tau2n = bnorm(:,1) * (rmu * (g2yi(:,2) + g1yi(:,3)))
355 tau3n = bnorm(:,1) * (rmu * (g3yi(:,2) + g1yi(:,4)))
H A DadvLES.f27 & g1yi(npro,ndof), g2yi(npro,ndof),
180 g1yi = zero
184 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1)
185 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
186 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
187 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4)
202 Sij(:,1) = g1yi(:,2)
205 Sij(:,4) = (g2yi(:,2)+g1yi(:,3))*pt5
206 Sij(:,5) = (g3yi(:,2)+g1yi(:,4))*pt5
358 rLui(:,1) =(aci(:,1) + u1 * g1yi(:,2)
[all …]
H A De3b.f61 dimension g1yi(npro,ndof), g2yi(npro,ndof),