Related papers: Hyper Boris integrators for kinetic plasma simulat…
This paper presents an efficient Krylov subspace iterative solver for the three-dimensional (3D) Helmholtz equation with non-constant coefficients and absorbing boundary conditions, combining high-resolution compact schemes with low-order…
Micromagnetics depends on high-fidelity numerical methods for magnetization dynamics. This work proposes a third-order temporal accuracy scheme for the Landau-Lifshitz-Gilbert equation, addressing accuracy-efficiency trade-offs in existing…
Among the family of fourth-order time integration schemes, the two-stage Gauss--Legendre method, which is an implicit Runge--Kutta method based on collocation, is the only superconvergent. The computational cost of this implicit scheme for…
The efficient numerical solution of many kinetic models in plasma physics is impeded by the stiffness of these systems. Exponential integrators are attractive in this context as they remove the CFL condition induced by the linear part of…
Continuum Vlasov simulations can be utilized for highly accurate modelling of fully kinetic plasmas. Great progress has been made recently regarding the applicability of the method in realistic plasma configurations. However, a reduction of…
This paper serves to treat boundary conditions numerically with high order accuracy in order to match the two-stage fourth-order finite volume schemes for hyperbolic problems developed in [{\em J. Li and Z. Du, A two-stage fourth order…
In this work, we develop first-order (Hessian-free) and zero-order (derivative-free) implementations of the Cubically regularized Newton method for solving general non-convex optimization problems. For that, we employ finite difference…
We consider the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients (methods with real coefficients of order greater than two necessarily have negative…
The treatment of hypercubic lattice artifacts is essential for the calculation of non-perturbative renormalization constants of RI-MOM schemes. It has been shown that for the RI'-MOM scheme a large part of these artifacts can be calculated…
Abridged. Kinetic plasma simulations are nowadays commonly used to study a wealth of non-linear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in…
We present a general, high-order, fully explicit relaxation scheme which can be applied to any system of nonlinear hyperbolic conservation laws in multiple dimensions. The scheme consists of two steps. In a first (relaxation) step, the…
A novel linear-algebraic algorithm, multiple Arnoldi method, was developed in an interdisciplinary study between physics and applied mathematics and realized one-hundred-million-atom (100-nm-scale) electronic state calculations on the K…
Scientific computing is an essential tool for scientific discovery and engineering design, and its computational cost is always a main concern in practice. To accelerate scientific computing, it is a promising approach to use machine…
We propose a new numerical approach to compute nonclassical solutions to hyperbolic conservation laws. The class of finite difference schemes presented here is fully conservative and keep nonclassical shock waves as sharp interfaces,…
The electromagnetic behavior of plasmonic structures can be predicted after discretizing and solving a linear system of equations, derived from a continuous surface integral equation (SIE) and the appropriate boundary conditions, using a…
Blood flow, dam or ship construction and numerous other problems in biomedical and general engineering involve incompressible flows interacting with elastic structures. Such interactions heavily influence the deformation and stress states…
While quantum computing provides an exponential advantage in solving system of linear equations, there is little work to solve system of nonlinear equations with quantum computing. We propose quantum Newton's method (QNM) for solving…
This paper presents a fast algorithm to solve a spectral estimation problem for two-dimensional random fields. The latter is formulated as a convex optimization problem with the Itakura-Saito pseudodistance as the objective function subject…
Reduced-order simulation is an emerging method for accelerating physical simulations with high DOFs, and recently developed neural-network-based methods with nonlinear subspaces have been proven effective in diverse applications as more…
We describe a number of recently developed techniques for improving the performance of large-scale nuclear configuration interaction calculations on high performance parallel computers. We show the benefit of using a preconditioned block…