Skip to content

Commit

Permalink
FDS Source: Issue #13235. Allow printout of Y+ and U_tau
Browse files Browse the repository at this point in the history
  • Loading branch information
mcgratta committed Aug 13, 2024
1 parent 38e9b7c commit 1e40e25
Show file tree
Hide file tree
Showing 3 changed files with 5 additions and 1 deletion.
1 change: 1 addition & 0 deletions Source/cons.f90
Original file line number Diff line number Diff line change
Expand Up @@ -268,6 +268,7 @@ MODULE GLOBAL_CONSTANTS
LOGICAL :: CHECK_BOUNDARY_ONE_D_ARRAYS=.FALSE. !< Flag that indicates that ONE_D array dimensions need to be checked
LOGICAL :: TENSOR_DIFFUSIVITY=.FALSE. !< If true, use experimental tensor diffusivity model for spec and tmp
LOGICAL :: TEST_CHAR_MASS_TRANSFER_MODEL=.FALSE. !< Experimental flag to test mass transfer resistence in char model
LOGICAL :: OUTPUT_WALL_QUANTITIES=.FALSE. !< Flag to force call to WALL_MODEL

INTEGER, ALLOCATABLE, DIMENSION(:) :: CHANGE_TIME_STEP_INDEX !< Flag to indicate if a mesh needs to change time step
INTEGER, ALLOCATABLE, DIMENSION(:) :: SETUP_PRESSURE_ZONES_INDEX !< Flag to indicate if a mesh needs to keep searching for ZONEs
Expand Down
2 changes: 2 additions & 0 deletions Source/read.f90
Original file line number Diff line number Diff line change
Expand Up @@ -15835,6 +15835,8 @@ SUBROUTINE GET_QUANTITY_INDEX(SMOKEVIEW_LABEL,SMOKEVIEW_BAR_LABEL,OUTPUT_INDEX,O
ENDIF
ENDIF

IF (QUANTITY=='VISCOUS WALL UNITS' .OR. QUANTITY=='FRICTION VELOCITY') OUTPUT_WALL_QUANTITIES = .TRUE.

! Assigne MATL_INDEX when MATL_ID is specified

IF (MATL_ID/='null') THEN
Expand Down
3 changes: 2 additions & 1 deletion Source/wall.f90
Original file line number Diff line number Diff line change
Expand Up @@ -111,7 +111,8 @@ SUBROUTINE WALL_BC(T,DT,NM)
CALL CALCULATE_RHO_D_F(B1,BC,WALL_INDEX=IW)
ENDIF

IF (WC%BOUNDARY_TYPE==SOLID_BOUNDARY .AND. (ANY(SPECIES_MIXTURE%CONDENSATION_SMIX_INDEX>0).OR.DEPOSITION)) THEN
IF (WC%BOUNDARY_TYPE==SOLID_BOUNDARY .AND. &
(ANY(SPECIES_MIXTURE%CONDENSATION_SMIX_INDEX>0) .OR. DEPOSITION .OR. OUTPUT_WALL_QUANTITIES)) THEN
CALL WALL_MODEL(SLIP_COEF,B2%U_TAU,B2%Y_PLUS,MU_DNS(BC%IIG,BC%JJG,BC%KKG)/RHO(BC%IIG,BC%JJG,BC%KKG),SF%ROUGHNESS,&
0.5_EB/B1%RDN,B1%U_TANG)
ENDIF
Expand Down

0 comments on commit 1e40e25

Please sign in to comment.