cut-cell-buffer_mod.F

Go to the documentation of this file.

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/.
!||====================================================================
!||    i22edge_mod          ../common_source/modules/interfaces/cut-cell-buffer_mod.F
!||--- called by ------------------------------------------------------
!||    anim_nodal_p_elems   ../engine/source/output/anim/generate/anim_nodal_p_elems.F
!||    h3d_velvecc22        ../engine/source/output/h3d/h3d_results/h3d_velvecc22.F
!||    h3d_velvecz22        ../engine/source/output/h3d/h3d_results/h3d_velvecz22.F
!||    i22buce              ../engine/source/interfaces/intsort/i22buce.F
!||    i22bufbric_mod       ../common_source/modules/interfaces/cut-cell-search_mod.F
!||    i22cor3              ../engine/source/interfaces/int22/i22cor3.F
!||    i22get_prev_data     ../engine/source/interfaces/int22/i22get_prev_data.F
!||    i22ident             ../engine/source/interfaces/int22/i22ident.F
!||    i22intersect         ../engine/source/interfaces/int22/i22intersect.F
!||    i22subvol            ../engine/source/interfaces/int22/i22subvol.F
!||    i22tri_mod           ../common_source/modules/interfaces/cut-cell-search_mod.F
!||    i22trivox            ../engine/source/interfaces/intsort/i22trivox.F
!||    int22listcombi       ../engine/source/interfaces/int22/int22ListCombi.F
!||    nodal_schlieren      ../engine/source/output/anim/generate/nodal_schlieren.F
!||    nodald               ../engine/source/output/anim/generate/nodald.F
!||    nodalt               ../engine/source/output/anim/generate/nodalt.F
!||    nodalvfrac           ../engine/source/output/anim/generate/nodalvfrac.F
!||    nodalzvol            ../engine/source/output/anim/generate/nodalzvol.F
!||    resol                ../engine/source/engine/resol.F
!||    velvecc22            ../engine/source/output/anim/generate/velvec.F
!||    velvecz22            ../engine/source/output/anim/generate/velvecz22.F
!||====================================================================
      MODULE i22edge_mod
 
#include "my_real.inc"
 
        INTEGER,ALLOCATABLE,DIMENSION(:)       :: nbold
        INTEGER, ALLOCATABLE, DIMENSION(:,:)   :: list_b_old
        INTEGER,DIMENSION(:,:),ALLOCATABLE     :: iiad22
        INTEGER,DIMENSION(:)  ,ALLOCATABLE     :: ngb, idb, nelb, idlocb
        my_real,DIMENSION(:)  ,ALLOCATABLE     :: v22max_l, dx22min_l    !kinematic time step
 
        my_real,ALLOCATABLE,DIMENSION(:,:,:)   :: uvarl           ! multi-material merging
        my_real,ALLOCATABLE,DIMENSION(:,:)     :: eint_l          ! mono material only
        my_real,ALLOCATABLE,DIMENSION(:,:)     :: rho_l           ! mono material only
        my_real,ALLOCATABLE,DIMENSION(:,:,:)   :: mom_l           ! mono material only
        my_real,ALLOCATABLE,DIMENSION(:,:,:)   :: sig_l           ! mono material only
        my_real,ALLOCATABLE,DIMENSION(:,:,:)   :: vold_l          ! mono material only
        my_real,ALLOCATABLE,DIMENSION(:,:,:)   :: supercellvol_l  ! geometrical (both multi and mono material)
        INTEGER,ALLOCATABLE,DIMENSION(:,:,:)   :: unlinked_cells_l
        INTEGER,ALLOCATABLE,DIMENSION(:)       :: n_unlinked_l
        INTEGER,ALLOCATABLE,DIMENSION(:)       :: imergel
 
        TYPE edge_entity
          INTEGER          ,DIMENSION(2)       :: node
          INTEGER                              :: nbcut                   !number of intersection points on edge
          INTEGER                              :: cutshell(2)
          my_real,          DIMENSION(2)       :: cutcoor                 !local coordinates 1D on edge N1->N2
          my_real,          DIMENSION(3,2)     :: cutpoint                !coordinates x,y,z
          my_real,          DIMENSION(3,2)     :: cutvel                  !velocity of intersection point
          my_real,          DIMENSION(3)       :: vector
          my_real                              :: len
        END TYPE edge_entity
 
        TYPE cut_plane
          my_real                              :: n(3)                    !normal vector
          my_real                              :: b(3)                    !basis
          my_real                              :: scut(1)                 !Cut Area
          my_real                              :: seff(1)                 !Effective Section = Wetted Surface : Cut Aera - Holes         (computed by suming lagrangian projection)
          my_real                              :: p(3,8)                  !Cut Surface Polygon (coordinates)
          my_real                              :: vel(3)
          INTEGER                              :: np                      !Cut Surface Polygon (number of points) ! NP<=6 | NP<=8 for additional Scut (closed surface hypothesis)
        END TYPE
 
        TYPE list_adj
          INTEGER                              :: num                     !number of adjacent elements (level2)
          INTEGER                              :: num_inv(48)             !to define bijection with adjacent cut cell
          INTEGER                              :: iv(48)                  !internal ID
          INTEGER                              :: ib(48)                  !cut cell id
          INTEGER                              :: icell(48)               !secnd cell id
          INTEGER                              :: secndface(48)           !secnd face ID to pickup correct flux
        END TYPE list_adj
 
        TYPE list_secnd
          my_real                              :: vol_unmerged
          INTEGER                              :: num
          INTEGER                              :: numsecndnodes
          INTEGER                              :: fm(24)                  !Main face common face
          INTEGER                              :: fv(24)                  !Secnd corresponding face
          INTEGER                              :: iv(24)                  !brick id which contains secnd cell
          INTEGER                              :: ibv(24)                 !cut cell address in cut cell buffer
          INTEGER                              :: icellv(24)              !secnd cell id
          my_real                              :: vol(24)                 !secnd cell volume
          my_real                              :: surf_v(24)              !Area of linking face
          INTEGER                              :: numnod_cell(24)         !noeuds sommets
          INTEGER                              :: listnodid(24,8)         !Global nodes ID (max8) describing the cell, for each of 9 cells
        END TYPE list_secnd
 
        TYPE polyface_entity
          INTEGER                              :: numpoint                !gives number of points (node+intersec.pts) on a given cell face   (%NumPOINT_Face(6,9))
          INTEGER                              :: nadjcell                !Number of adjacent cell for a given cell 1:9 and a given face 1:6  (%NAdjCell(9,6) )
          INTEGER                              :: adjacent_cell(5)        !index3: IADJ (adjacent cell id 1:%NAdjCell)
          my_real                              :: adjacent_flux(5)        !1 adjacent cell = 1 flux (plus pratique pour maillage non conforme)
          my_real                              :: adjacent_upwflux(5)
          my_real                              :: center(3)               !Centroides des facettes de polyedres
          my_real                              :: surf                    !face0=intersection surface
          my_real                              :: vel(3)                  !velocity on polyhedra faces. Face 0 is for further implementation for porosity. Computed in i22subvol to take into account lagrangian vel.
          my_real                              :: f_face(3)               !FVM force on faces (and also velocity storage)
          my_real                              :: u_n
          my_real                              :: w(3)
        END TYPE polyface_entity
 
        TYPE polyface_entity_2
          INTEGER                              :: numpoint                !gives number of points (node+intersec.pts) on a given cell face   (%NumPOINT_Face(6,9))
          INTEGER                              :: nadjcell                !Number of adjacent cell for a given cell 1:9 and a given face 1:6  (%NAdjCell(9,6) )
          !INTEGER                              :: Adjacent_Cell(5)       !index3: IADJ (adjacent cell id 1:%NAdjCell)
          !my_real                              :: Adjacent_FLUX(5)       !1 adjacent cell = 1 flux (plus pratique pour maillage non conforme)
          !my_real                              :: Adjacent_upwFLUX(5)
          my_real                              :: center(3)               !polyhedra Centroid
          my_real                              :: surf                    !face0=intersection surface
          my_real                              :: f_face(3)               !FVM force on faces (and also velocity storage)
          my_real                              :: u_n(9)       !%U_N(index<=4) face0 is cut surface : can be the union of max. 4 cut polygons
          !possible optimization : 1 single U_N per polyhedron ipoly <9
          my_real                              :: w(3)
 
        END TYPE polyface_entity_2
 
        TYPE mat_entity
          my_real                              :: rho
          my_real                              :: rhoe
          my_real                              :: rhou(3)
          my_real                              :: sig(6)
          my_real                              :: ssp
          my_real, ALLOCATABLE,DIMENSION(:)    :: uvar  !law37:size=5   !law51:size=N0PHAS+4*NVPHAS   !0:otherwise
        END TYPE mat_entity
 
        TYPE poly_entity
          INTEGER                              :: numnod                  !vertexes (%NumNOD_Cell(9))
          INTEGER                              :: numpoint                !total number of intersection point which describes the volume (%NumPOINT_Cell(9))
          INTEGER                              :: listnodid(8)            !Global nodes ID (max8) describing the cell, for each of 9 cells  (%ListNodID(8,9))
          INTEGER                              :: ismain                  !=1 main =0:secnd  (index is icell)   (%IsMain(9))
          INTEGER                              :: whereismain(4)          !locate main cell to link with index1 : ICELL,  index2: (1)face, (2):icellm, (3):IDm, (4):IBm !(%WhereIsMain(9,4))
          INTEGER                              :: id_free_node            !post-treatment  (!ID_FREE_NODE(9) )
          my_real                              :: adjacent_flu1           ! Adjacent_FLU1(9)
          my_real                              :: ddvol                   !d/dT (DELTA V)
          my_real                              :: ddvol_upw
          my_real                              :: phi                     !data for convection (mass or energy)
          my_real                              :: dphi                    !convected increment (mass or energy)
          my_real                              :: vfracm(4)               !main cell VFRAC
          my_real                              :: vnew                    !subvolumes (if exists)
          my_real                              :: vold
          my_real                              :: old_vnew
          my_real                              :: cellcenter(3)           !centroid of control volume
          my_real                              :: dvol(0:2)               !index1:0:previous,1=real,2=prediction   index2 :0-9=icell (0=supercell)
          my_real                              :: cutface
          my_real                              :: dm
          TYPE(polyface_entity_2)              :: face0
          TYPE(polyface_entity)                :: face(1:6)
          !TYPE(MAT_ENTITY)                     :: MAT
        END TYPE poly_entity
 
        TYPE node_entity
          INTEGER                              :: old_whichcell           ! [1,8] -> [1,9]
          INTEGER                              :: whichcell               !node local id in[1,8] -> [1,9] = local cut cell id.
          INTEGER                              :: nodwasmain              !=1 belonged to the main cell, 0 otherwise. (index is brick node)
          INTEGER                              :: wherewasmain            !index=node : result=face
        END TYPE node_entity
 
        TYPE brick_entity
          INTEGER                              :: id
          INTEGER                              :: ng
          INTEGER                              :: idloc
          INTEGER                              :: icode
          INTEGER                              :: old_icode
          INTEGER                              :: idble
          INTEGER                              :: nbits
          INTEGER                              :: npqts
          INTEGER                              :: nbcut                  !number of cut planes
          INTEGER                              :: mainid
          INTEGER                              :: wascut
          INTEGER                              :: newinbuffer
          INTEGER                              :: oldmainstrongnode
          INTEGER                              :: mlw
          INTEGER                              :: itask
          INTEGER                              :: secid_cell(8)           !cut local id in [1,8] -> [1,52] = sec type. Each possible cut is numbered from 1 to 52 ; 8 nodes => 8 possible cut cells
          INTEGER                              :: seg_add_lft             !EL FIRST : first lagrangian face in CAND_E
          INTEGER                              :: seg_add_llt             !EL LAST  : last lagrangian face in CAND_E
          INTEGER                              :: adjacent_brick(6,5)     !index2: 1:IV, 2:NGV, 3:IDLOCV, 4:IBV, 5:IFV
          INTEGER                              :: mergetarget(3,5)        !index1 (1)=Face (2)=MCELL, (3)=IB
          INTEGER                              :: ntarget                 !number of merged location
          INTEGER                              :: closedsurf(6)           !tells if brick face has a double cut (only for multiple cut, in this case cell 9 has no adjacent cells, only for contact forces to get pressure differential)
          INTEGER                              :: ismergetarget
          my_real                              :: n(6,3)                  !unitary normal vector on each 6 faces
          my_real                              :: vnew_scell              !SuperCell volume (if exists)
          my_real                              :: vold_scell              !volume main (if exists)    : need for epsdot in sigeps51, merge/demerge, link switching , srho3
          my_real                              :: scellcenter(3)          !SuperCell centroid
          my_real                              :: face_brick(6)
          my_real                              :: uncutvol                !brick volume. the one in gbuf can up extended if secnd is merged.
          my_real                              :: fcell(3)                !FVM force on centroid
          my_real                              :: fext_cell(3)
          my_real                              :: sig(0:6)                !0 is pressure 1:6 is tensor
          my_real                              :: rhoc
          my_real                              :: mach
          my_real                              :: dvol
          my_real                              :: poly9wonodes(6,2)       !Poly9woNodes(1:6,1) : tells if polygona has no nodes
          !Poly9woNodes(1:6,2) : tells about adjacent cell  0=closed surface  1=communication with adjacent cell
          TYPE(poly_entity)                    :: poly(9)
          TYPE(node_entity)                    :: node(8)
          TYPE(list_adj)                       :: adj_elems               !to be cleaned ?
          TYPE(cut_plane)                      :: pcut(16)
          TYPE(list_secnd)                     :: secndlist               !list of secnd cells linked to the given main cell
          TYPE(edge_entity)                    :: Edge(12)                !list of edges for velocity or coordinates
          CHARACTER*14                         :: SECTYPE(8)
          TYPE(mat_entity)                     :: bakmat
        END TYPE brick_entity
 
 
        TYPE(list_secnd),DIMENSION(:,:),ALLOCATABLE    :: old_secndlist
 
 
      END MODULE