fenicsR13: A Tensorial Mixed Finite Element Solver for the Linear R13 Equations Using the FEniCS Computing Platform

Abstract

We present a mixed finite element solver for the linearized R13 equations of non-equilibrium gas dynamics. The Python implementation builds upon the software tools provided by the FEniCS computing platform. We describe a new tensorial approach utilizing the extension capabilities of FEniCS's Unified Form Language (UFL) to define required differential operators for tensors above second degree. The presented solver serves as an example for implementing tensorial variational formulations in FEniCS, for which the documentation and literature seem to be very sparse. Using the software abstraction levels provided by the UFL allows an almost one-to-one correspondence between the underlying mathematics and the resulting source code. Test cases support the correctness of the proposed method using validation with exact solutions. To justify the usage of extended gas flow models, we discuss typical application cases involving rarefaction effects. We provide the documented and validated solver publicly.

arXiv comment29 pages, 13 figures, 8 listings, 3 table. Submitted to ACM Transactions on Mathematical Software (TOMS)

References in corpus (9)

• DOLFIN: Automated Finite Element Computing
• Firedrake: automating the finite element method by composing abstractions
• A compiler for variational forms
• Galerkin Least-Squares Stabilization in Ice Sheet Modeling - Accuracy, Robustness, and Comparison to other Techniques
• Regularized 13-Moment Equations for Inverse Power Law Models
• fenicsR13: A Tensorial Mixed Finite Element Solver for the Linear R13 Equations Using the FEniCS Computing Platform
• Moment Method for the Boltzmann Equation of Reactive Quaternary Gaseous Mixture
• Spline Moment Models for the one-dimensional Boltzmann-BGK equation
• Fast convergence and asymptotic preserving of the General Synthetic Iterative Scheme

Cited by in corpus (3)

• fenicsR13: A Tensorial Mixed Finite Element Solver for the Linear R13 Equations Using the FEniCS Computing Platform
• A generalized fundamental solution technique for the regularized 13-moment system in rarefied gas flows
• A Quasi-Optimal Factorization Preconditioner for Periodic Schrödinger Eigenstates in Anisotropically Expanding Domains