Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shape Matching”

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 73 records · Page 4

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Imaging the effect of dipole tilt on magetotail boundaries

Eight years of Interplanetary Monitoring Platform (IMP) 8, four years of International Sun Earth Explorer (ISEE) 2, and one year of IMP 7 magnetometer data have been combined to produce an 'image' of the average magnetic field for a YZ cross section (aberrated GSM) of the magnetotail at a downtail distance of 25 R(sub E). The shape of the neutral sheet and magnetopause boundaries can be observed directly from the images. A fitting function that qualitatively matches the observed boundary shape can then be chosen. This approach improves on previous fits to possible unsuitable functional forms specified independently of the data. In addition, as a refinement of previous studies, we have corrected for varying solar wind dynamic pressure and the effects of tail flaring. We find the magnetopause is displaced above the XY plane with increasing dipole tilt. The neutral sheet is found to curve slightly more than the model of Fairlfield (1980) during times of large dipole tilt and near the flanks appears to differ substantially from the neutral sheet shape given by the analytic model of Voigt (1984), the more recent neutral sheet model of Dandouras (1988) based on the Voigt model, and the semi-empirical model of Tsyganenko (1989).

Hammond, C. Max↗

Build-Up Approach to Updating the Mock Quiet Spike Beam Model

When a new aircraft is designed or a modification is done to an existing aircraft, the aeroelastic properties of the aircraft should be examined to ensure the aircraft is flight worthy. Evaluating the aeroelastic properties of a new or modified aircraft can include performing a variety of analyses, such as modal and flutter analyses. In order to produce accurate results from these analyses, it is imperative to work with finite element models (FEM) that have been validated by or correlated to ground vibration test (GVT) data, Updating an analytical model using measured data is a challenge in the area of structural dynamics. The analytical model update process encompasses a series of optimizations that match analytical frequencies and mode shapes to the measured modal characteristics of structure. In the past, the method used to update a model to test data was "trial and error." This is an inefficient method - running a modal analysis, comparing the analytical results to the GVT data, manually modifying one or more structural parameters (mass, CG, inertia, area, etc.), rerunning the analysis, and comparing the new analytical modal characteristics to the GVT modal data. If the match is close enough (close enough defined by analyst's updating requirements), then the updating process is completed. If the match does not meet updating-requirements, then the parameters are changed again and the process is repeated. Clearly, this manual optimization process is highly inefficient for large FEM's and/or a large number of structural parameters. NASA Dryden Flight Research Center (DFRC) has developed, in-house, a Mode Matching Code that automates the above-mentioned optimization process, DFRC's in-house Mode Matching Code reads mode shapes and frequencies acquired from GVT to create the target model. It also reads the current analytical model, as we11 as the design variables and their upper and lower limits. It performs a modal analysis on this model and modifies it to create an updated model that has similar mode shapes and frequencies as those of the target model. The Mode Matching Code output frequencies and modal assurance criteria (MAC) values that allow for the quantified comparison of the updated model versus the target model. A recent application of this code is the F453 supersonic flight testing platform, NASA DFRC possesses a modified F-15B that is used as a test bed aircraft for supersonic flight experiments. Traditionally, the finite element model of the test article is generated. A GVT is done on the test article ta validate and update its FEM. This FEM is then mated to the F-15B model, which was correlated to GVT data in fall of 2004, A GVT is conducted with the test article mated to the aircraft, and this mated F-15B/ test article FEM is correlated to this final GVT.

Herrera, Claudia Y.↗

A Three-Dimensional Coupled Internal/External Simulation of a Film-Cooled Turbine Vane

A three-dimensional Navier-Stokes simulation has been performed for a realistic film-cooled turbine vane using the LeRC-HT code. The simulation includes the flow regions inside the coolant plena and film cooling holes in addition to the external flow. The vane is the subject of an upcoming NASA Glenn Research Center experiment and has both circular cross-section and shaped film cooling holes. This complex geometry is modeled using a multi-block grid which accurately discretizes the actual vane geometry including shaped holes. The simulation matches operating conditions for the planned experiment and assumes periodicity in the spanwise direction on the scale of one pitch of the film cooling hole pattern. Two computations were performed for different isothermal wall temperatures, allowing independent determination of heat transfer coefficients and film effectiveness values. The results indicate separate localized regions of high heat transfer coefficient values, while the shaped holes provide a reduction in heat flux through both parameters. Hole exit data indicate rather simple skewed profiles for the round holes, but complex profiles for the shaped holes with mass fluxes skewed strongly toward their leading edges.

Heidmann, James D.↗

Increased Efficiency LED

In a semiconductor light emitting diode (LED) a large proportion of the light produced will be internally reflected when it strikes the end face of the LED at a semiconductor-air interface. This is due to the high index of refraction of most LED semiconductor materials. This problem may be partially overcome by modification of the shape of the LED by a reverse taper of the sides of the LED. Light is redirected by the taper to strike the interface at an angle closer to normal. This allows light to exit the LED that would be totally internally reflected in an untapered LED. The novelty of the present invention lies in the tapering of the sides of the LED for increased through transmission of light. Prior art devices have made use surface modifications such as roughening, addition of guiding layers, use of index matching materials and hemispherical shaping of the emitting end of the LED. Current technology LEDs have a transmission efficiency of only a few percent. An increase in this efficiency would be of great value to any technology that makes use of LEDs for light generation. The present invention is related to LAR 1 5050-SB, Collection of Light from Optical Fibers of NA greater than 1.

Egalon, Claudio Oliviera↗

Collisionless tearing in a field-reversed sheet pinch assuming nonparallel propagation

The problem of collisionless linear tearing is examined assuming a wave vector with a component normal to the equilibrium field. The geometry is defined and the general form of the linear dispersion equation is calculated. The linear theory results when k is parallel to B are reviewed, and Ampere's law is calculated for the external adiabatic region when k times B does not equal zero, using two-fluid theory. A solution is obtained for the approximate form of the perturbed currents and vector potential assuming quasi-parallel k. The resonant current contributions within the singular layer are calculated, obtaining an estimate of the dispersion equation. The form of the adiabatic currents within the singular layer is calculated, showing that an x-z current system persists even in the limit k perpendicular to B goes to zero. Finally, the perturbed vector potential solutions across the singular layer are matched to obtain the shape of the complete eigenfunction.

Quest, K. B.↗

The optical structure of the central core in the peculiar supernova remnant CTB 80

Deep and high resolution H-alpha + forbidden N II, O III, and red continuum interference-filter images of the bright central radio core of the supernova remnant CTB 80 are presented. These photographs reveal a surprisingly filamentary structure not visible from previous imaging and indicate considerable morphological differences between forbidden O III and H-alpha + forbidden N II emission regions. The H-alpha + forbidden N II filaments appear elongated in the E-W direction like that of the radio emission plateau surrounding the central core. The forbidden O III emission is much smaller in angular size but matches well the position, shape and size of the core's 6 cm radio emission. Also detected are the two possible optical counterparts reported by Blair and Schild (1984) to the remnants' X-ray point source and it is noted that both may have coincident H-alpha + forbidden N II emission.

Fesen, R. A.↗

Generic hypersonic vehicle performance model

An integrated computational model of a generic hypersonic vehicle was developed for the purpose of determining the vehicle's performance characteristics, which include the lift, drag, thrust, and moment acting on the vehicle at specified altitude, flight condition, and vehicular configuration. The lift, drag, thrust, and moment are developed for the body fixed coordinate system. These forces and moments arise from both aerodynamic and propulsive sources. SCRAMjet engine performance characteristics, such as fuel flow rate, can also be determined. The vehicle is assumed to be a lifting body with a single aerodynamic control surface. The body shape and control surface location are arbitrary and must be defined. The aerodynamics are calculated using either 2-dimensional Newtonian or modified Newtonian theory and approximate high-Mach-number Prandtl-Meyer expansion theory. Skin-friction drag was also accounted for. The skin-friction drag coefficient is a function of the freestream Mach number. The data for the skin-friction drag coefficient values were taken from NASA Technical Memorandum 102610. The modeling of the vehicle's SCRAMjet engine is based on quasi 1-dimensional gas dynamics for the engine diffuser, nozzle, and the combustor with heat addition. The engine has three variable inputs for control: the engine inlet diffuser area ratio, the total temperature rise through the combustor due to combustion of the fuel, and the engine internal expansion nozzle area ratio. The pressure distribution over the vehicle's lower aft body surface, which acts as an external nozzle, is calculated using a combination of quasi 1-dimensional gas dynamic theory and Newtonian or modified Newtonian theory. The exhaust plume shape is determined by matching the pressure inside the plume, calculated from the gas dynamic equations, with the freestream pressure, calculated from Newtonian or Modified Newtonian theory. In this manner, the pressure distribution along the vehicle after body expansion surface is then determined. The aerodynamic modeling, the engine modeling, and the exhaust plume analysis are described in more detail. A description of the computer code used to perform the above calculations is given and an input/output example is then given. The computer code is available on a Macintosh floppy disk.

Chavez, Frank R.↗

Optical, IUE, and ROSAT observations of the eclipsing nova-like variable V347 Puppis (LB 1800)

Using time-resolved optical spectroscopy and UBVRI and high-speed photometry obtained at Mount Stromlo Observatory, Mount John University Observatory, and the South African Astronomical Observatory; International Ultraviolet Explorer (IUE) ultraviolet spectroscopy; and Roentgen Satellite (ROSAT) survey X-ray fluxes, we present a study of the accretion disk, hot spot, and emission line regions in the bright eclipsing nova-like variable V347 Pup (LB 1800). In the optical and UV, V347 Pup is a strong emission line source with a continuum spectrum which is remarkably red for a high-M cataclysmic variable. Consistent with its high inclination, we interpret the continuum spectrum as the superposition of the spectrum of the cool (T(sub eff) approximately 7000 K) outer edge and the hot (T(sub eff) approximately 100,000 K) inner regions of a self-eclipsed accretion disk. For the assumed parameters, the model matches the level and shape of the observed spectrum for an inclination of approximately 88 and a distance of approximately 300 pc. The prominent hump in the optical and UV light curves just before eclipse manifests the presence of the hot spot where the accretion stream strikes the edge of the disk. The wavelength dependence of the amplitude of the hump is best modeled by a spot having an effective temperature of approximately 25,000 K and an area of approximately 3 x 10(exp 18) sq cm if the spot radiates like a blackbody, or an effective temperatue of approximately 14,000 K and an area of approximately 3 x 10(exp 19) sq cm if it radiates with a stellar spectrum. In either case, the hot spot produces only one-tenth of the predicted luminosity for the assumed mass-transfer rate of 10(exp -8) solar mass/yr. Either the hot spot is 'buried' in the edge of the accretion disk, or a significant fraction of its luminosity is radiated away in lines. The difference in azimuth between the peak of the hump and the dynamically expected location of the hot spot suggests that the spot's emitting surface is rotated forward by approximately 36 deg relative to the edge of the disk.

Mauche, Christopher W.↗

Excitation of solar p-modes

We investigate the rates at which energy is supplied to individual p-modes as a function of their frequencies nu and angular degrees l. The observationally determined rates are compared with those calculated on the hypothesis that the modes are stochastically excited by turbulent convection. The observationally determined excitation rate is assumed to be equal to the product of the mode's energy E and its (radian) line width Gamma. We obtain E from the mode's mean square surface velocity with the aid of its velocity eigenfuction. We assume that Gamma measures the mode's energy decay rate, even though quasi-elastic scattering may dominate true absorption. At fixed l, E(Gamma) arises as nu(exp 7) at low nu, reaches a peak at nu approximately equal 3.5 mHz, and then declines as nu(exp 4.4) at higher nu . At fixed nu, E(Gamma) exhibits a slow decline with increasing l. To calculate energy input rates, P(sub alpha), we rely on the mixing-length model of turbulent convection. We find entropy fluctuations to be about an order of magnitude more effective than the Reynolds stress in exciting p-modes . The calculated P(sub alpha) mimic the nu(exp 7) dependence of E(Gamma) at low nu and the nu(exp -4.4) dependence at high nu. The break of 11.4 powers in the nu-dependence of E(Gamma) across its peak is attributed to a combination of (1) the reflection of high-frequency acoustic waves just below the photosphere where the scale height drops precipitously and (2) the absence of energy-bearing eddies with short enough correlation times to excite high-frequency modes. Two parameters associated with the eddy correlation time are required to match the location and shape of the break. The appropriate values of these parameters, while not unnatural, are poorly constrained by theory. The calculated P(sub alpha) can also be made to fit the magnitude of E(Gamma) with a reasonable value for the eddy aspect ratio. Our resutls suggest a possible explanation for the decline of mode energy with increasing l at fixed nu. Entropy fluctuations couple to changes in volume associated with the oscillation mode. These decrease with decreasing n at fixed nu, becoming almost zero for the f-mode.

Goldreich, Peter↗

On Flowfield Periodicity in the NASA Transonic Flutter Cascade: Experimental Study - Part 1

An extensive study to improve flow uniformity and periodicity in the NASA Transonic Flutter Cascade is presented here. The results are reported in two independent parts dealing with the experimental approach and the analytical approach. The first part, the Experimental Study, focuses first on the data sets acquired in this facility in the past and explains several discrepancies, particularly the questions of actual flow incidence and cascade back pressure levels. Next, available means for control and modifications of the cascade flowfield, boundary layer bleed and tailboard settings are presented in detail. This is followed by experimental data sets acquired in modified test facility configurations that were based on analytical predictions of the cascade flowfield. Finally, several important conclusions about improving the cascade flowfield uniformity and blade load periodicity are summarized. The important conclusions are: (1) boundary layer bleed does not improve the cascade flow periodicity; (2) tunnel wall contours must be carefully matched to the expected shape of cascade streamlines; (3) actual flow incidence for each cascade configuration rather must be measured instead of relying on the tunnel geometry; and (4) the current cascade configuration exhibits a very high blade load uniformity over six blades from blade #2 to blade #7, and the facility is now ready for unsteady pressure data acquisition.

Lepicovsky, J.↗

The Micromechanics of the Moving Contact Line

A transient moving contact line is investigated experimentally. The dynamic interface shape between 20 and 800 microns from the contact line is compared with theory. A novel experiment is devised, in which the contact line is set into motion by electrically altering the solid-liquid surface tension gamma(sub SL). The contact line motion simulates that of spontaneous wetting along a vertical plate with a maximum capillary number Ca approx. = 4 x 10(exp -2). The images of the dynamic meniscus are analyzed as a funtion of Ca. For comparison, the steady-state hydrodynamic equation based on the creeping flow model in a wedge geometry and the three-region uniform perturbation expansion of Cox (1986) is adopted. The interface shape is well depicted by the uniform solutions for Ca <= 10(exp -3). However, for Ca > 10(exp -3), the uniform solution over-predicts the viscous bending. This over-prediction can be accounted for by modifying the slip coefficient within the intermediate solution. With this correction, the measured interface shape is seen to match the theoretical prediction for all capillary numbers. The amount of slip needed to fit the measurements does not scale with the capillary number.

Lichter, Seth↗

Simulations of inducers at Low-Flow Off-Design Conditions

The ability to accurately model details of inlet back flow for inducers operating at low-flow, off-design conditions is evaluated. A sub-scale version of a three- bladed liquid hydrogen inducer tested in water with detailed velocity and pressure measurements is used as a numerical test bed. Under low-flow, off-design conditions the length of the separation zone as well as the swirl velocity magnitude was under predicted with a standard k-E model. When the turbulent viscosity coefficient was reduced good comparison was obtained at all the flow conditions examined with both the magnitude and shape of the profile matching well with the experimental data taken half a diameter upstream of the leading edge. The velocity profiles and incidence angles at the leading edge itself were less sensitive to the back flow length predictions indicating that single-phase performance predictions may be well predicted even if the details of flow separation modeled are incorrect. However, for cavitating flow situations the prediction of the correct swirl in the back flow and the pressure depression in the core becomes critical since it leads to vapor formation. The simulations have been performed using the CRUNCH CFD@ code that has a generalized multi-element unstructured framework and an advanced multi-phase formulation for cryogenic fluids. The framework has been validated rigorously for predictions of temperature and pressure depression in cryogenic fluid cavities and has also been shown to predict the cavitation breakdown point for inducers at design conditions.

Hosangadi, Ashvin↗

Simulations of a Liquid Hydrogen Inducer at Low-Flow Off-Design Flow Conditions

The ability to accurately model details of inlet back flow for inducers operating a t low-flow, off-design conditions is evaluated. A sub-scale version of a three-bladed liquid hydrogen inducer tested in water with detailed velocity and pressure measurements is used as a numerical test bed. Under low-flow, off-design conditions the length of the separation zone as well as the swirl velocity magnitude was under predicted with a standard k-E model. When the turbulent viscosity coefficient was reduced good comparison was obtained a t all the flow conditions examined with both the magnitude and shape of the profile matching well with the experimental data taken half a diameter upstream of the leading edge. The velocity profiles and incidence angles a t the leading edge itself were less sensitive to the back flow length predictions indicating that single-phase performance predictions may be well predicted even if the details of flow separation modeled are incorrect. However, for cavitating flow situations the prediction of the correct swirl in the back flow and the pressure depression in the core becomes critical since it leads to vapor formation. The simulations have been performed using the CRUNCH CFD(Registered Trademark) code that has a generalized multi-element unstructured framework and a n advanced multi-phase formulation for cryogenic fluids. The framework has been validated rigorously for predictions of temperature and pressure depression in cryogenic fluid cavities and has also been shown to predict the cavitation breakdown point for inducers a t design conditions.

Hosangadi, A.↗

Advanced Colloids Experiment (ACE) Science Overview

The Advanced Colloids Experiment is being conducted on the International Space Station (ISS) using the Light Microscopy Module (LMM) in the Fluids Integrated Rack (FIR). Work to date will be discussed and future plans and opportunities will be highlighted. The LMM is a microscope facility designed to allow scientists to process, manipulate, and characterize colloidal samples in micro-gravity where the absence of gravitational settling and particle jamming enables scientists to study such things as:a.The role that disordered and ordered-packing of spheres play in the phase diagram and equation of state of hard sphere systems,b.crystal nucleation and growth, growth instabilities, and the glass transition, c.gelation and phase separation of colloid polymer mixtures,d.crystallization of colloidal binary alloys,e.competition between crystallization and phase separation,f.effects of anisotropy and specific interactions on packing, aggregation, frustration and crystallization,g.effects of specific reversible and irreversible interactions mediated in the first case by hybridization of complementary DNA strands attached to separate colloidal particles,h.Lock and key interactions between colloids with dimples and spheres which match the size and shape of the dimples,i.finding the phase diagrams of isotropic and interacting particles,j.new techniques for complex self-assembly including scenarios for self-replication, k.critical Casimir forces,l.biology (real and model systems) in microgravity,m.etc. By adding additional microscopy capabilities to the existing LMM, NASA will increase the tools available for scientists that fly experiments on the ISS enabling scientists to observe directly what is happening at the particle level. Presently, theories are needed to bridge the gap between what is being observed (at a macroscopic level when photographing samples) with what is happening at a particle (or microscopic) level. What is happening at a microscopic level will be directly accessible with the availability of the Light Microscopy Module (LMM) on ISS. To meet these goals, the ACE experiment is being built-up in stages, with the availability of confocal microscopy being the ultimate objective. Supported by NASAs Physical Sciences Research Program, ESAESTEC, and the authors respective governments.

Microgravity Biology Research↗

The JWST/NIRCam Coronagraph: Mask Design and Fabrication

The NIRCam instrument on the James Webb Space Telescope will provide coronagraphic imaging from lambda =1-5 microns of high contrast sources such as extrasolar planets and circumstellar disks. A Lyot coronagraph with a variety of circular and wedge-shaped occulting masks and matching Lyot pupil stops will be implemented. The occulters approximate grayscale transmission profiles using halftone binary patterns comprising wavelength-sized metal dots on anti-reflection coated sapphire substrates. The mask patterns are being created in the Micro Devices Laboratory at the Jet Propulsion Laboratory using electron beam lithography. Samples of these occulters have been successfully evaluated in a coronagraphic testbed. In a separate process, the complex apertures that form the Lyot stops will be deposited onto optical wedges. The NIRCam coronagraph flight components are expected to be completed this year.

Coronagraph↗

Acceptable Tolerances for Matching Icing Similarity Parameters in Scaling Applications

This paper reviews past work and presents new data to evaluate how changes in similarity parameters affect ice shapes and how closely scale values of the parameters should match reference values. Experimental ice shapes presented are from tests by various researchers in the NASA Glenn Icing Research Tunnel. The parameters reviewed are the modified inertia parameter (which determines the stagnation collection efficiency), accumulation parameter, freezing fraction, Reynolds number, and Weber number. It was demonstrated that a good match of scale and reference ice shapes could sometimes be achieved even when values of the modified inertia parameter did not match precisely. Consequently, there can be some flexibility in setting scale droplet size, which is the test condition determined from the modified inertia parameter. A recommended guideline is that the modified inertia parameter be chosen so that the scale stagnation collection efficiency is within 10 percent of the reference value. The scale accumulation parameter and freezing fraction should also be within 10 percent of their reference values. The Weber number based on droplet size and water properties appears to be a more important scaling parameter than one based on model size and air properties. Scale values of both the Reynolds and Weber numbers need to be in the range of 60 to 160 percent of the corresponding reference values. The effects of variations in other similarity parameters have yet to be established.

Anderson, David N.↗

Experimental performance of a microwave cavity plasma disk ion source

The detailed description and performance characteristics of a microwave ion source are presented. This ion source utilizes an internally tuned, single-mode (or selective multimode) cylindrical cavity applicator to focus and match microwave energy into a disk-shaped discharge zone. The combination of mode focus control and variable, internal cavity matching allows the efficient operation and beam extraction over a wide range of pressures, powers and gaseous inputs. Experimental measurements of ion beam current versus accelerating voltage and input microwave power in xenon and oxygen gas are presented. Ion source specific energy and mass utilization versus experimental variables are also determined. The experimental performance demonstrates the ability of this ion source to extract an ion beam with a well matched, stable, and continuous operation over a wide range of input gases, low pressures, and over input gas flow rates in excess of 100 to a few sccm. Double Langmuir probe measurements in xenon gas indicate high degrees of ionization, and electron and ion concentrations in excess of 100 critical densities in the microwave discharge zone. This ion source has many potential uses such as spacecraft electric propulsion, material ion beam processing, and neutral beam ion sources.

Root, J.↗