Skip to content

Commit

Permalink
FDS Verification: add 3 component pine_wood_TGA cases to dataplot and…
Browse files Browse the repository at this point in the history 
… tune parameters to better match Anca-Couce data
  • Loading branch information
@rmcdermo
rmcdermo committed Dec 23, 2024
1 parent 5c8ba22 commit 4c16b55
Show file tree
Hide file tree
Showing 7 changed files with 117 additions and 29 deletions.
20 changes: 10 additions & 10 deletions Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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'
Expand All @@ -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'/

Expand Down
38 changes: 19 additions & 19 deletions Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds
View file
Original file line number Diff line number Diff line change
@@ -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'
Expand All @@ -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'
Expand All @@ -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'
Expand All @@ -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
Expand All @@ -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'/

Expand Down
4 changes: 4 additions & 0 deletions Utilities/Matlab/FDS_verification_dataplot_inputs.csv
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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
3 changes: 3 additions & 0 deletions Verification/FDS_Cases.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
27 changes: 27 additions & 0 deletions Verification/Pyrolysis/pine_wood_TGA_exp03_3C.fds
View file
Original file line number Diff line number Diff line change
@@ -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 /
27 changes: 27 additions & 0 deletions Verification/Pyrolysis/pine_wood_TGA_exp16_3C.fds
View file
Original file line number Diff line number Diff line change
@@ -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 /
27 changes: 27 additions & 0 deletions Verification/Pyrolysis/pine_wood_TGA_exp17_3C.fds
View file
Original file line number Diff line number Diff line change
@@ -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 /

0 comments on commit 4c16b55

Please sign in to comment.