• ### Improving the IMEX method with a residual balanced decomposition(1705.04870)

Feb. 8, 2019 math.NA, cs.NA
In numerical time-integration with implicit-explicit (IMEX) methods, a within-step adaptable decomposition called residual balanced decomposition is introduced. With this decomposition, the requirement of a small enough residual in the iterative solver can be removed, consequently, this allows to exchange stability for efficiency. This decomposition transfers any residual occurring in the implicit equation of the implicit-step into the explicit part of the decomposition. By balancing the residual, the accuracy of the local truncation error of the time-stepping method becomes independent from the accuracy by which the implicit equation is solved. In order to balance the residual, the original IMEX decomposition is adjusted after the iterative solver has been stopped. For this to work, the traditional IMEX time-stepping algorithm needs to be changed. We call this new method the shortcut-IMEX (SIMEX). SIMEX can gain computational efficiency by exploring the trade-off between the computational effort placed in solving the implicit equation and the size of the numerically stable time-step. Typically, increasing the number of solver iterations increases the largest stable step-size. Both multi-step and Runge-Kutta (RK) methods are suitable for use with SIMEX. Here, we show the efficiency of a SIMEX-RK method in overcoming parabolic stiffness by applying it to a nonlinear reaction-advection-diffusion equation. In order to define a stability region for SIMEX, a region in the complex plane is depicted by applying SIMEX to a suitable PDE model containing diffusion and dispersion. A myriad of stability regions can be reached by changing the RK tableau and the solver.
• ### Electromagnetic instability of the Thomson Problem(cond-mat/0506616)

The classical Thomson problem of $n$ charged particles confined to the surface of a sphere of radius $a$ is analyzed within the Darwin approximation of electrodynamics. For $n<n_c(a)$ the ground state corresponds to a hexagonal Wigner crystal with a number of topological defects. However, if $n>n_c(a)$ the Wigner lattice is unstable with respect to small perturbations and the ground state becomes spontaneously magnetized for finite $n$.