diff --git a/Utilities/Input_Libraries/MATL/cellulose_1C_MATL.fds b/Utilities/Input_Libraries/MATL/cellulose_1C_MATL.fds new file mode 100644 index 00000000000..a5607a48a3d --- /dev/null +++ b/Utilities/Input_Libraries/MATL/cellulose_1C_MATL.fds @@ -0,0 +1,107 @@ +&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 = 'CELLULOSE' + FYI = 'Kashiwagi and Nambu CNF 1992, Table 1' + DENSITY = 360. + CONDUCTIVITY_RAMP = 'k dry wood' + SPECIFIC_HEAT_RAMP = 'c_v dry wood' + N_REACTIONS = 2 + N_S = 1.8,1.3 + N_T = 0,0 + N_O2 = 0.0,0.5 + A = 1.803e+15,4.276e+11 + E = 2.2E+05,1.6E+05 + MATL_ID(1,1) = 'CHAR' + NU_MATL(1,1) = 0.24 + SPEC_ID(1:4,1) = 'WATER VAPOR','CARBON MONOXIDE','CARBON DIOXIDE','FUEL VAPOR' + NU_SPEC(1:4,1) = 0.5,0.01,0.03,0.22 + MATL_ID(1,2) = 'CHAR' + NU_MATL(1,2) = 0.21 + SPEC_ID(1:5,2) = 'OXYGEN','WATER VAPOR','CARBON MONOXIDE','CARBON DIOXIDE','FUEL VAPOR' + NU_SPEC(1:5,2) = -0.41,0.8,0.08,0.24,0.08 + HEAT_OF_REACTION(1:2) = 570.,-5700. + ADJUST_H = F/ + +&MATL ID = 'CHAR' + DENSITY = 75. + CONDUCTIVITY_RAMP = 'k + k_r char' + SPECIFIC_HEAT_RAMP = 'c_v char' + N_REACTIONS = 1 + N_S = 1 + N_O2 = 0.78 + A = 5.667e+09 + E = 1.6E+05 + MATL_ID(1,1) = 'ASH' + NU_MATL(1,1) = 0.03 + SPEC_ID(1:4,1) = 'OXYGEN','WATER VAPOR','CARBON MONOXIDE','CARBON DIOXIDE' + NU_SPEC(1:4,1) = -1.65,0.3,0.5,1.82 + HEAT_OF_REACTION(1) = -25000. + ADJUST_H = F + GAS_DIFFUSION_DEPTH = 0.001/ + +&MATL ID = 'ASH' + 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 / + +! 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 / +&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.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 / +&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.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 / + + + + + + + + + + + diff --git a/Verification/FDS_Cases.sh b/Verification/FDS_Cases.sh index 89e14ca63db..88a29c5601d 100755 --- a/Verification/FDS_Cases.sh +++ b/Verification/FDS_Cases.sh @@ -424,6 +424,9 @@ $QFDS -d Pyrolysis cable_11_jacket_mcc.fds $QFDS -d Pyrolysis cable_23_jacket_mcc.fds $QFDS -d Pyrolysis cable_701_jacket_mcc.fds $QFDS -d Pyrolysis cell_burn_away.fds +$QFDS -d Pyrolysis cellulose_TGA_1KPM_Air.fds +$QFDS -d Pyrolysis cellulose_TGA_3KPM_Air.fds +$QFDS -d Pyrolysis cellulose_TGA_5KPM_Air.fds $QFDS -d Pyrolysis birch_tga_1step_2.fds $QFDS -d Pyrolysis birch_tga_1step_20.fds $QFDS -d Pyrolysis enthalpy.fds diff --git a/Verification/Pyrolysis/cellulose_TGA_1KPM_Air.fds b/Verification/Pyrolysis/cellulose_TGA_1KPM_Air.fds new file mode 100644 index 00000000000..73fdfee49be --- /dev/null +++ b/Verification/Pyrolysis/cellulose_TGA_1KPM_Air.fds @@ -0,0 +1,27 @@ +&HEAD CHID='cellulose_TGA_1KPM_Air', TITLE='Kashiwagi and Nambu CNF 1992'/ + +&MESH IJK=10,10,10,XB=-.5,.5,-.5,.5,0,1/ + +&TIME T_END=1. / + +&SPEC ID='NITROGEN', BACKGROUND=T/ +&SPEC ID='OXYGEN', MASS_FRACTION_0=0.23/ +&SPEC ID='WATER VAPOR'/ +&SPEC ID='CARBON MONOXIDE'/ +&SPEC ID='CARBON DIOXIDE'/ +&SPEC ID='FUEL VAPOR'/ + +&SURF ID = 'FUEL' + MATL_ID = 'CELLULOSE' + SURFACE_VOLUME_RATIO = 7590 + LENGTH = 0.1 + GEOMETRY = 'CYLINDRICAL' + TGA_ANALYSIS = T + TGA_HEATING_RATE = 1./ + +&PART ID='FUEL', SAMPLING_FACTOR=1, SURF_ID='FUEL' / +&INIT PART_ID='FUEL', XYZ=0.0,0.0,0.3, N_PARTICLES=1/ + +&CATF OTHER_FILES='../../Utilities/Input_Libraries/MATL/cellulose_1C_MATL.fds'/ + +&TAIL / diff --git a/Verification/Pyrolysis/cellulose_TGA_3KPM_Air.fds b/Verification/Pyrolysis/cellulose_TGA_3KPM_Air.fds new file mode 100644 index 00000000000..e80f1ad3bae --- /dev/null +++ b/Verification/Pyrolysis/cellulose_TGA_3KPM_Air.fds @@ -0,0 +1,27 @@ +&HEAD CHID='cellulose_TGA_3KPM_Air', TITLE='Kashiwagi and Nambu CNF 1992'/ + +&MESH IJK=10,10,10,XB=-.5,.5,-.5,.5,0,1/ + +&TIME T_END=1. / + +&SPEC ID='NITROGEN', BACKGROUND=T/ +&SPEC ID='OXYGEN', MASS_FRACTION_0=0.23/ +&SPEC ID='WATER VAPOR'/ +&SPEC ID='CARBON MONOXIDE'/ +&SPEC ID='CARBON DIOXIDE'/ +&SPEC ID='FUEL VAPOR'/ + +&SURF ID = 'FUEL' + MATL_ID = 'CELLULOSE' + SURFACE_VOLUME_RATIO = 7590 + LENGTH = 0.1 + GEOMETRY = 'CYLINDRICAL' + TGA_ANALYSIS = T + TGA_HEATING_RATE = 3./ + +&PART ID='FUEL', SAMPLING_FACTOR=1, SURF_ID='FUEL' / +&INIT PART_ID='FUEL', XYZ=0.0,0.0,0.3, N_PARTICLES=1/ + +&CATF OTHER_FILES='../../Utilities/Input_Libraries/MATL/cellulose_1C_MATL.fds'/ + +&TAIL / diff --git a/Verification/Pyrolysis/cellulose_TGA_5KPM_Air.fds b/Verification/Pyrolysis/cellulose_TGA_5KPM_Air.fds new file mode 100644 index 00000000000..f2969cb93b2 --- /dev/null +++ b/Verification/Pyrolysis/cellulose_TGA_5KPM_Air.fds @@ -0,0 +1,27 @@ +&HEAD CHID='cellulose_TGA_5KPM_Air', TITLE='Kashiwagi and Nambu CNF 1992'/ + +&MESH IJK=10,10,10,XB=-.5,.5,-.5,.5,0,1/ + +&TIME T_END=1. / + +&SPEC ID='NITROGEN', BACKGROUND=T/ +&SPEC ID='OXYGEN', MASS_FRACTION_0=0.23/ +&SPEC ID='WATER VAPOR'/ +&SPEC ID='CARBON MONOXIDE'/ +&SPEC ID='CARBON DIOXIDE'/ +&SPEC ID='FUEL VAPOR'/ + +&SURF ID = 'FUEL' + MATL_ID = 'CELLULOSE' + SURFACE_VOLUME_RATIO = 7590 + LENGTH = 0.1 + GEOMETRY = 'CYLINDRICAL' + TGA_ANALYSIS = T + TGA_HEATING_RATE = 5./ + +&PART ID='FUEL', SAMPLING_FACTOR=1, SURF_ID='FUEL' / +&INIT PART_ID='FUEL', XYZ=0.0,0.0,0.3, N_PARTICLES=1/ + +&CATF OTHER_FILES='../../Utilities/Input_Libraries/MATL/cellulose_1C_MATL.fds'/ + +&TAIL /