Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “literature”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Long-term multi-wavelength study of 1ES 0647+250

The BL Lac object 1ES 0647+250 is one of the few distant γ-ray emitting blazars detected at very high energies (VHEs; &100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during a period of low activity in the years 2009-2011 as well as during three flaring activities in the years 2014, 2019, and 2020, with the highest VHE flux in the last epoch. An extensive multi-instrument data set was collected as part of several coordinated observing campaigns over these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluated the variability and correlation of the emission in the different energy bands with the fractional variability and the Z transformed discrete correlation function, as well as its spectral evolution in X-rays and γ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimated its distance in an indirect manner through a comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data. We find significant long-term variability, especially in X-rays and VHE γ rays. Furthermore, significant (3-4σ) correlations were found between the radio, optical, and high-energy (HE) γ-ray fluxes, with the radio emission delayed by about ~400 days with respect to the optical and γ-ray bands. The spectral analysis reveals a harder-when-brighter trend during the non-flaring state in the X-ray domain. However, no clear patterns were observed for either the enhanced states or the HE (30 MeV < E < 100 GeV) and VHE γ-ray emission of the source. The indirect estimation of the redshift yielded a value of z = 0.45 ± 0.05, which is compatible with some of the values reported in the literature. The SEDs related to the low-activity state and the three flaring states of 1ES 0647+250 can be described reasonably well with the both one-component and two-component leptonic scenarios. However, the long-term correlations indicate the need for an additional radio-producing region located about 3.6 pc downstream from the gamma-ray producing region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

In-depth analysis of solar models with high-metallicity abundances and updated opacity tables

Context. As a result of the high-quality constraints available for the Sun, we are able to carry out detailed combined analyses using neutrino, spectroscopic, and helioseismic observations. These studies lay the ground for future improvements of the key physical components of solar and stellar models because ingredients such as the equation of state, the radiative opacities, or the prescriptions for macroscopic transport processes of chemicals are then used to study other stars in the Universe. Aims. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss the effect on their properties of changes in the micro and macro physical ingredients. Methods. We carried out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged 3D models that were claimed to resolve the solar modelling problem. We compared these models to helioseismic and neutrino constraints. Results. The properties of the solar models are significantly affected by the use of the recent OPLIB opacity tables and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those for which the OP opacities were used. We show that a modification of the temperature gradient just below the base of the convective zone is required to remove the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a localised increase in the opacity of a few percent. Conclusions. We conclude that the existing degeneracies and issues in solar modelling are not removed by using an increase in the solar metallicity, in contradiction to what has been suggested in the recent literature. Therefore, standard solar models cannot be used as an argument for a high-metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

All-electron ab initio Bethe-Salpeter equation approach to neutral excitations in molecules with numeric atom-centered orbitals

The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as “Thiel’s set.” Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted “tier2+aug2” is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.

Liu, Chi↗
DOE OSTI2020

Comparing GGA, GGA+ U , and meta-GGA functionals for redox-dependent binding at open metal sites in metal–organic frameworks

Metal–organic frameworks (MOFs) with open metal sites have been widely investigated for the selective adsorption of small molecules via redox mechanisms where charge transfer can take place between the binding site and the adsorbate of interest. Quantum-chemical screening methods based on density functional theory have emerged as a promising route to accelerate the discovery of MOFs with enhanced binding affinities toward various adsorbates. However, the success of this approach is linked to the accuracy of the underlying density functional approximations (DFAs). In this work, we compare commonly used generalized gradient approximation (GGA), GGA+U, and meta-GGA exchange-correlation functionals in modeling redox-dependent binding at open metal sites in MOFs using O 2 and N 2 as representative small molecules. We find that the self-interaction error inherent to the widely used Perdew, Burke, and Ernzerhof (PBE) GGA predicts metal sites that are artificially redox-active, as evidenced by their strong binding affinities, short metal–adsorbate bond distances, and large degree of charge transfer. The incorporation of metal-specific, empirical Hubbard U corrections based on the transition metal oxide literature systematically reduces the redox activity of the open metal sites, often improving agreement with experiment. Additionally, the binding behavior shifts from strong chemisorption to weaker physisorption as a function of U. The M06-L meta-GGA typically predicts binding energies between those of PBE-D3(BJ) and PBE-D3(BJ)+U when using empirically derived U values from the transition metal oxide literature. Despite the strong sensitivity of the binding affinities toward a given DFA, the GGA, GGA+U, and meta-GGA approaches often yield the same qualitative trends and structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Data-driven materials research enabled by natural language processing and information extraction

Given the emergence of data science and machine learning throughout all aspects of society, but particularly in the scientific domain, there is increased importance placed on obtaining data. Data in materials science are particularly heterogeneous, based on the significant range in materials classes that are explored and the variety of materials properties that are of interest. This leads to data that range many orders of magnitude, and these data may manifest as numerical text or image-based information, which requires quantitative interpretation. The ability to automatically consume and codify the scientific literature across domains - enabled by techniques adapted from the field of natural language processing - therefore has immense potential to unlock and generate the rich datasets necessary for data science and machine learning. This review focuses on the progress and practices of natural language processing and text mining of materials science literature and highlights opportunities for extracting additional information beyond text contained in figures and tables in articles. Here, we discuss and provide examples for several reasons for the pursuit of natural language processing for materials, including data compilation, hypothesis development, and understanding the trends within and across fields. Current and emerging natural language processing methods along with their applications to materials science are detailed. We, then, discuss natural language processing and data challenges within the materials science domain where future directions may prove valuable.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experimental observations of exploding bridgewire detonator function

Exploding bridgewire detonators are an industry standard technology used for over 75 years and valued for their precise timing and safety characteristics. Despite widespread use, their functional mechanism remains controversial with both shock and non-shock mechanisms attributed. In this work, we re-examine bridgewire detonator function with a suite of modern diagnostics and compare these observations with the existing literature. Traditional detonator observations consisted of voltage applied to the bridgewire and time dependent current, integral response measurements such as case motion and more recently Schlieren imaging of the detonator surface. In this work, we add visible light emission, X-ray transmission, proton radiography, and temperature measurements during detonator function in addition to voltage, current, and function times. The addition of in-situ observations of light emission, temperature, and density gives us new insight into the mechanisms of explosive bridgewire detonator function. We see a distinct separation in time, location, symmetry, and velocity of bridgewire output and detonation onset. During the time between bridgewire burst and the initiation of detonation, we observe a temperature ramp in the input pellet. In this paper, we present the suite of measurements and comparisons with the literature on integral response measurements.

42 ENGINEERING↗
DOE OSTI2021

An automated multi-component gas adsorption system (MC GAS)

The knowledge gap on adsorption of complex mixtures in the literature relative to single component data represents a persistent obstacle to developing accurate process models for adsorption separations. The collection of mixed gas adsorption data is an imminent need for improved understanding of the behavior of adsorbent systems in these diverse adsorption applications. Current approaches to understanding mixture adsorption using predictive theories based on pure component adsorption experiments often fail to capture the behavior of more complex, non-ideal systems. In this work, we present an automated volumetric instrument for the measurement of mixed gas adsorption isotherms. This instrument was validated by comparison to other in-house instruments and data available in the literature, and the binary adsorption measurements were found to be thermodynamically consistent. The automation of this instrument allows for rapid collection of high-quality mixture adsorption data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

The role of pre-existing heterogeneities in materials under shock and spall

There has been a challenge for many decades to understand how heterogeneities influence the behavior of materials under shock loading, eventually leading to spall formation and failure. Experimental, analytical, and computational techniques have matured to the point where systematic studies of materials with complex microstructures under shock loading and the associated failure mechanisms are feasible. This is enabled by more accurate diagnostics as well as characterization methods. As interest in complex materials grows, understanding and predicting the role of heterogeneities in determining the dynamic behavior becomes crucial. Early computational studies, hydrocodes, in particular, historically preclude any irregularities in the form of defects and impurities in the material microstructure for the sake of simplification and to retain the hydrodynamic conservation equations. Contemporary computational methods, notably molecular dynamics simulations, can overcome this limitation by incorporating inhomogeneities albeit at a much lower length and time scale. This review discusses literature that has focused on investigating the role of various imperfections in the shock and spall behavior, emphasizing mainly heterogeneities such as second-phase particles, inclusions, and voids under both shock compression and release. Pre-existing defects are found in most engineering materials, ranging from thermodynamically necessary vacancies, to interstitial and dislocation, to microstructural features such as inclusions, second phase particles, voids, grain boundaries, and triple junctions. This literature review explores the interaction of these heterogeneities under shock loading during compression and release. Systematic characterization of material heterogeneities before and after shock loading, along with direct measurements of Hugoniot elastic limit and spall strength, allows for more generalized theories to be formulated. Further, continuous improvement toward time-resolved, in situ experimental data strengthens the ability to elucidate upon results gathered from simulations and analytical models, thus improving the overall ability to understand and predict how materials behave under dynamic loading.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Moiré engineering in van der Waals heterostructures

Isolated atomic planes can be assembled into a multilayer van der Waals (vdW) heterostructure in a precisely chosen sequence. These heterostructures feature moiré patterns if the constituent 2D material layers are stacked in an incommensurable way, due to a lattice mismatch or twist. This design-by-stacking has opened up the promising area of moiré engineering, a term that can be understood in two different perspectives, namely, (i) structural—engineering a moiré pattern by introducing twist, relative strain, or defects that affect the commensurability of the layers and (ii) functional—exploiting a moiré pattern to find and tune resulting physical properties of a vdW heterostructure. The latter meaning, referring to the application of a moiré pattern, is seen in the literature in the specific context of the observation of correlated electronic states and unconventional superconductivity in twisted bilayer graphene. The former meaning, referring to the design of the moiré pattern itself, is present in the literature but less commonly discussed or less understood. The underlying link between these two perspectives lies in the deformation field of the moiré superlattice. In this Perspective, we describe a path from designing a moiré pattern to employing the moiré pattern to tune physical properties of a vdW heterostructure. We also discuss the concept of moiré engineering in the context of twistronics, strain engineering, and defect engineering in vdW heterostructures. Although twistronics is always associated with moiré superlattices, strain and defect engineering are often not. Here, in this work, we demonstrate how strain and defect engineering can be understood within the context of moiré engineering. Adopting this perspective, we note that moiré engineering creates a compelling opportunity to design and develop multiscale electronic devices.

2D materials↗
DOE OSTI2023

Advanced data analysis in inertial confinement fusion and high energy density physics

Bayesian analysis enables flexible and rigorous definition of statistical model assumptions with well-characterized propagation of uncertainties and resulting inferences for single-shot, repeated, or even cross-platform data. This approach has a strong history of application to a variety of problems in physical sciences ranging from inference of particle mass from multi-source high-energy particle data to analysis of black-hole characteristics from gravitational wave observations. The recent adoption of Bayesian statistics for analysis and design of high-energy density physics (HEDP) and inertial confinement fusion (ICF) experiments has provided invaluable gains in expert understanding and experiment performance. In this Review, we discuss the basic theory and practical application of the Bayesian statistics framework. We highlight a variety of studies from the HEDP and ICF literature, demonstrating the power of this technique. Due to the computational complexity of multi-physics models needed to analyze HEDP and ICF experiments, Bayesian inference is often not computationally tractable. Two sections are devoted to a review of statistical approximations, efficient inference algorithms, and data-driven methods, such as deep-learning and dimensionality reduction, which play a significant role in enabling use of the Bayesian framework. We provide additional discussion of various applications of Bayesian and machine learning methods that appear to be sparse in the HEDP and ICF literature constituting possible next steps for the community. We conclude by highlighting community needs, the resolution of which will improve trust in data-driven methods that have proven critical for accelerating the design and discovery cycle in many application areas.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electrical Analysis of Pulsed Laser Annealed Poly-Si: Ga/SiOx Passivating Contacts

In this contribution, we examine Ga hyper-doped poly-Si/SiOx contacts realized by pulsed laser melting (PLM). Here, we use Ga as a novel p-type dopant and B as a conventional dopant to induce non-equilibrium doping using an excimer laser. We perform simulations to visualize the maximum melt depth profiles within the poly-Si, with a goal of distributing dopants close to the tunneling oxide, but at the same time preserving the passivation. Hall measurements show that sheet resistance for B is lower than Ga due to its higher solid solubility limit in Si. After comparing the Hall active dopant concentration with the chemical concentration obtained by SIMS measurement, we show nearly 100% activation B activation reaching 10^21 cm-3, while only ~20% activation for Ga. Nevertheless, we achieve active doping concentrations of Ga in poly-Si six times higher than its solid solubility limit in Si (~10^19 cm-3). We compare our Hall mobilities with values in the literature for c-Si and show that B mobilities in laser-treated poly-Si are close to that of the literature value for B, while Ga mobilities are lower, possibly due to additional scattering channels within grain boundaries and deformed lattice. We also compare our results on PLM samples with conventional furnace annealed samples, and we show much higher percent activation and mobilities. Previously, we showed a low contact resistivity of 35.5 +/- 2.4 m..omega..cm2. Here, we further confirm this result by scanning spreading resistance microscopy and Kelvin force nanoprobe microscopy. We demonstrate that our poly-Si: Ga/nCz contact exhibits large drift and diffusion currents under normal cell operating voltage, which widens the laser processing window for a good metal/poly-Si/c-Si contact.

high-efficiency solar cells↗
DOE OSTI2024

The bond energy of UN + : Guided ion beam studies of the reactions of U + with N 2 and NO

In this work, a guided ion beam tandem mass spectrometer was used to study the reactions of U + with N 2 and NO. Reaction cross sections were measured over a wide range of energy for both systems. In each reaction, UN + is formed by an endothermic process, thereby enabling the direct measurement of the threshold energy and determination of the UN + bond dissociation energy. For the reaction of U + + N 2 , a threshold energy (E 0 ) of 4.02 ± 0.11 eV was measured, leading to D 0 (UN + ) = 5.73 ± 0.11 eV. The reaction of U + + NO yields UO + through an exothermic, barrierless process that proceeds with 94 ± 23% efficiency at the lowest energy. Analysis of the endothermic UN + cross section in this reaction provides E 0 = 0.72 ± 0.11 eV and, therefore, D 0 (UN + ) = 5.78 ± 0.11 eV. Averaging the values obtained from both reactions, we report D 0 (UN + ) = 5.76 ± 0.13 eV as our best value (uncertainty of two standard deviations). Combined with precise literature values for the ionization energies of U and UN, we also derive D 0 (UN) = 5.86 ± 0.13 eV. Both bond dissociation energies agree well with high-level theoretical treatments in the literature. The formation of UN + in reaction of U + with NO also exhibits a considerable increase in reaction probability above ~3 eV. Theory suggests that this may be consistent with the formation of UN + in excited quintet spin states, which we hypothesize are dynamically favored because the number of 5f electrons in reactants and products is conserved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Hiding in plain sight: The prevalence and impact of trions and Fermi polarons in transient absorption spectroscopy experiments of 2D semiconductors

Transient absorption (TA) spectroscopy is one of the most popular experimental methods to measure the excited state lifetimes and charge carrier recombination mechanisms in two dimensional (2D) semiconductors. This fundamental information is essential for designing and optimizing the next generation of ultrathin and lightweight 2D semiconductor-based optoelectronic devices. However, the interpretation of TA spectroscopy data varies across the community. The community lacks a unifying physical explanation for how and why experimental variables such as incident light intensity, sample-substrate interactions, and/or applied bias affect TA spectral data. This Perspective (1) compares the physical chemistry TA literature to nanomaterial physics literature from a historical perspective, (2) reviews multiple physical explanations that the TA community developed to explain spectral features and experimental trends, (3) provides a unifying explanation for how and why trions—and, more generally, Fermi polarons—contribute to TA spectra, and (4) quantifies the extent to which various physical interpretations and data analysis procedures yield different timescales and mechanisms for the same set of experimental results. We highlight the importance of considering trions/Fermi polarons in TA measurements and their implications for advancing our understanding of 2D material properties.

2D materials↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗