firemodels · rmcdermo · Dec 23, 2024 · Dec 23, 2024
diff --git a/Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds b/Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds
@@ -5,10 +5,10 @@
       CONDUCTIVITY_RAMP  = 'k dry wood'
       SPECIFIC_HEAT_RAMP = 'c_v dry wood'
       N_REACTIONS        = 2
-      N_S                = 0.87,0.63
-      N_O2               = 0.0,0.72
-      A                  = 4.70E+6,1.45E+10
-      E                  = 1.05E+05,1.27E+05
+      N_S                = 0.50,0.63
+      N_O2               = 0.00,0.72
+      A                  = 1.00E+07,1.00E+13
+      E                  = 1.00E+05,1.58E+05
       MATL_ID(1,1)       = 'CHAR'
       NU_MATL(1,1)       = 0.31
       MATL_ID(1,2)       = 'CHAR'
@@ -17,25 +17,25 @@
       NU_SPEC(1,1)       = 0.69
       SPEC_ID(1:2,2)     = 'OXYGEN','FUEL VAPOR'
       NU_SPEC(1:2,2)     = -0.1, 0.79
-      HEAT_OF_REACTION   = 416. /
+      HEAT_OF_REACTION(1:2) = 416./
 
 &MATL ID                 = 'CHAR'
-      DENSITY            = 73.
+      DENSITY            = 112.     ! 112/360 = 0.31 to match yield from AC Tab 7
       CONDUCTIVITY_RAMP  = 'k + k_r char'
       SPECIFIC_HEAT_RAMP = 'c_v char'
       N_REACTIONS        = 1
       N_S                = 0.56
       N_O2               = 0.68
-      A                  = 8.31E+07
-      E                  = 1.24E+05
+      A                  = 7.00E+06
+      E                  = 1.09E+05
       MATL_ID            = 'ASH'
       NU_MATL            = 0.02
       SPEC_ID            = 'OXYGEN','CARBON DIOXIDE'
       NU_SPEC            = -2.67,3.65
-      HEAT_OF_REACTION   = -32000. /
+      HEAT_OF_REACTION   = -32000./
 
 &MATL ID                 = 'ASH'
-      DENSITY            = 5.7
+      DENSITY            = 2.0
       CONDUCTIVITY_RAMP  = 'k + k_r ash'
       SPECIFIC_HEAT_RAMP = 'c_v ash'/
 

diff --git a/Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds b/Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds
@@ -1,13 +1,13 @@
 
-&MATL ID                 = 'PINE 1'
+&MATL ID                 = 'PINE 1' ! peak of anaerobic curve (T_ref approx 360 C)
       DENSITY            = 360.
       CONDUCTIVITY_RAMP  = 'k dry wood'
       SPECIFIC_HEAT_RAMP = 'c_v dry wood'
       N_REACTIONS        = 2
-      N_S                = 0.56,0.3
-      N_O2               = 0.0,0.61
-      A                  = 4.9E+09,8.9E+09
-      E                  = 1.46E+05,1.16E+05
+      N_S                = 0.56,0.50
+      N_O2               = 0.00,0.61
+      A                  = 1.00E+11,1.00E+16
+      E                  = 1.50E+05,1.89E+05
       MATL_ID(1,1)       = 'CHAR'
       NU_MATL(1,1)       = 0.25
       MATL_ID(1,2)       = 'CHAR'
@@ -16,17 +16,17 @@
       NU_SPEC(1,1)       = 0.75
       SPEC_ID(1:2,2)     = 'OXYGEN','FUEL VAPOR'
       NU_SPEC(1:2,2)     = -0.1, 0.85
-      HEAT_OF_REACTION   = 416. /
+      HEAT_OF_REACTION(1:2) = 416./
 
-&MATL ID                 = 'PINE 2'
+&MATL ID                 = 'PINE 2' ! first hump of anaerobic curve
       DENSITY            = 360.
       CONDUCTIVITY_RAMP  = 'k dry wood'
       SPECIFIC_HEAT_RAMP = 'c_v dry wood'
       N_REACTIONS        = 2
       N_S                = 1,1
-      N_O2               = 0.0,0.49
-      A                  = 5.0E+10,2.0E+04
-      E                  = 1.44E+05,7.5E+04
+      N_O2               = 0.00,0.49
+      A                  = 1.00E+15,1.00E+05
+      E                  = 1.93E+05,8.50E+04
       MATL_ID(1,1)       = 'CHAR'
       NU_MATL(1,1)       = 0.25
       MATL_ID(1,2)       = 'CHAR'
@@ -35,17 +35,17 @@
       NU_SPEC(1,1)       = 0.75
       SPEC_ID(1:2,2)     = 'OXYGEN','FUEL VAPOR'
       NU_SPEC(1:2,2)     = -0.1, 0.85
-      HEAT_OF_REACTION   = 416. /
+      HEAT_OF_REACTION(1:2) = 416./ AC Tab 8
 
-&MATL ID                 = 'PINE 3'
+&MATL ID                 = 'PINE 3' ! broad third peak (tail of anaerobic curve)
       DENSITY            = 360.
       CONDUCTIVITY_RAMP  = 'k dry wood'
       SPECIFIC_HEAT_RAMP = 'c_v dry wood'
       N_REACTIONS        = 2
       N_S                = 1.25,5.67
-      N_O2               = 0.0,0.66
-      A                  = 2.9E+11,2.6
-      E                  = 1.64E+05,1.64E+05
+      N_O2               = 0.00,0.66
+      A                  = 1.00,0.00
+      E                  = 4.00E+04,1.67E+05
       MATL_ID(1,1)       = 'CHAR'
       NU_MATL(1,1)       = 0.25
       MATL_ID(1,2)       = 'CHAR'
@@ -54,16 +54,16 @@
       NU_SPEC(1,1)       = 0.75
       SPEC_ID(1:2,2)     = 'OXYGEN','FUEL VAPOR'
       NU_SPEC(1:2,2)     = -0.1, 0.85
-      HEAT_OF_REACTION   = 416. /
+      HEAT_OF_REACTION(1:2) = 416./
 
 &MATL ID                 = 'CHAR'
-      DENSITY            = 73.
+      DENSITY            = 90.
       CONDUCTIVITY_RAMP  = 'k + k_r char'
       SPECIFIC_HEAT_RAMP = 'c_v char'
       N_REACTIONS        = 1
       N_S                = 0.56
       N_O2               = 0.68
-      A                  = 7.55E+07
+      A                  = 8.31E+07
       E                  = 1.24E+05
       MATL_ID            = 'ASH'
       NU_MATL            = 0.02
@@ -72,7 +72,7 @@
       HEAT_OF_REACTION   = -32000. /
 
 &MATL ID                 = 'ASH'
-      DENSITY            = 5.7
+      DENSITY            = 2.0
       CONDUCTIVITY_RAMP  = 'k + k_r ash'
       SPECIFIC_HEAT_RAMP = 'c_v ash'/
 

diff --git a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
@@ -728,6 +728,10 @@ f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_1C_cat_git.txt,Pyrolysis/pine_wood
 d,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 13,Exp (2.5 K/min),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp13_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (2.5 K/min),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,20.5% O_2 (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,3.20E-03,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_rate,N/A,end,0,pine wood TGA,kd,k,TeX
 f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 14,Exp (5 K/min),r^,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp14_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (5 K/min),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,20.5% O_2 (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,3.20E-03,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_rate,N/A,end,0,pine wood TGA,kd,k,TeX
 f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 15,Exp (10 K/min),b*,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp15_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (10 K/min),b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,20.5% O_2 (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,3.20E-03,1,no,0.05 0.90,East,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_rate,N/A,end,0,pine wood TGA,kd,k,TeX
+d,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 14,Exp (20.5 % O_2),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp14_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (20.5 % O_2),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,5 K/min (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,2.00E-03,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_oxygen,N/A,end,0,pine wood TGA,kd,k,TeX
+f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 16,Exp (8.2 % O_2),r^,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp16_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (8.2 % O_2),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,5 K/min (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,2.00E-03,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_oxygen,N/A,end,0,pine wood TGA,kd,k,TeX
+f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 17,Exp (4.3 % O_2),b*,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp17_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (4.3 % O_2),b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,5 K/min (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,2.00E-03,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_oxygen,N/A,end,0,pine wood TGA,kd,k,TeX
+f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_3C_cat_git.txt,Pyrolysis/pine_wood_TGA.csv,3,4,Temp,MLR 3,Exp (0 % O_2),gs,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/pine_wood_TGA_exp03_3C_cat_tga.csv,2,3,Temp,Total MLR,FDS (0 % O_2),g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,5 K/min (pine\_wood\_TGA\_3C),Temperature (°C),Normalized Mass Loss Rate (1/s),200,550,1,0,2.00E-03,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/pine_wood_TGA_3C_oxygen,N/A,end,0,pine wood TGA,kd,k,TeX
 s,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
 g,sphere_leak,Complex_Geometry/sphere_leak_git.txt,Complex_Geometry/sphere_leak.csv,1,2,Time,Pressure,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Complex_Geometry/sphere_leak_devc.csv,2,3,Time,Pressure,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Rise (sphere\_leak),Time (s),Pressure (Pa),0,100,1,0,5000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/sphere_leak,Relative Error,max,0.05,Pressure Effects,k+,k,TeX
 d,sphere_radiate,Complex_Geometry/sphere_radiate_git.txt,Complex_Geometry/sphere_radiate.csv,1,2,Time,HF,Exact,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Complex_Geometry/sphere_radiate_devc.csv,2,3,Time,HF1,FDS,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat Flux (sphere\_radiate),Time (s),Heat Flux (kW/m²),0,0.01,1,0,8,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/sphere_radiate,Relative Error,max,0.07,Radiation,bs,b,TeX
diff --git a/Verification/FDS_Cases.sh b/Verification/FDS_Cases.sh
@@ -434,14 +434,17 @@ $QFDS -d Pyrolysis part_baking_soda_420K.fds
 $QFDS -d Pyrolysis part_baking_soda_450K.fds
 $QFDS -d Pyrolysis part_baking_soda_500K.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp03_1C.fds
+$QFDS -d Pyrolysis pine_wood_TGA_exp03_3C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp13_1C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp13_3C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp14_1C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp14_3C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp15_1C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp15_3C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp16_1C.fds
+$QFDS -d Pyrolysis pine_wood_TGA_exp16_3C.fds
 $QFDS -d Pyrolysis pine_wood_TGA_exp17_1C.fds
+$QFDS -d Pyrolysis pine_wood_TGA_exp17_3C.fds
 $QFDS -d Pyrolysis pyrolysis_1.fds
 $QFDS -d Pyrolysis pyrolysis_2.fds
 $QFDS -d Pyrolysis specified_hrr.fds

diff --git a/Verification/Pyrolysis/pine_wood_TGA_exp03_3C.fds b/Verification/Pyrolysis/pine_wood_TGA_exp03_3C.fds
@@ -0,0 +1,27 @@
+&HEAD CHID='pine_wood_TGA_exp03_3C'/
+
+&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.0/
+&SPEC ID='CARBON DIOXIDE'/
+&SPEC ID='FUEL VAPOR'/
+&SPEC ID='WATER VAPOR'/
+
+&SURF ID = 'FUEL'
+      MATL_ID(1,1:3) = 'PINE 1','PINE 2','PINE 3'
+      MATL_MASS_FRACTION(1,1:3) = 0.55,0.1,0.35
+      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/pine_wood_3C_MATL.fds'/
+
+&TAIL /
diff --git a/Verification/Pyrolysis/pine_wood_TGA_exp16_3C.fds b/Verification/Pyrolysis/pine_wood_TGA_exp16_3C.fds
@@ -0,0 +1,27 @@
+&HEAD CHID='pine_wood_TGA_exp16_3C'/
+
+&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.093/
+&SPEC ID='CARBON DIOXIDE'/
+&SPEC ID='FUEL VAPOR'/
+&SPEC ID='WATER VAPOR'/
+
+&SURF ID = 'FUEL'
+      MATL_ID(1,1:3) = 'PINE 1','PINE 2','PINE 3'
+      MATL_MASS_FRACTION(1,1:3) = 0.55,0.1,0.35
+      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/pine_wood_3C_MATL.fds'/
+
+&TAIL /
diff --git a/Verification/Pyrolysis/pine_wood_TGA_exp17_3C.fds b/Verification/Pyrolysis/pine_wood_TGA_exp17_3C.fds
@@ -0,0 +1,27 @@
+&HEAD CHID='pine_wood_TGA_exp17_3C'/
+
+&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.049/
+&SPEC ID='CARBON DIOXIDE'/
+&SPEC ID='FUEL VAPOR'/
+&SPEC ID='WATER VAPOR'/
+
+&SURF ID = 'FUEL'
+      MATL_ID(1,1:3) = 'PINE 1','PINE 2','PINE 3'
+      MATL_MASS_FRACTION(1,1:3) = 0.55,0.1,0.35
+      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/pine_wood_3C_MATL.fds'/
+
+&TAIL /