Related papers: Simulating Pitch Angle Scattering Using An Explici…
We develop an Explicitly Solvable Energy-Conserving (ESEC) algorithm for the Stochastic Differential Equation (SDE) describing the pitch-angle scattering process in magnetized plasmas. The Cayley transform is used to calculate both the…
A class of explicit numerical schemes is developed to solve for the relativistic dynamics and spin of particles in electromagnetic fields, using the Lorentz-BMT equation formulated in the Clifford algebra representation of Baylis. It is…
Particles and fields are standard components in numerical simulations like transport simulations in nuclear physics and have very well understood dynamics. Still, a common problem is the interaction between particles and fields due to their…
Due to manufacturing defects or wear and tear, industrial components may have uncertainties. In order to evaluate the performance of machined components, it is crucial to quantify the uncertainty of the scattering surface. This brings up an…
Thermal decay rate over an edge-shaped barrier at high dissipation is studied numerically through the computer modeling. Two sorts of the stochastic Langevin type equations are applied: (i) the Langevin equations for the coordinate and…
We discuss the electromagnetic scattering and radiation problems of multilayered spheres, reviewing the history of the Lorentz-Mie theory and the numerical stability issues encountered in handling multilayered spheres. By combining…
In this work, the lattice Boltzmann method (LBM) is assessed as a time-domain numerical approach for electromagnetic wave scattering. Owing to its explicit formulation and suitability for parallel computation on structured grids, LBM…
In this paper, we propose and analyze a positivity-preserving, energy stable numerical scheme for certain type reaction-diffusion systems involving the Law of Mass Action with the detailed balance condition. The numerical scheme is…
We investigate the use of stochastic methods for zero energy quantum scattering based on a path integral approach. With the application to the scattering of a projectile from a nuclear many body target in mind, we use the potential…
We present a method for constructing numerical schemes with up to 3rd strong convergence order for solution of a class of stochastic differential equations, including equations of the Langevin type. The construction proceeds in two stages.…
Langevin simulation provides an effective way to study collisional effects in beams by reducing the six-dimensional Fokker-Planck equation to a group of stochastic ordinary differential equations. These resulting equations usually have…
Time evolution and scattering simulation in phenomenological models are of great interest for testing and validating the potential for near-term quantum computers to simulate quantum field theories. Here, we simulate one-particle…
The Langevin dynamics is a diffusion process extensively used, in particular in molecular dynamics simulations, to sample Gibbs measures. Some alternatives based on (piecewise deterministic) kinetic velocity jump processes have gained…
In this work, we aim at constructing numerical schemes, that are as efficient as possible in terms of cost and conservation of invariants, for the Vlasov--Fokker--Planck system coupled with Poisson or Amp\`ere equation. Splitting methods…
We propose a novel iterative numerical method to solve the three-dimensional inverse obstacle scattering problem of recovering the shape of the obstacle from far-field measurements. To address the inherent ill-posed nature of the inverse…
The present work introduces a simple, yet effective particle coalescing procedure for two-dimensional SPH simulations with spatially varying resolution. In addition to the regular conservation properties of former algorithms concerning the…
This paper develops an efficient numerical method for the inverse scattering problem of a time-harmonic plane wave incident on a perfectly reflecting random periodic structure. The method is based on a novel combination of the Monte Carlo…
Lorentz invariant structure-preserving algorithms possess reference-independent secular stability, which is vital for simulating relativistic multi-scale dynamical processes. The splitting method has been widely used to construct…
Time-centered, hence second-order, methods for integrating the relativistic momentum of charged particles in an electromagnetic field are derived. A new method is found by averaging the momentum before use in the magnetic rotation term, and…
In this article, we propose a numerical approach to solve quantum mechanical scattering problems, using phase function method, by considering neutron-proton interaction as an example. The nonlinear phase equation, obtained from the…