Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “probability density 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.

699 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Alternative Perspectives on Risk

The goal of the commercial air transport system is to provide air transportation to the flying public at an acceptable cost with minimal risk. in an ideal situation these three goals would support each other. In fact, it is sometimes the case that the goals conflict: getting passengers to their destinations on time may conflict with fixing a minor mechanical malfunction that may or may not impact safety; flying a route that will avoid turbulence, thereby providing passengers with a more comfortable ride, may consume more fuel; managing traffic density may mean aircraft are delayed or must use an approach that will result in a long taxi to their gates, costing time and fuel. Various players in the system--pilots, dispatchers, controllers, as well as managers in the airline carriers and traffic management system--make decisions every day that involve trade-offs of benefits and costs. The prospect of revisions in the air traffic management system, with shifts in responsibilities from controllers to users, including airline operations center personnel and pilots, means that individuals may be performing either new jobs or old jobs under new guidance. It will be essential to know how the various players (a) perceive the risks and benefits associated with the decisions they will make under the old and new control structures, and (b) how much risk they are willing to accept in making decisions. Risk is here defined as the probability and magnitude of negative events (after Slovic, 1987). Of primary interest are risks associated with traffic, weather, and operational factors such as schedule, fuel consumption, and passenger service. Previous research has documented differences between groups in perceptions of risks associated with both everyday and aviation related situations. Risk perception varies as a function of familiarity with the situation, degree to which one is potentially affected by the risk, the level of control one has over the situation, and one's level of experience and responsibility in the situation. In our presentation we will consider several factors that may influence differences in risk perception, risk tolerance and risk management among the three major categories of participants in the aviation system (pilots, dispatchers, and controllers). Primary factors that may affect risk attitudes (the collective term we will use for the three components of risk) include: expertise or job-specific knowledge or training, personal involvement or vulnerability to consequences, goals, time horizon (imminence of consequences), span of control and type of control, and information or technology support. These will be considered in the context of five different types of risk (physical threat, economic, social, legal, and ethical). A study designed to study differences between pilots and controllers in their perceptions and responses to traffic risks under the present and a user-preferred control situation will be described. Future studies involving dispatchers' risk perceptions regarding various types of risk will be discussed.

Davison, Jeannie↗
NASA NTRS2013

NASA Tech Briefs, May 2013

Topics include: Test Waveform Applications for JPL STRS Operating Environment; Pneumatic Proboscis Heat-Flow Probe; Method to Measure Total Noise Temperature of a Wireless Receiver During Operation; Cursor Control Device Test Battery; Functional Near-Infrared Spectroscopy Signals Measure Neuronal Activity in the Cortex; ESD Test Apparatus for Soldering Irons; FPGA-Based X-Ray Detection and Measurement for an X-Ray Polarimeter; Sequential Probability Ratio Test for Spacecraft Collision Avoidance Maneuver Decisions; Silicon/Carbon Nanotube Photocathode for Splitting Water; Advanced Materials and Fabrication Techniques for the Orion Attitude Control Motor; Flight Hardware Packaging Design for Stringent EMC Radiated Emission Requirements; RF Reference Switch for Spaceflight Radiometer Calibration; An Offload NIC for NASA, NLR, and Grid Computing; Multi-Scale CNT-Based Reinforcing Polymer Matrix Composites for Lightweight Structures; Ceramic Adhesive and Methods for On-Orbit Repair of Re-Entry Vehicles; Self-Healing Nanocomposites for Reusable Composite Cryotanks; Pt-Ni and Pt-Co Catalyst Synthesis Route for Fuel Cell Applications; Aerogel-Based Multilayer Insulation with Micrometeoroid Protection; Manufacturing of Nanocomposite Carbon Fibers and Composite Cylinders; Optimized Radiator Geometries for Hot Lunar Thermal Environments; A Mission Concept: Re-Entry Hopper-Aero-Space-Craft System on-Mars (REARM-Mars); New Class of Flow Batteries for Terrestrial and Aerospace Energy Storage Applications; Reliability of CCGA 1152 and CCGA 1272 Interconnect Packages for Extreme Thermal Environments; Using a Blender to Assess the Microbial Density of Encapsulated Organisms; Mixed Integer Programming and Heuristic Scheduling for Space Communication; Video Altimeter and Obstruction Detector for an Aircraft; Control Software for Piezo Stepping Actuators; Galactic Cosmic Ray Event-Based Risk Model (GERM) Code; Sasquatch Footprint Tool; and Multi-User Space Link Extension (SLE) System.

Source record↗
NASA NTRS1994

Photoionization of the diffuse interstellar medium and galactic halo by OB associtations

Assuming smoothly varying H I distributions in te Galactic disk, we have calculated the geometry of diffuse II regions due to OB associations in the Galactic plane. Near the solar circle, OB associations with a Lyman continuum (Lyc) photon luminosity Psi(sub Lyc) = 3.3 x 10(exp 7) cm(exp -2) s(exp -1), produce H II regions that are density bounded in the vertical direction (H II chimneys) allowing Lyc to escape the gaseous disk and penetrate into the Galactic halo. We provide analytic formulae for the Lyc escape fraction as functions of S(sub 0) O-star catalog of Garmany and a new Lyc stellar stellar Lyc stellar flux calibration, we find a production rate of Lyc photons by OB associations within 2.5 kpc of Psi(sub Lyc) = 3.3 x 10(exp 7) cm(exp -2) s(exp -1). Integrating the fraction of Lyc photons that escape the disk over our adopted luminosity function of OB associations, we estimate that approximately 7% of the ionizing photons, or Phi(sub Lyc) = 2.3 x 10(exp 6) cm(exp -2) s(exp -1), escape each side of the H I disk layer and penetrate the diffuse ionized medium ('Reynolds layer'). This flux is sufficient to explain the potoionization of this, although we have not constructed a model for the observed H-alpha emission and pulsar dispersion measures that is fully consistent with the absorption rate of Lyc in the H II layer. Since our quiescent model does not account for the effects of dynamic chimneys and superbubbles, which should enhance Lyc escape, we conclude the O stars are the probable source of ionizing radiation for the Reynolds layer. For a random distribution of OB associations throughout the disk, the Lyc flux is nearly uniform for heights Z is greater than approximately 0.8 kpc above the midplane.

Dove, James B.↗
NASA NTRS1998

HST Spatially Resolved Spectra of the Accretion Disc and Gas Stream of the Nova-Like Variable UX Ursae Majoris

Time-resolved eclipse spectroscopy of the nova-like variable UX UMa obtained with the Hubble Space Telescope/Faint Object Spectrograph (HST/FOS) on 1994 August and November is analysed with eclipse mapping techniques to produce spatially resolved spectra of its accretion disk and gas stream as a function of distance from the disk centre. The inner accretion disk is characterized by a blue continuum filled with absorption bands and lines, which cross over to emission with increasing disk radius, similar to that reported at optical wavelengths. The comparison of spatially resolved spectra at different azimuths reveals a significant asymmetry in the disk emission at ultraviolet (UV) wavelengths, with the disk side closest to the secondary star showing pronounced absorption by an 'iron curtain' and a Balmer jump in absorption. These results suggest the existence of an absorbing ring of cold gas whose density and/or vertical scale increase with disk radius. The spectrum of the infalling gas stream is noticeably different from the disc spectrum at the same radius suggesting that gas overflows through the impact point at the disk rim and continues along the stream trajectory, producing distinct emission down to 0.1 R(sub LI). The spectrum of the uneclipsed light shows prominent emission lines of Lyalpha, N v lambda1241, SiIV Lambda 1400, C IV Lambda 1550, HeII Lambda 1640, and MgII Lambda 2800, and a UV continuum rising towards longer wavelengths. The Balmer jump appears clearly in emission indicating that the uneclipsed light has an important contribution from optically thin gas. The lines and optically thin continuum emission are most probably emitted in a vertically extended disk chromosphere + wind. The radial temperature profiles of the continuum maps are well described by a steady-state disc model in the inner and intermediate disk regions (R greater than or equal to 0.3R(sub LI) ). There is evidence of an increase in the mass accretion rate from August to November (from V = 10 (exp -8.3 +/-0.1) to 10(exp -8.1 +/- 0.1 solar mass yr(exp -1)), in accordance with the observed increase in brightness. Since the UX UMA disc seems to be in a high mass accretion, high-viscosity regime in both epochs, this result suggests that the mass transfer rate of UX UMA varies substantially (approximately equal to 50 per cent) on time-scales of a few months. It is suggested that the reason for the discrepancies between the prediction of the standard disk model and observations is not an inadequate treatment of radiative transfer in the disc atmosphere, but rather the presence of addition important sources of light in the system besides the accretion disk (e.g., optically thin contiuum emission from the disk wind and possible absorption by circumstellar cool gas).

Baptista, Raymundo↗
NASA NTRS2004

Risk Assessment and Control through Countermeasure System Iplementation for Long-term Crew Exposure to Microgravity

Experience with the International Space Station (ISS) program demonstrates the degree to which engineering design and operational solutions must protect crewmembers from health risks due to long-term exposure to the microgravity environment. Risks to safety and health due to degradation in the microgravity environment include crew inability to complete emergency or nominal activities, increased risk of injury, and inability to complete safe return to the ground due to reduced strength or embrittled bones. These risks without controls slowly increase in probability for the length of the mission and become more significant for increasing mission durations. Countermeasures to microgravity include hardware systems that place a crewmember s body under elevated stress to produce an effect similar to daily exposure to gravity. The ISS countermeasure system is predominately composed of customized exercise machines. Historical treatment of microgravity countermeasure systems as medical research experiments unintentionally reduced the foreseen importance and therefore the capability of the systems to function in a long-term operational role. Long-term hazardous effects and steadily increasing operational risks due to non-functional countermeasure equipment require a more rigorous design approach and incorporation of redundancy into seemingly non- mission-critical hardware systems. Variations in the rate of health degradation and responsiveness to countermeasures among the crew population drastically increase the challenge for design requirements development and verification of the appropriate risk control strategy. The long-term nature of the hazards and severe limits on logistical re-supply mass, volume and frequency complicates assessment of hardware availability and verification of an adequate maintenance and sparing plan. Design achievement of medically defined performance requirements by microgravity countermeasure systems and incorporation of adequate failure tolerance significantly reduces these risks. Future implementation of on-site monitoring hardware for critical health parameters such as bone mineral density would allow greater responsiveness, efficiency, and optimized design of the countermeasures system.

Gernand, Jeremy M.↗
NASA NTRS1997

Improvement in Suppression of Pulsed Nd:YAG Laser Light With Iodine Absorption Cells for Filtered Rayleigh Scattering Measurements

Filtered Rayleigh scattering using iodine absorption cells is an effective technique for obtaining density, temperature, and velocity measurements in high speed confined flows. By tuning a single frequency laser to a strong iodine absorption line, stray scattered laser light can be greatly suppressed. For example, the minimum transmission predicted by an iodine absorption model calculation is less than 10(exp -5) at the 18788.44/cm line using a 200 mm absorption cell containing iodine vapor at 0.46 T. Measurements obtained by other researches using a CW Nd:YAG laser agree with the model calculations. However, measurements made by us and by others using Q-switched, injection-seeded, frequency doubled Nd:YAG lasers only show minimum transmission of about 3 x 10(exp -3). This greatly reduces the applicability of the filtered Rayleigh scattering technique using these lasers in experiments having large amounts of stray scattered laser light. The purposes of the present study are to characterize the spectrum of the excess light transmitted by the iodine cell and to make changes to the laser to reduce the transmitted laser light. Transmission data as a function of laser frequency for the iodine absorption line at 18788.44/cm are presented. A planar mirror Fabry-Perot interferometer was used to characterize the frequency spectrum of the light passed through the cell. Measurements taken with the laser tuned to the center of the iodine absorption line show the light transmitted through the iodine cell to have a component with a bandwidth of about 40 GHz. This is probably caused by other modes in the laser that exist in spite of the single frequency injection beam. A second broadband component was also observed, possibly caused by the laser flash lamps or by fluorescence. An intracavity etalon was installed in the laser oscillator cavity to suppress the 40 GHz component. Measurements taken with the etalon tuned to the injection frequency showed a reduction in the transmitted laser light. This improvement allows the iodine cell to block significantly more of the stray laser light in filtered Rayleigh scattering experiments. Examples are given of filtered Rayleigh scattering measurements showing the effect of the etalon on measurements taken in a Mach 3 flow in the NASA Lewis 4 inch by 10 inch supersonic wind tunnel.

Seasholtz, Richard G.↗
NASA NTRS2007

Titan Orbiter Aerorover Mission with Enceladus Science (TOAMES)

Cassini and Huygens have made exciting discoveries at Titan and Enceladus, and at the same time made us aware of how little we understand about these bodies. For example, the source, and/or recycling mechanism, of methane in Titan's atmosphere is still puzzling. Indeed, river beds (mostly dry) and lakes have been spotted, and occasional clouds have been seen, but the physics to explain the observations is still mostly lacking, since our "image" of Titan is still sketchy and quite incomplete. Enceladus, only -500 km in extent, is even more puzzling, with its fiery plumes of vapor, dust and ice emanating from its south polar region, "feeding" Saturn's E ring. Long term variability of magnetospheric plasma, neutral gas, E-ring ice grain density, radio emissions, and corotation of Saturn's planetary magnetic field in response to Enceladus plume activity are of great interest for Saturn system science. Both Titan and Enceladus are bodies of considerable astrobiological interest in view of high organic abundances at Titan and potential subsurface liquid water at Enceladus. We propose to develop a new mission to Titan and Enceladus, the Titan Orbiter Aerorover Mission with Enceladus Science (TOAMES), to address these questions using novel new technologies. TOAMES is a multi-faceted mission that starts with orbit insertion around Saturn using aerobraking with Titan's extended atmosphere. We then have an orbital tour around Saturn (for 1-2 years) and close encounters with Enceladus, before it goes into orbit around Titan (via aerocapture). During the early reconnaissance phase around Titan, perhaps 6 months long, the orbiter will use altimetry, radio science and remote sensing instruments to measure Titan's global topography, subsurface structure and atmospheric winds. This information will be used to determine where and when to release the Aerorover, so that it can navigate safely around Titan and identify prime sites for surface sampling and analysis. In situ instruments will sample the upper atmosphere which may provide the seed population for the complex organic chemistry on the surface. The Aerorover will probably use a "hot air" Montgolfier balloon concept using the waste heat from the MMRTG 1-2 kwatts. New technologies will need to be developed and miniaturization will be required to maintain functionality while controlling mass, power and cost. Duty cycling will be used. The Aerorover will have all the instruments needed to sample Titan's atmosphere and surface with possible methane lakes-rivers. It will e.g., use multi-spectral imagers and for last 6 months of mission, balloon payload will land on surface at predetermined site to take core samples of the surface and use seismometers to help probe the interior. All remote (and active) sensors on the orbiter will share a - 1 meter telescope, called MIDAS (Multiple Instrument Distributed Aperture Sensor). MIDAS observations in stable orbit at Titan can provide full global maps of Titan's surface and could additionally provide long term observations of the Saturn system including Enceladus for extended mission phases over many years, potentially for decades. Experience from the Hubble Space Telescope has shown strong public interest and commitment to exciting generational missions.

Sittler, Edward C.↗
NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗
NASA NTRS2024

Quantification of Uncertainty and Risk Sensitivity for Safety of Emerging Operations

The growing need to develop and deploy small unmanned aerial vehicles (sUAVs) for various applications in the airspace necessitates reliable tools to accurately predict the flight trajectories of the sUAVs. The knowledge of the predicted trajectories help decision makers anticipate potential conflict, assess the risk, and take appropriate risk mitigation actions. In addition, uncertainties in vehicle models, weather, and controller action further highlights the need for reliable prediction tools. In this project, the application of mixed sparse grid-based quadrature and generalized polynomial chaos(gPC) expansion method for uncertainty quantification and collision assessment in air traffic consisting of fixed-wing small unmanned aerial vehicles (sUAV) was studied. From the results obtained, it can be concluded that this provides a reliable framework to carry out quantitative conflict assessment in an unmanned air traffic, which when employed, can improve the functionalities of the unmanned traffic management system. It was observed that the results from the gPC expansion framework developed in the project can be utilized to conduct rapid probabilistic collision assessment for near real-time unmanned traffic management in the airspace. From the vehicle models, position updates, and wind-field data, a priori gPC based 3-σcon-fidence ellipses can provide estimates of potential conflict at some future instants. The computational costs scaled linearly when the uncertain inputs were fewer. Further, the largest allowable distribution of para-metric uncertainties that leads to the smallest risk of collision in traffic of small unmanned aerial vehicles could be calculated. The time of closest approach between two sUAVs can be established paving way for development of proactive mitigation strategies. The separation between the sUAVs was found to be most significantly affected by uncertainties in the maximum available thrusts, zero-lift drag coefficients, and wing planform areas of the sUAVs. The study of uncertain wind-fields indicated that a heterogeneous traffic mix resulted in an increased probability of conflict. Increased measurement update rate reduced the uncertain-ties in the trajectories of the vehicles, further reducing the probability of conflict but rapid updates of all vehicles in the airspace poses a stringent communication limitation. The gPC framework also provided the means to analyze vehicle impact (crash region) due to loss of control resulting from actuator failure in sUAS traffic, essentially to predict impact and crash zones for representative vehicles. The predicted regions when compared with non-participant density, provides a means to develop an early mitigation strategy, should the sUAV detect an imminent actuator failure.

Rajnish Bhusal↗
NASA NTRS2008

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Aerosols and especially their effect on clouds are one of the key components of the climate system and the hydrological cycle [Ramanathan et al., 2001]. Yet, the aerosol effect on clouds remains largely unknown and the processes involved not well understood. A recent report published by the National Academy of Science states "The greatest uncertainty about the aerosol climate forcing - indeed, the largest of all the uncertainties about global climate forcing - is probably the indirect effect of aerosols on clouds [NRC, 2001]." The aerosol effect on clouds is often categorized into the traditional "first indirect (i.e., Twomey)" effect on the cloud droplet sizes for a constant liquid water path [Twomey, 1977] and the "semi-direct" effect on cloud coverage [e.g., Ackerman et al ., 2001]." Enhanced aerosol concentrations can also suppress warm rain processes by producing a narrow droplet spectrum that inhibits collision and coalescence processes [e.g., Squires and Twomey, 1961; Warner and Twomey, 1967; Warner, 1968; Rosenfeld, 19991. The aerosol effect on precipitation processes, also known as the second type of aerosol indirect effect [Albrecht, 1989], is even more complex, especially for mixed-phase convective clouds. Table 1 summarizes the key observational studies identifying the microphysical properties, cloud characteristics, thermodynamics and dynamics associated with cloud systems from high-aerosol continental environments. For example, atmospheric aerosol concentrations can influence cloud droplet size distributions, warm-rain process, cold-rain process, cloud-top height, the depth of the mixed phase region, and occurrence of lightning. In addition, high aerosol concentrations in urban environments could affect precipitation variability by providing an enhanced source of cloud condensation nuclei (CCN). Hypotheses have been developed to explain the effect of urban regions on convection and precipitation [van den Heever and Cotton, 2007 and Shepherd, 2005]. Please see Tao et al. (2007) for more detailed description on aerosol impact on precipitation. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated.

Tao, Wei-Kuo↗
NASA NTRS2013

Enigma Version 12

Enigma Version 12 software combines model building, animation, and engineering visualization into one concise software package. Enigma employs a versatile user interface to allow average users access to even the most complex pieces of the application. Using Enigma eliminates the need to buy and learn several software packages to create an engineering visualization. Models can be created and/or modified within Enigma down to the polygon level. Textures and materials can be applied for additional realism. Within Enigma, these models can be combined to create systems of models that have a hierarchical relationship to one another, such as a robotic arm. Then these systems can be animated within the program or controlled by an external application programming interface (API). In addition, Enigma provides the ability to use plug-ins. Plugins allow the user to create custom code for a specific application and access the Enigma model and system data, but still use the Enigma drawing functionality. CAD files can be imported into Enigma and combined to create systems of computer graphics models that can be manipulated with constraints. An API is available so that an engineer can write a simulation and drive the computer graphics models with no knowledge of computer graphics. An animation editor allows an engineer to set up sequences of animations generated by simulations or by conceptual trajectories in order to record these to highquality media for presentation. Enigma Version 12 Lyndon B. Johnson Space Center, Houston, Texas 28 NASA Tech Briefs, September 2013 Planetary Protection Bioburden Analysis Program NASA's Jet Propulsion Laboratory, Pasadena, California This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware.

Shores, David↗
NASA NTRS1998

Angular Broadening: Effects of Nonzero, Spatially Varying Plasma Frequency Between the Source and Observer

Angular broadening of radiation due to scattering by density irregularities is usually described using geometric optics (GO) or the parabolic wave equation (PWE) with the assumptions that the radiation frequency f greatly exceeds the local plasma frequency f(sub p0) or that f(sub p0)/f is constant along the path. These assumptions are inappropriate for many solar system radio phenomena. Here the PWE and GO formalisms are extended to treat angular broadening in plasmas with nonzero, spatially varying ratios, f(sub p0)(z)/f < 1. The new PWE results show that the correlation function, scattered angular spectrum, and other quantities are modified by inclusion of a denominator factor [1 - f(sup 2, sub p0)(z prime)/f prime] inside the path integral over z prime, while the mean-square scattering angle depends on both f(sub p0)(z)/f at the observer and the foregoing factor inside the path integral. The PWE and GO predictions for are identical and involve equivalent assumptions. Previous GO and PWE results are recovered in the limits that f(sub p0)(z prime)/f is constant or zero. The new PWE and GO results will permit more accurate calculation of angular broadening for solar system and astrophysical sources. Moreover and importantly, due to the PWE and GO results for being identical, previous GO analyses of in solar system contexts are essentially correct, except for the neglect of or minor deficiencies in the treatment of nonzero, spatially varying f(sub p0)/f effects. The identical GO and PWE results for and the form of the PWE equation for the correlation function raise questions as to whether diffraction is unimportant for angular broadening (under the usual PWE conditions). Future direct comparisons of the PWE predictions with angular spectra calculated using existing GO ray-tracing codes should answer these questions. Diffraction effects are probably important when the medium and turbulence are not sufficiently homogeneous transverse to the central ray path. The implications are also discussed for studies of the 2-3.5 kHz radiation observed in the outer heliosphere.

Cairns, Iver H.↗
NASA NTRS2000

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗
NASA NTRS1996

Early Energetic Particle Irradiation of the HED Parent Body Regolith

Previous studies have shown that many individual grains within the dark phase of the Kapoeta howardite were irradiated with energetic particles while residing on the surface of the early HED regolith. Particle tracks in these grains vary in density by more than an order of magnitude and undoubtedly were formed by energetic heavy (Fe) ions associated with early solar flares. Early Irradiation of HED Regolith: Concentrations of excess Ne alone are not sufficient to decide between competing galactic and solar irradiation models. However, from recent studies of depth samples of oriented lunar rocks, we have shown that the cosmogenic 21-Ne/22-Ne ratio produced in feldspar differs substantially between Galactic Cosmic Radiation (GCR) and solar protons, and that this difference is exactly that predicted from cross-section data. Using Ne literature data and new isotopic data we obtained on acid-etched, separated feldspar from both the light and dark phases of Kapoeta, we derive 21-Ne/22-Ne = 0.80 for the recent GCR irradiation and 21-Ne/22-Ne = 0.68 for the early regolith irradiation. This derived ratio indicates that the early Ne production in the regolith occurred by both galactic and solar protons. If we adopt a likely one-component regolith model in which all grains were exposed to galactic protons but individual grains had variable exposure to solar protons, we estimate that this early GCR irradiation lasted for about 3-6 m.y. More complex two-component regolith models involving separate solar and galactic irradiation would permit this GCR age to be longer. Higher-energy solar protons would permit the GCR to be longer. Higher-energy solar protons would permit the GCR age to be shorter. Further, cosmogenic 126(Xe) in Kapoeta dark is no more than a factor of about 2 higher than that observed in Kapoeta light. Because 126(Xe) can only be formed by galactic protons and not solar protons, these data support a short GCR irradiation for the HED regolith. This would also be the maximum time peRiod for the solar irradiation. Various asteroidal regolith models, based on Monte Carlo modeling of impact rates as a function of size and on irradiation features of meteorites, predict surface exposure times of about 0.1 to 10 m.y., and depend on such factors as gravity, rock mechanical properties, and micrometeoroid flux. Because the depth at which solar Fe tracks are produced (is much less than 1 micrometer) is much less than the depth at which Solar Cosmic Rays (SCR) Ne is produced (about 1 cm), for a reasonably well-stirred HED regolith the "surface exposure time" for SCR 21-Ne production should be significantly longer than that for solar tracks and some other surface irradiation features. Enhanced Solar Proton Irradiation: For bulk samples of Kapoeta dark feldspar and a one-component regolith model, the derived ratio of 21-Ne/22-Ne = 0.68 implies that the early production ratio of SCR 21-Ne to GCR 21-Ne was about 0.5-1.5. This ratio is independent of any assumptions about the fraction of dark grains that are irradiated or of the variability in the degree of solar irradiation among grains. The 21-Ne SCR/GCR ratio indirectly derived from bulk Kapoeta pyroxene is somewhat larger, as is the ratio derived for simple two-component regolith models. Individual feldspar grains that were extensively solar irradiated would require even larger 21-Ne SCR/GCR production ratios. In contrast, the theoretical SCR/GCR production ratio for lunar feldspar with 0 g/CM2 shield ing is is less than or equal to 2, and the lowest ratio observed in near-surface samples of lunar anorthosites is less than or equal to 1. Considering the greater solar distance of Vesta (compared to the Moon), the likelihood that SCR 21-Ne was acquired under some shielding where production rates are lower, and the likelihood that the exposure time to galactic protons exceeded the exposure time to solar protons because of their very different penetration depths, the 21-Ne SCR/GCR production ratio on the HED parent body was probably < 0.1. The relatively large difference between the derived 21-Ne SCR/GCR ratio in Kapoeta dark feldspar and the estimated production ratio strongly indicates that the early solar irradiation involved a flux -20-50x the recent solar flux. This enhanced proton flux was probably associated with an overall greater solar activity in the first approximately 10(exp 7) to 10(exp 8) years of solar history.

Bogard, D. D.↗
NASA NTRS2013

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗