相关论文: Two-scale exponential integrators with uniform acc…
We present a new algorithm which is named the Dynamical Functional Particle Method, DFPM. It is based on the idea of formulating a finite dimensional damped dynamical system whose stationary points are the solution to the original…
Recent developments in high peak-power table-top laser systems reaching highly relativistic light intensities have led to significant advances in laser-driven particle acceleration schemes (mainly the laser wakefield acceleration, LWFA)…
Molecular dynamics (MD) simulations of complex electrochemical systems, such as ionic liquid supercapacitors, are increasingly including the constant potential method (CPM) to model conductive electrodes at specified potential difference,…
We propose a novel algorithm for the temporal integration of the resistive magnetohydrodynamics (MHD) equations. The approach is based on exponential Rosenbrock schemes in combination with Leja interpolation. It naturally preserves Gauss's…
An electrostatic, implicit particle-in-cell (PIC) model for collisionless, fully magnetized, paraxial plasma expansions in a magnetic nozzle is introduced with exact charge, energy, and magnetic moment conservation properties. The approach…
In this paper, we propose two novel fourth-order integrators that exhibit uniformly high accuracy and long-term near conservations for solving the nonlinear Dirac equation (NLDE) in the nonrelativistic regime. In this regime, the solution…
Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear, computationally intensive, and can be difficult to interpret, while often only…
We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based…
We present the standard electromagnetic Particle-in-Cell method, starting from the discrete approximation of derivatives on a uniform grid. The application to second-order, centered, finite-difference discretization of the equations of…
Many integral equation-based methods are available for problems of time-harmonic electromagnetic scattering from perfect electric conductors. Among the many challenges that arise in such calculations are the avoidance of spurious…
The two-dimensional extension of the one-dimensional PDM-Lagrangians and their nonlocal point transformation mappings into constant unit-mass exactly solvable Lagrangians is introduced. The conditions on the related two-dimensional…
In this paper, we formulate and analyse exponential integrations when applied to nonlinear Schr\"{o}dinger equations in a normal or highly oscillatory regime. A kind of exponential integrators with energy preservation, optimal convergence…
Maxwell equations describe the propagation of electromagnetic waves and are therefore fundamental to understanding many problems encountered in the study of antennas and electromagnetics. The aim of this paper is to propose and analyse an…
The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very…
We propose a parallel version of the cross interpolation algorithm and apply it to calculate high-dimensional integrals motivated by Ising model in quantum physics. In contrast to mainstream approaches, such as Monte Carlo and quasi Monte…
Across many plasma applications, the underlying phenomena and interactions among the involved processes are known to exhibit three-dimensional characteristics. Furthermore, the global properties and evolution of plasma systems are often…
Angular integrals arise in a wide range of perturbative quantum field theory calculations. In this work we investigate angular integrals with three denominators in $d=4-2\varepsilon$ dimensions. We derive integration-by-parts relations for…
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an…
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…
A common optimization problem in the areas of magnetized plasmas and fusion energy is the design of magnets to produce a given three-dimensional magnetic field distribution to high precision. When designing arrays of permanent magnets for…