OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
eflux2.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "vect01_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine eflux2 (pm, ixq, v, x, flux, flu1, fill, dfill, vol, ale_connect, alph)

Function/Subroutine Documentation

◆ eflux2()

subroutine eflux2 ( pm,
integer, dimension(nixq,numelq) ixq,
v,
x,
flux,
flu1,
fill,
dfill,
vol,
type(t_ale_connectivity), intent(in) ale_connect,
alph )

Definition at line 33 of file eflux2.F.

35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE ale_connectivity_mod
39 use element_mod , only : nixq
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C G l o b a l P a r a m e t e r s
46C-----------------------------------------------
47#include "mvsiz_p.inc"
48C-----------------------------------------------
49C C o m m o n B l o c k s
50C-----------------------------------------------
51#include "com01_c.inc"
52#include "com04_c.inc"
53#include "vect01_c.inc"
54#include "param_c.inc"
55C-----------------------------------------------
56C D u m m y A r g u m e n t s
57C-----------------------------------------------
58 INTEGER IXQ(NIXQ,NUMELQ)
59 my_real pm(npropm,nummat), v(3,numnod), x(3,numnod), flux(4,*), flu1(*), fill(*), dfill(*), vol(*), alph(*)
60 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
61C-----------------------------------------------
62C L o c a l V a r i a b l e s
63C-----------------------------------------------
64 INTEGER MAT(MVSIZ), NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ), I,II,J,IAD2
65 my_real
66 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
67 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
68 . n1y(mvsiz), n2y(mvsiz), n3y(mvsiz), n4y(mvsiz), n1z(mvsiz), n2z(mvsiz), n3z(mvsiz), n4z(mvsiz),
69 . flux1(mvsiz), flux2(mvsiz), flux3(mvsiz), flux4(mvsiz),
70 . vy1(mvsiz), vy2(mvsiz), vy3(mvsiz), vy4(mvsiz), vz1(mvsiz), vz2(mvsiz), vz3(mvsiz),
71 . vz4(mvsiz), vdy1(mvsiz), vdy2(mvsiz),
72 . vdy3(mvsiz), vdy4(mvsiz), vdz1(mvsiz), vdz2(mvsiz), vdz3(mvsiz), vdz4(mvsiz), upw(mvsiz),
73 . reduc,upwl(4,mvsiz)
74C-----------------------------------------------
75C S o u r c e L i n e s
76C-----------------------------------------------
77 mat(1) = 1
78 DO i=lft,llt
79 ii=i+nft
80 mat(i)=ixq(1,ii)
81 nc1(i)=ixq(2,ii)
82 nc2(i)=ixq(3,ii)
83 nc3(i)=ixq(4,ii)
84 nc4(i)=ixq(5,ii)
85
86 y1(i)=x(2,nc1(i))
87 z1(i)=x(3,nc1(i))
88
89 y2(i)=x(2,nc2(i))
90 z2(i)=x(3,nc2(i))
91
92 y3(i)=x(2,nc3(i))
93 z3(i)=x(3,nc3(i))
94
95 y4(i)=x(2,nc4(i))
96 z4(i)=x(3,nc4(i))
97
98 vdy1(i)=v(2,nc1(i))
99 vdz1(i)=v(3,nc1(i))
100
101 vdy2(i)=v(2,nc2(i))
102 vdz2(i)=v(3,nc2(i))
103
104 vdy3(i)=v(2,nc3(i))
105 vdz3(i)=v(3,nc3(i))
106
107 vdy4(i)=v(2,nc4(i))
108 vdz4(i)=v(3,nc4(i))
109 ENDDO
110C-----------------------------------------------
111C RELATIVE VELOCITY ON FACE
112C-----------------------------------------------
113 DO i=lft,llt
114 vy1(i)=half*(vdy1(i)+vdy2(i))
115 vy2(i)=half*(vdy2(i)+vdy3(i))
116 vy3(i)=half*(vdy3(i)+vdy4(i))
117 vy4(i)=half*(vdy4(i)+vdy1(i))
118 vz1(i)=half*(vdz1(i)+vdz2(i))
119 vz2(i)=half*(vdz2(i)+vdz3(i))
120 vz3(i)=half*(vdz3(i)+vdz4(i))
121 vz4(i)=half*(vdz4(i)+vdz1(i))
122 ENDDO
123C------------------------------------------
124C NORMAL VECTORS
125C------------------------------------------
126 DO i=lft,llt
127 n1y(i)= (z2(i)-z1(i))
128 n1z(i)=-(y2(i)-y1(i))
129
130 n2y(i)= (z3(i)-z2(i))
131 n2z(i)=-(y3(i)-y2(i))
132
133 n3y(i)= (z4(i)-z3(i))
134 n3z(i)=-(y4(i)-y3(i))
135
136 n4y(i)= (z1(i)-z4(i))
137 n4z(i)=-(y1(i)-y4(i))
138 ENDDO
139 IF(n2d == 1)THEN
140 DO i=lft,llt
141 n1y(i) = n1y(i)*(y1(i)+y2(i))*half
142 n1z(i) = n1z(i)*(y1(i)+y2(i))*half
143 n2y(i) = n2y(i)*(y2(i)+y3(i))*half
144 n2z(i) = n2z(i)*(y2(i)+y3(i))*half
145 n3y(i) = n3y(i)*(y3(i)+y4(i))*half
146 n3z(i) = n3z(i)*(y3(i)+y4(i))*half
147 n4y(i) = n4y(i)*(y1(i)+y4(i))*half
148 n4z(i) = n4z(i)*(y1(i)+y4(i))*half
149 ENDDO
150 ENDIF
151C--------------------
152C VOLUME FLUXES
153C--------------------
154 DO i=lft,llt
155 flux1(i)=(vy1(i)*n1y(i)+vz1(i)*n1z(i))
156 flux2(i)=(vy2(i)*n2y(i)+vz2(i)*n2z(i))
157 flux3(i)=(vy3(i)*n3y(i)+vz3(i)*n3z(i))
158 flux4(i)=(vy4(i)*n4y(i)+vz4(i)*n4z(i))
159 ENDDO
160C-----------------------------------------------
161C SPECIFIC CASE LAW51
162C-----------------------------------------------
163 IF(nint(pm(19,mat(1))) == 51)THEN
164 DO i=lft,llt
165 flux(1,i)=flux1(i)
166 flux(2,i)=flux2(i)
167 flux(3,i)=flux3(i)
168 flux(4,i)=flux4(i)
169 ENDDO
170 RETURN
171 ENDIF
172C-----------------------------------------------
173C UPWIND
174C-----------------------------------------------
175 DO j=1,4
176 DO i=lft,llt
177 upwl(j,i) = pm(16,mat(i))
178 ENDDO
179 ENDDO
180
181 !======================================================!
182 ! BOUNDARY FACE : no volume flux by default !
183 ! slip wall bc !
184 !======================================================!
185 DO i=lft,llt
186 iad2 = ale_connect%ee_connect%iad_connect(i + nft)
187 reduc=pm(92,mat(i))
188
189 ii=ale_connect%ee_connect%connected(iad2 + 1 - 1)
190 IF(ii == 0)THEN
191 flux1(i)=flux1(i)*reduc
192 ENDIF
193
194 ii=ale_connect%ee_connect%connected(iad2 + 2 - 1)
195 IF(ii == 0)THEN
196 flux2(i)=flux2(i)*reduc
197 ENDIF
198
199 ii=ale_connect%ee_connect%connected(iad2 + 3 - 1)
200 IF(ii == 0)THEN
201 flux3(i)=flux3(i)*reduc
202 ENDIF
203
204 ii=ale_connect%ee_connect%connected(iad2 + 4 - 1)
205 IF(ii == 0)THEN
206 flux4(i)=flux4(i)*reduc
207 ENDIF
208
209 ENDDO
210C-----------------------------------------------
211C SPECIFIC TREATMENT FOR LAW20 (OBSOLETE)
212C-----------------------------------------------
213 IF(jmult/=0)THEN
214 DO i=lft,llt
215 upw(i) =pm(16,mat(i))
216 ENDDO
217 CALL amulf2(fill,dfill,flux,flu1,vol,ale_connect,alph,
218 . flux1, flux2, flux3, flux4, upw,
219 . nc1, nc2, nc3, nc4)
220 ELSE
221 DO i=lft,llt
222 flux(1,i)= flux1(i) - upwl(1,i) * abs(flux1(i))
223 flux(2,i)= flux2(i) - upwl(2,i) * abs(flux2(i))
224 flux(3,i)= flux3(i) - upwl(3,i) * abs(flux3(i))
225 flux(4,i)= flux4(i) - upwl(4,i) * abs(flux4(i))
226 flu1(i) = flux1(i) + upwl(1,i) * abs(flux1(i))
227 . + flux2(i) + upwl(2,i) * abs(flux2(i))
228 . + flux3(i) + upwl(3,i) * abs(flux3(i))
229 . + flux4(i) + upwl(4,i) * abs(flux4(i))
230 ENDDO
231 ENDIF
232
233 RETURN
amulf2
subroutine amulf2(fill, dfill, flux, flu1, vol, ale_connect, alpv, flux1, flux2, flux3, flux4, upw, nc1, nc2, nc3, nc4)
Definition amulf2.F:34
my_real
#define my_real
Definition cppsort.cpp:32
ale_connectivity_mod
Definition ale_connectivity_mod.F:225