Search NASA⌕ Search

SEARCH · Search NASA

Results for “NEUTRAL ELEMENT”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Extreme ultraviolet spectroscopy and photometry of VV Puppis during a high accretion state

We determine the physical properties of the accretion region of the AM Her-type binary VV Puppis using extreme ultraviolet (EUV) medium-resolution spectroscopy and photometry obtained with the Extreme Ultraviolet Explorer (EUVE) observatory. The EUV continuum from VV Pup was detected in the wavelength range from 75 to 135 A and was simultaneously recorded with the Deep Survey/Spectrometer (DS/S) imaging telescope, allowing for the extraction of an accurate light curve. VV Pup appeared to have entered a high-accretion state just prior to the pointed EUVE observations. We use the EUV light curve to infer the diameter of the accretion region (d = 220 km) assuming a hemispherical geometry and a radius of 9000 km for the white dwarf. We perform a model atmosphere analysis and, based on the light curve properties and assuming a distance of 145 pc, we derive an effective temperature of the accretion region in the range 270,000 is less than T(sub eff) is less than 360,000 K and a neutral hydrogen column density in the local interstellar medium of n(sub H) = 1.9 - 3.7 x 10(exp 19)/sq cm. The total EUV/soft X-ray energy radiated by the accretion region is approximately 3.5 x 10(exp 32) ergs/s. Our results provide a first verification of past suggestions that deep heating of the white dwarf surface produces the soft X-ray flux from the polars. We present a possible detection of O VI absortion features, and we suggest that extensive EUVE observations targeting high-accretion events may result in oxygen and heavier element abundance determination in the accretion region.

Vennes, Stephane↗
DOE OSTI2020

Laser-Induced Plasma as a Method for the Metallic Materials Hardness Estimation: An Alternative Approach

The application of an alternative laser-induced breakdown spectroscopy (LIBS) method based on transversely excited atmospheric (TEA) CO2 laser was investigated for the first time for estimating the hardness of metallic materials. The human eye-safe (TEA) CO2 laser, operating at 10.6 µm, was used for plasma generation. The LIBS spectra were recorded by employing a cost-effective CCD camera for the time-integrated and spatial resolved measurements. The cast iron and aluminum alloys samples with different hardness have been tested. The ratio between the magnesium ionic and neutral lines in LIBS spectra was applied for estimating the material hardness. In addition, the hardness of all samples included in this study was determined using the conventional method for material hardness determination, i.e. Vickers hardness test. The linear dependence of magnesium lines intensity ratio on the sample’s hardness was obtained for both kinds of materials. Profilometric measurements were used to verify that the newly-introduced method is almost nondestructive for the investigated metals. The results obtained in the present work confirmed the potential of using this original LIBS system not only for the fast elemental analysis but also for the direct estimation of the hardness of metals and alloys. © 2020, Springer Science+Business Media, LLC, part of Springer Nature.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Deriving the intermediate spectra and photocycle kinetics from time-resolved difference spectra of bacteriorhodopsin. The simpler case of the recombinant D96N protein

The bacteriorhodopsin photocycle contains more than five spectrally distinct intermediates, and the complexity of their interconversions has precluded a rigorous solution of the kinetics. A representation of the photocycle of mutated D96N bacteriorhodopsin near neutral pH was given earlier (Varo, G., and J. K. Lanyi. 1991. Biochemistry. 30:5008-5015) as BRhv-->K<==>L<==>M1-->M2--> BR. Here we have reduced a set of time-resolved difference spectra for this simpler system to three base spectra, each assumed to consist of an unknown mixture of the pure K, L, and M difference spectra represented by a 3 x 3 matrix of concentration values between 0 and 1. After generating all allowed sets of spectra for K, L, and M (i.e., M1 + M2) at a 1:50 resolution of the matrix elements, invalid solutions were eliminated progressively in a search based on what is expected, empirically and from the theory of polyene excited states, for rhodopsin spectra. Significantly, the average matrix values changed little after the first and simplest of the search criteria that disallowed negative absorptions and more than one maximum for the M intermediate. We conclude from the statistics that during the search the solutions strongly converged into a narrow region of the multidimensional space of the concentration matrix. The data at three temperatures between 5 and 25 degrees C yielded a single set of spectra for K, L, and M; their fits are consistent with the earlier derived photocycle model for the D96N protein.

NASA Program Exobiology↗
NASA NTRS2014

Artificial Intelligence Based Control Power Optimization on Tailless Aircraft

Traditional methods of control allocation optimization have shown difficulties in exploiting the full potential of controlling large arrays of control devices on innovative air vehicles. Artificial neutral networks are inspired by biological nervous systems and neurocomputing has successfully been applied to a variety of complex optimization problems. This project investigates the potential of applying neurocomputing to the control allocation optimization problem of Hybrid Wing Body (HWB) aircraft concepts to minimize control power, hinge moments, and actuator forces, while keeping system weights within acceptable limits. The main objective of this project is to develop a proof-of-concept process suitable to demonstrate the potential of using neurocomputing for optimizing actuation power for aircraft featuring multiple independently actuated control surfaces. A Nastran aeroservoelastic finite element model is used to generate a learning database of hinge moment and actuation power characteristics for an array of flight conditions and control surface deflections. An artificial neural network incorporating a genetic algorithm then uses this training data to perform control allocation optimization for the investigated aircraft configuration. The phase I project showed that optimization results for the sum of required hinge moments are improved by more than 12% over the best Nastran solution by using the neural network optimization process.

Gern, Frank↗
DOE OSTI2026

Optimal Hairpin Winding Configuration for EV Traction Motors to Enhance Active CMV Cancellation Capability of NPL.X Inverter

Active common-mode voltage (CMV) cancellation in neutral-point-less (NPL.X) inverter is an effective approach for reducing common-mode electromagnetic interference and alleviating voltage stress associated with high-speed switching in EV traction drives. However, the existing hairpin winding configuration often exhibits impedance imbalance due to slot leakage inductance and manufacturing constraints, limiting the effectiveness of active CMV cancellation. This paper investigates symmetric hairpin winding configurations for a dual-three-phase EV traction machine to improve impedance symmetry. Considering the manufacturing constraints of six-layer hairpin windings, four feasible symmetric winding configurations are analyzed. Two-dimensional finite-element analysis is performed to extract the impedance parameters of each configuration. The extracted impedance parameters are used to quantify the impedance mismatch between complementary phase pairs and evaluate its influence on active CMV cancellation. The results show that the fully interleaved winding configuration achieves perfectly balanced phase impedances with 0 % mismatch, whereas the mismatch in the other configurations ranges from approximately 2.4 % to 0.8 %. Consequently, the resulting total CMV overshoot is significantly suppressed to near zero voltage under high-voltage and high switching frequency operating conditions, demonstrating the importance of winding symmetry for maximizing the active CMV cancellation capability of the NPL.X inverter.

Lee, Kangbeen [Purdue Univ., West Lafayette, IN (U↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
NASA NTRS2015

Pressure Profiles in a Loop Heat Pipe Under Gravity Influence

During the operation of a loop heat pipe (LHP), the viscous flow induces pressure drops in various elements of the loop. The total pressure drop is equal to the sum of pressure drops in vapor grooves, vapor line, condenser, liquid line and primary wick, and is sustained by menisci at liquid and vapor interfaces on the outer surface of the primary wick in the evaporator. The menisci will curve naturally so that the resulting capillary pressure matches the total pressure drop. In ground testing, an additional gravitational pressure head may be present and must be included in the total pressure drop when LHP components are placed in a non-planar configuration. Under gravity-neutral and anti-gravity conditions, the fluid circulation in the LHP is driven solely by the capillary force. With gravity assist, however, the flow circulation can be driven by the combination of capillary and gravitational forces, or by the gravitational force alone. For a gravity-assist LHP at a given elevation between the horizontal condenser and evaporator, there exists a threshold heat load below which the LHP operation is gravity driven and above which the LHP operation is capillary force and gravity co-driven. The gravitational pressure head can have profound effects on the LHP operation, and such effects depend on the elevation, evaporator heat load, and condenser sink temperature. This paper presents a theoretical study on LHP operations under gravity neutral, anti-gravity, and gravity-assist modes using pressure diagrams to help understand the underlying physical processes. Effects of the condenser configuration on the gravitational pressure head and LHP operation are also discussed.

0000↗
NASA NTRS2015

Pressure Profiles in a Loop Heat Pipe under Gravity Influence

During the operation of a loop heat pipe (LHP), the viscous flow induces pressure drops in various elements of the loop. The total pressure drop is equal to the sum of pressure drops in vapor grooves, vapor line, condenser, liquid line and primary wick, and is sustained by menisci at liquid and vapor interfaces on the outer surface of the primary wick in the evaporator. The menisci will curve naturally so that the resulting capillary pressure matches the total pressure drop. In ground testing, an additional gravitational pressure head may be present and must be included in the total pressure drop when LHP components are placed in a non-planar configuration. Under gravity-neutral and anti-gravity conditions, the fluid circulation in the LHP is driven solely by the capillary force. With gravity assist, however, the flow circulation can be driven by the combination of capillary and gravitational forces, or by the gravitational force alone. For a gravity-assist LHP at a given elevation between the horizontal condenser and evaporator, there exists a threshold heat load below which the LHP operation is gravity driven and above which the LHP operation is capillary force and gravity co-driven. The gravitational pressure head can have profound effects on the LHP operation, and such effects depend on the elevation, evaporator heat load, and condenser sink temperature. This paper presents a theoretical study on LHP operations under gravity-neutral, anti-gravity, and gravity-assist modes using pressure diagrams to help understand the underlying physical processes. Effects of the condenser configuration on the gravitational pressure head and LHP operation are also discussed.

Thermal Control Systems↗
NASA NTRS1999

MEMS Micropropulsion Activities at JPL

A status of MEMS-based micropropulsion activities conducted at JPL will be given. These activities include work conducted on the so called Vaporizing Liquid Micro-Thruster (VLM) which recently underwent proof-of-concept testing, demonstrating the ability to vaporize water propellant at 2 W and 2 V. Micro-ion engine technologies, such m field emitter arrays and micro-grids are being studied. Focus in the field emitter area is on arrays able to survive in thruster plumes and micro-ion engine plasmas to serve as neutralizers aW engine cathodes. Integrated, batch-fabricated Ion repeller grid structures are being studied as well as different emitter tip materials are being investigated to meet these goals. A micro-isolation valve is being studied to isolate microspacecraft feed system during long interplanetary cruises, avoiding leakage and prolonging lifetime and reliability of such systems. This concept relies on the melting of a thin silicon barrier. Burst pressure values as high as 2,900 psig were obtained for these valves and power requirements to melt barriers ranging between 10 - 50 microns in thickness, as determined through thermal finite element calculations, varied between 10 - 30 W to be applied over a duration of merely 0.5 ms.

Mueller, Juergen↗
NASA NTRS2003

Understanding the Wind/Ejecta of ETA Carinae with the STIS

The mysterious Luminous Blue Variable, Eta Carinae, and its ejecta, the Homunculus, are yielding much new and exciting information to the Space Telescope Imaging Spectrograph (STIS). With 0.05 in. pixels and R=5000, we have followed changes in the Central Source and the surrounding ejecta over the past 5.5 year period in the spectral region 1650-10300A. Neutral, partially-ionized and highly-excited nebular regions have been resolved. With 0.025 in., R=30000 and R=120000 between 1180 and 3160A, we are following the changes in the Central Source, the wind structure and at least 20 absorption components, thought to be within the wall of the Humunculus. The Central Source proves to be slightly extended at HST angular resolutions. Theoretical models of the nebular line fluxes are providing new insight to the sources of excitation. Photo-excitation, line and continuum pumping, plus shocks, all play a role. Laboratory measurements of Fe-like elements and improved atomic data bases prove necessary for proper interpretation. Increased angular resolution by larger-aperture space telescopes and improved optics and detectors will lead to new understanding of massive stars, their wind structures and history of ejecta.

Gull, Theodore R.↗
DOE OSTI2026

Numerical-heating effects in atmospheric pressure streamer discharges simulated with a PIC code

Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.

Nikic, Dejan [University of New Mexico, Albuquerqu↗
DOE OSTI2021

A Graphics Processing Unit (GPU) Approach to Large Eddy Simulation (LES) for Transport and Contaminant Dispersion

Recent advances in the development of large eddy simulation (LES) atmospheric models with corresponding atmospheric transport and dispersion (AT&D) modeling capabilities have made it possible to simulate short, time-averaged, single realizations of pollutant dispersion at the spatial and temporal resolution necessary for common atmospheric dispersion needs, such as designing air sampling networks, assessing pollutant sensor system performance, and characterizing the impact of airborne materials on human health. The high computational burden required to form an ensemble of single-realization dispersion solutions using an LES and coupled AT&D model has, until recently, limited its use to a few proof-of-concept studies. An example of an LES model that can meet the temporal and spatial resolution and computational requirements of these applications is the joint outdoor-indoor urban large eddy simulation (JOULES). A key enabling element within JOULES is the computationally efficient graphics processing unit (GPU)-based LES, which is on the order of 150 times faster than if the LES contaminant dispersion simulations were executed on a central processing unit (CPU) computing platform. JOULES is capable of resolving the turbulence components at a suitable scale for both open terrain and urban landscapes, e.g., owing to varying environmental conditions and a diverse building topology. In this paper, we describe the JOULES modeling system, prior efforts to validate the accuracy of its meteorological simulations, and current results from an evaluation that uses ensembles of dispersion solutions for unstable, neutral, and stable static stability conditions in an open terrain environment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Development of Individually Addressable Micro-Mirror-Arrays for Space Applications

We have been developing a 32 x 32 prototype array of individually addressable Micro-Mirrors capable of operating at cryogenic temperature for Earth and Space Science applications. Micro-Mirror-Array technology has the potential to revolutionize imaging and spectroscopy systems for NASA's missions of the 21st century. They can be used as programmable slits for the Next Generation Space Telescope, as smart sensors for a steerable spectrometer, as neutral density filters for bright scene attenuation etc. The, entire fabrication process is carried out in the Detector Development Laboratory at NASA, GSFC. The fabrication process is low temperature compatible and involves integration of conventional CMOS technology and surface micro-machining used in MEMS. Aluminum is used as the mirror material and is built on a silicon substrate containing the CMOS address circuit. The mirrors are 100 microns x l00 microns in area and deflect by +/- 10 deg induced by electrostatic actuation between two parallel plate capacitors. A pair of thin aluminum torsion straps allow the mirrors to tilt. Finite-element-analysis and closed form solutions using electrostatic and mechanical torque for mirror operation were developed and the results were compared with laboratory performance. The results agree well both at room temperature and at cryogenic temperature. The development demonstrates the first cryogenic operation of two-dimensional Micro-Mirrors with bi-state operation. Larger arrays will be developed meeting requirements for different science applications. Theoretical analysis, fabrication process, laboratory test results and different science applications will be described in detail.

Dutta, Sanghamitra B.↗
DOE OSTI2023

Elements of H-mode pedestal structure

Abstract This paper reviews current understanding of key physics elements that control the H-mode pedestal structure, which exists at the boundary of magnetically confined plasmas. The structure of interest is the width, height and gradient of temperature, density and pressure profiles in the pedestal. Emphasis is placed on understanding obtained from combined experimental, theoretical and simulation work and on results observed on multiple machines. Pedestal profiles are determined by the self-consistent interaction of sources, transport and magnetohydrodynamic limits. The heat source is primarily from heat deposited in the core and flowing to the pedestal. This source is computed from modeling of experimental data and is generally well understood. Neutrals at the periphery of the plasma provide the dominant particle source in current machines. This source has a complex spatial structure, is very difficult to measure and is poorly understood. For typical H-mode operation, the achievable pedestal pressure is limited by repetitive, transient magnetohydrodynamic instabilities. First principles models of peeling–ballooning modes are generally able to explain the observed limits. In some regimes, instability occurs below the predicted limits and these remain unexplained. Several mechanisms have been identified as plausible sources of heat transport. These include neoclassical processes for ion heat transport and several turbulent processes, driven by the steep pedestal gradients, as sources of electron and ion heat transport. Reduced models have successfully predicted the pedestal or density at the pedestal top. Firming up understanding of heat and particle transport remains a primary challenge for developing more complete predictive pedestal models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A Combined Computational and Mathematical Analysis of Interconnect Fatigue Potential in Photovoltaic Modules

A finite element model of a 60-cell monocrystalline silicon glass-polymer photovoltaic module was simulated with ±1.0 kPa and ±2.4 kPa loads applied to the glass to calculate the deformation under load. Cell-to-cell displacements were used to approximate interconnect strain and stress. A mathematical fatigue cycle life relation was fitted to data for the interconnect material (copper), to generate a life prediction at each interconnect location based on the local stress means, reversal extents, and amplitudes. Interconnect stress was found to be significantly asymmetric about zero despite symmetric positive and negative module loads due to laminate thickness offsets about the neutral plane and the effects of module framing. Cycle life results indicated that interconnect fatigue failure was unlikely to occur over a 30-year lifetime of conservative wind and snow load cycles since the typical cell design feature of leaving some unconstrained length between the cell edge and first solder pad increases the effective gauge length and decreases the stress levels below the material endurance limit. Follow-up analyses found that 3.6 mm and 6.4 mm were the minimum unconstrained lengths required to survive the assumed lifetime of wind and snow cycles, respectively, confirming that typical industrial module constructions with 8–15 mm unconstrained lengths should survive conservatively. Notably, large magnitude, low-cycle snow loading was consistently the limiting factor requiring a longer unconstrained interconnect length. Finally, insights and workflows from this study inform module interconnection design limits for survival against mechanical fatigue in deployment environments.

14 SOLAR ENERGY↗
DOE OSTI2022

The strangeness program at GlueX

The GlueX experiment located at Jefferson Lab studies the spectrum of hadrons using photoproduction on a LH2 target in a wide variety of final states. With its detector system capable of measuring neutral and charged final state particles over almost the full solid angle, and very good particle iden- tification capabilities, GlueX can measure many different hadrons containing strangeness. A linearly polarized photon beam allows the measurement of po- larization observables, which contain information about the production mecha- nisms involved in generating strange particles in photoproduction. In addition, GlueX can perform precise cross-section measurements, which help to study the spectrum of strange hadrons. In this presentation, the GlueX experiment is introduced, and recent progress of its strangeness program is discussed. We present recent results on ?0 beam asymmetries, ?(1520) spin-density matrix elements and ongoing studies of the ?(1405) lineshape. We also present our recent progress on measurements of ?? ? and ?(?) photoproduction. Also, future prospects for strangeness measurements at GlueX are discussed.

Pauli, Peter↗
DOE OSTI2026

Structural and mechanistic insights into protective non-neutralizing antibodies targeting Crimean-Congo hemorrhagic fever virus nucleocapsid protein

Abstract Crimean-Congo Hemorrhagic Fever Virus (CCHFV) is a tick-borne virus endemic to Africa, Asia, and expanding regions within Europe. With mortality rates approaching 40%, rising incidence, and no currently approved countermeasures, CCHFV is recognized as a priority public health threat. CCHFV nucleocapsid protein (NP) has long been a key target for diagnostics. Recently, NP-specific humoral responses have also been correlated with protection conferred by protective vaccines candidates. Additionally, the first non-neutralizing monoclonal antibody (mAb) 9D5 demonstrated protective efficacy against CCHFV challenge, underscoring NP as a viable antiviral target. Here, nine anti-NP mAb were utilized to identify four antigenic sites on NP as well as localize these sites to the head or stalk domains. These mAb also revealed variable levels of in vivo protection, independent from whether the epitope site is located in the head or stalk regions. Additionally, three X-ray crystallography structures were obtained that included CCHFV NP from strain Afg09-2990 in complex with the most potent mAb (9D5). This, along with additional structures of two unbound NPs, revealed structural elements critical for mAb-9D5 broad-spectrum protective characteristics. These findings provide a path towards the rapid identification of broadly protective anti-NP mAb countermeasures.

Moresco, Vanessa↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗