Skip to content

FDS Validation: updates to Kashiwagi Gasification input files #14013

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

Merged
rmcdermo merged 3 commits into from
Jan 8, 2025
Merged

FDS Validation: updates to Kashiwagi Gasification input files #14013

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
c8472ac
FDS Validation: remove old Kashiwagi input files
rmcdermo Jan 8, 2025
acb9239
FDS Validation: update Kashiwagi input files
rmcdermo Jan 8, 2025
649e275
FDS Validation: rename new input files
rmcdermo Jan 8, 2025
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
91 changes: 56 additions & 35 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/MATL.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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 /
Expand All @@ -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 /
Expand All @@ -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 /



Expand Down
37 changes: 23 additions & 14 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_25_1C.fds
View file
Open in desktop
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'/
Expand All @@ -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'/
Expand All @@ -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/
Expand Down
80 changes: 0 additions & 80 deletions Validation/Kashiwagi_Gasification/FDS_Input_Files/pine_10p5O2_25_1C_NEW.fds
View file
Open in desktop

This file was deleted.

Loading
Loading