OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps46.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "scr06_c.inc"
#include "scr14_c.inc"
#include "com01_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine sigeps46 (nel, nuparam, nuvar, time, timestep, uparam, rho0, rho, volume, eint, epspxx, epspyy, epspzz, epspxy, epspyz, epspzx, depsxx, depsyy, depszz, depsxy, depsyz, depszx, sigoxx, sigoyy, sigozz, sigoxy, sigoyz, sigozx, signxx, signyy, signzz, signxy, signyz, signzx, sigvxx, sigvyy, sigvzz, sigvxy, sigvyz, sigvzx, soundsp, viscmax, uvar, off, ngl, ipt, ipm, mat, wxx, wyy, wzz, mtn, deltax, aire)

Function/Subroutine Documentation

◆ sigeps46()

subroutine sigeps46 ( integer nel,
integer nuparam,
integer nuvar,
time,
timestep,
uparam,
rho0,
rho,
volume,
eint,
epspxx,
epspyy,
epspzz,
epspxy,
epspyz,
epspzx,
depsxx,
depsyy,
depszz,
depsxy,
depsyz,
depszx,
sigoxx,
sigoyy,
sigozz,
sigoxy,
sigoyz,
sigozx,
signxx,
signyy,
signzz,
signxy,
signyz,
signzx,
sigvxx,
sigvyy,
sigvzz,
sigvxy,
sigvyz,
sigvzx,
soundsp,
viscmax,
uvar,
off,
integer, dimension(nel) ngl,
integer ipt,
integer, dimension(npropmi,*) ipm,
integer, dimension(nel) mat,
wxx,
wyy,
wzz,
integer mtn,
deltax,
aire )

Definition at line 30 of file sigeps46.F.

42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE ale_mod
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50#include "mvsiz_p.inc"
51#include "param_c.inc"
52#include "scr06_c.inc"
53#include "scr14_c.inc"
54#include "com01_c.inc"
55C-----------------------------------------------
56C I N P U T A r g u m e n t s
57C-----------------------------------------------
58 INTEGER NEL, NUPARAM, NUVAR,MAT(NEL),NGL(NEL),IPM(NPROPMI,*),IPT
59 my_real time,timestep,uparam(*),
60 . rho(nel),rho0(nel),volume(nel),eint(nel),
61 . epspxx(nel),epspyy(nel),epspzz(nel),
62 . epspxy(nel),epspyz(nel),epspzx(nel),
63 . depsxx(nel),depsyy(nel),depszz(nel),
64 . depsxy(nel),depsyz(nel),depszx(nel),
65 . sigoxx(nel),sigoyy(nel),sigozz(nel),
66 . sigoxy(nel),sigoyz(nel),sigozx(nel),
67 . deltax(nel),aire(nel)
68C-----------------------------------------------
69C O U T P U T A r g u m e n t s
70C-----------------------------------------------
71 my_real
72 . signxx(nel),signyy(nel),signzz(nel),
73 . signxy(nel),signyz(nel),signzx(nel),
74 . sigvxx(nel),sigvyy(nel),sigvzz(nel),
75 . sigvxy(nel),sigvyz(nel),sigvzx(nel),
76 . soundsp(nel),viscmax(nel)
77C-----------------------------------------------
78C I N P U T O U T P U T A r g u m e n t s
79C-----------------------------------------------
80 my_real uvar(nel,nuvar), off(nel),wxx(nel),wyy(nel),wzz(nel)
81C-----------------------------------------------
82C VARIABLES FOR FUNCTION INTERPOLATION
83C-----------------------------------------------
84 INTEGER MTN
85C-----------------------------------------------
86C L o c a l V a r i a b l e s
87C-----------------------------------------------
88 INTEGER I,J,ISGS,IADBUF,K
89 my_real
90 . vis,c1,smag2,
91 . dav(mvsiz),dxx(mvsiz),dyy(mvsiz),dzz(mvsiz),
92 . dxy(mvsiz),dyz(mvsiz),dzx(mvsiz),ss(mvsiz),
93 . ca,vis2(mvsiz),vis2n,viscpression,
94 . vis1(mvsiz),delta(mvsiz),fac,kapa,e,c,vis1o,vis1n,y2p,ypm
95C-----------------------------------------------
96C S o u r c e L i n e s
97C-----------------------------------------------
98C CONSTANT
99C
100 ypm=onzep225
101 kapa=zep4187
102 e=ninep793
103C
104C INITIALIZATION
105C
106 IF(time==zero .AND. ale%GLOBAL%INCOMP==1)THEN
107 DO i=1,nel
108 sigoxx(i)=third*(sigoxx(i)+sigoyy(i)+sigozz(i))
109 sigoyy(i)=sigoxx(i)
110 sigozz(i)=sigoxx(i)
111 ENDDO
112 ENDIF
113C
114 IF(time==zero .AND. mtn==47)THEN
115 DO i=1,nel
116 uvar(i,1)=ypm
117 ENDDO
118 ENDIF
119C-----------------------------------------------
120C VISCOUS FLOW (L.E.S)
121C-----------------------------------------------
122 iadbuf = ipm(7,mat(1))-1
123 vis = uparam(iadbuf+1)
124 c1 = uparam(iadbuf+2)
125 isgs = int(uparam(iadbuf+3))
126 smag2= uparam(iadbuf+4)
127 ca = uparam(iadbuf+5)
128C
129C PRESSURE
130C
131 IF(ale%GLOBAL%INCOMP==1)THEN
132 DO i=1,nel
133 signxx(i)=sigoxx(i)+c1*(depsxx(i)+depsyy(i)+depszz(i))
134 signyy(i)=signxx(i)
135 signzz(i)=signxx(i)
136 signxy(i)=zero
137 signyz(i)=zero
138 signzx(i)=zero
139 soundsp(i) = sqrt(c1/rho(i))
140 ENDDO
141 ELSE
142 DO i=1,nel
143 signxx(i)=c1*(one -rho(i)/rho0(i))
144 signyy(i)=signxx(i)
145 signzz(i)=signxx(i)
146 signxy(i)=zero
147 signyz(i)=zero
148 signzx(i)=zero
149 soundsp(i) = sqrt(c1/rho(i))
150 ENDDO
151 ENDIF
152C
153C VISCOUS STRESS
154C TURBULENT VISCOSITY IN SUC CELL
155C
156 DO i=1,nel
157 dav(i)=third*(epspxx(i)+epspyy(i)+epspzz(i))
158 dxy(i)=half*epspxy(i)
159 dyz(i)=half*epspyz(i)
160 dzx(i)=half*epspzx(i)
161 ss(i)=sqrt(two*(epspxx(i)**2+epspyy(i)**2+epspzz(i)**2+dxy(i)**2+dyz(i)**2+dzx(i)**2))
162 dxx(i)=epspxx(i)-dav(i)
163 dyy(i)=epspyy(i)-dav(i)
164 dzz(i)=epspzz(i)-dav(i)
165 ENDDO
166
167 IF(mtn==46)THEN
168 IF(n2d==0)THEN
169 IF(isgs==3)THEN
170 DO i=1,nel
171 delta(i)=deltax(i)**2
172 ENDDO
173 ELSE
174 DO i=1,nel
175 delta(i)=volume(i)**two_third
176 ENDDO
177 ENDIF
178 ELSE
179 IF(isgs==3)THEN
180 DO i=1,nel
181 delta(i)=deltax(i)**2
182 ENDDO
183 ELSE
184 DO i=1,nel
185 delta(i)=aire(i)
186 ENDDO
187 ENDIF
188 ENDIF
189 DO i=1,nel
190 vis1(i)=vis+smag2*ss(i)*delta(i)
191 vis2(i)=ca*vis1(i)
192 uvar(i,1)=vis1(i)/vis
193 ENDDO
194 ELSE
195
196C BOUNDARY LAYER
197 DO i=1,nel
198 IF(isgs/=0)THEN
199 c=ss(i)*kapa*deltax(i)*deltax(i)/vis
200 y2p=uvar(i,1)
201
202C ITERATION FOR Y2+ ET EQUIVALENT VISCOSITY (SAME AS K-EPSILON WITH DOMAIN CLOSURE)
203
204 y2p=c/log(e*y2p)
205 y2p=max(y2p,ypm)
206 vis1n=kapa*y2p/log(e*y2p)
207 vis2n=zero
208 IF(isgs>=2)vis2n=kapa*y2p
209 uvar(i,1)=y2p
210 ELSE
211 vis1n=one
212 vis2n=zero
213 uvar(i,1)=one
214 ENDIF
215 vis1(i)=vis*vis1n
216 vis2(i)=vis*vis2n
217 ENDDO
218 ENDIF
219C
220 DO i=1,nel
221 vis1(i)=rho(i)*vis1(i)
222 vis2(i)=three*rho(i)*vis2(i)
223 viscpression=vis2(i)*dav(i)
224 viscmax(i) = vis1(i)+half*vis2(i)
225 vis1(i)=2.*vis1(i)
226 sigvxx(i)=vis1(i)*dxx(i)+viscpression
227 sigvyy(i)=vis1(i)*dyy(i)+viscpression
228 sigvzz(i)=vis1(i)*dzz(i)+viscpression
229 sigvxy(i)=vis1(i)*dxy(i)
230 sigvyz(i)=vis1(i)*dyz(i)
231 sigvzx(i)=vis1(i)*dzx(i)
232 ENDDO
233C
234C VORTICITY STORAGE
235C
236 IF((anim_e(10)==1 .OR. anim_se(10)==1) .AND. timestep/=0.)THEN
237 fac=four/timestep
238 IF(n2d==0)THEN
239 DO i=1,nel
240 uvar(i,2)=fac*sqrt(wxx(i)**2+wyy(i)**2+wzz(i)**2)
241 ENDDO
242 ELSE
243 DO i=1,nel
244 uvar(i,2)=fac*wzz(i)
245 ENDDO
246 ENDIF
247 ENDIF
248C
249 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
ale_mod
Definition ale_mod.F:149
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:249