Related papers: A Scheme to Numerically Evolve Data for the Confor…
This article is concerned with the derivation of numerical reconstruction schemes for the inverse moving source problem on determining source profiles in (time-fractional) evolution equations. As a continuation of the theoretical result on…
The numerical solution of partial differential equations (PDEs) is challenging because of the need to resolve spatiotemporal features over wide length and timescales. Often, it is computationally intractable to resolve the finest features…
We propose in this paper efficient first/second-order time-stepping schemes for the evolutional Navier-Stokes-Nernst-Planck-Poisson equations. The proposed schemes are constructed using an auxiliary variable reformulation and sophisticated…
This paper is concerned with the adaptive numerical treatment of stochastic partial differential equations. Our method of choice is Rothe's method. We use the implicit Euler scheme for the time discretization. Consequently, in each step, an…
An unsteady problem is considered for a space-fractional equation in a bounded domain. A first-order evolutionary equation involves a fractional power of an elliptic operator of second order. Finite element approximation in space is…
We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real-time using higher-order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this…
We propose a general strategy to discretize the Dyson series without applying direct numerical quadrature to high-dimensional integrals, and extend this framework to open quantum systems. The resulting discretization can also be interpreted…
An unsteady problem is considered for a space-fractional equation in a bounded domain. A first-order evolutionary equation involves the square root of an elliptic operator of second order. Finite element approximation in space is employed.…
The numerical solution of time-dependent radiative transfer problems is challenging, both, due to the high dimension as well as the anisotropic structure of the underlying integro-partial differential equation. In this paper we propose a…
Convergence of a full discretization of a second order stochastic evolution equation with nonlinear damping is shown and thus existence of a solution is established. The discretization scheme combines an implicit time stepping scheme with…
We propose a new numerical scheme of evolution for the Einstein equations using the discrete variational derivative method (DVDM). We derive the discrete evolution equation of the constraint using this scheme and show the constraint…
Initial data are the starting point for any numerical simulation. In the case of numerical relativity, Einstein's equations constrain our choices of these initial data. We will examine several of the formalisms used for specifying Cauchy…
Numerical schemes for Einstein's vacuum equation are developed. Einstein's equation in harmonic gauge is second order symmetric hyperbolic. It is discretized in four-dimensional spacetime by Finite Differences, Finite Elements, and Interior…
This is the first in a series of papers on the construction and validation of a three-dimensional code for general relativistic hydrodynamics, and its application to general relativistic astrophysics. This paper studies the consistency and…
A general framework for the numerical approximation of evolution problems is presented that allows to preserve exactly an underlying Hamiltonian- or gradient structure. The approach relies on rewriting the evolution problem in a particular…
Numerous Optimization Algorithms have a time-varying update rule thanks to, for instance, a changing step size, momentum parameter or, Hessian approximation. In this paper, we apply unrolled or automatic differentiation to a time-varying…
Hamiltonian systems are known to conserve the Hamiltonian function, which describes the energy evolution over time. Obtaining a numerical spatio-temporal scheme that accurately preserves the discretized Hamiltonian function is often a…
We review recent efforts to re-formulate the Einstein equations for fully relativistic numerical simulations. The so-called numerical relativity (computational simulations in general relativity) is a promising research field matching with…
We propose a computer-assisted approach to studying the effective continuum behavior of spatially discrete evolution equations. The advantage of the approach is that the "coarse model" (the continuum, effective equation) need not be…
We compare the performance of several discretizations of the simple pendulum equation in a series of numerical experiments. The stress is put on the long-time behaviour. We choose for the comparison numerical schemes which preserve the…