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

07/12/2020
by   Lambert Theisen, et al.
0

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.

READ FULL TEXT

page 1

page 17

page 18

page 20

research
08/10/2022

Mixed finite element methods for the ferrofluid model with magnetization paralleled to the magnetic field

Mixed finite element methods are considered for a ferrofluid flow model ...
research
02/01/2018

Slate: extending Firedrake's domain-specific abstraction to hybridized solvers for geoscience and beyond

Within the finite element community, discontinuous Galerkin (DG) and mix...
research
02/01/2018

A domain-specific language for the hybridization and static condensation of finite element methods

In this paper, we introduce a domain-specific language (DSL) for concise...
research
06/06/2022

A High Order Stabilized Solver for the Volume Averaged Navier-Stokes Equations

The Volume-Averaged Navier-Stokes equations are used to study fluid flow...
research
09/14/2022

An Automated Process for 2D and 3D Finite Element Overclosure and Gap Adjustment using Radial Basis Function Networks

In biomechanics, geometries representing complicated organic structures ...
research
11/23/2021

A Conservative Finite Element Solver for MHD Kinematics equations: Vector Potential method and Constraint Preconditioning

A new conservative finite element solver for the three-dimensional stead...

Please sign up or login with your details

Forgot password? Click here to reset