Related papers: Hyper Boris integrators for kinetic plasma simulat…
We present a 3D finite element solver for the nonlinear Poisson-Nernst-Planck (PNP) equations for electrodiffusion, coupled to the Stokes system of fluid dynamics. The model serves as a building block for the simulation of macromolecule…
In a companion paper (hep-th/0512317), we have presented an approximation scheme to solve the Non Perturbative Renormalization Group equations that allows the calculation of the $n$-point functions for arbitrary values of the external…
Modeling and simulation of High Power Microwave (HPM) breakdown, a multiscale phenomenon, is computationally expensive and requires solving Maxwell's equations (EM solver) coupled with a plasma continuity equation (plasma solver). In this…
To approximate solutions of complex nonlinear partial differential equations remains a computational challenge, especially for sets of equations relevant in industry, such as Euler or Navier-Stokes equations. Even the most sophisticated…
We present a novel methodology for the numerical solution of problems of diffraction by infinitely thin screens in three dimensional space. Our approach relies on new integral formulations as well as associated high-order quadrature rules.…
We propose a new type of multilevel method for solving eigenvalue problems based on Newton iteration. With the proposed iteration method, solving eigenvalue problem on the finest finite element space is replaced by solving a small scale…
We present sixth- and eighth-order Hermite integrators for astrophysical $N$-body simulations, which use the derivatives of accelerations up to second order ({\it snap}) and third order ({\it crackle}). These schemes do not require previous…
We present a new solver for nonlinear parabolic problems that is L-stable and achieves high order accuracy in space and time. The solver is built by first constructing a single-dimensional heat equation solver that uses fast O(N)…
In [Z. Hu, R. Li, and Z. Qiao. Acceleration for microflow simulations of high-order moment models by using lower-order model correction. J. Comput. Phys., 327:225-244, 2016], it has been successfully demonstrated that using lower-order…
We develop a highly accurate analytic approximation for small-scale non-cold relic perturbations by solving the collisionless Boltzmann equation in the quasi-stationary regime. The approximation is implemented in CLASSIER (CLASS Integral…
The accurate prediction of occurrence and strength of kinetic instabilities in plasmas remains a significant challenge in nuclear fusion research. To accurately capture the plasmas dynamics one is required to solve the Vlasov equation for…
Superlevel partitioning is combined with a simple relaxation procedure to construct an iterative technique for solving equations of statistical equilibrium. In treating an $N$-level model atom, the technique avoids the $N^{3}$ scaling in…
It is demonstrated that the inherent structure of Boris pusher may result in the accumulation of errors in numerical integration when estimating the works of electric fields on particles. This, in turn, leads to an incorrect estimations of…
The simulation of turbulence in the boundary region of a tokamak is crucial for understanding and optimizing the performance of fusion reactors. In this work, the use of low-rank linear algebra techniques is shown to enhance the efficiency…
We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled to hydrodynamics equations. It is meant to be applied to simulations of core-collapse supernovae. We handle…
Particle-in-cell codes are now standard tools for studying ultra-intense laser-plasma interactions. Motivated by direct laser acceleration of electrons in sub-critical plasmas, we examine temporal resolution requirements that must be…
We describe an algorithm for the numerical solution of second order linear differential equations in the highly-oscillatory regime. It is founded on the recent observation that the solutions of equations of this type can be accurately…
In this paper we present a technique for constructing robust solvers for stiff algebraic source terms, such as those typically used for modelling relaxation processes in hyperbolic systems of partial differential equations describing…
In this paper, we present an iterative three-point method with memory based on the family of King's methods to solve nonlinear equations. This proposed method has eighth order convergence and costs only four function evaluations per…
We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler…