Merge remote-tracking branch 'firemodels/master' into ht3d

mcgratta · Sep 11, 2023 · 2401d2b · 2401d2b
2 parents 0d8179f + f76ecc6
commit 2401d2b
Show file tree

Hide file tree

Showing 7 changed files with 157 additions and 108 deletions.
diff --git a/.github/workflows/linux.yml b/.github/workflows/linux.yml
@@ -31,7 +31,7 @@ jobs:
     runs-on: [ubuntu-latest]
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
 
       # install oneapi components from apt repository based on
       # oneapi-ci/scripts/setup_apt_repo_linux.sh
@@ -84,7 +84,7 @@ jobs:
     runs-on: [ubuntu-latest]
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - name: install openmpi
       run: |
         sudo apt-get update

diff --git a/.github/workflows/osx.yml b/.github/workflows/osx.yml
@@ -46,7 +46,7 @@ jobs:
         shell: bash
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
 
       # install oneapi components from web installer based on
       # oneapi-ci/scripts/install_macos.sh
@@ -141,7 +141,7 @@ jobs:
         shell: bash
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
     - name: install openmpi
       run: |
         brew install open-mpi

diff --git a/.github/workflows/windows.yml b/.github/workflows/windows.yml
@@ -44,7 +44,7 @@ jobs:
         shell: cmd
 
     steps:
-    - uses: actions/checkout@v3
+    - uses: actions/checkout@v4
 
       # install oneapi components from web installer based on
       # oneapi-ci/scripts/install_windows.bat

diff --git a/Manuals/FDS_User_Guide/FDS_User_Guide.tex b/Manuals/FDS_User_Guide/FDS_User_Guide.tex
@@ -9210,7 +9210,7 @@ \subsection{Spatially-Integrated Outputs}
 \begin{lstlisting}
 &DEVC XB=..., QUANTITY='TEMPERATURE', ID='maxT', SPATIAL_STATISTIC='MAX' /
 \end{lstlisting}
-causes FDS to write out the maximum gas phase temperature over the volume bounded by {\ct XB}. Other {\ct SPATIAL\_STATISTIC}'s are discussed below.  Some are appropriate for gas phase output quantities, some for solid phase, and some for both. Note that if {\ct XB} is used for a point device without a {\ct SPATIAL\_STATISTIC}, then FDS will define {\ct XYZ} to be the center of the volume defined by {\ct XB}.
+causes FDS to write out the maximum gas phase temperature over the volume bounded by {\ct XB}. Other {\ct SPATIAL\_STATISTIC}'s are discussed below.  Some are appropriate for gas phase output quantities, some for solid phase, and some for both. Note that if {\ct XB} is used for a point device without a {\ct SPATIAL\_STATISTIC}, then FDS will define {\ct XYZ} to be the center of the volume defined by {\ct XB}. A {\ct SPATIAL\_STATISTIC} can only be used for a gas phase {\ct QUANTITY} in Table~\ref{info:gasoutputquantities} with both D and S in the file type column or for a solid phase {\ct QUANTITY}  in Table~\ref{info:solidoutputquantities} with both D and B in the file type column.
 
 For solid phase output quantities, like heat fluxes and surface temperatures, the specification of a {\ct SURF\_ID} along with the appropriate statistic limits the calculation to only those surfaces. You can further limit the search by using the sextuplet of coordinates {\ct XB} to force FDS to only compute statistics for surface cells within the given volume. Be careful to account for the fact that the solid surface might shift to conform to the underlying numerical grid. Note that you do not (and should not) specify an orientation via the parameter {\ct IOR} when using a spatial statistic. {\ct IOR} is only needed to find a specific point on the solid surface.
 

diff --git a/Validation/NIST_Polymers/FDS_Input_Files/NIST_Gasification_Test_R3.fds b/Validation/NIST_Polymers/FDS_Input_Files/NIST_Gasification_Test_R3.fds
@@ -1,48 +1,65 @@
-&HEAD CHID='NIST_Gasification_Test_R3', TITLE='Black PMMA at 50 kW/m2 in Gasification Apparatus' /
+&HEAD CHID='NIST_Gasification_Test_R3', TITLE='Black PMMA at 50 kW/m2 in Gasification Apparatus Test R3' /
 
-&MESH IJK=3,3,4, XB=-0.15,0.15,-0.15,0.15,0.0,0.4 /
+&MESH XB=0,1,0,1,0,1,IJK=5,1,4/
 
-&TIME T_END=500., WALL_INCREMENT=1, DT=0.1 /
+&MISC SOLID_PHASE_ONLY=T,TMPA=20.45/
+&VENT XB=0.0,0.2,0,1,0,0,SURF_ID='SAMPLE1'/
 
-&MISC Y_O2_INFTY=0.001 /
+&TIME T_END=500,DT=0.01/
 
-&MATL ID                = 'Kaowool PM Insulation'
-      DENSITY           = 256.
-      SPECIFIC_HEAT     = 1.07
-      CONDUCTIVITY_RAMP = 'k_kaowool' /
+&CATF OTHER_FILES='pmma_properties.txt' /
 
-&RAMP ID='k_kaowool', T= 260., F=0.0576 /
-&RAMP ID='k_kaowool', T= 538., F=0.085  /
-&RAMP ID='k_kaowool', T= 816., F=0.125  /
-&RAMP ID='k_kaowool', T=1093., F=0.183  /
+&MATL ID='KPM',
+      DENSITY = 256.,
+      SPECIFIC_HEAT_RAMP = 'KPM_CP',
+      CONDUCTIVITY_RAMP = 'KPM_K',
+      EMISSIVITY = 0.92/
 
-&CATF OTHER_FILES='pmma_properties.txt' /
+&MATL ID='AL',
+      DENSITY = 2700.,
+      SPECIFIC_HEAT = 0.89,
+      CONDUCTIVITY = 100.,
+      EMISSIVITY = 0.3/
+
+! The PMMA was bonded to the the Kaowoll PM (KPM) using a combustible adhesive.
+! The mass of adhesive was treated as if it were PMMA and the sample thickess adjusted accordingly.
+&SURF ID='SAMPLE1',
+      HEAT_TRANSFER_COEFFICIENT=13.,
+      COLOR='WHITE',
+      MATL_ID='BLACK PMMA','KPM','AL',
+      THICKNESS=0.006091,0.02852,0.00635,
+      STRETCH_FACTOR=1,1,1,
+      EXTERNAL_FLUX=50/
 
-&SURF ID='PMMA SLAB'
-      COLOR='BLACK'
-      STRETCH_FACTOR=1.
-      HEAT_TRANSFER_COEFFICIENT=13.
-      EXTERNAL_FLUX=50 
-      BACKING='INSULATED'
-      MATL_ID='BLACK PMMA','Kaowool PM Insulation'
-      THICKNESS=0.0059,0.0254 /
+&RAMP ID='KPM_CP', T=26.85, F=0.7653/
+&RAMP ID='KPM_CP', T=226.85, F=1.0211/
+&RAMP ID='KPM_CP', T=526.85, F=1.2128/
+&RAMP ID='KPM_CP', T=626.85, F=1.1711/
+&RAMP ID='KPM_CP', T=1326.85, F=1.2834/
 
-&VENT XB=-0.05,0.05,-0.05,0.05,0.0,0.0, SURF_ID = 'PMMA SLAB' /
+&RAMP ID='KPM_K', T=260, F=0.0576/
+&RAMP ID='KPM_K', T=538, F=0.085/
+&RAMP ID='KPM_K', T=816, F=0.125/
+&RAMP ID='KPM_K', T=1093, F=0.183/
 
-&VENT MB='XMIN', SURF_ID='OPEN' /
-&VENT MB='XMAX', SURF_ID='OPEN' /
-&VENT MB='YMIN', SURF_ID='OPEN' /
-&VENT MB='YMAX', SURF_ID='OPEN' /
-&VENT MB='ZMAX', SURF_ID='OPEN' /
 
-&DUMP DT_DEVC=1. /
+&DEVC ID='Mass', QUANTITY='CONTROL VALUE',CTRL_ID='MP', UNITS='[g]' /
+&DEVC ID='MLR', XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='BURNING RATE', CONVERSION_FACTOR=1000, UNITS='[g/m2/s]' /
+&DEVC ID='Back Surface Temperature', XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='INSIDE WALL TEMPERATURE',DEPTH=-0.03487, CONVERSION_ADDEND=273, UNITS='[K]'/
+&DEVC ID='Top Surface Temperature',  XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='BLACK PMMA',ID='SD_B',IOR=3/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='PMMA_INT',ID='SD_I',IOR=3/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='PMMA_CHAR',ID='SD_C',IOR=3/
 
-&DEVC ID='Mass', XB=-0.05,0.05,-0.05,0.05,0.00,0.01, IOR=3, QUANTITY='SURFACE DENSITY', SPATIAL_STATISTIC='SURFACE INTEGRAL'
-      CONVERSION_FACTOR=1000, UNITS='[g]' /
-&DEVC ID='MLR', XB=-0.05,0.05,-0.05,0.05,0.00,0.01, IOR=3, QUANTITY='MASS FLUX', SPEC_ID='MMA', CONVERSION_FACTOR=1000, UNITS='[g/m2/s]' /
-&DEVC ID='Back Surface Temperature', XYZ=0.0,0.0,0.0, IOR=3, QUANTITY='BACK WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
-&DEVC ID='Top Surface Temperature',  XYZ=0.0,0.0,0.0, IOR=3, QUANTITY='WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
+! Sums the virgin, intermediate, and char masses
+&CTRL ID='SDsum',FUNCTION_TYPE='SUM',INPUT_ID='SD_B','SD_I','SD_C'/
 
-&TAIL / 
+! The functions below remove that portion of the PMMA mass that is the adhesive to give
+! SDs - PMMA surface density in kg/m2
+! SDn - Normalized surface density
+! MP - Mass of PMMA in g
+&CTRL ID='SDs',FUNCTION_TYPE='SUBTRACT',INPUT_ID='SDsum','CONSTANT',CONSTANT=0.40821/
+&CTRL ID='SDn',FUNCTION_TYPE='DIVIDE',INPUT_ID='SDs','CONSTANT',CONSTANT=7.3697/
+&CTRL ID='MP',FUNCTION_TYPE='MULTIPLY',INPUT_ID='SDn','CONSTANT',CONSTANT=28.254/
 
 
diff --git a/Validation/NIST_Polymers/FDS_Input_Files/NIST_Gasification_Test_R4.fds b/Validation/NIST_Polymers/FDS_Input_Files/NIST_Gasification_Test_R4.fds
@@ -1,48 +1,64 @@
-&HEAD CHID='NIST_Gasification_Test_R4', TITLE='Black PMMA at 50 kW/m2 in Gasification Apparatus' /
+&HEAD CHID='NIST_Gasification_Test_R4', TITLE='Black PMMA at 50 kW/m2 in Gasification Apparatus Test R4' /
 
-&MESH IJK=3,3,4, XB=-0.15,0.15,-0.15,0.15,0.0,0.4 /
+&MESH XB=0,1,0,1,0,1,IJK=5,1,4/
 
-&TIME T_END=500., WALL_INCREMENT=1, DT=0.1 /
+&MISC SOLID_PHASE_ONLY=T,TMPA=20.45/
+&VENT XB=0.0,0.2,0,1,0,0,SURF_ID='SAMPLE1'/
 
-&MISC Y_O2_INFTY=0.001 /
+&TIME T_END=500,DT=0.01/
 
-&MATL ID                = 'Kaowool PM Insulation'
-      DENSITY           = 256.
-      SPECIFIC_HEAT     = 1.07
-      CONDUCTIVITY_RAMP = 'k_kaowool' /
+&CATF OTHER_FILES='pmma_properties.txt' /
 
-&RAMP ID='k_kaowool', T= 260., F=0.0576 /
-&RAMP ID='k_kaowool', T= 538., F=0.085  /
-&RAMP ID='k_kaowool', T= 816., F=0.125  /
-&RAMP ID='k_kaowool', T=1093., F=0.183  /
+&MATL ID='KPM',
+      DENSITY = 256.,
+      SPECIFIC_HEAT_RAMP = 'KPM_CP',
+      CONDUCTIVITY_RAMP = 'KPM_K',
+      EMISSIVITY = 0.92/
 
-&CATF OTHER_FILES='pmma_properties.txt' /
+&MATL ID='AL',
+      DENSITY = 2700.,
+      SPECIFIC_HEAT = 0.89,
+      CONDUCTIVITY = 100.,
+      EMISSIVITY = 0.3/
 
-&SURF ID='PMMA SLAB'
-      COLOR='BLACK'
-      STRETCH_FACTOR=1.
-      HEAT_TRANSFER_COEFFICIENT=13.
-      EXTERNAL_FLUX=50 
-      BACKING='INSULATED'
-      MATL_ID='BLACK PMMA','Kaowool PM Insulation'
-      THICKNESS=0.00565,0.0254 /
+! The PMMA was bonded to the the Kaowoll PM (KPM) using a combustible adhesive.
+! The mass of adhesive was treated as if it were PMMA and the sample thickess adjusted accordingly.
+&SURF ID='SAMPLE1',
+      HEAT_TRANSFER_COEFFICIENT=13.,
+      COLOR='WHITE',
+      MATL_ID='BLACK PMMA','KPM','AL',
+      THICKNESS=0.006077,0.02852,0.00635,
+      STRETCH_FACTOR=1,1,1,
+      EXTERNAL_FLUX=50/
 
-&VENT XB=-0.05,0.05,-0.05,0.05,0.0,0.0, SURF_ID = 'PMMA SLAB' /
+&RAMP ID='KPM_CP', T=26.85, F=0.7653/
+&RAMP ID='KPM_CP', T=226.85, F=1.0211/
+&RAMP ID='KPM_CP', T=526.85, F=1.2128/
+&RAMP ID='KPM_CP', T=626.85, F=1.1711/
+&RAMP ID='KPM_CP', T=1326.85, F=1.2834/
 
-&VENT MB='XMIN', SURF_ID='OPEN' /
-&VENT MB='XMAX', SURF_ID='OPEN' /
-&VENT MB='YMIN', SURF_ID='OPEN' /
-&VENT MB='YMAX', SURF_ID='OPEN' /
-&VENT MB='ZMAX', SURF_ID='OPEN' /
+&RAMP ID='KPM_K', T=260, F=0.0576/
+&RAMP ID='KPM_K', T=538, F=0.085/
+&RAMP ID='KPM_K', T=816, F=0.125/
+&RAMP ID='KPM_K', T=1093, F=0.183/
 
-&DUMP DT_DEVC=1. /
+&DEVC ID='Mass', QUANTITY='CONTROL VALUE',CTRL_ID='MP', UNITS='[g]' /
+&DEVC ID='MLR', XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='BURNING RATE', CONVERSION_FACTOR=1000, UNITS='[g/m2/s]' /
+&DEVC ID='Back Surface Temperature', XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='INSIDE WALL TEMPERATURE',DEPTH=-0.03487, CONVERSION_ADDEND=273, UNITS='[K]'/
+&DEVC ID='Top Surface Temperature',  XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='BLACK PMMA',ID='SD_B',IOR=3/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='PMMA_INT',ID='SD_I',IOR=3/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='PMMA_CHAR',ID='SD_C',IOR=3/
 
-&DEVC ID='Mass', XB=-0.05,0.05,-0.05,0.05,0.00,0.01, IOR=3, QUANTITY='SURFACE DENSITY', SPATIAL_STATISTIC='SURFACE INTEGRAL'
-      CONVERSION_FACTOR=1000, UNITS='[g]' /
-&DEVC ID='MLR', XB=-0.05,0.05,-0.05,0.05,0.00,0.01, IOR=3, QUANTITY='MASS FLUX', SPEC_ID='MMA', CONVERSION_FACTOR=1000, UNITS='[g/m2/s]' /
-&DEVC ID='Back Surface Temperature', XYZ=0.0,0.0,0.0, IOR=3, QUANTITY='BACK WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
-&DEVC ID='Top Surface Temperature',  XYZ=0.0,0.0,0.0, IOR=3, QUANTITY='WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
+! Sums the virgin, intermediate, and char masses
+&CTRL ID='SDsum',FUNCTION_TYPE='SUM',INPUT_ID='SD_B','SD_I','SD_C'/
 
-&TAIL / 
+! The functions below remove that portion of the PMMA mass that is the adhesive to give
+! SDs - PMMA surface density in kg/m2
+! SDn - Normalized surface density
+! MP - Mass of PMMA in g
+&CTRL ID='SDs',FUNCTION_TYPE='SUBTRACT',INPUT_ID='SDsum','CONSTANT',CONSTANT=0.66487/
+&CTRL ID='SDn',FUNCTION_TYPE='DIVIDE',INPUT_ID='SDs','CONSTANT',CONSTANT=7.3535/
+&CTRL ID='MP',FUNCTION_TYPE='MULTIPLY',INPUT_ID='SDn','CONSTANT',CONSTANT=28.192/
 
 
diff --git a/Validation/NIST_Polymers/FDS_Input_Files/NIST_Gasification_Test_R5.fds b/Validation/NIST_Polymers/FDS_Input_Files/NIST_Gasification_Test_R5.fds
@@ -1,48 +1,64 @@
-&HEAD CHID='NIST_Gasification_Test_R5', TITLE='Black PMMA at 50 kW/m2 in Gasification Apparatus' /
+&HEAD CHID='NIST_Gasification_Test_R5', TITLE='Black PMMA at 50 kW/m2 in Gasification Apparatus Test R5' /
 
-&MESH IJK=3,3,4, XB=-0.15,0.15,-0.15,0.15,0.0,0.4 /
+&MESH XB=0,1,0,1,0,1,IJK=5,1,4/
 
-&TIME T_END=500., WALL_INCREMENT=1, DT=0.1 /
+&MISC SOLID_PHASE_ONLY=T,TMPA=20.45/
+&VENT XB=0.0,0.2,0,1,0,0,SURF_ID='SAMPLE1'/
 
-&MISC Y_O2_INFTY=0.001 /
+&TIME T_END=500,DT=0.01/
 
-&MATL ID                = 'Kaowool PM Insulation'
-      DENSITY           = 256.
-      SPECIFIC_HEAT     = 1.07
-      CONDUCTIVITY_RAMP = 'k_kaowool' /
+&CATF OTHER_FILES='pmma_properties.txt' /
 
-&RAMP ID='k_kaowool', T= 260., F=0.0576 /
-&RAMP ID='k_kaowool', T= 538., F=0.085  /
-&RAMP ID='k_kaowool', T= 816., F=0.125  /
-&RAMP ID='k_kaowool', T=1093., F=0.183  /
+&MATL ID='KPM',
+      DENSITY = 256.,
+      SPECIFIC_HEAT_RAMP = 'KPM_CP',
+      CONDUCTIVITY_RAMP = 'KPM_K',
+      EMISSIVITY = 0.92/
 
-&CATF OTHER_FILES='pmma_properties.txt' /
+&MATL ID='AL',
+      DENSITY = 2700.,
+      SPECIFIC_HEAT = 0.89,
+      CONDUCTIVITY = 100.,
+      EMISSIVITY = 0.3/
 
-&SURF ID='PMMA SLAB'
-      COLOR='BLACK'
-      STRETCH_FACTOR=1.
-      HEAT_TRANSFER_COEFFICIENT=13.
-      EXTERNAL_FLUX=50 
-      BACKING='INSULATED'
-      MATL_ID='BLACK PMMA','Kaowool PM Insulation'
-      THICKNESS=0.0062,0.0254 /
+! The PMMA was bonded to the the Kaowoll PM (KPM) using a combustible adhesive.
+! The mass of adhesive was treated as if it were PMMA and the sample thickess adjusted accordingly.
+&SURF ID='SAMPLE1',
+      HEAT_TRANSFER_COEFFICIENT=13.,
+      COLOR='WHITE',
+      MATL_ID='BLACK PMMA','KPM','AL',
+      THICKNESS=0.006523,0.02852,0.00635,
+      STRETCH_FACTOR=1,1,1,
+      EXTERNAL_FLUX=50/
 
-&VENT XB=-0.05,0.05,-0.05,0.05,0.0,0.0, SURF_ID = 'PMMA SLAB' /
+&RAMP ID='KPM_CP', T=26.85, F=0.7653/
+&RAMP ID='KPM_CP', T=226.85, F=1.0211/
+&RAMP ID='KPM_CP', T=526.85, F=1.2128/
+&RAMP ID='KPM_CP', T=626.85, F=1.1711/
+&RAMP ID='KPM_CP', T=1326.85, F=1.2834/
 
-&VENT MB='XMIN', SURF_ID='OPEN' /
-&VENT MB='XMAX', SURF_ID='OPEN' /
-&VENT MB='YMIN', SURF_ID='OPEN' /
-&VENT MB='YMAX', SURF_ID='OPEN' /
-&VENT MB='ZMAX', SURF_ID='OPEN' /
+&RAMP ID='KPM_K', T=260, F=0.0576/
+&RAMP ID='KPM_K', T=538, F=0.085/
+&RAMP ID='KPM_K', T=816, F=0.125/
+&RAMP ID='KPM_K', T=1093, F=0.183/
 
-&DUMP DT_DEVC=1. /
+&DEVC ID='Mass', QUANTITY='CONTROL VALUE',CTRL_ID='MP', UNITS='[g]' /
+&DEVC ID='MLR', XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='BURNING RATE', CONVERSION_FACTOR=1000, UNITS='[g/m2/s]' /
+&DEVC ID='Back Surface Temperature', XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='INSIDE WALL TEMPERATURE',DEPTH=-0.03487, CONVERSION_ADDEND=273, UNITS='[K]'/
+&DEVC ID='Top Surface Temperature',  XYZ=0.1,0.5,0.0, IOR=3, QUANTITY='WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='BLACK PMMA',ID='SD_B',IOR=3/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='PMMA_INT',ID='SD_I',IOR=3/
+&DEVC XYZ=0.1,0.5,0.0,QUANTITY='SURFACE DENSITY',MATL_ID='PMMA_CHAR',ID='SD_C',IOR=3/
 
-&DEVC ID='Mass', XB=-0.05,0.05,-0.05,0.05,0.00,0.01, IOR=3, QUANTITY='SURFACE DENSITY', SPATIAL_STATISTIC='SURFACE INTEGRAL'
-      CONVERSION_FACTOR=1000, UNITS='[g]' /
-&DEVC ID='MLR', XB=-0.05,0.05,-0.05,0.05,0.00,0.01, IOR=3, QUANTITY='MASS FLUX', SPEC_ID='MMA', CONVERSION_FACTOR=1000, UNITS='[g/m2/s]' /
-&DEVC ID='Back Surface Temperature', XYZ=0.0,0.0,0.0, IOR=3, QUANTITY='BACK WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
-&DEVC ID='Top Surface Temperature',  XYZ=0.0,0.0,0.0, IOR=3, QUANTITY='WALL TEMPERATURE', CONVERSION_ADDEND=273, UNITS='[K]'/
+! Sums the virgin, intermediate, and char masses
+&CTRL ID='SDsum',FUNCTION_TYPE='SUM',INPUT_ID='SD_B','SD_I','SD_C'/
 
-&TAIL / 
+! The functions below remove that portion of the PMMA mass that is the adhesive to give
+! SDs - PMMA surface density in kg/m2
+! SDn - Normalized surface density
+! MP - Mass of PMMA in g
+&CTRL ID='SDs',FUNCTION_TYPE='SUBTRACT',INPUT_ID='SDsum','CONSTANT',CONSTANT=0.56704/
+&CTRL ID='SDn',FUNCTION_TYPE='DIVIDE',INPUT_ID='SDs','CONSTANT',CONSTANT=7.8930/
+&CTRL ID='MP',FUNCTION_TYPE='MULTIPLY',INPUT_ID='SDn','CONSTANT',CONSTANT=30.289/