Source author record

Maxime Breden

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

8 published item(s)

preprint2026arXiv

A simple rigorous integrator for semilinear parabolic PDEs

Simulations of the dynamics generated by partial differential equations (PDEs) provide approximate, numerical solutions to initial value problems. Such simulations are ubiquitous in scientific computing, but the correctness of the results is usually not guaranteed. We propose a new method for the rigorous integration of parabolic PDEs, i.e., the derivation of rigorous and explicit error bounds between the numerically obtained approximate solution and the exact one, which is then proven to exist over the entire time interval considered. These guaranteed error bounds are obtained a posteriori, using a fixed point reformulation based on a piece-wise in time constant approximation of the linearization around the numerical solution. Our setup leads to relatively simple-to-understand estimates, which has several advantages. Most critically, it allows us to optimize various aspects of the proof, and in particular to provide an adaptive time-stepping strategy. In case the solution converges to a stable hyperbolic equilibrium, we are also able to prove this convergence, applying our rigorous integrator with a final, infinitely long timestep. We showcase the ability of our method to rigorously integrate over relatively long time intervals, and to capture non-trivial dynamics, via examples on the Swift--Hohenberg equation, the Ohta--Kawasaki equation and the Kuramoto--Sivashinsky equation. We expect that the simplicity and efficiency of the approach will enable generalization to a wide variety of other parabolic PDEs, as well as applications to boundary value problems.

preprint2026arXiv

Constructive proofs for some semilinear PDEs on $H^2(e^{|x|^2/4},\mathbb{R}^d)$

We develop computer-assisted tools to study semilinear equations of the form \begin{equation*} -Δu -\frac{x}{2}\cdot \nabla{u}= f(x,u,\nabla u) ,\quad x\in\mathbb{R}^d. \end{equation*} Such equations appear naturally in several contexts, and in particular when looking for self-similar solutions of parabolic PDEs. We develop a general methodology, allowing us not only to prove the existence of solutions, but also to describe them very precisely. We introduce a spectral approach based on an eigenbasis of $\mathcal{L}:= -Δ-\frac{x}{2}\cdot \nabla$ in spherical coordinates, together with a quadrature rule allowing to deal with nonlinearities, in order to get accurate approximate solutions. We then use a Newton-Kantorovich argument, in an appropriate weighted Sobolev space, to prove the existence of a nearby exact solution. We apply our approach to nonlinear heat equations, to nonlinear Schrödinger equations and to a generalised viscous Burgers equation, and obtain both radial and non-radial self-similar profiles.

preprint2022arXiv

A posteriori validation of generalized polynomial chaos expansions

Generalized polynomial chaos expansions are a powerful tool to study differential equations with random coefficients, allowing in particular to efficiently approximate random invariant sets associated to such equations. In this work, we use ideas from validated numerics in order to obtain rigorous a posteriori error estimates together with existence results about gPC expansions of random invariant sets. This approach also provides a new framework for conducting validated continuation, i.e. for rigorously computing isolated branches of solutions in parameter-dependent systems, which generalizes in a straightforward way to multi-parameter continuation. We illustrate the proposed methodology by rigorously computing random invariant periodic orbits in the Lorenz system, as well as branches and 2-dimensional manifolds of steady states of the Swift-Hohenberg equation.

preprint2021arXiv

Computer-assisted proofs for some nonlinear diffusion problems

In the last three decades, powerful computer-assisted techniques have been developed in order to validate a posteriori numerical solutions of semilinear elliptic problems of the form $Δu +f(u,\nabla u) = 0$. By studying a well chosen fixed point problem defined around the numerical solution, these techniques make it possible to prove the existence of a solution in an explicit (and usually small) neighborhood the numerical solution. In this work, we develop a similar approach for a broader class of systems, including nonlinear diffusion terms of the form $ΔΦ(u)$. In particular, this enables us to obtain new results about steady states of a cross-diffusion system from population dynamics: the (non-triangular) SKT model. We also revisit the idea of automatic differentiation in the context of computer-assisted proof, and propose an alternative approach based on differential-algebraic equations.

preprint2015arXiv

Computation of maximal local (un)stable manifold patches by the parameterization method

In this work we develop some automatic procedures for computing high order polynomial expansions of local (un)stable manifolds for equilibria of differential equations. Our method incorporates validated truncation error bounds, and maximizes the size of the image of the polynomial approximation relative to some specified constraints. More precisely we use that the manifold computations depend heavily on the scalings of the eigenvectors: indeed we study the precise effects of these scalings on the estimates which determine the validated error bounds. This relationship between the eigenvector scalings and the error estimates plays a central role in our automatic procedures. In order to illustrate the utility of these methods we present several applications, including visualization of invariant manifolds in the Lorenz and FitzHugh-Nagumo systems and an automatic continuation scheme for (un)stable manifolds in a suspension bridge problem. In the present work we treat explicitly the case where the eigenvalues satisfy a certain non-resonance condition.

preprint2015arXiv

Global bifurcation diagram of steady states of systems of PDEs via rigorous numerics: a 3-component reaction-diffusion system

In this paper, we use rigorous numerics to compute several global smooth branches of steady states for a system of three reaction-diffusion PDEs introduced by Iida et al. [J. Math. Biol., {\bf 53}, 617--641 (2006)] to study the effect of cross-diffusion in competitive interactions. An explicit and mathematically rigorous construction of a global bifurcation diagram is done, except in small neighborhoods of the bifurcations. The proposed method, even though influenced by the work of van den Berg et al. [Math. Comp., {\bf 79}, 1565--1584 (2010)], introduces new analytic estimates, a new {\em gluing-free} approach for the construction of global smooth branches and provides a detailed analysis of the choice of the parameters to be made in order to maximize the chances of performing successfully the computational proofs.

preprint2015arXiv

Rigorous numerics for nonlinear operators with tridiagonal dominant linear part

We present a method designed for computing solutions of infinite dimensional non linear operators $f(x) = 0$ with a tridiagonal dominant linear part. We recast the operator equation into an equivalent Newton-like equation $x = T(x) = x - Af(x)$, where $A$ is an approximate inverse of the derivative $Df(\overline x)$ at an approximate solution $\overline x$. We present rigorous computer-assisted calculations showing that $T$ is a contraction near $\overline x$, thus yielding the existence of a solution. Since $Df(\overline x)$ does not have an asymptotically diagonal dominant structure, the computation of $A$ is not straightforward. This paper provides ideas for computing $A$, and proposes a new rigorous method for proving existence of solutions of nonlinear operators with tridiagonal dominant linear part.

preprint2015arXiv

Smoothness of moments of the solutions of discrete coagulation equations with diffusion

In this paper, we establish smoothness of moments of the solutions of discrete coagulation-diffusion systems. As key assumptions, we suppose that the coagulation coefficients grow at most sub-linearly and that the diffusion coefficients converge towards a strictly positive limit (those conditions also imply the existence of global weak solutions and the absence of gelation).