Related papers: A New Time-Dependent Finite Difference Method for …
Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different…
A systematic and robust laser speckle contrast imaging (LSCI) method and procedure is presented, covering the LSCI system calibration, static scattering removal, and measurement noise estimation and correction to obtain a true flow speckle…
A particle scheme for scalar conservation laws in one space dimension is presented. Particles representing the solution are moved according to their characteristic velocities. Particle interaction is resolved locally, satisfying exact…
Scale-resolving simulations of high Reynolds number incompressible flows are often limited by the Courant-Friedrichs-Lewy (CFL) stability restriction imposed by explicit time-stepping schemes, resulting in small time step sizes and long…
The DistFlow model accurately represents power flows in distribution systems, but the model's nonlinearities result in computational challenges for many applications. Accordingly, a linear approximation known as \mbox{LinDistFlow} (and its…
We present a new algorithm for the discretization of the Vlasov-Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a…
We present exact, analytic and simple solutions of relativistic perfect fluid hydrodynamics. The solutions allow us to calculate the rapidity distribution of the particles produced at the freeze-out, and fit them to the measured rapidity…
We examine the problem of particle acceleration at a relativistic shocks assuming pitch-angle scattering and using a Hartree-Fock method to approximate the associated eigenfunctions. This leads to a simple transcendental equation…
The Stochastic Weighted Particle Method (SWPM) of Rjasanow and Wagner is a generalization of the Direct Simulation Monte Carlo method for computing the probability density function of the velocities of a system of interacting particles for…
We present a new family of relativistic lattice kinetic schemes for the efficient simulation of relativistic flows in both strongly-interacting (fluid) and weakly-interacting (rarefied gas) regimes. The method can also deal with both…
In order to prevent velocity, pressure, and temperature spikes at material discontinuities occurring when the interface-capturing schemes inconsistently simulate compressible multi-material flows(when the specific heats ratio is…
In this paper, we investigate the recently developed lattice Boltzmann model for relativistic hydrodynamics. To this purpose, we perform simulations of shock waves in quark-gluon plasma in the low and high viscosities regime, using three…
The acceleration time scale for the process of first-order Fermi acceleration in relativistic shock waves with oblique magnetic field configurations is investigated by the method of Monte Carlo particle simulations. We demonstrate the…
Explicit numerical finite difference schemes for partial differential equations are well known to be easy to implement but they are particularly problematic for solving equations whose solutions admit shocks, blowups and discontinuities.…
This paper presents a tensor-recovery method to solve probabilistic power flow problems. Our approach generates a high-dimensional and sparse generalized polynomial-chaos expansion that provides useful statistical information. The result…
A novel explicit method to model Lorentz linear dispersive media with finite difference method are presented. The method shows an explicit method without any modification to the Leap-Frogging scheme. The polarizations of the Lorentz media…
In this paper, a novel non-intrusive probabilistic power flow (PPF) analysis method based on the low-rank approximation (LRA) is proposed, which can accurately and efficiently estimate the probabilistic characteristics (e.g., mean,…
Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of…
In plasma edge simulations, kinetic Monte Carlo (MC) is often used to simulate neutral particles and estimate source terms. For large-sized reactors, like ITER and DEMO, high particle collision rates lead to a substantial computational cost…
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…