| /phasta/phSolver/compressible/ |
| H A D | e3visc.f | 3 & rmu, rlm, rlm2mu, 50 dimension rmu(npro), rlm(npro), 70 stiff(:, 3, 3) = rmu 71 stiff(:, 4, 4) = rmu 73 stiff(:, 5, 3) = rmu * u2 74 stiff(:, 5, 4) = rmu * u3 80 stiff(:, 3, 7) = rmu 81 stiff(:, 5, 7) = rmu * u2 87 stiff(:, 4,12) = rmu 88 stiff(:, 5,12) = rmu * u3 [all …]
|
| H A D | getdiff.f | 2 & rmu, rlm, rlm2mu, con, shp, argument 42 & rmu(npro), rlm(npro), 74 rmu = datmat(1,2,2) + (datmat(1,2,1)-datmat(1,2,2)) 91 rmu(:)= datmat(1,2,1) 93 rmu(:)=fmax*datmat(1,2,1) 95 rmu(:)= datmat(1,2,1)*( 113 rmu(:)=datmat(1,2,1) 115 rmu(:)=fmax*datmat(1,2,1) 117 rmu(:)=datmat(1,2,1) 127 rmu = datmat(1,2,1) [all …]
|
| H A D | e3ql.f | 37 & u3(npro), rmu(npro), 88 & rmu, rlm, rlm2mu, con, shape, 98 qdi(:,2) = rmu * g1yi(:,3) 99 & + rmu * g2yi(:,2) 100 qdi(:,3) = rmu * g1yi(:,4) 101 & + rmu * g3yi(:,2) 102 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3) 103 & + rmu * u3 * g1yi(:,4) 104 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 105 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) [all …]
|
| H A D | e3q.f | 38 & u3(npro), rmu(npro), 93 & rmu, rlm, rlm2mu, con, shape, 107 qdi(:,2) = rmu * g1yi(:,3) 108 & + rmu * g2yi(:,2) 109 qdi(:,3) = rmu * g1yi(:,4) 110 & + rmu * g3yi(:,2) 111 qdi(:,4) = rlm2mu * u1 * g1yi(:,2) + rmu * u2 * g1yi(:,3) 112 & + rmu * u3 * g1yi(:,4) 113 & + rmu * u2 * g2yi(:,2) + rlm * u1 * g2yi(:,3) 114 & + rmu * u3 * g3yi(:,2) + rlm * u1 * g3yi(:,4) [all …]
|
| H A D | e3b.f | 84 dimension rmu(npro), rlm(npro), 191 & rmu, rlm, rlm2mu, con, shape, 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))) 205 & + bnorm(:,3) * (rmu *(g3yi(:,3) + g2yi(:,4))) 206 tau3n = bnorm(:,1) * (rmu *(g3yi(:,2) + g1yi(:,4))) 207 & + bnorm(:,2) * (rmu *(g3yi(:,3) + g2yi(:,4))) 449 dimension rmu(npro), rlm(npro), 573 call getDiffSclr (T, cp, rmu, [all …]
|
| H A D | e3.f | 51 & acl(npro,nshl,ndof), rmu(npro), 140 & rmu, rlm, rlm2mu, 179 & rmu, rlm, rlm2mu, 202 & rho, rmu, cp, 271 & rmu, rlm2mu, con) 280 & rmu, rlm2mu, con) 369 & rmu(npro), con(npro), 435 & rho, rmu, con, 445 & call e3sourceSclr (Sclr, rho, rmu, 481 & rmu, Sclr, rho, [all …]
|
| H A D | e3bdg.f | 6 & rmu, rlm2mu, con) argument 22 & rmu(npro), con(npro) 604 A0(:,3,3) = A0(:,3,3) + rmu !rmu 605 A0(:,4,4) = A0(:,4,4) + rmu !rmu 607 A0(:,5,3) = A0(:,5,3) + rmu * u2 608 A0(:,5,4) = A0(:,5,4) + rmu * u3 863 A0(:,2,2) = A0(:,2,2) + rmu !rmu 865 A0(:,4,4) = A0(:,4,4) + rmu !rmu 866 A0(:,5,2) = A0(:,5,2) + rmu * u1 868 A0(:,5,4) = A0(:,5,4) + rmu * u3 [all …]
|
| H A D | e3bdg_nd.f | 6 & rmu, rlm2mu, con) argument 22 & rmu(npro), con(npro), 567 A0(:,3,3) = A0(:,3,3) + rmu !rmu 568 A0(:,4,4) = A0(:,4,4) + rmu !rmu 570 A0(:,5,3) = A0(:,5,3) + rmu * u2 571 A0(:,5,4) = A0(:,5,4) + rmu * u3 782 A0(:,2,2) = A0(:,2,2) + rmu !rmu 784 A0(:,4,4) = A0(:,4,4) + rmu !rmu 785 A0(:,5,2) = A0(:,5,2) + rmu * u1 787 A0(:,5,4) = A0(:,5,4) + rmu * u3 [all …]
|
| H A D | e3ivar.f | 11 & rmu, rlm, rlm2mu, argument 80 & rmu(npro), rlm(npro), 230 & rmu, rlm, rlm2mu, con, shape, 496 & rho, rmu, con, argument 573 & rmu(npro), con(npro), 650 call getDiffSclr(T, cp, rmu,
|
| H A D | e3tau.f | 1 subroutine e3tau (rho, cp, rmu, argument 48 dimension rmu(npro), cv(npro), 101 fact=pt5*rho*uh1/rmu 135 fact=h2o2u*h2o2u*rk*celt/rmu 138 tau(:,3)=min(tau(:,3),fact*rmu/con)*temper 166 & +fff*rmu**2*(gijd(:,1)**2 + gijd(:,3)**2 + gijd(:,6)**2 + 284 subroutine e3tau_nd (rho, cp, rmu, T, argument 333 dimension rmu(npro), cv(npro), 745 subroutine e3tau_nd_modal (rho, cp, rmu, T, 795 dimension rmu(npro), cv(npro), [all …]
|
| H A D | e3ls.f | 2 & rho, rmu, cp, argument 190 call e3tau (rho, cp, rmu, 208 call e3tau_nd (rho, cp, rmu, T, 776 & rho, rmu, rTLSt, argument 833 & acti(npro), rmu(npro), 891 call e3tauSclr(rho, rmu, A0t,
|
| /phasta/phSolver/incompressible/ |
| H A D | getdiff.f | 1 subroutine getDiff( dwl,yl, shape, xmudmi, xl, rmu, rho) argument 9 real*8 yl(npro,nshl,ndof), rmu(npro), xmudmi(npro,ngauss), 38 rmu = datmat(1,2,1) 65 rmu = datmat(1,2,2) + (datmat(1,2,1)-datmat(1,2,2))*prop_blend 86 rmu = rmu + xmudmi(:,intp) 89 call AddEddyViscSA(yl, shape, rmu) 92 call AddEddyViscKE(yl, dwl, shape, rho, sigmaInv, rmu) 97 call EviscDESIC (xl,rmu,xmudmi) 191 subroutine AddEddyViscSA(yl,shape,rmu) 201 & rmu [all …]
|
| H A D | getfld.f | 1 subroutine getFld (T, cp, rmu, rlm, rlm2mu, argument 34 & rmu(npro), rlm(npro), 42 rmu = datmat(1,2,1) 48 rmu = datmat(1,2,1) * (T/datmat(2,2,1))*sqrt(T/datmat(2,2,1)) 57 rlm = -pt66 * rmu 61 rlm = (datmat(1,3,1) - pt66) * rmu 68 rlm2mu = rlm + two * rmu
|
| H A D | e3q.f | 35 & rmu(npro) 87 call getdiff(dwl, yl, shape, xmudmi, xl, rmu, tmp) 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)) 98 qdi(:,5) = two * rmu * g2yi(:,3) 99 qdi(:,8) = rmu * (g2yi(:,4) + g3yi(:,3)) 103 qdi(:,3) = rmu * (g1yi(:,4) + g3yi(:,2)) 104 qdi(:,6)= rmu * (g2yi(:,4) + g3yi(:,3)) [all …]
|
| H A D | e3ql.f | 35 & rmu(npro), 66 call getdiff(dwl, yl, shape, xmudmi, xl,rmu, tmp) 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)) 77 qdi(:,5) = two * rmu * g2yi(:,3) 78 qdi(:,8) = rmu * (g2yi(:,4) + g3yi(:,3)) 82 qdi(:,3) = rmu * (g1yi(:,4) + g3yi(:,2)) 83 qdi(:,6)= rmu * (g2yi(:,4) + g3yi(:,3)) [all …]
|
| H A D | e3stab.f | 3 & rmu, tauC, tauM, argument 35 & rmu(npro), uBar(npro,3), unorm(npro) 63 rnu=rmu*rhoinv 130 tauC=tauM*(one+tauM*rmu*rmu) 169 2 + fff * rmu** 2 190 tauC=tauM*(one+tauM*rmu*rmu) 229 2 + fff * rmu** 2 270 tauC=pt5*tauM*min(one,pt5*tauM/rmu)*taucfct 294 tauM =min(dts,min(fact,fact*fact*unorm*pt33/rmu)) 330 & rLui, rmu, uBar ) argument [all …]
|
| H A D | e3res.f | 4 & rLui, rmu, rho, argument 44 & rLui(npro,nsd), rmu(npro), rho(npro), 103 tmp1 = rmu * ( g2yi(:,2) + g1yi(:,3) ) 104 tmp2 = rmu * ( g3yi(:,3) + g2yi(:,4) ) 105 tmp3 = rmu * ( g1yi(:,4) + g3yi(:,2) ) 125 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp 129 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp 133 rGNa(:,3,3) = two * rmu * g3yi(:,4) + tmp 139 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp - u1(:)*u1(:)*rho(:) 143 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp - u2(:)*u2(:)*rho(:) [all …]
|
| H A D | e3.f | 60 & rmu(npro), tauC(npro), 94 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho) 112 & rmu, tauC, tauM, 120 & rLui, rmu, rho, 131 & rLui, rmu,
|
| H A D | e3lhs.f | 3 & rLui, rmu, argument 37 & rLui(npro,3), rmu(npro), 62 rmu = tlW * rmu 110 t2(:,1) = rmu * shg(:,b,1) 111 t2(:,2) = rmu * shg(:,b,2) 112 t2(:,3) = rmu * shg(:,b,3)
|
| H A D | e3bvar.f | 5 & u1, u2, u3, rmu, argument 51 dimension yl(npro,nshl,ndof), rmu(npro), 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))) 353 & + bnorm(:,2) * two * rmu * g2yi(:,3) 354 & + bnorm(:,3) * (rmu * (g3yi(:,3) + g2yi(:,4))) 355 tau3n = bnorm(:,1) * (rmu * (g3yi(:,2) + g1yi(:,4))) 356 & + bnorm(:,2) * (rmu * (g3yi(:,3) + g2yi(:,4))) [all …]
|
| H A D | e3b.f | 70 dimension rou(npro), rmu(npro) 121 call getdiff(dwl, yl, shape, xmudmi, xlb, rmu, rho) 131 & rmu, unm,
|
| H A D | e3ivar.f | 346 & rmu(npro) local 451 rmu = datmat(1,2,1) 465 call getdiff(dwl, yl, shpfun, xmudmi, xl, rmu, rho) 466 call e3uBar(rho, src, dxidx, rLui, rmu, uBar)
|
| /phasta/phSolver/common/ |
| H A D | hierarchic.f | 369 rmu = datmat(1,2,1) 371 yvl(:,intp,6 ) = two * rmu * gradV(:,1,1) 372 yvl(:,intp,7 ) = two * rmu * gradV(:,2,2) 373 yvl(:,intp,8 ) = two * rmu * gradV(:,3,3) 375 yvl(:,intp,9 ) = rmu * ( gradV(:,1,2) + gradV(:,2,1) ) 376 yvl(:,intp,10) = rmu * ( gradV(:,1,3) + gradV(:,3,1) ) 377 yvl(:,intp,11) = rmu * ( gradV(:,2,3) + gradV(:,3,2) )
|
| H A D | scatnu.f | 22 rmu=datmat(1,2,1) 23 xmudmi=min(xmudmi,1000.0*rmu) !don't let it get larger than 1000 mu
|
| H A D | elm3keps.f | 121 real*8 rmu(npro), local 160 rmu(:)=datmat(1,2,1) 174 nuInv = rho(e)/rmu(e)
|