multi_bilan.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/.
!||====================================================================
!||    multi_bilan           ../engine/source/multifluid/multi_bilan.f
!||--- called by ------------------------------------------------------
!||    multi_globalize       ../engine/source/multifluid/multi_globalize.F
!||    multi_update_global   ../engine/source/multifluid/multi_update_global.f
!||--- uses       -----------------------------------------------------
!||    multi_fvm_mod         ../common_source/modules/ale/multi_fvm_mod.F90
!||====================================================================
      SUBROUTINE multi_bilan(PARTSAV, VOL, MULTI_FVM,
     .     NEL, NFT, IPARTS, GRESAV, IGRTH, GRTH)
C-----------------------------------------------
C     M o d u l e s
C-----------------------------------------------
      USE multi_fvm_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real, INTENT(INOUT) :: partsav(npsav, *), gresav(*)
      my_real, INTENT(IN) :: vol(*)
      TYPE(multi_fvm_struct), INTENT(IN) :: MULTI_FVM
      INTEGER, INTENT(IN) :: NEL, NFT, IPARTS(*), IGRTH(*), GRTH(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, M
      my_real :: ei_sav(nel), ek_sav(nel), 
     .     xm_sav(nel), ym_sav(nel), zm_sav(nel), rbidon(1), va2(nel),
     .     mas(nel), off(nel), vx(nel), vy(nel), vz(nel)
      
      rbidon(1) = zero
      off(1:nel) = one
      DO i = 1, nel
         mas(i) = vol(i) * multi_fvm%RHO(i + nft)
         vx(i) = multi_fvm%VEL(1, i + nft)
         vy(i) = multi_fvm%VEL(2, i + nft)
         vz(i) = multi_fvm%VEL(3, i + nft)
         va2(i) = vx(i)**2 + vy(i)**2 + vz(i)**2
      ENDDO
      m = iparts(1)
      DO i=1,nel
         partsav(1, m) = partsav(1, m) + multi_fvm%EINT(i + nft) * vol(i)
         partsav(2, m) = partsav(2, m) + half * mas(i) * va2(i)
         partsav(6, m) = partsav(6, m) + mas(i)
         partsav(3, m) = partsav(3, m) + mas(i) * vx(i)
         partsav(4, m) = partsav(4, m) + mas(i) * vy(i)
         partsav(5, m) = partsav(5, m) + mas(i) * vz(i)
      ENDDO
 
      END SUBROUTINE multi_bilan