Merge pull request #15030 from cxp484/FireX

FireX: Merge with firemodels/master

Author: cxp484

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -10562,7 +10562,15 @@ \subsection{Thermocouples}
 \label{fig:thermocouple_time_constant}
 \end{figure}
 
+\subsection{Bi-Directional Probe}
+\label{info:bidir_probe}
 
+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}.
 
 \subsection{Volume Flow}
 \label{info:volume_flow}
@@ -11730,6 +11738,7 @@ \section{Device, Control, and Other Miscellaneous Output Quantities}
 \endhead
 \ct{ACTUATED SPRINKLERS}                       & Section~\ref{info:TIMING}                     &                & D            \\ \hline
 \ct{ASPIRATION}                                & Section~\ref{info:aspiration_detector}        & \%/m           & D            \\ \hline
+\ct{BI-DIRECTIONAL PROBE}                      & Section~\ref{info:bidir_probe}                & m/s            & D            \\ \hline
 \ct{CHAMBER OBSCURATION}                       & Section~\ref{info:smoke_detector}             & \%/m           & D            \\ \hline
 \ct{CELL INDEX I}                              & Mesh cell index in x                          &                & D,S          \\ \hline
 \ct{CELL INDEX J}                              & Mesh cell index in y                          &                & D,S          \\ \hline
@@ -13176,6 +13185,7 @@ \section{\texorpdfstring{{\tt PROP}}{PROP} (Device Properties)}
 \ct{ALPHA_E}                          & Real          & Section~\ref{info:smoke_detector}         &                       & 0.                   \\ \hline
 \ct{BETA_C}                           & Real          & Section~\ref{info:smoke_detector}         &                       & 1.                   \\ \hline
 \ct{BETA_E}                           & Real          & Section~\ref{info:smoke_detector}         &                       & 1.                   \\ \hline
+\ct{CALIBRATION_CONSTANT}             & Real          & Section~\ref{info:bidir_probe}            &                       & 0.93                 \\ \hline
 \ct{CHARACTERISTIC_VELOCITY}          & Real          & Section~\ref{info:pressure_coefficient}   & m/s                   & 1.                   \\ \hline
 \ct{C_FACTOR}                         & Real          & Section~\ref{info:sprinklers}             & (m/s)$^{1/2}$         & 0.                   \\ \hline
 \ct{DENSITY}                           & Real          & Section~\ref{info:THERMOCOUPLE}           & kg/m$^3$              & 8908.                \\ \hline
@@ -13212,6 +13222,7 @@ \section{\texorpdfstring{{\tt PROP}}{PROP} (Device Properties)}
 \ct{PDPA_RADIUS}                      & Real          & Section~\ref{PDPA}                        & m                     & 0.                   \\ \hline
 \ct{PDPA_START}                       & Real          & Section~\ref{PDPA}                        & s                     & 0.                   \\ \hline
 \ct{PRESSURE_RAMP}                    & Character     & Section~\ref{info:sprinklers}             &                       &                      \\ \hline
+\ct{PROBE_DIAMETER}                    & Real          & Section~\ref{info:bidir_probe}           & m                     &  0.0254              \\ \hline
 \ct{P0}                                & Real          & Section~\ref{info:velocity_patch}         & m/s                   &  0.                  \\ \hline
 \ct{PX(3)}                             & Real          & Section~\ref{info:velocity_patch}         & m/s                   &  0.                  \\ \hline
 \ct{PXX(3,3)}                          & Real          & Section~\ref{info:velocity_patch}         & m/s                   &  0.                  \\ \hline
@@ -13227,6 +13238,7 @@ \section{\texorpdfstring{{\tt PROP}}{PROP} (Device Properties)}
 \ct{SPRAY_PATTERN_MU}                & Real          & Section~\ref{info:sprinklers}             & degrees               & 0.                   \\ \hline
 \ct{SPRAY_PATTERN_SHAPE}             & Character     & Section~\ref{info:sprinklers}             &                       & \ct{'GAUSSIAN'}     \\ \hline
 \ct{SPRAY_PATTERN_TABLE}             & Character     & Section~\ref{info:sprinklers}             &                       &                      \\ \hline
+\ct{TC}                               & Logical       & Section~\ref{info:bidir_probe}            &                       & \ct{T}               \\ \hline
 \ct{TIME_CONSTANT}                    & Real          & Section~\ref{info:THERMOCOUPLE}           & s                     &                      \\ \hline
 \ct{VELOCITY_COMPONENT}               & Integer       & Section~\ref{info:velocity_patch}         &                       &                      \\ \hline
 \ct{VIEW_ANGLE}                       & Real          & Section~\ref{info:heat_flux}              & degrees               & 180.                 \\ \hline
@@ -14449,6 +14461,7 @@ \chapter{Error Codes}
 885 \> \ct{DEVC ... must have coordinates, even if not a point quantity.} \> Section~\ref{info:DEVC} \\
 886 \> \ct{DEVC ... must have either an output QUANTITY or PROP_ID.} \> Section~\ref{info:DEVC} \\
 887 \> \ct{DEVC ... must have an ORIENTATION.} \> Section~\ref{info:heat_flux} \\
+    \>                                         \> or Section~\ref{info:bidir_probe} \\
 888 \> \ct{DEVC ... TEMPORAL_STATISTIC is not recognized.} \> Section~\ref{info:statistics} \\
 889 \> \ct{DEVC ... SPATIAL_STATISTIC is not recognized.} \> Section~\ref{info:statistics} \\
 890 \> \ct{DEVC ... STATISTICS is not recognized.} \> Section~\ref{info:statistics} \\

Source/data.f90

@@ -608,6 +608,12 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(113)%UNITS = 'kW/m2'
 OUTPUT_QUANTITY(113)%SHORT_NAME  = 'hflux'
 
+! Model of a Bi-Directional Probe
+
+OUTPUT_QUANTITY(114)%NAME  = 'BI-DIRECTIONAL PROBE'
+OUTPUT_QUANTITY(114)%UNITS = 'm/s'
+OUTPUT_QUANTITY(114)%SHORT_NAME  = 'bdp'
+
 ! Special Outputs for Partially Stirred Batch Reactor Model
 
 OUTPUT_QUANTITY(130)%NAME  = 'EXTINCTION'

Source/devc.f90

@@ -11,13 +11,13 @@ MODULE DEVICE_VARIABLES
 
 TYPE PROPERTY_TYPE
    REAL(EB) :: DENSITY,DIAMETER,EMISSIVITY,HEAT_TRANSFER_COEFFICIENT,SPECIFIC_HEAT,RTI,TIME_CONSTANT, &
-               ACTIVATION_TEMPERATURE,ACTIVATION_OBSCURATION, &
+               ACTIVATION_TEMPERATURE,ACTIVATION_OBSCURATION, CALIBRATION_CONSTANT,&
                ALPHA_E,ALPHA_C,BETA_E,BETA_C,CHARACTERISTIC_VELOCITY,PARTICLE_VELOCITY,MASS_FLOW_RATE,FLOW_RATE,FLOW_TAU, &
                GAUGE_EMISSIVITY,GAUGE_TEMPERATURE,INITIAL_TEMPERATURE,K_FACTOR,C_FACTOR,OPERATING_PRESSURE,OFFSET,&
-               SPRAY_ANGLE(2,2),P0=0._EB,PX(3)=0._EB,PXX(3,3)=0._EB,VIEW_ANGLE
+               SPRAY_ANGLE(2,2),P0=0._EB,PX(3)=0._EB,PXX(3,3)=0._EB,VIEW_ANGLE,PROBE_DIAMETER
    INTEGER  :: PDPA_M=0,PDPA_N=0,N_SMOKEVIEW_PARAMETERS=0,N_SMOKEVIEW_IDS=0,N_INSERT,I_VEL=0,PARTICLES_PER_SECOND
    LOGICAL  :: PDPA_INTEGRATE=.TRUE.,PDPA_NORMALIZE=.TRUE.,HISTOGRAM_NORMALIZE=.TRUE.,HISTOGRAM=.FALSE., &
-               HISTOGRAM_CUMULATIVE=.FALSE.,SPARK=.FALSE.
+               HISTOGRAM_CUMULATIVE=.FALSE.,SPARK=.FALSE.,TC=.TRUE.
    REAL(EB) :: PDPA_START=0._EB,PDPA_END=1.E6_EB,PDPA_RADIUS=0.1_EB
    REAL(EB), ALLOCATABLE, DIMENSION(:) :: TABLE_ROW, V_FACTOR
    INTEGER  :: PART_INDEX=-1,FLOW_RAMP_INDEX,SPRAY_PATTERN_INDEX,Z_INDEX=-999,Y_INDEX=-999,PRESSURE_RAMP_INDEX

Source/dump.f90

@@ -8494,13 +8494,13 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
 INTEGER, INTENT(IN) :: II,JJ,KK,IND,IND2,NM,Y_INDEX,Z_INDEX,ELEM_INDX,PART_INDEX,VELO_INDEX,PIPE_INDEX,PROP_INDEX,REAC_INDEX, &
                        MATL_INDEX
 INTEGER, INTENT(IN), OPTIONAL :: ICC_IN,JCC_IN
-REAL(EB) :: H_TC,TMP_TC,RE_D,NUSSELT,VEL,K_G,MU_G,&
+REAL(EB) :: H_TC,TMP_TC,RE_D,NUSSELT,VEL,K_G,MU_G,COSTHETA,FAC,DP,&
             Q_SUM,TMP_G,UU,VV,WW,VEL2,Y_MF_INT,PATHLENGTH,EXT_COEF,MASS_EXT_COEF,ZZ_FUEL,ZZ_OX,&
             VELSR,WATER_VOL_FRAC,RHS,DT_C,DT_E,T_RATIO,Y_E_LAG, H_G,H_G_SUM,CPBAR,CP,ZZ_GET(1:N_TRACKED_SPECIES),RCON,&
             EXPON,Y_SPECIES,MEC,Y_SPECIES2,Y_H2O,R_Y_H2O,R_DN,SGN,Y_ALL(N_SPECIES),H_S,D_Z_N(0:I_MAX_TEMP),&
             DISSIPATION_RATE,S11,S22,S33,S12,S13,S23,DUDX,DUDY,DUDZ,DVDX,DVDY,DVDZ,DWDX,DWDY,DWDZ,ONTHDIV,SS,ETA,DELTA,R_DX2,&
             UVW,UODX,VODY,WODZ,XHAT,ZHAT,BBF,GAMMA_LOC,VC,VOL,PHI,GAS_PHASE_OUTPUT_CC,&
-            GAS_PHASE_OUTPUT_CFA,CFACE_AREA,VELOCITY_COMPONENT(1:3),ATOTV(1:3),TMP_F,R_D,MW
+            GAS_PHASE_OUTPUT_CFA,CFACE_AREA,VELOCITY_COMPONENT(1:3),ATOTV(1:3),TMP_F,R_D,MW,RHO_AIR,PROBE_TMP
 INTEGER :: N,I,J,K,NN,IL,III,JJJ,KKK,IP,JP,KP,FED_ACTIVITY,IP1,JP1,KP1,IM1,JM1,KM1,IIM1,JJM1,KKM1,NR,NS,RAM,&
            ICC,JCC,NCELL,AXIS,ICF,NFACE,JCF,JCC_LO,JCC_HI,PDPA_FORMULA,IC
 REAL(FB) :: RN
@@ -9048,6 +9048,30 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
       ENDIF
       GAS_PHASE_OUTPUT_RES = (GAS_PHASE_OUTPUT_RES - K_G*(TMP(IP,JP,KP)-TMP(II,JJ,KK))*R_DN)*0.001
 
+   CASE(114) ! BI-DIRECTIONAL PROBE
+      IF (PY%TC) THEN
+         PROBE_TMP = GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,110,IND2,Y_INDEX,Z_INDEX,ELEM_INDX,PART_INDEX,VELO_INDEX,PIPE_INDEX,&
+                                           PROP_INDEX,REAC_INDEX,MATL_INDEX,ICC_IN,JCC_IN) + TMPM
+      ELSE
+         PROBE_TMP = TMP(II,JJ,KK)
+      ENDIF
+      UU  = U(II,JJ,KK)
+      VV  = V(II,JJ,KK)
+      WW  = W(II,JJ,KK)
+      VEL2 = UU**2+VV**2+WW**2
+      DP = 0.5_EB*VEL2*RHO(II,JJ,KK)
+      COSTHETA = (UU*ORIENTATION_VECTOR(1,DV%ORIENTATION_INDEX)+VV*ORIENTATION_VECTOR(2,DV%ORIENTATION_INDEX)+&
+                  WW*ORIENTATION_VECTOR(3,DV%ORIENTATION_INDEX))/SQRT(VEL2)
+      FAC = -2.308_EB*ABS(COSTHETA)**3 + 2.533_EB*ABS(COSTHETA)**2 + 0.7847_EB*ABS(COSTHETA)
+      VEL = FAC*SQRT(VEL2)
+      ZZ_GET(1:N_TRACKED_SPECIES) = ZZ(II,JJ,KK,1:N_TRACKED_SPECIES)
+      CALL GET_VISCOSITY(ZZ_GET,MU_G,TMP(II,JJ,KK))
+      RE_D = MIN(3800._EB,MAX(40._EB,RHO(II,JJ,KK)*VEL*PY%PROBE_DIAMETER/MU_G))
+      FAC = 1.533_EB-0.001366_EB*RE_D+0.000001688_EB*RE_D**2-0.0000000009706_EB*RE_D**3+&
+             0.0000000000002555_EB*RE_D**4-2.484E-17_EB*RE_D**5
+      RHO_AIR = 350.9736_EB/PROBE_TMP !350 is 0.0288 101325/ 8.314472
+      GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*VEL*PY%CALIBRATION_CONSTANT*FAC
+
    CASE(130) ! EXTINCTION
       ZZ_GET(1:N_TRACKED_SPECIES) = ZZ(II,JJ,KK,1:N_TRACKED_SPECIES)
       ZZ_FUEL = 0._EB
@@ -9724,11 +9748,11 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
    CASE(562) ! BFZ
       GAS_PHASE_OUTPUT_RES = BFZ(II,JJ,KK)
    CASE(563) ! BFX MINUS
-      GAS_PHASE_OUTPUT_RES = BFX(II-1,JJ,KK)
+      GAS_PHASE_OUTPUT_RES = BFX(MAX(0,II-1),JJ,KK)
    CASE(564) ! BFY MINUS
-      GAS_PHASE_OUTPUT_RES = BFY(II,JJ-1,KK)
+      GAS_PHASE_OUTPUT_RES = BFY(II,MAX(0,JJ-1),KK)
    CASE(565) ! BFZ MINUS
-      GAS_PHASE_OUTPUT_RES = BFZ(II,JJ,KK-1)
+      GAS_PHASE_OUTPUT_RES = BFZ(II,JJ,MAX(0,KK-1))
 
  END SELECT IND_SELECT
 

Source/read.f90

@@ -6595,30 +6595,30 @@ SUBROUTINE READ_PROP
 USE DEVICE_VARIABLES
 USE PHYSICAL_FUNCTIONS, ONLY : SPRAY_ANGLE_DISTRIBUTION
 USE MATH_FUNCTIONS, ONLY: GET_TABLE_INDEX
-REAL(EB) :: ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E, &
+REAL(EB) :: ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E,CALIBRATION_CONSTANT, &
             HEAT_TRANSFER_COEFFICIENT,DIAMETER,DENSITY,SPECIFIC_HEAT, &
             C_FACTOR,CHARACTERISTIC_VELOCITY,ORIFICE_DIAMETER,EMISSIVITY, &
             PARTICLE_VELOCITY,FLOW_RATE,FLOW_TAU,GAUGE_EMISSIVITY,GAUGE_TEMPERATURE,INITIAL_TEMPERATURE,K_FACTOR,&
             LENGTH,SPRAY_ANGLE(2,2),OFFSET,OPERATING_PRESSURE,RTI,PDPA_START,PDPA_END,PDPA_RADIUS,MASS_FLOW_RATE,&
-            SPRAY_PATTERN_MU,SPRAY_PATTERN_BETA,HISTOGRAM_LIMITS(2),P0,PX(3),PXX(3,3),TIME_CONSTANT,VIEW_ANGLE
+            SPRAY_PATTERN_MU,SPRAY_PATTERN_BETA,HISTOGRAM_LIMITS(2),P0,PX(3),PXX(3,3),TIME_CONSTANT,VIEW_ANGLE,PROBE_DIAMETER
 INTEGER ::I,N,NN,PDPA_M,PDPA_N,PARTICLES_PER_SECOND,VELOCITY_COMPONENT,HISTOGRAM_NBINS,FED_ACTIVITY
-LOGICAL :: PDPA_INTEGRATE,PDPA_NORMALIZE,HISTOGRAM_NORMALIZE,HISTOGRAM,HISTOGRAM_CUMULATIVE,SPARK
+LOGICAL :: PDPA_INTEGRATE,PDPA_NORMALIZE,HISTOGRAM_NORMALIZE,HISTOGRAM,HISTOGRAM_CUMULATIVE,SPARK,TC
 CHARACTER(LABEL_LENGTH) :: SMOKEVIEW_ID(SMOKEVIEW_OBJECTS_DIMENSION),QUANTITY='null',PART_ID='null',FLOW_RAMP='null', &
                  SPRAY_PATTERN_TABLE='null',SPEC_ID='null',&
                  PRESSURE_RAMP='null',SMOKEVIEW_PARAMETERS(SMOKEVIEW_OBJECTS_DIMENSION), &
                  SPRAY_PATTERN_SHAPE='GAUSSIAN'
 TYPE (PROPERTY_TYPE), POINTER :: PY
 
-NAMELIST /PROP/ ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E, &
+NAMELIST /PROP/ ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E,CALIBRATION_CONSTANT,&
                 CHARACTERISTIC_VELOCITY,C_FACTOR,DENSITY,DIAMETER,EMISSIVITY,FED_ACTIVITY,FLOW_RAMP,FLOW_RATE,FLOW_TAU, &
                 GAUGE_EMISSIVITY,GAUGE_TEMPERATURE,HEAT_TRANSFER_COEFFICIENT,HISTOGRAM,HISTOGRAM_CUMULATIVE, &
                 HISTOGRAM_LIMITS,HISTOGRAM_NBINS,HISTOGRAM_NORMALIZE,ID, &
                 INITIAL_TEMPERATURE,K_FACTOR,LENGTH,MASS_FLOW_RATE,OFFSET,OPERATING_PRESSURE,ORIFICE_DIAMETER,P0,&
                 PARTICLES_PER_SECOND,PARTICLE_VELOCITY,PART_ID,PDPA_END,&
                 PDPA_INTEGRATE,PDPA_M,PDPA_N,PDPA_NORMALIZE,PDPA_RADIUS,&
-                PDPA_START,PRESSURE_RAMP,PX,PXX,QUANTITY,RTI,SMOKEVIEW_ID,SMOKEVIEW_PARAMETERS,SPARK,&
+                PDPA_START,PRESSURE_RAMP,PROBE_DIAMETER,PX,PXX,QUANTITY,RTI,SMOKEVIEW_ID,SMOKEVIEW_PARAMETERS,SPARK,&
                 SPEC_ID,SPECIFIC_HEAT,SPRAY_ANGLE,&
-                SPRAY_PATTERN_BETA,SPRAY_PATTERN_MU,SPRAY_PATTERN_SHAPE,SPRAY_PATTERN_TABLE,TIME_CONSTANT,VELOCITY_COMPONENT,&
+                SPRAY_PATTERN_BETA,SPRAY_PATTERN_MU,SPRAY_PATTERN_SHAPE,SPRAY_PATTERN_TABLE,TC,TIME_CONSTANT,VELOCITY_COMPONENT,&
                 VIEW_ANGLE
 
 ! Count the PROP lines in the input file. Note how many of these are cables.
@@ -6661,6 +6661,7 @@ SUBROUTINE READ_PROP
    PY%ALPHA_E                  = ALPHA_E
    PY%BETA_C                   = BETA_C
    PY%BETA_E                   = BETA_E
+   PY%CALIBRATION_CONSTANT     = CALIBRATION_CONSTANT
    PY%DENSITY                  = DENSITY
    PY%DIAMETER                 = DIAMETER
    PY%EMISSIVITY               = EMISSIVITY
@@ -6676,6 +6677,7 @@ SUBROUTINE READ_PROP
    PY%OFFSET                   = OFFSET
    PY%OPERATING_PRESSURE       = OPERATING_PRESSURE
    PY%PART_ID                  = PART_ID
+   PY%PROBE_DIAMETER           = PROBE_DIAMETER
    PY%QUANTITY                 = QUANTITY
    IF (PY%PART_ID/='null' .AND. PY%QUANTITY == 'null' ) PY%QUANTITY = 'NOZZLE FLOW RATE'
    PY%RTI                      = RTI
@@ -6734,6 +6736,7 @@ SUBROUTINE READ_PROP
       CALL SPRAY_ANGLE_DISTRIBUTION(PY%SPRAY_LON,PY%SPRAY_LAT,PY%SPRAY_LON_CDF,PY%SPRAY_LAT_CDF, &
                                     SPRAY_PATTERN_BETA,SPRAY_PATTERN_MU,PY%SPRAY_ANGLE,SPRAY_PATTERN_SHAPE,NDC2)
    ENDIF
+   PY%TC = TC
    PY%TIME_CONSTANT = TIME_CONSTANT
 
    ! PDPA model
@@ -6903,6 +6906,7 @@ SUBROUTINE SET_PROP_DEFAULTS
 ALPHA_E                  = 0.0_EB
 BETA_C                   = -1.0_EB
 BETA_E                   = -1.0_EB
+CALIBRATION_CONSTANT     = 0.93_EB
 DENSITY                  = 8908._EB    ! kg/m3 (Nickel)
 DIAMETER                 = 0.001       ! m
 EMISSIVITY               = 0.85_EB
@@ -6942,6 +6946,7 @@ SUBROUTINE SET_PROP_DEFAULTS
 HISTOGRAM_NORMALIZE      = .TRUE.
 PRESSURE_RAMP            = 'null'
 P0                       = -1.E10_EB
+PROBE_DIAMETER           = 0.0254_EB
 PX                       = 0._EB
 PXX                      = 0._EB
 QUANTITY                 = 'null'
@@ -6958,6 +6963,7 @@ SUBROUTINE SET_PROP_DEFAULTS
 SPRAY_PATTERN_SHAPE      = 'GAUSSIAN'
 SPRAY_PATTERN_MU         = -1._EB
 SPRAY_PATTERN_BETA       = 5.0_EB
+TC                       = .TRUE.
 TIME_CONSTANT            = -1._EB
 FED_ACTIVITY             = 2 ! light work
 VELOCITY_COMPONENT       = 0
@@ -14209,7 +14215,8 @@ SUBROUTINE READ_DEVC
           DV%QUANTITY(1)=='GAUGE HEAT FLUX GAS' .OR. &
           DV%QUANTITY(1)=='RADIANCE' .OR. &
           DV%QUANTITY(1)=='ADIABATIC SURFACE TEMPERATURE GAS' .OR. &
-          DV%QUANTITY(1)=='RADIOMETER GAS') THEN
+          DV%QUANTITY(1)=='RADIOMETER GAS' .OR. &
+          DV%QUANTITY(1)=='BI-DIRECTIONAL PROBE') THEN
          IF (DV%ORIENTATION_INDEX==0) THEN
             WRITE(MESSAGE,'(3A)')  'ERROR(887): DEVC ',TRIM(ID),' must have an ORIENTATION.'
             CALL SHUTDOWN(MESSAGE) ; RETURN

Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128.fds

@@ -26,6 +26,9 @@
 
 &SLCF PBY=0.0, QUANTITY='VELOCITY', VECTOR=.TRUE./
 &SLCF PBY=0.0, QUANTITY='MASS FRACTION', SPEC_ID='TRACER', CELL_CENTERED=T/
+&SLCF PBY=0.0, QUANTITY='BFX', CELL_CENTERED=T/
+&SLCF PBY=0.0, QUANTITY='BFY', CELL_CENTERED=T/
+&SLCF PBY=0.0, QUANTITY='BFZ', CELL_CENTERED=T/
 
 &DEVC XB=-0.00390625,-0.99609375,0,0,0.00390625,0.99609375, POINTS=128, TIME_HISTORY=T, TIME_AVERAGED=F,
       QUANTITY='MASS FRACTION', SPEC_ID='TRACER', ID='Y_TRACER'/

Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128_1mesh.fds

@@ -25,6 +25,9 @@
 
 &SLCF PBY=0.0, QUANTITY='VELOCITY', VECTOR=.TRUE./
 &SLCF PBY=0.0, QUANTITY='MASS FRACTION', SPEC_ID='TRACER', CELL_CENTERED=T/
+&SLCF PBY=0.0, QUANTITY='BFX', CELL_CENTERED=T/
+&SLCF PBY=0.0, QUANTITY='BFY', CELL_CENTERED=T/
+&SLCF PBY=0.0, QUANTITY='BFZ', CELL_CENTERED=T/
 
 &DEVC XB=-0.00390625,-0.99609375,0,0,0.00390625,0.99609375, POINTS=128, TIME_HISTORY=T, TIME_AVERAGED=F,
       QUANTITY='MASS FRACTION', SPEC_ID='TRACER', ID='Y_TRACER'/