Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two-dimensional electron gas”

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 145 records · Page 8

X-point and divertor filament dynamics from gas puff imaging on TCV

Abstract A new gas puff imaging diagnostic has been installed on the TCV tokamak, providing two-dimensional insights into scrape-off-layer (SOL) turbulence dynamics above, at and below the magnetic X-point. A detailed study in L-mode, attached, lower single-null discharges shows that statistical properties have little poloidal variations, while vast differences are present in the 2D behaviour of intermittent filaments. Strongly elongated filaments, just above the X-point and in the divertor far-SOL, show a good consistency in shape and dynamics with field-line tracing from filaments at the outboard midplane, highlighting their connection. In the near-SOL of the outer divertor leg, short-lived, high frequency and more circular (diameter ∼15 sound Larmor radii) filaments are observed. These divertor-localised filaments appear born radially at the position of maximum density and display a radially outward motion with velocity ≈400 m s −1 that is comparable to radial velocities of upstream-connected filaments. Conversely, in these discharges ( B → × ∇ B pointing away from the divertor), these divertor filaments’ poloidal velocities differ strongly from those of upstream-connected filaments. The importance of divertor-localised filaments upon radial transport and profile broadening is explored using filament statistics and in situ kinetic profile measurements along the divertor leg. This suggests that these filaments contribute significantly to electron density profile broadening in the divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Engineering Nanostructured Interfaces of Hexagonal Boron Nitride-Based Materials for Enhanced Catalysis

Hexagonal boron nitrides (h-BNs) are attractive two-dimensional (2D) nanomaterials that consist of alternating B and N atoms and layered honeycomb-like structures similar to graphene. They have exhibited unique properties and promising application potentials in the field of energy storage and transformation. Recent advances in utilizing h-BN as a metal-free catalyst in the oxidative dehydrogenation of propane have triggered broad interests in exploring h-BN in catalysis. However, h-BN-based materials as robust nanocatalysts in heterogeneous catalysis are still underexplored because of the limited methodologies capable of affording h-BN with controllable crystallinity, abundant porosity, high purity, and defect engineering, which played important roles in tuning their catalytic performance. Here, in this account, our recent progress in addressing the above issues will be highlighted, including the synthesis of high-quality h-BN-based nanomaterials via both bottom-up and top-down pathways and their catalytic utilization as metal-free catalysts or as supports to tune the interfacial electronic properties on the metal nanoparticles (NPs). First, we will focus on the large-scale fabrication of h-BN nanosheets (h-BNNSs) with high crystallinity, improved surface area, satisfactory purity, and tunable defects. h-BN derived from the traditional approaches using boron trioxide and urea as the starting materials generally contains carbon/oxygen impurities and has low crystallinity. Several new strategies were developed to address the issues. Using bulk h-BN as the precursor via gas exfoliation in liquid nitrogen, single- or few-layered h-BNNS with abundant defects could be generated. Amorphous h-BN precursors could be converted to h-BN nanosheets with high crystallinity assisted by a magnesium metallic flux via a successive dissolution/precipitation/crystallization procedure. The as-fabricated h-BNNS featured high crystallinity and purity as well as abundant porosity. An ionothermal metathesis procedure was developed using inorganic molten salts (NaNH 2 and NaBH 4 ) as the precursors. The h-BN scaffolds could be produced on a large scale with high yield, and the as-afforded materials possessed high purity and crystallinity. Second, utilization of the as-prepared h-BN library as metal-free catalysts in dehydrogenation and hydrogenation reactions will be summarized, in which they exhibited enhanced catalytic activity over the counterparts from the previous synthesis method. Third, the interface modulation between metal NPs with the as-prepared defects’ abundant h-BN support will be highlighted. The h-BN-based strong metal–support interaction (SMSI) nanocatalysts were constructed without involving reducible metal oxides via the ionothermal procedure we developed by deploying specific inorganic metal salts, acting as robust nanocatalysts in CO oxidation. Under conditions simulated for practical exhaust systems, promising catalytic efficiency together with high thermal stability and sintering resistance was achieved. Across all of these examples, unique insights into structures, defects, and interfaces that emerge from in-depth characterization through microscopy, spectroscopy, and diffraction will be highlighted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sheath formation around a dielectric droplet in a He atmospheric pressure plasma

Interactions at the interface between atmospheric pressure plasmas and liquids are being investigated to address applications ranging from nanoparticle synthesis to decontamination and fertilizer production. Furthermore, many of these applications involve activation of droplets wherein the droplet is fully immersed in the plasma and synergistically interacts with the plasma. To better understand these interactions, two-dimensional modeling of radio frequency (RF) glow discharges at atmospheric pressure operated in He with an embedded lossy dielectric droplet (tens of microns in size) was performed. The properties of the sheath that forms around the droplet were investigated over the RF cycle. The electric field in the bulk plasma polarizes the dielectric droplet while the electron drift in the external electric field is shadowed by the droplet. The interaction between the bulk and sheath electric fields produces a maximum in E/N (electric field/gas number density) at the equator on one side of the droplet where the bulk and sheath fields are aligned in the same direction and a minimum along the opposite equator. Due to resistive heating, the electron temperature T e is maximum 45° above and below the equator of the droplet where power deposition per electron is the highest. Although the droplet is, on the average, negatively charged, the charge density on the droplet is positive on the poles and negative on the equator, as the electron motion is primarily due to diffusion at the poles but due to drift at the equator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing and controlling oxygen impurity diffusion in h -BN semi-bulk crystals

Combining its unique features of ultrawide bandgap (UWBG) and two-dimensional nature, h-BN has been explored for emerging applications such as deep ultraviolet optoelectronic devices and single photon emitters. One of the unusual applications of h-BN is for solid-state neutron detectors by utilizing the property of high thermal neutron capture cross section of B-10 as well as its UWBG properties. Although a record high detection efficiency of 59% has been attained by h-BN detectors, the understanding/minimization of defects and impurities is still needed to further advance the h-BN material and detector technologies. We report metal organic chemical vapor deposition growth and oxygen (O) impurity diffusion in thick h-BN. The diffusion coefficient (D) of O impurities has been measured via the evolution of an oxygen related emission with the etching depth, providing a value of D of ∼ 2 × 10−13 cm2/s at 1450 °C and supporting the interpretation that oxygen in h-BN is a substitutional donor. A multiple-buffer-layer approach was employed to mitigate to a certain degree the issue of oxygen diffusion from sapphire substrate during growth. It was demonstrated that the performance of h-BN neutron detectors fabricated from the wafer incorporating multiple buffer layers was significantly improved, as manifested by the enhanced thermal neutron detection efficiency. The advancement of the crystal growth technology of h-BN semi-bulk crystals creates applications in optoelectronic and power electronic devices utilizing the UWBG semiconductor properties of h-BN, while high efficiency h-BN neutron detectors have the potential to supplant the traditional He-3 gas detectors in various application areas by offering the obvious advantages of UWBG semiconductor technologies.

Physics↗

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗

Plasma assisted NH 3 /H 2 /air ignition in nanosecond discharges with non-equilibrium energy transfer

Ammonia (NH 3 ), with its high energy density and easiness to store and transport as a hydrogen carrier, has become a promising alternative green fuel. However, its adoption in power generation is hindered by challenges such as low burning velocity, slow low-temperature oxidation, high NO x emissions, and ignition difficulty. Here, this work computationally investigates the effects of non-equilibrium energy transfer by nanosecond discharges on NH 3 ignition and flame propagation in an NH 3 /H 2 /air flow at 700 K and 1 atm. The simulation results demonstrate that NH 3 /air mixtures require a large ignition energy due to their large critical ignition radius. It is shown that adding 30 % hydrogen significantly reduces the critical ignition radius and minimum ignition energy. Two-dimensional modeling further shows a non-monotonic dependence of ignition kernel volume on the applied voltage and reduced electric field. The optimum ignition enhancement occurs at 200 Td where the generation of electronically excited species and radicals including N 2 (B), O( 1 D) and OH becomes most efficient. Higher voltages divert electron energy toward ionization, which makes it less effective for NH 3 ignition. The study also identifies an optimal electrode gap size for a given pulse energy. Smaller gap sizes increase deposited energy density, raising temperature and radical concentrations. However, excessive reduction of the gap distance reduces flame propagation speed due to the flame stretch effect in rich mixtures with the effective Lewis number greater than unity. A nonlinear relationship between pulse repetition frequency and ignition kernel volume is observed in a nanosecond pulsed high frequency discharge (NPHFD). An optimal frequency range of 200 kHz to 2 MHz is found when two pulses are used. In addition, an optimal number of pulses exists for each pulse repetition frequency, with higher frequencies requiring more pulses to maximize the overlap region. These findings provide critical insights on developing controlled plasma discharge techniques for efficient NH 3 ignition in reactive flows within internal combustion engines and gas turbines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A dual dynamic shutter system for accelerating ion irradiation sample throughput via lateral gas implantation gradients

Ion irradiation for material performance testing is limited due to its serial nature, which allows for only one value of the implantation (appm) versus dose (dpa) parameter space to be explored for each ion and experiment at a time. While ion irradiation can accelerate the process by up to three orders of magnitude compared to neutron irradiation experiments, the sample throughput for ion irradiation remains relatively low. To address these limitations, a novel capability has been developed at the Michigan Ion Beam Laboratory (MIBL), enabling for the creation of single- and two-dimensional lateral ion implantation gradients using recently installed motorized-controlled ion-beam shutters. This advancement can generate a wide scope of the two-dimensional (H+, He2+) implantation parameter space within a single sample. Integration of this new capability now allows for dual- and triple-ion beam experiments to be performed with full user control over not only the ion implantation depth, but also laterally across the sample by imposing ion implantation concentration gradients, thus providing researchers with a high-throughput means for material testing under various irradiation conditions. Furthermore, recent improvements in MIBL's microbeam ion-beam analysis (IBA) target station now allow for probing these concentration gradients in irradiated alloys with exceptional spatial resolution, down to 10 µm. These two approaches promise to significantly improve ion irradiation capabilities and increase the sample throughput by several orders of magnitude. The application of the shutter technique plus the subsequent microbeam characterization of the imposed implantation gradients are showcased by two proof-of-principle ion-irradiated experiments, one performed on single-crystal Si and the other on the fusion-candidate alloy F82H-IEA. These advancements mark a substantial leap in ion-beam technology, offering researchers a robust, high-throughput method to efficiently investigate candidate alloys with high technological readiness for both advanced fission and fusion reactor applications, in a time- and cost-effective manner.

36 - MATERIALS SCIENCE↗

Visualizing Size-Dependent Dynamics of CeO 2-δ {100}-Supported CoO x Nanoparticles Under CO 2 Hydrogenation Conditions

Carbon dioxide is a major greenhouse gas. In order to optimize processes focused on its chemical valorization, one needs detailed information about the effects of CO 2 and/or CO 2 /H 2 mixtures on the structure and morphology of metal/oxide catalysts. Here, in this study, the evolution of a catalyst with cobalt supported on CeO 2 -cube nanostructures under CO 2 hydrogenation conditions was investigated by using a set of in situ characterization techniques (X-ray absorption fine structure, X-ray diffraction, diffuse reflectance infrared Fourier transform spectroscopy, and environmental transmission electron microscopy (TEM)). The {100} facets of the ceria support displayed an unexpectedly high stability due to strong interactions with the aggregates of cobalt oxide. A significant influence of interfacial bonding between CoO x and CeO 2-δ {100} is evident through a clear preference in the orientation of CoO x nanoparticles (NPs) with respect to the substrate. For initially reduced Co/CeO 2 -cube nanostructures, a kinetically controlled oxidation of cobalt upon the introduction of CO 2 was observed during the early stages of CO 2 hydrogenation. Environmental TEM revealed the size-dependent morphological behavior of cobalt oxide NPs due to strong interactions with the CeO 2 {100} surface. When the environment was switched from H 2 to a mixture of H 2 and CO 2 at 250 °C, small CoO x NPs (in the largest dimension < 2.5 nm) rapidly transformed from a pyramidal three-dimensional (3D) form to a planar, monatomic layer attached to the concurrently oxidized CeO 2-δ {100} surface. This maximizes the number of sites available for the binding of CO 2 or reaction intermediates. The shape transformation reflected the oxophilic character of cobalt and strong metal–support interactions. The removal of CO 2 from the gas phase led to a reduction of the cobalt oxide NPs by hydrogen and a reversible two-dimensional → 3D transformation. In contrast, no significant morphological changes, apart from further oxidation, were observed for big CoO x NPs (in the largest dimension > 3 nm). These trends are not seen for nanoparticles of noble metals. The observed morphological and structural changes in the small CoO x NPs affected the stability of reaction intermediates and modified the selectivity of the CoO x /CeO 2 catalyst system for methane production.

36 MATERIALS SCIENCE↗

Range-separated hybrid functionals for mixed dimensional heterojunctions: Application to phthalocyanines/MoS 2

We analyze the electronic structure and level alignment of transition-metal phthalocyanine (MPc) molecules adsorbed on two-dimensional MoS 2 employing density functional theory (DFT) calculations. We develop a procedure for multi-objective optimal tuning of parameters of range-separated hybrid functionals in these mixed-dimensional systems. Using this procedure, which leads to the asymptotically-correct exchange-correlation potential between molecule and two-dimensional material, we obtain electronic structures consistent with experimental photoemission results for both energy level alignment and electronic bandgaps, representing a significant advance compared to standard DFT methods. We elucidate the MoS 2 valence resonance with the transition-metal phthalocyanine non-frontier 3d orbitals and its dependence on the transition metal atomic number. Based on our calculations, we derive parameter-free, model self-energy corrections that quantitatively accounts for the effects of the heterogeneous dielectric environment on the electronic structure of these mixed-dimensional heterojunctions.

2D materials↗

Physics insights from a large-scale 2D UEDGE simulation database for detachment control in KSTAR

A large-scale database of two-dimensional UEDGE simulations has been developed to study detachment physics in KSTAR and to support surrogate models for control applications. Nearly 70 000 steady-state solutions were generated, systematically scanning upstream density, input power, plasma current, impurity fraction, and anomalous transport coefficients, with magnetic and electric drifts across the magnetic field included. The database identifies robust detachment indicators, with strike-point electron temperature at detachment onset consistently T e,target ~ 3-4 eV, largely insensitive to upstream conditions. Scaling relations reveal weaker impurity sensitivity than one-dimensional models and show that heat flux widths follow Eich’s scaling only for uniform, low D and χ. Distinctive in–out divertor asymmetries are observed in KSTAR, differing qualitatively from DIII-D. Complementary time-dependent simulations quantify plasma response to gas puffing, with delays of 5 - 15 ms at the outer strike point and ∼ 40 ms for the low-magnetic-field-side radiation front. These dynamics are well captured by first-order-plus-dead-time models and are consistent with experimentally observed detachment-control behavior in KSTAR (Gupta et al 2025 Plasma Phys. Control. Fusion (submitted)).

Physics - Plasma physics↗

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗

Two-Dimensional Thomson Scattering in Laser-Produced Plasmas

We present two-dimensional (2D) optical Thomson scattering measurements of electron density and temperature in laser-produced plasmas. The novel instrument directly measures ne(x,y) and Te(x,y) in two dimensions over large spatial regions (cm2) with sub-mm spatial resolution, by automatically translating the scattering volume while the plasma is produced repeatedly by irradiating a solid target with a high-repetition-rate laser beam (10 J, ∼1012 W/cm2, 1 Hz). In this paper, we describe the design and motorized auto-alignment of the instrument and the computerized algorithm that autonomously fits the spectral distribution function to the tens-of-thousands of measured scattering spectra, and captures the transition from the collective to the non-collective regime with distance from the target. As an example, we present the first 2D scattering measurements in laser-driven shock waves in ambient nitrogen gas at a pressure of 0.13 mbar.

Zhang, Haiping (ORCID:0009000228133852)↗

A two component model of the diurnal variations in the thermospheric composition

A self-consistent, two-dimensional, two-component diffusion model is presented to describe the diurnal component of thermosphere dynamics. In the regions where oxygen and helium are the minor constituents, diffusion is shown to increase the amplitudes of these constituents and to advance the times of their diurnal maxima by several hours with respect to the gas temperature and the major species. Both effects are in basic agreement with satellite observations. The diffusion process is relatively insignificant for temperature, total mass density, and wind fields. The magnitude of the temperature-density phase delay is discussed with respect to its dependence on the electron density distribution.

Mayr, H. G.↗

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗

Electron Paramagnetic Resonance Imaging of the Spatial Distribution of Free Radicals in PMR-15 Polyimide Resins

Prior studies have shown that free radicals generated by heating polyimides above 300 C are stable at room temperature and are involved in thermo-oxidative degradation in the presence of oxygen gas. Electron Paramagnetic Resonance Imaging (EPRI) is a technique to determine the spatial distribution of free radicals. X-band (9.5 GHz) EPR images of PMR-15 polyimide were obtained with a spatial resolution of about 0.18 mm along a 2 mm dimension of the sample. In a polyimide sample that was not thermocycled, the radical distribution was uniform along the 2 mm dimension of the sample. For a polyimide sample that was exposed to thermocycling in air for 300 one-hour cycles at 335 C, one-dimensional EPRI showed a higher concentration of free radicals in the surface layers than in the bulk sample. A spectral-spatial two-dimensional image showed that the EPR lineshape of the surface layer remained the same as that of the bulk. These EPRI results suggest that the thermo-oxidative degradation of PMR-15 resin involves free radicals present in the oxygen-rich surface layer.

Ahn, Myong K.↗

Electron Paramagnetic Resonance Imaging of the Spatial Distribution of Free Radicals in PMR-15 Polyimide Resins

Prior studies have shown that free radicals generated by heating polyimides above 300 C are stable at room temperature and are involved in thermo-oxidative degradation in the presence of oxygen gas. Electron paramagnetic resonance imaging (EPRI) is a technique to determine the spatial distribution of free radicals. X-band (9.5 GHz) EPR images of PMR-15 polyimide were obtained with a spatial resolution of approximately 0.18 mm along a 2-mm dimension of the sample. In a polyimide sample that was not thermocycled, the radical distribution was uniform along the 2-mm dimension of the sample. For a polyimide sample that was exposed to thermocycling in air for 300 1-h cycles at 335 C, one-dimensional EPRI showed a higher concentration of free radicals in the surface layers than in the bulk sample. A spectral-spatial two-dimensional image showed that the EPR lineshape of the surface layer remained the same as that of the bulk. These EPRI results suggest that the thermo-oxidative degradation of PMR-15 resin involves free radicals present in the oxygen-rich surface layer.

Ahn, Myong K.↗

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

Wear Mechanisms in Electron Sources for Ion Propulsion, 1: Neutralizer Hollow Cathode

Upon the completion of two long-duration life tests of a 30-cm ion engine, the orifice channel of the neutralizer hollow cathode was eroded away to as much as twice its original diameter. Whereas the neutralizer cathode orifice opened significantly, no noticeable erosion of the discharge cathode orifice was observed. Noquantitative explanation of these erosion trends has been established since the completion of the two life tests. A two-dimensional model of the partially ionized gas inside these devices has been developed and applied to the neutralizer hollow cathode. The numerical simulations show that the main mechanism responsible for the channel erosion is sputtering by Xe+. These ions are accelerated by the sheath along the channel and bombard the surface with kinetic energy/charge of about 17 V at the beginning of cathode life. The density of the ions inside the neutralizer orifice is computed to be as high as 2.1 x 10(sup 22) m(sup -3). Because of the 3.5-times larger diameter of the discharge cathode orifice, the ion density inside the orifice is more than 40 times lower and the sheath drop 7 V lower compared with the values in the neutralizer. At these conditions, Xe+ can cause no significant sputtering of the surface.

ionized gas↗