OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
ale51_antidiff2.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_antidiff2 ../engine/source/ale/ale51/ale51_antidiff2.F
25!||--- called by ------------------------------------------------------
26!|| afluxt ../engine/source/ale/ale51/afluxt.F
27!||--- uses -----------------------------------------------------
28!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
29!|| ale_mod ../common_source/modules/ale/ale_mod.F
30!|| element_mod ../common_source/modules/elements/element_mod.F90
31!|| segvar_mod ../engine/share/modules/segvar_mod.F
32!||====================================================================
33 SUBROUTINE ale51_antidiff2(FLUX,ALE_CONNECT,ALPH,VOL,IXQ,FLUX_VOIS,
34 . N4_VOIS,ITAB,ITRIMAT,
35 . SEGVAR)
36C-----------------------------------------------
37C D e c r i p t i o n
38C-----------------------------------------------
39C
40C ** Warning, this formulation (for volume fractions) using \eta_3 parameter from /UPWIND option is not conservative.
41C It was marked as OBSOLETE. Forget it and use MUSCL instead
42C
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE segvar_mod
47 USE ale_connectivity_mod
48 USE ale_mod
49 use element_mod , only : nixq
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "spmd_c.inc"
58#include "com04_c.inc"
59#include "com08_c.inc"
60#include "vect01_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER IXQ(NIXQ,NUMELQ),N4_VOIS(NUMELQ+NQVOIS,4),ITAB(NUMNOD)
65 INTEGER,INTENT(IN) :: ITRIMAT
66 my_real flux(4,*),alph(*),vol(*),flux_vois(numelq+nqvois,6)
67 TYPE(t_segvar),INTENT(IN) :: SEGVAR
68 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
69C-----------------------------------------------
70C L o c a l V a r i a b l e s
71C-----------------------------------------------
72 INTEGER I,II,K,JV(6),KV(6),IAD2,IAD3
73 my_real
74 . vol0,av0,uav0,alphi,ualphi,aaa,ff(6),udt,phi0
75C-----------------------------------------------
76C S o u r c e L i n e s
77C-----------------------------------------------
78 jv = -huge(jv(1))
79 kv = -huge(kv(1))
80 ff = -huge(ff(1))
81C-----------------------------------------------
82C WET SURFACE
83C-----------------------------------------------
84 IF(dt1 > zero)THEN
85 udt = one/dt1
86 ELSE
87 udt = zero
88 ENDIF
89C
90 DO i=lft,llt
91 ii = i + nft
92 iad2 = ale_connect%ee_connect%iad_connect(ii)
93 vol0 = vol(i)*udt
94 av0 = alph(ii) * vol0
95 uav0 = vol0 - av0
96 alphi = zero
97 ualphi = zero
98 phi0 = zero
99C-----------------------------------------------
100C adjacent face & total outgoing volume flux
101C-----------------------------------------------
102 DO k=1,4
103 IF(flux(k,ii) > zero)THEN
104 jv(k) = ale_connect%ee_connect%connected(iad2 + k - 1)
105 kv(k) = k
106 IF(jv(k) == 0)THEN
107 jv(k) = ii
108 ff(k) = alph(ii)*flux(k,ii)
109 ELSEIF(jv(k) < 0)THEN
110 !case <0 : ebcs (-IVOIS is then segment identifier)
111 ff(k) = segvar%PHASE_ALPHA(itrimat,-jv(k)) *flux(k,ii)
112 ELSEIF(jv(k) <= numelq)THEN
113 iad3 = ale_connect%ee_connect%iad_connect(jv(k))
114 IF(ale_connect%ee_connect%connected(iad3 + 1 - 1) == ii) kv(k) = 1
115 IF(ale_connect%ee_connect%connected(iad3 + 2 - 1) == ii) kv(k) = 2
116 IF(ale_connect%ee_connect%connected(iad3 + 3 - 1) == ii) kv(k) = 3
117 IF(ale_connect%ee_connect%connected(iad3 + 4 - 1) == ii) kv(k) = 4
118 ff(k) = alph(jv(k))*flux(k,ii)
119 ELSE
120 !spmd
121 ff(k) = alph(jv(k))*flux(k,ii)
122 ENDIF
123c outgoing volume flux (estimation)
124 alphi = alphi + ff(k)
125c initial outgoing flux
126 phi0 = phi0 + flux(k,ii)
127 ENDIF
128 ENDDO
129c outgoing rest
130 ualphi = phi0 - alphi
131C-----------------------------------------------
132C outgoing volume flux by face
133C-----------------------------------------------
134 IF(alphi > av0.AND.av0 > zero)THEN
135C-----------------------------------------------
136C outgoing volume flux > volume
137C-----------------------------------------------
138 aaa = av0 / alphi
139 DO k=1,4
140 IF(flux(k,ii) > zero)THEN
141 ff(k) = ff(k) * aaa
142 ENDIF
143 ENDDO
144 ELSEIF(ualphi > uav0.AND.uav0 > zero)THEN
145
146C-----------------------------------------------
147C outgoing rest > available volume
148C-----------------------------------------------
149 aaa = uav0/ualphi
150 DO k=1,4
151 IF(flux(k,ii) > zero)THEN
152 ff(k) = flux(k,ii) + (ff(k)-flux(k,ii))*aaa
153 ENDIF
154 ENDDO
155c ELSE
156 ENDIF
157C-----------------------------------------------
158C outgoing flux
159C-----------------------------------------------
160 DO k=1,4
161 IF(flux(k,ii) > zero)THEN
162 ff(k) = 0.5* (ff(k)*(1.-ale%UPWIND%UPWSM)
163 . + alph(ii)*flux(k,ii)*(1.+ale%UPWIND%UPWSM))
164 flux(k,ii) = ff(k)
165 IF(jv(k) < 0)THEN
166 !nothing to do
167 ELSEIF(jv(k) <= numelq)THEN
168 flux(kv(k),jv(k)) = -flux(k,ii)
169 ELSE
170 flux_vois(ii,k) = flux(k,ii)
171 n4_vois(ii,1) = itab(ixq(2,ii))
172 n4_vois(ii,2) = itab(ixq(3,ii))
173 n4_vois(ii,3) = itab(ixq(4,ii))
174 n4_vois(ii,4) = itab(ixq(5,ii))
175 ENDIF
176 ENDIF
177 ENDDO
178C-----------------------------------------------
179C incoming volume flux from EBCS
180C-----------------------------------------------
181 IF(nsegflu > 0)THEN
182 iad2 = ale_connect%ee_connect%iad_connect(ii)
183 DO k=1,4
184 IF(flux(k,ii) < zero .AND. ale_connect%ee_connect%connected(iad2 + k - 1) < 0)THEN
185 flux(k,ii) = segvar%PHASE_ALPHA(itrimat,-ale_connect%ee_connect%connected(iad2 + k - 1))*flux(k,ii)
186 ENDIF
187 ENDDO
188 ENDIF
189C-------------
190 enddo! next I
191C-----------------------------------------------
192 RETURN
193 END
194C
ale51_antidiff2
subroutine ale51_antidiff2(flux, ale_connect, alph, vol, ixq, flux_vois, n4_vois, itab, itrimat, segvar)
Definition ale51_antidiff2.F:36
my_real
#define my_real
Definition cppsort.cpp:32
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
ale_mod
Definition ale_mod.F:150
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:253
segvar_mod
Definition segvar_mod.F:63
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:256
segvar_mod::t_segvar
Definition segvar_mod.F:71