FDS Verification: Add bi-dir verification case.

Author: drjfloyd

Manuals/Bibliography/FDS_general.bib

@@ -4068,6 +4068,17 @@ @ARTICLE{Mawhinney:FT2012
   year         = {2012}
 }
 
+@ARTICLE{McCaffrey:1976,
+	author      = {B.J. McCaffrey and G. Heskestad},
+	title       = {{A Robust Bidirectional Low-Velocity Probe
+	for Flame and Fire Application}},
+	journal     = {Combustion and Flame},
+	volume      = {26},
+	number      = {},
+	pages       = {125-127},
+	year        = {1976},
+}
+
 @TECHREPORT{McCaffrey:NBSIR_79-1910,
   author       = {McCaffrey, B.J.},
   title        = {{Purely Buoyant Diffusion Flames: Some Experimental Results}},

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -10550,7 +10550,7 @@ \subsection{Thermocouples}
    \tau = \frac{D_{\rm eff} \, \rho_{\rm eff} \, c_{\rm eff} }{6 \, h}  \quad ; \quad  h=\frac{k \, \NU}{D_{\rm eff}} \quad ; \quad \NU = 2 + 0.6 \, \RE^{1/2} \, \PR^{1/3}  \quad ; \quad
    \RE = \frac{\rho \|\bu\| D_{\rm eff}}{\mu}  \label{TC_tau}
 \ee
-For a given value of the time constant, $\tau$, and effective thermal properties, Eq.~(\ref{TC_tau}) can be solved implicitly for the effective diameter. The \ct{TIME_CONSTANT} is specified on a \ct{PROP} line which is identified by the \ct{DEVC} line using a \ct{PROP_ID}. If you specify the \ct{TIME_CONSTANT}, you can still specify the effective \ct{EMISSIVITY}, \ct{DENSITY}, \ct{SPECIFIC_HEAT}, and \ct{HEAT_TRANSFER_COEFFICIENT} as well, but not the \ct{DIAMETER} because this will be calculated automatically. In most cases, it is sufficient to simply specify the \ct{TIME_CONSTANT}.
+For a given value of the time constant, $\tau$, and effective thermal properties, Eq.~(\ref{TC_tau}) can be solved implicitly for the effective diameter. The \ct{TIME_CONSTANT} is specified on a \ct{PROP} line which is identified by the \ct{DEVC} line using a \ct{PROP_ID}. If you specify the \ct{TIME_CONSTANT}, you can still specify the effective \ct{EMISSIVITY}, \ct{DENSITY}, \ct{SPECIFIC_HEAT} (or \ct{SPECIFIC_HEAT_RAMP}), and \ct{HEAT_TRANSFER_COEFFICIENT} as well, but not the \ct{DIAMETER} because this will be calculated automatically. In most cases, it is sufficient to simply specify the \ct{TIME_CONSTANT}.
 
 Figure~\ref{fig:thermocouple_time_constant} shows the results of a simple test case called \ct{thermocouple_time_constant} whose input file is in the \ct{Heat_Transfer} samples folder. Three thermocouples with given time constants of 0.5~s, 3.0~s, and 8.0~s are suddenly subjected to a 20~m/s air stream of 30~$^\circ$C. It is expected that each TC should reach a temperature of 26.32~$^\circ$C at their given time constants.
 \begin{figure}[!ht]
@@ -10562,13 +10562,24 @@ \subsection{Thermocouples}
 
 \subsection{Bi-Directional Probe}
 \label{info:bidir_probe}
+\label{bi_dir}
 
 The output quantity \ct{BI-DIRECTIONAL PROBE} is the velocity of a modeled bi-directional probe. A bi-directional probe uses the following equation:
 \be
 C \sqrt{\frac{2 \Delta P}{\rho}}
 \label{BDP}
 \ee
-where $C$ is a calibration constant (default value is 0.93), $\Delta P$ is the pressure difference across the probe, and $\rho$ is the gas density at the probe. In a typical experiment, the gas density is computed assuming standard pressure (101325 Pa), the molecular weight of air (28.8 g/mol), and the temperature as measured by a thermocouple near the probe. Bi-directional probes have biases due to both the Reynolds number (based on the probe diameter) of the flow and the angle of the flow with respect to the probe axis. This model accounts for those sensitivities and the impact of density differences from varied molecular weight at the probe. The orientation of the probe can be specified with either \ct{IOR} or \ct{ORIENTATION} on \ct{DEVC}. A probe with \ct{IOR}=-1 would have a positive velocity output when the flow is in the negative x direction. Parameters for the probe can be specified with a \ct{PROP_ID} on the \ct {DEVC}. The calibration constant and the probe diameter (default of 0.0254 m) can be set respectively with \ct{CALIBRATION_CONSTANT} and \ct{PROBE_DIAMETER} on \ct{PROP}. If the probe temperature is an aspirated thermocouple or other measurement not sensitive to the radiative environment, then set \ct{TC=F} on \ct{PROP}.
+where $C$ is a calibration constant (default value is 0.93), $\Delta P$ is the pressure difference across the probe, and $\rho$ is the gas density at the probe. In a typical experiment, the gas density is computed assuming standard pressure (101325 Pa), the molecular weight of air (28.8 g/mol), and the temperature as measured by a thermocouple near the probe.
+
+Bi-directional probes have biases due to both the Reynolds number (based on the probe diameter) of the flow and the angle of the flow with respect to the probe axis~\cite{McCaffrey:1976}. At low Reynolds number a probe will measure a higher effective velocity. As the angle of the flow vector with the axis increases, the effective velocity at first increases up to an angle of 30$^\circ$ due to a low pressure region forming downstream of the probe, and then decreases reaching no measured flow at an angle of 90$^\circ$. This model accounts for these sensitivities and the impact of density differences from varied molecular weight at the probe. The orientation of the probe can be specified with either \ct{IOR} or \ct{ORIENTATION} on \ct{DEVC}. A probe with \ct{IOR}=-1 would have a positive velocity output when the flow is in the -x direction. Parameters for the probe can be specified with a \ct{PROP_ID} on the \ct {DEVC}. The calibration constant (default of 0.93) and the probe diameter (default of 0.0254 m) can be set respectively with \ct{CALIBRATION_CONSTANT} and \ct{PROBE_DIAMETER} on \ct{PROP}. If the probe temperature is an aspirated thermocouple or other measurement not sensitive to the radiative environment, then set \ct{TC=F} on \ct{PROP}. Thermocouple specific properties for a bi-directional probe, see Section~\ref{info:THERMOCOUPLE}, should be set with the same \ct{PROP} as for the probe.
+
+Figure~\ref{bi_dir_fig} shows the results of a bi-directional probe with varying angle to a 1~m/s flow and varying flow speed. 
+\begin{figure}[ht]
+	\includegraphics[width=3in]{SCRIPT_FIGURES/bi_dir}
+	\includegraphics[width=3in]{SCRIPT_FIGURES/bi_dir_2}
+	\caption[Results of the \ct{bi_dir} test case]{Measured velocities using a bi-directional probe. (Left) Varying angle for a 1~m/s flow. (Right) Probe axis aligned flow with varied flow speed.}
+	\label{bi_dir_fig}
+\end{figure}
 
 \subsection{Volume Flow}
 \label{info:volume_flow}

Utilities/Matlab/FDS_verification_dataplot_inputs.csv

@@ -36,6 +36,8 @@ d,back_wall_test_2,Heat_Transfer/back_wall_test_2_git.txt,Heat_Transfer/back_wal
 d,back_wall_test_2,Heat_Transfer/back_wall_test_2_git.txt,Heat_Transfer/back_wall_test_2_devc.csv,2,3,Time,T_F_Target|T_B_Target,Exact Front|Exact Back,ko|ro,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Heat_Transfer/back_wall_test_2_devc.csv,2,3,Time,T_F3|T_B3,FDS Front|FDS Back,k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Rotated GEOM Surface Temp. (back\_wall\_test\_2),Time (s),Temperature (°C),0,20,1,0,1000,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/back_wall_test_2_temp_geom_rotated,Relative Error,end,0.02,Heat Transfer,r^,r,TeX
 d,beam_detector,Detectors/beam_detector_git.txt,Detectors/beam_detector.csv,1,2,Time,Obs|Obs|Obs,Exact|Exact|Exact,k-|k-|k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Detectors/beam_detector_devc.csv,2,3,Time,beam_1|beam_2|beam_3,FDS (beam\_1)|FDS (beam\_2)|FDS (beam\_3),ko|ro|go,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Path Obscuration (beam\_detector),Time (s),Path Obscuration (%),0,1,1,0,110,1,no,0.25 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/beam_detector_obs,Relative Error,end,0.01,Aerosols,r^,r,TeX
 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,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

Verification/Controls/bi_dir.csv

@@ -0,0 +1,5 @@
+Time,V1_0,V1_180,V1_30,V1_55,Vp1_0,V10_0
+5,1.03,-1.03,1.1,0.86,12.7,10.2
+10,1.03,-1.03,1.1,0.86,12.7,10.2
+15,1.03,-1.03,1.1,0.86,12.7,10.2
+20,1.03,-1.03,1.1,0.86,12.7,10.2

Verification/Controls/bi_dir.fds

@@ -0,0 +1,31 @@
+&HEAD CHID='bi_dir',TITLE='Bi-directional probe test'/
+
+&MISC STRATIFICATION=F/
+&RADI RADIATION=F/
+
+&TIME T_END=20/
+
+&DUMP NFRAMES=20/
+
+&MESH XB=0.0,1.0,0,1,0,1,IJK=10,10,10/
+&MESH XB=1.1,2.1,0,1,0,1,IJK=10,10,10/
+&MESH XB=2.2,3.2,0,1,0,1,IJK=10,10,10/
+
+&SURF ID='SLIP',NO_SLIP=T,DEFAULT=T/
+&SURF ID='0p1',VEL=-0.1,COLOR='RED'/
+&SURF ID='1p0',VEL=-1,COLOR='GREEN'/
+&SURF ID='10p0',VEL=-10,COLOR='BLUE'/
+
+&VENT PBX=0.0,SURF_ID='0p1'/
+&VENT PBX=1.1,SURF_ID='1p0'/
+&VENT PBX=2.2,SURF_ID='10p0'/
+
+&VENT MB='XMAX',SURF_ID='OPEN'/
+
+&DEVC XYZ=0.5,0.5,0.5,QUANTITY='BI-DIRECTIONAL PROBE',ORIENTATION=1,0,0,ID='Vp1_0',CONVERSION_FACTOR=100/
+&DEVC XYZ=1.6,0.5,0.5,QUANTITY='BI-DIRECTIONAL PROBE',ORIENTATION=1,0,0,ID='V1_0'/
+&DEVC XYZ=1.6,0.5,0.5,QUANTITY='BI-DIRECTIONAL PROBE',ORIENTATION=-1,0,0,ID='V1_180'/
+&DEVC XYZ=1.6,0.5,0.5,QUANTITY='BI-DIRECTIONAL PROBE',ORIENTATION=0.8660,0.5000,0,ID='V1_30'/
+&DEVC XYZ=1.6,0.5,0.5,QUANTITY='BI-DIRECTIONAL PROBE',ORIENTATION=0.5736,0.8192,0,ID='V1_55'/
+&DEVC XYZ=2.7,0.5,0.5,QUANTITY='BI-DIRECTIONAL PROBE',ORIENTATION=1,0,0,ID='V10_0'/
+

Verification/FDS_Cases.sh

@@ -177,6 +177,7 @@ $QFDS -p 5 -d Complex_Geometry geom_stretched_grid.fds
 $QFDS -p 3 -d Complex_Geometry thin_object_mass.fds
 
 $QFDS -d Controls activate_vents.fds
+$QFDS -d Controls bi_dir.fds
 $QFDS -d Controls control_test.fds
 $QFDS -d Controls control_test_2.fds
 $QFDS -d Controls create_remove.fds