Score: 0

Pressure-robust enriched Galerkin finite element methods for coupled Navier-Stokes and heat equations

Published: December 18, 2025 | arXiv ID: 2512.16716v1

By: Sanjeeb Poudel, Sanghyun Lee, Lin Mu

Potential Business Impact:

Simulates tricky fluid and heat flow in any shape.

Business Areas:
Advanced Materials Manufacturing, Science and Engineering

We propose a pressure-robust enriched Galerkin (EG) finite element method for the incompressible Navier-Stokes and heat equations in the Boussinesq regime. For the Navier-Stokes equations, the EG formulation combines continuous Lagrange elements with a discontinuous enrichment vector per element in the velocity space and a piecewise constant pressure space, and it can be implemented efficiently within standard finite element frameworks. To enforce pressure robustness, we construct velocity reconstruction operators that map the discrete EG velocity field into exactly divergence-free, H(div)-conforming fields. In particular, we develop reconstructions based on Arbogast-Correa (AC) mixed finite element spaces on quadrilateral meshes and demonstrate that the resulting schemes remain stable and accurate even on highly distorted grids. The nonlinearity of the coupled Navier-Stokes-Boussinesq system is treated with several iterative strategies, including Picard iterations and Anderson-accelerated iterations; our numerical study shows that Anderson acceleration yields robust and efficient convergence for high Rayleigh number flows within the proposed framework. The performance of the method is assessed on a set of benchmark problems and application-driven test cases. These numerical experiments highlight the potential of pressure-robust EG methods as flexible and accurate tools for coupled flow and heat transport in complex geometries.

Country of Origin
🇺🇸 United States

Page Count
22 pages

Category
Computer Science:
Computational Engineering, Finance, and Science