Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theoretical models”

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 415 records · Page 23

The sensitivity of a model of the ambient stratosphere to empirical inputs

The usefulness of a general application of sensitivity studies to the modeling of the ambient stratosphere using the GSFC steady-state photochemistry model is examined. The comparison of theoretical model results with field measurements of trace species is used to assess the applicability of the model to the stratosphere, to help determine the eddy diffusion profile, and to validate mechanisms within the photochemistry of the model. Specific examples of each are discussed and the results of a comprehensive sensitivity study used to elucidate certain points about the comparison are analyzed.

Butler, D. M.↗

Structure of the solar chromosphere. III - Models of the EUV brightness components of the quiet-sun

The described investigation is concerned with the solution of the non-LTE optically thick transfer equations for hydrogen, carbon, and other constituents to determine semiempirical models for six components of the quiet solar chromosphere. For a given temperature-height distribution, the solution is obtained of the equations of statistical equilibrium, radiative transfer for lines and continua, and hydrostatic equilibrium to find the ionization and excitation conditions for each atomic constituent. The emergent spectrum is calculated, and a trial and error approach is used to adjust the temperature distribution so that the emergent spectrum is in best agreement with the observed one. The relationship between semiempirical models determined in this way and theoretical models based on radiative equilibrium is discussed by Avrett (1977). Harvard Skylab EUV observations are used to determine models for a number of quiet-sun regions.

Vernazza, J. E.↗

Dendritic solidification of alloys in low gravity

The effect of low gravity during solidification on primary and secondary dendrite spacings is investigated. Dendritic spacing data for iron-carbon-type alloys and superalloys which were directionally solidified during aircraft low-g trajectories are examined. The low-g data data are compared with data for metal model materials and binary alloys, and theoretical models to explain the dendritic spacing are considered.

Curreri, P. A.↗

Polarization and wavelength diversities of Gulf Stream fronts imaged by AIRSAR

During the 1990 Gulf Stream Experiment, NASA/JPL AIRSAR imaged the north edge of the Gulf Stream near the coast of Virginia. Simultaneous in-situ measurements of currents, temperatures, salinities, etc. were made for several crossings of the north edge by the R/V Cape Henlopen. Measurements identified two fronts with shearing and converging flows. The polarimetric SAR images from the fronts showed two bright linear features. One of them corresponds to the temperature front, which separated the warm Gulf Stream water to the south from a cool, freshwater filament to the north. The other line, located about 8 km north of the temperature front, is believed to correspond to the velocity front between the filament and the slope water. At these fronts, wave-current interactions produced narrow bands of steep and breaking waves manifesting higher radar returns in polarimetric SAR images. In general, our AIRSAR imagery shows that the signal-to-clutter ratio of radar cross sections for the temperature front is higher than that of the velocity front. In this paper, we study the polarization and wavelength diversities of radar response of these two fronts using the P-, L-, and C-Band Polarimetric SAR data. The north-south flight path of the AIRSAR crossed the temperature front several times and provided valuable data for analysis. Three individual passes are investigated. We found that for the temperature front, the cross-pol (HV) responses are much higher than co-pol responses (VV and HH), and that P-Band HV has the highest signal to clutter ratio. For the velocity front, the ratio is the strongest in P-Band VV, and it is indistinguishable for all polarizations in C-Band. The radar cross sections for all three polarization (HH, HV, and VV) and for all three bands are modelled using an ocean wave model and a composite Bragg scattering model. In our initial investigations, the theoretical model agrees qualitatively with the AIRSAR observations.

Lee, J. S.↗

Studies of Particle Acceleration, Transport and Radiation in Impulsive Phase of Solar Flares

Solar activity and its most prominent aspect, the solar flares, have considerable influence on terrestrial and space weather. Solar flares also provide a suitable laboratory for the investigation of many plasma and high energy processes important in the magnetosphere of the Earth and many other space and astrophysical situations. Hence, progress in understanding of flares will have considerable scientific and societal impact. The primary goal of this grant is the understanding of two of the most important problems of solar flare physics, namely the determination of the energy release mechanism and how this energy accelerates particles. This is done through comparison of the observations with theoretical models, starting from observations and gradually proceeding to theoretically more complex situations as the lower foundations of our understanding are secured. It is generally agreed that the source of the flare energy is the annihilation of magnetic fields by the reconnection process. Exactly how this energy is released or how it is dissipated remains controversial. Moreover, the exact mechanism of the acceleration of the particles is still a matter of debate. Data from many spacecrafts and ground based instruments obtained over the past decades have given us some clues. Theoretical analyses of these data have led to the standard thick target model (STT) where most of the released energy goes into an (assumed) power law spectrum of accelerated particles, and where all the observed radiations are the consequence of the interaction of these particles with the flare plasma. However, some theoretical arguments, and more importantly some new observations, have led us to believe that the above picture is not complete. It appears that plasma turbulence plays a more prominent role than suspected previously, and that it is the most likely agent for accelerating particles. The model we have developed is based on production of a high level of plasma waves and turbulence in the reconnection region above a flare loop. This turbulence accelerates particles stochastically which radiate some of their energy in this region but carry most of their energy to the footpoints of the loop, where they lose all their energy and radiate bulk of the observed radiation as in the traditional thick target model. In the past we have worked on various aspects of this model. We have evaluated the interaction rates of the plasma waves with electrons and ions, developed theoretical frameworks for the acceleration, transport and radiative processes, and produced numerical codes for the investigation of these processes. The goal of this grant has been further development and testing of this new paradigm, with emphases on the relative acceleration of electrons and ions and on a comprehensive investigation of the turbulence generation, cascade, and damping processes. We review several pieces of important evidence that we have uncovered indicating the crucial roles of turbulence, in and we describe accomplishments during the past two years of this grant.

Petrosian, Vahe↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Drozda, Tomasz G.↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Tomasz G Drozda↗

A model for interphase chromosomes and evaluation of radiation-induced aberrations

We have developed a theoretical model for evaluating radiation-induced chromosomal exchanges by explicitly taking into account interphase (G(0)/G(1)) chromosome structure, nuclear organization of chromosomes, the production of double-strand breaks (DSBs), and the subsequent rejoinings in a faithful or unfaithful manner. Each of the 46 chromosomes for human lymphocytes (40 chromosomes for mouse lymphocytes) is modeled as a random polymer inside a spherical volume. The chromosome spheres are packed randomly inside a spherical nucleus with an allowed overlap controlled by a parameter Omega. The rejoining of DSBs is determined by a Monte Carlo procedure using a Gaussian proximity function with an interaction range parameter sigma. Values of Omega and sigma have been found which yield calculated results of interchromosomal aberration frequencies that agree with a wide range of experimental data. Our preferred solution is one with an interaction range of 0.5 microm coupled with a relatively small overlap parameter of 0.675 microm, which more or less confirms previous estimates. We have used our model with these parameter values and with resolution or detectability limits to calculate yields of translocations and dicentrics for human lymphocytes exposed to low-LET radiation that agree with experiments in the dose range 0.09 to 4 Gy. Five different experimental data sets have been compared with the theoretical results. Essentially all of the experimental data fall between theoretical curves corresponding to resolution limits of 1 Mbp and 20 Mbp, which may reflect the fact that different investigators use different limits for sensitivity or detectability. Translocation yields for mouse lymphocytes have also been calculated and are in good agreement with experimental data from 1 cGy to 10 cGy. There is also good agreement with recent data on complex aberrations. Our model is expected to be applicable to both low- and high-LET radiation, and we include a sample prediction of the yield of interchromosomal rejoining in the dose range 0.22 Gy to 2 Gy of 1000 MeV/nucleon iron particles. This dose range corresponds to average particle traversals per nucleus ranging from 1.0 to 9.12.

Non-NASA Center↗

Beam Heating of Samples: Modeling and Verification

Energy absorbed from the X-ray beam by the sample requires cooling by forced convection (i.e. cryostream) to minimize temperature increase and the damage caused to the sample by the X-ray heating. In this presentation we will first review the current theoretical models and recent studies in the literature, which predict the sample temperature rise for a given set of beam parameters. It should be noted that a common weakness of these previous studies is that none of them provide actual experimental confirmation. This situation is now remedied in our investigation where the problem of x-ray sample heating is taken up once more. We have theoretically investigated, and at the same time, in addition to the numerical computations, performed experiments to validate the predictions. We have modeled, analyzed and experimentally tested the temperature rise of a 1 mm diameter glass sphere (sample surrogate) exposed to an intense synchrotron X-ray beam, while it is being cooled in a uniform flow of nitrogen gas. The heat transfer, including external convection and internal heat conduction was theoretically modeled using CFD to predict the temperature variation in the sphere during cooling and while it was subjected to an undulator (ID sector 19) X-ray beam at the APS. The surface temperature of the sphere during the X-ray beam heating was measured using the infrared camera measurement technique described in a previous talk. The temperatures from the numerical predictions and experimental measurements are compared and discussed. Additional results are reported for the two different sphere sizes and for two different supporting pin orientations.

Kazmierczak, Michael↗

Performance model of a recirculating stack nickel hydrogen cell

A theoretical model of the nickel hydrogen battery cell has been utilized to describe the chemical and physical changes during charge and overcharge in a recirculating stack nickel hydrogen cell. In particular, the movement of gas and electrolyte have been examined as a function of the amount of electrolyte put into the cell stack during cell activation, and as a function of flooding in regions of the gas screen in this cell design. Additionally, a two-dimensional variation on this model has been utilized to describe the effects of non-uniform loading in the nickel-electrode on the movement of gas and electrolyte within the recirculating stack nickel hydrogen cell. The type of nonuniform loading that has been examined here is that associated with higher than average loading near the surface of the sintered nickel electrode, a condition present to some degree in many nickel electrodes made by electrochemical impregnation methods. The effects of high surface loading were examined primarily under conditions of overcharge, since the movement of gas and electrolyte in the overcharging condition was typically where the greatest effects of non-uniform loading were found. The results indicate that significant changes in the capillary forces between cell components occur as the percentage of free volume in the stack filled by electrolyte becomes very high. These changes create large gradients in gas-filled space and oxygen concentrations near the boundary between the separator and the hydrogen electrode when the electrolyte fill is much greater than about 95 percent of the stack free volume. At lower electrolyte fill levels, these gaseous and electrolyte gradients become less extreme, and shift through the separator towards the nickel electrode. Similarly, flooding of areas in the gas screen cause higher concentrations of oxygen gas to approach the platinum/hydrogen electrode that is opposite the back side of the nickel electrode. These results illustrate the need for appropriate pore size distributions, and the maintenance of both convective electrolyte and gas flow paths through the stack, if the recirculating stack nickel hydrogen cell design is to work properly.

Zimmerman, Albert H.↗

Two-dimensional wind-tunnel interference from measurements on two contours

This paper describes how wall-induced velocities near a model in a two-dimensional wind tunnel can be estimated from upwash distributions measured along two contours surrounding a model. The method is applicable to flows that can be represented by linear theory. It was derived by applying the Schwarz Integral Formula separately to the two contours and by exploiting the free-air relationship between upwashes along the contours. Advantages of the method are that only one flow quantity need by measured and no representation of the model is required. A weakness of the method is that it assumes streamwise interference velocity vanishes far upstream of the model. This method was applied to a simple theoretical model of flow in a solid-wall wind tunnel. The theoretical interference velocities and the velocities computed using the method were in excellent agreement. The method was then used to analyze experimental data acquired during adaptive-wall experiments at Ames Research Center. This analysis confirmed that the wall adjustments reduced wall-induced velocities near the model.

Schairer, E. T.↗

Photodissociation Regions

The interstellar medium of galaxies is the reservoir out of which stars are born and into which stars inject newly created elements as they age. The physical properties of the interstellar medium are governed in part by the radiation emitted by these stars. Far-ultraviolet (6 eV< hNu < 13.6 eV) photons from massive stars dominate the heating and influence the chemistry of the neutral atomic gas and much of the molecular gas in galaxies. Predominantly neutral regions of the interstellar medium in which the heating and chemistry are regulated by far ultraviolet photons are termed Photodissociation Regions (PDRs). These regions are the origin of most of the non-stellar infrared (IR) and the millimeter and submillimeter CO emission from galaxies. The importance of PDRs has become increasingly apparent with the advances in IR and submillimeter astronomy. The IR emission from PDRs includes fine structure lines of C, C(+) and O; rovibrational lines of H2; rotational lines of CO; broad mid-IR features of polycyclic aromatic hydrocarbons; and a luminous underlying IR continuum from interstellar dust. The transition of H to H2 and C(+) to CO occurs within PDRs. Comparison of observations with theoretical models of PDRs enables one to determine the density and temperature structure, the elemental abundances, the level of ionization, and the radiation field. PDR models have been applied to interstellar clouds near massive stars, planetary nebulae, red giant outflows, photoevaporating planetary disks around newly formed stars, diffuse clouds, the neutral intercloud medium, and molecular clouds in the interstellar radiation field-in summary, much of the interstellar medium in galaxies. Theoretical PDR models explain the observed correlations of the [CII] 158, micrometers with the CO J=1-0 emission, the CO J=1-0 luminosity with the interstellar molecular mass, and the [CII] 158 micrometers plus [OI] 63 micrometers luminosity with the IR continuum luminosity. On a more global scale, PDR models predict the existence of two stable neutral phases of the interstellar medium, elucidate the formation and destruction of star-forming molecular clouds, and suggest radiation-induced feedback mechanisms that may regulate star formation rates and the column density of gas through giant molecular clouds.

Hollenbach, David J.↗

Photodissociation Regions in the Interstellar Medium of Galaxies

The interstellar medium of galaxies is the reservoir out of which stars are born and into which stars inject newly created elements as they age. The physical properties of the interstellar medium are governed in part by the radiation emitted by these stars. Far-ultraviolet (6 eV less than h(nu) less than 13.6 eV) photons from massive stars dominate the heating and influence the chemistry of the neutral atomic gas and much of the molecular gas in galaxies. Predominantly neutral regions of the interstellar medium in which the heating and chemistry are regulated by far ultraviolet photons are termed Photo-Dissociation Regions (PDRs). These regions are the origin of most of the non-stellar infrared (IR) and the millimeter and submillimeter CO emission from galaxies. The importance of PDRs has become increasingly apparent with advances in IR and submillimeter astronomy. The IR emission from PDRs includes fine structure lines of C, C+, and O; rovibrational lines of H2, rotational lines of CO; broad middle features of polycyclic aromatic hydrocarbons; and a luminous underlying IR continuum from interstellar dust. The transition of H to H2 and C+ to CO occurs within PDRs. Comparison of observations with theoretical models of PDRs enables one to determine the density and temperature structure, the elemental abundances, the level of ionization, and the radiation field. PDR models have been applied to interstellar clouds near massive stars, planetary nebulae, red giant outflows, photoevaporating planetary disks around newly formed stars, diffuse clouds, the neutral intercloud medium, and molecular clouds in the interstellar radiation field-in summary, much of the interstellar medium in galaxies. Theoretical PDR models explain the observed correlations of the [CII] 158 microns with the COJ = 1-0 emission, the COJ = 1-0 luminosity with the interstellar molecular mass, and the [CII] 158 microns plus [OI] 63 microns luminosity with the IR continuum luminosity. On a more global scale, MR models predict the existence of two stable neutral phases of the interstellar medium, elucidate the formation and destruction of star-forming molecular clouds, and suggest radiation-induced feedback mechanisms that may regulate star formation rates and the column density of gas through giant molecular clouds.

Hollenbach, David J.↗

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗

Development of global model for atmospheric backscatter at CO2 wavelengths

The improvement of an understanding of the variation of the aerosol backscattering at 10.6 micron within the free troposphere and the development model to describe this was undertaken. The analysis combines theoretical modeling with the results contained within three independent data sets. The data sets are obtained by the SAGE I/SAM II satellite experiments, the GAMETAG flight series and by direct backscatter measurements. The theoretical work includes use of a bimodal, two component aerosol model, and the study of the microphysical and associated optical changes occurring within an aerosol plume. A consistent picture is obtained, which describes the variation of the aerosol backscattering function in the free troposphere with altitude, latitude, and season. Most data are available and greatest consistency is found inside the Northern Hemisphere.

Kent, G. S.↗

Closed coronal structures. III - Comparison of static models with X-ray, EUV, and radio observations

Numerical models of static coronal loops in energy balance are compared with high spatial resolution observations of extreme ultraviolet lines, broad-band X-ray emission, and interferometric observations at 2.8 cm of a solar active region. Difficulties of using scaling laws to test static models of coronal loops are reviewed. The theoretical model used for the comparison is summarized; the detailed X-ray, EUV, and microwave observations of the selected active region are presented; and the comparison of the model with the observations is performed. It is shown that simple static models with conductive flux vanishing at the loop base reproduce satisfactorily the observed properties in the upper portion of loop structures from compact, high-pressure loops in the core of the region to more extended, fainter loops and to large-scale loops interconnecting different active regions. Effects of changing loop parameters are investigated, and it is argued, that in contrast to the present approach, scaling laws cannot be used to discriminate between different static energy balance models. Some discrepancy is found between model predictions and observations for the lower sections of loop structures. Possible causes of the discrepancy are discussed.

Pallavicini, R.↗

Wet model of Saturn's ionosphere: Water from the rings

Current theoretical models of Saturn's ionosphere are difficult to reconcile with the ionospheric electron density profiles obtained from the Pioneer and Voyager radio occultation observations and the large diurnal variation of maximum ionospheric electron density deduced from studies of Saturn lightning discharges. A model of Saturn's ionosphere is proposed in which water plays a major role as a minor constituent present by virtue of downward diffusion from an external source. This model of the Saturn ionosphere is a classical 'F2' type layer resulting from the photodissociative production of H(+) from H2 and rapid chemical loss due to a series of charge exchange reactions with water. A planet-wide influx of about 4x10 to the 7th power molecules/sec/sq cm of water from the rings is consistent with the observed ionospheric electron densities and estimates of influx due to micrometeoride bombardment of the rings. An enhanced influx of water occurs at latitudes (-38 deg, +44 deg) magnetically connected to the inner edge of Saturn's B ring which results from an electromagnetic erosion process contributing substantially to the (local) upper atmosphere water content. Present day influx at these latitudes is possibly as large as 2x10 to the 9th power molecules/sec/sq cm.

Connerney, J. E. P.↗