Related papers: Highly Parallel Demagnetization Field Calculation …
The basic idea of fast Monte Carlo (MC) simulations is to perform particle-based MC simulations with the excluded-volume interactions modeled by "soft" repulsive potentials that allow particle overlapping. This gives much faster system…
The plasma edge flow, situated at the intricate boundary between plasma and neutral particles, plays a pivotal role in the design of nuclear fusion devices such as divertors and pumps. Traditional numerical simulation methods, such as the…
This article introduces a highly parallel algorithm for molecular dynamics simulations with short-range forces on single node multi- and many-core systems. The algorithm is designed to achieve high parallel speedups for strongly…
Accurately measuring magnetic field strength in the interstellar medium, including giant molecular clouds (GMCs), remains a significant challenge. We present a machine learning approach using Denoising Diffusion Probabilistic Models (DDPMs)…
We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms…
Fast multipole methods have O(N) complexity, are compute bound, and require very little synchronization, which makes them a favorable algorithm on next-generation supercomputers. Their most common application is to accelerate N-body…
The finite amplitude method (FAM) is a very efficient approach for solving the fully self-consistent random-phase approximation (RPA) equations. We use FAM to rederive the RPA matrices for general Skyrme-like functionals, calculate the…
The validation, verification, and uncertainty quantification of computationally expensive theoretical models of quantum many-body systems require the construction of fast and accurate emulators. In this work, we develop emulators for…
The oversampling multiscale finite element method (MsFEM) is one of the most popular methods for simulating composite materials and flows in porous media which may have many scales. But the method may be inapplicable or inefficient in some…
We propose extensions and improvements of the statistical analysis of distributed multipoles (SADM) algorithm put forth by Chipot et al. in [6] for the derivation of distributed atomic multipoles from the quantum-mechanical electrostatic…
The magnetic quadrupole structure formation during the interaction of two ultra-short high power laser pulses with a collisionless plasma is demonstrated with 2.5-dimensional particle-in-cell simulations. The subsequent expansion of the…
This paper presents the first parallel implementation of the novel "Interpolated Factored Green Function" (IFGF) method introduced recently for the accelerated evaluation of discrete integral operators arising in wave scattering and other…
We extend the foil winding homogenization method to magnetic field conforming formulations. We first propose a full magnetic field foil winding formulation by analogy with magnetic flux density conforming formulations. We then introduce the…
A fundamental task in particle-in-cell (PIC) simulations of plasma physics is solving for charged particle motion in electromagnetic fields. This problem is especially challenging when the plasma is strongly magnetized due to numerical…
The FHP algorithm allows to obtain the relativistic multipole moments of a vacuum stationary axisymmetric solution in terms of coefficients which appear in the expansion of its Ernst's potential on the symmetry axis. First of all, we will…
X-ray pulse-profile modeling of millisecond pulsars offers a direct route to measuring neutron star masses and radii, thereby constraining the dense-matter equation of state. However, standard analyses typically rely on \emph{ad hoc}…
We develop a multiscale hybrid scheme for simulations of soft condensed matter systems, which allows one to treat the system at the particle level in selected regions of space, and at the continuum level elsewhere. It is derived…
An efficient parallelization approach to simulate optical properties of ensembles of quantum emitters in realistic electromagnetic environments is considered. It relies on balancing computing load of utilized processors and is built into…
We present an indirect higher order boundary element method utilising NURBS mappings for exact geometry representation and an interpolation-based fast multipole method for compression and reduction of computational complexity, to counteract…
Different numerical approaches for the stray-field calculation in the context of micromagnetic simulations are investigated. We compare finite difference based fast Fourier transform methods, tensor grid methods and the finite-element…