FDS Source: Allow for VEL_T with VEL=0

Author: mcgratta

Source/velo.f90

@@ -1816,7 +1816,7 @@ SUBROUTINE VELOCITY_BC(T,NM,APPLY_TO_ESTIMATED_VARIABLES)
 REAL(EB) :: MUA,TSI,WGT,T_NOW,RAMP_T,OMW,MU_WALL,RHO_WALL,SLIP_COEF,VEL_T, &
             UUP(2),UUM(2),DXX(2),MU_DUIDXJ(-2:2),DUIDXJ(-2:2),PROFILE_FACTOR,VEL_GAS,VEL_GHOST, &
             MU_DUIDXJ_USE(2),DUIDXJ_USE(2),VEL_EDDY,U_TAU,Y_PLUS,U_NORM, &
-            DRAG_FACTOR,HT_SCALE_FACTOR,VEG_HT,VEL_MF
+            DRAG_FACTOR,HT_SCALE_FACTOR,VEG_HT,VEL_N
 INTEGER :: NOM(2),IIO(2),JJO(2),KKO(2),IE,II,JJ,KK,IEC,IOR,IWM,IWP,ICMM,ICMP,ICPM,ICPP,ICD,ICDO,IVL,I_SGN, &
            VELOCITY_BC_INDEX,IIGM,JJGM,KKGM,IIGP,JJGP,KKGP,SURF_INDEXM,SURF_INDEXP,ITMP,ICD_SGN,ICDO_SGN, &
            BOUNDARY_TYPE_M,BOUNDARY_TYPE_P,IS,IS2,IWPI,IWMI,VENT_INDEX
@@ -2267,29 +2267,33 @@ SUBROUTINE VELOCITY_BC(T,NM,APPLY_TO_ESTIMATED_VARIABLES)
 
             ELSE
 
-               PROFILE_FACTOR = 1._EB
-               IF (ABS(SF%VEL)>0._EB) THEN
+               VEL_N = 0.5_EB*(WCM_B1%U_NORMAL_S+WCP_B1%U_NORMAL_S)
+
+               IF (ABS(SF%VEL)>0._EB .OR. VEL_N==0._EB) THEN  ! User-specified normal velocity or no normal velocity at all
                   IF (ABS(SF%T_IGN-T_BEGIN)<=SPACING(SF%T_IGN) .AND. SF%RAMP(TIME_VELO)%INDEX>=1) THEN
                      TSI = T
                   ELSE
                      TSI=T-SF%T_IGN
                   ENDIF
+                  PROFILE_FACTOR = 1._EB
                   RAMP_T = EVALUATE_RAMP(TSI,SF%RAMP(TIME_VELO)%INDEX,TAU=SF%RAMP(TIME_VELO)%TAU)
                   IF (SF%VEL < 0._EB) THEN
                      IF (SF%RAMP(VELO_PROF_Z)%INDEX>0) PROFILE_FACTOR = EVALUATE_RAMP(ZC(KK),SF%RAMP(VELO_PROF_Z)%INDEX)
                      IF (IEC==1 .OR. (IEC==2 .AND. ICD==2)) VEL_T = RAMP_T*(PROFILE_FACTOR*(SF%VEL_T(2) + VEL_EDDY))
                      IF (IEC==3 .OR. (IEC==2 .AND. ICD==1)) VEL_T = RAMP_T*(PROFILE_FACTOR*(SF%VEL_T(1) + VEL_EDDY))
-                  ELSE
+                  ELSEIF (SF%VEL > 0._EB) THEN
                      IF (SF%PROFILE/=0) PROFILE_FACTOR = ABS(0.5_EB*(WCM_B1%U_NORMAL_0+WCP_B1%U_NORMAL_0)/SF%VEL)
                      IF (SF%RAMP(VELO_PROF_Z)%INDEX>0) PROFILE_FACTOR = EVALUATE_RAMP(ZC(KK),SF%RAMP(VELO_PROF_Z)%INDEX)
                      IF (IEC==1 .OR. (IEC==2 .AND. ICD==2)) VEL_T = RAMP_T*PROFILE_FACTOR*VEL_EDDY
                      IF (IEC==3 .OR. (IEC==2 .AND. ICD==1)) VEL_T = RAMP_T*PROFILE_FACTOR*VEL_EDDY
+                  ELSE  ! User-specified VEL_T but with VEL=0
+                     IF (IEC==1 .OR. (IEC==2 .AND. ICD==2)) VEL_T = RAMP_T*SF%VEL_T(2)
+                     IF (IEC==3 .OR. (IEC==2 .AND. ICD==1)) VEL_T = RAMP_T*SF%VEL_T(1)
                   ENDIF
-               ELSE ! Mass flux BC
-                  VEL_MF = 0.5_EB*(WCM_B1%U_NORMAL_S+WCP_B1%U_NORMAL_S)
-                  IF (VEL_MF < 0._EB) THEN
-                     IF (IEC==1 .OR. (IEC==2 .AND. ICD==2)) VEL_T = -SF%VEL_T(2)*VEL_MF
-                     IF (IEC==3 .OR. (IEC==2 .AND. ICD==1)) VEL_T = -SF%VEL_T(1)*VEL_MF
+               ELSE ! VEL_N is due to something else besides a user-specified VEL, like a MASS_FLUX BC
+                  IF (VEL_N < 0._EB) THEN
+                     IF (IEC==1 .OR. (IEC==2 .AND. ICD==2)) VEL_T = -SF%VEL_T(2)*VEL_N
+                     IF (IEC==3 .OR. (IEC==2 .AND. ICD==1)) VEL_T = -SF%VEL_T(1)*VEL_N
                   ENDIF
                ENDIF
 

Verification/Flowfields/lid_driven_cavity_2d.fds

@@ -7,7 +7,7 @@
 &SPEC ID='LJ AIR', VISCOSITY=1., BACKGROUND=.TRUE./
 
 &SURF ID='bot', COLOR='BLUE' /
-&SURF ID='top', COLOR='RED', VEL=-0, VEL_T(1)=1., TAU_V=0. /
+&SURF ID='top', COLOR='RED', VEL_T(1)=1., TAU_V=0. /
 
 &VENT PBX=0, SURF_ID='PERIODIC'/
 &VENT PBX=1, SURF_ID='PERIODIC'/

Verification/Flowfields/uniform_shear_dudy.fds

@@ -7,7 +7,7 @@
 &SPEC ID='LJ AIR', VISCOSITY=1., BACKGROUND=.TRUE./
 
 &SURF ID='bot', COLOR='BLUE' /
-&SURF ID='top', COLOR='RED', VEL=-0, VEL_T(1)=1., TAU_V=0. /
+&SURF ID='top', COLOR='RED', VEL_T(1)=1., TAU_V=0. /
 
 &VENT PBX=0, SURF_ID='PERIODIC'/
 &VENT PBX=1, SURF_ID='PERIODIC'/

Verification/Flowfields/uniform_shear_dudz.fds

@@ -7,7 +7,7 @@
 &SPEC ID='LJ AIR', VISCOSITY=1., BACKGROUND=.TRUE./
 
 &SURF ID='bot', COLOR='BLUE' /
-&SURF ID='top', COLOR='RED', VEL=-0, VEL_T(2)=1., TAU_V=0. /
+&SURF ID='top', COLOR='RED', VEL_T(2)=1., TAU_V=0. /
 
 &VENT PBX=0, SURF_ID='MIRROR'/
 &VENT PBX=1, SURF_ID='MIRROR'/