Ken Mattsson

Stable and Accurate Artificial Dissipation

Stability for nonlinear convection problems using centered difference schemes require the addition of artificial dissipation. In this paper we present dissipation operators that preserve both stability and accuracy for high order finite difference approximations of initial boundary value problems.

Erratum to: Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients

Finite difference operators approximating second derivatives with variable coefficients and satisfying a summation-by-parts rule have been derived for the second-, fourth- and sixth-order case by using the symbolic mathematics software Maple. The operators are based on the same norms as the corresponding approximations of the first derivative, which makes the construction of stable approximations to general multi-dimensional hyperbolic-parabolic problems straightforward.

Boundary Procedures for Summation-by-Parts Operators

Four different methods of imposing boundary conditions for the linear advection-diffusion equation and a linear hyperbolic system are considered. The methods are analyzed using the energy method and the Laplace transform technique. Numerical calculations are done, considering in particular the case when the initial data and boundary data are inconsistent.

Stable and Accurate Interpolation Operators for High-Order Multiblock Finite Difference Methods

Block-to-block interface interpolation operators are constructed for several common high-order finite difference discretizations. In contrast to conventional interpolation operators, these new interpolation operators maintain the strict stability, accuracy, and conservation properties of the base scheme even when nonconforming grids or dissimilar operators are used in adjoining blocks. The stability properties of the new operators are verified using eigenvalue analysis, and the accuracy properties are verified using numerical simulations of the Euler equations in two spatial dimensions.

Stable and accurate wave-propagation in discontinuous media

A time stable discretization is derived for the second-order wave equation with discontinuous coefficients. The discontinuity corresponds to inhomogeneity in the underlying medium and is treated by splitting the domain. Each (homogeneous) sub domain is discretized using narrow-diagonal summation by parts operators and, then, patched to its neighbors by using a penalty method, leading to fully explicit time integration. This discretization yields a time stable and efficient scheme. The analysis is verified by numerical simulations in one-dimension using high-order finite difference discretizations, and in three-dimensions using an unstructured finite volume discretization.

Steady-State Computations Using Summation-by-Parts Operators

This paper concerns energy stability on curvilinear grids and its impact on steady-state calulations. We have done computations for the Euler equations using fifth order summation-by-parts block and diagonal norm schemes. By imposing the boundary conditions weakly we obtain a fifth order energy-stable scheme. The calculations indicate the significance of energy stability in order to obtain convergence to steady state. Furthermore, the difference operators are improved such that faster convergence to steady state are obtained.

Stable Boundary Treatment for the Wave Equation on Second-Order Form

A stable and accurate boundary treatment is derived for the second-order wave equation. The domain is discretized using narrow-diagonal summation by parts operators and the boundary conditions are imposed using a penalty method, leading to fully explicit time integration. This discretization yields a stable and efficient scheme. The analysis is verified by numerical simulations in one-dimension using high-order finite difference discretizations, and in three-dimensions using an unstructured finite volume discretization.

High order finite difference methods for wave propagation in discontinuous media

High order finite difference approximations are derived for the second order wave equation with discontinuous coefficients, on rectangular geometries. The discontinuity is treated by splitting the domain at the discontinuities in a multi block fashion. Each sub-domain is discretized with compact second derivative summation by parts operators and the blocks are patched together to a global domain using the projection method. This guarantees a conservative, strictly stable and high order accurate scheme. The analysis is verified by numerical simulations in one and two spatial dimensions.

Boundary conditions for a divergence free velocity-pressure formulation of the Navier-Stokes equations

New sets of boundary conditions for the velocity-pressure formulation of the incompressible Navier-Stokes equations are derived. The boundary conditions have the same form on both inflow and outflow boundaries and lead to a divergence free solution. Moreover, the specific form of the boundary condition makes it possible derive a symmetric positive definite equation system for the internal pressure. Numerical experiments support the theoretical conclusions.

Optimal diagonal-norm SBP operators

Optimal boundary closures are derived for first derivative, finite difference operators of order 2, 4, 6 and 8. The closures are based on a diagonal-norm summation-by-parts (SBP) framework, thereby guaranteeing linear stability on piecewise curvilinear multi-block grids and entropy stability for nonlinear equations that support a convex extension. The new closures are developed by enriching conventional approaches with additional boundary closure stencils and non-equidistant grid distributions at the domain boundaries. Greatly improved accuracy is achieved near the boundaries, as compared with traditional diagonal-norm operators of the same order. The superior accuracy of the new optimal diagonal-norm SBP operators is demonstrated for linear hyperbolic systems in one dimension and for the nonlinear compressible Euler equations in two dimensions.