Search

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]

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