Natural convection from the vertical plate surface
Vertical steel plate divides two air environments with different temperatures.
Plane problem of steady-state heat transfer.
Height of the plate l = 1000 mm, thickness of the plate d = 2 mm.
Plate surface area is 1 m2.
temperature of the air below the plate T1 = 20°C,
temperature of the air above the plate T2 = -10°C
thermal conductivity of the plate λplate = 40 W/(m·K),
thermal conductivity of the ambient air λair = 0.027 W/(m·K).
Calculate the temperatures of the right and the left surfaces of the plate and the heat flux passing through the surfaces.
The heat flux flows from hot to cold areas (in our case from the left surface of the plate then through the plate toward air on the right).
Heat flux in the plate is regulated by the plate thermal conductivity. The heat exchange between air and the plate is regulated by natural convection. Convection coefficients depend on many factors such as environment temperature, shape and orientation of the surfaces.
For the vertical plate an average convection coefficients could be calculated using the similarity theory of heat transfer. We estimate that the temperatures of the upper and the lower surfaces of the plate equal to the average temperature of the plate which is Tplate = (T1+T2)/2 = 5 °C for calculation of the convection coefficients. The model is the thermal conductivity problem in the plate with the boundary conditions of the natural convection on the lower and the upper surfaces. The procedure of the solution is:
- Prandtl number (Pr), Grashof number (Gr), Relay number (Ra) and Nusselt number (Nu) are calculated by the known formulae*:
Pr = μ·C / λair, where μ - dynamical viscosity of the air, C - thermal capacity of the air;
Gr = L3ρ2g·ΔT·β / μ2, where L - characteristic size (in our case it is the plate height l), ρ - air density, g - acceleration of gravity, ΔT - temperature difference, β - coefficient of the air thermal expansion;
Ra = Gr·Pr;
Nu = [0.825 + 0.387·Ra1/6/(1 + (0.492/Pr)9/16)8/27]2.
Parameters of the air (ρ, C, μ, β) are calculated for the average temperatures:
below plate: Tair1 = (T1 + Tplate)/2 = (20+5)/2 = 12.5 °C;
C1 = 1000 J/(kg·K), μ1 = 1.87·10-5 N·s/m2, ρ1 = 1.25 kg/m3, β1 = 0.003501 1/K;
above plate: Tair2 = (T2 + Tplate)/2 = (-10+5)/2 = -2.5 °C**.
C2 = 1000 J/(kg·K), μ2 = 1.87·10-5 N·s/m2, ρ2 = 1.32 kg/m3, β2 = 0.003695 1/K.
Numerical values of the parameters:
ΔT1 = 20 - 5 = 15 °C, ΔT2 = 5 - (-10) = 15 °C;
Pr1 = 0.69; Pr2=0.69;
Gr1 = 2.30·109, Gr2 = 2.71·109;
Ra1 = 1.59·108, Ra2 = 1.88·108;
Nu1 = 141, Nu2 = 149.
Calculations are automated using the spreadsheet prepared by Harlan Bengtson and Lamar Stonecypher. Spreadsheet for the vertical plane convection is attached.
- Using the numerical values, we have
α = Nu·λair /L
After the substitution of the values, numerical values could be obtained
α1 = 3.82 W/(m2·K) for the left surface,
α2 = 4.02 W/(m2·K) for the right surface.
- QuickField problem is solved and the resulting temperatures of the surfaces and the heat flux through the surface is calculated.
Temperature of the left surface of the plate: T1 = 277.77 K (4.62 °C).
Temperature of the right surface of the plate: T2 = 277.76 K (4.61 °C).
Heat flux: F = 58.756 W (plate surface area is 1 m2).
Download simulation files and parameter calculation spreadsheet.
*CIBSE Guide C: Reference Data (2010). Butterworth-Heinemann, ISBN: 0750653604. Table 3.5.
**Remark: Since the formulas of the similarity theory relate to the big plates, in our case formula above gives the average value of the Nusselt number for each of the surfaces. This approach is very approximate by its nature, so there is no reason to make iterations to calculate the accurate Nusselt number related to actual surfaces temperatures.