Merge remote-tracking branch 'firemodels/master' into FireX

Author: cxp484

Manuals/Bibliography/FDS_general.bib

@@ -3230,7 +3230,18 @@ @inproceedings{Hu:2005
 	booktitle = {Fire Safety Science---Proceedings of the Eighth International Symposium},
 	pages = {1193-1204},
 	publisher = {International Association for Fire Safety Science},
-	year = {2005},
+	year = {2005}
+}
+
+@inproceedings{Hu:2007,
+        title = {{Towards large eddy simulations of flame extinction and carbon monoxide emission in compartment fires}},
+        author = {Z. Hu and Y. Utiskul and J.G. Quintiere and A. Trouve},
+        doi = {10.1016/j.proci.2006.08.053},
+        booktitle = {Proceedings of the Combustion Institute},
+        pages = {2437-2545},
+        volume = {31},
+        publisher = {Elsevier},
+        year = {2007}
 }
 
 @article{Huang:1963,

Manuals/FDS_Validation_Guide/Experiment_Chapter.tex

@@ -3677,6 +3677,33 @@ \section{USN High Bay Hangar Experiments}
 
 
 
+\section{Utiskul Compartment}
+\label{Utiskul_Compartment_Description}
+
+Hu and Utiskul conducted experiments and analysis aimed at characterizing the dynamics of compartment fires under poorly ventilated conditions~\cite{Hu:2005,Hu:2007}. The study considers four cases corresponding to different values of the global equivalence ratio. The experimental utilized a 40~cm by 40~cm by 40~cm compartment with two horizontal vents of equal size located at the top and bottom of one of the compartment walls. The vent width varied between 2~cm and 40~cm; the vent height varied between 1~cm and 3~cm. The fire was fueled by a circular heptane pan located at the center of the compartment floor. The compartment was instrumented with sensors and probes that provided the time history of the fuel mass loss rate, local temperatures, surface heat fluxes, and floor and ceiling-level concentrations of O$_2$, CO$_2$ and CO. Table~\ref{Hu_Summary} summarizes the four experiments.
+
+\begin{table}[h!]
+\centering
+\caption[Summary of the Hu/Utiskul experiments]{Summary of the Hu/Utiskul experiments.}
+\begin{tabular}{|c|c|c|c|}
+\hline
+Case    & Vent Dimensions (cm)  & Pan Diameter (cm)  & Fire Regime                                           \\ \hline \hline
+1       & 3 $\times$ 40         & 9.5                & steady; well-ventilated; flame stabilized on burner   \\ \hline
+2       & 1 $\times$ 40         & 19                 & steady; under-ventilated; flame stabilized at vents   \\ \hline
+3       & 3 $\times$ 10         & 9.5                & unsteady; under-ventilated; oscillating flame         \\ \hline
+4       & 1 $\times$ 2          & 9.5                & unsteady; under-ventilated; flame extinction          \\ \hline
+\end{tabular}
+\label{Hu_Summary}
+\end{table}
+
+\subsubsection{Modeling Notes}
+
+In the simulations of the experiments, it is assumed that the combustion of n-heptane can be simplified to two fast reactions, the first converting fuel to CO and soot, and the second converting CO and soot to CO$_2$. By default, 2/3 of the carbon in the fuel is converted to CO in the first step, the remaining 1/3 to soot. The heats of combustion for the reactions are calculated directly from the heats of formation of the individual molecules. The measured fuel mass loss rate was specified in the simulations.
+
+In the simulations, the model geometry included the compartment interior plus a comparable volume at the exterior to allow for a natural flow into and out of the compartment.
+
+
+
 \section{UWO Wind Tunnel Experiments}
 \label{UWO_Wind_Tunnel_Description}
 

Manuals/FDS_Validation_Guide/Species_Chapter.tex

@@ -1873,6 +1873,13 @@ \subsection{NIST Pool Fires}
 \clearpage
 
 \subsection{Utiskul Compartment}
+\label{sec:Utiskul_Compartment}
+
+A summary of the experiments is given in Sec.~\ref{Utiskul_Compartments_Description}. Figures~\ref{Utiskul_1} through \ref{Utiskul_4} display comparisons of measured and predicted compartment and ceiling jet temperatures, pressure, heat flux, and species concentration. The fuel mass loss rate is specified, not predicted, by the model. 
+
+The species concentrations are measured at extraction points 2~cm above the floor and 2~cm below the ceiling, 10~cm from the vents along the compartment centerline. The ceiling jet temperatures are taken 2~cm below the ceiling along the centerline, at locations 2.5~cm from the upper vent (TC2), 2.5~cm from the back wall (TC13), and at the ceiling center (TC17). The pressures are taken 1.5~cm above the floor and below the ceiling, 2.5~cm from the side wall, and 5~cm from the front wall, where the vents are located.
+
+Heat fluxes are measured at the center of the fuel pan looking up (HF1), 5~cm from the back wall looking up (HF2), halfway up the side wall aimed horizontally at the opposite wall (HF3), and at the center of the ceiling looking down (HF4).
 
 \begin{figure}[p]
 \begin{tabular*}{\textwidth}{l@{\extracolsep{\fill}}r}

Source/cons.f90

@@ -499,7 +499,7 @@ MODULE GLOBAL_CONSTANTS
 
 LOGICAL, ALLOCATABLE, DIMENSION(:) :: RADIATION_COMPLETED  !< Indicates that the radiation field is completely updated
 
-LOGICAL :: RANDOMIZE_RADIATION_DIRECTIONS=.FALSE.          !< If TRUE, randomly rotate solid angles each full RTE solve
+LOGICAL :: RANDOMIZE_RADIATION_DIRECTIONS=.FALSE.           !< If TRUE, randomly rotate solid angles each full RTE solve
 LOGICAL :: ALLOW_RANDOM_RADIATION_ROTATION=.FALSE.         !< A derived variable, to block random rotation fo cyl and 2d cases
 
 INTEGER :: NUMBER_SPECTRAL_BANDS=0                         !< Number of wavelength bands for rad solver (1 for gray gas)
@@ -830,6 +830,8 @@ MODULE GLOBAL_CONSTANTS
 ! VENT array
 REAL(EB), ALLOCATABLE, DIMENSION(:) :: VENT_TOTAL_AREA  !< Array holding grid-snapped areas for all vents
 
+LOGICAL :: HAS_MIRROR_BOUNDARY=.FALSE. !< To decide if mirror boundary exist for the setup.
+
 END MODULE GLOBAL_CONSTANTS
 
 

Source/main.f90

@@ -298,11 +298,10 @@ PROGRAM FDS
 DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
    CALL INITIALIZE_RADIATION_EXCHANGE(NM)
 ENDDO
-
+CALL CHECK_FOR_MIRROR_BOUNDARY()
 IF (MY_RANK==0 .AND. VERBOSE) CALL VERBOSE_PRINTOUT('Completed INITIALIZE_RADIATION_EXCHANGE')
 
 ! Initialize persistent MPI sends and receives and allocate buffer arrays.
-
 IF(RADIATION .AND. .NOT. ALLOW_RANDOM_RADIATION_ROTATION) THEN
    CALL CALCULATE_DIRECTION_COEFFICIENTS()
    CALL ALLOCATE_RADIAITON_RECV_PKG()

Source/radi.f90

@@ -2763,7 +2763,8 @@ MODULE RAD
 IMPLICIT NONE (TYPE,EXTERNAL)
 PRIVATE
 
-PUBLIC INIT_RADIATION,COMPUTE_RADIATION,BLACKBODY_FRACTION,CALCULATE_DIRECTION_COEFFICIENTS,INTERPOLATE_IL
+PUBLIC INIT_RADIATION,COMPUTE_RADIATION,BLACKBODY_FRACTION,CALCULATE_DIRECTION_COEFFICIENTS,INTERPOLATE_IL,&
+       CHECK_FOR_MIRROR_BOUNDARY
 
 REAL(EB) :: TYY_FAC
 INTEGER :: N_KAPPA_T=44                             !< Number of temperature points in absorption coefficient look-up table
@@ -3396,18 +3397,40 @@ SUBROUTINE CALCULATE_FVM_ANGLES()
       ENDDO
    ENDDO
 ENDDO
+ 
 
 END SUBROUTINE CALCULATE_FVM_ANGLES
 
+!> ! Determine whether a mirror boundary exists.
+SUBROUTINE CHECK_FOR_MIRROR_BOUNDARY()
+LOGICAL :: HAS_MIRROR_BOUNDARY_LOCAL
+INTEGER :: NM,IW,IERR
+TYPE(WALL_TYPE),POINTER :: WC
+HAS_MIRROR_BOUNDARY_LOCAL=.FALSE.
+MESH_LOOP: DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
+   CALL POINT_TO_MESH(NM)
+   DO IW=1,N_EXTERNAL_WALL_CELLS+N_INTERNAL_WALL_CELLS
+      WC=>WALL(IW)
+      IF (WC%BOUNDARY_TYPE==MIRROR_BOUNDARY) THEN
+         HAS_MIRROR_BOUNDARY_LOCAL=.TRUE.
+         EXIT MESH_LOOP
+      ENDIF
+   ENDDO
+ENDDO MESH_LOOP
+CALL MPI_ALLREDUCE(HAS_MIRROR_BOUNDARY_LOCAL,HAS_MIRROR_BOUNDARY,1,MPI_LOGICAL,MPI_LOR,MPI_COMM_WORLD,IERR)
+END SUBROUTINE CHECK_FOR_MIRROR_BOUNDARY
+
 
 
 !> \Calculate direction coefficients for with and without random rotations
 !> as described in the FDS Tech. Ref. Guide Vol. 1 Sec. 6.2.2.
 SUBROUTINE CALCULATE_DIRECTION_COEFFICIENTS()
 USE COMP_FUNCTIONS, ONLY : CURRENT_TIME
 
-INTEGER  :: NRA,N,IO,I,IERR
+INTEGER  :: NRA,N,NN,NI,K,IO,I,IERR
 REAL(EB) :: TNOW,AXIS(3),MERIDIAN(3), AZIMUTH(3),E1(3),E2(3),REF(3),PSI,DLO,MAGTMP,RNDNUM(3)
+REAL(EB) :: PROJ_SUM_X, PROJ_SUM_Y, PROJ_SUM_Z
+REAL(EB) :: REFLECTED_DIRECTION(3), DOT_MAX, DOT_TMP
 
 TNOW=CURRENT_TIME()
 
@@ -3485,6 +3508,24 @@ SUBROUTINE CALCULATE_DIRECTION_COEFFICIENTS()
    ENDDO
 ENDIF
 
+! Normalize such that sum of DLX/DLY/DLZ is PI.
+IF (ALLOW_RANDOM_RADIATION_ROTATION) THEN
+
+   PROJ_SUM_X = 0.5_EB*SUM(ABS(DLX(1:NRA)))
+   PROJ_SUM_Y = 0.5_EB*SUM(ABS(DLY(1:NRA)))
+   PROJ_SUM_Z = 0.5_EB*SUM(ABS(DLZ(1:NRA)))
+
+   IF (PROJ_SUM_X > TWO_EPSILON_EB) &
+      DLX(1:NRA) = PI/PROJ_SUM_X*DLX(1:NRA)
+
+   IF (PROJ_SUM_Y > TWO_EPSILON_EB) &
+      DLY(1:NRA) = PI/PROJ_SUM_Y*DLY(1:NRA)
+
+   IF (PROJ_SUM_Z > TWO_EPSILON_EB) &
+      DLZ(1:NRA) = PI/PROJ_SUM_Z*DLZ(1:NRA)
+
+ENDIF
+
 ! Set (wall normal)*(angle vector) value
 DO N = 1,NRA
    DLN( 0,N) = 0._EB !prevent undefined variable errors
@@ -3496,6 +3537,27 @@ SUBROUTINE CALCULATE_DIRECTION_COEFFICIENTS()
    DLN( 3,N) =  DLZ(N)
 ENDDO
 
+
+! Calculate mirror directions
+IF (HAS_MIRROR_BOUNDARY .AND. ALLOW_RANDOM_RADIATION_ROTATION) THEN
+   DO N = 1,NRA
+      DO K = 1,3
+         REFLECTED_DIRECTION = DLANG(:,N)
+         REFLECTED_DIRECTION(K) = -REFLECTED_DIRECTION(K)
+         DOT_MAX = -HUGE(1._EB)
+         NI = 1
+         DO NN = 1,NRA
+            DOT_TMP = DOT_PRODUCT(REFLECTED_DIRECTION,DLANG(:,NN))
+            IF (DOT_TMP > DOT_MAX) THEN
+               DOT_MAX = DOT_TMP
+               NI = NN
+            ENDIF
+         ENDDO
+         DLM(N,K) = NI
+      ENDDO
+   ENDDO
+ENDIF
+
 ! In axially symmetric case, each angle represents two symmetric angles. So weight the intensities by two.
 WEIGH_CYL = 1._EB
 IF (CYLINDRICAL) THEN

Validation/Process_All_Output.sh

@@ -194,6 +194,7 @@ PROCESS USCG_HAI
 PROCESS USFS_Catchpole
 PROCESS USFS_Corsica
 PROCESS USN_Hangars
+PROCESS Utiskul_Compartment
 PROCESS UWO_Wind_Tunnel
 PROCESS Vettori_Flat_Ceiling
 PROCESS Vettori_Sloped_Ceiling

Validation/Run_Serial.sh

@@ -39,6 +39,7 @@ cd Ulster_SBI; ./Run_All.sh $OPTIONS; cd ..
 cd UMD_Polymers; ./Run_All.sh $OPTIONS; cd ..
 cd UMD_Burning_Rate_Emulator; ./Run_All.sh $OPTIONS; cd ..
 cd USCG_HAI; ./Run_All.sh $OPTIONS; cd ..
+cd Utiskul_Compartment; ./Run_All.sh $OPTIONS; cd ..
 cd Vettori_Flat_Ceiling; ./Run_All.sh $OPTIONS; cd ..
 cd Vettori_Sloped_Ceiling; ./Run_All.sh $OPTIONS; cd ..
 cd VTT_Sprays; ./Run_All.sh $OPTIONS; cd ..

Validation/Utiskul_Compartment/FDS_Input_Files/common_properties.txt

@@ -0,0 +1,78 @@
+
+&MESH IJK=60,40,40, XB=-0.2,0.4,-0.2,0.2,0.0,0.4/
+
+&DUMP DT_DEVC=5., DT_HRR=5. /
+
+&REAC FUEL = 'N-HEPTANE', N_SIMPLE_CHEMISTRY_REACTIONS=2 /
+
+&SURF ID        = 'Kaowool Type M Board'
+      MATL_ID   = 'CERAMIC FIBER'
+      COLOR     = 'GRAY'
+      THICKNESS = 0.025
+      DEFAULT   = T /
+
+&MATL ID                 = 'CERAMIC FIBER'
+      FYI                = 'Morgan Advanced Ceramics'
+      CONDUCTIVITY_RAMP  = 'k_fiber'
+      SPECIFIC_HEAT_RAMP = 'c_fiber'
+      DENSITY            =  272.
+      EMISSIVITY         = 0.90 /
+
+&RAMP ID='k_fiber', T= 260., F=0.068 /
+&RAMP ID='k_fiber', T= 538., F=0.098 /
+&RAMP ID='k_fiber', T= 816., F=0.141 /
+&RAMP ID='k_fiber', T=1093., F=0.219 /
+
+&RAMP ID='c_fiber', T= 100., F=1.04  /
+&RAMP ID='c_fiber', T= 550., F=1.14  /
+&RAMP ID='c_fiber', T= 982., F=1.09  /
+&RAMP ID='c_fiber', T=1000., F=1.05  /
+
+&SLCF PBY=0., QUANTITY='TEMPERATURE', CELL_CENTERED=.TRUE./
+&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='N-HEPTANE', CELL_CENTERED=.TRUE./
+&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=.TRUE./
+&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', CELL_CENTERED=.TRUE./
+&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', CELL_CENTERED=.TRUE./
+&SLCF PBY=0., QUANTITY='HRRPUV', CELL_CENTERED=.TRUE./
+&SLCF PBY=0., QUANTITY='VELOCITY', VECTOR=.TRUE./
+
+&DEVC XB=0.20,0.20,-0.2,0.2,0.37,0.4, QUANTITY='MASS FLUX X', SPATIAL_STATISTIC='AREA INTEGRAL', ID='M1'/
+&DEVC XB=0.20,0.20,-0.2,0.2,0.0,0.03, QUANTITY='MASS FLUX X', SPATIAL_STATISTIC='AREA INTEGRAL', ID='M2'/
+
+&DEVC XYZ= 0.15 0.175 0.385, QUANTITY='PRESSURE', ID='P1'/
+&DEVC XYZ= 0.15 0.175 0.015, QUANTITY='PRESSURE', ID='P2'/
+
+&DEVC XB=-0.0475,0.0475,-0.0475,0.0475,0.00,0.00, QUANTITY='BURNING RATE', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='B'/
+
+&DEVC XYZ= 0.1 0.0 0.38, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', ID='G1'/
+&DEVC XYZ= 0.1 0.0 0.02, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', ID='G2'/
+&DEVC XYZ= 0.1 0.0 0.38, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', ID='G3'/
+&DEVC XYZ= 0.1 0.0 0.02, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', ID='G4'/
+&DEVC XYZ= 0.1 0.0 0.38, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', ID='G5'/
+&DEVC XYZ= 0.1 0.0 0.02, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', ID='G6'/
+
+&DEVC XYZ=  0.200 0.0 0.385, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC1'/
+&DEVC XYZ=  0.175 0.0 0.380, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC2'/
+&DEVC XYZ=  0.175 0.0 0.340, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC3'/
+&DEVC XYZ=  0.175 0.0 0.300, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC4'/
+&DEVC XYZ=  0.175 0.0 0.260, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC5'/
+&DEVC XYZ=  0.175 0.0 0.220, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC6'/
+&DEVC XYZ=  0.175 0.0 0.180, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC7'/
+&DEVC XYZ=  0.175 0.0 0.140, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC8'/
+&DEVC XYZ=  0.175 0.0 0.100, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC9'/
+&DEVC XYZ=  0.175 0.0 0.060, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC10'/
+&DEVC XYZ=  0.175 0.0 0.020, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC11'/
+&DEVC XYZ=  0.200 0.0 0.015, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC12'/
+&DEVC XYZ=  -0.175 0.0 0.380, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC13'/
+&DEVC XYZ=  -0.175 0.0 0.260, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC14'/
+&DEVC XYZ=  -0.175 0.0 0.140, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC15'/
+&DEVC XYZ=  -0.175 0.0 0.020, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC16'/
+&DEVC XYZ=  0.000 0.0 0.380, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC17'/
+&DEVC XYZ=  0.000 0.0 0.260, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC18'/
+&DEVC XYZ=  0.000 0.0 0.140, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC19'/
+
+&DEVC XYZ=  0.00 0.00 0.00, IOR=3, QUANTITY= 'GAUGE HEAT FLUX', ID='HF1'/
+&DEVC XYZ= -0.15 0.00 0.00, IOR=3, QUANTITY= 'GAUGE HEAT FLUX', ID='HF2'/
+&DEVC XYZ= -0.20 0.00 0.20, IOR=1, QUANTITY= 'GAUGE HEAT FLUX', ID='HF3'/
+&DEVC XYZ=  0.00 0.00 0.40, IOR=-3, QUANTITY= 'GAUGE HEAT FLUX', ID='HF4'/
+

Validation/Utiskul_Compartment/FDS_Input_Files/hu_1p.fds

@@ -1,13 +1,7 @@
-&HEAD CHID='hu_1p',TITLE='Hu et al. IAFSS (2005) compartment fire Case 1 prescribed mass loss rate'/
-
-&MESH IJK=60,40,40, XB=-0.2,0.4,-0.2,0.2,0.0,0.4/
+&HEAD CHID='hu_1p', TITLE='Hu et al. IAFSS (2005) compartment fire Case 1 prescribed mass loss rate'/
 
 &TIME T_END=350./
 
-&DUMP DT_DEVC=5., DT_HRR=5. /
-
-&REAC FUEL = 'N-HEPTANE', N_SIMPLE_CHEMISTRY_REACTIONS=2 /
-
 &VENT XB=0.40,0.40,-0.20,0.20,0.00,0.4, SURF_ID='OPEN' /
 &OBST XB=0.20,0.22,-0.20,0.20,0.00,0.4 /
 &HOLE XB=0.19,0.23,-0.20,0.20,0.37,0.4 /
@@ -17,76 +11,13 @@
 
 &SURF ID='FIRE', SPEC_ID(1)='N-HEPTANE', MASS_FLUX(1)=0.02, RAMP_MF(1)='case1'/
 
-&RAMP ID='case1' T=0, F=0.0/
-&RAMP ID='case1' T=10, F=0.6/
-&RAMP ID='case1' T=70, F=1.0/
-&RAMP ID='case1' T=250, F=1.0/
+&RAMP ID='case1' T=  0, F=0.0 /
+&RAMP ID='case1' T= 10, F=0.6 /
+&RAMP ID='case1' T= 70, F=1.0 /
+&RAMP ID='case1' T=250, F=1.0 /
 &RAMP ID='case1' T=280, F=0.75/
-&RAMP ID='case1' T=350, F=0.0/
-
-&SURF ID        = 'INSULATION'
-      MATL_ID   = 'CERAMIC FIBER'
-      COLOR     = 'GRAY'
-      THICKNESS = 0.020
-      BACKING   = 'EXPOSED'
-      DEFAULT   = .TRUE. /
-
-&MATL ID                 = 'CERAMIC FIBER'
-      CONDUCTIVITY_RAMP  = 'k_fiber'
-      SPECIFIC_HEAT      = 1.04
-      DENSITY            =  128.
-      EMISSIVITY         = 0.97 /
-&RAMP ID='k_fiber',T= 20.,F=0.09 /
-&RAMP ID='k_fiber',T=300.,F=0.09 /
-&RAMP ID='k_fiber',T=600.,F=0.17 /
-&RAMP ID='k_fiber',T=900.,F=0.25 /
-
-&SLCF PBY=0., QUANTITY='TEMPERATURE', CELL_CENTERED=.TRUE./
-&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='N-HEPTANE', CELL_CENTERED=.TRUE./
-&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=.TRUE./
-&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', CELL_CENTERED=.TRUE./
-&SLCF PBY=0., QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', CELL_CENTERED=.TRUE./
-&SLCF PBY=0., QUANTITY='HRRPUV', CELL_CENTERED=.TRUE./
-&SLCF PBY=0., QUANTITY='VELOCITY', VECTOR=.TRUE./
-
-&DEVC XB=0.20,0.20,-0.2,0.2,0.37,0.4, QUANTITY='MASS FLUX X', SPATIAL_STATISTIC='AREA INTEGRAL', ID='M1'/
-&DEVC XB=0.20,0.20,-0.2,0.2,0.0,0.03, QUANTITY='MASS FLUX X', SPATIAL_STATISTIC='AREA INTEGRAL', ID='M2'/
-
-&DEVC XYZ= 0.15 0.175 0.385, QUANTITY='PRESSURE', ID='P1'/
-&DEVC XYZ= 0.15 0.175 0.015, QUANTITY='PRESSURE', ID='P2'/
-
-&DEVC XB=-0.0475,0.0475,-0.0475,0.0475,0.00,0.00, QUANTITY='BURNING RATE', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='B'/
-
-&DEVC XYZ= 0.1 0.0 0.38, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', ID='G1'/
-&DEVC XYZ= 0.1 0.0 0.02, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', ID='G2'/
-&DEVC XYZ= 0.1 0.0 0.38, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', ID='G3'/
-&DEVC XYZ= 0.1 0.0 0.02, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', ID='G4'/
-&DEVC XYZ= 0.1 0.0 0.38, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', ID='G5'/
-&DEVC XYZ= 0.1 0.0 0.02, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', ID='G6'/
-
-&DEVC XYZ=  0.200 0.0 0.385, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC1'/
-&DEVC XYZ=  0.175 0.0 0.380, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC2'/
-&DEVC XYZ=  0.175 0.0 0.340, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC3'/
-&DEVC XYZ=  0.175 0.0 0.300, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC4'/
-&DEVC XYZ=  0.175 0.0 0.260, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC5'/
-&DEVC XYZ=  0.175 0.0 0.220, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC6'/
-&DEVC XYZ=  0.175 0.0 0.180, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC7'/
-&DEVC XYZ=  0.175 0.0 0.140, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC8'/
-&DEVC XYZ=  0.175 0.0 0.100, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC9'/
-&DEVC XYZ=  0.175 0.0 0.060, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC10'/
-&DEVC XYZ=  0.175 0.0 0.020, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC11'/
-&DEVC XYZ=  0.200 0.0 0.015, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC12'/
-&DEVC XYZ=  -0.175 0.0 0.380, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC13'/
-&DEVC XYZ=  -0.175 0.0 0.260, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC14'/
-&DEVC XYZ=  -0.175 0.0 0.140, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC15'/
-&DEVC XYZ=  -0.175 0.0 0.020, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC16'/
-&DEVC XYZ=  0.000 0.0 0.380, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC17'/
-&DEVC XYZ=  0.000 0.0 0.260, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC18'/
-&DEVC XYZ=  0.000 0.0 0.140, QUANTITY='THERMOCOUPLE', RELATIVE=.TRUE., ID='TC19'/
+&RAMP ID='case1' T=350, F=0.0 /
 
-&DEVC XYZ=  0.00 0.00 0.00, IOR=3, QUANTITY= 'GAUGE HEAT FLUX', ID='HF1'/
-&DEVC XYZ= -0.15 0.00 0.00, IOR=3, QUANTITY= 'GAUGE HEAT FLUX', ID='HF2'/
-&DEVC XYZ= -0.20 0.00 0.20, IOR=1, QUANTITY= 'GAUGE HEAT FLUX', ID='HF3'/
-&DEVC XYZ=  0.00 0.00 0.40, IOR=-3, QUANTITY= 'GAUGE HEAT FLUX', ID='HF4'/
+&CATF OTHER_FILES='common_properties.txt' /
 
 &TAIL /