Search

Total variation diminishing multigrid methods have been developed for first order accurate discretizations of hyperbolic conservation laws. This technique is based on a so-called upwind biased residual interpolation and allows for algorithms devoid of spurious numerical oscillations in the transient phase. In this paper, we justify the introduction of such prolongation and restriction operators by rewriting the algorithm in a matrix-vector notation. This perspective sheds new light on multigrid procedures for hyperbolic problems and provides a direct extension for high order accurate difference approximations. The new multigrid procedure is presented, advantages and disadvantages are discussed and numerical experiments are performed. Funding agencies: Linkoping University; VINNOVA, the Swedish Governmental Agency for Innovation SystemsVinnova [2013-01209]

Multigrid schemes for high order finite difference methods on summation-by-parts form are studied by comparing the effect of different interpolation operators. By using the standard linear prolongation and restriction operators, the Galerkin condition leads to inaccurate coarse grid discretizations. In this paper, an alternative class of interpolation operators that bypass this issue and preserve the summation-by-parts property on each grid level is considered. Clear improvements of the convergence rate for relevant model problems are achieved. Funding agencies: VINNOVA, the Swedish Governmental Agency for Innovation Systems [2013-01209]

Diagonal norm finite difference based time integration methods in summation-by-parts form are investigated. The second, fourth, and sixth order accurate discretizations are proven to have eigenvalues with strictly positive real parts. This leads to provably invertible fully discrete approximations of initial boundary value problems. Our findings also allow us to conclude that the Runge--Kutta methods based on second, fourth, and sixth order summation-by-parts finite difference time discretizations automatically satisfy previously unreported stability properties. The procedure outlined in this article can be extended to even higher order summation-by-parts approximations with repeating stencil. Funding agencies: VINNOVA, the Swedish Governmental Agency for Innovation SystemsVinnova [2013-01209]

Many physical phenomena can be described mathematically by means of partial differential equations. These mathematical formulations are said to be well-posed if a unique solution, bounded by the gi ...

We present a modified formulation of the dual time-stepping technique which makes use of two derivatives in pseudo-time. This new technique retains and improves the convergence properties to the stationary solution. When compared with the conventional dual time-stepping, the method with two derivatives reduces the stiffness of the problem and requires fewer iterations for full convergence to steady-state. In the current formulation, these positive effects require that an approximation of the square root of the spatial operator is available and inexpensive. Funding agencies: Linkoping University; Swedish Governmental Agency for Innovation SystemsVinnova [2013-01209]; VINNOVAVinnova

Fully-implicit discrete formulations in summation-by-parts form for initial-boundary value problems must be invertible in order to provide well functioning procedures. We prove that, under mild assumptions, pseudo-spectral collocation methods for the time derivative lead to invertible discrete systems when energy-stable spatial discretizations are used.

We discuss conservative and stable numerical approximations in summation-by-parts form for linear hyperbolic problems with variable coefficients. An extended setting, where the boundary or interface may or may not be included in the grid, is considered. We prove that conservative and stable formulations for variable coefficient problems require a boundary and interface conforming grid and exact numerical mimicking of integration-by-parts. Finally, we comment on how the conclusions from the linear analysis carry over to the nonlinear setting. Funding agencies: VINNOVA [2013-01209]

Localizing a radiant source is a widespread problem to many scientific and technological research areas. E.g. localization based on range measurements stays at the core of technologies like radar, sonar and wireless sensors networks. In this manuscript we study in depth the model for source localization based on range measurements obtained from the source signal, from the point of view of algebraic geometry. In the case of three receivers, we find unexpected connections between this problem and the geometry of Kummer's and Cayley's surfaces. Our work gives new insights also on the localization based on range differences., 38 pages, 18 figures