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

Mickaël D. Chekroun

Mickaël D. Chekroun appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.DSmath.APmath-phmath.MPnlin.CDmath.PRmath.CAmath.NAmath.OCNumerical Analysisphysics.ao-phphysics.data-an

Catalog footprint

What is connected

15works
12topics
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

15 published item(s)

preprint2022arXiv

Deep spectral computations in linear and nonlinear diffusion problems

We propose a flexible machine-learning framework for solving eigenvalue problems of diffusion operators in moderately large dimension. We improve on existing Neural Networks (NNs) eigensolvers by demonstrating our approach ability to compute (i) eigensolutions for non-self adjoint operators with small diffusion (ii) eigenpairs located deep within the spectrum (iii) computing several eigenmodes at once (iv) handling nonlinear eigenvalue problems. To do so, we adopt a variational approach consisting of minimizing a natural cost functional involving Rayleigh quotients, by means of simple adiabatic technics and multivalued feedforward neural parametrisation of the solutions. Compelling successes are reported for a 10-dimensional eigenvalue problem corresponding to a Kolmogorov operator associated with a mixing Stepanov flow. We moreover show that the approach allows for providing accurate eigensolutions for a 5-D Schrödinger operator having $32$ metastable states. In addition, we address the so-called Gelfand superlinear problem having exponential nonlinearities, in dimension $4$, and for nontrivial domains exhibiting cavities. In particular, we obtain NN-approximations of high-energy solutions approaching singular ones. We stress that each of these results are obtained using small-size neural networks in situations where classical methods are hopeless due to the curse of dimensionality. This work brings new perspectives for the study of Ruelle-Pollicot resonances, dimension reduction, nonlinear eigenvalue problems, and the study of metastability when the dynamics has no potential.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Transitions in Stochastic Non-equilibrium Systems: Efficient Reduction and Analysis

A central challenge in physics is to describe non-equilibrium systems driven by randomness, such as a randomly growing interface, or fluids subject to random fluctuations that account e.g. for local stresses and heat fluxes not related to the velocity and temperature gradients. For deterministic systems with infinitely many degrees of freedom, normal form and center manifold theory have shown a prodigious efficiency to often completely characterize how the onset of linear instability translates into the emergence of nonlinear patterns. However, in presence of random fluctuations, this reduction procedure is seriously challenged due to large excursions caused by the noise, and the approach needs to be revisited. We present an alternative framework to cope with these difficulties exploiting the approximation theory of stochastic invariant manifolds and energy estimates measuring the defect of parameterization of the high-modes. To operate for fluid problems, these error estimates are derived under assumptions regarding dissipation effects brought by the high-modes that suitably counterbalance the loss of regularity due to the nonlinear terms. The approach enables us to predict, from the reduced equations, the occurrence in large probability of a stochastic analogue to the pitchfork bifurcation, as long as the noise's intensity and the eigenvalue's magnitude of the mildly unstable mode scale accordingly. Our parameterization formulas involve non-Markovian coefficients, which depend explicitly on the history of the noise path that drives the SPDE dynamics, and their memory content is self-consistently determined by the intensity of the random force and its interaction through the SPDE's nonlinear terms. Applications to a stochastic Rayleigh-Benard problem are detailed, for which conditions for a stochastic pitchfork bifurcation (in large probability) to occur, are clarified.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal management of harvested population at the edge of extinction

Optimal control of harvested population at the edge of extinction in an unprotected area, is considered. The underlying population dynamics is governed by a Kolmogorov-Petrovsky-Piskunov equation with a harvesting term and space-dependent coefficients while the control consists of transporting individuals from a natural reserve. The nonlinear optimal control problem is approximated by means of a Galerkin scheme. Convergence result about the optimal controlled solutions and error estimates between the corresponding optimal controls, are derived. For certain parameter regimes, nearly optimal solutions are calculated from a simple logistic ordinary differential equation (ODE) with a harvesting term, obtained as a Galerkin approximation of the original partial differential equation (PDE) model. A critical allowable fraction $\underlineα$ of the reserve's population is inferred from the reduced logistic ODE with a harvesting term. This estimate obtained from the reduced model allows us to distinguish sharply between survival and extinction for the full PDE itself, and thus to declare whether a control strategy leads to success or failure for the corresponding rescue operation while ensuring survival in the reserve's population. In dynamical terms, this result illustrates that although continuous dependence on the forcing may hold on finite-time intervals, a high sensitivity in the system's response may occur in the asymptotic time. We believe that this work, by its generality, establishes bridges interesting to explore between optimal control problems of ODEs with a harvesting term and their PDE counterpart.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ruelle-Pollicott Resonances of Stochastic Systems in Reduced State Space. Part II: Stochastic Hopf Bifurcation

The spectrum of the generator (Kolmogorov operator) of a diffusion process, referred to as the Ruelle-Pollicott (RP) spectrum, provides a detailed characterization of correlation functions and power spectra of stochastic systems via decomposition formulas in terms of RP resonances. Stochastic analysis techniques relying on the theory of Markov semigroups for the study of the RP spectrum and a rigorous reduction method is presented in Part I. This framework is here applied to study a stochastic Hopf bifurcation in view of characterizing the statistical properties of nonlinear oscillators perturbed by noise, depending on their stability. In light of the Hörmander theorem, it is first shown that the geometry of the unperturbed limit cycle, in particular its isochrons, is essential to understand the effect of noise and the phenomenon of phase diffusion. In addition, it is shown that the spectrum has a spectral gap, even at the bifurcation point, and that correlations decay exponentially fast. Explicit small-noise expansions of the RP eigenvalues and eigenfunctions are then obtained, away from the bifurcation point, based on the knowledge of the linearized deterministic dynamics and the characteristics of the noise. These formulas allow one to understand how the interaction of the noise with the deterministic dynamics affect the decay of correlations. Numerical results complement the study of the RP spectrum at the bifurcation, revealing useful scaling laws. The analysis of the Markov semigroup for stochastic bifurcations is thus promising in providing a complementary approach to the more geometric random dynamical system approach. This approach is not limited to low-dimensional systems and the reduction method presented in part I is applied to a stochastic model relevant to climate dynamics in part III.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Variational approach to closure of nonlinear dynamical systems: Autonomous case

A general, variational approach to derive low-order reduced systems for nonlinear systems subject to an autonomous forcing, is introduced. The approach is based on the concept of optimal parameterizing manifold (PM) that substitutes the more classical notion of slow manifold or invariant manifold when breakdown of slaving occurs. An optimal PM provides the manifold that describes the average motion of the neglected scales as a function of the resolved scales and allows, in principle, for determining the best vector field of the reduced state space that describes e.g. the dynamics' slow motion. The underlying optimal parameterizations are approximated by dynamically-based formulas derived analytically from the original equations. These formulas are contingent upon the determination of only a few (scalar) parameters obtained from minimization of cost functionals, depending on training dataset collected from direct numerical simulation. In practice, a training period of length comparable to a characteristic recurrence or decorrelation time of the dynamics, is sufficient for the efficient derivation of optimized parameterizations. Applications to the closure of low-order models of Atmospheric Primitive Equations and Rayleigh-Bénard convection are then discussed. The approach is finally illustrated --- in the context of the Kuramoto-Sivashinsky turbulence --- as providing efficient closures without slaving for a cutoff scale $k_c$ placed within the inertial range and the reduced state space is just spanned by the unstable modes, without inclusion of any stable modes whatsoever. The underlying optimal PMs obtained by our variational approach are far from slaving and allow for remedying the excessive backscatter transfer of energy to the low modes encountered by classical invariant manifold approximations in their standard forms when the latter are used at this cutoff wavelength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low-Dimensional Galerkin Approximations of Nonlinear Delay Differential Equations

This article revisits the approximation problem of systems of nonlinear delay differential equations (DDEs) by a set of ordinary differential equations (ODEs). We work in Hilbert spaces endowed with a natural inner product including a point mass, and introduce polynomials orthogonal with respect to such an inner product that live in the domain of the linear operator associated with the underlying DDE. These polynomials are then used to design a general Galerkin scheme for which we derive rigorous convergence results and show that it can be numerically implemented via simple analytic formulas. The scheme so obtained is applied to three nonlinear DDEs, two autonomous and one forced: (i) a simple DDE with distributed delays whose solutions recall Brownian motion; (ii) a DDE with a discrete delay that exhibits bimodal and chaotic dynamics; and (iii) a periodically forced DDE with two discrete delays arising in climate dynamics. In all three cases, the Galerkin scheme introduced in this article provides a good approximation by low-dimensional ODE systems of the DDE's strange attractor, as well as of the statistical features that characterize its nonlinear dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Data-driven non-Markovian closure models

This paper has two interrelated foci: (i) obtaining stable and efficient data-driven closure models by using a multivariate time series of partial observations from a large-dimensional system; and (ii) comparing these closure models with the optimal closures predicted by the Mori-Zwanzig (MZ) formalism of statistical physics. Multilayer stochastic models (MSMs) are introduced as both a generalization and a time-continuous limit of existing multilevel, regression-based approaches to closure in a data-driven setting; these approaches include empirical model reduction (EMR), as well as more recent multi-layer modeling. It is shown that the multilayer structure of MSMs can provide a natural Markov approximation to the generalized Langevin equation (GLE) of the MZ formalism. A simple correlation-based stopping criterion for an EMR-MSM model is derived to assess how well it approximates the GLE solution. Sufficient conditions are derived on the structure of the nonlinear cross-interactions between the constitutive layers of a given MSM to guarantee the existence of a global random attractor. This existence ensures that no blow-up can occur for a very broad class of MSM applications, a class that includes non-polynomial predictors and nonlinearities that do not necessarily preserve quadratic energy invariants. The approach is applied to conceptual nonlinear models borrowed from climate dynamics and population dynamics. In both cases, it is shown that the resulting closure models are able to capture the main statistical features of the dynamics, even in presence of weak time-scale separation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Finite-Horizon Parameterizing Manifolds, and Applications to Suboptimal Control of Nonlinear Parabolic PDEs

This article proposes a new approach based on finite-horizon parameterizing manifolds (PMs) for the design of low-dimensional suboptimal controllers to optimal control problems of nonlinear partial differential equations (PDEs) of parabolic type. Given a finite horizon $[0,T]$ and a low-mode truncation of the PDE, a PM provides an approximate parameterization of the uncontrolled high modes by the controlled low ones so that the unexplained high-mode energy is reduced, in an $L^2$-sense, when this parameterization is applied. Analytic formulas of such PMs are derived by application of the method of pullback approximation of the high-modes (Chekroun, Liu, and Wang, 2013, http://arxiv.org/pdf/1310.3896v1.pdf). These formulas allow for an effective derivation of reduced ODE systems, aimed to model the evolution of the low-mode truncation of the controlled state variable, where the high-mode part is approximated by the PM function applied to the low modes. A priori error estimates between the resulting PM-based low-dimensional suboptimal controller $u_R^\ast$ and the optimal controller $u^*$ are derived. These estimates demonstrate that the closeness of $u_R^\ast$ to $u^*$ is mainly conditioned on two factors: (i) the parameterization defect of a given PM, associated respectively with $u_R^\ast$ and $u^*$; and (ii) the energy kept in the high modes of the PDE solution either driven by $u_R^\ast$ or $u^*$ itself. The practical performances of such PM-based suboptimal controllers are numerically assessed for various optimal control problems associated with a Burgers-type equation. The numerical results show that a PM-based reduced system allows for the design of suboptimal controllers with good performances provided that the associated parameterization defects and energy kept in the high modes are small enough, in agreement with the rigorous results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical weather prediction in two dimensions with topography, using a finite volume method

We aim to study a finite volume scheme to solve the two dimensional inviscid primitive equations of the atmosphere with humidity and saturation, in presence of topography and subject to physically plausible boundary conditions to the system of equations. In that respect, a version of a projection method is introduced to enforce the compatibility condition on the horizontal velocity field, which comes from the boundary conditions. The resulting scheme allows for a significant reduction of the errors near the topography when compared to more standard finite volume schemes. In the numerical simulations, we first present the associated good convergence results that are satisfied by the solutions simulated by our scheme when compared to particular analytic solutions. We then report on numerical experiments using realistic parameters. Finally, the effects of a random small-scale forcing on the velocity equation is numerically investigated. The numerical results show that such a forcing is responsible for recurrent large-scale patterns to emerge in the temperature and velocity fields.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-Markovian Reduced Systems for Stochastic Partial Differential Equations: The Additive Noise Case

This article proposes for stochastic partial differential equations (SPDEs) driven by additive noise, a novel approach for the approximate parameterizations of the ``small'' scales by the ``large'' ones, along with the derivaton of the corresponding reduced systems. This is accomplished by seeking for stochastic parameterizing manifolds (PMs) introduced in a previous work by the authors(*), which are random manifolds aiming to provide --- in a mean square sense --- such approximate parameterizations. Backward-forward systems are designed to give access to such PMs as pullback limits depending through the nonlinear terms on the time-history of the dynamics of the low modes when the latter is simply approximated by its stochastic linear component. It is shown that the corresponding pullback limits can be efficiently determined, leading in turn to an operational procedure for the derivation of non-Markovian reduced systems able to achieve good modeling performances in practice. This is illustrated on a stochastic Burgers-type equation, where it is shown that the corresponding non-Markovian features of these reduced systems play a key role to reach such performances.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymptotics of the Coleman-Gurtin model

This paper is concerned with the integrodifferential equation $$\partial_t u-Δu -\int_0^\infty κ(s)Δu(t-s)\,\d s + φ(u)=f$$ arising in the Coleman-Gurtin's theory of heat conduction with hereditary memory, in presence of a nonlinearity $φ$ of critical growth. Rephrasing the equation within the history space framework, we prove the existence of global and exponential attractors of optimal regularity and finite fractal dimension for the related solution semigroup, acting both on the basic weak-energy space and on a more regular phase space.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Averaging of time-periodic systems without a small parameter

In this article, we present a new approach to averaging in non-Hamiltonian systems with periodic forcing. The results here do not depend on the existence of a small parameter. In fact, we show that our averaging method fits into an appropriate nonlinear equivalence problem, and that this problem can be solved formally by using the Lie transform framework to linearize it. According to this approach, we derive formal coordinate transformations associated with both first-order and higher-order averaging, which result in more manageable formulae than the classical ones. Using these transformations, it is possible to correct the solution of an averaged system by recovering the oscillatory components of the original non-averaged system. In this framework, the inverse transformations are also defined explicitly by formal series; they allow the estimation of appropriate initial data for each higher-order averaged system, respecting the equivalence relation. Finally, we show how these methods can be used for identifying and computing periodic solutions for a very large class of nonlinear systems with time-periodic forcing. We test the validity of our approach by analyzing both the first-order and the second-order averaged system for a problem in atmospheric chemistry.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Climate dynamics and fluid mechanics: Natural variability and related uncertainties

The purpose of this review-and-research paper is twofold: (i) to review the role played in climate dynamics by fluid-dynamical models; and (ii) to contribute to the understanding and reduction of the uncertainties in future climate-change projections. To illustrate the first point, we focus on the large-scale, wind-driven flow of the mid-latitude oceans which contribute in a crucial way to Earth's climate, and to changes therein. We study the low-frequency variability (LFV) of the wind-driven, double-gyre circulation in mid-latitude ocean basins, via the bifurcation sequence that leads from steady states through periodic solutions and on to the chaotic, irregular flows documented in the observations. This sequence involves local, pitchfork and Hopf bifurcations, as well as global, homoclinic ones. The natural climate variability induced by the LFV of the ocean circulation is but one of the causes of uncertainties in climate projections. Another major cause of such uncertainties could reside in the structural instability in the topological sense, of the equations governing climate dynamics, including but not restricted to those of atmospheric and ocean dynamics. We propose a novel approach to understand, and possibly reduce, these uncertainties, based on the concepts and methods of random dynamical systems theory. As a very first step, we study the effect of noise on the topological classes of the Arnol'd family of circle maps, a paradigmatic model of frequency locking as occurring in the nonlinear interactions between the El Nino-Southern Oscillations (ENSO) and the seasonal cycle. It is shown that the maps' fine-grained resonant landscape is smoothed by the noise, thus permitting their coarse-grained classification. This result is consistent with stabilizing effects of stochastic parametrization obtained in modeling of ENSO phenomenon via some general circulation models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Models of population dynamics under the influence of external perturbations: mathematical results

In this Note, we describe the stationary equilibria and the asymptotic behaviour of an heterogeneous logistic reaction-diffusion equation under the influence of autonomous or time-periodic forcing terms. We show that the study of the asymptotic behaviour in the time-periodic forcing case can be reduced to the autonomous one, the last one being described in function of the "size" of the external perturbation. Our results can be interpreted in terms of maximal sustainable yields from populations. We briefly discuss this last aspect through a numerical computation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On population resilience to external perturbations

We study a spatially explicit harvesting model in periodic or bounded environments. The model is governed by a parabolic equation with a spatially dependent nonlinearity of Kolmogorov--Petrovsky--Piskunov type, and a negative external forcing term $-δ$. Using sub- and supersolution methods and the characterization of the first eigenvalue of some linear elliptic operators, we obtain existence and nonexistence results as well as results on the number of stationary solutions. We also characterize the asymptotic behavior of the evolution equation as a function of the forcing term amplitude. In particular, we define two critical values $δ^*$ and $δ_2$ such that, if $δ$ is smaller than $δ^*$, the population density converges to a "significant" state, which is everywhere above a certain small threshold, whereas if $δ$ is larger than $δ_2$, the population density converges to a "remnant" state, everywhere below this small threshold. Our results are shown to be useful for studying the relationships between environmental fragmentation and maximum sustainable yield from populations. We present numerical results in the case of stochastic environments.

0 citations0 reviewsNo code linkedOpen access