FDS6.9.1 and FDS6.10.0 #14408
Description
Hi, I'm getting different results using two versions of the 2D GEOM test. It's strange that it doesn't behave as well in FDS 6.10.0 as it did in the previous version, please help me out, thanks!
&HEAD CHID='10AULMAT'/
&TIME T_END=4.0, WALL_INCREMENT=1/
&DUMP COLUMN_DUMP_LIMIT=.TRUE., DT_RESTART=1.0, NFRAMES=4000, MASS_FILE=.TRUE., VELOCITY_ERROR_FILE=.TRUE./
&MISC GVEC=0.0,0.0,0.0, TMPA=25.0, CFL_MIN=0.6, CFL_MAX=0.8, SIMULATION_MODE='LES', VN_MIN=0.6, VN_MAX=0.8/
&PRES SOLVER=ULMAT /
&MESH ID='Mesh01', IJK=350,1,350, XB=-0.07,0.07,-2.0E-4,2.0E-4,-0.07,0.07/
&SPEC ID='NITROGEN', BACKGROUND=.TRUE./
&SPEC ID='OXYGEN', MASS_FRACTION_0=0.233/
&SPEC ID='CARBON DIOXIDE'/
&SPEC ID='ETHYLENE'/
&SPEC ID='WATER VAPOR'/
&REAC ID='Ethylene 1step',
FUEL='ETHYLENE',
A=2e12,
E=125520,
SPEC_ID_NU='ETHYLENE', 'OXYGEN', 'CARBON DIOXIDE', 'WATER VAPOR',
NU=-1,-3,2,2,
SPEC_ID_N_S='ETHYLENE' ,'OXYGEN',
N_S=0.1,1.65/
&DEVC ID='net heat flux[s]', QUANTITY='NET HEAT FLUX', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='rad heat flux[s]', QUANTITY='RADIATIVE HEAT FLUX', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='incident heat flux[s]', QUANTITY='INCIDENT HEAT FLUX', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='gas temp near the wall[s]', QUANTITY='GAS TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='wall temp[s]', QUANTITY='WALL TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='inside wall temp 0.00005[s]', DEPTH=5.0E-5, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='inside wall temp 0.0001[s]', DEPTH=1.0E-4, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='core temp 0.00015[s]', DEPTH=1.5E-4, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='core temp 0.0002[s]', DEPTH=2.0E-4, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='core temp 0.00025[s]', DEPTH=2.5E-4, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='HTC[s]', QUANTITY='HEAT TRANSFER COEFFICIENT', XYZ=0.0,0.0,4.0E-4, IOR=3/
&DEVC ID='gas temperature[g]', QUANTITY='TEMPERATURE', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='Gas density[g]', QUANTITY='DENSITY', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='Rad loss[g]', QUANTITY='RADIATION LOSS', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='rad heat flux gas[g]', QUANTITY='RADIATIVE HEAT FLUX GAS', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='VEL[g]', QUANTITY='VELOCITY', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='u-vel[g]', QUANTITY='U-VELOCITY', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='w-vel[g]', QUANTITY='W-VELOCITY', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='C2H4 mass frac[g]', QUANTITY='MASS FRACTION', SPEC_ID='ETHYLENE', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='CO2 mass frac[g]', QUANTITY='MASS FRACTION', SPEC_ID='CARBON DIOXIDE', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='H2O mass frac[g]', QUANTITY='MASS FRACTION', SPEC_ID='WATER VAPOR', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='O2 mass frac[g]', QUANTITY='MASS FRACTION', SPEC_ID='OXYGEN', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='N2 mass frac[g]', QUANTITY='MASS FRACTION', SPEC_ID='NITROGEN', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='C2H4 vol frac[g]', QUANTITY='VOLUME FRACTION', SPEC_ID='ETHYLENE', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='CO2 vol frac[g]', QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='H2O vol frac[g]', QUANTITY='VOLUME FRACTION', SPEC_ID='WATER VAPOR', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='O2 vol frac[g]', QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&DEVC ID='N2 vol frac[g]', QUANTITY='VOLUME FRACTION', SPEC_ID='NITROGEN', XYZ=0.0,0.0,4.1E-4, ORIENTATION=0.0,0.0,1.0/
&MATL ID='PE',
SPECIFIC_HEAT=2.3,
CONDUCTIVITY=0.34,
DENSITY=920.0,
EMISSIVITY=0.3,
N_REACTIONS=1,
HEAT_OF_REACTION=140.6,
MATL_ID(1,1)='moltenPE',
NU_MATL(1,1)=1.0,
N_S=0.0,
REFERENCE_TEMPERATURE=118.0,
HEATING_RATE=10.0,
PYROLYSIS_RANGE=5.0/
&MATL ID='moltenPE',
SPECIFIC_HEAT=2.5,
CONDUCTIVITY=0.25,
DENSITY=760.0,
EMISSIVITY=0.1,
N_REACTIONS=1,
HEAT_OF_REACTION=962.0,
SPEC_ID(1,1)='ETHYLENE',
NU_SPEC(1,1)=1.0,
A=5.0E+18,
E=3.0E+5/
&MATL ID='NiCr',
SPECIFIC_HEAT=0.48139,
CONDUCTIVITY=17.4,
DENSITY=8410.0,
EMISSIVITY=0.7/
&MATL ID='Stainless',
SPECIFIC_HEAT=0.5,
CONDUCTIVITY=16.2,
DENSITY=7850.0/
&SURF ID='Wire',
RGB=146,202,166,
MATL_ID(1,1)='PE',
MATL_ID(2,1)='NiCr',
MATL_MASS_FRACTION(1,1)=1.0,
MATL_MASS_FRACTION(2,1)=1.0,
THICKNESS(1:2)=1.5E-4,2.5E-4,
GEOMETRY='CYLINDRICAL',
INNER_RADIUS=0.0,
INTERNAL_HEAT_SOURCE=0.,2802548/
&GEOM ID='cylinder',
SURF_ID='Wire',
CYLINDER_LENGTH=0.0004,
CYLINDER_RADIUS=4E-4,
CYLINDER_ORIGIN=0.,0.,0.,
CYLINDER_AXIS=0.,1.,0.,
CYLINDER_NSEG_AXIS=1,
CYLINDER_NSEG_THETA=16/
&BNDF QUANTITY='ADIABATIC SURFACE TEMPERATURE'/
&BNDF QUANTITY='HEAT TRANSFER COEFFICIENT'/
&BNDF QUANTITY='INCIDENT HEAT FLUX'/
&BNDF QUANTITY='NET HEAT FLUX'/
&BNDF QUANTITY='RADIATIVE HEAT FLUX'/
&BNDF QUANTITY='WALL TEMPERATURE'/
&SLCF QUANTITY='TEMPERATURE', CELL_CENTERED=.TRUE., ID='Slice', PBY=0.0/
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='OXYGEN', CELL_CENTERED=.TRUE., ID='Slice01', PBY=0.0/
&SLCF QUANTITY='DENSITY', CELL_CENTERED=.TRUE., ID='Slice02', PBY=0.0/
&SLCF QUANTITY='HRRPUV', CELL_CENTERED=.TRUE., ID='Slice03', PBY=0.0/
&SLCF QUANTITY='PRESSURE', CELL_CENTERED=.TRUE., ID='Slice04', PBY=0.0/
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='ETHYLENE', CELL_CENTERED=.TRUE., ID='Slice05', PBY=0.0/
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', CELL_CENTERED=.TRUE., ID='Slice06', PBY=0.0/
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='WATER VAPOR', CELL_CENTERED=.TRUE., ID='Slice07', PBY=0.0/
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='NITROGEN', CELL_CENTERED=.TRUE., ID='Slice08', PBY=0.0/
&SLCF QUANTITY='VELOCITY', VECTOR=.TRUE., ID='Slice09', PBY=0.0/
&TAIL /
The gas distribution and velocity vector lines around the cylinder are more symmetric in FDS 6.9.1, but very asymmetric in FDS 6.10.1.