Open AccessAn improved method of Nusselt number calculation
Author(s) -
HoLiu Phyllis,
Hager Bradford H.,
Raefsky Arthur
Publication year - 1987
Publication title -
geophysical journal of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0016-8009
DOI - 10.1111/j.1365-246x.1987.tb01375.x
Subject(s) - nusselt number , dimensionless quantity , finite element method , heat flux , smoothing , mathematics , context (archaeology) , bilinear interpolation , mathematical analysis , mechanics , geometry , heat transfer , physics , thermodynamics , reynolds number , geology , paleontology , statistics , turbulence
Summary Accurate computation of heat flux is a primary motivation for convection calculations, yet heat fluxes have generally been calculated to substantially less than optimal accuracy, given the accompanying temperature and velocity fields. We analyse this problem in the context of the finite element method, using bilinear shape functions to interpolate the nodal values of velocity and temperature throughout the medium. Our method, however, is general and also applies to finite difference and higher order finite elements. Consider (dimensionless) vertical heat flux q z = wT – dT/dz, with w vertical velocity and T temperature. The first term is quadratic within (bilinear) elements and continuous at nodes, while the second term is constant within elements and discontinuous at nodes. Usual practice is to evaluate q z at element centres. We instead assign to each element the volume average of q z over that element. Element values are then projected to nodes using a modification of a technique used for pressure smoothing (Hughes et al. ) Our approach reduces the calculated depth variation in horizontally‐averaged flux by more than a factor of 10 and shows more rapid convergence of Nusselt number Nu as a function of grid size. It also eliminates the artificial diffusion of heat out of supposed insulating side walls.