Skip to content

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

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
marcosvanella merged 1 commit into from
Dec 12, 2025
+8 −5
Merged

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

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 8 additions & 5 deletions Source/part.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand All @@ -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:)
Expand Down
Loading