From 4c16b551c1e05f6ce180735e550029a3b25ba438 Mon Sep 17 00:00:00 2001 From: Randy McDermott Date: Mon, 23 Dec 2024 17:22:30 -0500 Subject: [PATCH] FDS Verification: add 3 component pine_wood_TGA cases to dataplot and tune parameters to better match Anca-Couce data --- .../MATL/pine_wood_1C_MATL.fds | 20 +++++----- .../MATL/pine_wood_3C_MATL.fds | 38 +++++++++---------- .../FDS_verification_dataplot_inputs.csv | 4 ++ Verification/FDS_Cases.sh | 3 ++ .../Pyrolysis/pine_wood_TGA_exp03_3C.fds | 27 +++++++++++++ .../Pyrolysis/pine_wood_TGA_exp16_3C.fds | 27 +++++++++++++ .../Pyrolysis/pine_wood_TGA_exp17_3C.fds | 27 +++++++++++++ 7 files changed, 117 insertions(+), 29 deletions(-) create mode 100644 Verification/Pyrolysis/pine_wood_TGA_exp03_3C.fds create mode 100644 Verification/Pyrolysis/pine_wood_TGA_exp16_3C.fds create mode 100644 Verification/Pyrolysis/pine_wood_TGA_exp17_3C.fds diff --git a/Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds b/Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds index 6e103ea58d3..e4788436a6f 100644 --- 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 index 9972daea7a6..ab6e74ce714 100644 --- 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 index 77a8604deb1..b0636a029be 100644 --- 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 index d9c20a06828..89e14ca63db 100755 --- a/Verification/FDS_Cases.sh +++ b/Verification/FDS_Cases.sh @@ -434,6 +434,7 @@ $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 @@ -441,7 +442,9 @@ $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 new file mode 100644 index 00000000000..a7cf376aaa6 --- /dev/null +++ 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 new file mode 100644 index 00000000000..ae4485156cf --- /dev/null +++ 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 new file mode 100644 index 00000000000..e3c125c82c8 --- /dev/null +++ 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 /