Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

FDS Validation: store new input files just in case #13931

Merged
merged 1 commit into from
Dec 20, 2024
Merged

FDS Validation: store new input files just in case #13931

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
103 changes: 103 additions & 0 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/MATL.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 /











80 changes: 80 additions & 0 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_25_1C_NEW.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 /
80 changes: 80 additions & 0 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_40_1C_NEW.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 /
80 changes: 80 additions & 0 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_21O2_25_1C_NEW.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -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 /
Loading
Loading