Merge pull request #15296 from mcgratta/master

FDS Verification: Remove birch_tga.m

Author: mcgratta

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

@@ -38,6 +38,8 @@ d,beam_detector,Detectors/beam_detector_git.txt,Detectors/beam_detector.csv,1,2,
 d,beam_detector,Detectors/beam_detector_git.txt,Detectors/beam_detector_ss_check.csv,1,2,L,value,Exact,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Detectors/beam_detector_ss_smv.csv,1,2,L,value,Smokeview,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Smoke Obscuration (beam\_detector),Distance (m),Pixel Value,0,10,1,0,255,1,no,0.25 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/beam_detector,N/A,end,0.02,Aerosols,r^,r,TeX
 d,bi_dir,Controls/bi_dir_git.txt,Controls/bi_dir.csv,1,2,Time,V1_0|V1_180|V1_30|V1_55,Expected 0^\circ|Expected 180^\circ|Expected 30^\circ|Expected 55^\circ,ko|ro|bo|go,0,100000,,5,100000,-1.00E+09,1.00E+09,0,Controls/bi_dir_devc.csv,2,3,Time,V1_0|V1_180|V1_30|V1_55,FDS 0^\circ|FDS 180^\circ|FDS 30^\circ|FDS 55^\circ,k-|r-|b-|g-,0,100000,,5,100000,-1.00E+09,1.00E+09,0,Measured Velocity (bi\_dir),Time (s),Velocity (m/s),0,20,1,-1.5,1.5,1,no,0.25 0.92,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/bi_dir,Relative Error,mean,0.02,Controls,kd,k,TeX
 d,bi_dir,Controls/bi_dir_git.txt,Controls/bi_dir.csv,1,2,Time,Vp1_0|V10_0,Expected 0.1 m/s \times 100 |Expected  10 m/s,ko|ro,0,100000,,5,100000,-1.00E+09,1.00E+09,0,Controls/bi_dir_devc.csv,2,3,Time,Vp1_0|V10_0,FDS 0.1 m/s \times 100 |FDS  10 m/s,k-|r-,0,100000,,5,100000,-1.00E+09,1.00E+09,0,Measured Velocity (bi\_dir),Time (s),Velocity (m/s),0,20,1,0,15,1,no,0.25 0.92,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/bi_dir_2,Relative Error,mean,0.02,Controls,kd,k,TeX
+d,birch_tga,Pyrolysis/birch_tga_1step_2_git.txt,Pyrolysis/birch_tga_2_exp.csv,1,2,Temp,Mass,Exp (2 K/min),b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/birch_tga_1step_2_devc.csv,2,3,T_wall,Mass Fraction,FDS (2 K/min),b--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,TGA Test (birch_tga),Temperature ($^\circ$C),Mass Fraction,0,800,1,0,1.1,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/birch_tga,Relative Error,end,0.01,Pyrolysis,kd,k,LaTeX
+f,birch_tga,Pyrolysis/birch_tga_1step_2_git.txt,Pyrolysis/birch_tga_20_exp.csv,1,2,Temp,Mass,Exp (20 K/min),r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolysis/birch_tga_1step_20_devc.csv,2,3,T_wall,Mass Fraction,FDS (20 K/min),r--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,TGA Test (birch_tga),Temperature ($^\circ$C),Mass Fraction,0,800,1,0,1.1,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/birch_tga,Relative Error,end,0.01,Pyrolysis,kd,k,LaTeX
 d,bound_test_1,Species/bound_test_1_git.txt,Species/bound_test_1.csv,1,2,Time,F1|A1|P1|F2|A2|P2|Sum,F1|A1|P1|F2|A2|P2|Sum,ko|ro|bo|go|mo|co|k+,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Species/bound_test_1_devc.csv,2,3,Time,XF1|XA1|XP1|XF2|XA2|XP2|SUM,FDS F1|FDS A1|FDS P1|FDS F2|FDS A2|FDS P2|FDS SUM,k-|r-|b-|g-|m-|c-|k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Boundedness Test (bound\_test\_1),Time (s),Volume Fraction,0,1,1,0,1.05,1,no,0.05 0.90,EastOutside,,1.3,linear,FDS_Verification_Guide/SCRIPT_FIGURES/bound_test_1,Absolute Error,end,0.001,Species,kd,k,TeX
 d,bound_test_2,Species/bound_test_2_git.txt,Species/bound_test_2.csv,1,2,Time,F1|F2|A1|A2|P1|P2|P3|Sum,F1|F2|A1|A2|P1|P2|P3|Sum,ko|ro|bo|go|mo|co|yo|k+,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Species/bound_test_2_devc.csv,2,3,Time,XF1|XF2|XA1|XA2|XP1|XP2|XP3|SUM,FDS F1|FDS F2|FDS A1|FDS A2|FDS P1|FDS P2|FDS P3|FDS SUM,k-|r-|b-|g-|m-|c-|y-|k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Boundedness Test (bound\_test\_2),Time (s),Volume Fraction,0,1,1,0,1.05,1,no,0.05 0.90,EastOutside,,1.3,linear,FDS_Verification_Guide/SCRIPT_FIGURES/bound_test_2,Absolute Error,end,0.001,Species,kd,k,TeX
 d,box_burn_away1,Fires/box_burn_away1_git.txt,Fires/box_burn_away.csv,1,2,Time,Mass (kg),Ideal,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Fires/box_burn_away1_devc.csv,2,3,Time,Mass fuel,FDS (fuel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pyrolyzed Mass (box\_burn\_away1),Time (s),Mass (kg),0,30,1,0,1.5,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/box_burn_away1,Relative Error,end,0.02,Fires,r*,r,TeX

Utilities/Matlab/FDS_verification_script.m

@@ -65,8 +65,6 @@
 disp('turb_model...');                    turb_model
 disp('jet_decay...');                     jet_decay
 disp('wall_model...');                    wall_model
-disp('pyrolysis...');                     pyrolysis
-disp('birch_tga...');                     birch_tga
 disp('pulsating...');                     pulsating
 disp('compression_wave...');              compression_wave
 disp('soborot_mass_transport...');        soborot_mass_transport

Utilities/Matlab/scripts/birch_tga.m

@@ -1,6 +1,6 @@
 &HEAD CHID='birch_tga_1step_2', TITLE='TGA TEST of evaporation' / 
 
-&MESH IJK=3,1,4, XB=-2,2,-0.5,0.5,0,1 / 
+&MESH IJK=3,1,4, XB=-3,3,-0.5,0.5,0,1 / 
 
 &TIME T_END=24000., WALL_INCREMENT=1, DT=1.0 /
 
@@ -74,8 +74,8 @@
 4800	= 10 K/min
 2400	= 20 K/min
 
-&DEVC XYZ = 0.0,0.0,0.0,IOR = 3, QUANTITY = 'SURFACE DENSITY' /
-&DEVC XYZ = 0.0,0.0,1.0,IOR = -3, QUANTITY = 'WALL TEMPERATURE' /
-&DEVC XYZ = 0.0,0.0,0.0,IOR = 3, QUANTITY = 'WALL TEMPERATURE' /
+&DEVC ID='Mass Fraction', XYZ=0.0,0.0,0.0, IOR= 3, QUANTITY='NORMALIZED MASS' /
+&DEVC ID='T_wall',        XYZ=0.0,0.0,1.0, IOR=-3, QUANTITY='WALL TEMPERATURE' /
+&DEVC ID='T_sample',      XYZ=0.0,0.0,0.0, IOR= 3, QUANTITY='WALL TEMPERATURE' /
 
 &TAIL /

Utilities/Matlab/scripts/pyrolysis.m

@@ -1,6 +1,6 @@
 &HEAD CHID='birch_tga_1step_20', TITLE='TGA TEST of evaporation' / 
 
-&MESH IJK=3,1,4, XB=-2,2,-0.5,0.5,0,1 /  
+&MESH IJK=3,1,4, XB=-3,3,-0.5,0.5,0,1 /  
 
 &TIME T_END=2400., WALL_INCREMENT=1, DT=1.0 /  
 
@@ -74,8 +74,8 @@
 4800	= 10 K/min
 2400	= 20 K/min
 
-&DEVC XYZ = 0.0,0.0,0.0,IOR = 3, QUANTITY = 'SURFACE DENSITY' /
-&DEVC XYZ = 0.0,0.0,1.0,IOR = -3, QUANTITY = 'WALL TEMPERATURE' /
-&DEVC XYZ = 0.0,0.0,0.0,IOR = 3, QUANTITY = 'WALL TEMPERATURE' /
+&DEVC ID='Mass Fraction', XYZ=0.0,0.0,0.0, IOR= 3, QUANTITY='NORMALIZED MASS' /
+&DEVC ID='T_wall',        XYZ=0.0,0.0,1.0, IOR=-3, QUANTITY='WALL TEMPERATURE' /
+&DEVC ID='T_sample',      XYZ=0.0,0.0,0.0, IOR= 3, QUANTITY='WALL TEMPERATURE' /
 
 &TAIL /