Merge pull request #15714 from marcosvanella/master

FDS Source : Add VEL_PART as option for particle ejection velocity in…

Author: marcosvanella

Source/part.f90

@@ -857,9 +857,12 @@ SUBROUTINE PARTICLE_FACE_INSERT(WALL_INDEX,CFACE_INDEX)
                BC%Z = CFA_Z + SF%EMBER_GENERATION_HEIGHT(1) + &
                       (SF%EMBER_GENERATION_HEIGHT(2)-SF%EMBER_GENERATION_HEIGHT(1))*REAL(RN3,EB)
             ENDIF
-            LP%U = DOT_PRODUCT(CFA_BC%NVEC,(/-B1%U_NORMAL,SF%VEL_T(1),SF%VEL_T(2)/))
-            LP%V = DOT_PRODUCT(CFA_BC%NVEC,(/SF%VEL_T(1),-B1%U_NORMAL,SF%VEL_T(2)/))
-            LP%W = DOT_PRODUCT(CFA_BC%NVEC,(/SF%VEL_T(1),SF%VEL_T(2),-B1%U_NORMAL/))
+            VEL_PART = B1%U_NORMAL
+            IF (.NOT.LPC%STATIC .AND. SF%VEL_PART >-999999._EB) VEL_PART = SF%VEL_PART
+            LP%U = DOT_PRODUCT(CFA_BC%NVEC,(/-VEL_PART,SF%VEL_T(1),SF%VEL_T(2)/))
+            LP%V = DOT_PRODUCT(CFA_BC%NVEC,(/SF%VEL_T(1),-VEL_PART,SF%VEL_T(2)/))
+            LP%W = DOT_PRODUCT(CFA_BC%NVEC,(/SF%VEL_T(1),SF%VEL_T(2),-VEL_PART/))
+
          ENDIF WALL_OR_CFACE_IF_2
 
          ! Update idicies in case offset puts location in a different cell
@@ -3608,7 +3611,7 @@ SUBROUTINE PARTICLE_MASS_ENERGY_TRANSFER(T,DT,NM)
 IF (.NOT. PARTICLES_EXISTED .OR. ((.NOT. LIQUID_DROPLETS) .AND. (.NOT. SOLID_PARTICLES))) THEN
    ! No particles ever existed in the mesh or only MASSLESS defined, then there is nothing to do.
    T_USED(8)=T_USED(8)+CURRENT_TIME()-TNOW
-   RETURN   
+   RETURN
 ENDIF
 
 IF (LIQUID_DROPLETS) THEN
@@ -3632,7 +3635,7 @@ SUBROUTINE PARTICLE_MASS_ENERGY_TRANSFER(T,DT,NM)
       B2 => BOUNDARY_PROP2(CFA%B2_INDEX)
       B2%WORK1 = B1%TMP_F
    ENDDO
-   
+
    RHO_INTERIM => WORK1 ; RHO_INTERIM(0:IBP1,0:JBP1,0:KBP1) = RHO(0:IBP1,0:JBP1,0:KBP1)
    TMP_INTERIM => WORK2 ; TMP_INTERIM(0:IBP1,0:JBP1,0:KBP1) = TMP(0:IBP1,0:JBP1,0:KBP1)
    ZZ_INTERIM  => SWORK1 ; ZZ_INTERIM(0:IBP1,0:JBP1,0:KBP1,1:) = ZZ(0:IBP1,0:JBP1,0:KBP1,1:)