An asymptotic preserving scheme based on a new formulation for NLS in the semiclassical limit

Christophe Besse (1, 2) , Rémi Carles (3) , Florian Méhats (4, 5)

1 LPP - Laboratoire Paul Painlevé - UMR 8524
2 SIMPAF - SImulations and Modeling for PArticles and Fluids
3 I3M - Institut de Mathématiques et de Modélisation de Montpellier
4 IPSO - Invariant Preserving SOlvers
5 IRMAR - Institut de Recherche Mathématique de Rennes

We consider the semiclassical limit for the nonlinear Schrodinger equation. We introduce a phase/amplitude representation given by a system similar to the hydrodynamical formulation, whose novelty consists in including some asymptotically vanishing viscosity. We prove that the system is always locally well-posed in a class of Sobolev spaces, and globally well-posed for a fixed positive Planck constant in the one-dimensional case. We propose a second order numerical scheme which is asymptotic preserving. Before singularities appear in the limiting Euler equation, we recover the quadratic physical observables as well as the wave function with mesh size and time step independent of the Planck constant. This approach is also well suited to the linear Schrodinger equation.

Domaines: Analyse numérique [math.NA]; Equations aux dérivées partielles [math.AP]

Christophe Besse, Rémi Carles, Florian Méhats. An asymptotic preserving scheme based on a new formulation for NLS in the semiclassical limit. Multiscale Modeling and Simulation: A SIAM Interdisciplinary Journal, 2013, 11 (4), pp.1228-1260. ⟨10.1137/120899017⟩. ⟨hal-00752011⟩