Merge pull request #13931 from rmcdermo/master

FDS Validation: store new input files just in case
firemodels · Dec 20, 2024 · b1a1407 · b1a1407
2 parents 2628a82 + d90053c
commit b1a1407
Show file tree

Hide file tree

Showing 7 changed files with 583 additions and 0 deletions.
diff --git a/Validation/Kashiwagi_Gasification/FDS_Input_Files/MATL.fds b/Validation/Kashiwagi_Gasification/FDS_Input_Files/MATL.fds
@@ -0,0 +1,103 @@
+&MATL ID                 = 'MOISTURE'
+      DENSITY            = 1000.
+      CONDUCTIVITY       = 0.062
+      SPECIFIC_HEAT      = 4.184
+      N_REACTIONS        = 1
+      N_S                = 1.0
+      N_T                = 0.
+      A                  = 4.29E+03
+      E                  = 4.38E+04
+      SPEC_ID            = 'WATER VAPOR'
+      NU_SPEC            = 1.0
+      HEAT_OF_REACTION   = 2410./
+
+&MATL ID                 = 'PINE'
+      FYI                = 'Properties from Lautenberger and Pello CNF 2009'
+      DENSITY            = 360.
+      CONDUCTIVITY_RAMP  = 'k dry wood'
+      SPECIFIC_HEAT_RAMP = 'c_v dry wood'
+      N_REACTIONS        = 2
+      N_S                = 1.0,1.0
+      N_T                = 0.0,0.0
+      N_O2               = 0.0,1.0
+      A                  = 2.5E+04,0
+      E                  = 0.8E+05,1.0E+05 !!! sensitive !!!
+      MATL_ID(1,1)       = 'CHAR'
+      NU_MATL(1,1)       = 0.3
+      MATL_ID(1,2)       = 'CHAR'
+      NU_MATL(1,2)       = 0.3
+      SPEC_ID(1,1)       = 'FUEL VAPOR'
+      NU_SPEC(1,1)       = 0.7
+      SPEC_ID(1:2,2)     = 'OXYGEN','FUEL VAPOR'
+      NU_SPEC(1:2,2)     = -0.1, 0.8
+      HEAT_OF_REACTION(1:2) = 600.,-900./
+
+&MATL ID                 = 'CHAR'
+      DENSITY            = 90.
+      CONDUCTIVITY_RAMP  = 'k + k_r char'
+      SPECIFIC_HEAT_RAMP = 'c_v char'
+      N_REACTIONS        = 1
+      N_S                = 1.00 ! 0.56
+      N_O2               = 1.00 ! 0.68
+      A                  = 8.0E+10
+      E                  = 1.0E+05
+      MATL_ID            = 'ASH'
+      NU_MATL            = 0.06
+      SPEC_ID            = 'OXYGEN','CARBON DIOXIDE'
+      NU_SPEC            = -2.67,3.61
+      HEAT_OF_REACTION   = -11000.
+      GAS_DIFFUSION_DEPTH(1) = 0.01/
+
+&MATL ID                 = 'ASH'
+      DENSITY            = 5.7
+      CONDUCTIVITY_RAMP  = 'k + k_r ash'
+      SPECIFIC_HEAT_RAMP = 'c_v ash'/
+
+&RAMP ID='k dry wood', T=  20.0, F=0.173549 /
+&RAMP ID='k dry wood', T= 300.0, F=0.258492 /
+&RAMP ID='k dry wood', T= 600.0, F=0.331944 /
+&RAMP ID='k dry wood', T= 900.0, F=0.395609 /
+&RAMP ID='k dry wood', T=1200.0, F=0.452913 /
+
+&RAMP ID='c_v dry wood', T=  20.0, F=1.638272 /
+&RAMP ID='c_v dry wood', T= 300.0, F=2.550559 /
+&RAMP ID='c_v dry wood', T= 600.0, F=3.367617 /
+&RAMP ID='c_v dry wood', T= 900.0, F=4.092512 /
+&RAMP ID='c_v dry wood', T=1200.0, F=4.756277 /
+
+# includes effective radiative conductivity
+&RAMP ID='k + k_r char', T=  20.0, F=0.069043 /
+&RAMP ID='k + k_r char', T= 300.0, F=0.121335 /
+&RAMP ID='k + k_r char', T= 600.0, F=0.227944 /
+&RAMP ID='k + k_r char', T= 900.0, F=0.419637 /
+&RAMP ID='k + k_r char', T=1200.0, F=0.727922 /
+
+&RAMP ID='c_v char', T=  20.0, F=1.210882 /
+&RAMP ID='c_v char', T= 300.0, F=1.463983 /
+&RAMP ID='c_v char', T= 600.0, F=1.649246 /
+&RAMP ID='c_v char', T= 900.0, F=1.793039 /
+&RAMP ID='c_v char', T=1200.0, F=1.912405 /
+
+# includes effective radiative conductivity
+&RAMP ID='k + k_r ash', T=  20.0, F=0.066647 /
+&RAMP ID='k + k_r ash', T= 300.0, F=0.141158 /
+&RAMP ID='k + k_r ash', T= 600.0, F=0.326017 /
+&RAMP ID='k + k_r ash', T= 900.0, F=0.679571 /
+&RAMP ID='k + k_r ash', T=1200.0, F=1.261558 /
+
+&RAMP ID='c_v ash', T=  20.0, F=1.234783 /
+&RAMP ID='c_v ash', T= 300.0, F=1.525267 /
+&RAMP ID='c_v ash', T= 600.0, F=1.741593 /
+&RAMP ID='c_v ash', T= 900.0, F=1.911420 /
+&RAMP ID='c_v ash', T=1200.0, F=2.053577 /
+
+
+
+
+
+
+
+
+
+
+
diff --git a/Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_25_1C_NEW.fds b/Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_25_1C_NEW.fds
@@ -0,0 +1,80 @@
+&HEAD CHID='pine_10p5O2_25_1C', TITLE='Kashiwagi et al. Combust. Flame (1987), white pine gasification, 10.5 vol% O2, 25 kW/m2'/
+
+&MESH IJK=8,8,8, XB=-0.028,0.048,-0.038,0.038,-0.038,0.038/
+
+&TIME T_END=900., WALL_INCREMENT=1/
+
+&DUMP DT_DEVC=10., DT_PROF=10., DT_RESTART=10./
+
+&MISC TMPA=30./
+
+&RADI NUMBER_RADIATION_ANGLES=16, TIME_STEP_INCREMENT=1, ANGLE_INCREMENT=1/
+
+&VENT DB='ZMIN', SURF_ID='inflow'/
+&VENT DB='ZMAX', SURF_ID='OPEN'/
+
+&SURF ID='inflow', VEL=-0.01, COLOR='BLUE'/
+&SURF ID='sample', COLOR='BROWN'/
+&SURF ID='exposed sample', CONVECTION_LENGTH_SCALE=0.019, BLOWING=T, STRETCH_FACTOR=1,
+      THICKNESS=0.038, BACKING='EXPOSED', MATL_ID(1,1:2)='PINE','MOISTURE', MATL_MASS_FRACTION(1,1:2)=0.95,0.05,
+      COLOR='GRAY 10', EXTERNAL_FLUX=25./
+&SURF ID='default wall', EMISSIVITY=0., HEAT_TRANSFER_COEFFICIENT=0., COLOR='INVISIBLE', FREE_SLIP=T, DEFAULT=T/
+
+&OBST XB=0,0.038,-0.019,0.019,-0.019,0.019, SURF_ID6='exposed sample','sample','sample','sample','sample','sample'/
+
+&SPEC ID='NITROGEN', BACKGROUND=T/
+&SPEC ID='OXYGEN', MASS_FRACTION_0=0.115/ 10.5 vol % O2 to mass fraction
+&SPEC ID='CARBON DIOXIDE'/
+&SPEC ID='FUEL VAPOR'/
+&SPEC ID='WATER VAPOR'/
+
+! simple material for setup
+! &MATL ID='PINE', DENSITY=360., SPECIFIC_HEAT=1.15, CONDUCTIVITY=0.3 /
+
+! full material model with kinetics
+&CATF OTHER_FILES='MATL.fds'/
+
+&BNDF QUANTITY='MASS FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='MASS FLUX', SPEC_ID='OXYGEN', CELL_CENTERED=T/
+&BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='BACK WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='TOTAL HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='RADIATIVE HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='CONVECTIVE HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='WALL THICKNESS', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE DENSITY', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE OXYGEN MASS FRACTION', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE OXYGEN ITERATIONS', CELL_CENTERED=T/
+&BNDF QUANTITY='HEAT TRANSFER COEFFICIENT', CELL_CENTERED=T/
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='mf', CONVERSION_FACTOR=69252, UNITS='mg/cm2 s'/
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='OXYGEN',         SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='norm mf O2',  CONVERSION_FACTOR=-25937, UNITS='mg/cm2 s'/ 69252/(-2.67)
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='CARBON DIOXIDE', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='norm mf CO2', CONVERSION_FACTOR= 18973, UNITS='mg/cm2 s'/ 69252/( 3.65)
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='OXYGEN',         SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='NORM MASS O2 CONSUMED',  CONVERSION_FACTOR=-0.374531835205993/ -1/2.67
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='CARBON DIOXIDE', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='NORM MASS CO2 PRODUCED', CONVERSION_FACTOR= 0.273972602739726/  1/3.65
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='SURFACE DENSITY', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='MASS', SMOOTHING_TIME=10/
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='SURFACE DENSITY', MATL_ID='ASH', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='INSTANT VALUE', ID='NORM MASS ASH PRODUCED', CONVERSION_FACTOR= 50./ 1/0.02
+
+&CTRL ID='PID WT LOSS', FUNCTION_TYPE='PID', INPUT_ID='MASS', PROPORTIONAL_GAIN=0, DIFFERENTIAL_GAIN=1, INTEGRAL_GAIN=0/
+&DEVC XYZ=0,0,0, ID='wt loss mf', QUANTITY='CONTROL VALUE', CTRL_ID='PID WT LOSS', CONVERSION_FACTOR=-69252, UNITS='mg/cm2 s'/
+
+&DEVC XYZ=-0.001,0,0, QUANTITY='WALL TEMPERATURE', IOR=-1, ID='T0mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.005, IOR=-1, ID='T5mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.010, IOR=-1, ID='T10mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.015, IOR=-1, ID='T15mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='BACK WALL TEMPERATURE', IOR=-1, ID='T38mm'/
+
+&PROF ID='Tx' , XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', IOR=-1 /
+&PROF ID='MOISTUREx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='MOISTURE', IOR=-1 /
+&PROF ID='PINEx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='PINE', IOR=-1 /
+&PROF ID='CHARx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='CHAR', IOR=-1 /
+&PROF ID='ASHx' , XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='ASH' , IOR=-1 /
+&PROF ID='O2x'  , XYZ=-0.001,0,0, QUANTITY='SOLID OXYGEN MASS FRACTION', MATL_ID='CHAR' , IOR=-1 /
+
+&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
+&SLCF PBY=0, QUANTITY='TEMPERATURE', CELL_CENTERED=T/
+&SLCF PBY=0, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=T/
+
+&TAIL /
diff --git a/Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_40_1C_NEW.fds b/Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_40_1C_NEW.fds
@@ -0,0 +1,80 @@
+&HEAD CHID='pine_10p5O2_40_1C', TITLE='Kashiwagi et al. Combust. Flame (1987), white pine gasification, 10.5 vol% O2, 40 kW/m2'/
+
+&MESH IJK=8,8,8, XB=-0.028,0.048,-0.038,0.038,-0.038,0.038/
+
+&TIME T_END=900., WALL_INCREMENT=1/
+
+&DUMP DT_DEVC=10., DT_PROF=10., DT_RESTART=10./
+
+&MISC TMPA=30./
+
+&RADI NUMBER_RADIATION_ANGLES=16, TIME_STEP_INCREMENT=1, ANGLE_INCREMENT=1/
+
+&VENT DB='ZMIN', SURF_ID='inflow'/
+&VENT DB='ZMAX', SURF_ID='OPEN'/
+
+&SURF ID='inflow', VEL=-0.01, COLOR='BLUE'/
+&SURF ID='sample', COLOR='BROWN'/
+&SURF ID='exposed sample', CONVECTION_LENGTH_SCALE=0.019, BLOWING=T, STRETCH_FACTOR=1,
+      THICKNESS=0.038, BACKING='EXPOSED', MATL_ID(1,1:2)='PINE','MOISTURE', MATL_MASS_FRACTION(1,1:2)=0.95,0.05,
+      COLOR='GRAY 10', EXTERNAL_FLUX=40./
+&SURF ID='default wall', EMISSIVITY=0., HEAT_TRANSFER_COEFFICIENT=0., COLOR='INVISIBLE', FREE_SLIP=T, DEFAULT=T/
+
+&OBST XB=0,0.038,-0.019,0.019,-0.019,0.019, SURF_ID6='exposed sample','sample','sample','sample','sample','sample'/
+
+&SPEC ID='NITROGEN', BACKGROUND=T/
+&SPEC ID='OXYGEN', MASS_FRACTION_0=0.115/ 10.5 vol % O2 to mass fraction
+&SPEC ID='CARBON DIOXIDE'/
+&SPEC ID='FUEL VAPOR'/
+&SPEC ID='WATER VAPOR'/
+
+! simple material for setup
+! &MATL ID='PINE', DENSITY=360., SPECIFIC_HEAT=1.15, CONDUCTIVITY=0.3 /
+
+! full material model with kinetics
+&CATF OTHER_FILES='MATL.fds'/
+
+&BNDF QUANTITY='MASS FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='MASS FLUX', SPEC_ID='OXYGEN', CELL_CENTERED=T/
+&BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='BACK WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='TOTAL HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='RADIATIVE HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='CONVECTIVE HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='WALL THICKNESS', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE DENSITY', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE OXYGEN MASS FRACTION', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE OXYGEN ITERATIONS', CELL_CENTERED=T/
+&BNDF QUANTITY='HEAT TRANSFER COEFFICIENT', CELL_CENTERED=T/
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='mf', CONVERSION_FACTOR=69252, UNITS='mg/cm2 s'/
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='OXYGEN',         SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='norm mf O2',  CONVERSION_FACTOR=-25937, UNITS='mg/cm2 s'/ 69252/(-2.67)
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='CARBON DIOXIDE', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='norm mf CO2', CONVERSION_FACTOR= 18973, UNITS='mg/cm2 s'/ 69252/( 3.65)
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='OXYGEN',         SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='NORM MASS O2 CONSUMED',  CONVERSION_FACTOR=-0.374531835205993/ -1/2.67
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='CARBON DIOXIDE', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='NORM MASS CO2 PRODUCED', CONVERSION_FACTOR= 0.273972602739726/  1/3.65
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='SURFACE DENSITY', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='MASS', SMOOTHING_TIME=10/
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='SURFACE DENSITY', MATL_ID='ASH', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='INSTANT VALUE', ID='NORM MASS ASH PRODUCED', CONVERSION_FACTOR= 50./ 1/0.02
+
+&CTRL ID='PID WT LOSS', FUNCTION_TYPE='PID', INPUT_ID='MASS', PROPORTIONAL_GAIN=0, DIFFERENTIAL_GAIN=1, INTEGRAL_GAIN=0/
+&DEVC XYZ=0,0,0, ID='wt loss mf', QUANTITY='CONTROL VALUE', CTRL_ID='PID WT LOSS', CONVERSION_FACTOR=-69252, UNITS='mg/cm2 s'/
+
+&DEVC XYZ=-0.001,0,0, QUANTITY='WALL TEMPERATURE', IOR=-1, ID='T0mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.005, IOR=-1, ID='T5mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.010, IOR=-1, ID='T10mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.015, IOR=-1, ID='T15mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='BACK WALL TEMPERATURE', IOR=-1, ID='T38mm'/
+
+&PROF ID='Tx' , XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', IOR=-1 /
+&PROF ID='MOISTUREx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='MOISTURE', IOR=-1 /
+&PROF ID='PINEx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='PINE', IOR=-1 /
+&PROF ID='CHARx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='CHAR', IOR=-1 /
+&PROF ID='ASHx' , XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='ASH' , IOR=-1 /
+&PROF ID='O2x'  , XYZ=-0.001,0,0, QUANTITY='SOLID OXYGEN MASS FRACTION', MATL_ID='CHAR' , IOR=-1 /
+
+&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
+&SLCF PBY=0, QUANTITY='TEMPERATURE', CELL_CENTERED=T/
+&SLCF PBY=0, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=T/
+
+&TAIL /
diff --git a/Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_21O2_25_1C_NEW.fds b/Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_21O2_25_1C_NEW.fds
@@ -0,0 +1,80 @@
+&HEAD CHID='pine_21O2_25_1C', TITLE='Kashiwagi et al. Combust. Flame (1987), white pine gasification, 21 vol% O2, 25 kW/m2'/
+
+&MESH IJK=8,8,8, XB=-0.028,0.048,-0.038,0.038,-0.038,0.038/
+
+&TIME T_END=900., WALL_INCREMENT=1/
+
+&DUMP DT_DEVC=10., DT_PROF=10., DT_RESTART=10./
+
+&MISC TMPA=30./
+
+&RADI NUMBER_RADIATION_ANGLES=16, TIME_STEP_INCREMENT=1, ANGLE_INCREMENT=1/
+
+&VENT DB='ZMIN', SURF_ID='inflow'/
+&VENT DB='ZMAX', SURF_ID='OPEN'/
+
+&SURF ID='inflow', VEL=-0.01, COLOR='BLUE'/
+&SURF ID='sample', COLOR='BROWN'/
+&SURF ID='exposed sample', CONVECTION_LENGTH_SCALE=0.019, BLOWING=T, STRETCH_FACTOR=1,
+      THICKNESS=0.038, BACKING='EXPOSED', MATL_ID(1,1:2)='PINE','MOISTURE', MATL_MASS_FRACTION(1,1:2)=0.95,0.05,
+      COLOR='GRAY 10', EXTERNAL_FLUX=25./
+&SURF ID='default wall', EMISSIVITY=0., HEAT_TRANSFER_COEFFICIENT=0., COLOR='INVISIBLE', FREE_SLIP=T, DEFAULT=T/
+
+&OBST XB=0,0.038,-0.019,0.019,-0.019,0.019, SURF_ID6='exposed sample','sample','sample','sample','sample','sample'/
+
+&SPEC ID='NITROGEN', BACKGROUND=T/
+&SPEC ID='OXYGEN', MASS_FRACTION_0=0.232/ 21 vol % O2 to mass fraction
+&SPEC ID='CARBON DIOXIDE'/
+&SPEC ID='FUEL VAPOR'/
+&SPEC ID='WATER VAPOR'/
+
+! simple material for setup
+! &MATL ID='PINE', DENSITY=360., SPECIFIC_HEAT=1.15, CONDUCTIVITY=0.3 /
+
+! full material model with kinetics
+&CATF OTHER_FILES='MATL.fds'/
+
+&BNDF QUANTITY='MASS FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='MASS FLUX', SPEC_ID='OXYGEN', CELL_CENTERED=T/
+&BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='BACK WALL TEMPERATURE', CELL_CENTERED=T/
+&BNDF QUANTITY='TOTAL HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='RADIATIVE HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='CONVECTIVE HEAT FLUX', CELL_CENTERED=T/
+&BNDF QUANTITY='WALL THICKNESS', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE DENSITY', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE OXYGEN MASS FRACTION', CELL_CENTERED=T/
+&BNDF QUANTITY='SURFACE OXYGEN ITERATIONS', CELL_CENTERED=T/
+&BNDF QUANTITY='HEAT TRANSFER COEFFICIENT', CELL_CENTERED=T/
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='mf', CONVERSION_FACTOR=69252, UNITS='mg/cm2 s'/
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='OXYGEN',         SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='norm mf O2',  CONVERSION_FACTOR=-25937, UNITS='mg/cm2 s'/ 69252/(-2.67)
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='CARBON DIOXIDE', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='norm mf CO2', CONVERSION_FACTOR= 18973, UNITS='mg/cm2 s'/ 69252/( 3.65)
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='OXYGEN',         SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='NORM MASS O2 CONSUMED',  CONVERSION_FACTOR=-0.374531835205993/ -1/2.67
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='MASS FLUX', SPEC_ID='CARBON DIOXIDE', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='TIME INTEGRAL', ID='NORM MASS CO2 PRODUCED', CONVERSION_FACTOR= 0.273972602739726/  1/3.65
+
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='SURFACE DENSITY', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='MASS', SMOOTHING_TIME=10/
+&DEVC XB=-0.028,0.048,-0.038,0.038,-0.038,0.038, QUANTITY='SURFACE DENSITY', MATL_ID='ASH', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', TEMPORAL_STATISTIC='INSTANT VALUE', ID='NORM MASS ASH PRODUCED', CONVERSION_FACTOR= 50./ 1/0.02
+
+&CTRL ID='PID WT LOSS', FUNCTION_TYPE='PID', INPUT_ID='MASS', PROPORTIONAL_GAIN=0, DIFFERENTIAL_GAIN=1, INTEGRAL_GAIN=0/
+&DEVC XYZ=0,0,0, ID='wt loss mf', QUANTITY='CONTROL VALUE', CTRL_ID='PID WT LOSS', CONVERSION_FACTOR=-69252, UNITS='mg/cm2 s'/
+
+&DEVC XYZ=-0.001,0,0, QUANTITY='WALL TEMPERATURE', IOR=-1, ID='T0mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.005, IOR=-1, ID='T5mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.010, IOR=-1, ID='T10mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', DEPTH=0.015, IOR=-1, ID='T15mm'/
+&DEVC XYZ=-0.001,0,0, QUANTITY='BACK WALL TEMPERATURE', IOR=-1, ID='T38mm'/
+
+&PROF ID='Tx' , XYZ=-0.001,0,0, QUANTITY='INSIDE WALL TEMPERATURE', IOR=-1 /
+&PROF ID='MOISTUREx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='MOISTURE', IOR=-1 /
+&PROF ID='PINEx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='PINE', IOR=-1 /
+&PROF ID='CHARx', XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='CHAR', IOR=-1 /
+&PROF ID='ASHx' , XYZ=-0.001,0,0, QUANTITY='SOLID DENSITY', MATL_ID='ASH' , IOR=-1 /
+&PROF ID='O2x'  , XYZ=-0.001,0,0, QUANTITY='SOLID OXYGEN MASS FRACTION', MATL_ID='CHAR' , IOR=-1 /
+
+&SLCF PBY=0, QUANTITY='VELOCITY', VECTOR=T/
+&SLCF PBY=0, QUANTITY='TEMPERATURE', CELL_CENTERED=T/
+&SLCF PBY=0, QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=T/
+
+&TAIL /