Skip to content

Commit e029e8a

Browse files
mcgrattamcgratta
authored
Merge pull request #14098 from mcgratta/master
FDS Source: Issue #13516. LAYER_DIVIDE_DEPTH instead of PYROLYSIS DEPTH
2 parents 1836b87 + 89ca851 commit e029e8a

File tree

6 files changed

+13
-13
lines changed

6 files changed

+13
-13
lines changed

‎Manuals/FDS_User_Guide/FDS_User_Guide.tex‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -2875,9 +2875,9 @@ \subsection{Solid Phase Gas Transport}
28752875

28762876
Suppose, for example, that the solid is composed of a layer of insulation on top of a layer of plastic on top of a layer of steel. The steel is impermeable. By setting \ct{LAYER_DIVIDE=2.0} to the \ct{SURF} line whose first layer is the insulation directs the vapors generated by the insulation and plastic to be driven out of the exterior surface of the insulation. Similarly, for the \ct{SURF} line that is applied to the steel, specifying \ct{LAYER_DIVIDE=0.0} would indicate that no fuel vapors are to escape the steel surface. Note that in this instance, the sum of the values of \ct{LAYER_DIVIDE} is not equal to the number of layers, but this is not a problem because the layer of steel does not generate any gases.
28772877

2878-
The solid phase output \ct{QUANTITY='PYROLYSIS DEPTH'} can check if the \ct{LAYER_DIVIDE} is being applied correctly. For example, the input line
2878+
The solid phase output \ct{QUANTITY='LAYER DIVIDE DEPTH'} can check if the \ct{LAYER_DIVIDE} is being applied correctly. For example, the input line
28792879
\begin{lstlisting}
2880-
&BNDF QUANTITY='PYROLYSIS DEPTH', CELL_CENTERED=T /
2880+
&BNDF QUANTITY='LAYER DIVIDE DEPTH', CELL_CENTERED=T /
28812881
\end{lstlisting}
28822882
shows the distance from the surface where the pyrolyates are directed towards that surface, as opposed to the opposite side of the solid.
28832883

@@ -11492,6 +11492,7 @@ \section{Solid Phase Output Quantities}
1149211492
\ct{INCIDENT HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
1149311493
\ct{INSIDE WALL TEMPERATURE} & Section~\ref{info:DEPTH} & $^\circ$C & D,Pr \\ \hline
1149411494
\ct{INSIDE WALL DEPTH} & Section~\ref{info:DEPTH} & m & D,Pr \\ \hline
11495+
\ct{LAYER DIVIDE DEPTH} & Section~\ref{info:LAYER_DIVIDE} & m & B,D \\ \hline
1149511496
\ct{MASS FLUX}$^{1,4}$ & Section~\ref{info:wallflux} & \unit{kg/(m^2.s)} & B,D \\ \hline
1149611497
\ct{MASS FLUX WALL}$^1$ & Section~\ref{info:wallflux} & \unit{kg/(m^2.s)} & B,D \\ \hline
1149711498
\ct{MPUA}$^2$ & Section~\ref{bucket_test_1} & kg/m$^2$ & B,D \\ \hline
@@ -11503,7 +11504,6 @@ \section{Solid Phase Output Quantities}
1150311504
\ct{NORMALIZED MASS LOSS RATE}$^4$ & Section~\ref{info:material_components} & 1/s & D \\ \hline
1150411505
\ct{OXIDATIVE HRRPUA} & Section~\ref{sec:reaction_rates} & \unit{kW/m^2} & B,D \\ \hline
1150511506
\ct{PRESSURE COEFFICIENT} & Section~\ref{info:pressure_coefficient} & & B,D \\ \hline
11506-
\ct{PYROLYSIS DEPTH} & Section~\ref{info:LAYER_DIVIDE} & m & B,D \\ \hline
1150711507
\ct{RADIATIVE HEAT FLUX} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
1150811508
\ct{RADIOMETER} & Section~\ref{info:heat_flux} & \unit{kW/m^2} & B,D \\ \hline
1150911509
\ct{REFERENCE_HEAT_FLUX} & Section~\ref{info:scaled_burning} & \unit{kW/m^2} & B,D \\ \hline

‎Source/data.f90‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1509,9 +1509,9 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
15091509
OUTPUT_QUANTITY(-47)%INSIDE_SOLID = .TRUE.
15101510
OUTPUT_QUANTITY(-47)%BNDF_APPROPRIATE = .FALSE.
15111511

1512-
OUTPUT_QUANTITY(-48)%NAME= 'PYROLYSIS DEPTH'
1512+
OUTPUT_QUANTITY(-48)%NAME= 'LAYER DIVIDE DEPTH'
15131513
OUTPUT_QUANTITY(-48)%UNITS= 'm'
1514-
OUTPUT_QUANTITY(-48)%SHORT_NAME= 'p-depth'
1514+
OUTPUT_QUANTITY(-48)%SHORT_NAME= 'ld-depth'
15151515

15161516
OUTPUT_QUANTITY(-51)%NAME = 'ENTHALPY FLUX WALL'
15171517
OUTPUT_QUANTITY(-51)%UNITS= 'kW/m2'

‎Source/dump.f90‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -8982,9 +8982,9 @@ REAL(EB) FUNCTION SOLID_PHASE_OUTPUT(INDX,Y_INDEX,Z_INDEX,PART_INDEX,OPT_WALL_IN
89828982
SOLID_PHASE_OUTPUT = 0.5_EB*( ONE_D%X(I_DEPTH-1) + ONE_D%X(I_DEPTH) )
89838983
ENDIF
89848984

8985-
CASE(48) ! PYROLYSIS DEPTH
8985+
CASE(48) ! LAYER DIVIDE DEPTH
89868986
IF (SF%THERMAL_BC_INDEX==THERMALLY_THICK) THEN
8987-
SOLID_PHASE_OUTPUT = ONE_D%PYROLYSIS_DEPTH
8987+
SOLID_PHASE_OUTPUT = ONE_D%LAYER_DIVIDE_DEPTH
89888988
ELSE
89898989
SOLID_PHASE_OUTPUT = 0._EB
89908990
ENDIF
@@ -9853,7 +9853,7 @@ SUBROUTINE DUMP_PROF(T,NM)
98539853
IF (NWP==0) CYCLE PROF_LOOP
98549854
CALL GET_WALL_NODE_WEIGHTS(NWP,ONE_D%N_LAYERS,ONE_D%N_LAYER_CELLS,ONE_D%LAYER_THICKNESS,SF%GEOMETRY, &
98559855
ONE_D%X(0:NWP),SF%LAYER_DIVIDE,DX_S(1:NWP),RDX_S(0:NWP+1),RDXN_S(0:NWP),DX_WGT_S(0:NWP),DXF,DXB,LAYER_INDEX,MF_FRAC,&
9856-
SF%INNER_RADIUS,ONE_D%PYROLYSIS_DEPTH)
9856+
SF%INNER_RADIUS,ONE_D%LAYER_DIVIDE_DEPTH)
98579857
ELSE
98589858
NWP = SF%N_CELLS_INI
98599859
IF (NWP==0) CYCLE PROF_LOOP

‎Source/func.f90‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -3979,7 +3979,7 @@ SUBROUTINE PACK_BOUNDARY_ONE_D(NM,IC,RC,LC,OS,OD_INDEX,UNPACK_IT,COUNT_ONLY)
39793979
IC=IC+1 ; IF (.NOT.COUNT_ONLY) CALL EQUATE(OS%INTEGERS(IC) , ONE_D%MATL_INDEX(NN) , UNPACK_IT)
39803980
ENDDO
39813981

3982-
RC=RC+1 ; IF (.NOT.COUNT_ONLY) CALL EQUATE(OS%REALS(RC) ,ONE_D%PYROLYSIS_DEPTH , UNPACK_IT)
3982+
RC=RC+1 ; IF (.NOT.COUNT_ONLY) CALL EQUATE(OS%REALS(RC) ,ONE_D%LAYER_DIVIDE_DEPTH , UNPACK_IT)
39833983
I1=RC+1 ; RC=I1+ONE_D%N_MATL-1 ; IF (.NOT.COUNT_ONLY) CALL EQUATE(OS%REALS(I1:RC),ONE_D%M_DOT_S_PP(1:RC-I1+1) , UNPACK_IT)
39843984
I1=RC+1 ; RC=I1+ONE_D%N_CELLS_MAX ; IF (.NOT.COUNT_ONLY) CALL EQUATE(OS%REALS(I1:RC),ONE_D%X(0:RC-I1) , UNPACK_IT)
39853985
I1=RC+1 ; RC=I1+ONE_D%N_CELLS_OLD-1 ; IF (.NOT.COUNT_ONLY) CALL EQUATE(OS%REALS(I1:RC),ONE_D%DX_OLD(1:RC-I1+1) , UNPACK_IT)

‎Source/type.f90‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -228,7 +228,7 @@ MODULE TYPES
228228
INTEGER, ALLOCATABLE, DIMENSION(:) :: RAMP_IHS_INDEX !< (1:N_LAYERS) RAMP index for HEAT_SOURCE
229229
INTEGER, ALLOCATABLE, DIMENSION(:) :: MATL_INDEX !< (1:N_MATL) Number of materials
230230

231-
REAL(EB) :: PYROLYSIS_DEPTH=0._EB !< Distance from surface that divides the direction of evolved pyrolyzates (m)
231+
REAL(EB) :: LAYER_DIVIDE_DEPTH=0._EB !< Distance from surface that divides the direction of evolved pyrolyzates (m)
232232

233233
REAL(EB), ALLOCATABLE, DIMENSION(:) :: M_DOT_S_PP !< (1:N_MATL) Mass production rate of solid species
234234
REAL(EB), ALLOCATABLE, DIMENSION(:) :: X !< (0:NWP) Depth (m), \f$ x_{{\rm s},i} \f$

‎Source/wall.f90‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -2049,7 +2049,7 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
20492049
NWP = SUM(ONE_D%N_LAYER_CELLS(1:ONE_D%N_LAYERS))
20502050
CALL GET_WALL_NODE_WEIGHTS(NWP,ONE_D%N_LAYERS,ONE_D%N_LAYER_CELLS(1:ONE_D%N_LAYERS),ONE_D%LAYER_THICKNESS,SF%GEOMETRY, &
20512051
ONE_D%X(0:NWP),LAYER_DIVIDE,DX_S(1:NWP),RDX_S(0:NWP+1),RDXN_S(0:NWP),DX_WGT_S(0:NWP),DXF,DXB,&
2052-
LAYER_INDEX(0:NWP+1),MF_FRAC(1:NWP),SF%INNER_RADIUS,ONE_D%PYROLYSIS_DEPTH)
2052+
LAYER_INDEX(0:NWP+1),MF_FRAC(1:NWP),SF%INNER_RADIUS,ONE_D%LAYER_DIVIDE_DEPTH)
20532053
ELSE COMPUTE_GRID
20542054
NWP = SF%N_CELLS_INI
20552055
DXF = SF%DXF
@@ -2476,7 +2476,7 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
24762476
CALL GET_WALL_NODE_WEIGHTS(NWP,ONE_D%N_LAYERS,ONE_D%N_LAYER_CELLS(1:ONE_D%N_LAYERS), &
24772477
ONE_D%LAYER_THICKNESS(1:ONE_D%N_LAYERS),SF%GEOMETRY,X_S_NEW(0:NWP),LAYER_DIVIDE,DX_S(1:NWP), &
24782478
RDX_S(0:NWP+1),RDXN_S(0:NWP),DX_WGT_S(0:NWP),DXF,DXB,LAYER_INDEX(0:NWP+1),MF_FRAC(1:NWP),&
2479-
SF%INNER_RADIUS,ONE_D%PYROLYSIS_DEPTH)
2479+
SF%INNER_RADIUS,ONE_D%LAYER_DIVIDE_DEPTH)
24802480
ENDIF
24812481

24822482
N_LAYER_CELLS_NEW(1:ONE_D%N_LAYERS) = ONE_D%N_LAYER_CELLS(1:ONE_D%N_LAYERS)
@@ -2625,7 +2625,7 @@ SUBROUTINE SOLID_HEAT_TRANSFER(NM,T,DT_BC,PARTICLE_INDEX,WALL_INDEX,CFACE_INDEX,
26252625
ONE_D%LAYER_THICKNESS(1:ONE_D%N_LAYERS))
26262626
CALL GET_WALL_NODE_WEIGHTS(NWP_NEW,ONE_D%N_LAYERS,N_LAYER_CELLS_NEW(1:ONE_D%N_LAYERS),ONE_D%LAYER_THICKNESS,SF%GEOMETRY, &
26272627
X_S_NEW(0:NWP_NEW),LAYER_DIVIDE,DX_S(1:NWP_NEW),RDX_S(0:NWP_NEW+1),RDXN_S(0:NWP_NEW),&
2628-
DX_WGT_S(0:NWP_NEW),DXF,DXB,LAYER_INDEX(0:NWP_NEW+1),MF_FRAC(1:NWP_NEW),SF%INNER_RADIUS,ONE_D%PYROLYSIS_DEPTH)
2628+
DX_WGT_S(0:NWP_NEW),DXF,DXB,LAYER_INDEX(0:NWP_NEW+1),MF_FRAC(1:NWP_NEW),SF%INNER_RADIUS,ONE_D%LAYER_DIVIDE_DEPTH)
26292629
IF(NWP_NEW < NWP) ONE_D%DX_OLD(NWP_NEW:NWP) = 0._EB ! Zero out old values if needed
26302630
ONE_D%DX_OLD(1:NWP_NEW) = DX_S(1:NWP_NEW)
26312631
ONE_D%LAYER_THICKNESS_OLD(1:ONE_D%N_LAYERS) = ONE_D%LAYER_THICKNESS(1:ONE_D%N_LAYERS)

0 commit comments

Comments
 (0)