𝔖 Bobbio Scriptorium
✦   LIBER   ✦

Implicit and semi-implicit schemes: Algorithms

✍ Scribed by R. Keppens; G. Tóth; M.A. Botchev; A. Van Der Ploeg


Publisher
John Wiley and Sons
Year
1999
Tongue
English
Weight
214 KB
Volume
30
Category
Article
ISSN
0271-2091

No coin nor oath required. For personal study only.

✦ Synopsis


This study formulates general guidelines to extend an explicit code with a great variety of implicit and semi-implicit time integration schemes. The discussion is based on their specific implementation in the Versatile Advection Code, which is a general purpose software package for solving systems of non-linear hyperbolic (and/or parabolic) partial differential equations, using standard high resolution shock capturing schemes. For all combinations of explicit high resolution schemes with implicit and semi-implicit treatments, it is shown how second-order spatial and temporal accuracy for the smooth part of the solutions can be maintained. Strategies to obtain steady state and time accurate solutions implicitly are discussed. The implicit and semi-implicit schemes require the solution of large linear systems containing the Jacobian matrix. The Jacobian matrix itself is calculated numerically to ensure the generality of this implementation. Three options are discussed in terms of applicability, storage requirements and computational efficiency. One option is the easily implemented matrix-free approach, but the Jacobian matrix can also be calculated by using a general grid masking algorithm, or by an efficient implementation for a specific Lax-Friedrich-type total variation diminishing (TVD) spatial discretization. The choice of the linear solver depends on the dimensionality of the problem. In one dimension, a direct block tridiagonal solver can be applied, while in more than one spatial dimension, a conjugate gradient (CG)-type iterative solver is used. For advection-dominated problems, preconditioning is needed to accelerate the convergence of the iterative schemes. The modified block incomplete LU-preconditioner is implemented, which performs very well. Examples from two-dimensional hydrodynamic and magnetohydrodynamic computations are given. They model transonic stellar outflow and recover the complex magnetohydrodynamic bow shock flow in the switch-on regime found in De Sterck et al. [Phys. Plasmas, 5, 4015 (1998)].


📜 SIMILAR VOLUMES


Accuracy and stability of semi-implicit
✍ Gjesdal, Thor ;Teigland, Rune 📂 Article 📅 1998 🏛 John Wiley and Sons 🌐 English ⚖ 134 KB

This paper discusses implementation strategies for second-order ®nite dierence discretizations of advection. Purely explicit and implicit methods both have disadvantages, and we consider semi-implicit schemes in which the ¯ux is split into a primary implicit part and a secondary explicit correction.

Multidomain implicit numerical scheme
✍ A. Povitsky; M. Wolfshtein 📂 Article 📅 1997 🏛 John Wiley and Sons 🌐 English ⚖ 270 KB

A multidomain method for the solution of elliptic CFD problems with an ADI scheme is described. Two methods of treatment of internal boundary conditions for ADI functions are discussed, namely an explicit and a semiimplicit method. Stability conditions for the proposed methods are derived theoretica