Replies: 2 comments 2 replies
-
Create a steel box with 1 m sides. Add a hot inner box 0.8 m on a side centered inside the box, and then record &HEAD CHID='A', TITLE='Check heat transfer through interpolated boundaries' / &MESH IJK=21,20,20, XB=-0.50, 0.55,-0.5, 0.5,-0.5, 0.5 / &MISC GVEC=0.,0.,0., NOISE=.FALSE., STRATIFICATION=.FALSE. / &TIME T_END=3600., DT=0.5 / &SPEC ID='NITROGEN', BACKGROUND=.TRUE. / &MATL ID = 'STEEL' &SURF ID='HOT', TMP_FRONT=500., EMISSIVITY=1., COLOR='RED' / &SURF ID = 'STEEL PLATE' &OBST XB=-0.40, 0.40,-0.40, 0.40,-0.40, 0.40, SURF_ID='HOT' / &OBST XB=-0.55,-0.50,-0.50, 0.50,-0.50, 0.50 / &VENT MB='XMIN', SURF_ID='OPEN' / &BNDF QUANTITY='WALL TEMPERATURE', CELL_CENTERED=.TRUE. / &DEVC XYZ=0.49,0.00,0.00, QUANTITY='RADIATIVE HEAT FLUX', ID='Flux_in', IOR=-1 / &DEVC XYZ=-.56,0.00,0.00, QUANTITY='WALL TEMPERATURE', ID='Temp_s_o', IOR=-1 / &SLCF PBY=0., QUANTITY='INTEGRATED INTENSITY', CELL_CENTERED=.TRUE. / &DEVC ID='QC_O',XYZ=-.56,0.00,0.00,QUANTITY='CONVECTIVE HEAT FLUX',IOR=-1/ &DEVC ID='QC_I',XYZ=-.49,0.00,0.00,QUANTITY='CONVECTIVE HEAT FLUX',IOR=1/ &TAIL / |
Beta Was this translation helpful? Give feedback.
-
This is just too much. It is difficult to untangle all this. If you believe that the code has a bug, submit a simple case to the Issue Tracker. No one has the time to sort out all these cases and answer all your questions. |
Beta Was this translation helpful? Give feedback.
-
I am trying to investigate heat transfer, including convection and radiation, to a door near a heat source.
Therefore, I made some modifications to the sample file (back_wall_test.fds) to output the temperature of the door.
The modifications at this point are as follows.
•Initialized the EMISSIVITY of the steel (EMISSIVITY=1 -> EMISSIVITY=0.9).
•Removed ”HEAT_TRANSFER_COEFFICIENT=0” from STEEL PLATE and HOT.
•Temperature output devices were added to the wall surface.
•Temperature SLCFs were added at each axis origin.
Comparing the results of this modified sample file (A.fds) and the two cases (g981.fds,EMISSIVITY0.fds) created based on the modified sample file (A.fds), I found the following unknowns in each.
(1) Initialization of GVEC (changing g from 0 to 9.81). (file:g981.fds)
I want to consider combustion in the upper layer, so I want to make gravity work and create a high temperature layer. So I removed "GVEC=0,0,0" and initialized it.
From Equation 8.3 in the user guide, if g=9.81, h should be higher than when g=0. Therefore, I thought that more heat would be transferred to the outside if GVEC was initialized.
However, as shown in the figure below, the result is that less heat is transferred to the outside than when g=0.
Q1: Could you tell me what caused the loss of heat transfer when g was changed from 0 to 9.81?
(2) I set EMISSIVITY=0 and tried to see if only the convection term could be output. (file:EMISSIVITY0.fds)
"back_wall_test.fds" outputs only the radiative heat flux term. I tried to see if it is possible to output only the convective heat flux term.
From Equation 22.33 in the User's Guide, I thought that if EMISSIVITY = 0, then εs = 0 and only the radiative heat flux term could be set to 0. The result is shown in the graph below. Where QC_O is the convective heat flux incident on the heater-side wall, QC_I is the convective heat flux incident on the opposite wall, QR_O is the radiant heat flux incident on the heater-side wall, and QR_I is the radiant heat flux incident on the opposite wall.
In fact, when I set EMISSIVITY=0 and output the results, I confirmed that RADIATIVE_HEAT_FLUX is set to 0. On the other hand, the value of CONVECTIVE_HEAT_FLUX was different from that when EMISSIVITY=0.9.
Q2:Does the value of EMISSIVITY also affect the convective heat flux term?
Q3: Could you please tell me why the value of CONVECTIVE_HEAT_FLUX was different for EMISSIVITY=0 and EMISSIVITY=0.9?
I appreciate your clarification on the above three questions.
Beta Was this translation helpful? Give feedback.
All reactions