Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “image coding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogeneities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogenities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ICME and In-Situ Process Monitoring for Rapid Qualification of Components Made by Laser-based Powder Bed Additive Manufacturing Processes for Nuclear Structural Applications

Additive manufacturing (AM) through selective powder bed melting of successive layers is being considered as disruptive technology for rapid production of low-cost nuclear reactor internal components with complex geometries. However, there is a potential for non-uniform distribution of physical features, such as porosity or microstructural differences, due to variations of temperature across a component build. These heterogenities make the deployment of traditional material qualification and non-destructive evaluation of AM components difficult. This research project explored the feasibility of using in-situ process monitoring methods and integrated computational materials engineering (ICME) principles as an alternate qualification methodology and approach. This project included six individual tasks: (i) design of artifacts relevant to the nuclear power industry, (ii) evaluation of laser processing and in-situ measurements, (iii) computational modeling, (iv) ex-situ microstructural characterization, (v) evaluation of scaling the methodology for large-scale structures, and (vi) development of a data package to codes and standards organizations. As a part of this research, laser powder bed fusion of metals (316L, Alloy 718, and Ti6Al4V) was explored. This document provides an overview of the research performed over three years and serves as the final U.S. Department of Energy report for this project. In the third year, the project focused on the following activities: (a) in-situ infrared and optical imaging of each layer during processing of representative stainless steel geometries; (b) development of methodologies to analyze the data; and (c) application of ICME methodologies and other experiments to estimate the relevance of defects and microstructure to tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Design of source coders and joint source/channel coders for noisy channels

A theory behind a proposed joint source/channel coding approach is developed and a variable rate design approach which provides substantial improvement over current joint source/channel coder designs is obtained. The Rice algorithm as applied to the output of the Gamma Ray Detector of the Mars Orbiter is evaluated. An alternative algorithm is obtained which outperforms the Rice both in terms of data compression and noisy channel performance. A high-fidelity low-rate image compression algorithm is developed which provides almost distortionless compression of high resolution images.

Sayood, Khalid↗

AXAF Detector Backgrounds Produced By Cosmic Ray Protons

One of the science instruments on the Advanced X-ray Astrophysics Facility (AXAF), planned for launch in 1998 into a highly elliptical (10,000 km x 140,000 km) orbit, is a microchannel plate High Resolution Camera (HRC). This detector is designed to provide imaging and spectroscopic observations of x-rays emitted by stellar sources in the 0.1 to 10 keV energy range. Described here are analyses made to determine the expected time-dependent detector background from prompt and delayed (activation) radiation initiated by galactic cosmic-ray (GCR) proton interactions in the spacecraft and payload. Numerical simulations were made using the coupled set of Monte Carlo radiation transport codes, analysis software, and data bases shown. The major codes are HETC for nucleon-meson transport, EGS for simulating electromagnetic cascades, and MORSE for low-energy (less than 15 MeV) neutron transport. The simulation follows the transport history of photons in the energy range from - 100 GeV down to approx. 0.1 keV due to gamma-ray sources from neutral pion decay, high-energy (spallation) collisions, and low-energy neutron inelastic scattering and capture reactions. Also included is radioisotope production and the tracking of gamma-rays, electrons, and positrons from induced radioactivity.

Armstrong, T. W.↗

Demonstration and Automation of Reflected Target Optical Measurement for Heliostats

Accurate optical surfaces are a primary driver of concentrated solar power plant performance. Errors in pointing and tracking mirrors, the canting of individual mirror facets, and the surface slope of the mirror itself can be caused by errors during assembly, transportation, wind loading, gravity, and many other sources. The tools that exist to measure these error sources today largely rely on fringe deflectometry (SOFAST, QDec, others), or photogrammetry with targets attached to the mirror surface. Since 2022, NREL has been developing a measurement method called the Reflected Target Non-intrusive Assessment (ReTNA) system. This system differs from most established methods in that we perform deflectometry with a pattern of coded targets, identified in space with photogrammetry. Reflected target systems have several advantages over traditional fringe deflectometry systems. Firstly, they can be operated in bright or ambient lighting, a challenge for fringe systems that use a projector and screen. Reflected target systems also can use a much lighter and less expensive target than projector-based systems. Lastly, 2D slope measurement can be solved from a single image, which leads to several advantages for accommodating faster measurements and smaller sized targets. These advantages make ReTNA particularly well-suited for applications where there are space or lighting constraints, like performing heliostat quality assurance on an assembly line. It's also useful when a lightweight, flexible system is needed, like for heliostat developers to quickly measure a new heliostat design at different orientations, to observe gravitational effects on the mirror surface shape. In the last year, significant improvements were made to this tool to make it more useful for these applications. These improvements were focused around validation of the ReTNA measurement system, and automation of the setup and measurement process. First, we present an improved ReTNA layout, for use on the heliostat assembly line. Next, we detail the various changes to the ReTNA software and computer vision methods to automate data collection in this new setup, and lessons learned from this process. The goal with this new setup is to perform a full heliostat surface characterization without removing the mirror from the assembly line. Lastly, we share results from several ReTNA validation studies undertaken over the last year. These include repeated ReTNA measurement on demonstration mirror facets, comparisons with other optical measurement tools, and some studies aimed at quantifying the uncertainty of ReTNA measurement under various constraints (mirror-target spacing, camera resolution, etc.). These results are compared with 2024 HelioCon performance targets, and our planned next steps for the ReTNA measurement system are presented.

CSP↗

The Role of Ontologies in Schema-based Program Synthesis

Program synthesis is the process of automatically deriving executable code from (non-executable) high-level specifications. It is more flexible and powerful than conventional code generation techniques that simply translate algorithmic specifications into lower-level code or only create code skeletons from structural specifications (such as UML class diagrams). Key to building a successful synthesis system is specializing to an appropriate application domain. The AUTOBAYES and AUTOFILTER systems, under development at NASA Ames, operate in the two domains of data analysis and state estimation, respectively. The central concept of both systems is the schema, a representation of reusable computational knowledge. This can take various forms, including high-level algorithm templates, code optimizations, datatype refinements, or architectural information. A schema also contains applicability conditions that are used to determine when it can be applied safely. These conditions can refer to the initial specification, to intermediate results, or to elements of the partially-instantiated code. Schema-based synthesis uses AI technology to recursively apply schemas to gradually refine a specification into executable code. This process proceeds in two main phases. A front-end gradually transforms the problem specification into a program represented in an abstract intermediate code. A backend then compiles this further down into a concrete target programming language of choice. A core engine applies schemas on the initial problem specification, then uses the output of those schemas as the input for other schemas, until the full implementation is generated. Since there might be different schemas that implement different solutions to the same problem this process can generate an entire solution tree. AUTOBAYES and AUTOFILTER have reached the level of maturity where they enable users to solve interesting application problems, e.g., the analysis of Hubble Space Telescope images. They are large (in total around 100kLoC Prolog), knowledge intensive systems that employ complex symbolic reasoning to generate a wide range of non-trivial programs for complex application do- mains. Their schemas can have complex interactions, which make it hard to change them in isolation or even understand what an existing schema actually does. Adding more capabilities by increasing the number of schemas will only worsen this situation, ultimately leading to the entropy death of the synthesis system. The root came of this problem is that the domain knowledge is scattered throughout the entire system and only represented implicitly in the schema implementations. In our current work, we are addressing this problem by making explicit the knowledge from Merent parts of the synthesis system. Here; we discuss how Gruber's definition of an ontology as an explicit specification of a conceptualization matches our efforts in identifying and explicating the domain-specific concepts. We outline the dual role ontologies play in schema-based synthesis and argue that they address different audiences and serve different purposes. Their first role is descriptive: they serve as explicit documentation, and help to understand the internal structure of the system. Their second role is prescriptive: they provide the formal basis against which the other parts of the system (e.g., schemas) can be checked. Their final role is referential: ontologies also provide semantically meaningful "hooks" which allow schemas and tools to access the internal state of the program derivation process (e.g., fragments of the generated code) in domain-specific rather than language-specific terms, and thus to modify it in a controlled fashion. For discussion purposes we use AUTOLINEAR, a small synthesis system we are currently experimenting with, which can generate code for solving a system of linear equations, Az = b.

Bures, Tomas↗

Imaging of the gravitational lens system PG 1115+080 with the Hubble Space Telescope

This paper is the first of a series presenting observations of gravitational lenses and lens candidates, taken with the Wide Field/Planetary Camera of the HST. We have resolved the gravitational lens system PG 1115+080 into four point sources and a red, extended object that is presumably the lens galaxy; we present accurate relative intensities, colors, and positions of the four images, and lower accuracy intensity and position of the lens galaxy, all at the epoch 1991.2. Comparison with earlier data shows no compelling evidence for relative intensity variations between the QSO components having so far been observed. The new data agree with earlier conclusions that the system is rather simple, and can be produced by the single observed galaxy. The absence of asymmetry in the HST images implies that the emitting region of the quasar itself has an angular radius smaller than about 10 milliarcsec (100 pc for H0 = 50, q0 = 0.5).

Kristian, Jerome↗

Derivation of scaled surface reflectances from AVIRIS data

A method for retrieving 'scaled surface reflectances' assuming horizontal surfaces having Lambertian reflectances from spectral data collected by Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is presented here. In this method, the integrated water vapor amount on a pixel by pixel basis is derived from the 0.94 micron and 1.14 micron water vapor absorption features. The transmission spectra of H2O, CO2, O3, N2O, CO, CH4, and O2 in the 0.4-2.5 micron region are simulated. The scattering effect due to atmospheric molecules and aerosols is modeled with the 5S computer code. The AVIRIS radiances are divided by solar irradiances above the atmosphere to obtain the apparent reflectances. The scaled surface reflectances are derived from the apparent reflectances using the simulated atmospheric gaseous transmittances and the simulated molecular and aerosol scattering data. The scaled surface reflectances differ from the real surface reflectances by a multiplicative factor. In order to convert the scaled surface reflectances into real surface reflectances, the slopes and aspects of the surfaces must be known.

Gao, Bo-Cai↗

Modelling of a space-based Ronchi ruling experiment for high-precision astrometry

An end-to-end modelling program for an astrometric telescope employing a Ronchi ruling has been developed. The program models the aberrated images formed anywhere in the field-of-view. It then determines apparent centroids by simulating the motion of a Ronchi ruling across the field. Photo-electron statistics are included. A 6-term plate-constant model is used to determine the apparent motion of a target star within a reference frame as the object is re-observed through both ideal and perturbed optics. The modelling code is accurate at the submicroarcsecond level. optics. The modelling code is accurate at the submicroarcsecond level.

Shaklan, S.↗

Testing of Haar-Like Feature in Region of Interest Detection for Automated Target Recognition (ATR) System

The objectives of this project were to develop a ROI (Region of Interest) detector using Haar-like feature similar to the face detection in Intel's OpenCV library, implement it in Matlab code, and test the performance of the new ROI detector against the existing ROI detector that uses Optimal Trade-off Maximum Average Correlation Height filter (OTMACH). The ROI detector included 3 parts: 1, Automated Haar-like feature selection in finding a small set of the most relevant Haar-like features for detecting ROIs that contained a target. 2, Having the small set of Haar-like features from the last step, a neural network needed to be trained to recognize ROIs with targets by taking the Haar-like features as inputs. 3, using the trained neural network from the last step, a filtering method needed to be developed to process the neural network responses into a small set of regions of interests. This needed to be coded in Matlab. All the 3 parts needed to be coded in Matlab. The parameters in the detector needed to be trained by machine learning and tested with specific datasets. Since OpenCV library and Haar-like feature were not available in Matlab, the Haar-like feature calculation needed to be implemented in Matlab. The codes for Adaptive Boosting and max/min filters in Matlab could to be found from the Internet but needed to be integrated to serve the purpose of this project. The performance of the new detector was tested by comparing the accuracy and the speed of the new detector against the existing OTMACH detector. The speed was referred as the average speed to find the regions of interests in an image. The accuracy was measured by the number of false positives (false alarms) at the same detection rate between the two detectors.

neural network↗

Synthetic Observations of Helioseismic and Magnetic Imager

Synthetic observations of the FeI line (6173A), which corresponds to the operational wavelength of Helioseismic and Magnetic Imager (HMI) on board of the NASA Solar Dynamics Observatory (SDO). Synthetic observations represent 4D maps for four angular distances from the solar disk center (0deg), 30deg, 45deg, and 60deg at 30degrees latitude in range +/- 0.75 Angstrom from the reference line of 6173.341 Angstrom, number of points across the spectral line is 256, spatial resolution of 100km, and temporal resolution of seconds. Synthetic data set obtained using SPINOR code (Frutiger et al. 2000). Synthetic observations computed from hydrodynamic simulations, which cover 25Mm of upper layers of the convection zone and 1Mm of atmosphere, the horizontal size of the domain of 80x80 Mm. Simulation results obtained the “StellarBox” code (Wray et al., 2015, 2018) with imposed rotation rate at 30deg latitude using F-plane approximation.

Irina Kitiashvili↗

The Atmospheric Response to High Nonthermal Electron-beam Fluxes in Solar Flares. II. Hydrogen-broadening Predictions for Solar Flare Observations with the Daniel K. Inouye Solar Telescope

Redshifted components of chromospheric emission lines in the hard X-ray impulsive phase of solar flares have recently been studied through their 30 s evolution with the high resolution of the Interface Region Imaging Spectrograph. Radiative-hydrodynamic flare models show that these redshifts are generally reproduced by electron-beam-generated chromospheric condensations. The models produce large ambient electron densities, and the pressure broadening of the hydrogen Balmer series should be readily detected in observations. To accurately interpret the upcoming spectral data of flares with the DKIST, we incorporate nonideal, nonadiabatic line-broadening profiles of hydrogen into the RADYN code. These improvements allow time-dependent predictions for the extreme Balmer line wing enhancements in solar flares. We study two chromospheric condensation models, which cover a range of electron-beam fluxes (1 − 5 × 10^(11) erg s^(−1) cm^(−2)) and ambient electron densities (1 − 60 × 10^(13) cm^(−3)) in the flare chromosphere. Both models produce broadening and redshift variations within 10 s of the onset of beam heating. In the chromospheric condensations, there is enhanced spectral broadening due to large optical depths at Hα, Hβ, and Hγ, while the much lower optical depth of the Balmer series H12−H16 provides a translucent window into the smaller electron densities in the beam-heated layers below the condensation. The wavelength ranges of typical DKIST/ViSP spectra of solar flares will be sufficient to test the predictions of extreme hydrogen wing broadening and accurately constrain large densities in chromospheric condensations.

Adam F. Kowalski↗

In-flight calibration of the spectral and radiometric characteristics of AVIRIS in 1991

On 7 Mar. 1991, an in-flight calibration experiment was held at the Ivanpah Playa in southeastern California for the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) imaging spectrometer. Five AVIRIS overflights were acquired of a calibration target designated on the Ivanpah Playa surface. At the time of the overflights, the reflectance of the calibration target was measured with a field spectrometer. In addition, the atmospheric optical depths and water vapor abundance were measured from a radiometer station adjacent to the calibration target. These in-situ measurements were used to constrain the MODTRAN radiative transfer code to model the upwelling spectral radiance incident to the sensor aperture during the overflights. Analyses of this modeled radiance in conjunction with the laboratory-calibrated radiance were used to determine the spectral and radiometric calibration of AVIRIS while in flight.

Green, Robert O.↗

Measurements of Heat Transfer, Flow, and Pressures in a Simulated Turbine Blade Internal Cooling Passage

An experimental study was made to obtain quantitative information on heat transfer, flow, and pressure distribution in a branched duct test section that had several significant features of an internal cooling passage of a turbine blade. The objective of this study was to generate a set of experimental data that could be used for validation of computer codes that would be used to model internal cooling. Surface heat transfer coefficients and entrance flow conditions were measured at nominal entrance Reynolds numbers of 45,000, 335,000, and 726,000. Heat transfer data were obtained by using a steady-state technique in which an Inconel heater sheet is attached to the surface and coated with liquid crystals. Visual and quantitative flow-field data from particle image velocimetry measurements for a plane at midchannel height for a Reynolds number of 45,000 were also obtained. The flow was seeded with polystyrene particles and illuminated by a laser light sheet. Pressure distribution measurements were made both on the surface with discrete holes and in the flow field with a total pressure probe. The flow-field measurements yielded flow-field velocities at selected locations. A relatively new method, pressure sensitive paint, was also used to measure surface pressure distribution. The pressure paint data obtained at Reynolds numbers of 335,000 and 726,000 compared well with the more standard method of measuring pressures by using discrete holes.

Russell, Louis M.↗

Experimental Results from the Thermal Energy Storage-2 (TES-2) Flight Experiment

Thermal Energy Storage-2 (TES-2) is a flight experiment that flew on the Space Shuttle Endeavour (STS-72), in January 1996. TES-2 originally flew with TES-1 as part of the OAST-2 Hitchhiker payload on the Space Shuttle Columbia (STS-62) in early 1994. The two experiments, TES-1 and TES-2 were identical except for the fluoride salts to be characterized. TES-1 provided data on lithium fluoride (LiF), TES-2 provided data on a fluoride eutectic (LiF/CaF2). Each experiment was a complex autonomous payload in a Get-Away-Special payload canister. TES-1 operated flawlessly for 22 hr. Results were reported in a paper entitled, Effect of Microgravity on Materials Undergoing Melting and Freezing-The TES Experiment, by David Namkoong et al. A software failure in TES-2 caused its shutdown after 4 sec of operation. TES-1 and 2 were the first experiments in a four experiment suite designed to provide data for understanding the long duration microgravity behavior of thermal energy storage salts that undergo repeated melting and freezing. Such data have never been obtained before and have direct application for the development of space-based solar dynamic (SD) power systems. These power systems will store energy in a thermal energy salt such as lithium fluoride or a eutectic of lithium fluoride/calcium difluoride. The stored energy is extracted during the shade portion of the orbit. This enables the solar dynamic power system to provide constant electrical power over the entire orbit. Analytical computer codes were developed for predicting performance of a space-based solar dynamic power system. Experimental verification of the analytical predictions were needed prior to using the analytical results for future space power design applications. The four TES flight experiments were to be used to obtain the needed experimental data. This paper will address the flight results from the first and second experiments, TES-1 and 2, in comparison to the predicted results from the Thermal Energy Storage Simulation (TESSIM) analytical computer code. An analysis of the TES-2 data was conducted by Cleveland State University Professor, Mounir Ibrahim. TESSIM validation was based on two types of results; temperature history of various points on the containment vessel and TES material distribution within the vessel upon return from flight. The TESSIM prediction showed close comparison with the flight data. Distribution of the TES material within the vessel was obtained by a tomography imaging process. The frozen TES material was concentrated toward the colder end of the canister. The TESSIM prediction indicated a similar pattern. With agreement between TESSIM and the flight data, a computerized representation was produced to show the movement and behavior of the void during the entire melting and freezing cycles.

Tolbert, Carol↗

Physical Properties of Sub-Galactic Clumps at 0.5 ≤ Z ≤ 1.5 in the UVUDF

We present an investigation of clumpy galaxies in the Hubble Ultra Deep Field at 0.5 equal to or less than z equal to or less than 1.5 in the rest-frame far-ultraviolet (FUV) using Hubble Space Telescope Wide Field Camera 3 broadband imaging in F225W, F275W, and F336W. An analysis of 1404 galaxies yields 209 galaxies that host 403 kpc scale clumps. These host galaxies appear to be typical star-forming galaxies, with an average of 2 clumps per galaxy and reaching a maximum of 8 clumps. We measure the photometry of the clumps and determine the mass, age, and star formation rates (SFR) using the spectral energy distribution fitting code FAST (Fitting and Assessment of Synthetic Templates). We find that clumps make an average contribution of 19% to the total rest-frame FUV flux of their host galaxy. Individually, clumps contribute a median of 5% to the host galaxy SFR and an average of approximately 4% to the host galaxy mass, with total clump contributions to the host galaxy stellar mass ranging widely from lower than 1% up to 93%. Clumps in the outskirts of galaxies are typically younger, with higher SFRs, than clumps in the inner regions. The results are consistent with clump migration theories in which clumps form through violent gravitational instabilities in gas-rich turbulent disks, eventually migrate toward the center of the galaxies, and coalesce into the bulge.

galaxies: structure↗

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗