Skip to main content

Viktor Linders: Iterative methods and structure preservation

Time: Thu 2023-03-30 14.00 - 15.00

Location: KTH, 3424, Lindstedsvägen 25

Participating: Viktor Linders (Lund)


Consistency, conservation, stability and accuracy are essential ingredients in the design of successful discretizations of partial differential equations. Considerable effort has been put into ensuring these properties for spatial and temporal discretizations. In the context of implicit schemes, iterative methods are used to approximate solutions to the systems of (nonlinear) equations that arise in each time step. However, iterative methods have not been designed with these concepts in mind.

In this talk we will discuss questions surrounding the preservation of structures by iterative methods and the potential impact that this could have on the quality of numerical solutions. We will see that not all iterative methods are conservative, and among those that are, many are inconsistent. In addition, iterative methods violate entropy estimates that the discretization satisfies, and can convert a dissipative scheme into an anti-dissipative one. This has a profound impact on the errors in numerical simulations. We will examine several solutions to these issues, and see that enforcing the preservation of the underlying structure can result in better numerical solutions, even when considerably larger tolerances are used.