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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

A Canonical Response in Rainfall Characteristics to Global Warming: Projections by IPCC CMIP5 Models

Changes in rainfall characteristics induced by global warming are examined based on probability distribution function (PDF) analysis, from outputs of 14 IPCC (Intergovernmental Panel on Climate Change), CMIP (5th Coupled Model Intercomparison Project) models under various scenarios of increased CO2 emissions. Results show that collectively CMIP5 models project a robust and consistent global and regional rainfall response to CO2 warming. Globally, the models show a 1-3% increase in rainfall per degree rise in temperature, with a canonical response featuring large increase (100-250 %) in frequency of occurrence of very heavy rain, a reduction (5-10%) of moderate rain, and an increase (10-15%) of light rain events. Regionally, even though details vary among models, a majority of the models (>10 out of 14) project a consistent large scale response with more heavy rain events in climatologically wet regions, most pronounced in the Pacific ITCZ and the Asian monsoon. Moderate rain events are found to decrease over extensive regions of the subtropical and extratropical oceans, but increases over the extratropical land regions, and the Southern Oceans. The spatial distribution of light rain resembles that of moderate rain, but mostly with opposite polarity. The majority of the models also show increase in the number of dry events (absence or only trace amount of rain) over subtropical and tropical land regions in both hemispheres. These results suggest that rainfall characteristics are changing and that increased extreme rainfall events and droughts occurrences are connected, as a consequent of a global adjustment of the large scale circulation to global warming.

Lau, William K. M.↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗
NASA NTRS2007

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗
NASA NTRS2011

Characterizing Vegetation Cover in Global Savannas with an Annual Foliage Clumping Index Derived from the MODIS BRDF Product

The global savanna biome is characterized by enormous diversity in the physiognomy and spatial structure of the vegetation. The foliage clumping index can be calculated from bidirectional reflectance distribution function (BRDF) data. It measures the response of the darkspot reflectance to increased shadow associated with clumped vegetation and is related to leaf area index. Clumping index theoretically declines with increasing woody cover until the tree canopy begins to become uniform. In this study, clumping index is calculated for Moderate Resolution Imaging Spectroradiometer BRDF data for the Australian tropical savanna, the tropical savannas of South America, and the tropical savannas of east, west and southern Africa and compared with site-based measurements of tree canopy cover, and with area-based classifications of land cover. There were differences in sensitivity of clumping index between red and near-infrared reflectance channels, and between savanna systems with markedly different woody vegetation physiognomy. Clumping index was broadly related to foliage cover from historical site data in Australia and in West Africa and Kenya, but not in Southern Africa nor with detailed site-based demographic data in the cerrado of Brazil. However, clumping index decreased with proportion of woody cover in land cover datasets for east Africa, Australia and the Colombian Llanos. There was overlap in the range of clumping index values for forest, cerrado and campo land covers in Brazil. Clumping index was generally negatively correlated with percentage tree cover from the MODIS Vegetation Continuous Fields product, but regional differences in the relationship were evident. There were large differences in the frequency distributions of clumping index from savanna, woody savanna and grassland land cover classes between global ecoregions. The clumping index shows differing sensitivity to savanna woody cover for red and NIR reflectance, and requires regional calibration for application as a universal indicator.

Hill, Michael J.↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗
NASA NTRS1994

Effect of B(sub y) on neutral line ridges and dynamical source ordering

In this paper we study the effect of a uniform magnetic field B(sub y) in the cross-tail y direction on the dynamics and distribution of energetic ions in a current sheet model including a neutral line. Martin and Speiser (1988) have previously shown that a 'ridge' in the velocity space distribution function is a remote neutral line signature in a two-dimensional field without B(sub y). Our results show that as B(sub y) is increased, using nominal tail parameters, there is very little change in the ridge signature for small B(sub y) (up to about 4 times B(sub z)). For intermediate values (up to about 10 times times B(sub Y)) the ridge becomes observable further from the neutral line, while close to the X line the ridge is strongly modified. For large B(sub y) (of the order of B(sub x)) the ridge evolves into the new neutral line signatures which depend strongly on whether the observations point is above or below the current sheet. These results are used to estimate B(sub y) less than 2 nT in the Active Magentospheric Particle Traces Explorers/Ion Release Module (AMPTE/IRM) event modeled by Speiser and Martin (1994), which is consistent with onboard magnetometer measurements. We further show that the cause of the new structures appearing at large B(sub Y) is the same source seperation effect that produced the original ridge and that all these signatures are relatively insensitive to a small parallel electric field. Finally, we point our a dynamical ordering which occurs at large B(sub y): Using high-resolution simulations, we find a common boundary separating regions of positive and negative initial x, y, and z position, as well as initial pitch angle. This boundary is purely dynamical in nature and is independent of modeled source distributions.

Martin, R. F., Jr.↗
NASA NTRSDate not supplied

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗
NASA NTRSDate not supplied

Overview of NASA's Ocean Color Instrument Solar Calibration Architecture, Pre-Launch Tests and Preliminary On-Orbit Results

Launched in February 2024, the PACE mission represents NASA’s next investment in ocean biology, clouds, and aerosol data records. A key feature of PACE is the inclusion of an advanced satellite radiometer known as the Ocean Color Instrument (OCI), a global mapping radiometer that combines multispectral and hyperspectral remote sensing. Like its predecessors, OCI will provide two day global coverage of TOA radiances. Unlike its predecessors, OCI will cover a spectral range from 340nm to 2260nm. Below 900nm, OCI will include two spectrometers that continuously span the ultraviolet to 600nm and 600nm to near-infrared spectral regions to provide hyperspectral radiances sampled every 2.5 nm, with a bandwidth of 5 nm for each channel. Wavelengths above 900nm are measured in seven discrete multispectral bands of varying bandwidths, six of which are at similar wavelengths to those on heritage missions to support both atmospheric and ocean color applications. Nominal spatial resolution is similar to the SeaWiFS instrument with 1050 m at nadir. As for SeaWiFS, the pixel size increases due to a ~20 degree tilt and as a function of scan angle. Variations in the radiometric sensitivity of each OCI channel over time will be monitored by solar diffuser measurements for short term instrument gain adjustments and independent lunar measurements for trend adjustments of long time periods, similar to the approach employed for the VIIRS instrument [4]. The OCI flight-unit was built at NASA’s Goddard Space Flight Center. At the time of this writing, OCI has completed on-orbit commissioning activities and normal science operations have begun. A key aspect of the OCI architecture is the capability to trend absolute and relative calibration changes over the course of mission life with solar calibration. Every 24 hours, the PACE spacecraft performs an inertial hold as the ground track nears the North Pole which orients a Quasi-Volume Diffuser (QVD) mounted on OCI towards the sun. By knowing the irradiance of the sun and the reflectivity of the target, the absolute radiance at the input to the OCI aperture can be computed as OCI scans the target. The allowable absolute uncertainty budget for each solar calibration measurement is 1.6% 1-sigma below 900nm at beginning of life (BOL) and the allowable relative uncertainty budget is ~0.26% 1-sigma. The Solar Calibration Assembly (SCA) consists of three targets selectable via a single mechanism which also opens a door. The targets consist of a Daily Bright Target (DBT), Monthly Bright Target (MBT), and Daily Dim Target (DDT). The bright targets are quartz QVDs with the monthly target being used to track the degradation of the daily target. The dim target is used to track CCD linearity using Progressive Time-Delay Integration (PTDI). A composite baffle is attached to the SCA housing aperture to block Earth shine and stray light from the spacecraft. The SCA assembly is mounted to a view port ~90° from OCI nadir. This paper provides an overview of driving solar calibration requirements, error-budgets and early trade studies which drove the solar calibration assembly (SCA) architecture and on-orbit maneuver. Measurements of the diffuser Bidirectional Reflectance Distribution Function (BRDF) at TNO, Netherlands and GSFC are briefly described. Optical modelling and test results at the sub-system and instrument level are included. Finally, preliminary measurements on-orbit are compared to pre-launch predictions.

Joseph J Knuble↗
NASA NTRSDate not supplied

Overview of NASA's Ocean Color Instrument Solar Calibration Architecture, Pre-Launch Tests and Preliminary On-Orbit Results

Launched in February 2024, the PACE mission represents NASA’s next investment in ocean biology, clouds, and aerosol data records. A key feature of PACE is the inclusion of an advanced satellite radiometer known as the Ocean Color Instrument (OCI), a global mapping radiometer that combines multispectral and hyperspectral remote sensing. Like its predecessors, OCI will provide two day global coverage of TOA radiances. Unlike its predecessors, OCI will cover a spectral range from 340nm to 2260nm. Below 900nm, OCI will include two spectrometers that continuously span the ultraviolet to 600nm and 600nm to near-infrared spectral regions to provide hyperspectral radiances sampled every 2.5 nm, with a bandwidth of 5 nm for each channel. Wavelengths above 900nm are measured in seven discrete multispectral bands of varying bandwidths, six of which are at similar wavelengths to those on heritage missions to support both atmospheric and ocean color applications. Nominal spatial resolution is similar to the SeaWiFS instrument with 1050 m at nadir. As for SeaWiFS, the pixel size increases due to a ~20 degree tilt and as a function of scan angle. Variations in the radiometric sensitivity of each OCI channel over time will be monitored by solar diffuser measurements for short term instrument gain adjustments and independent lunar measurements for trend adjustments of long time periods, similar to the approach employed for the VIIRS instrument [4]. The OCI flight-unit was built at NASA’s Goddard Space Flight Center. At the time of this writing, OCI has completed on-orbit commissioning activities and normal science operations have begun. A key aspect of the OCI architecture is the capability to trend absolute and relative calibration changes over the course of mission life with solar calibration. Every 24 hours, the PACE spacecraft performs an inertial hold as the ground track nears the North Pole which orients a Quasi-Volume Diffuser (QVD) mounted on OCI towards the sun. By knowing the irradiance of the sun and the reflectivity of the target, the absolute radiance at the input to the OCI aperture can be computed as OCI scans the target. The allowable absolute uncertainty budget for each solar calibration measurement is 1.6% 1-sigma below 900nm at beginning of life (BOL) and the allowable relative uncertainty budget is ~0.26% 1-sigma. The Solar Calibration Assembly (SCA) consists of three targets selectable via a single mechanism which also opens a door. The targets consist of a Daily Bright Target (DBT), Monthly Bright Target (MBT), and Daily Dim Target (DDT). The bright targets are quartz QVDs with the monthly target being used to track the degradation of the daily target. The dim target is used to track CCD linearity using Progressive Time-Delay Integration (PTDI). A composite baffle is attached to the SCA housing aperture to block Earth shine and stray light from the spacecraft. The SCA assembly is mounted to a view port ~90° from OCI nadir. This paper provides an overview of driving solar calibration requirements, error-budgets and early trade studies which drove the solar calibration assembly (SCA) architecture and on-orbit maneuver. Measurements of the diffuser Bidirectional Reflectance Distribution Function (BRDF) at TNO, Netherlands and GSFC are briefly described. Optical modelling and test results at the sub-system and instrument level are included. Finally, preliminary measurements on-orbit are compared to pre-launch predictions.

ocean color↗
NASA NTRS1968

Independent doubly truncated gamma variables

Density and distribution functions of the sum of independent variables, each having a truncated gamma density function, were derived for use in the measurement of complex physical phenomena.

Lavender, D. E.↗