OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
ale51_upwind3.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| ale51_upwind3 ../engine/source/ale/ale51/ale51_upwind3.F
25!||--- called by ------------------------------------------------------
26!|| afluxt ../engine/source/ale/ale51/afluxt.F
27!|| ale51_finish ../engine/source/ale/ale51/ale51_finish.F
28!|| ale51_init ../engine/source/ale/ale51/ale51_init.F
29!||--- uses -----------------------------------------------------
30!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
31!|| element_mod ../common_source/modules/elements/element_mod.F90
32!|| i22tri_mod ../common_source/modules/interfaces/cut-cell-search_mod.F
33!||====================================================================
34 SUBROUTINE ale51_upwind3(PM ,IXS ,FLUX ,FLU1 ,ALE_CONNECT,
35 . ITRIMAT,DDVOL,QMV ,IFLG ,
36 . NV46)
37C-----------------------------------------------
38C D e s c r i p t i o n
39C-----------------------------------------------
40C This subroutines computes from direct fluxes :
41C -1- FLUX(1:6,*), FLU1(*) : Decomposition Factors from Biased Upwind Method
42C -2- QMV(7:12) :
43C -3- DDVOL : D/DV . D/DT . VOL = DV/DT
44C
45C-----------------------------------------------
46C M o d u l e s
47C-----------------------------------------------
48 USE i22tri_mod
49 USE ale_connectivity_mod
50 use element_mod , only : nixs
51C-----------------------------------------------
52C I m p l i c i t T y p e s
53C-----------------------------------------------
54#include "implicit_f.inc"
55C-----------------------------------------------
56C G l o b a l P a r a m e t e r s
57C-----------------------------------------------
58#include "mvsiz_p.inc"
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "com01_c.inc"
63#include "com04_c.inc"
64#include "vect01_c.inc"
65#include "param_c.inc"
66C-----------------------------------------------
67C D u m m y A r g u m e n t s
68C-----------------------------------------------
69 INTEGER IXS(NIXS,NUMELS), ITRIMAT,IFLG,NV46
70 my_real PM(NPROPM,NUMMAT), FLUX(NV46,*), FLU1(*),DDVOL(*)
71 my_real,TARGET :: qmv(12,*)
72 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER I,J,II,MAT(MVSIZ),IAD2
77 my_real REDUC, UPWL(NV46,MVSIZ),FLUX1(MVSIZ), FLUX2(MVSIZ), FLUX3(MVSIZ), FLUX4(MVSIZ),FLUX5(MVSIZ), FLUX6(MVSIZ)
78 my_real, TARGET :: qmvl(7:12,mvsiz)
79 my_real, DIMENSION(:) ,POINTER :: pqmv7,pqmv8,pqmv9,pqmv10,pqmv11,pqmv12
80 LOGICAL :: debug_outp
81C-----------------------------------------------
82C S o u r c e L i n e s
83C-----------------------------------------------
84 DO i=lft,llt
85 flux1(i) = flux(1,i)
86 flux2(i) = flux(2,i)
87 flux3(i) = flux(3,i)
88 flux4(i) = flux(4,i)
89 ENDDO
90 IF(nv46==6)THEN
91 DO i=lft,llt
92 flux5(i) = flux(5,i)
93 flux6(i) = flux(6,i)
94 ENDDO
95 ENDIF
96C-----------------------------------------------
97 DO i=lft,llt
98 ii=i+nft
99 mat(i)=ixs(1,ii)
100 ENDDO
101C-----------------------------------------------
102C UPWIND
103C-----------------------------------------------
104 DO j=1,nv46
105 DO i=lft,llt
106 upwl(j,i)=pm(16,mat(i))
107 ENDDO
108 ENDDO
109 !======================================================!
110 ! BOUNDARY FACE : no volume flux by default !
111 ! slip wall bc !
112 !======================================================!
113 DO i=lft,llt
114 iad2 = ale_connect%ee_connect%iad_connect(i + nft)
115 reduc=pm(92,mat(i))
116 !---face1---!
117 ii=ale_connect%ee_connect%connected(iad2 + 1 - 1)
118 IF(ii == 0)THEN
119 flux1(i)=flux1(i)*reduc
120 ENDIF
121 !---face2---!
122 ii=ale_connect%ee_connect%connected(iad2 + 2 - 1)
123 IF(ii == 0)THEN
124 flux2(i)=flux2(i)*reduc
125 ENDIF
126 !---face3---!
127 ii=ale_connect%ee_connect%connected(iad2 + 3 - 1)
128 IF(ii == 0)THEN
129 flux3(i)=flux3(i)*reduc
130 ENDIF
131 !---face4---!
132 ii=ale_connect%ee_connect%connected(iad2 + 4 - 1)
133 IF(ii == 0)THEN
134 flux4(i)=flux4(i)*reduc
135 ENDIF
136 IF(nv46==6)THEN
137 !---face5---!
138 ii=ale_connect%ee_connect%connected(iad2 + 5 - 1)
139 IF(ii == 0)THEN
140 flux5(i)=flux5(i)*reduc
141 ENDIF
142 !---face6---!
143 ii=ale_connect%ee_connect%connected(iad2 + 6 - 1)
144 IF(ii == 0)THEN
145 flux6(i)=flux6(i)*reduc
146 ENDIF
147 ENDIF
148 ! IF(NV46==6)THEN
149 ENDDO
150
151 DO i=lft,llt
152 flux(1,i) = flux1(i)-upwl(1,i)*abs(flux1(i))
153 flux(2,i) = flux2(i)-upwl(2,i)*abs(flux2(i))
154 flux(3,i) = flux3(i)-upwl(3,i)*abs(flux3(i))
155 flux(4,i) = flux4(i)-upwl(4,i)*abs(flux4(i))
156 ENDDO
157 IF(nv46==6)THEN
158 DO i=lft,llt
159 flux(5,i) = flux5(i)-upwl(5,i)*abs(flux5(i))
160 flux(6,i) = flux6(i)-upwl(6,i)*abs(flux6(i))
161 ENDDO
162 ENDIF
163
164 IF(iflg == 1)THEN
165 pqmv7 => qmv(07,lft:llt)
166 pqmv8 => qmv(08,lft:llt)
167 pqmv9 => qmv(09,lft:llt)
168 pqmv10 => qmv(10,lft:llt)
169 pqmv11 => qmv(11,lft:llt)
170 pqmv12 => qmv(12,lft:llt)
171 ELSE
172 pqmv7 => qmvl(07,lft:llt)
173 pqmv8 => qmvl(08,lft:llt)
174 pqmv9 => qmvl(09,lft:llt)
175 pqmv10 => qmvl(10,lft:llt)
176 pqmv11 => qmvl(11,lft:llt)
177 pqmv12 => qmvl(12,lft:llt)
178 ENDIF
179 !MOMENTUM ADVECTION NEEDED FOR AMOMT3.F
180
181 IF(nv46==6)THEN
182 DO i=lft,llt
183 pqmv7(i) = flux1(i)+upwl(1,i)*abs(flux1(i))
184 pqmv8(i) = flux2(i)+upwl(2,i)*abs(flux2(i))
185 pqmv9(i) = flux3(i)+upwl(3,i)*abs(flux3(i))
186 pqmv10(i) = flux4(i)+upwl(4,i)*abs(flux4(i))
187 pqmv11(i) = flux5(i)+upwl(5,i)*abs(flux5(i))
188 pqmv12(i) = flux6(i)+upwl(6,i)*abs(flux6(i))
189 flu1(i) = pqmv7(i) + pqmv8(i) + pqmv9(i) + pqmv10(i) +pqmv11(i) + pqmv12(i)
190 ENDDO
191 IF(itrimat > 0)THEN
192 DO i=lft,llt
193 ddvol(i)=half*( flu1(i)+ flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i)+flux(5,i)+flux(6,i) )
194 ENDDO
195 ENDIF
196 ELSE
197 DO i=lft,llt
198 pqmv7(i) = flux1(i)+upwl(1,i)*abs(flux1(i))
199 pqmv8(i) = flux2(i)+upwl(2,i)*abs(flux2(i))
200 pqmv9(i) = flux3(i)+upwl(3,i)*abs(flux3(i))
201 pqmv10(i) = flux4(i)+upwl(4,i)*abs(flux4(i))
202 flu1(i) = pqmv7(i) + pqmv8(i) + pqmv9(i) + pqmv10(i)
203 ENDDO
204 IF(itrimat > 0)THEN
205 DO i=lft,llt
206 ddvol(i)=half*( flu1(i)+ flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i) )
207 ENDDO
208 ENDIF
209 ENDIF
210
211 !INTERFACE 22 ONLY - OUTPUT---------------! OBSOLETE
212 debug_outp = .false.
213 if(ibug22_upwind /= 0)then
214 if(ibug22_upwind > 0)then
215 do i=lft,llt
216 if(ixs(11,i+nft) == ibug22_upwind)debug_outp=.true.
217 enddo
218 elseif(ibug22_upwind==-1)then
219 debug_outp = .true.
220 endif
221 if(((itrimat /= ibug22_itrimat).and.(ibug22_itrimat /= -1)))debug_outp=.false.
222 endif
223!#!include "lockon.inc"
224 if(debug_outp)then
225 print *, " |----ale51_upwind3.F-----|"
226 print *, " | THREAD INFORMATION |"
227 print *, " |------------------------|"
228 print *, " NCYCLE=", ncycle
229 do i=lft,llt
230 if(ibug22_upwind/=ixs(11,i+nft).and.ibug22_upwind/=-1)cycle
231 !if (tag22(i)==zero)print *," UNCUT"!cycle
232 print *, " brique =", ixs(11,nft+i)
233 print *, " itrimat =", itrimat
234 write (*,fmt='(A,6E26.14)')" Flux(1:6) =", flux(1:6,i)
235 write (*,fmt='(A,1E26.14)')" Flu1 =", flu1(i)
236 print *, " ------------------------"
237 enddo
238!#!include "lockoff.inc"
239 endif
240 !-----------------------------------------!
241
242C-----------------------------------------------
243 RETURN
244 END
245C
ale51_upwind3
subroutine ale51_upwind3(pm, ixs, flux, flu1, ale_connect, itrimat, ddvol, qmv, iflg, nv46)
Definition ale51_upwind3.F:37
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
i22tri_mod
Definition cut-cell-search_mod.F:212
i22tri_mod::ibug22_upwind
integer ibug22_upwind
Definition cut-cell-search_mod.F:403
i22tri_mod::ibug22_itrimat
integer ibug22_itrimat
Definition cut-cell-search_mod.F:393
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:256