270 research outputs found
Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator
In this paper, we describe a new class of fast solvers for separable elliptic
partial differential equations in cylindrical coordinates with
free-space radiation conditions. By combining integral equation methods in the
radial variable with Fourier methods in and , we show that
high-order accuracy can be achieved in both the governing potential and its
derivatives. A weak singularity arises in the Fourier transform with respect to
that is handled with special purpose quadratures. We show how these solvers
can be applied to the evaluation of the Coulomb collision operator in kinetic
models of ionized gases.Comment: 20 pages, 5 figure
Landau Collision Integral Solver with Adaptive Mesh Refinement on Emerging Architectures
The Landau collision integral is an accurate model for the small-angle
dominated Coulomb collisions in fusion plasmas. We investigate a high order
accurate, fully conservative, finite element discretization of the nonlinear
multi-species Landau integral with adaptive mesh refinement using the PETSc
library (www.mcs.anl.gov/petsc). We develop algorithms and techniques to
efficiently utilize emerging architectures with an approach that minimizes
memory usage and movement and is suitable for vector processing. The Landau
collision integral is vectorized with Intel AVX-512 intrinsics and the solver
sustains as much as 22% of the theoretical peak flop rate of the Second
Generation Intel Xeon Phi, Knights Landing, processor
Detailed Numerical Simulation of Multi-Dimensional Plasma Assisted Combustion
Interaction between flames and plasmas are the guiding thread of this work. Nanosecond Repetitively Pulsed (NRP) discharges are non-thermal plasmas which have shown interesting features for combustion control. They can interact with flames not only through heat, but also chemically by producing active species. In this work, fully-coupled plasma assited combustion simulations are targeted. To achieve this goal, plasma discharge capabilities are built in the low temperature plasma code, AVIP. The corresponding numerical methods, as well as validation cases regarding each set of equations, are first presented. To simulate plasma discharges, the coupled drift-diffusion equations and the Poisson equation are considered. AVIP is coupled to the AVBP code which solves the reactive Navier-Stokes equations to describe combustion phenomena. In a second part, we start by constructing and validating a fully-detailed chemistry for methane-air mixtures in zero-dimensional reactors before reducing it for multi dimensional simulations. The multi-dimensional streamer simulation capabilities of the code are then assessed using simple chemistries. All the validated parts of the code come together in a fully detailed simulation of ignition using NRP discharges. We finish by discussing phenomenological models built upon the knowledge that we gained from fully-detailed simulations. In a last part, finally, attempt to solve the Poisson and generalized Poisson equations using neural networks, which have a potential for speedup compared to classical linear solvers, is carried out
Simulating radiative shocks in nozzle shock tubes
We use the recently developed Center for Radiative Shock Hydrodynamics
(CRASH) code to numerically simulate laser-driven radiative shock experiments.
These shocks are launched by an ablated beryllium disk and are driven down
xenon-filled plastic tubes. The simulations are initialized by the
two-dimensional version of the Lagrangian Hyades code which is used to evaluate
the laser energy deposition during the first 1.1ns. The later times are
calculated with the CRASH code. This code solves for the multi-material
hydrodynamics with separate electron and ion temperatures on an Eulerian
block-adaptive-mesh and includes a multi-group flux-limited radiation diffusion
and electron thermal heat conduction. The goal of the present paper is to
demonstrate the capability to simulate radiative shocks of essentially
three-dimensional experimental configurations, such as circular and elliptical
nozzles. We show that the compound shock structure of the primary and wall
shock is captured and verify that the shock properties are consistent with
order-of-magnitude estimates. The produced synthetic radiographs can be used
for comparison with future nozzle experiments at high-energy-density laser
facilities.Comment: submitted to High Energy Density Physic
- …