FDS Source: add TANGENTIAL VELOCITY output for B1%U_TANG

Author: ericvmueller

Manuals/FDS_User_Guide/FDS_User_Guide.tex

@@ -2271,7 +2271,7 @@ \subsubsection{Default Convective Heat Transfer Model}
 \end{center}
 \end{table}
 
-The length scale, $L$, is specified by \ct{CONVECTION_LENGTH_SCALE} on the \ct{SURF} line. By default, it is 1~m for plates, and the diameter of a sphere or cylinder.
+The length scale, $L$, is specified by \ct{CONVECTION_LENGTH_SCALE} on the \ct{SURF} line. By default, it is 1~m for plates, and the diameter of a sphere or cylinder. The velocity used to compute Re can be output with the solid phase quantity \ct{'TANGENTIAL VELOCITY'}.
 
 \subsubsection{Output for Convective Heat Transfer Regime}
 \label{info:convection_regime}
@@ -11737,6 +11737,7 @@ \section{Solid Phase Output Quantities}
 \ct{SUBSTEPS}                                  & Section~\ref{info:solid_phase_stability}      &                    & B,D          \\ \hline
 \ct{SURFACE DENSITY}$^4$                       & Section~\ref{info:material_components}        & kg/m$^2$           & B,D          \\ \hline
 \ct{SURFACE DEPOSITION}$^1$                    & Section~\ref{info:deposition}                 & kg/m$^2$           & B,D          \\ \hline
+\ct{TANGENTIAL VELOCITY}                       & Section~\ref{info:convection}                 & m/s                & B,D          \\ \hline
 \ct{TOTAL MASS FLUX WALL}$^1$                  & Section~\ref{info:wallflux}                   & kg/s/m$^2$         & B,D          \\ \hline
 \ct{VELOCITY ERROR}                            & Section~\ref{info:PRES}                       & m/s                & B,D          \\ \hline
 \ct{VISCOUS STRESS WALL}                       & Section~\ref{info:distributed_forces}         & Pa                 & B,D          \\ \hline

Source/data.f90

@@ -1752,6 +1752,10 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 OUTPUT_QUANTITY(-100)%SPEC_ID_REQUIRED = .FALSE.
 OUTPUT_QUANTITY(-100)%BNDF_APPROPRIATE = .TRUE.
 
+OUTPUT_QUANTITY(-101)%NAME = 'TANGENTIAL VELOCITY'
+OUTPUT_QUANTITY(-101)%UNITS= 'm/s'
+OUTPUT_QUANTITY(-101)%SHORT_NAME = 'u_tang'
+
 END SUBROUTINE DEFINE_OUTPUT_QUANTITIES
 
 

Source/dump.f90

@@ -9538,6 +9538,9 @@ REAL(EB) FUNCTION SOLID_PHASE_OUTPUT(INDX,Y_INDEX,Z_INDEX,PART_INDEX,OPT_WALL_IN
    CASE(100) ! CONDENSATION HEAT FLUX
       SOLID_PHASE_OUTPUT = B1%Q_CONDENSE * 0.001_EB
 
+   CASE(101) ! TANGENTIAL VELOCITY
+      IF (ASSOCIATED(B1)) SOLID_PHASE_OUTPUT = B1%U_TANG
+
 END SELECT SOLID_PHASE_SELECT
 
 END FUNCTION SOLID_PHASE_OUTPUT