Skip to content

FDS Source: Initialize AIT zone temp Issue #15175 #15182

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
drjfloyd merged 2 commits into from
Sep 23, 2025
Merged

FDS Source: Initialize AIT zone temp Issue #15175 #15182

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
c4a789b
FDS Source: Initialize AIT zone temp Issue #15175
drjfloyd Sep 23, 2025
e88f8b0
FDS Source: remove debug
drjfloyd Sep 23, 2025
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
29 changes: 15 additions & 14 deletions Source/read.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4790,6 +4790,7 @@ SUBROUTINE READ_REAC
RN%AIT_EXCLUSION_ZONE(IZ)%CTRL_ID = AIT_EXCLUSION_ZONE_CTRL_ID(IZ)
ENDIF
ENDDO

RN%E = E*1000._EB
RN%E_IN = E
RN%EQUATION = EQUATION
Expand Down Expand Up @@ -5046,6 +5047,7 @@ SUBROUTINE SET_REAC_DEFAULTS
A = -1._EB
AUTO_IGNITION_TEMPERATURE = -TMPM
AIT_EXCLUSION_ZONE = -1.E6_EB
AIT_EXCLUSION_ZONE_TEMPERATURE = -TMPM
AIT_EXCLUSION_ZONE_CTRL_ID = 'null'
AIT_EXCLUSION_ZONE_DEVC_ID = 'null'
C = 0._EB
Expand Down Expand Up @@ -5542,13 +5544,13 @@ SUBROUTINE PROC_REAC_1
ELSEIF(TRIM(SPECIES_MIXTURE(NS)%ID) == "O2" .OR. TRIM(SPECIES_MIXTURE(NS)%ID) == "OXYGEN") THEN
I_O2 = NS
ELSEIF(TRIM(SPECIES_MIXTURE(NS)%ID) == "CO2" .OR. TRIM(SPECIES_MIXTURE(NS)%ID) == "CARBON DIOXIDE") THEN
I_CO2 = NS
I_CO2 = NS
ELSEIF(TRIM(SPECIES_MIXTURE(NS)%ID) == "H2O" .OR. TRIM(SPECIES_MIXTURE(NS)%ID) == "WATER VAPOR") THEN
I_H2O = NS
ENDIF
I_H2O = NS
ENDIF
ENDDO

ENDIF
ENDIF

END SUBROUTINE PROC_REAC_1

Expand Down Expand Up @@ -6605,7 +6607,7 @@ SUBROUTINE READ_PROP
PY%SPECIFIC_HEAT(J) = 451.5_EB+0.001_EB*(REAL(MAX(293,MIN(1493,J)),EB)-293._EB)*(601.0_EB-451.5_EB)
ENDDO
ENDIF
PY%SPECIFIC_HEAT = PY%SPECIFIC_HEAT/TIME_SHRINK_FACTOR
PY%SPECIFIC_HEAT = PY%SPECIFIC_HEAT/TIME_SHRINK_FACTOR
IF (SPECIFIC_HEAT_RAMP /= 'null') CALL GET_RAMP_INDEX(SPECIFIC_HEAT_RAMP,'TEMPERATURE',PY%SPECIFIC_HEAT_RAMP_INDEX)
PY%C_FACTOR = C_FACTOR
PY%CHARACTERISTIC_VELOCITY = CHARACTERISTIC_VELOCITY
Expand Down Expand Up @@ -6852,7 +6854,7 @@ SUBROUTINE SET_PROP_DEFAULTS
DIAMETER = 0.001 ! m
EMISSIVITY = 0.85_EB
HEAT_TRANSFER_COEFFICIENT= -1._EB ! W/m2/K
SPECIFIC_HEAT = -1._EB ! kJ/kg/K
SPECIFIC_HEAT = -1._EB ! kJ/kg/K
SPECIFIC_HEAT_RAMP = 'null'
C_FACTOR = 0.0_EB
CHARACTERISTIC_VELOCITY = 1.0_EB ! m/s
Expand Down Expand Up @@ -7019,7 +7021,7 @@ SUBROUTINE PROC_PROP
PY%V_FACTOR = PY%PARTICLE_VELOCITY/SQRT(PY%OPERATING_PRESSURE)
ENDIF
ENDIF

! Check units of specific heat

IF (PY%SPECIFIC_HEAT_RAMP_INDEX > 0) THEN
Expand Down Expand Up @@ -8857,8 +8859,8 @@ SUBROUTINE READ_SURF(QUICK_READ)
IF (DELAMINATION_TMP(NL) > 0._EB) SF%DELAMINATION_TMP(NL) = DELAMINATION_TMP(NL) + TMPM
IF (DELAMINATION_DENSITY(NL) > 0._EB) SF%DELAMINATION_DENSITY(NL) = DELAMINATION_DENSITY(NL)
ENDDO


! Adjust min and max temperature and convert C to K

IF (TMP_FRONT >= -TMPM) TMPMIN = MIN(TMPMIN,TMP_FRONT+TMPM)
Expand Down Expand Up @@ -13700,7 +13702,7 @@ SUBROUTINE READ_DEVC
IF (PROP_ID/='null') THEN
CALL GET_PROPERTY_INDEX(PROP_INDEX,'DEVC',PROP_ID)
IF (PROPERTY(PROP_INDEX)%IGNITION_ZONE) XB_XYZ= .FALSE.
ENDIF
ENDIF
IF (POINTS==1 .AND. SPATIAL_STATISTIC=='null' .AND. ALL(XB>-1E6_EB) .AND. ALL(XYZ<=-1E6_EB) .AND. XB_XYZ) THEN
XYZ(1) = 0.5_EB*(XB(1)+XB(2))
XYZ(2) = 0.5_EB*(XB(3)+XB(4))
Expand Down Expand Up @@ -14929,7 +14931,7 @@ SUBROUTINE PROC_DEVC
WRITE(MESSAGE,'(3A)') 'ERROR(899): DEVC ',TRIM(DV%ID),' should not contain a + sign.'
CALL SHUTDOWN(MESSAGE) ; RETURN
ENDIF

! Check for HVAC outputs with no HVAC inputs

IF ((DV%DUCT_ID/='null' .OR. DV%NODE_ID(1)/='null') .AND. .NOT. HVAC_SOLVE) THEN
Expand Down Expand Up @@ -15006,7 +15008,7 @@ SUBROUTINE PROC_DEVC
ENDIF
IGNITION_ZONES(N_IGNITION_ZONES)%DEVC_INDEX = N
IGNITION_ZONES(N_IGNITION_ZONES)%DEVC_ID = DV%ID
ENDIF
ENDIF
ENDIF

IF(TRIM(ODE_SOLVER)=='CVODE' .AND. N_IGNITION_ZONES >= 1) THEN
Expand Down Expand Up @@ -15228,7 +15230,6 @@ SUBROUTINE PROC_DEVC
2._EB*PY%DENSITY*PY%SPECIFIC_HEAT(293)*PY%DIAMETER
PY%DIAMETER = PY%DIAMETER - F/DFDD
TOL = ABS(F/DFDD)
WRITE(*,*)'C',TRIM(PY%ID),PY%DIAMETER
ENDDO
ENDIF
ENDIF
Expand Down Expand Up @@ -15526,7 +15527,7 @@ SUBROUTINE READ_PROF
CALL SHUTDOWN(MESSAGE) ; RETURN
ENDIF
ENDIF

IF (INIT_ID/='null') THEN
INIT_DO: DO I=1,N_INIT
IF (INIT_ID==INITIALIZATION(I)%ID) THEN
Expand Down
Loading