2022/07/26 by Ionut Danaila, Francky Luddens, Danaila, Ionut +3
Engineering · #Advanced Numerical Methods in Computational Mathematics #Computational Fluid Dynamics and Aerodynamics #Elasticity and Material Modeling #FOS: Mathematics #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.2207.12732
openalex publication_date 2022/07/26 · openalex created_date 2022/07/28 · openalex updated_date 2026/07/28
The purpose of this article is to study the convergence of a low order finite element approximation for a natural convection problem. We prove that the discretization based on P1 polynomials for every variable (velocity, pressure and temperature) is well-posed if used with a penalty term in the divergence equation, to compensate the loss of an inf-sup condition. With mild assumptions on the pressure regularity, we recover convergence for the Navier-Stokes-Boussinesq system, provided the penalty term is chosen in accordance with the mesh size. We express conditions to obtain optimal order of convergence. We illustrate theoretical convergence results with extensive examples. The computational cost that can be saved by this approach is also assessed.