powder__burn_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/.

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

!||    powder_burn_mod   ../common_source/eos/powder_burn.F

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

!||    eosmain           ../common_source/eos/eosmain.F

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

       MODULE powder_burn_mod

       IMPLICIT NONE

       CONTAINS

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

!||    powder_burn        ../common_source/eos/powder_burn.F

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

!||    eosmain            ../common_source/eos/eosmain.F

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

!||    constant_mod       ../common_source/modules/constant_mod.F

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

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


      SUBROUTINE powder_burn    ( NPROPM, NUMMAT ,

     1                            IFLAG , NEL    ,PM    ,OFF       , EINT  ,MU          ,MU2 ,

     2                            ESPE  , DVOL   ,DF    ,VNEW      , MAT   ,PSH         ,

     3                            PNEW  , DPDM   ,DPDE  ,MAT_PARAM ,

     4                            VAREOS, NVAREOS,DT1   ,RHO0      ,BFRAC )

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

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

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

C This is a template to introduce numerical solving of Deflagration EoS (experimental) based on LAW105 Implementation

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

!! \details STAGGERED SCHEME IS EXECUTED IN TWO PASSES IN EOSMAIN : IFLG=0 THEN IFLG=1

!! \details COLLOCATED SCHEME IS DOING A SINGLE PASS : IFLG=2

!! \details

!! \details  STAGGERED SCHEME

!! \details     EOSMAIN / IFLG = 0 : DERIVATIVE CALCULATION FOR SOUND SPEED ESTIMATION c[n+1] REQUIRED FOR PSEUDO-VISCOSITY (DPDE:partial derivative, DPDM:total derivative)

!! \details     MQVISCB            : PSEUDO-VISCOSITY Q[n+1]

!! \details     MEINT              : INTERNAL ENERGY INTEGRATION FOR E[n+1] : FIRST PART USING P[n], Q[n], and Q[n+1] CONTRIBUTIONS

!! \details     EOSMAIN / IFLG = 1 : UPDATE P[n+1], T[N+1]

!! \details                          INTERNAL ENERGY INTEGRATION FOR E[n+1] : LAST PART USING P[n+1] CONTRIBUTION

!! \details                            (second order integration dE = -P.dV where P = 0.5(P[n+1] + P[n]) )

!! \details  COLLOCATED SCHEME

!! \details     EOSMAIN / IFLG = 2 : SINGLE PASS FOR P[n+1] AND DERIVATIVES

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

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

C   M o d u l e s

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

        USE matparam_def_mod, ONLY : matparam_struct_

        use constant_mod , only : zero, one, zep07, three100, em04

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

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

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

        implicit none

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

C   I n c l u d e   F i l e s

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

#include      "my_real.inc"

#include      "comlock.inc"

#include      "com06_c.inc"

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

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

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

        INTEGER,INTENT(IN) :: NVAREOS,NUMMAT,NPROPM

        INTEGER MAT(NEL), IFLAG, NEL

        my_real,INTENT(IN) :: PM(NPROPM,NUMMAT),OFF(NEL),MU(NEL),MU2(NEL),DVOL(NEL),DF(NEL),VNEW(NEL),RHO0(NEL)

        my_real,INTENT(INOUT) :: PNEW(NEL),DPDM(NEL),DPDE(NEL),EINT(NEL),ESPE(NEL)

        my_real,INTENT(INOUT) :: vareos(nel,nvareos),bfrac(nel)

        my_real,INTENT(IN) :: dt1

        TYPE(matparam_struct_), INTENT(IN) :: MAT_PARAM !< material data structure

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

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

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

        INTEGER I

        my_real :: bulk,P0,PSH(NEL),DD,EG,Gr,CC,ALPHA,FSCALE_B,FSCALE_P,FSCALE_G,FSCALE_RHO, C1,C2

        my_real :: TMP, TMP2

        my_real :: compac

        my_real :: total_bfrac

        my_real :: dpdm_gas, dpdm_powder

        my_real :: mass,ps,pg,rho_s, rho_g,pold

        my_real :: rho(nel)

        INTEGER :: funcb,funcg

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

C   E x t e r n a l   F u n c t i o n s

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

        my_real,EXTERNAL :: finter

        ! EXTERNAL FINTER

        !      Y = FINTER(IFUNC(J),X,NPF,TF,DYDX)

        !      Y       : y = f(x)

        !      X       : x

        !      DYDX    : f'(x) = dy/dx

        !      IFUNC(J): FUNCTION INDEX

        !            J : FIRST(J=1), SECOND(J=2) .. FUNCTION USED FOR THIS LAW

        !      NPF,TF  : FUNCTION PARAMETER

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

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

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


        ! --- STAGGERED SCHEME IFLG=0

        IF(iflag == 0) THEN

          bulk       = mat_param%EOS%UPARAM(01)

          p0         = mat_param%EOS%UPARAM(02)

          psh(1:nel) = mat_param%EOS%UPARAM(03)

          dd         = mat_param%EOS%UPARAM(04)

          eg         = mat_param%EOS%UPARAM(05)

          gr         = mat_param%EOS%UPARAM(06)

          cc         = mat_param%EOS%UPARAM(07)

          alpha      = mat_param%EOS%UPARAM(08)

          fscale_b   = mat_param%EOS%UPARAM(09)

          fscale_p   = mat_param%EOS%UPARAM(10)

          fscale_g   = mat_param%EOS%UPARAM(11)

          fscale_rho = mat_param%EOS%UPARAM(12)

          c1         = mat_param%EOS%UPARAM(13)

          c2         = mat_param%EOS%UPARAM(14)

          funcb      = mat_param%EOS%FUNC(1)

          funcg      = mat_param%EOS%FUNC(2)

          compac = one - zep07  !1-0.07 = 0.93

          total_bfrac = zero

          IF(dt1 == zero)THEN

            DO i=1,nel

              espe(i) = eg

              eint(i) = eg*rho0(i)*vnew(i)

              vareos(i,1) = p0                 !PS

              vareos(i,2) = zero               !PG

              vareos(i,3) = rho0(i)/compac     !RHO_S

              vareos(i,4) = zero               !RHO_G

              vareos(i,5) = p0                 !POLD

              vareos(i,6) = zero               !F(t_old)

              vareos(i,7) = rho0(i)*vnew(i)    !Mass0

              dpdm(i) = bulk

              dpde(i) = zero

            ENDDO

          ENDIF

          DO i=1,nel

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

            ! INIT.                          !

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

            rho(i) = rho0(i) * (one + mu(i))

            mass = rho(i)*vnew(i)

            espe = eint(i)/mass

            ps   =vareos(i,1)

            pg   =vareos(i,2)

            rho_s=vareos(i,3)

            rho_g=vareos(i,4)

            pold =vareos(i,5)

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

            ! SOUND SPEED DERIVATIVE                                             !

            !  DPDM ! dP/dmu total derivative                                    !

            !  dPdE ! partial derivative (where E=Eint/V0 , and rho0V0=rho.V)    !

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

            tmp = (one+mu(i))*rho0(i)/dd

            tmp2 = (one+mu(i)) ; tmp2=tmp2*tmp2

            dpdm_gas = eg*exp(tmp)*(one+tmp) + pg/tmp2*(one+mu(i))*exp(tmp)

            dpdm_powder = bulk

            dpdm(i) = total_bfrac * dpdm_gas + (one-total_bfrac)*dpdm_powder   !total derivative of mixture

            dpde(i) = total_bfrac * (one+mu(i)*tmp)                            !partial derivative of mixture

          ENDDO


        ELSEIF(iflag == 1) THEN


        ELSEIF (iflag == 2) THEN

          DO i=1, nel

            IF (vnew(i) > zero) THEN

              pnew(i) = zero

              dpdm(i) = zero !total derivative

              dpde(i) = zero !partial derivative

            ENDIF

          ENDDO


        ENDIF

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

        RETURN


      END SUBROUTINE powder_burn

      END MODULE powder_burn_mod

constant_mod
Definition constant_mod.F:355

powder_burn_mod
Definition powder_burn.F:28

powder_burn_mod::powder_burn
subroutine powder_burn(npropm, nummat, iflag, nel, pm, off, eint, mu, mu2, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde, mat_param, vareos, nvareos, dt1, rho0, bfrac)
Definition powder_burn.F:44