Search

The recently developed collisional radiative (CR) models for neutral xenon (Xe I) and singly charged xenon ion (Xe II) described by Chaplin et al. in the paper "Insights from Collisional-Radiative Models of Neutral and Singly-Ionized Xenon in Hall Thrusters" also submitted to the 2021 AIAA Propulsion and Energy Forum are more thoroughly investigated. Comparisons between the existing and recently developed JPL Xe I CR models are done and reasons for discrepancies between the two models are explored. An experiment to validate the Xe II CR model is carried out in the anode region of the UCLA Plasma interactions facility. The Xe II CR model is applied to measured xenon ion emission line intensities to determine electron temperature and density. Spatial resolution of these optical measurements is obtained through an Abel inversion for direct comparison to Langmuir probe results. Although the measured electron densities were on the order of $10^{17}$ m$^{-3}$ and the Xe II CR model is most suitable for use with densities above $10^{18}$ m$^{-3}$ due to the strong dependence of the line ratios on electron temperature in that regime, initial results from the validation experiment show reasonable agreements between the optical and Langmuir probe results. The 441 nm to 605 nm line ratio appears to have the most promise for non-intrusive electron temperature measurements with the Xe II CR model. Shortcomings of the Xe II CR model revealed in the validation experiment are to be explored further.

Wear measurements of the inner pole cover in HERMeS have shown as much as ten times higher erosion rates than past numerical simulations predicted. Also, measurements of ion velocity distribution functions (IVDF) suggest that ion heating occurs near the chamfered regions of the acceleration channel and around the inner pole surface. In our companion paper we present theoretical evidence that the plasma conditions in this region can excite instabilities in the lower hybrid frequency. Here, we provide further evidence of their presence through comparisons between numerical simulations and measurements. After including in the 2-D axisymmetric hydrodynamic-PIC code Hall2De an idealized model of the ion temperature based on a scaling with the lower hybrid frequency, we find the agreement between the computed and measured IVDF is significantly improved. In addition, the computed erosion rates at the inner pole surface are now in close agreement with wear test measurements.

Excessive electron back streaming is one of the primary life-limiting failure modes for gridded ion thrusters. Physics-based modeling of the optics grid erosion and electron back streaming margins is augmented with statistical uncertainty quantification techniques to generate a life expectancy distribution of this particular failure mechanism for the NEXT ion thruster. Generation of distributions instead of a single point estimate provide a more comprehensive picture of thruster failure probabilities and life expectation for various mission applications.

The Psyche mission will employ Stationary Plasma Thrusters (SPT-140) for primary propulsion in space and will be the first to use Hall thrusters beyond lunar orbit. One of the longest standing risks in the implementation of this EP technology on missions like Psyche is our limited understanding of the effects of ground test facilities on the performance and wear of the thruster.

Numerical simulations have been conducted to assess the erosion rates at the pole covers of the magnetically shielded miniature (MaSMi) Hall thruster. MaSMi is part of the Ascendant Sub-kW Transcelestial Electric Propulsion System (ASTRAEUS), whose design requirements demand a xenon throughput of 100 kg, thus enabling deep space missions using SmallSats. The results of our simulations show that the proposed thickness of the graphite pole covers will not be completely eroded at end of life. We find that the erosion rates are larger near the corners of the covers, which are directly exposed to the plasma in the acceleration channel.

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

UNKNOWN

UNKNOWN

Accelerator (accel) grid sputtering by ions formed through charge-exchange (CEX) reactions between beam ions and residual neutral gas is a critical life-limiting mechanism for gridded ion thrusters. The three-dimensional ion optics code CEX3D is designed to simulate this grid erosion for a single beamlet, with a particular emphasis on non-axisymmetric features such as the "pits and grooves" erosion commonly observed on the accel grid downstream face in two-grid thrusters. The treatment of CEX ions in the code was recently upgraded with a new particle-in-cell (PIC) module to account for the influence of these ions' space charge on the electrostatic potential downstream of the grids. In order to achieve reasonable computation times while resolving the Debye length near the grids and avoiding gross violations of the Courant-Friedrichs-Lewy (CFL) condition, macroparticle velocities in the PIC calculation are limited through a rescaling procedure that preserves ion trajectories and space charge density. The code accounts for beam divergence, finite momentum transfer in CEX collisions, and radial losses of CEX ions from the beam; these effects are important for determining the CEX ion flux to the accel grid because the calculated potential downstream of the grids can become very flat. The upgraded code has been used to simulate operation of NASA's Evolutionary Xenon Thruster (NEXT) during the 51 kHr Long Duration Test - a selection of results is presented and compared with experimental data.

We have developed a physics-based model based on a pseudo-particle description of the electron cyclotron drift instability. A key improvement of the model with respect to previous work is that linear theory is not applied in the event of wave saturation and deviations of electrons or ions from a Maxwellian distribution function. In the acceleration region, the anomalous collision frequency is computed as the minimum value necessary to prevent the electron drift velocity from exceeding the thermal velocity. A functional based on the electron equilibration time is defined to control the transition from high to low resistivity regions. The model was previously applied to a single Hall thruster at its nominal operating condition, showing promising results that captured accurately the location of the thruster’s acceleration region. In this paper, we extend the use of this firstprinciples models to two additional thrusters, also considering multiple operating conditions for each of them. Numerical results are compared to experimental measurements obtained with non-invasive laser induced fluorescence. In general, the agreement between experiments and simulations is good. The model is able to predict the location of the acceleration region for all cases. We observe however that fine details, such as changes in the plasma potential gradient within the acceleration regions, are not captured. The model is also insensitive to changes in the magnetic field strength while experiments show that small shifts in location (of less than 5% of the acceleration channel length) occur. We plan to address the weaknesses of our method with the help of physical insight gained from kinetic simulations of the acceleration region.

Azimuthal direction 1-D and azimuthal and axial direction 2-D Particle in Cell (PIC) simulations of the Hall thruster acceleration and near-plume regions are presented. The axial extent of the 2-D computational domain is based on the ExB drift parameter, 𝑫𝑬𝒙𝑩, the number of times around the channel an electron ExB drifts before being scattered. With only high energy main beam ions, the simulations produced a broad potential profile in the acceleration region. The addition of a small amount of ionization in the acceleration domain at a rate consistent with Hall thruster parameters causes the potential profile to develop two distinct regions. In the upstream portion of the computational domain, the potential profile is very steep, typical of Hall thruster acceleration regions. In the downstream portion, the potential profile is almost flat, like potentials in Hall thruster near plume regions. This bifurcation happens naturally; there is no other change in code parameters between the two results. Without the addition of a small amount of ionization in the acceleration region, the axial electron transport is very low. With ionization, the calculated ratio of electron current to ion current is the same order of those measured in Hall thrusters. This suggests that in regions with lower energy electrons and ions, plasma waves are effective in allowing crossfield transport, and axial electric fields are small. In the acceleration region, the plasma waves are less effective at scattering the energetic electrons, and high axial electric fields are needed to produce the observed cross-field electron transport.

We present a computational model based on the hypothesis that the higher-than-expected electron mobility across magnetic field lines observed in Hall thrusters is due to the growth and eventual saturation of the electron cyclotron drift instability. The key improvement of this model with respect to previous work is that it decouples the saturation of the instability from the correlation that exists between mobility and magnitude of the wave perturbations in linear theory, thereby enabling first-principles simulations. The model has been incorporated in the 2-D (r-z) multi-fluids code Hall2De and simulations of the H6 thruster at 300 V, 20 A are numerically stable and achieve a steady-state solution at a computational cost that is not significantly higher than previous Hall2De simulations. The computed location of the acceleration region is found to be within 10% of the length of the acceleration channel compared to that inferred by experiments. The simulations also capture well the plasma gradients along the channel centerline of this thruster. To further establish the validity of this model, we plan to conduct simulations of other thrusters and/or operating conditions.

Several recent studies have shown that a shift of the red-edge inflection near 0.70 micron in vegetation reflectance spectra is an indicator of metal stress, partially attributable to changes in chlorophyll concentration. This 'red-edge shift', however, is difficult to detect and has been reported both toward longer (red) and shorter (blue) wavelengths. Our work demonstrates that direct measurement of the depth and width of the chlorophyll absorption band at 0.67 micron using digital feature extraction and absorption band characterization procedures developed for the analysis of mineral spectra is a more consistent indicator of metal stress. Additionally, the magnitude of these parameters is generally greater than that of the red edge shift and thus should be more amenable to detection and mapping using field and aircraft spectrometers.

The results of studies that have calculated projections of food deficits for less developed countries (LDC's) toward the year 2000 are compared. Included in the comparison are: (1) reports entitled "The Global 2000 Report to the President,""The Future of the World Economy," and "The Model of International Relations in Agriculture"; (2) work done at Iowa State University, the International Food Policy Research Institute, and the Food and Agriculture Organization of the United Nations; and (3) work sponsored by the Club of Rome. Almost without exception, these studies forsee a future for the food economies of LDC's characterized by rising deficits of consumption over production and/or increasing real food prices. These projects, most of which were initiated during the food crisis of 1973/74, are found to be at variance with observed trends in the post-war era and events occurring since the time they were conducted. A preponderance of attention on supply side factors, at the expense of a thorough investigation of transitory demand side influences, was characteristic of most of the studies. This incomplete diagnosis of the causes of the crisis contributed to a bias in the projections. (Author/RM)

To see what effects practice--frequent writing without teacher correction--would have on the writing of students in grades one and two in one small elementary school, writing samples were collected from all pupils at each grade level in October and May. Five factors were examined: general impression, length, vocabulary, spelling, and syntactic development. The students' teachers admitted that they had not been providing frequent opportunities for their pupils to write, but were willing to try this approach. Results showed significant improvement at both grade levels in all categories measured. The most interesting result was the startling variation in the amount of writing that individual children produced from day to day. It varied from only a few words one day to a hundred or more the next. The results indicate that, contrary to popular belief, many first grade students can write well at the beginning of the school year. Since children in first grade are already writing and are all learning to write correctly, the role of the teacher is not to teach beginning writing; rather it is to nurture beginning writers. (Tables of findings and writing samples are included.) (HOD)

Research program provides better understanding of pyrotechnic phenomenon for design purposes. Evaluating potential for damage to spacecraft by activation of pyrotechnic mechanisms, pyrotechnic-shock tests conducted on three configurations: pin pullers on orthogonal double Hopkinson bar arrangement; pin pullers on mockup of Halogen Occultation Experiment (HALOE) structure; and section of separation joint on single Hopkinson bar. Strains and accelerations measured. Strains converted to output stresses, forces, and moments. Acceleration shock-response spectra obtained for both acceleration and force signals. Results of research useful to designers in making comparison and evaluation tests before committing to costly spacecraft hardware.

To evaluate the potential for spacecraft damage caused by activation of pyrotechnic mechanisms, pyrotechnic shock tests were conducted on three configurations: pin-pullers on an orthogonal double Hopkinson bar arrangement, pin-pullers on a mock-up of the HALOE (Halogen Occultation Experiment) structure, and a section of separation joint on a single Hopkinson bar. Strains and accelerations were measured. The strains were converted to output stresses, forces and moments. Acceleration shock response spectra are presented for both acceleration- and force-time signals. The devices were attached to the Hopkinson bars with adaptors, which are typical of attachments used in practice. In order to predict the effect of the adaptors, finite element models of the adaptor-bar combinations were analyzed and results are presented with regard to grid size, time interval, and frequency range required to predict the response to pulses of 10 and 100 microsecond duration.

Two of the 43 Shuttle Imaging Radar (SIR-B) experiments carried out from the 41-G shuttle mission in 1984 involved a 2600-km swath across the Canadian Shield, with the objectives of studying the structure of province boundaries and developing techniques for the geologic use of orbital radar. Despite degraded single incidence angle imagery resulting from system problems, valuable experience has been obtained with data over a test site near Bancroft, Ontario. It has been found that even subdued glaciated topography can be effectively imaged, variations in backscatter being caused by variations in local incidence angle rather than shadowing. It has been demonstrated that small incidence angles are more sensitive to topography than large angles. Backscatter is extremely sensitive to look direction, topographic features nearly normal to the illumination being highlighted, and those nearly parallel to it being suppressed. It is concluded that orbital radar can provide a valuable tool for geologic studies of the Canadian Shield and similar areas, if suitable look angles and at least two look directions can be utilized for each area.

The effects of ice formations on the section lift, drag, and pitching-moment coefficients of an unswept NACA 65A004 airfoil section of 6-foot chord were studied.. The magnitude of the aerodynamic penalties was primarily a function of the shape and size of the ice formation near the leading edge of the airfoil. The exact size and shape of the ice formations were determined photographically and found to be complex functions of the operating and icing conditions. In general, icing of the airfoil at angles of attack less than 40 caused large increases in section drag coefficients (as much as 350 percent in 8 minutes of heavy glaze icing), reductions in section lift coefficients (up to 13 percent), and changes in the pitching-moment coefficient from diving toward climbing moments. At angles of attack greater than 40 the aerodynamic characteristics depended mainly on the ice type. The section drag coefficients generally were reduced by the addition of rime ice (by as much as 45 percent in 8 minutes of icing). In glaze icing, however, the drag increased at these angles of attack. The section lift coefficients were variably affected by rime-ice formations; however, in glaze icing, lift increases at high angles of attack amounted to as much as 9 percent for an icing time of 8 minutes. Pitching-moment-coefficient changes in icing conditions were somewhat erratic and depended on the icing condition. Rotation of the iced airfoil to angles of attack other than that at which icing occurred caused sufficiently large changes in the pitching-moment coefficient that, in flight, rapid corrections in trim might be required in order to avoid a hazardous situation.

An empirical relation has been obtained by which the change in drag coefficient caused by ice formations on an unswept NACA 65AO04 airfoil section can be determined from the following icing and operating conditions: icing time, airspeed, air total temperature, liquid-water content, cloud droplet impingement efficiencies, airfoil chord length, and angles of attack. The correlation was obtained by use of measured ice heights and ice angles. These measurements were obtained from a variety of ice formations, which were carefully photographed, cross-sectioned, and weighed. Ice weights increased at a constant rate with icing time in a rime icing condition and at progressively increasing rates in glaze icing conditions. Initial rates of ice collection agreed reasonably well with values predicted from droplet impingement data. Experimental droplet impingement rates obtained on this airfoil section agreed with previous theoretical calculations for angles of attack of 40 or less. Disagreement at higher angles of attack was attributed to flow separation from the upper surface of the experimental airfoil model.

Heating requirements for satisfactory cyclic de-icing over a wide range of icing and operating conditions have been determined for a gas-heated, 36deg swept airfoil of 6.9-foot chord with a partial-span leading-edge slat. Comparisons of heating requirements and effectiveness were made between the slatted and unslatted portions of the airfoil. Studies were also made comparing cyclic de-icing with continuous anti-icing, and cycll.cde-icing systems with and without leading-edge ice-free parting strips. De-icing heat requirements were approximately the same with either heated or unheated parting strips because of the aerodynamic effects of the 36deg sweep angle and the spanwise saw-tooth profile of leading-edge glaze-ice deposits. Cyclic de-icing heat-source requirements were found to be one-fourth or less of the heat requirements for complete anti-icing. The primary factors that affected the performance of the cyclic de-icing heating system were ambient air temperature, heat distribution, and thermal lag.

Studies were made to determine the effect of ice formations on the section drag of a 6.9-foot-chord 36 degree swept NACA 63A-009 airfoil with partial-span leading-edge slat. In general, the icing of a thin swept airfoil will result in greater aerodynamic penalties than for a thick unswept airfoil. Glaze-ice formations at the leading edge of the airfoil caused large increases in section drag even at liquid-water content of 0.39 gram per cubic meter. The use of an ice-free parting strip in the stagnation region caused a negligible change in drag compared with a completely unheated airfoil. Cyclic de-icing when properly applied caused the drag to decrease almost to the bare-airfoil drag value.

The effects of primary and runback ice formations on the section drag of a 36 deg swept NACA 63A-009 airfoil section with a partial-span leading-edge slat were studied over a range of angles of attack from 2 to 8 deg and airspeeds up to 260 miles per hour for icing conditions with liquid-water contents ranging from 0.39 to 1.23 grams per cubic meter and datum air temperatures from 10 to 25 F. The results with slat retracted showed that glaze-ice formations caused large and rapid increases in section drag coefficient and that the rate of change in section drag coefficient for the swept 63A-009 airfoil was about 2-1 times that for an unswept 651-212 airfoil. Removal of the primary ice formations by cyclic de-icing caused the drag to return almost to the bare-airfoil drag value. A comprehensive study of the slat icing and de-icing characteristics was prevented by limitations of the heating system and wake interference caused by the slat tracks and hot-gas supply duct to the slat. In general, the studies showed that icing on a thin swept airfoil will result in more detrimental aerodynamic characteristics than on a thick unswept airfoil.

Several methods of cyclic de-icing of a gas-heated airfoil were investigated to determine ice-removal characteristics and heating requirements. The cyclic de-icing system with a spanwise ice-free parting strip in the stagnation region and a constant-temperature gas-supply duct gave the quickest and most reliable ice removal. Heating requirements for the several methods of cyclic de-icing are compared, and the savings over continuous ice prevention are shown. Data are presented to show the relation of surface temperature, rate of surface heating, and heating time to the removal of ice.

Investigations have been made in flight and in wind tunnels to determine which components of turbojet installations are most critical in icing conditions, and to evaluate several methods of icing protection. From these studies, the requirements necessary for adequate icing protection and the consequent penalties on engine performance can be estimated. Because investigations have indicated that the compressor-inlet screen constitutes the greatest icing hazard and is difficult to protect, complete removal or retraction of the screen upon encountering an icing condition is recommended. In the absence of the screen, the inlet guide vanes of an axial-flow-type turbojet engine constitute the greatest danger to engine operation in an icing condition; a centrifugal-type engine, on the other hand, is relatively unsusceptible to icing once the screen has been removed. Of the three icing-protection systems investigated, surface heating, hot-gas bleedback, and inertia-separation inlets, only the first two offer an acceptable solution to the problem of engine icing protection. Surface heating, either by gas heating or electrical means, appears to be the most acceptable icing-protection method with regard to performance losses. Hot-gas bleedback, although causing undesirable thrust losses, offers an easy means of obtaining icing protection for some installations. The final choice of an icing-protection system depends, however, on the supply of heated gas and electrical power available and on the allowable performance and. weight penalties associated with each system.

A two-dimensional inlet-guide-vane cascade was investigated to determine the effects of ice formations on the pressure losses across the guide vanes and to evaluate the heated gas flow and temperature required to prevent Icing at various conditions. A gas flow of approximately 0.4 percent of the inlet-air flow was necessary for anti-icing a hollow guide-vane stage at an inlet-gas temperature of 500 F under the following icing conditions: air velocity, 280 miles per hour; water content, 0.9 gram per cubic meter; and Inlet-air static temperature, 00 F. Also presented are the anti-icing gas flows required with modifications of the hollow Internal gas passage, which show heatinput savings greater than 50 percent.

The effect of modifying the gas passage of hollow metal airfoils by the additIon of internal fins and partitions was experimentally investigated and comparisons were made among a basic unfinned airfoil section and two airfoil designs having metal fins attached at the leading edge of the internal gas passage. An analysis considering the effects of heat conduction in the airfoil metal was made to determine the internal modification effectiveness that may be obtained in gas-heated components, such as turbojet-inlet guide vanes, support struts, hollow propeller blades, arid. thin wings. Over a wide range of heated-gas flow and tunnel-air velocity, the increase In surface-heating rates with internal finning was marked (up to 3.5 times), with the greatest increase occurring at the leading edge where anti-icing heat requirements are most critical. Variations in the amount and the location of internal finning and. partitioning provided. control over the local rates of surface heat transfer and permitted efficient anti-icing utilization of the gas-stream heat content.