OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
nodal_schlieren.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!|| nodal_schlieren ../engine/source/output/anim/generate/nodal_schlieren.F
25!||--- called by ------------------------------------------------------
26!|| genani ../engine/source/output/anim/generate/genani.F
27!|| h3d_nodal_scalar ../engine/source/output/h3d/h3d_results/h3d_nodal_scalar.F
28!||--- calls -----------------------------------------------------
29!|| agrad3 ../engine/source/ale/ale3d/agrad3.F
30!||--- uses -----------------------------------------------------
31!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
32!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
33!|| element_mod ../common_source/modules/elements/element_mod.F90
34!|| i22bufbric_mod ../common_source/modules/interfaces/cut-cell-search_mod.f
35!|| i22edge_mod ../common_source/modules/interfaces/cut-cell-buffer_mod.F
36!|| i22tri_mod ../common_source/modules/interfaces/cut-cell-search_mod.F
37!|| initbuf_mod ../engine/share/resol/initbuf.F
38!||====================================================================
39 SUBROUTINE nodal_schlieren(WA4, X,IXS, IXQ, ITAB,IPARG ,IBID,ELBUF_TAB,ALE_CONNECTIVITY)
40C-----------------------------------------------
41C D e s c r i p t i o n
42C-----------------------------------------------
43C This subroutine outputs data for schlieren.
44C schlieren is eta = exp (-C ||grad(rho)||)
45C C is a constant which help user to adjust "brightness"
46C RADIOSS outputs density gradient which is recuired to output schlieren.
47C 'C' cosntant must be tuned during post-treatment then
48C is it introduced with HV result math.
49C-----------------------------------------------
50C P r e - C o n d i t i o n s
51C-----------------------------------------------
52C IALEL > 0
53C where IALEL =IPARG(7,NG)+IPARG(11,NG)
54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE initbuf_mod
58 USE elbufdef_mod
59 USE i22bufbric_mod
60 USE i22edge_mod
61 USE i22tri_mod
62 USE ale_connectivity_mod
63 use element_mod , only : nixs,nixq
64C-----------------------------------------------
65C I m p l i c i t T y p e s
66C-----------------------------------------------
67#include "implicit_f.inc"
68C-----------------------------------------------
69C C o m m o n B l o c k s
70C-----------------------------------------------
71#include "com01_c.inc"
72#include "com04_c.inc"
73#include "vect01_c.inc"
74#include "param_c.inc"
75#include "mvsiz_p.inc"
76C-----------------------------------------------
77C D u m m y A r g u m e n t s
78C-----------------------------------------------
79 INTEGER IXQ(NIXQ,*),IXS(NIXS,*),ITAB(*),IPARG(NPARG,*)
80 REAL WA4(*)
81 my_real :: x(3,*)
82 INTEGER :: IBID
83 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
84 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECTIVITY
85C-----------------------------------------------
86C L o c a l V a r i a b l e s
87C-----------------------------------------------
88 INTEGER I, ITYP, NEL,
89 . IR1, IR2, J, JJ, NNO_L, NNI_L, II3, II4, JJJ, NNI,
90 . IALEL,NNOD,IPOS,NGv,IDLOCv,K, IAD2
91 INTEGER IV(6), IE
92 INTEGER NG
93 TYPE(g_bufel_) ,POINTER :: GBUF
94 my_real, ALLOCATABLE, DIMENSION(:) :: count_vol
95 my_real d,v, dphi(mvsiz)
96
97 my_real :: grad(6,mvsiz)
98C-----------------------------------------------
99C D e s c r i p t i o n
100C-----------------------------------------------
101C This subroutine write nodal shadowgraph (schlieren)
102C-----------------------------------------------
103C S o u r c e L i n e s
104C-----------------------------------------------
105
106 wa4(1:numnod) = zero
107
108 RETURN
109
110 !TODO : 2D,3D, QUAD, HEXA,TETRA
111
112
113
114
115
116 nnod = nixs-3 !8-node brick or 4-node quad
117
118 !---------------------------------------------------------!
119 ! ALE STANDARD FORMULATION : 3D !
120 !---------------------------------------------------------!
121 ALLOCATE(count_vol(numnod))
122 count_vol(:) = 0
123 DO ng = 1, ngroup
124 nel =iparg(2,ng)
125 nft =iparg(3,ng)
126 ityp =iparg(5,ng)
127 ialel =iparg(7,ng)+iparg(11,ng)
128 IF(ityp/=1 .AND. ityp/=2)cycle
129 IF(ialel==0)cycle
130 gbuf => elbuf_tab(ng)%GBUF
131 DO i=1,nel
132 j = i+nft
133! PHI(J) = GBUF%RHO(I)
134 ENDDO
135 ENDDO
136
137 DO ng = 1, ngroup
138 nel =iparg(2,ng)
139 nft =iparg(3,ng)
140 ityp =iparg(5,ng)
141 ialel =iparg(7,ng)+iparg(11,ng)
142 IF(ityp/=1 .AND. ityp/=2)cycle
143 IF(ialel==0)cycle
144 gbuf => elbuf_tab(ng)%GBUF
145 DO i=1,nel
146 lft = 1
147 llt = nel
148 CALL agrad3(
149 1 ixs, x, ale_connectivity,grad)
150 ie =nft+i
151 iad2 = ale_connectivity%ee_connect%iad_connect(ie)
152 iv(1)=ale_connectivity%ee_connect%connected(iad2 + 1 - 1)
153 iv(2)=ale_connectivity%ee_connect%connected(iad2 + 2 - 1)
154 iv(3)=ale_connectivity%ee_connect%connected(iad2 + 3 - 1)
155 iv(4)=ale_connectivity%ee_connect%connected(iad2 + 4 - 1)
156 iv(5)=ale_connectivity%ee_connect%connected(iad2 + 5 - 1)
157 iv(6)=ale_connectivity%ee_connect%connected(iad2 + 6 - 1)
158 IF(iv(1)<=0)iv(1)=ie
159 IF(iv(2)<=0)iv(2)=ie
160 IF(iv(3)<=0)iv(3)=ie
161 IF(iv(4)<=0)iv(4)=ie
162 IF(iv(5)<=0)iv(5)=ie
163 IF(iv(6)<=0)iv(6)=ie
164 dphi(i) = zero
165! . (PHI(IV(1))-PHI(IE))*GRAD(1,I)
166! . +(PHI(IV(2))-PHI(IE))*GRAD(2,I)
167! . +(PHI(IV(3))-PHI(IE))*GRAD(3,I)
168! . +(PHI(IV(4))-PHI(IE))*GRAD(4,I)
169! . +(PHI(IV(5))-PHI(IE))*GRAD(5,I)
170! . +(PHI(IV(6))-PHI(IE))*GRAD(6,I)
171 DO j=2,nnod+1
172 jj=ixs(j,nft+i)
173 k = j-1
174 wa4(jj) = wa4(jj)+ dphi(i)
175 count_vol(jj) = count_vol(jj) + 1
176 ENDDO
177 ENDDO
178 enddo!next NG
179
180 !applying weight factor
181 DO i=1,numnod
182 IF(count_vol(i)/=zero)THEN
183 wa4(i)=wa4(i)/count_vol(i)
184 ENDIF
185 ENDDO
186 DEALLOCATE(count_vol)
187
188 END
189
my_real
#define my_real
Definition cppsort.cpp:32
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
i22bufbric_mod
Definition cut-cell-search_mod.F:95
i22edge_mod
Definition cut-cell-buffer_mod.F:48
i22tri_mod
Definition cut-cell-search_mod.F:212
initbuf_mod
Definition initbuf.F:160
nodal_schlieren
subroutine nodal_schlieren(wa4, x, ixs, ixq, itab, iparg, ibid, elbuf_tab, ale_connectivity)
Definition nodal_schlieren.F:40
agrad3
subroutine agrad3(ixs, x, ale_connectivity, grad, nel)
Definition agrad3.F:30
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:256