We report on a high-fidelity, spectral/hp element algorithm developed for the direct numerical simulation of thermal convection problems. We consider the incompressible Navier-Stokes (NS) and advection-diffusion equations coupled through a thermal body-forcing term. The flow is driven by a prescribed flowrate forcing with explicit treatment of the nonlinear advection terms. The explicit treatment of the body-force term also decouples both the NS and the advection-diffusion equations. The problem is then temporally discretized using an implicit-explicit scheme in conjunction with a velocity-correction splitting scheme to decouple the velocity and pressure fields in the momentum equation. Although not unique, this type of discretization has not been widely applied to thermal convection problems and we therefore provide a comprehensive overview of the algorithm and implementation which is available through the open-source package Nektar++. After verifying the algorithm on a number of illustrative problems we then apply the code to investigate flow in a channel with uniform or streamwise sinusoidal lower wall, in addition to a patterned sinusoidal heating. We verify the solver against previously published two-dimensional results. Finally, for the first time we consider a three-dimensional problem with a streamwise sinusoidal lower wall and sinusoidal heating which, for the chosen parameter, leads to the unusual dynamics of an initially unsteady two-dimensional instability leading to a steady three-dimensional nonlinear saturated state.
K E Y W O R D Sincompressible Navier-Stokes, spectral/hp element method, thermal convection
INTRODUCTIONThermal convection intuitively plays an important role in many engineering processes and natural phenomena such as atmospheric circulation, daily weather events, and the generation of the Sun's magnetic field. 1 Techniques used for the thermal management of power plants, engines, turbo-machinery, and even computer design (that we use to simulate fluid flow problems) often employ the principle of thermal convection. Despite a plethora of applications, simulations of thermal convection flow problems using the spectral/hp element method are rather limited. This method combines the geometric flexibility of classical h-type (mesh-size h) finite element or finite volume methods with theThis is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.