Skip to content

Commit

Permalink
Merge pull request firemodels#14013 from rmcdermo/master
Browse files Browse the repository at this point in the history
FDS Validation: updates to Kashiwagi Gasification input files
rmcdermo authored Jan 8, 2025
2 parents 3951a31 + 649e275 commit 01501b4
Showing 20 changed files with 200 additions and 1,062 deletions.
91 changes: 56 additions & 35 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/MATL.fds
Original file line number Diff line number Diff line change
@@ -12,58 +12,67 @@
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 !!!
N_O2 = 0.0,0.7
A = 1.0E+06,1.1E+05
E = 1.05E+05,8.8E+04
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
MATL_ID(1,2) = 'CHAR'
NU_MATL(1,2) = 0.3
SPEC_ID(1:2,2) = 'OXYGEN','FUEL VAPOR'
NU_SPEC(1:2,2) = -0.1, 0.8
HEAT_OF_REACTION(1:2) = 600.,-900./
NU_SPEC(1:2,2) = -0.5,1.2
HEAT_OF_REACTION(1:2) = 200.,-5700.
ADJUST_H = F/

&MATL ID = 'CHAR'
DENSITY = 90.
DENSITY = 108.
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/
N_S = 1
N_O2 = 0.8
A = 5.0e+08
E = 1.4E+05
MATL_ID(1,1) = 'ASH'
NU_MATL(1,1) = 0.01
SPEC_ID(1:2,1) = 'OXYGEN','CARBON DIOXIDE'
NU_SPEC(1:2,1) = -2.67,3.66
HEAT_OF_REACTION(1) = -20000.
ADJUST_H = F
GAS_DIFFUSION_DEPTH = 0.001/

&MATL ID = 'ASH'
DENSITY = 5.7
DENSITY = 5.0
CONDUCTIVITY_RAMP = 'k + k_r ash'
SPECIFIC_HEAT_RAMP = 'c_v ash'/

! from Lautenberger
&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 /
! &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 /

! 1.15 from Kashiwagi paper
&RAMP ID='c_v dry wood', T= 20.0, F=2.5 /
&RAMP ID='c_v dry wood', T= 300.0, F=2.5 /
&RAMP ID='c_v dry wood', T= 600.0, F=2.5 /
&RAMP ID='c_v dry wood', T= 900.0, F=2.5 /
&RAMP ID='c_v dry wood', T=1200.0, F=2.5 /

# includes effective radiative conductivity
&RAMP ID='k + k_r char', T= 20.0, F=0.069043 /
@@ -72,11 +81,17 @@
&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 /
! &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 /

&RAMP ID='c_v char', T= 20.0, F=1.5 /
&RAMP ID='c_v char', T= 300.0, F=1.5 /
&RAMP ID='c_v char', T= 600.0, F=1.5 /
&RAMP ID='c_v char', T= 900.0, F=1.5 /
&RAMP ID='c_v char', T=1200.0, F=1.5 /

# includes effective radiative conductivity
&RAMP ID='k + k_r ash', T= 20.0, F=0.066647 /
@@ -85,11 +100,17 @@
&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 /
! &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 /

&RAMP ID='c_v ash', T= 20.0, F=1.5 /
&RAMP ID='c_v ash', T= 300.0, F=1.5 /
&RAMP ID='c_v ash', T= 600.0, F=1.5 /
&RAMP ID='c_v ash', T= 900.0, F=1.5 /
&RAMP ID='c_v ash', T=1200.0, F=1.5 /



Original file line number Diff line number Diff line change
@@ -1,27 +1,30 @@
&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=16,16,16, XB=-0.03,0.05,-0.04,0.04,-0.04,0.04/
&MESH IJK=16,16,16, 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./

&MISC TEST_CHAR_MASS_TRANSFER_MODEL=T/
&MISC TMPA=30., TEST_NEW_CHAR_MODEL=T, BNDF_DEFAULT=F/

&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.005, COLOR='BLUE'/
&SURF ID='sample', STRETCH_FACTOR=1, THICKNESS=0.038, BACKING='EXPOSED', MATL_ID='PINE', MOISTURE_FRACTION=0.0526, COLOR='BROWN'/
&SURF ID='exposed sample', STRETCH_FACTOR=1, THICKNESS=0.038, BACKING='EXPOSED', MATL_ID='PINE', MOISTURE_FRACTION=0.0526, COLOR='GRAY 10', EXTERNAL_FLUX=25./
&SURF ID='sample', COLOR='BROWN'/
&SURF ID='exposed sample', CONVECTION_LENGTH_SCALE=0.038, 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'/
&OBST XB=0,0.038,-0.019,0.019,-0.019,0.019, SURF_ID6='exposed sample','sample','sample','sample','sample','sample', BNDF_OBST=T/

&SPEC ID='NITROGEN', BACKGROUND=T/
&SPEC ID='OXYGEN', MASS_FRACTION_0=0.115/ 10.5 vol % O2 to mass fraction
&SPEC ID='CARBON MONOXIDE'/
&SPEC ID='CARBON DIOXIDE'/
&SPEC ID='FUEL VAPOR'/
&SPEC ID='WATER VAPOR'/
@@ -30,27 +33,31 @@
! &MATL ID='PINE', DENSITY=360., SPECIFIC_HEAT=1.15, CONDUCTIVITY=0.3 /

! full material model with kinetics
&CATF OTHER_FILES='../../../Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds'/
&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='SUBSTEPS', 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/
&BNDF QUANTITY='BLOWING CORRECTION', CELL_CENTERED=T/

&DEVC XB=-0.03,0.05,-0.04,0.04,-0.04,0.04, QUANTITY='MASS FLUX', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='mf', CONVERSION_FACTOR=69252, UNITS='mg/cm2 s'/
&DEVC XB=-0.03,0.05,-0.04,0.04,-0.04,0.04, 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.03,0.05,-0.04,0.04,-0.04,0.04, 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', 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.03,0.05,-0.04,0.04,-0.04,0.04, 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.03,0.05,-0.04,0.04,-0.04,0.04, 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='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.03,0.05,-0.04,0.04,-0.04,0.04, QUANTITY='SURFACE DENSITY', SURF_ID='exposed sample', SPATIAL_STATISTIC='SURFACE INTEGRAL', ID='MASS', SMOOTHING_TIME=10/
&DEVC XB=-0.03,0.05,-0.04,0.04,-0.04,0.04, 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
&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'/
@@ -62,9 +69,11 @@
&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/

This file was deleted.

Loading

0 comments on commit 01501b4

Please sign in to comment.