Skip to content

Commit fffb5bc

Browse files
mcgrattamcgratta
committed
FDS Source: Issue #13696. Remove INITIAL_SPEED
1 parent 95782e5 commit fffb5bc

File tree

7 files changed

+9
-22
lines changed

7 files changed

+9
-22
lines changed

‎Manuals/FDS_User_Guide/FDS_User_Guide.tex‎

Lines changed: 2 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -6835,7 +6835,7 @@ \subsection{Pressure Gradient Force}
68356835

68366836
The {\ct PRESSURE\_GRADIENT\_FORCE} is input on the {\ct WIND} line in units of Pa/m, for example,
68376837
\begin{lstlisting}
6838-
&WIND PRESSURE_GRADIENT_FORCE=0.01, DIRECTION=315, INITIAL_SPEED=5 /
6838+
&WIND PRESSURE_GRADIENT_FORCE=0.01, DIRECTION=315 /
68396839
\end{lstlisting}
68406840
The wind {\ct DIRECTION} follows the usual meteorological convention---a northerly wind has direction of 0$^\circ$ and blows from north to south, or in the negative $y$ direction in the FDS coordinate system. An easterly wind has a direction of 90$^\circ$ and blows from east to west, or in the negative $x$ direction. You can vary the direction of the wind in time only, not in space:
68416841
\begin{lstlisting}
@@ -6844,7 +6844,7 @@ \subsection{Pressure Gradient Force}
68446844
&RAMP ID='dir', T= 600, F=330 /
68456845
&RAMP ID='dir', T=1200, F=350 /
68466846
\end{lstlisting}
6847-
Here, {\ct T} is time in seconds and {\ct F} is the wind direction in degrees. The {\ct INITIAL\_SPEED} and {\ct INITIAL\_DIRECTION} set the horizontal wind components at the start of the simulation. This is provided as a convenience because it may take on the order of hours of simulation time to slowly increase the wind speed via the pressure gradient force alone. The {\ct INITIAL\_SPEED} {\ct INITIAL\_DIRECTION} only set the initial wind speed and direction, but has no longer term effect. You may vary the {\ct PRESSURE\_GRADIENT\_FORCE} in time using {\ct RAMP\_PGF\_T}.
6847+
Here, {\ct T} is time in seconds and {\ct F} is the wind direction in degrees. You may vary the {\ct PRESSURE\_GRADIENT\_FORCE} in time using {\ct RAMP\_PGF\_T}.
68486848

68496849
There are some applications, like tunnels, where the pressure gradient force is a convenient way to introduce an air flow due to some external force. In such cases, it is sometimes more convenient to use the vector {\ct FORCE\_VECTOR(1:3)} along with the corresponding time ramps {\ct RAMP\_FVX\_T}, {\ct RAMP\_FVY\_T}, and {\ct RAMP\_FVZ\_T} to control each individual component. The {\ct PRESSURE\_GRADIENT\_FORCE} is simply the magnitude of the {\ct FORCE\_VECTOR} with the same units.
68506850

@@ -13675,8 +13675,6 @@ \section{\texorpdfstring{{\tt WIND}}{WIND} (Wind and Atmospheric Parameters)}
1367513675
{\ct FORCE\_VECTOR(3)} & Real & Section~\ref{info:force_vector} & Pa/m & 0. \\ \hline
1367613676
{\ct GEOSTROPHIC\_WIND(2)} & Real & Section~\ref{info:geostrophic_wind} & m/s & \\ \hline
1367713677
{\ct GROUND\_LEVEL} & Real & Section~\ref{info:stratification} & m & 0. \\ \hline
13678-
{\ct INITIAL\_DIRECTION} & Real & Section~\ref{info:force_vector} & deg & 270. \\ \hline
13679-
{\ct INITIAL\_SPEED} & Real & Section~\ref{info:force_vector} & m/s & 0. \\ \hline
1368013678
{\ct L} & Real & Section~\ref{info:WIND} & m & 0 \\ \hline
1368113679
{\ct LAPSE\_RATE} & Real & Section~\ref{info:stratification} & $^\circ$C/m & 0 \\ \hline
1368213680
{\ct LATITUDE} & Real & Section~\ref{info:coriolis_force} & degrees & \\ \hline

‎Source/cons.f90‎

Lines changed: 0 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -308,8 +308,6 @@ MODULE GLOBAL_CONSTANTS
308308
REAL(EB) :: U0 !< Wind speed in the \f$ x \f$ direction (m/s)
309309
REAL(EB) :: V0 !< Wind speed in the \f$ y \f$ direction (m/s)
310310
REAL(EB) :: W0 !< Wind speed in the \f$ z \f$ direction (m/s)
311-
REAL(EB) :: INITIAL_SPEED=-1._EB !< Initial wind speed (m/s) which is assumed to die off
312-
REAL(EB) :: INITIAL_DIRECTION=270._EB !< Initial wind direction (deg)
313311
REAL(EB) :: GVEC(3) !< Gravity vector (m/s2)
314312
REAL(EB) :: FVEC(3)=0._EB !< Force vector (N/m3)
315313
REAL(EB) :: OVEC(3)=0._EB !< Coriolis vector (1/s)

‎Source/init.f90‎

Lines changed: 3 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -363,20 +363,11 @@ SUBROUTINE INITIALIZE_MESH_VARIABLES_1(DT,NM)
363363
M%RHO(:,:,K) = M%RHO_0(K)
364364
M%RHOS(:,:,K) = M%RHO_0(K)
365365
M%TMP(:,:,K) = M%TMP_0(K)
366+
M%U(:,:,K) = M%U_WIND(K)
367+
M%V(:,:,K) = M%V_WIND(K)
368+
M%W(:,:,K) = M%W_WIND(K)
366369
ENDDO
367370

368-
IF (INITIAL_SPEED<0._EB) THEN
369-
DO K=0,M%KBP1
370-
M%U(:,:,K) = M%U_WIND(K)
371-
M%V(:,:,K) = M%V_WIND(K)
372-
M%W(:,:,K) = M%W_WIND(K)
373-
ENDDO
374-
ELSE
375-
M%U = -INITIAL_SPEED*SIN(INITIAL_DIRECTION*DEG2RAD)
376-
M%V = -INITIAL_SPEED*COS(INITIAL_DIRECTION*DEG2RAD)
377-
M%W = 0._EB
378-
ENDIF
379-
380371
! Custom velocity RAMPS (for verification)
381372

382373
IF (I_RAMP_UX>0) THEN

‎Source/read.f90‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -2077,7 +2077,7 @@ SUBROUTINE READ_WIND
20772077
TYPE(RESERVED_RAMPS_TYPE), POINTER :: RRP,RRPX
20782078
INTEGER, PARAMETER :: N_MO_PTS=51 ! number of Monin-Obukhov ramp points
20792079

2080-
NAMELIST /WIND/ CORIOLIS_VECTOR,DIRECTION,FORCE_VECTOR,FYI,GEOSTROPHIC_WIND,GROUND_LEVEL,INITIAL_DIRECTION,INITIAL_SPEED,&
2080+
NAMELIST /WIND/ CORIOLIS_VECTOR,DIRECTION,FORCE_VECTOR,FYI,GEOSTROPHIC_WIND,GROUND_LEVEL,&
20812081
L,LAPSE_RATE,LATITUDE,PRESSURE_GRADIENT_FORCE,RAMP_DIRECTION_T,RAMP_DIRECTION_Z,&
20822082
RAMP_PGF_T,RAMP_FVX_T,RAMP_FVY_T,RAMP_FVZ_T,RAMP_SPEED_T,RAMP_SPEED_Z,RAMP_TMP0_Z,&
20832083
SIGMA_THETA,SPEED,STRATIFICATION,TAU_THETA,THETA_STAR,TMP_REF,U_STAR,U0,&

‎Validation/Loughborough_Jet_Fires/FDS_Input_Files/jet1.fds‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -14,7 +14,7 @@
1414
&VENT MB='XMIN', SURF_ID='OPEN' /
1515
&VENT MB='XMAX', SURF_ID='OPEN' /
1616

17-
&WIND INITIAL_SPEED=6.3, INITIAL_DIRECTION=271., DIRECTION=271., PRESSURE_GRADIENT_FORCE=0.08 /
17+
&WIND DIRECTION=271., PRESSURE_GRADIENT_FORCE=0.08 /
1818

1919
&SURF ID='GROUND', ROUGHNESS=0.5 /
2020

‎Validation/Loughborough_Jet_Fires/FDS_Input_Files/jet2.fds‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -14,7 +14,7 @@
1414
&VENT MB='XMIN', SURF_ID='OPEN' /
1515
&VENT MB='XMAX', SURF_ID='OPEN' /
1616

17-
&WIND INITIAL_SPEED=6.2, INITIAL_DIRECTION=297., DIRECTION=297., PRESSURE_GRADIENT_FORCE=0.08 /
17+
&WIND DIRECTION=297., PRESSURE_GRADIENT_FORCE=0.08 /
1818

1919
&SURF ID='GROUND', ROUGHNESS=0.5 /
2020

‎Validation/Loughborough_Jet_Fires/FDS_Input_Files/jet3.fds‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -14,7 +14,7 @@
1414
&VENT MB='XMIN', SURF_ID='OPEN' /
1515
&VENT MB='XMAX', SURF_ID='OPEN' /
1616

17-
&WIND INITIAL_SPEED=3.6, INITIAL_DIRECTION=267., DIRECTION=267., PRESSURE_GRADIENT_FORCE=0.08 /
17+
&WIND DIRECTION=267., PRESSURE_GRADIENT_FORCE=0.08 /
1818

1919
&SURF ID='GROUND', ROUGHNESS=0.5 /
2020

0 commit comments

Comments
 (0)