1. A comparative numerical study of finite element methods resulting in mass conservation for Poisson's problem: Primal hybrid, mixed and their hybridized formulations.
- Author
-
Oliari, Victor B., Hancco Ancori, Ricardo J., and Devloo, Philippe R.B.
- Subjects
- *
FINITE element method , *CONSERVATION of mass , *APPROXIMATION error , *LINEAR systems , *CONDENSATION - Abstract
This paper presents a numerical comparison of finite-element methods resulting in local mass conservation at the element level for Poisson's problem, namely the primal hybrid and mixed methods. These formulations result in an indefinite system. Alternative formulations yielding a positive-definite system are obtained after hybridizing each method. The choice of approximation spaces yields methods with enhanced accuracy for the pressure variable, and results in systems with identical size and structure after static condensation. A regular pressure precision mixed formulation is also considered based on the classical RT k space. The simulations are accelerated using Open multi-processing (OMP) and Thread-Building Blocks (TBB) multithreading paradigms alongside either a coloring strategy or atomic operations ensuring a thread-safe execution. An additional parallel strategy is developed using C++ threads, which is based on the producer-consumer paradigm, and uses locks and semaphores as synchronization primitives. Numerical tests show the optimal parallel strategy for these finite-element formulations, and the computational performance of the methods are compared in terms of simulation time and approximation errors. Additional results are developed during the process. Numerical solvers often fail to find an accurate solution to the highly indefinite systems arising from finite-element formulations, and this paper documents a matrix ordering strategy to stabilize the resolution. A procedure to enable static condensation based on the introduction of piecewise constant functions that also fulfills Neumann's compatibility condition, and yet computes an average pressure per element is presented. • Numerical comparison of the mixed, primal hybrid, and other locally conservatives finite element methods. • Introduction of a twice hybridized primal hybrid method, yielding a positive definite matrix after static condensation. • Comparison of multi-threaded strategies and thread-safety mechanisms. • A numerical procedure to stabilize the linear system resolution for indefinite matrices is described. • A general procedure that enables static condensation and fulfills Neumann's compatibility condition. [ABSTRACT FROM AUTHOR]
- Published
- 2024
- Full Text
- View/download PDF