An efficient implementation of mass conserving characteristic-based schemes in 2D and 3D

08/16/2019
by   Hanz Martin Cheng, et al.
0

In this paper, we develop the ball-approximated characteristics (B-char) method, which is an algorithm for efficiently implementing characteristic-based schemes in 2D and 3D. Core to the implementation of numerical schemes is the evaluation of integrals, which in the context of characteristic-based schemes with piecewise constant approximations boils down to computing the intersections between two regions. In the literature, these regions are approximated by polytopes (polygons in 2D and polyhedra in 3D) and, due to this, the implementation in 3D is nontrivial. The main novelty in this paper is the approximation of the regions by balls, whose intersections are much cheaper to compute than those of polytopes. Of course, balls cannot fully tessellate a region, and hence some mass may be lost. We perform some adjustments, and also solve an optimisation problem, in order to yield a scheme that is both locally and globally mass conserving. This algorithm can achieve results that are similar to those obtained from an implementation which uses polytopal intersections, with a much cheaper computational cost.

READ FULL TEXT

page 15

page 17

page 18

page 19

page 20

research
07/27/2022

A mass conservative Eulerian-Lagrangian Runge-Kutta discontinuous Galerkin method for wave equations with large time stepping

We propose an Eulerian-Lagrangian (EL) Runge-Kutta (RK) discontinuous Ga...
research
07/06/2022

Spline-oriented inter/extrapolation-based multirate schemes of higher order

Multirate integration uses different time step sizes for different compo...
research
10/29/2014

Efficient optimisation of structures using tabu search

This paper presents a novel approach to the optimisation of structures u...
research
09/24/2021

Approximations of Piecewise Deterministic Markov Processes and their convergence properties

Piecewise deterministic Markov processes (PDMPs) are a class of stochast...
research
05/09/2020

High order mixed finite elements with mass lumping for elasticity on triangular grids

A family of conforming mixed finite elements with mass lumping on triang...
research
04/19/2021

Efficient implementation of characteristic-based schemes on unstructured triangular grids

Using characteristics to treat advection terms in time-dependent PDEs le...
research
11/04/2022

Smallest covering regions and highest density regions for discrete distributions

This paper examines the problem of computing a canonical smallest coveri...

Please sign up or login with your details

Forgot password? Click here to reset