Related papers: A note on the convergence of parametrised non-reso…
In this work we develop a method for computing mathematically rigorous enclosures of some one dimensional manifolds of heteroclinic orbits for nonlinear maps. Our method exploits a rigorous curve following argument build on high order…
We develop a validated numerical procedure for continuation of local stable/unstable manifold patches attached to equilibrium solutions of ordinary differential equations. The procedure has two steps. First we compute an accurate high order…
This paper develops a computational method for studying stable/unstable manifolds attached to periodic orbits of differential equations. The method uses high order Chebyshev-Taylor series approximations in conjunction with the…
Inverse problems use physical measurements along with a computational model to estimate the parameters or state of a system of interest. Errors in measurements and uncertainties in the computational model lead to inaccurate estimates. This…
The a posteriori error estimates are studied for a class of nonlinear stead-state Poisson-Nernst-Planck equations, which are a coupled system consisting of the Nernst-Planck equation and the Poisson equation. Both the global upper bounds…
We study the computation of local approximations of invariant manifolds of parabolic fixed points and parabolic periodic orbits of periodic vector fields. If the dimension of these manifolds is two or greater, in general, it is not possible…
We consider a family of parallel methods for constrained optimization based on projected gradient descents along individual coordinate directions. In the case of polyhedral feasible sets, local convergence towards a regular solution occurs…
One of the very few mathematically rigorous nonlinear model reduction methods is the restriction of a dynamical system to a low-dimensional, sufficiently smooth, attracting invariant manifold. Such manifolds are usually found using local…
The parametrisation method for invariant manifolds is a powerful technique for deriving reduced-order models in the context of nonlinear vibrating systems, allowing accurate computations of nonlinear normal modes. Thanks to arbitrary order…
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…
We consider the proximal gradient method on Riemannian manifolds for functions that are possibly not geodesically convex. Starting from the forward-backward-splitting, we define an intrinsic variant of the proximal gradient method that uses…
Based on the auxiliary subspace techniques, a hierarchical basis a posteriori error estimator is proposed for the Stokes problem in two and three dimensions. For the error estimator, we need to solve only two global diagonal linear systems…
It is known that the asymptotic invariant manifolds around an unstable periodic orbit in conservative systems can be represented by convergent series (Cherry 1926, Moser 1956, 1958, Giorgilli 2001). The unstable and stable manifolds…
In this work, we obtain an a-posteriori theorem for the existence of partly hyperbolic invariant tori in analytic Hamiltonian systems: autonomous, periodic, and quasi-periodic. The method of proof is based on the convergence of a KAM…
This work studies a posteriori error estimates and their use for time-dependent acoustic scattering problems, formulated as a time-dependent boundary integral equation based on a single-layer ansatz. The integral equation is discretized by…
In this paper we study the existence and regularity of stable manifolds associated to fixed points of parabolic type in the differentiable and analytic cases, using the parametrization method. The parametrization method relies on a suitable…
We present and implement an algorithm for computing the invariant circle and the corresponding stable manifolds for 2-dimensional maps. The algorithm is based on the parameterization method, and it is backed up by an a-posteriori theorem…
One approach to parametric and adaptive model reduction is via the interpolation of orthogonal bases, subspaces or positive definite system matrices. In all these cases, the sampled inputs stem from matrix sets that feature a geometric…
We develop new adaptive algorithms for temporal integration of nonlinear evolution equations on tensor manifolds. These algorithms, which we call step-truncation methods, are based on performing one time step with a conventional…
Nonuniform tubular neighborhoods of curves in Euclidean n-space are studied by using weighted distance functions and generalizing the normal exponential map. Different notions of injectivity radii are introduced to investigate singular but…