Related papers: Exponential Integrators for Resistive Magnetohydro…
We introduce a novel structure-preserving method in order to approximate the compressible ideal Magnetohydrodynamics (MHD) equations. This technique addresses the MHD equations using a non-divergence formulation, where the contributions of…
Leja points on a compact $K \subset \mathbb{C}$ are known to provide efficient points for interpolation, but their actual implementation can be computationally challenging. So-called pseudo Leja points are a more tractable solution, yet…
Hamiltonian Monte Carlo (HMC) improves the computational efficiency of the Metropolis algorithm by reducing its random walk behavior. Riemannian Manifold HMC (RMHMC) further improves HMC's performance by exploiting the geometric properties…
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square…
We present a novel gradient regularization to completely eliminate the magnetic divergence error in meshless magnetohydrodynamics (MHD), which offers a high spatial resolution and conservative advantage, due to its Lagrangian nature.…
We propose explicit symplectic integrators of molecular dynamics (MD) algorithms for rigid-body molecules in the canonical and isothermal-isobaric ensembles. We also present a symplectic algorithm in the constant normal pressure and lateral…
Directional interpolation is a fast and efficient compression technique for high-frequency Helmholtz boundary integral equations, but it requires a very large amount of storage in its original form. Algebraic recompression can significantly…
The effort to generate matrix exponentials and associated differentials, required to determine the time evolution of quantum systems, frequently constrains the evaluation of problems in quantum control theory, variational circuit…
In this paper, we establish the non-uniqueness of solutions to the ideal magnetohydrodynamics equations in any dimension greater than three by proving the existence of infinitely many compactly supported weak solutions. In particular, these…
In this paper, we consider the task of efficiently computing the numerical solution of evolutionary complex Ginzburg--Landau equations on Cartesian product domains with homogeneous Dirichlet/Neumann or periodic boundary conditions. To this…
One of the most accurate methods for solving the time-dependent Schr\"{o}dinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid…
Geometric integrators of the Schr\"{o}dinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to…
Modern simulation codes for general relativistic ideal magnetohydrodynamics are all facing a long standing technical problem given by the need to recover fundamental variables from those variables that are evolved in time. In the…
A variety of high-energy astrophysical phenomena are powered by the release -- via magnetic reconnection -- of the energy stored in oppositely directed fields. Single-fluid resistive magnetohydrodynamic (MHD) simulations with uniform…
In this paper, we introduce novel fast matrix inversion algorithms that leverage triangular decomposition and recurrent formalism, incorporating Strassen's fast matrix multiplication. Our research places particular emphasis on triangular…
Multiderivative time integrators have a long history of development for ordinary differential equations, and yet to date, only a small subset of these methods have been explored as a tool for solving partial differential equations (PDEs).…
The simulation of multi-body systems with frictional contacts is a fundamental tool for many fields, such as robotics, computer graphics, and mechanics. Hard frictional contacts are particularly troublesome to simulate because they make the…
This paper studies the nonlinear evolution of magnetic field turbulence in proximity of steady ideal MHD configurations characterized by a small electric current, a small plasma flow, and approximate flux surfaces, a physical setting that…
This paper is concerned with the development and testing of advanced time-stepping methods suited for the integration of time-accurate, real-world applications of computational fluid dynamics (CFD). The performance of several time…
By leveraging the natural geometry of a smooth probabilistic system, Hamiltonian Monte Carlo yields computationally efficient Markov Chain Monte Carlo estimation. At least provided that the algorithm is sufficiently well-tuned. In this…