Related papers: Simulating Pitch Angle Scattering Using An Explici…
Electromagnetic wave scattering from planar dielectric films deposited on one-dimensional, randomly rough, perfectly conducting substrates is studied by numerical simulations for both p- and s-polarization. The reduced Rayleigh equation,…
Transport and acceleration of charged particles in turbulent media is a topic of great interest in space physics and interstellar astrophysics. These processes are dominated by the scattering of particles off magnetic irregularities. The…
A formalism is introduced for the non-perturbative, purely numerical, solution of the reduced Rayleigh equation for the scattering of light from two-dimensional penetrable rough surfaces. As an example, we apply this formalism to study the…
In this article, we study the problem of sampling from distributions whose densities are not necessarily smooth nor logconcave. We propose a simple Langevin-based algorithm that does not rely on popular but computationally challenging…
We present an iterative algorithm to compute numerical approximations of the potential for the Schr\"odinger operator from scattering data. Four different types of scattering data are used as follows: fixed energy, fixed incident angle,…
A finite difference numerical scheme is proposed and analyzed for the Cahn-Hilliard-Stokes system with Flory-Huggins energy functional. A convex splitting is applied to the chemical potential, which in turns leads to the implicit treatment…
In this paper, two finite difference numerical schemes are proposed and analyzed for the droplet liquid film model, with a singular Leonard-Jones energy potential involved. Both first and second order accurate temporal algorithms are…
A numerical algorithm is presented for solving the direct scattering problems by the Modified Rayleigh Conjecture Method (MRC) introduced by A.G.Ramm. Some numerical examples are given. They show that the method is numerically efficient.
In the present manuscript, we consider the problem of dispersive wave simulation on a rotating globally spherical geometry. In this Part IV, we focus on numerical aspects while the model derivation was described in Part III. The algorithm…
The behavior of a Lattice Monte Carlo algorithm (if it is designed correctly) must approach that of the continuum system that it is designed to simulate as the time step and the mesh step tend to zero. However, we show for an algorithm for…
The Lorentz equations describe the motion of electrically charged particles in electric and magnetic fields and are used widely in plasma physics. The most popular numerical algorithm for solving them is the Boris method, a variant of the…
In electromagnetism, acoustics, and quantum mechanics, scattering problems can routinely be solved numerically by virtue of perfectly matched layers (PMLs) at simulation domain boundaries. Unfortunately, the same has not been possible for…
In this paper we study the problem of computing the effective diffusivity for a particle moving in chaotic and stochastic flows. In addition we numerically investigate the residual diffusion phenomenon in chaotic advection. The residual…
We introduce a new iterative method to recover a real compact supported potential of the Schr\"odinger operator from their fixed angle scattering data. The method combines a fixed point argument with a suitable approximation of the…
We formulate a smoothed-particle hydrodynamics numerical method, traditionally used for the Euler equations for fluid dynamics in the context of astrophysical simulations, to solve the non-linear Schrodinger equation in the Madelung…
We propose a novel formulation for parametric finite element methods to simulate surface diffusion of closed curves, which is also called as the curve diffusion. Several high-order temporal discretizations are proposed based on this new…
We propose a Langevin diffusion-based algorithm for non-convex optimization and sampling on a product manifold of spheres. Under a logarithmic Sobolev inequality, we establish a guarantee for finite iteration convergence to the Gibbs…
We propose a new simple and explicit numerical scheme for time-homogeneous stochastic differential equations. The scheme is based on sampling increments at each time step from a skew-symmetric probability distribution, with the level of…
This paper presents a robust numerical solution to the electromagnetic scattering problem involving multiple multi-layered cavities in both transverse magnetic and electric polarizations. A transparent boundary condition is introduced at…
A new energy and enstrophy conserving scheme is evaluated using a suite of test cases over the global spherical domain or bounded domains. The evaluation is organized around a set of pre-defined properties: accuracy of individual opeartors,…