Related papers: A simple solver for the two-fluid plasma model bas…
The Finite Difference Time Domain (FDTD) method is a widely used numerical technique for solving Maxwell's equations, particularly in computational electromagnetics and photonics. It enables accurate modeling of wave propagation in complex…
We prove existence and uniqueness of smooth solutions with large initial data for a system of equations modeling the interaction of short waves, governed by a nonlinear Schr\"odinger equation, and long waves, described by the equations of…
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO…
Studying time-dependent behavior in lasers is analytically difficult due to the saturating non-linearity inherent in the Maxwell-Bloch equations and numerically demanding because of the computational resources needed to discretize both time…
We revisit, both numerically and analytically, the finite-time blowup of the infinite-energy solution of 3D Euler equations of stagnation-point-type introduced by Gibbon et al. (1999). By employing the method of mapping to regular systems,…
We present a spectrally-accurate scheme to turn a boundary integral formulation for an elliptic PDE on a single unit cell geometry into one for the fully periodic problem. Applications include computing the effective permeability of…
Numerical methods for the Euler equations with a singular source are discussed in this paper. The stationary discontinuity induced by the singular source and its coupling with the convection of fluid presents challenges to numerical…
In this work, we detail the GPU-porting of an in-house pseudo-spectral solver tailored towards large-scale simulations of interface-resolved simulation of drop- and bubble-laden turbulent flows. The code relies on direct numerical…
We present a refinement of the Spectral Method by incorporating an optimization method into it and generalize it to two space dimensions. We then apply this Refined Spectral Method as an extremely accurate technique for finding the bound…
A novel 3-D higher-order finite-difference time-domain framework with complex frequency-shifted perfectly matched layer for the modeling of wave propagation in cold plasma is presented. Second- and fourth-order spatial approximations are…
We prove that the 3-D free-surface incompressible Euler equations with regular initial geometries and velocity fields have solutions which can form a finite-time "splash" (or "splat") singularity first introduced in [9], wherein the…
The Poisson equation is critical to get a self-consistent solution in plasma fluid simulations used for Hall effect thrusters and streamer discharges, since the Poisson solution appears as a source term of the unsteady nonlinear flow…
The Vlasov-Maxwell equations provide kinetic simulations of collisionless plasmas, but numerically solving them on classical computers is often impractical. This is due to the computational resource constraints imposed by the time evolution…
The 3-D exact analytical solutions of ideal two fluid plasma, single fluid plasma (MHD) and neutral fluid equations have been found using physically justifiable assumptions. Surprisingly these solutions satisfy all non-linearities in the…
In this article, we study the gas expansion problem by turning a sharp corner into vacuum for the two-dimensional pseudo-steady compressible Euler equations with a convex equation of state. This problem can be considered as interaction of a…
The accuracy and stability of implicit CFD codes are frequently impaired by the decoupling between variables, which can ultimately lead to numerical divergence. Coupled solvers, which solve all the governing equations simultaneously, have…
Two fundamental models in plasma physics are given by the Vlasov-Maxwell-Boltzmann system and the compressible Euler-Maxwell system which both capture the complex dynamics of plasmas under the self-consistent electromagnetic interactions at…
We develop a spectrally accurate numerical method to compute solutions of a model partial differential equation used in plasma physics to describe diffusion in velocity space due to Fokker-Planck collisions. The solution is represented as a…
For one dimensional or multidimensional compressible Euler system of polytropic gases, it is well known that the smooth solution will generally develop singularities in finite time. However, for three dimensional Chaplygin gases, due to the…
We describe a parallel algorithm for solving the time-independent 3d Schrodinger equation using the finite difference time domain (FDTD) method. We introduce an optimized parallelization scheme that reduces communication overhead between…