Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density modification”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

CALIPSO Lidar Calibration at 532 nm: Version 4 Nighttime Algorithm

Data products from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on board Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) were recently updated following the implementation of new (version 4) calibration algorithms for all of the level 1 attenuated backscatter measurements. In this work we present the motivation for and the implementation of the version 4 nighttime 532 nm parallel channel calibration. The nighttime 532 nm calibration is the most fundamental calibration of CALIOP data, since all of CALIOP’s other radiometric calibration procedures – i.e., the 532 nm daytime calibration and the 1064 nm calibrations during both nighttime and daytime – depend either directly or indirectly on the 532 nm nighttime calibration. The accuracy of the 532 nm nighttime calibration has been significantly improved by raising the molecular normalization altitude from 30-34 km to 36-39 km to substantially reduce stratospheric aerosol contamination. Due to the greatly reduced molecular number density and consequently reduced signal-to-noise ratio (SNR) at these higher altitudes, the signal is now averaged over a larger number of samples using data from multiple adjacent granules. As well, an enhanced strategy for filtering the radiation-induced noise from high energy particles was adopted. Further, the meteorological model used in the earlier versions has been replaced by the improved MERRA-2 model. An aerosol scattering ratio of 1.01 ± 0.01 is now explicitly used for the calibration altitude. These modifications lead to globally revised calibration coefficients which are, on average, 2-3% lower than in previous data releases. Further, the new calibration procedure is shown to eliminate biases at high altitudes that were present in earlier versions and consequently leads to an improved representation of stratospheric aerosols. Validation results using airborne lidar measurements are also presented. Biases relative to collocated measurements acquired by the Langley Research Center (LaRC) airborne high spectral resolution lidar (HSRL) are reduced from 3.6% ± 2.2% in the version 3 data set to 1.6% ± 2.4 % in the version 4 release.

Jayanta Kar↗
NASA NTRS2021

Tonlé Sap Food Security & Agriculture III: Evaluating Changes in Ecosystem Vitality and Freshwater Health in the Tonlé Sap Basin using Remotely Sensed Data

Tonlé Sap Lake, located in Cambodia, Southeast Asia, is one of the most productive inland fisheries in the world. With the unique reverse flow hydrology of the Tonlé Sap River, this freshwater system contains rich biodiversity and provides critical freshwater resources for the local community. Overfishing, stronger seasonality, drought, dam construction, forest fires, and untreated industrial domestic sewage threaten the ecosystem vitality and economic success of Tonlé Sap. In collaboration with Conservation International and the Cambodian Ministry of Water Resources and Meteorology’s Tonlé Sap Authority, we calibrated and finalized remotely-sensed proxies for sub-indicators of the Freshwater Health Index (FHI). We used NASA Earth observation data from Moderate Resolution Imaging Spectroradiometer (MODIS), and the Gravity Recovery and Climate Experiment (GRACE). These datasets were used in RS proxies and a Soil and Water Assessment Tool (SWAT) model that previous teams developed for sub-indicators of FHI. They included landcover, bank modification, and water quality metrics. It was determined that the ground water storage in Tonlé Sap Lake exhibited a slight declining trend over from April 1st, 2002 to February 3rd, 2017. To calculate chlorophyll-a concentration in the Tonlé Sap Lake, we used Copernicus Global Land Service mission data which provided a Trophic State Index of the lake. Between January 1st, 2000 and December 31st, 2020, the lake and surrounding region within the study area boundary displayed a slight decrease in vegetation density and consistently high chlorophyll-a concentrations. The SWAT model calculated nitrogen and phosphorus content measured in outlet points of the lake. We demonstrated that remotely sensed data is valuable for providing additional information for the FHI, but is not fully capable of replacing its in situ counterpart. We also concluded the water quantity is on a slowly declining trend within the basin.

Joey Lindsay↗
NASA NTRS1979

A fluid description of plasma double-layers

The space-charge double-layer that forms between two plasmas with different densities and thermal energies was investigated using three progressively realistic models which are treated by fluid theory, and take into account four species of particles: electrons and ions reflected by the double-layer, and electrons and ions transmitted through it. The two plasmas are assumed to be cold, and the self-consistent potential, electric field and space-charge distributions within the double-layer are determined. The effects of thermal velocities are taken into account for the reflected particles, and the modifications to the cold plasma solutions are established. Further modifications due to thermal velocities of the transmitted particles are examined. The applicability of a one dimensional fluid description, rather than plasma kinetic theory, is discussed. Theoretical predictions are compared with double layer potentials and lengths deduced from laboratory and space plasma experiments.

Levine, J. S.↗
NASA NTRS1986

Effects on optical systems from interactions with oxygen atoms in low earth orbits

Modifications of material surface properties due to interactions with ambient atomic oxygen have been observed on surfaces facing the orbital direction in low earth orbits. Some effects are very damaging to surface optical properties while some are more subtle and even beneficial. Most combustible materials are heavily etched, and some coatings, such as silver and osmium, are seriously degraded or removed as volatile oxides. The growth of oxide films on metals and semiconductors considered stable in dry air was measured. Material removal, surface roughness, reflectance, and optical densities are reported. Effects of temperature, contamination, and overcoatings are noted.

Peters, P. N.↗
NASA NTRS1983

A case study of the energy dissipation of the gravity wave field based on satellite altimeter measurements

Wave breaking is proposed as the primary energy dissipation mechanism for the gravity wave field. The energy dissipation rate is calculated based on the statistical model proposed by Longuet-Higgins (1969) with a modification of the breaking criterion incorporating the surface stress according to Phillips and Banner (1974). From this modified model, an analytic expression is found for the wave attenuation rate and the half-life time of the wave field which depend only on the significant slope of the wave field and the ratio of friction velocity to initial wave phase velocity. These expressions explain why the freshly generated wave field does not last long, but why swells are capable of propagating long distances without substantial change in energy density. It is shown that breaking is many orders of magnitude more effective in dissipating wave energy than the molecular viscosity, if the significant slope is higher than 0.01. Limited observational data from satellite and laboratory are used to compare with the analytic results, and show good agreement.

Huang, N. E.↗
NASA NTRS1986

The temperature structure in accretion flows onto massive protostars

Radiation transfer problems involved in the infall of dust and gas during star formation are studied. Dust properties are discussed, and modifications of spherical radiative transfer equations are presented that permit forward scattering by dust to be treated for the small size of the star relative to the inner radius of the shell. A procedure for deriving the stellar radiation field incident on the inner edge of the shell is developed. The temperature correction procedure of Cassinelli and Hartmann (1975) for extended stellar atmospheres is modified so that the multitemperature nature of the grains in the cloud may be derived. Temperature distributions for three schematic models in which the density is prespecified are discussed. Radiative acceleration of grains is addressed, showing that the proper mean opacity differs by a large factor from the Rosseland mean opacity that is commonly used. Emergent fluxes for the models are given.

Wolfire, Mark G.↗
NASA NTRS1989

Detailed design of a Ride Quality Augmentation System for commuter aircraft

The design of a Ride Quality Augmentation System (RQAS) for commuter aircraft is documented. The RQAS is designed for a Cessna 402B, an 8 passenger prop twin representative to this class of aircraft. The purpose of the RQAS is the reduction of vertical and lateral accelerations of the aircraft due to atmospheric turbulence by the application of active control. The detailed design of the hardware (the aircraft modifications, the Ride Quality Instrumentation System (RQIS), and the required computer software) is examined. The aircraft modifications, consisting of the dedicated control surfaces and the hydraulic actuation system, were designed at Cessna Aircraft by Kansas University-Flight Research Laboratory. The instrumentation system, which consist of the sensor package, the flight computer, a Data Acquisition System, and the pilot and test engineer control panels, was designed by NASA-Langley. The overall system design and the design of the software, both for flight control algorithms and ground system checkout are detailed. The system performance is predicted from linear simulation results and from power spectral densities of the aircraft response to a Dryden gust. The results indicate that both accelerations are possible.

Suikat, Reiner↗
NASA NTRS2019

Applications of Atomic Layer Deposition in the Modification of Carbon Nanotubes

Atomic Layer Deposition (ALD) is a cost-effective nanoadditive-manufacturing technique that allows for conformal coating of substrates with atomic control in a benign temperature and pressure environment. Using paired precursor gases, thin films can be deposited on flat or textured surfaces ranging from glass, polymers, aerogels, and metals. Through atomic layer control, where single layers of atoms are deposited, fabrication of metal transparent films, nano-laminates, and coatings of nano-channels and pores is achievable. Reaction mechanisms in ALD are normally self-limiting, allowing for atomically accurate control of nanometer (nm) thicknesses. Therefore, high uniformity and precise thickness control make ALD an attractive process for the creation of novel nano-scale devices. Decreases in resistivity and density of electrical wire are needed to improve the function of electronics, electric motors and cables. Such improvements may be accomplished by adding ballistically conducting, metallic carbon nanotubes (CNT) to Cu. ALD was used to coat multiple substrates including CNT with Cu in an effort to make CNT-Cu composites which is more conductive and less dense than Cu. In addition, ALD can be used to as a method to deposit the catalyst layer seed layer for CNT growth using Ni. The ALD of copper metallic films can follow multiple reaction pathways depending on the ALD precursors used. For this work the reaction pathway was deposition of copper oxide and then post process annealing in a hydrogen environment. Copper(II) diacetylacetonate (Cu(acac)2) and ozone are used as precursors for copper oxide. As-deposited copper oxide films prepared at 180°C resulted in a growth per cycle of 1.0 A/cycle with low film resistivity.

Dwivedi, Vivek H.↗
NASA NTRS1974

Electron capture in carbon dwarf supernovae

The rates of electron capture on heavier elements under the extreme conditions predicted for dwarf star supernovae have been computed, incorporating modifications that seem to be indicated by present experimental results. An estimate of the maximum possible value of such rates is also given. The distribution of nuclei in nuclear statistical equilibrium has been calculated for the range of expected supernovae conditions, including the effects of the temperature dependence of nuclear partition functions. These nuclide abundance distributions are then used to compute nuclear equilibrium thermodynamic properties. The effects of the electron capture on such equilibrium matter are discussed. In the context of the 'carbon detonation' supernova model, the dwarf central density required to ensure core collapse to a neutron star configuration is found to be slightly higher than that obtained by Bruenn (1972) with the electron capture rates of Hansen (1966).-

Mazurek, T. J.↗
NASA NTRS1991

Nonequilibrium ionization effects in asymmetrically heated loops

The effects of nonequilibrium ionization on magnetic loop models with a steady siphon flow that is driven by a nonuniform heating rate are investigated. The model developed by Mariska (1988) to explain the observed redshifts of transition region emission lines is examined, and the number densities of the ions of carbon and oxygen along the loop are computed, with and without the approximation of ionization equilibrium. Considerable deviations from equilibrium were found. In order to determine the consequences of these nonequilibrium effects on the characteristics of the EUV emission from the loop plasma, the profiles and wavelength positions of all the important emission lines due to carbon and oxygen were calculated. The calculations are in broad agreement with Mariska's conclusions, although they show a significant diminution of the Doppler shifts, as well as modifications to the line widths. It is concluded that the inclusion of nonequilibrium effects make it more difficult to reproduce the observed characteristics of the solar transition region by means of the asymmetric-heating models.

Spadaro, D.↗
NASA NTRS1994

Inertial currents in isotropic plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasmas, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MHD plasma. Solutions are developed by taking the MHD limit ot two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗
NASA NTRS1993

Inertial Currents in Isotropic Plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MED plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗
NASA NTRS1999

Structural Reliability Analysis and Optimization: Use of Approximations

This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different approximations including the higher-order reliability methods (HORM) for representing the failure surface. This report is divided into several parts to emphasize different segments of the structural reliability analysis and design. Broadly, it consists of mathematical foundations, methods and applications. Chapter I discusses the fundamental definitions of the probability theory, which are mostly available in standard text books. Probability density function descriptions relevant to this work are addressed. In Chapter 2, the concept and utility of function approximation are discussed for a general application in engineering analysis. Various forms of function representations and the latest developments in nonlinear adaptive approximations are presented with comparison studies. Research work accomplished in reliability analysis is presented in Chapter 3. First, the definition of safety index and most probable point of failure are introduced. Efficient ways of computing the safety index with a fewer number of iterations is emphasized. In chapter 4, the probability of failure prediction is presented using first-order, second-order and higher-order methods. System reliability methods are discussed in chapter 5. Chapter 6 presents optimization techniques for the modification and redistribution of structural sizes for improving the structural reliability. The report also contains several appendices on probability parameters.

Grandhi, Ramana V.↗
NASA NTRS1974

Three-dimensional tertiary structure of yeast phenylalanine transfer RNA

Results of an analysis and interpretation of a 3-A electron density map of yeast phenylalanine transfer RNA. Some earlier detailed assignments of nucleotide residues to electron density peaks are found to be in error, even though the overall tracing of the backbone conformation of yeast phenylalanine transfer RNA was generally correct. A new, more comprehensive interpretation is made which makes it possible to define the tertiary interactions in the molecule. The new interpretation makes it possible to visualize a number of tertiary interactions which not only explain the structural role of most of the bases which are constant in transfer RNAs, but also makes it possible to understand in a direct and simple fashion the chemical modification data on transfer RNA. In addition, this pattern of tertiary interactions provides a basis for understanding the general three-dimensional folding of all transfer RNA molecules.

Kim, S. H.↗
NASA NTRS1994

Comet coma sample return instrument

The sample collection technology and instrument concept for the Sample of Comet Coma Earth Return Mission (SOCCER) are described. The scientific goals of this Flyby Sample Return are to return to coma dust and volatile samples from a known comet source, which will permit accurate elemental and isotopic measurements for thousands of individual solid particles and volatiles, detailed analysis of the dust structure, morphology, and mineralogy of the intact samples, and identification of the biogenic elements or compounds in the solid and volatile samples. Having these intact samples, morphologic, petrographic, and phase structural features can be determined. Information on dust particle size, shape, and density can be ascertained by analyzing penetration holes and tracks in the capture medium. Time and spatial data of dust capture will provide understanding of the flux dynamics of the coma and the jets. Additional information will include the identification of cosmic ray tracks in the cometary grains, which can provide a particle's process history and perhaps even the age of the comet. The measurements will be made with the same equipment used for studying micrometeorites for decades past; hence, the results can be directly compared without extrapolation or modification. The data will provide a powerful and direct technique for comparing the cometary samples with all known types of meteorites and interplanetary dust. This sample collection system will provide the first sample return from a specifically identified primitive body and will allow, for the first time, a direct method of matching meteoritic materials captured on Earth with known parent bodies.

Albee, A. L.↗
NASA NTRS1994

Inertial currents in isotropic plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, magnetohyrodynamic (MHD) plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗
NASA NTRS2000

Case Study Analyses of the SUCCESS DC-8 Scanning Lidar Database

Under project SUCCESS (Subsonic Aircraft Contrail and Cloud Effects Special Study) funded by the Atmospheric Effects of Aviation Program, SRI International (SRI) developed an angular scanning backscatter lidar for operation on the NASA DC-8 research aircraft and deployed the scanning lidar during the SUCCESS field campaign. The primary purpose of the lidar was to generate real-time video displays of clouds and contrails above, ahead of, and below the DC-8 as a means to help position the aircraft for optimum cloud and contrail sampling by onboard in situ sensors, and to help extend the geometrical domain of the in situ sampling records. A large, relatively complex lidar database was collected and several data examples were processed to illustrate the value of the lidar data for interpreting the other data records collected during SUCCESS. These data examples were used to develop a journal publication for the special SUCCESS Geophysical Research Letters issue. The data examples justified data analyses of a larger part of the DC-8 lidar database and is the objective of the current study. Efficient processing of the SUCCESS DC-8 scanning lidar database required substantial effort to enhance hardware and software components of the data system that was used for the initial analyses. MATLAB instructions are used to generate altitude and distance color-coded lidar displays corrected for effects introduced by aircraft pitch and forward movement during an angular scan time interval. Onboard in situ sensor atmospheric measurements are propagated to distances ahead of the DC-8 using recorded aircraft velocity so that they can be plotted on the lidar displays for comparison with lidar remotely observed aerosol distributions. Resulting lidar and in situ sensor polar scan displays over extended sampling intervals are integrated into a time series movie format for 36 case studies. Contrails and clouds were detected to ranges of 15 km by the forward-viewing angular scanning lidar and were progressively mapped as the aircraft approached and penetrated them. Near aircraft lidar observations were much better correlated with in situ sensor observations than lidar observations at greater distances ahead of the aircraft. The major cause of this difference was thought to be the about 2 deg. offset of the lidar viewing direction from the flight direction. Contrail spatial distributions were not of the quality obtainable from ground-based lidar observations. This results because contrails tend to become horizontally stratified, vertical distance between angular lidar observations increases with increased distance from the aircraft, and erratic aircraft motions during an angular scan. The most useful lidar observations were made with lidar viewing directions of vertically upward or vertically downward. These provided real-time information on aircraft altitudes to achieve optimum in situ cloud and contrail sampling. At sampling altitudes, the forward viewing angular scanning observations were useful for fine-tuning the aircraft altitude for cloud and contrail penetration. Best information on cloud and contrail properties were obtained from vertically directed lidar observations as the aircraft performed a series of upward and downward penetrations of contrails. This operational mode was especially well suited for lidar and radiometric evaluation of cloud and contrail optical and radiative properties. The vertical viewing lidar detected ice crystals thought to be precipitating from an aircraft contrail and their scavenging by a cirrus cloud layer. The lidar display indicates that the crystals are effective for increasing cirrus cloud density. Vertical angular scanning observations can evaluate the sharp decrease in lidar backscatter for small off-vertical viewing directions that result from horizontally aligned ice crystals and perhaps can provide additional information on crystal shapes. The about 2 deg. offset of the lidar viewing direction from the flight direction is thought to have greatly degraded the forward-viewing angular scanning observations and this mode of operation was not fully evaluated. However, the reasoning for this capability remains valid and the angular scan presentations collected during this program justifies modification of the lidar pod for true forward direction lidar viewing during future cloud and contrail studies.

Uthe, Edward E.↗
NASA NTRS1977

Some optical and kinetic properties of the nearby interstellar gas

A model for resonant scattering of sunlight by the nearby interstellar gas is developed which takes into account the modification of a Maxwell-Boltzmann velocity distribution by a central force. While the formalism applies to both hydrogen and helium, the latter is emphasized. Theoretical line profiles and integrated emission rates for helium are compared with those of approximate models and with existing observations. As a result of this analysis, an upward revision of the bulk flow velocity of interstellar helium relative to the solar system, previously deduced from He 584-A observations, is required. Velocities of 15-20 km/s, temperatures of the order of 5000-10,000 K, and densities of the order of 0.008-0.02 per cu cm are sufficient to explain the (helium) observations. However, higher-precision measurements are required in order to place more stringent limits on the gas parameters. Calculations of the flux of interstellar helium atoms impinging on earth's upper atmosphere show that the resulting impact ionization is small compared with terrestrial processes.

Meier, R. R.↗