Merge pull request #13716 from marcosvanella/master

FDS Verification : Add obst_activate cases with error tests.

Author: marcosvanella

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

@@ -465,6 +465,8 @@ d,ns2d_8_nupt1,NS_Analytical_Solution/ns2d_8_nupt1_git.txt,NS_Analytical_Solutio
 d,ns2d_16_nupt1,NS_Analytical_Solution/ns2d_16_nupt1_git.txt,NS_Analytical_Solution/ns2d_16_nupt1_exact.csv,1,2,Time,u-vel,Analytical (u-vel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_16_nupt1_devc.csv,2,3,Time,UVEL,FDS (UVEL),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (ns2d\_16\_nupt1),Time (s),Velocity (m/s),0,7,1,0.3,2.4,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_16_nupt1,Convergent Series,end,0,NS Analytical Solution,kd,k,TeX
 d,ns2d_32_nupt1,NS_Analytical_Solution/ns2d_32_nupt1_git.txt,NS_Analytical_Solution/ns2d_32_nupt1_exact.csv,1,2,Time,u-vel,Analytical (u-vel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_32_nupt1_devc.csv,2,3,Time,UVEL,FDS (UVEL),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (ns2d\_32\_nupt1),Time (s),Velocity (m/s),0,7,1,0.3,2.4,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_32_nupt1,Convergent Series,end,0,NS Analytical Solution,kd,k,TeX
 d,ns2d_64_nupt1,NS_Analytical_Solution/ns2d_64_nupt1_git.txt,NS_Analytical_Solution/ns2d_64_nupt1_exact.csv,1,2,Time,u-vel,Analytical (u-vel),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,NS_Analytical_Solution/ns2d_64_nupt1_devc.csv,2,3,Time,UVEL,FDS (UVEL),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (ns2d\_64\_nupt1),Time (s),Velocity (m/s),0,7,1,0.3,2.4,1,no,0.05 0.90,NorthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/ns2d_64_nupt1,Relative Error,mean,0.01,NS Analytical Solution,kd,k,TeX
+d,obst_activation_ulmat,Pressure_Solver/obst_activation_ulmat_git.txt,Pressure_Solver/obst_activation_exact.csv,1,2,Time,TEMP,Analytical (TEMP),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/obst_activation_ulmat_devc.csv,2,3,Time,TEMP,FDS (ULMAT),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (obst\_activation),Time (s),Max Temperature (°C),0,2,1,19.5,20.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/obst_activation,Absolute Error,max,0.5,Pressure Solver,kd,k,TeX
+f,obst_activation_default,Pressure_Solver/obst_activation_default_git.txt,Pressure_Solver/obst_activation_exact.csv,1,2,Time,TEMP,Analytical (TEMP),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/obst_activation_default_devc.csv,2,3,Time,TEMP,FDS (FFT),ko--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (obst\_activation),Time (s),Max Temperature (°C),0,2,1,19.5,20.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/obst_activation,Absolute Error,max,0.5,Pressure Solver,kd,k,TeX
 d,obst_coarse_fine_interface,Pressure_Effects/obst_coarse_fine_interface_git.txt,Pressure_Effects/obst_coarse_fine_interface_exact.csv,1,2,Time,DP,Analytical (DP),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Effects/obst_coarse_fine_interface_devc.csv,2,3,Time,DP,FDS (DP),k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Background Pressure (obst\_coarse\_fine\_interface),Time (s),Pressure (Pa),0,10,1,-10,10,1,no,0.05 0.90,NorthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/obst_coarse_fine_interface,Absolute Error,end,2,Pressure Effects,kd,k,TeX
 d,opening_ulmat,Pressure_Solver/opening_ulmat_git.txt,Pressure_Solver/opening_pressure_error.csv,1,2,Time,Pressure Tolerance,Ideal (Pressure Tolerance),ko--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/opening_ulmat_devc.csv,2,3,Time,perr-max,FDS (p err max),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Error (opening\_ulmat),Time (s),Pressure Error (Pa),0,10,1,0,1.00E-06,1,no,0.05 0.90,SouthEast,,1,semilogy,FDS_User_Guide/SCRIPT_FIGURES/opening_ulmat,Absolute Error,tolerance,1.00E-10,Pressure Solver,k+,k,TeX
 d,parabolic_profile,Flowfields/parabolic_profile_git.txt,Flowfields/parabolic_profile.csv,1,2,Time,Pressure,Exact (Pressure),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Flowfields/parabolic_profile_devc.csv,2,3,Time,pres,FDS (pres),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (parabolic\_profile),Time (s),Pressure (Pa),0,60,1,0,2500,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/parabolic_profile,Relative Error,end,0.01,Pressure Effects,k+,k,TeX

Verification/FDS_Cases.sh

@@ -405,6 +405,7 @@ $QFDS -p 8 -d Pressure_Solver pressure_iteration3d_uglmat.fds
 $QFDS -p 16 -d Pressure_Solver random_obstructions_fft.fds
 $QFDS -p 4 -d Pressure_Solver stairwell.fds
 $QFDS -p 2 -d Pressure_Solver ulmat_2zone.fds
+$QFDS -p 4 -d Pressure_Solver obst_activation_default.fds
 $QFDS -p 4 -d Pressure_Solver obst_activation_ulmat.fds
 
 $QFDS -d Pyrolysis cable_11_insulation_mcc.fds

Verification/Pressure_Solver/obst_activation_default.fds

@@ -0,0 +1,49 @@
+&HEAD CHID='obst_activation_default', TITLE='Test of various obst activation situations for default FFT Poisson solver.' /
+
+&MESH IJK=16,1,16, XB=0.0,0.8,-0.001,0.001,0.0,0.8 /
+&MESH IJK=16,1,16, XB=0.8,1.6,-0.001,0.001,0.0,0.8 /
+&MESH IJK=16,1,16, XB=1.6,2.4,-0.001,0.001,0.0,0.8 /
+&MESH IJK=16,1,16, XB=2.4,3.2,-0.001,0.001,0.0,0.8 /
+
+&TIME T_END=2.0 /
+
+&PRES CHECK_POISSON=T /
+
+&MISC SIMULATION_MODE='DNS', STRATIFICATION=.FALSE. /
+&RADI RADIATION=F /
+
+&SPEC ID='LJ AIR', VISCOSITY=1.2E-5, BACKGROUND=.TRUE. /
+
+&SURF ID='WALL', FREE_SLIP=.TRUE., DEFAULT=.TRUE. /
+&SURF ID='PLATE', NO_SLIP=.TRUE. /
+
+&DEVC ID='clock_1', XYZ=0,0,0, QUANTITY='TIME', SETPOINT=0.3, INITIAL_STATE=F /
+&DEVC ID='clock_2', XYZ=0,0,0, QUANTITY='TIME', SETPOINT=0.5, INITIAL_STATE=T /
+&DEVC ID='clock_3', XYZ=0,0,0, QUANTITY='TIME', SETPOINT=0.8, INITIAL_STATE=F /
+&DEVC ID='clock_4', XYZ=0,0,0, QUANTITY='TIME', SETPOINT=1.0, INITIAL_STATE=T /
+&DEVC ID='clock_5', XYZ=0,0,0, QUANTITY='TIME', SETPOINT=1.2, INITIAL_STATE=F /
+&DEVC ID='clock_6', XYZ=0,0,0, QUANTITY='TIME', SETPOINT=1.5, INITIAL_STATE=T /
+
+&OBST XB=0.1,0.2,-0.001,0.001,0.3,0.5, SURF_ID='PLATE', DEVC_ID='clock_1' /
+&OBST XB=1.4,1.6,-0.001,0.001,0.2,0.4, SURF_ID='PLATE', DEVC_ID='clock_2' /
+&OBST XB=2.4,2.5,-0.001,0.001,0.4,0.6, SURF_ID='PLATE', DEVC_ID='clock_3' /
+&OBST XB=3.0,3.2,-0.001,0.001,0.3,0.5, SURF_ID='PLATE', DEVC_ID='clock_4' /
+&OBST XB=0.6,1.0,-0.001,0.001,0.5,0.7, SURF_ID='PLATE', DEVC_ID='clock_5' /
+&OBST XB=1.9,2.1,-0.001,0.001,0.0,0.1, SURF_ID='PLATE', DEVC_ID='clock_6' /
+
+&SURF ID='BLOW', VEL=-0.5 /
+
+&VENT PBX=0.0, SURF_ID='BLOW' /
+&VENT PBX=3.2, SURF_ID='OPEN' /
+
+&SLCF PBY=0.,QUANTITY='PRESSURE', VECTOR=T /
+&SLCF PBY=0.,QUANTITY='DIVERGENCE', CELL_CENTERED=T /
+&SLCF PBY=0.,QUANTITY='ZONE PRESSURE SOLVER TYPE', CELL_CENTERED=T /
+
+&DEVC ID='TEMP', XB=0.0,3.2,-0.001,0.001,0.0,0.8, QUANTITY='TEMPERATURE', SPATIAL_STATISTIC='MAX', TEMPORAL_STATISTIC='INSTANT VALUE'/
+
+&DUMP DT_DEVC=0.5, SIG_FIGS=6 /
+  
+&TAIL /
+
+

Verification/Pressure_Solver/obst_activation_exact.csv

@@ -0,0 +1,6 @@
+Time,TEMP
+0.0,20.0
+0.5,20.0
+1.0,20.0
+1.5,20.0
+2.0,20.0

Verification/Pressure_Solver/obst_activation_ulmat.fds

@@ -1,4 +1,4 @@
-&HEAD CHID='obst_activation_ulmat', TITLE='Test of various obst activation situations for inseparable ULMAT' /
+&HEAD CHID='obst_activation_ulmat', TITLE='Test of various obst activation situations for inseparable ULMAT.' /
 
 &MESH IJK=16,1,16, XB=0.0,0.8,-0.001,0.001,0.0,0.8 /
 &MESH IJK=16,1,16, XB=0.8,1.6,-0.001,0.001,0.0,0.8 /
@@ -10,6 +10,7 @@
 &PRES CHECK_POISSON=T, SOLVER = 'ULMAT' /
 
 &MISC SIMULATION_MODE='DNS', STRATIFICATION=.FALSE. /
+&RADI RADIATION=F /
 
 &SPEC ID='LJ AIR', VISCOSITY=1.2E-5, BACKGROUND=.TRUE. /
 
@@ -39,6 +40,10 @@
 &SLCF PBY=0.,QUANTITY='DIVERGENCE', CELL_CENTERED=T /
 &SLCF PBY=0.,QUANTITY='ZONE PRESSURE SOLVER TYPE', CELL_CENTERED=T /
 
+&DEVC ID='TEMP', XB=0.0,3.2,-0.001,0.001,0.0,0.8, QUANTITY='TEMPERATURE', SPATIAL_STATISTIC='MAX', TEMPORAL_STATISTIC='INSTANT VALUE'/
+
+&DUMP DT_DEVC=0.5, SIG_FIGS=6 /
+  
 &TAIL /