diff --git a/Manuals/FDS_User_Guide/FDS_User_Guide.tex b/Manuals/FDS_User_Guide/FDS_User_Guide.tex index 7b8be93228e..383006ab3df 100644 --- a/Manuals/FDS_User_Guide/FDS_User_Guide.tex +++ b/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} \\ diff --git a/Source/data.f90 b/Source/data.f90 index eccbe2b8c18..b46fc0e5b66 100644 --- a/Source/data.f90 +++ b/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' diff --git a/Source/devc.f90 b/Source/devc.f90 index cfafe8670c4..e7de967e69c 100644 --- a/Source/devc.f90 +++ b/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 diff --git a/Source/dump.f90 b/Source/dump.f90 index 584cd6ce18e..5eafda355dd 100644 --- a/Source/dump.f90 +++ b/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 diff --git a/Source/read.f90 b/Source/read.f90 index c30fc959927..f792194b30a 100644 --- a/Source/read.f90 +++ b/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 diff --git a/Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128.fds b/Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128.fds index d81149bf637..9568df5d5d0 100644 --- a/Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128.fds +++ b/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'/ diff --git a/Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128_1mesh.fds b/Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128_1mesh.fds index df6788acec7..6181438e904 100644 --- a/Verification/Scalar_Analytical_Solution/soborot_superbee_square_wave_128_1mesh.fds +++ b/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'/