alefvm__main_8F_source.html

Copyright>        OpenRadioss

Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.

Copyright>

Copyright>        This program is free software: you can redistribute it and/or modify

Copyright>        it under the terms of the GNU Affero General Public License as published by

Copyright>        the Free Software Foundation, either version 3 of the License, or

Copyright>        (at your option) any later version.

Copyright>

Copyright>        This program is distributed in the hope that it will be useful,

Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of

Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the

Copyright>        GNU Affero General Public License for more details.

Copyright>

Copyright>        You should have received a copy of the GNU Affero General Public License

Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.

Copyright>

Copyright>

Copyright>        Commercial Alternative: Altair Radioss Software

Copyright>

Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss

Copyright>        software under a commercial license.  Contact Altair to discuss further if the

Copyright>        commercial version may interest you: https://www.altair.com/radioss/.

!||====================================================================

!||    alefvm_main            ../engine/source/ale/alefvm/alefvm_main.F

!||--- called by ------------------------------------------------------

!||    resol                  ../engine/source/engine/resol.F

!||--- calls      -----------------------------------------------------

!||    alefvm_scheme          ../engine/source/ale/alefvm/alefvm_scheme.F

!||    alefvm_sfint3          ../engine/source/ale/alefvm/alefvm_sfint3.F

!||    alefvm_sfint3_int22    ../engine/source/ale/alefvm/alefvm_sfint3_int22.F

!||    initbuf                ../engine/share/resol/initbuf.F

!||    my_barrier             ../engine/source/system/machine.F

!||--- uses       -----------------------------------------------------

!||    ale_connectivity_mod   ../common_source/modules/ale/ale_connectivity_mod.F

!||    alefvm_mod             ../common_source/modules/ale/alefvm_mod.F

!||    elbufdef_mod           ../common_source/modules/mat_elem/elbufdef_mod.F90

!||    i22tri_mod             ../common_source/modules/interfaces/cut-cell-search_mod.F

!||    initbuf_mod            ../engine/share/resol/initbuf.F

!||    intbufdef_mod          ../common_source/modules/interfaces/intbufdef_mod.F90

!||    segvar_mod             ../engine/share/modules/segvar_mod.F

!||====================================================================


      SUBROUTINE alefvm_main(

     1         X          , V       ,

     2         ELBUF_TAB  , VR      ,

     3         ALE_CONNECT, IPARG   , IXS    ,

     4         NALE       ,

     5         ITASK      , NODFT   , NODLT  ,IPM, NV46 ,MSNF     )

C-----------------------------------------------

C   D e s c r i p t i o n

C-----------------------------------------------

C 'alefvm' is related to a collocated scheme (built from FVM and based on Godunov scheme)

C  which was temporarily introduced for experimental option /INTER/TYPE22 (FSI coupling with cut cell method)

C This cut cell method is not completed, abandoned, and is not an official option.

C There is no other use for this scheme which is automatically enabled when /INTER/TYPE22 is defined (INT22>0 => IALEFVM=1).

C

C This subroutine is treating an uncut cell.

C-----------------------------------------------

C   M o d u l e s

C-----------------------------------------------

      USE initbuf_mod

      USE elbufdef_mod

      USE intbufdef_mod

      USE alefvm_mod

      USE i22tri_mod

      USE segvar_mod

      USE ale_connectivity_mod

C-----------------------------------------------

C   I m p l i c i t   T y p e s

C-----------------------------------------------

#include      "implicit_f.inc"

#include      "comlock.inc"

C-----------------------------------------------

C   C o m m o n   B l o c k s

C-----------------------------------------------

#include      "com01_c.inc"

#include      "com04_c.inc"

#include      "param_c.inc"

#include      "vect01_c.inc"

#include      "scr17_c.inc"

#include      "task_c.inc"

#include      "inter22.inc"

C-----------------------------------------------

C   D u m m y   A r g u m e n t s

C-----------------------------------------------

      INTEGER IXS(NIXS,*),IPARG(NPARG,NGROUP),NALE(*),NODFT,ITASK,NV46,NODLT,IPM(NPROPMI,NUMMAT)


      my_real X(3,*),V(3,*),VR(3,*),MSNF(*)

      TYPE(ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB

      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER N, NG, NF1,

     .   isolnod,   nel,

     .   istra

      INTEGER IPLA

      INTEGER IALEFVM_FLG, IMAT,NSG,NVC

      my_real

     .   bid


      TYPE(g_bufel_) ,POINTER :: GBUF

      TYPE(l_bufel_)  ,POINTER  :: LBUF


      INTEGER :: NIN,NBF,NBL,tNB

C-----------------------------------------------

C   P r e - C o n d i t i o n

C-----------------------------------------------

      IF(alefvm_param%IEnabled==0)RETURN

C-----------------------------------------------

C   S o u r c e   L i n e s

C-----------------------------------------------


      !---------------------------------------------------------!

      ! 2.  INIT/PARAMETERS                                     !

      !---------------------------------------------------------!

      nin            = 1

      bid            = zero

      !---------------------------------------------------------!

      ! 3.  INIT/MULTITHREADING                                 !

      !---------------------------------------------------------!

      nbf = 1+itask*nb/nthread

      nbl = (itask+1)*nb/nthread

      nbl = min(nbl,nb)

      tnb = nbl-nbf+1


      CALL my_barrier


      IF(int22/=0)THEN

         CALL alefvm_sfint3_int22(ixs, nv46, itask, nbf, nbl, nin)

      ENDIF


      !---------------------------------------------

      ! ALEFVM : FINITE VOLUME FOR MOMENTUM

      !---------------------------------------------


       alefvm_buffer%VERTEX(1:4,1:numnod) = zero ! for alefvm_expand_momentum.F => vel on nodes


       IF(alefvm_param%Ienabled > 0) THEN

         CALL my_barrier


         DO ng=itask+1,ngroup,nthread

           CALL initbuf(

     1                   iparg    ,ng      ,

     2                   mtn     ,nel     ,nft     ,iad     ,ity     ,

     3                   npt     ,jale    ,ismstr  ,jeul    ,jtur    ,

     4                   jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,

     5                   nvaux   ,jpor    ,jcvt    ,jclose  ,ipla    ,

     6                   irep    ,iint    ,igtyp   ,israt   ,isrot   ,

     7                   icsen   ,isorth  ,isorthg ,ifailure,jsms

     .                  )

           gbuf => elbuf_tab(ng)%GBUF

           IF (iparg(8,ng) == 1) cycle

           IF (jlag == 1 .OR. ity>2) cycle

           nsg           = iparg(10,ng)

           nvc           = iparg(19,ng)

           isolnod       = iparg(28,ng)

           istra         = iparg(44,ng)

           jsph          = 0

           isph2sol      = 0

           ipartsph      = iparg(69,ng)

           lft           = 1

           llt           = nel

           nf1           = nft+1

           IF(ity == 1 .AND. isolnod == 4)THEN

             !CALL S4FORC3()


           ELSEIF(ity == 1 .AND. isolnod /= 4)THEN

             lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)

             imat        = ixs(1,nft+lft) !all elem in group does have same material id.

             ialefvm_flg = ipm(251,imat)

             IF (ialefvm_flg <= 1) cycle

             CALL alefvm_sfint3(

     1                          ixs     , nv46       ,  ale_connect   , ialefvm_flg,

     2                          ipm     , iparg      ,  ng      ,

     3                          x       , gbuf%TAG22 ,nel   )

             CALL alefvm_scheme(

     1                          ixs     , ialefvm_flg,

     2                          gbuf%MOM, gbuf%VOL,  gbuf%RHO,

     3                          ipm     , gbuf%TAG22 ,

     4                          lbuf%SSP,gbuf%SIG ,  nel     )

           ELSEIF (ity == 2 .AND. jmult == 0) THEN

             !CALL QFORC2()

           ELSEIF(ity == 2 .AND. jmult /= 0)THEN

             !CALL BFORC2()

           ENDIF

         enddo!next NG

       endif!IF (ILAW11 /= 0) THEN


       CALL my_barrier


       DO n=nodft,nodlt

         IF(alefvm_buffer%VERTEX(4,n) == zero .OR. nale(n)==0)cycle

         v(1,n)      = alefvm_buffer%VERTEX(1,n)

         v(2,n)      = alefvm_buffer%VERTEX(2,n)

         v(3,n)      = alefvm_buffer%VERTEX(3,n)

       ENDDO


#include "vectorize.inc"

        !---TRANSLATIONS---!

       DO n=nodft,nodlt

         IF(msnf(n)<=zero)cycle

         IF(nale(n)==0)cycle

         v(1,n)    = v(1,n) / msnf(n)

         v(2,n)    = v(2,n) / msnf(n)

         v(3,n)    = v(3,n) / msnf(n)

       ENDDO

       IF (iroddl/=0) THEN

#include "vectorize.inc"

         !---TROTATIONS---!

         DO n=nodft,nodlt

           IF(nale(n)==0)cycle

           vr(1,n) = zero

           vr(2,n) = zero

           vr(3,n) = zero

         ENDDO

       ENDIF


       CALL my_barrier

       alefvm_buffer%VERTEX(4,1:numnod) = zero ! was needed for alefvm_expand_momentum.F reset here for gravity on next cycle


      RETURN


      END

alefvm_main
subroutine alefvm_main(x, v, elbuf_tab, vr, ale_connect, iparg, ixs, nale, itask, nodft, nodlt, ipm, nv46, msnf)
Definition alefvm_main.F:48

alefvm_scheme
subroutine alefvm_scheme(ixs, ialefvm_flg, mom, vol, rho, ipm, iad22, ssp, sig, nel)
Definition alefvm_scheme.F:39

alefvm_sfint3
subroutine alefvm_sfint3(ixs, nv46, ale_connect, ialefvm_flg, ipm, iparg, ng, x, iad22, nel)
Definition alefvm_sfint3.F:37

alefvm_sfint3_int22
subroutine alefvm_sfint3_int22(ixs, nv46, itask, nbf, nbl, nin)
Definition alefvm_sfint3_int22.F:35

min
#define min(a, b)
Definition macros.h:20

ale_connectivity_mod
Definition ale_connectivity_mod.F:223

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

i22tri_mod
Definition cut-cell-search_mod.F:212

i22tri_mod::nb
integer nb
Definition cut-cell-search_mod.F:487

initbuf_mod
Definition initbuf.F:160

initbuf_mod::initbuf
subroutine initbuf(iparg, ng, mtn, llt, nft, iad, ity, npt, jale, ismstr, jeul, jtur, jthe, jlag, jmult, jhbe, jivf, mid, jpor, jcvt, jclose, jpla, irep, iint, igtyp, israt, isrot, icsen, isorth, isorthg, ifailure, jsms)
Definition initbuf.F:261

segvar_mod
Definition segvar_mod.F:61

my_barrier
subroutine my_barrier
Definition machine.F:31