Related papers: Stable 3-level leapfrog integration in numerical r…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…
Intrusive Uncertainty Quantification methods such as stochastic Galerkin are gaining popularity, whereas the classical stochastic Galerkin approach is not ensured to preserve hyperbolicity of the underlying hyperbolic system. We apply a…
This work is concerned with practical stabilization of nonlinear systems by means of inf-convolution-based sample-and-hold control. It is a fairly general stabilization technique based on a generic non-smooth control Lyapunov function (CLF)…
This paper concerns models and convergence principles for dealing with stochasticity in a wide range of algorithms arising in nonlinear analysis and optimization in Hilbert spaces. It proposes a flexible geometric framework within which…
In this paper, we put forth distributed algorithms for solving loosely coupled unconstrained and constrained optimization problems. Such problems are usually solved using algorithms that are based on a combination of decomposition and first…
Linear Parameter-Varying (LPV) systems with jumps and piecewise differentiable parameters is a class of hybrid LPV systems for which no tailored stability analysis and stabilization conditions have been obtained so far. We fill this gap…
With the advent of massive data sets much of the computational science and engineering community has moved toward data-intensive approaches in regression and classification. However, these present significant challenges due to increasing…
This study concerns the efficiency of time-spectral methods for numerical solution of differential equations. It is found that the time-spectral method GWRM demonstrates insensitivity to stiffness and chaoticity due to the implicit nature…
Kinetic instabilities are one of the most challenging aspects in computational plasma physics. Accurately capturing their onset and evolution requires fine resolution of the high-dimensional distribution functions of each relevant species,…
This paper is concerned with a space-time adaptive numerical method for instationary porous media flows with nonlinear interaction between porosity and pressure, with focus on problems with discontinuous initial porosities. A convergent…
The stabilization of nonautonomous parabolic equations is achieved by feedback inputs tuning a finite number of actuators, where it is assumed that the input is subject to a time delay. To overcome destabilizing effects of the time delay,…
We consider entropy conservative and dissipative discretizations of nonlinear conservation laws with implicit time discretizations and investigate the influence of iterative methods used to solve the arising nonlinear equations. We show…
We propose an end-to-end quantum algorithm to simulate rapidly distorted turbulence via linear combination of Hamiltonian (LCHS). The algorithm comprises three primary stages: the efficient preparation of an initial turbulent state with a…
The Newton method is a powerful optimization algorithm, valued for its rapid local convergence and elegant geometric properties. However, its theoretical guarantees are usually limited to convex problems. In this work, we ask whether…
This paper studies multistep methods for the integration of reversible dynamical systems, with particular emphasis on the planar Kepler problem. It has previously been shown by Cano & Sanz-Serna that reversible linear multisteps for…
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
We present an efficient discontinuous Galerkin scheme for simulation of the incompressible Navier-Stokes equations including laminar and turbulent flow. We consider a semi-explicit high-order velocity-correction method for time integration…
In recent decades, there have been many attempts to construct symplectic integrators with variable time steps, with rather disappointing results. In this paper we identify the causes for this lack of performance, and find that they fall…
A time-stepping scheme with adaptivity in both the step size and the integration order is presented in the context of non-equilibrium dynamics described via Kadanoff-Baym equations. The accuracy and effectiveness of the algorithm are…
When training the parameters of a linear dynamical model, the gradient descent algorithm is likely to fail to converge if the squared-error loss is used as the training loss function. Restricting the parameter space to a smaller subset and…