BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Erwan Faou

Erwan Faou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APmath-phmath.MPmath.PRmath.CAmath.DSNumerical Analysisphysics.comp-phphysics.plasm-ph

Catalog footprint

What is connected

25works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

On the influence of the geometry on skin effect in electromagnetism

We consider the equations of electromagnetism set on a domain made of a dielectric and a conductor subdomain in a regime where the conductivity is large. Assuming smoothness for the dielectric--conductor interface, relying on recent works we prove that the solution of the Maxwell equations admits a multiscale asymptotic expansion with profile terms rapidly decaying inside the conductor. This skin effect is measured by introducing a skin depth function that turns out to depend on the mean curvature of the boundary of the conductor. We then confirm these asymptotic results by numerical experiments in various axisymmetric configurations. We also investigate numerically the case of a nonsmooth interface, namely a cylindrical conductor.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linearized wave turbulence convergence results for three-wave systems

We consider stochastic and deterministic three-wave semi-linear systems with bounded and almost continuous set of frequencies. Such systems can be obtained by considering nonlinear lattice dynamics or truncated partial differential equations on large periodic domains. We assume that the nonlinearity is small and that the noise is small or void and acting only in the angles of the Fourier modes (random phase forcing). We consider random initial data and assume that these systems possess natural invariant distributions corresponding to some Rayleigh-Jeans stationary solutions of the wave kinetic equation appearing in wave turbulence theory. We consider random initial modes drawn with probability laws that are perturbations of theses invariant distributions. In the stochastic case, we prove that in the asymptotic limit (small nonlinearity, continuous set of frequency and small noise), the renormalized fluctuations of the amplitudes of the Fourier modes converge in a weak sense towards the solution of the linearized wave kinetic equation around these Rayleigh-Jeans spectra. Moreover, we show that in absence of noise, the deterministic equation with the same random initial condition satisfies a generic Birkhoff reduction in a probabilistic sense, without kinetic description at least in some regime of parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On weakly turbulent solutions to the perturbed linear Harmonic oscillator

We introduce specific solutions to the linear harmonic oscillator, named bubbles. They form resonant families of invariant tori of the linear dynamics, with arbitrarily large Sobolev norms. We use these modulated bubbles of energy to construct a class of potentials which are real, smooth, time dependent and uniformly decaying to zero with respect to time, such that the corresponding perturbed quantum harmonic oscillator admits solutions which exhibit a logarithmic growth of Sobolev norms. The resonance mechanism is explicit in space variables and produces highly oscillatory solutions. We then give several recipes to construct similar examples using more specific tools based on the continuous resonant (CR) equation in dimension two.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence of a Normalized Gradient Algorithm for Computing Ground States

We consider the approximation of the ground state of the one-dimensional cubic nonlinear Schr{ö}dinger equation by a normalized gradient algorithm combined with linearly implicit time integrator, and finite difference space approximation. We show that this method, also called imaginary time evolution method in the physics literature, is con-vergent, and we provide error estimates: the algorithm converges exponentially towards a modified solitons that is a space discretization of the exact soliton, with error estimates depending on the discretization parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collision of almost parallel vortex filaments

We investigate the occurrence of collisions in the evolution of vortex filaments through a system introduced by Klein, Majda and Damodaran [KMD95] and Zakharov [Z88,Z99]. We first establish rigorously the existence of a pair of almost parallel vortex filaments, with opposite circulation, colliding at some point in finite time. The collision mechanism is based on the one of the self-similar solutions of the model, described in [BFM14]. In the second part of this paper we extend this construction to the case of an arbitrary number of filaments, with polygonial symmetry, that are perturbations of a configuration of parallel vortex filaments forming a polygon, with or without its center, rotating with constant angular velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-order Hamiltonian splitting for Vlasov-Poisson equations

We consider the Vlasov-Poisson equation in a Hamiltonian framework and derive new time splitting methods based on the decomposition of the Hamiltonian functional between the kinetic and electric energy. Assuming smoothness of the solutions, we study the order conditions of such methods. It appears that these conditions are of Runge-Kutta-Nystr{ö}m type. In the one dimensional case, the order conditions can be further simplified, and efficient methods of order 6 with a reduced number of stages can be constructed. In the general case, high-order methods can also be constructed using explicit computations of commutators. Numerical results are performed and show the benefit of using high-order splitting schemes in that context. Complete and self-contained proofs of convergence results and rigorous error estimates are also given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On numerical Landau damping for splitting methods applied to the Vlasov-HMF model

We consider time discretizations of the Vlasov-HMF (Hamiltonian Mean-Field) equation based on splitting methods between the linear and non-linear parts. We consider solutions starting in a small Sobolev neighborhood of a spatially homogeneous state satisfying a linearized stability criterion (Penrose criterion). We prove that the numerical solutions exhibit a scattering behavior to a modified state, which implies a nonlinear Landau damping effect with polynomial rate of damping. Moreover, we prove that the modified state is close to the continuous one and provide error estimates with respect to the time stepsize.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Collisions of vortex filament pairs

We consider the problem of collisions of vortex filaments for a model introduced by Klein, Majda and Damodaran, and Zakharov to describe the interaction of almost parallel vortex filaments in three-dimensional fluids. Since the results of Crow examples of collisions are searched as perturbations of antiparallel translating pairs of filaments, with initial perturbations related to the unstable mode of the linearized problem; most results are numerical calculations. In this article we first consider a related model for the evolution of pairs of filaments and we display another type of initial perturbation leading to collision in finite time. Moreover we give numerical evidence that it also leads to collision through the initial model. We finally study the self-similar solutions of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulations of Kinetic Electrostatic Electron Nonlinear (KEEN) Waves with Variable Velocity Resolution Grids and High-Order Time-Splitting

KEEN waves are nonlinear, non-stationary, self-organized asymptotic states in Vlasov plasmas outside the scope or purview of linear theory constructs such as electron plasma waves or ion acoustic waves. Nonlinear stationary mode theories such as those leading to BGK modes also do not apply. The range in velocity that is strongly perturbed by KEEN waves depends on the amplitude and duration of the ponderomotive force used to drive them. Smaller amplitude drives create highly localized structures attempting to coalesce into KEEN waves. These cases have much more chaotic and intricate time histories than strongly driven ones. The narrow range in which one must maintain adequate velocity resolution in the weakly driven cases challenges xed grid numerical schemes. What is missing there is the capability of resolving locally in velocity while maintaining a coarse grid outside the highly perturbed region of phase space. We here report on a new Semi-Lagrangian Vlasov-Poisson solver based on conservative non-uniform cubic splines in velocity that tackles this problem head on. An additional feature of our approach is the use of a new high-order time-splitting scheme which allows much longer simulations per computational e ort. This is needed for low amplitude runs which take a long time to set up KEEN waves, if they are able to do so at all. The new code's performance is compared to uniform grid simulations and the advantages quanti ed. The birth pains associated with KEEN waves which are weakly driven is captured in these simulations. These techniques allow the e cient simulation of KEEN waves in multiple dimensions which will be tackled next as well as generalizations to Vlasov-Maxwell codes which are essential to understanding the impact of KEEN waves in practice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of the Monte-Carlo error in a hybrid semi-lagrangian scheme

We consider Monte-Carlo discretizations of partial differential equations based on a combination of semi-lagrangian schemes and probabilistic representations of the solutions. We study the Monte-Carlo error in a simple case, and show that under an anti-CFL condition on the time-step $δt$ and on the mesh size $δx$ and for $N$ - the number of realizations - reasonably large, we control this error by a term of order $\mathcal{O}(\sqrt{δt /N})$. We also provide some numerical experiments to confirm the error estimate, and to expose some examples of equations which can be treated by the numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast weak-KAM integrators for separable Hamiltonian systems

We consider a numerical scheme for Hamilton-Jacobi equations based on a direct discretization of the Lax-Oleinik semi-group. We prove that this method is convergent with respect to the time and space stepsizes provided the solution is Lipschitz, and give an error estimate. Moreover, we prove that the numerical scheme is a geometric integrator satisfying a discrete weak-KAM theorem which allows to control its long time behavior. Taking advantage of a fast algorithm for computing min-plus convolutions based on the decomposition of the function into concave and convex parts, we show that the numerical scheme can be implemented in a very efficient way.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Plane wave stability of the split-step Fourier method for the nonlinear Schrödinger equation

Plane wave solutions to the cubic nonlinear Schrödinger equation on a torus have recently been shown to behave orbitally stable. Under generic perturbations of the initial data that are small in a high-order Sobolev norm, plane waves are stable over long times that extend to arbitrary negative powers of the smallness parameter. The present paper studies the question as to whether numerical discretizations by the split-step Fourier method inherit such a generic long-time stability property. This can indeed be shown under a condition of linear stability and a non-resonance condition. They can both be verified if the time step-size is restricted by a CFL condition in the case of a constant plane wave. The proof first uses a Hamiltonian reduction and transformation and then modulated Fourier expansions in time. It provides detailed insight into the structure of the numerical solution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Resonant time steps and instabilities in the numerical integration of Schrödinger equations

We consider the linear and non linear cubic Schrödinger equations with periodic boundary conditions, and their approximations by splitting methods. We prove that for a dense set of arbitrary small time steps, there exists numerical solutions leading to strong numerical instabilities preventing the energy conservation and regularity bounds obtained for the exact solution. We analyze rigorously these instabilities in the semi-discrete and fully discrete cases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The weakly nonlinear large box limit of the 2D cubic nonlinear Schrödinger equation

We consider the cubic nonlinear Schrödinger (NLS) equation set on a two dimensional box of size $L$ with periodic boundary conditions. By taking the large box limit $L \to \infty$ in the weakly nonlinear regime (characterized by smallness in the critical space), we derive a new equation set on $\R^2$ that approximates the dynamics of the frequency modes. This nonlinear equation turns out to be Hamiltonian and enjoys interesting symmetries, such as its invariance under Fourier transform, as well as several families of explicit solutions. A large part of this work is devoted to a rigorous approximation result that allows to project the long-time dynamics of the limit equation into that of the cubic NLS equation on a box of finite size.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analysis of exponential splitting methods for inhomogeneous parabolic equations

We analyze the convergence of the exponential Lie and exponential Strang splitting applied to inhomogeneous second-order parabolic equations with Dirichlet boundary conditions. A recent result on the smoothing properties of these methods allows us to prove sharp convergence results in the case of homogeneous Dirichlet boundary conditions. When no source term is present and natural regularity assumptions are imposed on the initial value, we show full-order convergence of both methods. For inhomogeneous equations, we prove full-order convergence for the exponential Lie splitting, whereas order reduction to 1.25 for the exponential Strang splitting. Furthermore, we give sufficient conditions on the inhomogeneity for full-order convergence of both methods. Moreover our theoretical convergence results explain the severe order reduction to 0.25 of splitting methods applied to problems involving inhomogeneous Dirichlet boundary conditions. We include numerical experiments to underline the sharpness of our theoretical convergence results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic preserving schemes for the Klein-Gordon equation in the non-relativistic limit regime

We consider the Klein-Gordon equation in the non-relativistic limit regime, i.e. the speed of light c tending to infinity. We construct an asymptotic expansion for the solution with respect to the small parameter depending on the inverse of the square of the speed of light. As the first terms of this asymptotic can easily be simulated our approach allows us to construct numerical algorithms that are robust with respect to the large parameter c producing high oscillations in the exact solution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Existence and stability of solitons for fully discrete approximations of the nonlinear Schrödinger equation

In this paper we study the long time behavior of a discrete approximation in time and space of the cubic nonlinear Schrödinger equation on the real line. More precisely, we consider a symplectic time splitting integrator applied to a discrete nonlinear Schrödinger equation with additional Dirichlet boundary conditions on a large interval. We give conditions ensuring the existence of a numerical soliton which is close in energy norm to the continuous soliton. Such result is valid under a CFL condition between the time and space stepsizes. Furthermore we prove that if the initial datum is symmetric and close to the continuous soliton, then the associated numerical solution remains close to the orbit of the continuous soliton for very long times.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasi-periodic solutions of the 2D Euler equation

We consider the two-dimensional Euler equation with periodic boundary conditions. We construct time quasi-periodic solutions of this equation made of localized travelling profiles with compact support propagating over a stationary state depending on only one variable. The direction of propagation is orthogonal to this variable, and the support is concentrated on flat strips of the stationary state. The frequencies of the solution are given by the locally constant velocities associated with the stationary state.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sobolev stability of plane wave solutions to the cubic nonlinear Schrödinger equation on a torus

It is shown that plane wave solutions to the cubic nonlinear Schrödinger equation on a torus behave orbitally stable under generic perturbations of the initial data that are small in a high-order Sobolev norm, over long times that extend to arbitrary negative powers of the smallness parameter. The perturbation stays small in the same Sobolev norm over such long times. The proof uses a Hamiltonian reduction and transformation and, alternatively, Birkhoff normal forms or modulated Fourier expansions in time.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weak backward error analysis for SDEs

We consider numerical approximations of stochastic differential equations by the Euler method. In the case where the SDE is elliptic or hypoelliptic, we show a weak backward error analysis result in the sense that the generator associated with the numerical solution coincides with the solution of a modified Kolmogorov equation up to high order terms with respect to the stepsize. This implies that every invariant measure of the numerical scheme is close to a modified invariant measure obtained by asymptotic expansion. Moreover, we prove that, up to negligible terms, the dynamic associated with the Euler scheme is exponentially mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Nekhoroshev type theorem for the nonlinear Schrödinger equation on the d-dimensional torus.

We prove a Nekhoroshev type theorem for the nonlinear Schrödinger equation $$ iu_t=-Δu+V\star u+\partial_{\bar u}g(u,\bar u)\, \quad x\in \T^d, $$ where $V$ is a typical smooth potential and $g$ is analytic in both variables. More precisely we prove that if the initial datum is analytic in a strip of width $ρ>0$ with a bound on this strip equals to $\eps$ then, if $\eps$ is small enough, the solution of the nonlinear Schrödinger equation above remains analytic in a strip of width $ρ/2$ and bounded on this strip by $C\eps$ during very long time of order $ \eps^{-α|\ln \eps|^β}$ for some constants $C> 0$, $α>0$ and $β<1$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reconciling alternate methods for the determination of charge distributions: A probabilistic approach to high-dimensional least-squares approximations

We propose extensions and improvements of the statistical analysis of distributed multipoles (SADM) algorithm put forth by Chipot et al. in [6] for the derivation of distributed atomic multipoles from the quantum-mechanical electrostatic potential. The method is mathematically extended to general least-squares problems and provides an alternative approximation method in cases where the original least-squares problem is computationally not tractable, either because of its ill-posedness or its high-dimensionality. The solution is approximated employing a Monte Carlo method that takes the average of a random variable defined as the solutions of random small least-squares problems drawn as subsystems of the original problem. The conditions that ensure convergence and consistency of the method are discussed, along with an analysis of the computational cost in specific instances.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Quasi invariant modified Sobolev norms for semi linear reversible PDEs

We consider a general class of infinite dimensional reversible differential systems. Assuming a non resonance condition on the linear frequencies, we construct for such systems almost invariant pseudo norms that are closed to Sobolev-like norms. This allows us to prove that if the Sobolev norm of index $s$ of the initial data $z_0$ is sufficiently small (of order $ε$) then the Sobolev norm of the solution is bounded by $2ε$ during very long time (of order $ε^{-r}$ with $r$ arbitrary). It turns out that this theorem applies to a large class of reversible semi linear PDEs including the non linear Schrödinger equation on the d-dimensional torus. We also apply our method to a system of coupled NLS equations which is reversible but not Hamiltonian. We also notice that for the same class of reversible systems we can prove a Birkhoff normal form theorem that in turn implies the same bounds on the Sobolev norms. Nevertheless the technics that we use to prove the existence of quasi invariant pseudo norms is much more simple and direct.

0 citations0 reviewsNo code linkedOpen access