Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distribution function”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Parameterizing Grid-Averaged Longwave Fluxes for Inhomogeneous Marine Boundary Layer Clouds

This paper examines the relative impacts on grid-averaged longwave flux transmittance (emittance) for Marine Boundary Layer (MBL) cloud fields arising from horizontal variability of optical depth tau and cloud sides, First, using fields of Landsat-inferred tau and a Monte Carlo photon transport algorithm, it is demonstrated that mean all-sky transmittances for 3D variable MBL clouds can be computed accurately by the conventional method of linearly weighting clear and cloudy transmittances by their respective sky fractions. Then, the approximations of decoupling cloud and radiative properties and assuming independent columns are shown to be adequate for computation of mean flux transmittance. Since real clouds have nonzero geometric thicknesses, cloud fractions A'(sub c) presented to isotropic beams usually exceed the more familiar vertically projected cloud fractions A(sub c). It is shown, however, that when A(sub c)less than or equal to 0.9, biases for all-sky transmittance stemming from use of A(sub c) as opposed to A'(sub c) are roughly 2-5 times smaller than, and opposite in sign to, biases due to neglect of horizontal variability of tau. By neglecting variable tau, all-sky transmittances are underestimated often by more than 0.1 for A(sub c) near 0.75 and this translates into relative errors that can exceed 40% (corresponding errors for all-sky emittance are about 20% for most values of A(sub c). Thus, priority should be given to development of General Circulation Model (GCM) parameterizations that account for the effects of horizontal variations in unresolved tau, effects of cloud sides are of secondary importance. On this note, an efficient stochastic model for computing grid-averaged cloudy-sky flux transmittances is furnished that assumes that distributions of tau, for regions comparable in size to GCM grid cells, can be described adequately by gamma distribution functions. While the plane-parallel, homogeneous model underestimates cloud transmittance by about an order of magnitude when 3D variable cloud transmittances are less than or equal to 0.2 and by approx. 20% to 100% otherwise, the stochastic model reduces these biases often by more than 80%.

Barker, Howard W.↗
NASA NTRS1998

Study of Atmospheric Nitric Oxide

We investigated the contribution of energetic nitrogen atoms to the production of nitric oxide in the thermosphere and their influence on the infrared emission spectrum. The nitric oxide molecules are important contributors to the cooling of the atmosphere. We first pointed out that in determining the energy distribution of the nitrogen atoms, it is important to take into account the thermal motion of the atmospheric gases. It had been ignored in all earlier studies. The source spectra are broadened considerably by the center of mass motion of the reactants. We worked out the consequences for the production of nitric oxide at night, using as sources of energetic N atoms, NO(+) + e yield N + O, N(D-2) + O yield N + O. The high energy tail is enhanced by orders of magnitude. We had earlier suggested (Sharma et al. 1993) that the reaction of energetic nitrogen atoms with O2 was responsible for the rotationally enhanced NO identified in the infrared spectrum. Our calculations provided quantitative confirmation of the suggestion. We proceeded to explore the validity of another approximation used in earlier analyses, the hard sphere approximation for the energy loss in elastic collisions. We carried out precise quantum mechanical calculations of the elastic 2 differential scattering of nitrogen atoms in collisions with oxygen atoms and showed that although the hard sphere approximation was nowhere of high precision, reasonable results could be obtained with an effective cross section of 6 x 10(exp 15)sq cm. We also initiated a program to include inelastic energy loss processes in the determination of the energy distribution function. We began a calculation of the rotation and vibrational excitation cross sections of molecular nitrogen and nitrogen atoms and developed a method for including inelastic energy loss as a function of scattering angle in the Boltzmann equation. A procedure for obtaining the solution of the Boltzman equation was worked out.

Dalgarno, A.↗
NASA NTRS1998

Development of a Probabilistic Dynamic Synthesis Method for the Analysis of Nondeterministic Structures

Accounting for the statistical geometric and material variability of structures in analysis has been a topic of considerable research for the last 30 years. The determination of quantifiable measures of statistical probability of a desired response variable, such as natural frequency, maximum displacement, or stress, to replace experience-based "safety factors" has been a primary goal of these studies. There are, however, several problems associated with their satisfactory application to realistic structures, such as bladed disks in turbomachinery. These include the accurate definition of the input random variables (rv's), the large size of the finite element models frequently used to simulate these structures, which makes even a single deterministic analysis expensive, and accurate generation of the cumulative distribution function (CDF) necessary to obtain the probability of the desired response variables. The research presented here applies a methodology called probabilistic dynamic synthesis (PDS) to solve these problems. The PDS method uses dynamic characteristics of substructures measured from modal test as the input rv's, rather than "primitive" rv's such as material or geometric uncertainties. These dynamic characteristics, which are the free-free eigenvalues, eigenvectors, and residual flexibility (RF), are readily measured and for many substructures, a reasonable sample set of these measurements can be obtained. The statistics for these rv's accurately account for the entire random character of the substructure. Using the RF method of component mode synthesis, these dynamic characteristics are used to generate reduced-size sample models of the substructures, which are then coupled to form system models. These sample models are used to obtain the CDF of the response variable by either applying Monte Carlo simulation or by generating data points for use in the response surface reliability method, which can perform the probabilistic analysis with an order of magnitude less computational effort. Both free- and forced-response analyses have been performed, and the results indicate that, while there is considerable room for improvement, the method produces usable and more representative solutions for the design of realistic structures with a substantial savings in computer time.

Brown, A. M.↗
NASA NTRS1998

The Use of Langmuir Probes in Non-Maxwellian Space Plasmas

Disturbance of the Maxwellian plasma may occur in the vicinity of a spacecraft due to photoemission, interactions between the spacecraft and thermospheric gases, or electron emissions from other devices on the spacecraft. Significant non-maxwellian plasma distributions may also occur in nature as a mixture of ionospheric and magnetospheric plasmas or secondaries produced by photoionization in the thermosphere or auroral precipitation. The general formulas for current collection (volt-ampere curves) by planar, cylindrical, and spherical Langmuir probes in isotropic and anisotropic non-maxwellian plasmas are examined. Examples are given of how one may identify and remove the non-maxwellian components in the Langmuir probe current to permit the ionospheric parameters to be determined. Theoretical volt-ampere curves presented for typical examples of non-maxwellian distributions include: two-temperature plasmas and a thermal plasma with an energetic electron beam. If the non-ionospheric electrons are Maxwellian at a temperature distinct from that of the ionosphere electrons, the volt-ampere curves can be fitted directly to obtain the temperatures and densities of both electron components without resorting to differenting the current. For an arbitrary isotropic distribution, the current for retarded particles is shown to be identical for the three geometries. For anisotropic distributions, the three probe geometries are not equally suited for measuring the ionospheric electron temperature and density or for determining the distribution function in the presence of non-maxwellian back-round electrons.

Hoegy, Walter R.↗
NASA NTRS1998

Probabilistic Modeling of Ceramic Matrix Composite Strength

Uncertainties associated with the primitive random variables such as manufacturing process (processing temperature, fiber volume ratio, void volume ratio), constituent properties (fiber, matrix and interface), and geometric parameters (ply thickness, interphase thickness) have been simulated to quantify the scatter in the first matrix cracking strength (FMCS) and the ultimate tensile strength of SCS-6/RBSN (SiC fiber (SCS-6) reinforced reaction-bonded silicon nitride composite) ceramic matrix composite laminate at room temperature. Cumulative probability distribution function for the FMCS and ultimate tensile strength at room temperature (RT) of (0)(sub 8), (0(sub 2)/90(sub 2), and (+/-45(sub 2))(sub S) laminates have been simulated and the sensitivity of primitive variables to the respective strengths have been quantified. Computationally predicted scatter of the strengths for a uniaxial laminate have been compared with those from limited experimental data. Also the experimental procedure used in the tests has been described briefly. Results show a very good agreement between the computational simulation and the experimental data. Dominating failure modes in (0)(sub 8), (0/90)(sub s) and (+/-45)(sub S) laminates have been identified. Results indicate that the first matrix cracking strength for the (0)(sub S), and (0/90)(sub S) laminates is sensitive to the thermal properties, modulus and strengths of both the fiber and matrix whereas the ultimate tensile strength is sensitive to the fiber strength and the fiber volume ratio. In the case of a (+/-45)(sub S), laminate, both the FMCS and the ultimate tensile strengths have a small scatter range and are sensitive to the fiber tensile strength as well as the fiber volume ratio.

Shan, Ashwin R.↗
NASA NTRS1997

A Model for Lower Hybrid Wave Excitation Compared with Observations by Viking

The mechanism of lower hybrid wave (LHW) excitation due to the O+ relative drift in a plasma subjected to low-frequency waves (LFWs) is used for analysis of Viking satellite data for events in the cusp/cleft region. In some cases, such a mechanism leads to LHW energy densities and ion distribution functions close to those observed, suggesting the proposed mechanism is a plausible candidate to explain certain classes of LHW generation events in space plasmas.

Khazanov, G. V.↗
NASA NTRS1998

Direct Acceleration of Pickup Ions at The Solar Wind Termination Shock: The Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

Ellison, Donald C.↗
NASA NTRS1998

A Gas-Kinetic Scheme for Reactive Flows

In this paper, the gas-kinetic BGK scheme for the compressible flow equations is extended to chemical reactive flow. The mass fraction of the unburnt gas is implemented into the gas kinetic equation by assigning a new internal degree of freedom to the particle distribution function. The new variable can be also used to describe fluid trajectory for the nonreactive flows. Due to the gas-kinetic BGK model, the current scheme basically solves the Navier-Stokes chemical reactive flow equations. Numerical tests validate the accuracy and robustness of the current kinetic method.

Lian,Youg-Sheng↗
NASA NTRS1996

Structure of the Magnetotail Current Sheet

An orbit tracing technique was used to generate current sheets for three magnetotail models. Groups of ions were followed to calculate the resulting cross-tail current. Several groups then were combined to produce a current sheet. The goal is a model in which the ions and associated electrons carry the electric current distribution needed to generate the magnetic field B in which ion orbits were traced. The region -20 R(sub E) less than x less than - 14 R(sub E) in geocentric solar magnetospheric coordinates was studied. Emphasis was placed on identifying the categories of ion orbits which contribute most to the cross-tail current and on gaining physical insight into the manner by which the ions carry the observed current distribution. Ions that were trapped near z = 0, ions that magnetically mirrored throughout the current sheet, and ions that mirrored near the Earth all were needed. The current sheet structure was determined primarily by ion magnetization currents. Electrons of the observed energies carried relatively little cross-tail current in these quiet time current sheets. Distribution functions were generated and integrated to evaluate fluid parameters. An earlier model in which B depended only on z produced a consistent current sheet, but it did not provide a realistic representation of the Earth's middle magnetotail. In the present study, B changed substantially in the x and z directions but only weakly in the y direction within our region of interest. Plasmas with three characteristic particle energies were used with each of the magnetic field models. A plasma was found for each model in which the density, average energy, cross-tail current, and bulk flow velocity agreed well with satellite observations.

Larson, Douglas J.↗
NASA NTRS1997

Kinetic-Energy Distribution of D(2p) Atoms from Analysis of the D Lyman-Alpha Line Profile

The kinetic-energy distribution of D(2p) atoms resulting from electron-impact dissociation of D2 has been measured. A high-resolution vacuum ultraviolet spectrometer was employed for the first measurement of the D Lyman-alpha (D L(alpha)) emission line profiles at 20- and 100-eV excitation energies. Analysis of the deconvoluted line profile of D L(alpha) at 100 eV reveals the existence of a narrow line central peak of 29+/-2 mA full width at half maximum and a broad pedestal wing structure about 190 mA wide. The wings of the line can be used to determine the fast atom distribution. The wings of D L(alpha) arise from dissociative excitation of a series of doubly excited states that cross the Franck-Condon region between 23 and 40 eV. The fast atom distribution at 100-eV electron impact energy spans the energy range from 1 to 10 eV with a peak value near 6 eV. Slow D(2p) atoms characterized by a distribution function with peak energy near 100 meV produce the central peak profile, which is nearly independent of the impact energy. The deconvoluted line profiles of the central peak at 20 eV for dissociative excitation of D2 and H2 are fitted with an analytical function for use in calibration of space flight instrumentation equipped with a D/H absorption cell. The kinetic-energy and line profile results are compared to similar measurements for H2. The absolute cross sections for the line center (slow atoms) and wings (fast atoms) and total emission line profile were measured from threshold to 400 eV. Analytical model coefficients are given for the energy dependence of the measured slow atom cross section.

Ciocca, M.↗
NASA NTRS1998

Auroral Microphysics Rocket (AMICIST) and Reflight of the Phaze Sounding Rocket

This grant was originally awarded for the flight of the AMICIST sounding rocket. However, upon launch failure of the PHAZE rocket, additional resources were placed in this grant to cover the launch of the PHAZE II rocket. AMICIST was successfully launched from the Poker Flat Range on February 24, 1995, and the PHAZE II was successfully launched also from the Poker Flat Range on February 10, 1997. The major objective of the AMICIST flight was to investigate the bursts of transverse ion acceleration occurring during aurora, commonly known as Lower Hybrid Solitary Structures, with the flight of two scientific payloads to unravel space from time effects. The data clearly showed that the structures of ion acceleration were a spatial phenomena having a scale size transverse to the local magnetic field of about 100 m. On the other hand, structures in the auroral electrons were generally observed on this flight to be temporal features that occurred at the same time at the two payloads separated by a few kilometers. The primary objective of the PHAZE rocket flight was to further study the temporal features of auroral electron precipitation having many electron detectors, some fixed in energy, that could study the distribution function of the auroral electrons on time scales of a fraction of a millisecond. A important paper on this topic has just been submitted for publication which uses the PHAZE data to show that the potential structure (electric field parallel to the magnetic field) that accelerates the auroral electrons within one Earth radius of the ground, is not static, but rather fluctuates with frequencies close to the local proton and hydrogen gyrofrequencies. The fluctuation appears to be an actual turning on and off of the electric field at these frequencies. When the potential is turned off, ambient electrons can enter the region and be accelerated along B when the potential is on creating field- aligned bursts which manifest themselves as flickering aurora seen from the ground. The on/off model can explain many of the dynamic features of the aurora which have been unexplained since man first looked up at the visual display. Much needs to be done yet to understand how ions couple into this process.

Arnoldy Roger L.↗
NASA NTRS1998

On Geomagnetism and Paleomagnetism

A statistical description of Earth's broad scale, core-source magnetic field has been developed and tested. The description features an expected, or mean, spatial magnetic power spectrum that is neither "flat" nor "while" at any depth, but is akin to spectra advanced by Stevenson and McLeod. This multipole spectrum describes the magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Natural variations of core multipole powers about their mean values are to be expected over geologic time and are described via trial probability distribution functions that neither require nor prohibit magnetic isotropy. The description is thus applicable to core-source dipole and low degree non-dipole fields despite axial dipole anisotropy. The description is combined with main field models of modem satellite and surface geomagnetic measurements to make testable predictions of: (1) the radius of Earth's core, (2) mean paleomagnetic field intensity, and (3) the mean rates and durations of both dipole power excursions and durable axial dipole reversals. The predicted core radius is 0.7% above the 3480 km seismologic value. The predicted root mean square paleointensity (35.6 mu T) and mean Virtual Axial Dipole Moment (about 6.2 lx 1022 Am(exp 2)) are within the range of various mean paleointensity estimates. The predicted mean rate of dipole power excursions, as defined by an absolute dipole moment <20% of the 1980 value, is 9.04/Myr and 14% less than obtained by analysis of a 4 Myr paleointensity record. The predicted mean rate of durable axial dipole reversals (2.26/Myr) is 2.3% more than established by the polarity time-scale for the past 84 Myr. The predicted mean duration of axial dipole reversals (5533 yr) is indistinguishable from an observational value. The accuracy of these predictions demonstrates the power and utility of the description, which is thought to merit further development and testing. It is suggested that strong stable stratification of Earth's uppermost outer core leads to a geologically long interval of no dipole reversals and a very nearly axisymmetric field outside the core. Statistical descriptions of other planetary magnetic fields are outlined.

Voorhies, Coerte V.↗
NASA NTRS1999

Aircraft NO(x) had no Unique Fingerprint on Sonex; Lightning Dominated Fresh NO(x) Sources

Key questions to which SONEX was directed were the following: Can aircraft corridors be detected? Is there a unique tracer for aircraft NO(x)? Can a "background" NO(x) (or NO(y) be defined? What fraction of NO(x) measured during SONEX was from aircraft? How representative was SONEX of the North Atlantic in 1997 and how typical of other years? We attempt to answer these questions through species-species correlations, probability distribution functions (PDFs), and meteorological history. There is not a unique aircraft tracer, largely due to the high variability of air mass origins and tracer ratios, which render "average" quantities meaningless. The greatest NO and NO(y) signals were associated with lightning and convective NO sources. Well-defined background CO, NO(y) and NO(y)/ozone ratio appear in subsets of two cross-track flights with subtropical origins and five flights with predominantly mid-latitude air. Forty percent of the observations on these 7 flights showed NO(y)/ozone to be above background, evidently due to unreacted NO(x). This NO(x) is a combination of aircraft, lightning and surface pollution injected by convection. The strongly subtropical signatures in SONEX observations, confirmed by pv (potential vorticity) values along flight tracks, argues for most of the unreacted NO(x) originating from lightning. Potential vorticity statistics along SONEX flight tracks in 1992-1998, and for the North Atlantic as a whole, show the SONEX meteorological environment to be representative of the North Atlantic flight corridor in the October-November period.

Thompson, A.↗
NASA NTRS1998

Probabilistic Evaluation of Fuselage-Type Composite Structures

A methodology is developed to simulate computationally the uncertain behavior of composite structures. The uncertain behavior includes buckling loads, natural frequencies, displacements, stress/strain, etc., which are the consequences of the random variation (scatter) of the primitive (independent random) variables in the constituent, plv. laminate and structural levels. This methodology is implemented in a computer code integrated probabilistic assessment of composite structures (IPACS). A fuselage-type composite structure is analyzed to demonstrate the code's capability . The probability distribution functions of the buckling loads, natural frequency, displacement, strain and stress are computed. The sensitivity of each primitive (independent random) variable to a given structural response is also identified from the analyses.

Shiao, Michael C.↗
NASA NTRS1997

Collisionless Plasma Modeling in an Arbitrary Potential Energy Distribution

A new technique for calculating a collisionless plasma along a field line is presented. The primary feature of the new model is that it can handle an arbitrary (including nonmonotonic) potential energy distribution. This was one of the limiting constraints on the existing models in this class, and these constraints are generalized for an arbitrary potential energy composition. The formulation for relating current density to the field-aligned potential as well as formulas for density, temperature and energy flux calculations are presented for several distribution functions, ranging from a bi-Lorentzian with a loss cone to an isotropic Maxwellian. A comparison of these results with previous models shows that the formulation reduces.to the earlier models under similar assumptions.

Liemohn, M. W.↗
NASA NTRS1998

A Simple Model of Cirrus Horizontal Inhomogeneity and Cloud Fraction

A simple model of horizontal inhomogeneity and cloud fraction in cirrus clouds has been formulated on the basis that all internal horizontal inhomogeneity in the ice mixing ratio is due to variations in the cloud depth, which are assumed to be Gaussian. The use of such a model was justified by the observed relationship between the normalized variability of the ice water mixing ratio (and extinction) and the normalized variability of cloud depth. Using radar cloud depth data as input, the model reproduced well the in-cloud ice water mixing ratio histograms obtained from horizontal runs during the FIRE2 cirrus campaign. For totally overcast cases the histograms were almost Gaussian, but changed as cloud fraction decreased to exponential distributions which peaked at the lowest nonzero ice value for cloud fractions below 90%. Cloud fractions predicted by the model were always within 28% of the observed value. The predicted average ice water mixing ratios were within 34% of the observed values. This model could be used in a GCM to produce the ice mixing ratio probability distribution function and to estimate cloud fraction. It only requires basic meteorological parameters, the depth of the saturated layer and the standard deviation of cloud depth as input.

Smith, Samantha A.↗
NASA NTRS1997

Banded Electron Structure Formation in the Inner Magnetosphere

Banded electron structures in energy-time spectrograms have been observed in the inner magnetosphere concurrent with a sudden relaxation of geomagnetic activity. In this study, the formation of these banded structures is considered with a global, bounce-averaged model of electron transport, and it is concluded that this structure is a natural occurrence when plasma sheet electrons are captured on closed drift paths near the Earth. These bands do not appear unless there is capture of plasma sheet electrons; convection along open drift paths making open pass around the Earth do not have time to develop this feature. The separation of high-energy bands from the injection population due to the preferential advection of the gradient-curvature drift creates spikes in the energy distribution, which overlap to form a series of bands in the energy spectrograms. The lowest band is the bulk of the injected population in the sub-key energy range. Using the Kp history for an observed banded structure event, a cloud of plasma sheet electrons is captured and the development of their distribution function is examined and discussed.

Liemohn, M. W.↗
NASA NTRS1998

Contamination Observed on the Passive Optical Sample Assembly (POSA)-1 Experiment

The Passive Optical Sample Assembly-1 (POSA-1), part of the Mir Environmental Effects Payload (MEEP), was designed to study the combined effects of contamination, atomic oxygen, ultraviolet radiation, vacuum, thermal cycling, and other constituents of the space environment on spacecraft materials. The MEEP program is a Phase 1 International Space Station Risk Mitigation Experiment. The payload was attached by EVA to the exterior of the Mir docking module during the Space Shuttle mission STS-76. It was removed during STS-86 after exposure to the Mir environment for 18 months. POSA-1 consists of nearly 400 samples of thermal control paints, chemical conversion coatings, mirrors, optics, and other candidate materials for International Space Station Alpha (ISSA). POSA-1 also flew state-of-the-art materials, such as electrically conductive thermal control coatings developed under contract by IITRI and AZ Technology, and passive instruments for monitoring the atomic oxygen and ultraviolet radiation dose to the experiment. Contamination was detected on the POSA-1 experiment. On the side facing space (-Z(sub B) direction using the Mir coordinate system), visible contamination was observed. The contamination was uniform as would occur from a slow photodeposition process. A very definite film was deposited on the optical samples. The deposition appears to have directionality with definite shadowing effects. On the side facing the main Mr core (+ZB direction using the Mir coordinate system), no visible contamination was noted. Pre- and post-flight characterization of candidate spacecraft materials is discussed, particularly the effect of contamination on total hemispherical spectral reflectance (250 - 2500 nm), infrared emittance, the bi-directional reflectance distribution function (BRDF), and vacuum ultraviolet spectral reflectance. Extensive photography under normal and black light was performed. Results of electron spectroscopy for chemical analysis (ESCA) and Fourier Transform Infrared Analysis (FTIR) on selected samples are also discussed.

Zwiener, James M.↗