Moving format

This format is used in VIGIE in the case of flow in a moving boundary, using non-structured meshes. This can also be used for reacting flows. Let's start by the definition of parameters and arrays used:

Note: The writing format can be modified.

Parameters:

• NS : Number of nodes.
• NT : Number of thetraedra.
• NFAC : Number of faces on boundary.
• NDIM : Space dimension.
• NELM : Number of points per volume.
• NE : Number of chimical species ($>= 1$).
• NVAR : Number of variables per node.

Arrays:

• COOR : Array of 3 point coordinates (X, Y, Z).
• NU : Array of indices of 4 vertices of tetraedra.
• LOGFAC : Array of logic of frontier faces.
• NSFAC : Array of indices of 3 vertices per faces.
• UA : Array of mesh nodes values.
  • UA (1 ,*)= $\rho$ : density
  • UA (2 ,*)= $\rho*U_{x}$ : density * ((X) speed component)
  • UA (3 ,*)= $\rho*U_{y}$ : density * ((Y) speed component)
  • UA (4 ,*)= $\rho*U_{z}$ : density * ((Z) speed component)
  • UA (5 ,*)= $E$ : Total energy minus the sum of the massic formation of enthalpie
  • UA (5+k,*)= $\rho_{k}$ : partial density of species $k$: k=1, NE-1
• OMEGA : Chimical reaction rate
• PRES : Pressure
• TEMP : Temperature

Moving format: An example:

C
        NDIM = 3
        NELM = 4
C
      REAL COOR(NDIM, NS),  UA(NVAR, NS)
      INTEGER  NU(4,NT), LOGFAC(NFAC), NSFAC(3,NFAC)
C
      INET1 = 60
      INET2 = 59
C      
      NE : read in datas
      NVAR = 4 + NE 
C
C      GEOMETRY DESCRIPTION FILE
C      -------------------------------------
      REWIND(INET1)
C
      WRITE( INET1, 112 ) NS, NT, NFAC
      WRITE( INET1, 115 ) ( ( COOR(I,IS), I= 1, NDIM ), IS= 1, NS )
      WRITE( INET1, 116 ) ( ( NU(K,JT), K= 1, NELM ), JT= 1, NT )
      WRITE( INET1, 116 ) ( LOGFAC(IFAC), IFAC= 1, NFAC )
      WRITE( INET1, 116 ) ( ( NSFAC(I,IFAC), I= 1, NELM - 1 ), IFAC= 1, NFAC )
C      
C      SOLUTION FILE AT TIME  T
C      --------------------------
C
      REWIND(INET2)
C
        WRITE( INET2, 115 ) ( (COOR(I,IS), I= 1, NDIM ), IS= 1, NS )
        WRITE( INET2, 115 ) ( (UA(I,IS), I= 1, NVAR ), IS= 1, NS )
        IF( NE .NE. 1 ) THEN
          READ(INET2,115) ( OMEGA(IS) ,  IS= 1,  NS )
          READ(INET2,115) (  PRES(IS) ,  IS= 1,  NS )
          READ(INET2,115) (  TEMP(IS) ,  IS= 1,  NS )
         ENDIF
        WRITE( INET2, 112 ) KT, T
C
 112  FORMAT( I6, E12.5 )
 113  FORMAT( 3(I6) )
 115  FORMAT( 5(E15.8,2x) )
 116  FORMAT( 10(I8,2x) )

Note: If the ``COOR'' are the same at time T, rewrite them.