Skip to content

FDS Source : CC_IBM, initialize pressure solver considering number of… #16127

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
Apr 7, 2026
+15 −5
Merged

FDS Source : CC_IBM, initialize pressure solver considering number of… #16127

Show file tree
Hide file tree
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
6 changes: 1 addition & 5 deletions Source/geom.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -23839,11 +23839,7 @@ SUBROUTINE READ_GEOM

CC_IBM = .TRUE.

! If unstructured projection defined set Pressure solver on unstructured grid.
IF (PRES_FLAG/=UGLMAT_FLAG) THEN
PRES_METHOD = 'ULMAT'
PRES_FLAG = ULMAT_FLAG
ENDIF
! For unstructured projection define Pressure solver on unstructured grid.
PRES_ON_WHOLE_DOMAIN = .FALSE.
IF (ABS(CCVOL_LINK-0.95_EB)<TWENTY_EPSILON_EB) CCVOL_LINK = DEFAULT_VOLFRAC_LINK
IF (CHECK_POISSON) GLMAT_VERBOSE=.TRUE.
Expand Down
12 changes: 12 additions & 0 deletions Source/init.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -996,6 +996,7 @@ SUBROUTINE INITIALIZE_MESH_VARIABLES_2(NM)
USE MEMORY_FUNCTIONS, ONLY: REALLOCATE_EDGE,REALLOCATE_REAL_ARRAY
USE PHYSICAL_FUNCTIONS, ONLY: GET_SPECIFIC_HEAT
USE CONTROL_VARIABLES
USE MISC_FUNCTIONS, ONLY: DEFINE_PRES_METHOD
INTEGER :: N,I,J,K,IPTS,JPTS,KPTS,N_EDGES_DIM,IW,IC,IERR,IPZ,IZERO,ICF,NSLICE
INTEGER, INTENT(IN) :: NM
REAL(EB) :: ZZ_GET(1:N_TRACKED_SPECIES),VC,RTRM,CP
Expand Down Expand Up @@ -1125,6 +1126,17 @@ SUBROUTINE INITIALIZE_MESH_VARIABLES_2(NM)
B1%PRESSURE_ZONE = M%PRESSURE_ZONE(BC%IIG,BC%JJG,BC%KKG)
ENDDO CFACE_LOOP_0

! For CC_IBM and unspecified pressure solver, select Pressure solver
! with number of pressure zones:
IF(CC_IBM .AND. PERIODIC_TEST/=7) THEN
IF (PRES_FLAG/=UGLMAT_FLAG .AND. N_ZONE>1) THEN
PRES_METHOD = 'ULMAT'
ELSEIF (PRES_FLAG/=ULMAT_FLAG .AND. N_ZONE<=1) THEN
PRES_METHOD = 'UGLMAT'
ENDIF
CALL DEFINE_PRES_METHOD(PRES_METHOD)
ENDIF

! If there is complex terrain using GEOM and Above Ground Level (AGL) slices,
! determine K index of gas phase quantities.

Expand Down
2 changes: 2 additions & 0 deletions Verification/Complex_Geometry/geom_stretched_grid.fds
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,6 +11,8 @@
&TIME T_END=1200/
&DUMP DT_DEVC=600/

&PRES SOLVER='ULMAT' /

&GEOM XB=-1200,-800,-1600,-1200,400,800/ GEOM on stretch mesh

&WIND SPEED=10.0, DIRECTION=206 /
Expand Down
Loading