This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I–ΔtL which maps the identity (no preconditioner) for Δt≪1 to Laplacian preconditioning for Δt≫1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for , away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

A cross-diffusion system modelling the information herding of individuals is analysed in a bounded domain with no-flux boundary conditions. The variables are the species' density and an influence function which modifies the information state of the individuals. The cross-diffusion term may stabilize or destabilize the system. Furthermore, it allows for a formal gradient-flow or entropy structure. Exploiting this structure, the global-in-time existence of weak solutions and the exponential decay to the constant steady state is proved in certain parameter regimes. This approach does not extend to all parameters. We investigate local bifurcations from homogeneous steady states analytically to determine whether this defines the validity boundary. This analysis shows that generically there is a gap in the parameter regime between the entropy approach validity and the first local bifurcation. Next, we use numerical continuation methods to track the bifurcating non-homogeneous steady states globally and to determine non-trivial stationary solutions related to herding behaviour. In summary, we find that the main boundaries in the parameter regime are given by the first local bifurcation point, the degeneracy of the diffusion matrix and a certain entropy decay validity condition. We study several parameter limits analytically as well as numerically, with a focus on the role of changing a linear damping parameter as well as a parameter controlling the cross-diffusion. We suggest that our paradigm of comparing bifurcation-generated obstructions to the parameter validity of global-functional methods could also be of relevance for many other models beyond the one studied here.

In this contribution, we introduce numerical continuation methods and bifurcation theory, techniques which find their roots in the study of dynamical systems, to the problem of tracing the parameter dependence of bound and resonant states of the quantum mechanical Schrödinger equation. We extend previous work on the subject [1] to systems of coupled equations.

Bound and resonant states of the Schrödinger equation can be determined through the poles of the S-matrix, a quantity that can be derived from the asymptotic form of the wave function. We introduce a regularization procedure that essentially transforms the S-matrix into its inverse and improves its smoothness properties, thus making it amenable to numerical continuation. This allows us to automate the process of tracking bound and resonant states when parameters in the Schrödinger equation are varied. We have applied this approach to a number of model problems with satisfying results.

We propose a fast local level set method for the inverse problem of gravimetry. The theoretical foundation for our approach is based on the following uniqueness result: if an open set D is star-shaped or x3-convex with respect to its center of gravity, then its exterior potential uniquely determines the open set D. To achieve this purpose constructively, the first challenge is how to parametrize this open set D as its boundary may have a variety of possible shapes. To describe those different shapes we propose to use a level-set function to parametrize the unknown boundary of this open set. The second challenge is how to deal with the issue of partial data as gravimetric measurements are only made on a part of a given reference domain Ω. To overcome this difficulty we propose a linear numerical continuation approach based on the single layer representation to find potentials on the boundary of some artificial domain containing the unknown set D. The third challenge is how to speed up the level set inversion process. Based on some features of the underlying inverse gravimetry problem such as the potential density being constant inside the unknown domain, we propose a novel numerical approach which is able to take advantage of these features so that the computational speed is accelerated by an order of magnitude. We carry out numerical experiments for both two- and three-dimensional cases to demonstrate the effectiveness of the new algorithm.