Related papers: High-order geometric integrators for the variation…
Sequential numerical methods for integrating initial value problems (IVPs) can be prohibitively expensive when high numerical accuracy is required over the entire interval of integration. One remedy is to integrate in a parallel fashion,…
We consider the systematic numerical approximation of Biot's quasistatic model for the consolidation of a poroelastic medium. Various discretization schemes have been analysed for this problem and inf-sup stable finite elements have been…
We develop a hybrid spatial discretization for the wave equation in second order form, based on high-order accurate finite difference methods and discontinuous Galerkin methods. The hybridization combines computational efficiency of finite…
This paper introduces an extension of the Morley element for approximating solutions to biharmonic equations. Traditionally limited to piecewise quadratic polynomials on triangular elements, the extension leverages weak Galerkin finite…
Numerical solving the Schr\"odinger equation with incommensurate potentials presents a great challenge since its solutions could be space-filling quasiperiodic structures without translational symmetry nor decay. In this paper, we propose…
Recently a new class of numerical integration methods -- ``mixed variable symplectic integrators'' -- has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of…
In this work, we develop energy-preserving iterative schemes for the (non-)linear systems arising in the Gauss integration of Poisson systems with quadratic Hamiltonian. Exploiting the relation between Gauss collocation integrators and…
Entropic optimal transport problems are regularized versions of optimal transport problems. These models play an increasingly important role in machine learning and generative modelling. For finite spaces, these problems are commonly solved…
The generalized-$\alpha$ method encompasses a wide range of time integrators. The method possesses high-frequency dissipation while minimizing unwanted low-frequency dissipation and the numerical dissipation can be controlled by the user.…
Exponential time differencing methods is a power tool for high-performance numerical simulation of computationally challenging problems in condensed matter physics, fluid dynamics, chemical and biological physics, where mathematical models…
A Gaussian resolution method for the computation of equilibrium density matrices rho(T) for a general multidimensional quantum problem is presented. The variational principle applied to the ``imaginary time'' Schroedinger equation provides…
Ehrenfest dynamics is a useful approximation for ab initio mixed quantum-classical molecular dynamics that can treat electronically nonadiabatic effects. Although a severe approximation to the exact solution of the molecular time-dependent…
In this paper, we propose high order numerical methods to solve a 2D advection diffusion equation, in the highly oscillatory regime. We use an integrator strategy that allows the construction of arbitrary high-order schemes {leading} to an…
Splitting methods are widely used for solving initial value problems (IVPs) due to their ability to simplify complicated evolutions into more manageable subproblems which can be solved efficiently and accurately. Traditionally, these…
The numerical simulation of the time-dependent Schr\"odinger equation for quantum systems is a very active research topic. Yet, resolving the solution sufficiently in space and time is challenging and mandates the use of modern…
We present a novel methodology for constructing arbitrarily high-order structure-preserving methods tailored for damped Hamiltonian systems. This method combines the idea of exponential integrator and energy-preserving collocation methods,…
We present a novel implicit scheme for the numerical solution of time-dependent conservation laws. The core idea of the presented method is to exploit and approximate the mixed spatial-temporal derivative of the solution that occurs…
In this paper we study the finite element approximation of systems of second-order nonlinear hyperbolic equations. The proposed numerical method combines a $hp$-version discontinuous Galerkin finite element approximation in the time…
In this paper we propose and investigate a general approach to constructing local energy-preserving algorithms which can be of arbitrarily high order in time for solving Hamiltonian PDEs. This approach is based on the temporal…
We construct a higher-order adaptive method for strong approximations of exit times of It\^o stochastic differential equations (SDE). The method employs a strong It\^o--Taylor scheme for simulating SDE paths, and adaptively decreases the…