摘要：This work presents an implementation of a stabilized finite element formulation for incompressible viscous flow coupled with the advective-diffusive temperature transport equation through the Boussinesq approximation approach. For solving of the incompressible Navier-Stokes equations we use the Streamline Upwind Petrov-Galerkin/Pressure Stabilized Petrov-Galerkin (SUPG/PSPG) formulation and for the advective transport equation the SUPG formulation is employed. A comparison of the computational performance between full Gaussian and reduced (centre of the element) integrations for the isoparametric 8-node hexahedron element is presented. A h-stabilization for both advective and viscous/diffusive terms is used to control the spurious hourglass modes introduced by under integrating the element. The implementation has been performed using the libMesh Finite Element Method (FEM) library (http://libmesh.sourceforge.net) which provides support for adaptive mesh refinement and coarsening (AMR/C) and parallel computations. A verification is made using the Kim-Moin problem and the computational performance is evaluated solving a natural convection problem in a parallel machine.
