SEARCH · Search NASA

Results for “geometrical parametrization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Sensitivity analysis of design parameters in RowWise borehole layout for ground heat exchangers

Ground source heat pump systems offer a promising pathway toward energy-efficient building heating and cooling. The performance and cost-effectiveness of these systems heavily depend on the design of the ground heat exchanger (GHE), particularly the spatial placement of boreholes used in large commercial buildings. Among various borefield layout strategies, the RowWise approach generates and optimizes borehole configurations within irregular polygonal land boundaries, providing land-use efficiency and installation flexibility. While multiple geometric design parameters constrain the RowWise layout, their influence on the system thermal performance and total drilling requirements remains unclear. Thus, this study presents a sensitivity analysis of key design parameters influencing the RowWise layout of vertical borehole GHEs, including the perimeter spacing ratio, borehole spacing, borehole field rotation angle, and borehole length. GHEDesigner and Ray Tune are employed to generate and assess different RowWise configurations. A series of parametric simulations was conducted to quantify the impact of each parameter on geometric distribution, economic cost, and computational speed. The results provide critical insights into the sensitivity and relative importance of different design parameters of RowWise method, offering practical guidance for designers aiming to optimize GHE layouts by balancing thermal efficiency, land constraints, and economic feasibility.

Xu, Dikai [Purdue University]↗

DOE OSTI2022

First Sagittarius A* Event Horizon Telescope Results. IV. Variability, Morphology, and Black Hole Mass

In this paper we quantify the temporal variability and image morphology of the horizon-scale emission from Sgr A*, as observed by the EHT in 2017 April at a wavelength of 1.3 mm. We find that the Sgr A* data exhibit variability that exceeds what can be explained by the uncertainties in the data or by the effects of interstellar scattering. The magnitude of this variability can be a substantial fraction of the correlated flux density, reaching ~100% on some baselines. Through an exploration of simple geometric source models, we demonstrate that ring-like morphologies provide better fits to the Sgr A* data than do other morphologies with comparable complexity. We develop two strategies for fitting static geometric ring models to the time-variable Sgr A* data; one strategy fits models to short segments of data over which the source is static and averages these independent fits, while the other fits models to the full data set using a parametric model for the structural variability power spectrum around the average source structure. Both geometric modeling and image-domain feature extraction techniques determine the ring diameter to be 51.8 ± 2.3 μas (68% credible intervals), with the ring thickness constrained to have an FWHM between ~30% and 50% of the ring diameter. To bring the diameter measurements to a common physical scale, we calibrate them using synthetic data generated from GRMHD simulations. This calibration constrains the angular size of the gravitational radius to be ${4.8}_{-0.7}^{+1.4}$ μas, which we combine with an independent distance measurement from maser parallaxes to determine the mass of Sgr A* to be ${4.0}_{-0.6}^{+1.1}\times {10}^{6}$ M⊙.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Improved galactic foreground removal for B-mode detection with clustering methods

ABSTRACT Characterizing the sub-mm Galactic emission has become increasingly critical especially in identifying and removing its polarized contribution from the one emitted by the cosmic microwave background (CMB). In this work, we present a parametric foreground removal performed on to sub-patches identified in the celestial sphere by means of spectral clustering. Our approach takes into account efficiently both the geometrical affinity and the similarity induced by the measurements and the accompanying errors. The optimal partition is then used to parametrically separate the Galactic emission encoding thermal dust and synchrotron from the CMB one applied on two nominal observations of forthcoming experiments from the ground and from the space. Moreover, the clustering is performed on tracers that are different from the data used for component separation, e.g. the spectral index maps of dust and synchrotron. Performing the parametric fit singularly on each of the clustering derived regions results in an overall improvement: both controlling the bias and the uncertainties in the CMB B-mode recovered maps. We finally apply this technique using the map of the number of clouds along the line of sight, $\mathcal {N}_c$, as estimated from H i emission data and perform parametric fitting on to patches derived by clustering on this map. We show that adopting the $\mathcal {N}_c$ map as a tracer for the patches related to the thermal dust emission, results in reducing the B-mode residuals post-component separation. The code is made publicly available https://github.com/giuspugl/fgcluster.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Bayesian optimization to design a novel x-ray shaping device

In radiation therapy, x-ray dose must be precisely sculpted to the tumor, while simultaneously avoiding surrounding organs at risk. This requires modulation of x-ray intensity in space and/or time. Typically, this is achieved using a multi leaf collimator (MLC) - a complex mechatronic device comprising over one hundred individually powered tungsten ‘leaves’ that move in or out of the radiation field as required. Here, an all-electronic x-ray collimation concept with no moving parts is presented, termed “SPHINX”: Scanning Pencil-beam High-speed Intensity-modulated X-ray source. SPHINX utilizes a spatially distributed bremsstrahlung target and collimator array in conjunction with magnetic scanning of a high energy electron beam to generate a plurality of small x-ray “beamlets.” A simulation framework was developed in Topas Monte Carlo incorporating a phase space electron source, transport through user defined magnetic fields, bremsstrahlung x-ray production, transport through a SPHINX collimator, and dose in water. This framework was completely parametric, meaning a simulation could be built and run for any supplied geometric parameters. This functionality was coupled with Bayesian optimization to find the best parameter set based on an objective function which included terms to maximize dose rate for a user defined beamlet width while constraining inter-channel cross talk and electron contamination. Designs for beamlet widths of 5, 7, and 10 mm 2 were generated. Each optimization was run for 300 iterations and took approximately 40 h on a 24-core computer. For the optimized 7-mm model, a simulation of all beamlets in water was carried out including a linear scanning magnet calibration simulation. Finally, a back-of-envelope dose rate formalism was developed and used to estimate dose rate under various conditions. The optimized 5–, 7–, and 10-mm models had beamlet widths of 5.1 , 7.2 , and 10.1 mm 2 and dose rates of 3574, 6351, and 10 015 Gy/C, respectively. The reduction in dose rate for smaller beamlet widths is a result of both increased collimation and source occlusion. For the simulation of all beamlets in water, the scanning magnet calibration reduced the offset between the collimator channels and beam centroids from 2.9 ±1.9 mm to 0.01 ±0.03 mm. A slight reduction in dose rate of approximately 2% per degree of scanning angle was observed. Based on a back-of-envelope dose rate formalism, SPHINX in conjunction with next-generation linear accelerators has the potential to achieve substantially higher dose rates than conventional MLC-based delivery, with delivery of an intensity modulated 100 x 100 mm 2 field achievable in 0.9 to 10.6 s depending on the beamlet widths used. Bayesian optimization was coupled with Monte Carlo modeling to generate SPHINX geometries for various beamlet widths. A complete Monte Carlo simulation for one of these designs was developed, including electron beam transport of all beamlets through scanning magnets, x-ray production and collimation, and dose in water. These results demonstrate that SPHINX is a promising candidate for sculpting radiation dose with no moving parts, and has the potential to vastly improve both the speed and robustness of radiotherapy delivery. A multi-beam SPHINX system may be a candidate for delivering magavoltage FLASH RT in humans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2021

Crossing Symmetric Dispersion Relations in Quantum Field Theories

For 2-2 scattering in quantum field theories, the usual fixed t dispersion relation exhibits only two-channel symmetry. This Letter considers a crossing symmetric dispersion relation, reviving certain old ideas from the 1970s. Rather than the fixed t dispersion relation, this needs a dispersion relation in a different variable z, which is related to the Mandelstam invariants s, t, u via a parametric cubic relation making the crossing symmetry in the complex z plane a geometric rotation. The resulting dispersion is manifestly three-channel crossing symmetric. We give simple derivations of certain known positivity conditions for effective field theories, including the null constraints, which lead to two sided bounds and derive a general set of new nonperturbative inequalities. We show how these inequalities enable us to locate the first massive string state from a low energy expansion of the four dilaton amplitude in type II string theory. We also show how a generalized (numerical) Froissart bound, valid for all energies, is obtained from this approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

Harp, Dylan R. (ORCID:0000000197778000)↗

DOE OSTI2023

Computational Fluid Dynamics Modeling of Solar Thermal Dry Reforming of Methane in a Parabolic Trough

Computational fluid dynamics simulations of solar-thermal dry reforming of methane using a parabolic trough configuration were performed. Parametric simulations of different combinations of gas flow rate, receiver tube emissivity, and geometric concentration ratio were conducted to determine configurations that could achieve the required catalyst temperatures of at least 700 °C to achieve high conversion of CH4 and CO2 to H2 and CO. Results showed that the concentration ratio of the parabolic trough collector had to be increased from ~70 to ~120 and the receiver-tube emissivity had to be reduced to ~0.2 to achieve bulk average catalyst temperatures of greater than 700 °C. Lower gas flow rates also reduced enthalpic heat losses and increased catalyst temperatures.

Ho, Clifford↗

DOE OSTI2021

Experiment and non-local crystal plasticity finite element study of nanoindentation on Al-8Ce-10Mg alloy

Cerium and magnesium strengthened aluminum alloys, or Al-Ce-Mg, is a recently developed alloy family that exhibits good mechanical properties at elevated temperatures (~300 °C). To examine the single-crystal properties of Al-Ce-Mg alloys, nanoindentation experiments are conducted in this study. A crystal plasticity finite element model (CPFEM) with the evolution of geometrically necessary dislocations (GNDs) is applied to simulate the indentation in individual grains. A parametric study is carried out to investigate the sensitivity of each crystal plasticity model parameter to the indentation behavior. The highly sensitive parameters are calibrated by matching the indentation load-depth curves, while the rest parameters are obtained from bulk polycrystal uniaxial tension tests. Overall, satisfactory matching between experiment and simulation is obtained for each individual grain. The calculated hardness, as determined from the experiment, shows the dependence on indent depth, which is captured by the GND model. Furthermore, the effect of grain orientation and neighboring grains to nanoindentation behavior have been discussed with the comparison between the simulation and experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Optimization of Conductive Fins to Minimize UO2 Fuel Temperature and Radial Temperature Gradient

To further the development of low-enriched uranium fuels, precedence has been placed on delivering the same amount of power while lowering the fuel temperature and radial temperature gradient. To address this, modeling efforts have resulted in a novel design featuring conductive fins of varying thermal conductivities and geometries inserted into the fuel matrix. These conductive inserts were not allowed to exceed 6% of the original fuel volume. This constraint was imposed due to other designs displacing 10% of fuel volume. A parametric study was performed that consisted of 2.56 million BISON simulations involving varying fin characteristics (i.e., fin thermal conductivity, number, and geometry) to determine the optimal geometric configuration for a desired amount of fuel volume displaced. The results from this study show that the thickness and length of each fin affect the fuel temperature and temperature gradient more than varying the number and thermal conductivity of the fins. The parametric study resulted in the development of an optimized combination to produce the lowest peak fuel temperature, lowest radial temperature gradient, and highest temperature reduction for the amount of original fuel volume displaced. The simulations presented in this work will eventually be compared with irradiation experiments of similar fuel designs at Idaho National Laboratory’s Advanced Test Reactor.

Paaren, Kyle M.↗

DOE OSTI2022

An additively-manufactured molten salt-to-supercritical carbon di-oxide primary heat exchanger for solar thermal power generation – Design and techno-economic performance

The design and techno-economic performance of a compact additively manufactured (AM) molten salt (MS)-to-supercritical carbon di-oxide (sCO2) primary heat exchanger (PHE) for solar thermal application is described. The PHE design consists of sCO2 flow through an array of microscale pin fins while the MS flows through mm-scale rectangular channels. Constraints imposed by AM using laser powder bed fusion method are considered in the design. Structural and fluid flow simulations are performed to arrive at a viable design of the core and headers. A simplified one-dimensional steady state model for the PHE is developed including the impact of surface roughness from the AM process. A process-based cost model is used to determine the tradeoff between thermofluidic design and manufacturing cost. A parametric study is performed using the thermo-fluidic and cost models to determine the set of geometrical and flow variables that result in high power density and low cost, while restricting the pressure drop on the sCO2 side to less than 2% of line pressure. Flow rates of MS and sCO2 were varied over heat capacity rate ratios ranging from 0.2 to 1. Results indicate that it is possible to design a low-pressure drop AM PHE with an effectiveness of 90% and a power density in excess of 10 MW/m3 (including headers). Fabrication of representative nickel superalloy specimens are shown to demonstrate that low-porosity parts with the requisite dimensional tolerance of PHE core can be generated.

14 SOLAR ENERGY↗

DOE OSTI2021

Process Robustness of Friction Stir Dovetailing of AA7099 to Steel with In Situ AA6061 Interlayer Linking

Rolled homogeneous armor (RHA) steel was joined to aluminum alloy AA7099-T7451 with an interlayer of AA6061-T651 by single-pass friction stir dovetailing (FSD). The evolution of the AA7099/AA6061 and AA6061/steel interfacial microstructures and mechanical properties of the joints were examined. The AA6061 interlayer prevented Zn migration from the AA7099 to the steel during FSD and inhibited formation of brittle, Zn-rich, Fe-Al intermetallic compounds within the joint. The flow of the dissimilar aluminum alloys at the interface was controlled with a partially threaded friction stir welding pin and a dovetail geometric configuration optimized in previous work. The process robustness was explored to establish parameter windows by elucidating parametric effects on the relationship between process conditions and the joint’s microstructure and mechanical properties. Ranges of process forge force (50–70 kN) and interface temperature (475°C–515°C) consistently yielded joints with adequate performance. The load-carrying capacity in lap shear tensile tests was 2025 ? 100 N/mm.

Friction Stir Dovetailing, FSW, AA7099, AA6061, Ro↗

DOE OSTI2020

Laminar drag reduction in microchannels with liquid infused textured surfaces

Liquid infused surfaces (LIS) achieve non-wetting properties through asperity structures containing pockets of a lubricating liquid rather than air. Although studies have demonstrated several potential applications of LIS, the effect of liquid infused texture on the flow of liquid in microchannels has not been systematically addressed. The present work focuses on the effect of liquid-infused texture geometry and the infused liquid properties on the flow characteristics of a bulk fluid in rectangular and round microchannel geometries. Considering four different textured geometries of transverse ridges, longitudinal ridges, in-lined square posts and staggered squares, a detailed parametric study is presented to assess the resulting drag reduction and friction factor, expressed in terms of the Poiseuille number, on the flow in microchannels. The effect of geometrical parameters is investigated for different asperity solid fraction (Φ s ), asperity height (h s ) and channel thickness (H) or tube diameter (D), while the effect of the infused fluid is studied in terms of a ratio of its viscosity relative to the viscosity of the flowing bulk fluid. Based on the results of the parametric study an extension of the conventional Moody diagram for non-wetting surfaces is developed. Finally, a design of experiments study is conducted to relate the drag reduction and Poiseuille number to the various governing parameters, using which conditions that maximize drag reduction are identified.

42 ENGINEERING↗

DOE OSTI2020

Parametric Design Study of a Power Electronics Package

The lifetime of a power electronics package is, to a large extent, determined by the reliability of its bonded interfaces under major loading conditions in an operational environment. Based on the application-level requirements of the package, bonded interface material is selected based on the results obtained from accelerated tests such as thermal cycling and power cycling. In addition to evaluating the reliability through accelerated tests, it is important to consider the impact of other component layers on the thermomechanical performance of the interface material, both from a material and geometric perspective. The co-efficient of thermal expansion (CTE) mismatch introduced by the use of different materials within a package and its structural design plays a critical role in determining the interface material reliability. In this paper, we present the results of a parametric modeling study of a power electronics package under thermal cycling in which the materials and geometric design of the different component layers were varied with respect to a baseline design to understand their impact on the reliability of the interface material. We chose the volume-averaged strain energy density per cycle computed at the corner region of the interface material as the metric for the reliability comparisons. Our results indicate that in addition to the CTE mismatch, the stiffness of individual component layers has a major impact on reliability. Among the different baseplates that we studied, aluminum silicon-carbide baseplates offered superior reliability over their copper and aluminum counterparts. We also found that the magnitude of the impact of stiffness variation—amongst the adjacent layers—on the reliability of the solder joint is dependent on the inherent CTE mismatch between the layers.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Parametric Design Study of a Power Electronics Package for Improving Solder Joint Reliability

47 OTHER INSTRUMENTATION↗

DOE OSTI2026

Biopolymer-Templated Titania Film Formation for Nanostructured Coatings Revealed by Machine Learning-Supported Time-Resolved Analysis

This study presents a machine learning approach to derive the film formation of biopolymer-templated titania nanostructures during spray deposition, in combination with in situ grazing-incidence small-angle X-ray scattering (GISAXS). A neural network trained on synthetic GISAXS data directly predicts domain-size distributions from experimental two-dimensional scattering patterns, capturing the full kinetics of nanostructure evolution with high temporal resolution. The predictions reveal hierarchical size distributions and periodic growth features, consistent with layer-by-layer spray deposition and validated by complementary scanning electron microscopy (SEM) imaging. Quantitative comparison with conventional parametric GISAXS fits shows good qualitative agreement, with systematic differences explained by domain-shape assumptions and resolved by applying a geometric scaling factor. Simulated SEM-like surfaces derived from neural network outputs reproduce the porous, foam-like nanoscale morphology observed experimentally, reinforcing the method’s credibility. This integrated approach enables real-time, nondestructive, statistically averaged monitoring of bulk nanostructure development in functional coatings, offering a scalable methodology to accelerate the characterization and process control of sustainably manufactured nanostructured titania films for energy-related applications such as photocatalysis and photovoltaics.

Heger, JulianEliah↗

DOE OSTI2024

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

Parametric study of optical transmission through plasmonic hole arrays modulated by the phase transition of vanadium dioxide

We have performed comprehensive electromagnetic simulations and preliminary experiments to explore the effects of geometrical and material parameters on the extraordinary optical transmission (EOT) through periodic arrays of subwavelength holes in a bilayer stack consisting of a gold or silver film atop a vanadium dioxide film (Au/Ag + VO 2 ), where the latter undergoes a semiconductor-to-metal phase transition. Using the finite-difference time-domain (FDTD) and finite-element methods (FEM), we vary iteratively the array periodicity, VO 2 film thickness and hole diameters, as well as the refractive index inside the VO 2 -layer holes and the VO 2 optical constants. For each variation, we compare the metallic-to-semiconducting ratios of the zero-order transmission ( T 00 ) peaks and find sharp maxima in these ratios within narrow parameter ranges. The maxima arise from Fabry-Perot and Fano-type resonances that minimize T 00 in the semiconducting phase of the perforated bilayers. At a fixed array period, the primary factors controlling the VO 2 -enabled EOT modulation are the VO 2 thickness, diameter of the VO 2 -layer holes, and absorption in the two VO 2 phases. Besides uncovering the origins of the higher metallic-phase T 00 , this study provides a protocol for optimizing the performance of the bilayer hole arrays for potential uses as dynamically tunable nano-optical devices.

Donev, Eugenii U. (ORCID:000000030555219X)↗

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↗