Permalink
Browse files

FDS Verification: Add additional test cases for MATL HoC

  • Loading branch information...
dr_jfloyd
dr_jfloyd committed Oct 10, 2018
1 parent be61004 commit 40ef1948972c83dbdbbb384c1c32d5c5c6449e95
Showing with 142 additions and 0 deletions.
  1. +2 −0 Verification/FDS_Cases.sh
  2. +63 −0 Verification/Fires/box_burn_away5.fds
  3. +77 −0 Verification/Fires/box_burn_away6.fds
@@ -133,6 +133,8 @@ $QFDS -d Fires box_burn_away1.fds
$QFDS -d Fires box_burn_away2.fds
$QFDS -d Fires box_burn_away3.fds
$QFDS -d Fires box_burn_away4.fds
$QFDS -d Fires box_burn_away5.fds
$QFDS -d Fires box_burn_away6.fds
$QFDS -d Fires box_burn_away_2D.fds
$QFDS -d Fires box_burn_away_2D_residue.fds
$QFDS -d Fires box_burn_away1_pyro3d_vs_pyro1d.fds
@@ -0,0 +1,63 @@
&HEAD CHID='box_burn_away5', TITLE='Test BURN_AWAY feature' /
The FOAM box is evaporated away by the high thermal radiation
from HOT surfaces. The mass of the box is 0.4^3 m3 * 20 kg/m3 = 1.28 kg.
This should be compared to the final value of fuel density volume integral,
computed by the first DEVC.
The gas species is mixture fraction fuel.
&MESH IJK=10,10,10 XB=0.0,1.0,0.0,1.0,0.0,1.0 /
&MISC Y_O2_INFTY=1.E-20/
&TIME T_END=30. DT = 0.01/
&MATL ID = 'FOAM'
HEAT_OF_REACTION = 800.
CONDUCTIVITY = 0.2
SPECIFIC_HEAT = 1.0
DENSITY = 20.
N_REACTIONS = 1
NU_SPEC = 1.
SPEC_ID = 'METHANE'
REFERENCE_TEMPERATURE= 200. /
&SURF ID = 'FOAM SLAB'
COLOR = 'TOMATO 3'
MATL_ID = 'FOAM'
THICKNESS = 0.1
BURN_AWAY = .TRUE.
BACKING = 'EXPOSED' /
&COMB N_SIMPLE_CHEMISTRY_REACTIONS=2/
&REAC FUEL='METHANE', AUTO_IGNITION_TEMPERATURE=15000. /
&DUMP SMOKE3D=.FALSE., MASS_FILE=.TRUE. /
&OBST XB=0.30,0.70,0.30,0.70,0.30,0.70, SURF_ID='FOAM SLAB', BULK_DENSITY=20. /
&SURF ID = 'HOT' TMP_FRONT = 1100. /
&VENT MB='XMIN', SURF_ID='HOT' /
&VENT MB='XMAX', SURF_ID='HOT' /
&VENT MB='YMIN', SURF_ID='HOT' /
&VENT MB='YMAX', SURF_ID='HOT' /
&VENT MB='ZMIN', SURF_ID='HOT' /
&VENT MB='ZMAX', SURF_ID='HOT' /
&BNDF QUANTITY='WALL TEMPERATURE' /
&BNDF QUANTITY='BURNING RATE' /
&BNDF QUANTITY='NET HEAT FLUX' /
&SLCF PBY=0.5, QUANTITY='TEMPERATURE' /
&SLCF PBY=0.5, QUANTITY='HRRPUV' /
&SLCF PBY=0.5, QUANTITY='INTEGRATED INTENSITY' /
&SLCF PBY=0.5, QUANTITY='MASS FRACTION',SPEC_ID='METHANE' /
&DEVC XB = 0,1,0,1,0,1, QUANTITY = 'DENSITY',SPEC_ID='METHANE',
STATISTICS = 'VOLUME INTEGRAL' ID = 'Mass fuel'/
&DEVC XB = 0,1,0,1,0,1, QUANTITY = 'HRR' ID = 'HRR' /
&TAIL /
@@ -0,0 +1,77 @@
&HEAD CHID='box_burn_away6', TITLE='Test BURN_AWAY feature' /
The FOAM box is evaporated away by the high thermal radiation
from HOT surfaces. The mass of the box is 0.4^3 m3 * 20 kg/m3 = 1.28 kg.
This should be compared to the final value of fuel density volume integral,
computed by the first DEVC.
The gas species is mixture fraction fuel.
&MESH IJK=10,10,10 XB=0.0,1.0,0.0,1.0,0.0,1.0 /
&TIME T_END=30. DT = 0.01/
&MATL ID = 'FOAM'
HEAT_OF_REACTION = 800.
CONDUCTIVITY = 0.2
SPECIFIC_HEAT = 1.0
DENSITY = 20.
N_REACTIONS = 1
NU_SPEC = 0.5,0.5
HEAT_OF_COMBUSTION = 50000,50000,
SPEC_ID = 'METHANE','PROPANE'
REFERENCE_TEMPERATURE= 200. /
&SURF ID = 'FOAM SLAB'
COLOR = 'TOMATO 3'
MATL_ID = 'FOAM'
THICKNESS = 0.1
BURN_AWAY = .TRUE.
BACKING = 'EXPOSED' /
&SPEC ID='NITROGEN',BACKGROUND=.TRUE./
&SPEC ID='OXYGEN'/
&SPEC ID='METHANE'/
&SPEC ID='PROPANE'/
&SPEC ID='CARBON DIOXIDE'/
&SPEC ID='WATER VAPOR'/
&REAC FUEL='METHANE', HEAT_OF_COMBUSTION=50000,
SPEC_ID_NU='METHANE','OXYGEN','CARBON DIOXIDE','WATER VAPOR',
NU=-1,-2,1,2/
&REAC FUEL='PROPANE', HEAT_OF_COMBUSTION=25000,
SPEC_ID_NU='PROPANE','OXYGEN','CARBON DIOXIDE','WATER VAPOR',
NU=-1,-5,3,4/
&DUMP SMOKE3D=.FALSE., MASS_FILE=.TRUE. /
&OBST XB=0.30,0.70,0.30,0.70,0.30,0.70, SURF_ID='FOAM SLAB', BULK_DENSITY=20. /
&SURF ID = 'HOT' TMP_FRONT = 1100. /
&VENT MB='XMIN', SURF_ID='HOT' /
&VENT MB='XMAX', SURF_ID='HOT' /
&VENT MB='YMIN', SURF_ID='HOT' /
&VENT MB='YMAX', SURF_ID='HOT' /
&VENT MB='ZMIN', SURF_ID='HOT' /
&VENT MB='ZMAX', SURF_ID='HOT' /
&BNDF QUANTITY='WALL TEMPERATURE' /
&BNDF QUANTITY='BURNING RATE' /
&BNDF QUANTITY='NET HEAT FLUX' /
&SLCF PBY=0.5, QUANTITY='TEMPERATURE' /
&SLCF PBY=0.5, QUANTITY='HRRPUV' /
&SLCF PBY=0.5, QUANTITY='INTEGRATED INTENSITY' /
&SLCF PBY=0.5, QUANTITY='MASS FRACTION',SPEC_ID='METHANE' /
&DEVC XB = 0,1,0,1,0,1, QUANTITY = 'DENSITY',SPEC_ID='METHANE',
STATISTICS = 'VOLUME INTEGRAL' ID = 'Mass fuel 1'/
&DEVC XB = 0,1,0,1,0,1, QUANTITY = 'DENSITY',SPEC_ID='PROPANE',
STATISTICS = 'VOLUME INTEGRAL' ID = 'Mass fuel 2'/
&DEVC XB = 0,1,0,1,0,1, QUANTITY = 'HRR' ID = 'HRR' /
&TAIL /

0 comments on commit 40ef194

Please sign in to comment.