OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
snorm3t.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!|| snorm3t ../engine/source/multifluid/snorm3t.F
25!||--- called by ------------------------------------------------------
26!|| multi_face_elem_data ../engine/source/multifluid/multi_face_data_elem.F
27!|| multi_fvm2fem ../engine/source/multifluid/multi_fvm2fem.f
28!||--- calls -----------------------------------------------------
29!|| arret ../engine/source/system/arret.F
30!||--- uses -----------------------------------------------------
31!|| element_mod ../common_source/modules/elements/element_mod.F90
32!||====================================================================
33 SUBROUTINE snorm3t (NEL, NFT, JALE, IXS, XGRID, WGRID,
34 . NORM, WFAC, SURF)
35 use element_mod , only : nixs
36C-----------------------------------------------
37C D e s c r i p t i o n
38C-----------------------------------------------
39C Computes normal vector to the faces of each element in a group
40C for a 3d solid element (hence Snorm3)
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C C o m m o n B l o c k s
47C-----------------------------------------------
48! NIXS
49C-----------------------------------------------
50C D u m m y A r g u m e n t s
51C-----------------------------------------------
52 INTEGER, INTENT(IN) :: NEL, NFT, JALE, IXS(NIXS, *)
53 my_real, INTENT(IN) ::
54 . xgrid(3, *), wgrid(*)
55 my_real, INTENT(OUT) :: wfac(3, 6, nel), surf(6, nel)
56 my_real, INTENT(OUT), TARGET :: norm(3, 6, nel)
57C-----------------------------------------------
58C L o c a l V a r i a b l e s
59C-----------------------------------------------
60 INTEGER :: II, NODE1, NODE2, NODE3, NODE4, KFACE
61 my_real
62 . x1(3), x2(3), x3(3), x4(3),
63 . w1(3), w2(3), w3(3), w4(3), xc(3), xf(3)
64 my_real, POINTER :: nx, ny, nz
65
66 DO ii = 1, nel
67C Nodes of the element locally stored in NODE* for memory access
68 node1 = ixs(2, ii + nft)
69 node2 = ixs(4, ii + nft)
70 node3 = ixs(7, ii + nft)
71 node4 = ixs(6, ii + nft)
72C Node coordinates
73 x1(1:3) = xgrid(1:3, node1)
74 x2(1:3) = xgrid(1:3, node2)
75 x3(1:3) = xgrid(1:3, node3)
76 x4(1:3) = xgrid(1:3, node4)
77
78 xc(1:3) = fourth * (x1(1:3) + x2(1:3) + x3(1:3) + x4(1:3))
79 IF (jale /= 0) THEN
80C Node grid velocities
81 w1(1:3) = wgrid(3 * (node1 - 1) + 1 : 3 * (node1 - 1) + 3)
82 w2(1:3) = wgrid(3 * (node2 - 1) + 1 : 3 * (node2 - 1) + 3)
83 w3(1:3) = wgrid(3 * (node3 - 1) + 1 : 3 * (node3 - 1) + 3)
84 w4(1:3) = wgrid(3 * (node4 - 1) + 1 : 3 * (node4 - 1) + 3)
85 ELSE
86! Euler
87 w1(1:3) = zero
88 w2(1:3) = zero
89 w3(1:3) = zero
90 w4(1:3) = zero
91 ENDIF
92C Face normal
93C Face 1 -> nodes 1 2 3
94 kface = 5
95 norm(1, kface, ii) = half * ((x3(2) - x1(2)) * (x2(3) - x1(3)) -
96 . (x3(3) - x1(3)) * (x2(2) - x1(2)))
97 norm(2, kface, ii) = half * ((x3(3) - x1(3)) * (x2(1) - x1(1)) -
98 . (x3(1) - x1(1)) * (x2(3) - x1(3)))
99 norm(3, kface, ii) = half * ((x3(1) - x1(1)) * (x2(2) - x1(2)) -
100 . (x3(2) - x1(2)) * (x2(1) - x1(1)))
101 nx => norm(1, kface, ii)
102 ny => norm(2, kface, ii)
103 nz => norm(3, kface, ii)
104 surf(kface, ii) = sqrt(nx * nx + ny * ny + nz * nz)
105 nx = -nx / surf(kface, ii)
106 ny = -ny / surf(kface, ii)
107 nz = -nz / surf(kface, ii)
108 xf(1:3) = third * (x1 + x2 + x3)
109 IF ((xf(1) - xc(1)) * nx + (xf(2) - xc(2)) * ny + (xf(3) - xc(3)) * nz < zero) THEN
110 CALL arret(2)
111 ENDIF
112C Face 2
113 kface = 6
114 norm(1, kface, ii) = half * ((x1(2) - x4(2)) * (x2(3) - x4(3)) -
115 . (x1(3) - x4(3)) * (x2(2) - x4(2)))
116 norm(2, kface, ii) = half * ((x1(3) - x4(3)) * (x2(1) - x4(1)) -
117 . (x1(1) - x4(1)) * (x2(3) - x4(3)))
118 norm(3, kface, ii) = half * ((x1(1) - x4(1)) * (x2(2) - x4(2)) -
119 . (x1(2) - x4(2)) * (x2(1) - x4(1)))
120 nx => norm(1, kface, ii)
121 ny => norm(2, kface, ii)
122 nz => norm(3, kface, ii)
123 surf(kface, ii) = sqrt(nx * nx + ny * ny + nz * nz)
124 nx = -nx / surf(kface, ii)
125 ny = -ny / surf(kface, ii)
126 nz = -nz / surf(kface, ii)
127 xf(1:3) = third * (x1 + x2 + x4)
128 IF ((xf(1) - xc(1)) * nx + (xf(2) - xc(2)) * ny + (xf(3) - xc(3)) * nz < zero) THEN
129 CALL arret(2)
130 ENDIF
131C Face 3
132 kface = 2
133 norm(1, kface, ii) = half * ((x2(2) - x4(2)) * (x3(3) - x4(3)) -
134 . (x2(3) - x4(3)) * (x3(2) - x4(2)))
135 norm(2, kface, ii) = half * ((x2(3) - x4(3)) * (x3(1) - x4(1)) -
136 . (x2(1) - x4(1)) * (x3(3) - x4(3)))
137 norm(3, kface, ii) = half * ((x2(1) - x4(1)) * (x3(2) - x4(2)) -
138 . (x2(2) - x4(2)) * (x3(1) - x4(1)))
139 nx => norm(1, kface, ii)
140 ny => norm(2, kface, ii)
141 nz => norm(3, kface, ii)
142 surf(kface, ii) = sqrt(nx * nx + ny * ny + nz * nz)
143 nx = -nx / surf(kface, ii)
144 ny = -ny / surf(kface, ii)
145 nz = -nz / surf(kface, ii)
146 xf(1:3) = third * (x4 + x2 + x3)
147 IF ((xf(1) - xc(1)) * nx + (xf(2) - xc(2)) * ny + (xf(3) - xc(3)) * nz < zero) THEN
148 CALL arret(2)
149 ENDIF
150C Face 4
151 kface = 4
152 norm(1, kface, ii) = half * ((x3(2) - x4(2)) * (x1(3) - x4(3)) -
153 . (x3(3) - x4(3)) * (x1(2) - x4(2)))
154 norm(2, kface, ii) = half * ((x3(3) - x4(3)) * (x1(1) - x4(1)) -
155 . (x3(1) - x4(1)) * (x1(3) - x4(3)))
156 norm(3, kface, ii) = half * ((x3(1) - x4(1)) * (x1(2) - x4(2)) -
157 . (x3(2) - x4(2)) * (x1(1) - x4(1)))
158 nx => norm(1, kface, ii)
159 ny => norm(2, kface, ii)
160 nz => norm(3, kface, ii)
161 surf(kface, ii) = sqrt(nx * nx + ny * ny + nz * nz)
162 nx = -nx / surf(kface, ii)
163 ny = -ny / surf(kface, ii)
164 nz = -nz / surf(kface, ii)
165 xf(1:3) = third * (x1 + x4 + x3)
166 IF ((xf(1) - xc(1)) * nx + (xf(2) - xc(2)) * ny + (xf(3) - xc(3)) * nz < zero) THEN
167 CALL arret(2)
168 ENDIF
169C Face grid velocity 1
170 wfac(1:3, 5, ii) = third * (w1(1:3) + w2(1:3) + w3(1:3))
171C Face grid velocity 2
172 wfac(1:3, 6, ii) = third * (w1(1:3) + w4(1:3) + w2(1:3))
173C Face grid velocity 3
174 wfac(1:3, 2, ii) = third * (w2(1:3) + w4(1:3) + w3(1:3))
175C Face grid velocity 4
176 wfac(1:3, 4, ii) = third * (w1(1:3) + w3(1:3) + w4(1:3))
177C Nullify quantities for KFACE = 1, 3
178 wfac(1:3, 1, ii) = zero
179 wfac(1:3, 3, ii) = zero
180 surf(1, ii) = zero
181 surf(3, ii) = zero
182 norm(1:3, 1, ii) = zero
183 norm(1:3, 3, ii) = zero
184 ENDDO
185 END SUBROUTINE snorm3t
my_real
#define my_real
Definition cppsort.cpp:32
norm
norm(diag(diag(diag(inv(mat))) -id.SOL), 2) % destroy mumps instance id.JOB
multi_fvm2fem
subroutine multi_fvm2fem(timestep, elbuf_tab, itask, ixs, ixq, iparg, xgrid, accele, vel, wgrid, ms, msnf, veul, stifn, fsky, iads, fskym, condn, condnsky, multi_fvm, nodadt_therm)
Definition multi_fvm2fem.F:49
snorm3t
subroutine snorm3t(nel, nft, jale, ixs, xgrid, wgrid, norm, wfac, surf)
Definition snorm3t.F:35
arret
subroutine arret(nn)
Definition arret.F:86