nick
unread,Jun 9, 2009, 11:15:07 AM6/9/09Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to FDS and Smokeview Discussions
Hello everyone,
I am quite curious about how the output methods in FDS (e.g., DEVC,
SLCF, and BNDF) works in FDS, especially how the data are
interpolated. So, I have come up with a very simple example, which I
have pasted below, consisting of 3x3 grid cells in X and Y direction
and 9 cells in Z direction. I have inserted a device in each of the
cell center to measure temperature and SLCF at Y=0, 0.1, 0.2 and 0.3
m. Now, I can have the cell center values as well as nodal values of
temperature for all the cells. And I am interested to know how the
SLCF produces the nodal value at any of the nodes from the surrounding
cell centered values. I have tried it on my own for time instant 90
seconds with fds2ascii_linux, taking 2 sec. time average and then
comparing it with cell center temperature value at the same time
instant but I could not figure out how slice file deduces the nodal
values from the cell center values. It would be great if an example
(at a single node) could be produced from this example. If needed, I
can provide all the result files from my simulations, so that you
don't have to run simulation again.
The second part of this study is to investigate the difference in
results produced by BNDF and SLCF & DEVC for solid surface
temperature. What I found is the solid temperature at the surface of
ceiling deduced by BNDF is far less than that determined by SLCF
(22.36 °C vs. 100 °C) while the temperature determined by a device
inside the solid yielded 20 °C. Can anyone explain why is there such a
discrepancy? I am also interested to know how BNDF determines the
solid surface temperatures. Any hint, link to any literature or
explanations will really help me understand FDS and will be highly
appreciated. Thank you.
here is the input file:
#######################################
Turbulence Investigation
Cs = 0.23, Domain size = 1.0 x 1.0 x 1.2 m, Grid size = 0.02,
radiation OFF
#######################################
&HEAD CHID='interpolation', TITLE='Turbulence model investigation,
with central blower, FDS version 5.2.0' /
######################################
Computational domain
######################################
&MESH IJK=3,3,9, XB=0.0,0.30,0.0,0.30,0.0,0.90 / 0.10 m (10 cm) sized
mesh in all directions; 10 cm floor and ceiling
&TIME T_END=100.0 / Total simulation time
&MISC SURF_DEFAULT='WALL',CSMAG =
0.23,RADIATION=.FALSE./,RESTART=.TRUE./ Default surface-Wall
&DUMP NFRAMES=400 / Formatted output
#####################################
Blower properties
#####################################
&SURF ID='SUPPLY', VEL=-0.80, COLOR='BLUE',TMP_FRONT=200 /
#####################################
Material and surface properties
#####################################
&MATL ID = 'CONC'
FYI = 'Properties completely fabricated'
SPECIFIC_HEAT = 1.5
CONDUCTIVITY = 1.5
DENSITY = 2400.0/ Concrete for walls and roof slab
MATL ID = 'PASSIVE'
FYI = 'Properties completely fabricated'
SPECIFIC_HEAT = 1.5
CONDUCTIVITY = 1.5
DENSITY = 2400.0/ Passive material for floor and burner
&SURF ID='WALL'
RGB=128,128,128
MATL_ID='CONC'
THICKNESS = 0.10/
&SURF ID='ADBATIC'
RGB=150,150,150
MATL_ID='CONC'
THICKNESS = 0.20
ADIABATIC=.TRUE./Adiabatic surface properties for passive material
&PART ID='smoke', MASSLESS=.TRUE., SAMPLING_FACTOR=1 /Particle
properties for smoke
########################################
Geometry of compartment
########################################
&OBST XB=0.00,0.30,0.00,0.30,0.70,0.90, SURF_ID='WALL' / Roof slab
&OBST XB=0.00,0.30,0.00,0.30,0.00,0.20, SURF_ID='ADBATIC' / Adiabatic
floor slab
&VENT XB=0.0,0.10,0.0,0.10,0.20,0.20, SURF_ID='SUPPLY' / Blower
geometry
###############################################
Boundary conditions at the domain extremities
###############################################
&VENT MB='ZMAX',SURF_ID='OPEN' /
&VENT MB='XMAX',SURF_ID='OPEN' /
&VENT MB='XMIN',SURF_ID='MIRROR' /
&VENT MB='YMAX',SURF_ID='OPEN' /
&VENT MB='YMIN',SURF_ID='MIRROR' /
####################################
Output quantities
####################################
&SLCF PBY=0.0, QUANTITY='TEMPERATURE',VECTOR=.TRUE./
&SLCF PBY=0.10, QUANTITY='TEMPERATURE',VECTOR=.TRUE./
&SLCF PBY=0.20, QUANTITY='TEMPERATURE',VECTOR=.TRUE./
&SLCF PBY=0.30, QUANTITY='TEMPERATURE',VECTOR=.TRUE./
&DEVC XYZ=0.05,0.05,0.45 ,QUANTITY='TEMPERATURE',ID='x1y1z1'/
&DEVC XYZ=0.15,0.05,0.45 ,QUANTITY='TEMPERATURE',ID='x2y1z1'/
&DEVC XYZ=0.25,0.05,0.45 ,QUANTITY='TEMPERATURE',ID='x3y1z1'/
&DEVC XYZ=0.05,0.05,0.55 ,QUANTITY='TEMPERATURE',ID='x1y1z2'/
&DEVC XYZ=0.15,0.05,0.55 ,QUANTITY='TEMPERATURE',ID='x2y1z2'/
&DEVC XYZ=0.25,0.05,0.55 ,QUANTITY='TEMPERATURE',ID='x3y1z2'/
&DEVC XYZ=0.05,0.05,0.65 ,QUANTITY='TEMPERATURE',ID='x1y1z3'/
&DEVC XYZ=0.15,0.05,0.65 ,QUANTITY='TEMPERATURE',ID='x2y1z3'/
&DEVC XYZ=0.25,0.05,0.65 ,QUANTITY='TEMPERATURE',ID='x3y1z3'/
&DEVC XYZ=0.05,0.05,0.75 ,QUANTITY='TEMPERATURE', ID='x1y1z4'/
&DEVC XYZ=0.15,0.05,0.75 ,QUANTITY='TEMPERATURE', ID='x2y1z4'/
&DEVC XYZ=0.25,0.05,0.75 ,QUANTITY='TEMPERATURE', ID='x3y1z4'/
&DEVC XYZ=0.05,0.05,0.85 ,QUANTITY='TEMPERATURE',ID='x1y1z5'/
&DEVC XYZ=0.15,0.05,0.85 ,QUANTITY='TEMPERATURE',ID='x2y1z5'/
&DEVC XYZ=0.25,0.05,0.85 ,QUANTITY='TEMPERATURE',ID='x3y1z5'/
=========================================================================================
&DEVC XYZ=0.05,0.15,0.45 ,QUANTITY='TEMPERATURE',ID='x1y2z1'/
&DEVC XYZ=0.15,0.15,0.45 ,QUANTITY='TEMPERATURE',ID='x2y2z1'/
&DEVC XYZ=0.25,0.15,0.45 ,QUANTITY='TEMPERATURE',ID='x3y2z1'/
&DEVC XYZ=0.05,0.15,0.55 ,QUANTITY='TEMPERATURE',ID='x1y2z2'/
&DEVC XYZ=0.15,0.15,0.55 ,QUANTITY='TEMPERATURE',ID='x2y2z2'/
&DEVC XYZ=0.25,0.15,0.55 ,QUANTITY='TEMPERATURE',ID='x3y2z3'/
&DEVC XYZ=0.05,0.15,0.65 ,QUANTITY='TEMPERATURE',ID='x1y1z3'/
&DEVC XYZ=0.15,0.15,0.65 ,QUANTITY='TEMPERATURE',ID='x2y1z3'/
&DEVC XYZ=0.25,0.15,0.65 ,QUANTITY='TEMPERATURE',ID='x3y1z3'/
&DEVC XYZ=0.05,0.15,0.75 ,QUANTITY='TEMPERATURE', ID='x1y2z4'/
&DEVC XYZ=0.15,0.15,0.75 ,QUANTITY='TEMPERATURE', ID='x2y2z4'/
&DEVC XYZ=0.25,0.15,0.75 ,QUANTITY='TEMPERATURE', ID='x3y2z4'/
&DEVC XYZ=0.05,0.15,0.85 ,QUANTITY='TEMPERATURE', ID='x1y2z5'/
&DEVC XYZ=0.15,0.15,0.85 ,QUANTITY='TEMPERATURE', ID='x2y2z5'/
&DEVC XYZ=0.25,0.15,0.85 ,QUANTITY='TEMPERATURE', ID='x3y2z5'/
==========================================================================================
&DEVC XYZ=0.05,0.25,0.45 ,QUANTITY='TEMPERATURE',ID='x1y3z1'/
&DEVC XYZ=0.15,0.25,0.45 ,QUANTITY='TEMPERATURE',ID='x2y3z1'/
&DEVC XYZ=0.25,0.25,0.45 ,QUANTITY='TEMPERATURE',ID='x3y3z1'/
&DEVC XYZ=0.05,0.25,0.55 ,QUANTITY='TEMPERATURE',ID='x1y3z2'/
&DEVC XYZ=0.15,0.25,0.55 ,QUANTITY='TEMPERATURE',ID='x2y3z2'/
&DEVC XYZ=0.25,0.25,0.55 ,QUANTITY='TEMPERATURE',ID='x3y3z3'/
&DEVC XYZ=0.05,0.25,0.65 ,QUANTITY='TEMPERATURE',ID='x1y1z3'/
&DEVC XYZ=0.15,0.25,0.65 ,QUANTITY='TEMPERATURE',ID='x2y1z3'/
&DEVC XYZ=0.25,0.25,0.65 ,QUANTITY='TEMPERATURE',ID='x3y1z3'/
&DEVC XYZ=0.05,0.25,0.75 ,QUANTITY='TEMPERATURE', ID='x1y3z4'/
&DEVC XYZ=0.15,0.25,0.75 ,QUANTITY='TEMPERATURE', ID='x2y3z4'/
&DEVC XYZ=0.25,0.25,0.75 ,QUANTITY='TEMPERATURE', ID='x3y3z4'/
&DEVC XYZ=0.05,0.25,0.85 ,QUANTITY='TEMPERATURE', ID='x1y3z5'/
&DEVC XYZ=0.15,0.25,0.85 ,QUANTITY='TEMPERATURE', ID='x2y3z5'/
&DEVC XYZ=0.25,0.25,0.85 ,QUANTITY='TEMPERATURE', ID='x3y3z5'/
===========================================================================================
&BNDF QUANTITY='HEAT TRANSFER COEFFICIENT'/
&BNDF QUANTITY='GAS TEMPERATURE'/
&BNDF QUANTITY='WALL TEMPERATURE'/
&TAIL/ End of file
Thank you and wish you all a nice day.
Best regards,
Nick