Skip to content

Commit 4c16b55

Browse files
rmcdermormcdermo
committed
FDS Verification: add 3 component pine_wood_TGA cases to dataplot and tune parameters to better match Anca-Couce data
1 parent 5c8ba22 commit 4c16b55

File tree

7 files changed

+117
-29
lines changed

7 files changed

+117
-29
lines changed

‎Utilities/Input_Libraries/MATL/pine_wood_1C_MATL.fds‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -5,10 +5,10 @@
55
CONDUCTIVITY_RAMP = 'k dry wood'
66
SPECIFIC_HEAT_RAMP = 'c_v dry wood'
77
N_REACTIONS = 2
8-
N_S = 0.87,0.63
9-
N_O2 = 0.0,0.72
10-
A = 4.70E+6,1.45E+10
11-
E = 1.05E+05,1.27E+05
8+
N_S = 0.50,0.63
9+
N_O2 = 0.00,0.72
10+
A = 1.00E+07,1.00E+13
11+
E = 1.00E+05,1.58E+05
1212
MATL_ID(1,1) = 'CHAR'
1313
NU_MATL(1,1) = 0.31
1414
MATL_ID(1,2) = 'CHAR'
@@ -17,25 +17,25 @@
1717
NU_SPEC(1,1) = 0.69
1818
SPEC_ID(1:2,2) = 'OXYGEN','FUEL VAPOR'
1919
NU_SPEC(1:2,2) = -0.1, 0.79
20-
HEAT_OF_REACTION = 416. /
20+
HEAT_OF_REACTION(1:2) = 416./
2121

2222
&MATL ID = 'CHAR'
23-
DENSITY = 73.
23+
DENSITY = 112. ! 112/360 = 0.31 to match yield from AC Tab 7
2424
CONDUCTIVITY_RAMP = 'k + k_r char'
2525
SPECIFIC_HEAT_RAMP = 'c_v char'
2626
N_REACTIONS = 1
2727
N_S = 0.56
2828
N_O2 = 0.68
29-
A = 8.31E+07
30-
E = 1.24E+05
29+
A = 7.00E+06
30+
E = 1.09E+05
3131
MATL_ID = 'ASH'
3232
NU_MATL = 0.02
3333
SPEC_ID = 'OXYGEN','CARBON DIOXIDE'
3434
NU_SPEC = -2.67,3.65
35-
HEAT_OF_REACTION = -32000. /
35+
HEAT_OF_REACTION = -32000./
3636

3737
&MATL ID = 'ASH'
38-
DENSITY = 5.7
38+
DENSITY = 2.0
3939
CONDUCTIVITY_RAMP = 'k + k_r ash'
4040
SPECIFIC_HEAT_RAMP = 'c_v ash'/
4141

‎Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds‎

Lines changed: 19 additions & 19 deletions
Original file line numberDiff line numberDiff line change
@@ -1,13 +1,13 @@
11

2-
&MATL ID = 'PINE 1'
2+
&MATL ID = 'PINE 1' ! peak of anaerobic curve (T_ref approx 360 C)
33
DENSITY = 360.
44
CONDUCTIVITY_RAMP = 'k dry wood'
55
SPECIFIC_HEAT_RAMP = 'c_v dry wood'
66
N_REACTIONS = 2
7-
N_S = 0.56,0.3
8-
N_O2 = 0.0,0.61
9-
A = 4.9E+09,8.9E+09
10-
E = 1.46E+05,1.16E+05
7+
N_S = 0.56,0.50
8+
N_O2 = 0.00,0.61
9+
A = 1.00E+11,1.00E+16
10+
E = 1.50E+05,1.89E+05
1111
MATL_ID(1,1) = 'CHAR'
1212
NU_MATL(1,1) = 0.25
1313
MATL_ID(1,2) = 'CHAR'
@@ -16,17 +16,17 @@
1616
NU_SPEC(1,1) = 0.75
1717
SPEC_ID(1:2,2) = 'OXYGEN','FUEL VAPOR'
1818
NU_SPEC(1:2,2) = -0.1, 0.85
19-
HEAT_OF_REACTION = 416. /
19+
HEAT_OF_REACTION(1:2) = 416./
2020

21-
&MATL ID = 'PINE 2'
21+
&MATL ID = 'PINE 2' ! first hump of anaerobic curve
2222
DENSITY = 360.
2323
CONDUCTIVITY_RAMP = 'k dry wood'
2424
SPECIFIC_HEAT_RAMP = 'c_v dry wood'
2525
N_REACTIONS = 2
2626
N_S = 1,1
27-
N_O2 = 0.0,0.49
28-
A = 5.0E+10,2.0E+04
29-
E = 1.44E+05,7.5E+04
27+
N_O2 = 0.00,0.49
28+
A = 1.00E+15,1.00E+05
29+
E = 1.93E+05,8.50E+04
3030
MATL_ID(1,1) = 'CHAR'
3131
NU_MATL(1,1) = 0.25
3232
MATL_ID(1,2) = 'CHAR'
@@ -35,17 +35,17 @@
3535
NU_SPEC(1,1) = 0.75
3636
SPEC_ID(1:2,2) = 'OXYGEN','FUEL VAPOR'
3737
NU_SPEC(1:2,2) = -0.1, 0.85
38-
HEAT_OF_REACTION = 416. /
38+
HEAT_OF_REACTION(1:2) = 416./ AC Tab 8
3939

40-
&MATL ID = 'PINE 3'
40+
&MATL ID = 'PINE 3' ! broad third peak (tail of anaerobic curve)
4141
DENSITY = 360.
4242
CONDUCTIVITY_RAMP = 'k dry wood'
4343
SPECIFIC_HEAT_RAMP = 'c_v dry wood'
4444
N_REACTIONS = 2
4545
N_S = 1.25,5.67
46-
N_O2 = 0.0,0.66
47-
A = 2.9E+11,2.6
48-
E = 1.64E+05,1.64E+05
46+
N_O2 = 0.00,0.66
47+
A = 1.00,0.00
48+
E = 4.00E+04,1.67E+05
4949
MATL_ID(1,1) = 'CHAR'
5050
NU_MATL(1,1) = 0.25
5151
MATL_ID(1,2) = 'CHAR'
@@ -54,16 +54,16 @@
5454
NU_SPEC(1,1) = 0.75
5555
SPEC_ID(1:2,2) = 'OXYGEN','FUEL VAPOR'
5656
NU_SPEC(1:2,2) = -0.1, 0.85
57-
HEAT_OF_REACTION = 416. /
57+
HEAT_OF_REACTION(1:2) = 416./
5858

5959
&MATL ID = 'CHAR'
60-
DENSITY = 73.
60+
DENSITY = 90.
6161
CONDUCTIVITY_RAMP = 'k + k_r char'
6262
SPECIFIC_HEAT_RAMP = 'c_v char'
6363
N_REACTIONS = 1
6464
N_S = 0.56
6565
N_O2 = 0.68
66-
A = 7.55E+07
66+
A = 8.31E+07
6767
E = 1.24E+05
6868
MATL_ID = 'ASH'
6969
NU_MATL = 0.02
@@ -72,7 +72,7 @@
7272
HEAT_OF_REACTION = -32000. /
7373

7474
&MATL ID = 'ASH'
75-
DENSITY = 5.7
75+
DENSITY = 2.0
7676
CONDUCTIVITY_RAMP = 'k + k_r ash'
7777
SPECIFIC_HEAT_RAMP = 'c_v ash'/
7878

‎Utilities/Matlab/FDS_verification_dataplot_inputs.csv‎

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -728,6 +728,10 @@ f,pine_wood_TGA,Pyrolysis/pine_wood_TGA_exp13_1C_cat_git.txt,Pyrolysis/pine_wood
728728
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
729729
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
730730
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
731+
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
732+
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
733+
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
734+
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
731735
s,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
732736
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
733737
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

‎Verification/FDS_Cases.sh‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -434,14 +434,17 @@ $QFDS -d Pyrolysis part_baking_soda_420K.fds
434434
$QFDS -d Pyrolysis part_baking_soda_450K.fds
435435
$QFDS -d Pyrolysis part_baking_soda_500K.fds
436436
$QFDS -d Pyrolysis pine_wood_TGA_exp03_1C.fds
437+
$QFDS -d Pyrolysis pine_wood_TGA_exp03_3C.fds
437438
$QFDS -d Pyrolysis pine_wood_TGA_exp13_1C.fds
438439
$QFDS -d Pyrolysis pine_wood_TGA_exp13_3C.fds
439440
$QFDS -d Pyrolysis pine_wood_TGA_exp14_1C.fds
440441
$QFDS -d Pyrolysis pine_wood_TGA_exp14_3C.fds
441442
$QFDS -d Pyrolysis pine_wood_TGA_exp15_1C.fds
442443
$QFDS -d Pyrolysis pine_wood_TGA_exp15_3C.fds
443444
$QFDS -d Pyrolysis pine_wood_TGA_exp16_1C.fds
445+
$QFDS -d Pyrolysis pine_wood_TGA_exp16_3C.fds
444446
$QFDS -d Pyrolysis pine_wood_TGA_exp17_1C.fds
447+
$QFDS -d Pyrolysis pine_wood_TGA_exp17_3C.fds
445448
$QFDS -d Pyrolysis pyrolysis_1.fds
446449
$QFDS -d Pyrolysis pyrolysis_2.fds
447450
$QFDS -d Pyrolysis specified_hrr.fds

‎Verification/Pyrolysis/pine_wood_TGA_exp03_3C.fds‎

Lines changed: 27 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,27 @@
1+
&HEAD CHID='pine_wood_TGA_exp03_3C'/
2+
3+
&MESH IJK=10,10,10,XB=-.5,.5,-.5,.5,0,1/
4+
5+
&TIME T_END=1. /
6+
7+
&SPEC ID='NITROGEN', BACKGROUND=T/
8+
&SPEC ID='OXYGEN', MASS_FRACTION_0=0.0/
9+
&SPEC ID='CARBON DIOXIDE'/
10+
&SPEC ID='FUEL VAPOR'/
11+
&SPEC ID='WATER VAPOR'/
12+
13+
&SURF ID = 'FUEL'
14+
MATL_ID(1,1:3) = 'PINE 1','PINE 2','PINE 3'
15+
MATL_MASS_FRACTION(1,1:3) = 0.55,0.1,0.35
16+
SURFACE_VOLUME_RATIO = 7590
17+
LENGTH = 0.1
18+
GEOMETRY = 'CYLINDRICAL'
19+
TGA_ANALYSIS = T
20+
TGA_HEATING_RATE = 5./
21+
22+
&PART ID='FUEL', SAMPLING_FACTOR=1, SURF_ID='FUEL' /
23+
&INIT PART_ID='FUEL', XYZ=0.0,0.0,0.3, N_PARTICLES=1/
24+
25+
&CATF OTHER_FILES='../../Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds'/
26+
27+
&TAIL /

‎Verification/Pyrolysis/pine_wood_TGA_exp16_3C.fds‎

Lines changed: 27 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,27 @@
1+
&HEAD CHID='pine_wood_TGA_exp16_3C'/
2+
3+
&MESH IJK=10,10,10,XB=-.5,.5,-.5,.5,0,1/
4+
5+
&TIME T_END=1. /
6+
7+
&SPEC ID='NITROGEN', BACKGROUND=T/
8+
&SPEC ID='OXYGEN', MASS_FRACTION_0=0.093/
9+
&SPEC ID='CARBON DIOXIDE'/
10+
&SPEC ID='FUEL VAPOR'/
11+
&SPEC ID='WATER VAPOR'/
12+
13+
&SURF ID = 'FUEL'
14+
MATL_ID(1,1:3) = 'PINE 1','PINE 2','PINE 3'
15+
MATL_MASS_FRACTION(1,1:3) = 0.55,0.1,0.35
16+
SURFACE_VOLUME_RATIO = 7590
17+
LENGTH = 0.1
18+
GEOMETRY = 'CYLINDRICAL'
19+
TGA_ANALYSIS = T
20+
TGA_HEATING_RATE = 5./
21+
22+
&PART ID='FUEL', SAMPLING_FACTOR=1, SURF_ID='FUEL' /
23+
&INIT PART_ID='FUEL', XYZ=0.0,0.0,0.3, N_PARTICLES=1/
24+
25+
&CATF OTHER_FILES='../../Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds'/
26+
27+
&TAIL /

‎Verification/Pyrolysis/pine_wood_TGA_exp17_3C.fds‎

Lines changed: 27 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,27 @@
1+
&HEAD CHID='pine_wood_TGA_exp17_3C'/
2+
3+
&MESH IJK=10,10,10,XB=-.5,.5,-.5,.5,0,1/
4+
5+
&TIME T_END=1. /
6+
7+
&SPEC ID='NITROGEN', BACKGROUND=T/
8+
&SPEC ID='OXYGEN', MASS_FRACTION_0=0.049/
9+
&SPEC ID='CARBON DIOXIDE'/
10+
&SPEC ID='FUEL VAPOR'/
11+
&SPEC ID='WATER VAPOR'/
12+
13+
&SURF ID = 'FUEL'
14+
MATL_ID(1,1:3) = 'PINE 1','PINE 2','PINE 3'
15+
MATL_MASS_FRACTION(1,1:3) = 0.55,0.1,0.35
16+
SURFACE_VOLUME_RATIO = 7590
17+
LENGTH = 0.1
18+
GEOMETRY = 'CYLINDRICAL'
19+
TGA_ANALYSIS = T
20+
TGA_HEATING_RATE = 5./
21+
22+
&PART ID='FUEL', SAMPLING_FACTOR=1, SURF_ID='FUEL' /
23+
&INIT PART_ID='FUEL', XYZ=0.0,0.0,0.3, N_PARTICLES=1/
24+
25+
&CATF OTHER_FILES='../../Utilities/Input_Libraries/MATL/pine_wood_3C_MATL.fds'/
26+
27+
&TAIL /

0 commit comments

Comments
 (0)