The implementation of a broad class of boundary conditions for non-linear hyperbolic systems

by   Edward W. G. Skevington, et al.

We propose methods that augment existing numerical schemes for the simulation of hyperbolic balance laws with Dirichlet boundary conditions to allow for the simulation of a broad class of differential algebraic conditions. Our approach is similar to that of Thompson (1987), where the boundary values were simulated by combining characteristic equations with the time derivative of the algebraic conditions, but differs in two important regards. Firstly, when the boundary is a characteristic of one of the fields Thompson's method can fail to produce reasonable values. We propose a method of combining the characteristic equations with extrapolation which ensures convergence. Secondly, the application of algebraic conditions can suffer from O(1) drift-off error, and we discuss projective time-stepping algorithms designed to converge for this type of system. Test problems for the shallow water equations are presented to demonstrate the result of simulating with and without the modifications discussed, illustrating their necessity for certain problems.



There are no comments yet.


page 1

page 2

page 3

page 4


A class of boundary conditions for time-discrete Green-Naghdi equations with bathymetry

This work is devoted to the structure of the time-discrete Green-Naghdi ...

Energy versus entropy estimates for nonlinear hyperbolic systems of equations

We compare and contrast information provided by the energy analysis of K...

Sharp stability for finite difference approximations of hyperbolic equations with boundary conditions

In this article, we consider a class of finite rank perturbations of Toe...

Computational homogenization of time-harmonic Maxwell's equations

In this paper we consider a numerical homogenization technique for curl-...

Weak Boundary Condition Enforcement for Linear Kirchhoff-Love Shells: Formulation, Error Analysis, and Verification

Stable and accurate modeling of thin shells requires proper enforcement ...

A linear system for pipe flow stability analysis allowing for boundary condition modifications

An accurate system to study the stability of pipe flow that ensures regu...

Neutron Transmission Strain Tomography for Non-Constant Strain-Free Lattice Spacing

Recently, several algorithms for strain tomography from energy-resolved ...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.