Source author record

A. S. Richardson

A. S. Richardson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph physics.comp-ph physics.chem-ph physics.flu-dyn

Catalog footprint

What is connected

5works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Simulating a pulsed power-driven plasma with ideal MHD

We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the MHD domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This, much more affordable ideal MHD representation, is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.

preprint2021arXiv

The rigid-beam model for simulating plasmas generated by intense electron beams

We introduce a simplified model of the electron-beam/plasma system to model the electrical breakdown caused by the inductive electric field created by a rapidly rising electron beam current. The rigid-beam model is a reduction to the problem geometry to cylindrical coordinated and simplifications to Maxwell's equations that are driven by a prescribed electron beam current density. The model is very convenient for comparing various reductions of the plasma dynamics and air chemistry equation while maintaining a good approximation to the overall magnitude of the beam-created electric field. The usefulness of this model is demonstrated by comparing results for two different fluid reductions of the plasma dynamics: one where the collision rates are computed from the local reduced electric field (E/p) and another where the collision rates are determined from the mean energy per particle. We find that the two methods give similar results at higher pressures where the energy relation rate is large but differs significantly at lower pressures where the characteristic inelastic energy loss time scale is comparable to or greater than the rise time of the electron beam current.

preprint2020arXiv

TurboPy: A Lightweight Python Framework for Computational Physics

Computational physics problems often have a common set of aspects to them that any particular numerical code will have to address. Because these aspects are common to many problems, having a framework already designed and ready to use will not only speed the development of new codes, but also enhance compatibility between codes. Some of the most common aspects of computational physics problems are: a grid, a clock which tracks the flow of the simulation, and a set of models describing the dynamics of various quantities on the grid. Having a framework that could deal with these basic aspects of the simulation in a common way could provide great value to computational scientists by solving various numerical and class design issues that routinely arise. This paper describes the newly developed computational framework that we have built for rapidly prototyping new physics codes. This framework, called turboPy, is a lightweight physics modeling framework based on the design of the particle-in-cell code turboWAVE. It implements a class (called Simulation) which drives the simulation and manages communication between physics modules, a class (called PhysicsModule) which handles the details of the dynamics of the various parts of the problem, and some additional classes such as a Grid class and a Diagnostic class to handle various ancillary issues that commonly arise.

preprint2011arXiv

Quasi-separatrix layers and three-dimensional reconnection diagnostics for line-tied tearing modes

In three-dimensional magnetic configurations for a plasma in which no closed field line or magnetic null exists, no magnetic reconnection can occur, by the strictest definition of reconnection. A finitely long pinch with line-tied boundary conditions, in which all the magnetic field lines start at one end of the system and proceed to the opposite end, is an example of such a system. Nevertheless, for a long system of this type, the physical behavior in resistive magnetohydrodynamics (MHD) essentially involves reconnection. This has been explained in terms comparing the geometric and tearing widths [1, 2]. The concept of a quasi-separatrix layer[3, 4] was developed for such systems. In this paper we study a model for a line-tied system in which the corresponding periodic system has an unstable tearing mode. We analyze this system in terms of two magnetic field line diagnostics, the squashing factor[3-5] and the electrostatic potential difference used in kinematic reconnection studies[6, 7]. We discuss the physical and geometric significance of these two diagnostics and compare them in the context of discerning tearing-like behavior in line-tied modes. [1] G. L. Delzanno and J. M. Finn. Physics of Plasmas, 15(3):032904, 2008. [2] Y.-M. Huang and E. G. Zweibel. Physics of Plasmas, 16(4):042102, 2009. [3] E. R. Priest and P. Démoulin. J. Geophys. Res., 100(A12):23443-23463, 1995. [4] P. Démoulin, J. C. Henoux, E. R. Priest, and C. H. Mandrini. Astron. Astrophys., 308:643-655, Apr. 1996. [5] V. S. Titov and G. Hornig. Advances in Space Research, 29(7):1087-1092, 2002. [6] Y. Lau and J. M. Finn. The Astrophysical Journal, 350:672-691, Feb. 1990. [7] Y. Lau and J. M. Finn. The Astrophysical Journal, 366:577-591, 1991.

preprint2010arXiv

The lower hybrid wave cutoff: A case study in eikonal methods

Eikonal, or ray tracing, methods are commonly used to estimate the propagation of radio frequency fields in plasmas. While the information gained from the rays is quite useful, an approximate solution for the fields would also be valuable, e.g., for comparison to full wave simulations. Such approximations are often difficult to perform numerically because of the special care which must be taken to correctly reconstruct the fields near reflection and focusing caustics. In this paper, we compare the standard eikonal method for approximating fields to a method based on the dynamics of wave packets. We compare the approximations resulting from these two methods to the analytical solution for a lower hybrid wave reflecting from a cutoff. The algorithm based on wave packets has the advantage that it can correctly deal with caustics, without any special treatment.