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

Go to the source code of this file.

Functions/Subroutines

subroutine alefvm_gravity (voln, ixs, rho)

Function/Subroutine Documentation

◆ alefvm_gravity()

subroutine alefvm_gravity ( voln,
integer, dimension(nixs,*) ixs,
rho )

Definition at line 31 of file alefvm_gravity.F.

32C-----------------------------------------------
33C D e s c r i p t i o n
34C-----------------------------------------------
35C 'alefvm' is related to a collocated scheme (built from FVM and based on Godunov scheme)
36C which was temporarily introduced for experimental option /INTER/TYPE22 (FSI coupling with cut cell method)
37C This cut cell method is not completed, abandoned, and is not an official option.
38C There is no other use for this scheme which is automatically enabled when /INTER/TYPE22 is defined (INT22>0 => IALEFVM=1).
39C
40C This subroutine is treating an uncut cell.
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE alefvm_mod
45 use element_mod , only : nixs
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C G l o b a l P a r a m e t e r s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "vect01_c.inc"
58#include "com01_c.inc"
59#include "param_c.inc"
60C-----------------------------------------------
61C D e s c r i p t i o n
62C-----------------------------------------------
63C This subroutines computes gravity forces for
64C finite volume scheme (IALEFVM==1)
65C-----------------------------------------------
66C D u m m y A r g u m e n t s
67C-----------------------------------------------
68 INTEGER :: IXS(NIXS,*)
69 my_real :: voln(mvsiz), rho(mvsiz)
70C-----------------------------------------------
71C L o c a l V a r i a b l e s
72C-----------------------------------------------
73 INTEGER :: I, J
74 INTEGER :: NVERTEX,INOd,II
75 my_real :: m_cell(mvsiz),accel(3,mvsiz)
76 LOGICAL :: debug_outp
77 INTEGER :: idbf,idbl
78C-----------------------------------------------
79C P r e - C o n d i t i o n s
80C-----------------------------------------------
81! IF(IALEFVM==0)RETURN
82! IF(NGRAV==0) RETURN !no there is Call_my barrier
83C-----------------------------------------------
84C S o u r c e L i n e s
85C-----------------------------------------------
86
87
88 !-------------------------------------------------------------!
89 ! GRAVITY COMPUTED FOR CELL CENTROIDS !
90 !-------------------------------------------------------------!
91 DO i=lft,llt
92 nvertex = 0
93 accel(1:3,i) = zero
94 !loop on vertexes
95 DO j=1,8
96 inod = ixs(1+j,i+nft)
97 IF(alefvm_buffer%VERTEX(4,inod)==zero) cycle
98 nvertex = nvertex + 1
99 accel(1,i) = accel(1,i) + alefvm_buffer%VERTEX(1,inod)
100 accel(2,i) = accel(2,i) + alefvm_buffer%VERTEX(2,inod)
101 accel(3,i) = accel(3,i) + alefvm_buffer%VERTEX(3,inod)
102 ENDDO
103 IF(nvertex>0)THEN
104 accel(1,i) = accel(1,i) / nvertex
105 accel(2,i) = accel(2,i) / nvertex
106 accel(3,i) = accel(3,i) / nvertex
107 ENDIF
108 enddo!next I
109
110 DO i=lft,llt
111 m_cell(i) = rho(i)*voln(i)
112 enddo!next I
113
114 DO i=lft,llt
115 ii = i + nft
116 alefvm_buffer%FCELL(1,ii) = alefvm_buffer%FCELL(1,ii) + m_cell(i) * accel(1,i)
117 alefvm_buffer%FCELL(2,ii) = alefvm_buffer%FCELL(2,ii) + m_cell(i) * accel(2,i)
118 alefvm_buffer%FCELL(3,ii) = alefvm_buffer%FCELL(3,ii) + m_cell(i) * accel(3,i)
119 enddo!next I
120
121 !DEBUG-OUTPUT---------------!
122 if(alefvm_param%IOUTP_GRAV /= 0)then
123 debug_outp = .false.
124 if(alefvm_param%IOUTP_GRAV>0)then
125 do i=lft,llt
126 ii = nft + i
127 if(ixs(11,ii)==alefvm_param%IOUTP_GRAV)THEN
128 debug_outp = .true.
129 idbf = i
130 idbl = i
131 EXIT
132 endif
133 enddo
134 elseif(alefvm_param%IOUTP_GRAV==-1)then
135 debug_outp=.true.
136 idbf = lft
137 idbl = llt
138 endif
139!#!include "lockon.inc"
140 if(debug_outp)then
141!#!include "lockon.inc"
142 print *, " |----alefvm_gravity.F----|"
143 print *, " | THREAD INFORMATION |"
144 print *, " |------------------------|"
145 print *, " NCYCLE =", ncycle
146 do i=idbf,idbl
147 ii = nft + i
148 print *, " brique=", ixs(11,nft+i)
149 write(*,fmt='(A24,1A26)') " ",
150 . "#--------- cell----------#"
151 write (*,fmt='(A,1E26.14)') " Rho =", rho(i)
152 write (*,fmt='(A,1E26.14)') " Vol =", voln(i)
153 write (*,fmt='(A,1E26.14)') " Mass =", m_cell(i)
154 write (*,fmt='(A,1E26.14)') " Accel-X =", accel(1,i)
155 write (*,fmt='(A,1E26.14)') " Accel-Y =", accel(2,i)
156 write (*,fmt='(A,1E26.14)') " Accel-Z =", accel(3,i)
157 write(*,fmt='(A24,8A26)') " ",
158 . "#--------- nod_1 ---------","#--------- nod_2 ---------",
159 . "#--------- nod_3 ---------","#--------- nod_4 ---------",
160 . "#--------- nod_5 ---------","#--------- nod_6 ---------",
161 . "#--------- nod_7 ---------","#--------- nod_8 --------#"
162 write (*,fmt='(A,8E26.14)') " acc-X =", alefvm_buffer%VERTEX(1,ixs(2:9,i))
163 write (*,fmt='(A,8E26.14)') " acc-Y =", alefvm_buffer%VERTEX(2,ixs(2:9,i))
164 write (*,fmt='(A,8E26.14)') " acc-Z =", alefvm_buffer%VERTEX(3,ixs(2:9,i))
165 print *, " "
166 enddo
167!#!include "lockoff.inc"
168 endif
169 endif
170 !-----------------------------------------!
171
172 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
alefvm_mod
Definition alefvm_mod.F:88
alefvm_mod::alefvm_buffer
type(alefvm_buffer_), target alefvm_buffer
Definition alefvm_mod.F:120
alefvm_mod::alefvm_param
type(alefvm_param_), target alefvm_param
Definition alefvm_mod.F:121