Output methods in FDS

1,076 views
Skip to first unread message

nick

unread,
Jun 9, 2009, 11:15:07 AM6/9/09
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

Kevin

unread,
Jun 9, 2009, 11:27:55 AM6/9/09
to FDS and Smokeview Discussions
You have spent alot of time on this, but it appears that you have not
read the Users Guide. SLCF is NOT appropriate for solid surfaces. It
is for the gas TEMPERATURE. The QUANTITY 'TEMPERATURE' is for the gas,
'WALL TEMPERATURE' is for the solid. There are going to be big
differences between solid and gas temperatures.
> ===========================================================================­==============
> &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'/
>
> ===========================================================================­================

nick

unread,
Jun 9, 2009, 11:49:50 AM6/9/09
to FDS and Smokeview Discussions
Thank you Kevin for your reply. I know I should not use slice for
solid surfaces but this is just for understanding what the program is
doing. But still, the difference in temperature between a device
inside the solid (INSIDE WALL TEMPERATURE) and BNDF is large. Can you
give me any hint how BNDF is computing the solid surface temperature?
Isn't the computation of solid surface temperature dependent on the
gas phase temperature of the cell adjacent to the solid? Thats why I
was expecting not a big difference between the gas temperature (from
device) and solid surface temperature (from BNDF).

Please do not take these questions negatively. These are only efforts
to understand FDS better. Thank you.


regards,
Nick

Kevin

unread,
Jun 9, 2009, 12:14:43 PM6/9/09
to FDS and Smokeview Discussions
This is the wrong way to analyze FDS. I suggest you read the FDS
Technical Reference Guide to better understand the mathematical model.
Then, set up simple test cases, like those in the Verification Suite,
to test the heat transfer routines. If you rely only on Smokeview
output, you might get confused by the interpretation of color or
QUANTITY name.
> > > Nick- Hide quoted text -
>
> - Show quoted text -
Reply all
Reply to author
Forward
0 new messages