Stiff directed lines in random media

We investigate the localization of stiff directed lines with bending energy by a short-range random potential. We apply perturbative arguments, Flory scaling arguments, a variational replica calculation, and functional renormalization to show that a stiff directed line in 1+d dimensions undergoes a localization transition with increasing disorder for d > 2=3. We demonstrate that this transition is accessible by numerical transfer matrix calculations in 1+1 dimensions and analyze the properties of the disorder dominated phase in detail. On the basis of the two-replica problem, we propose a relation between the localization of stiff directed lines in 1+d dimensions and of directed lines under tension in 1+3d dimensions, which is strongly supported by identical free energy distributions. This shows that pair interactions in the replicated Hamiltonian determine the nature of directed line localization transitions with consequences for the critical behavior of the Kardar-Parisi-Zhang (KPZ) equation. We support the proposed relation to directed lines via multifractal analysis revealing an analogous Anderson transition-like scenario and a matching correlation length exponent. Furthermore, we quantify how the persistence length of the stiff directed line is reduced by disorder.Comment: Accepted for publication by Physical Review

Very High Order \PNM Schemes on Unstructured Meshes for the Resistive Relativistic MHD Equations

In this paper we propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space dimensions. The nonlinear system under consideration is purely hyperbolic and contains a source term, the one for the evolution of the electric field, that becomes stiff for low values of the resistivity. For the spatial discretization we propose to use high order \PNM schemes as introduced in \cite{Dumbser2008} for hyperbolic conservation laws and a high order accurate unsplit time discretization is achieved using the element-local space-time discontinuous Galerkin approach proposed in \cite{DumbserEnauxToro} for one-dimensional balance laws with stiff source terms. The divergence free character of the magnetic field is accounted for through the divergence cleaning procedure of Dedner et al. \cite{Dedneretal}. To validate our high order method we first solve some numerical test cases for which exact analytical reference solutions are known and we also show numerical convergence studies in the stiff limit of the RRMHD equations using \PNM schemes from third to fifth order of accuracy in space and time. We also present some applications with shock waves such as a classical shock tube problem with different values for the conductivity as well as a relativistic MHD rotor problem and the relativistic equivalent of the Orszag-Tang vortex problem. We have verified that the proposed method can handle equally well the resistive regime and the stiff limit of ideal relativistic MHD. For these reasons it provides a powerful tool for relativistic astrophysical simulations involving the appearance of magnetic reconnection.Comment: 24 pages, 6 figures, submitted to JC

Numerical Methods for Singular Perturbation Problems

Consider the two-point boundary value problem for a stiff system of ordinary differential equations. An adaptive method to solve these problems even when turning points are present is discussed

Segregation by membrane rigidity in flowing binary suspensions of elastic capsules

Spatial segregation in the wall normal direction is investigated in suspensions containing a binary mixture of Neo-Hookean capsules subjected to pressure driven flow in a planar slit. The two components of the binary mixture have unequal membrane rigidities. The problem is studied numerically using an accelerated implementation of the boundary integral method. The effect of a variety of parameters was investigated, including the capillary number, rigidity ratio between the two species, volume fraction, confinement ratio, and the number fraction of the more floppy particle $X_f$ in the mixture. It was observed that in suspensions of pure species, the mean wall normal positions of the stiff and the floppy particles are comparable. In mixtures, however, the stiff particles were found to be increasingly displaced towards the walls with increasing $X_f$, while the floppy particles were found to increasingly accumulate near the centerline with decreasing $X_f$. The origins of this segregation is traced to the effect of the number fraction $X_f$ on the localization of the stiff and the floppy particles in the near wall region -- the probability of escape of a stiff particle from the near wall region to the interior is greatly reduced with increasing $X_f$, while the exact opposite trend is observed for a floppy particle with decreasing $X_f$. Simple model studies on heterogeneous pair collisions involving a stiff and a floppy particle mechanistically explain this observation. The key result in these studies is that the stiff particle experiences much larger cross-stream displacement in heterogeneous collisions than the floppy particle. A unified mechanism incorporating the wall-induced migration of deformable particles and the particle fluxes associated with heterogeneous and homogeneous pair collisions is presented.Comment: 19 Pages, 16 Figure

Extrapolation-based implicit-explicit general linear methods

For many systems of differential equations modeling problems in science and engineering, there are natural splittings of the right hand side into two parts, one non-stiff or mildly stiff, and the other one stiff. For such systems implicit-explicit (IMEX) integration combines an explicit scheme for the non-stiff part with an implicit scheme for the stiff part. In a recent series of papers two of the authors (Sandu and Zhang) have developed IMEX GLMs, a family of implicit-explicit schemes based on general linear methods. It has been shown that, due to their high stage order, IMEX GLMs require no additional coupling order conditions, and are not marred by order reduction. This work develops a new extrapolation-based approach to construct practical IMEX GLM pairs of high order. We look for methods with large absolute stability region, assuming that the implicit part of the method is A- or L-stable. We provide examples of IMEX GLMs with optimal stability properties. Their application to a two dimensional test problem confirms the theoretical findings