Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse modeling”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

On the Location of the gamma-Ray Outburst Emission in the BL Lacertae Object AO 0235 + 164 Through Observations Across the Electromagnetic Spectrum

We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array images at A = 7 mm with approx.0.15 milliarcsec resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7 mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary "core" and in the superluminal knot, both parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long timescales (months/ years), but the correspondence is poorer on shorter timescales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes.

Agudo, Ivan↗
NASA NTRS1996

SSME Condition Monitoring Using Neural Networks and Plume Spectral Signatures

For a variety of reasons, condition monitoring of the Space Shuttle Main Engine (SSME) has become an important concern for both ground tests and in-flight operation. The complexities of the SSME suggest that active, real-time condition monitoring should be performed to avoid large-scale or catastrophic failure of the engine. In 1986, the SSME became the subject of a plume emission spectroscopy project at NASA's Marshall Space Flight Center (MSFC). Since then, plume emission spectroscopy has recorded many nominal tests and the qualitative spectral features of the SSME plume are now well established. Significant discoveries made with both wide-band and narrow-band plume emission spectroscopy systems led MSFC to develop the Optical Plume Anomaly Detection (OPAD) system. The OPAD system is designed to provide condition monitoring of the SSME during ground-level testing. The operational health of the engine is achieved through the acquisition of spectrally resolved plume emissions and the subsequent identification of abnormal emission levels in the plume indicative of engine erosion or component failure. Eventually, OPAD, or a derivative of the technology, could find its way on to an actual space vehicle and provide in-flight engine condition monitoring. This technology step, however, will require miniaturized hardware capable of processing plume spectral data in real-time. An objective of OPAD condition monitoring is to determine how much of an element is present in the SSME plume. The basic premise is that by knowing the element and its concentration, this could be related back to the health of components within the engine. For example, an abnormal amount of silver in the plume might signify increased wear or deterioration of a particular bearing in the engine. Once an anomaly is identified, the engine could be shut down before catastrophic failure occurs. Currently, element concentrations in the plume are determined iteratively with the help of a non-linear computer code called SPECTRA, developed at the USAF Arnold Engineering Development Center. Ostensibly, the code produces intensity versus wavelength plots (i.e., spectra) when inputs such as element concentrations, reaction temperature, and reaction pressure are provided. However, in order to provide a higher-level analysis, element concentration is not specified explicitly as an input. Instead, two quantum variables, number density and broadening parameter, are used. Past experience with OPAD data analysis has revealed that the region of primary interest in any SSME plume spectrum lies in the wavelength band of 3300 A to 4330 A. Experience has also revealed that some elements, such as iron, cobalt and nickel, cause multiple peaks over the chosen wavelength range whereas other elements (magnesium, for example) have a few, relatively isolated peaks in the chosen wavelength range. Iteration with SPECTRA as a part of OPAD data analysis is an incredibly labor intensive task and not one to be performed by hand. What is really needed is the "inverse" of the computer code but the mathematical model for the inverse mapping is tenuous at best. However, building generalized models based upon known input/output mappings while ignoring details of the governing physical model is possible using neural networks. Thus the objective of the research project described herein was to quickly and accurately predict combustion temperature and element concentrations (i.e., number density and broadening parameter) from a given spectrum using a neural network. In other words, a neural network had to be developed that would provide a generalized "inverse" of the computer code SPECTRA.

Hopkins, Randall↗
NASA NTRS2007

NASA Tech Briefs, October 2007

Topics covered include; Wirelessly Interrogated Position or Displacement Sensors; Ka-Band Radar Terminal Descent Sensor; Metal/Metal Oxide Differential Electrode pH Sensors; Improved Sensing Coils for SQUIDs; Inductive Linear-Position Sensor/Limit-Sensor Units; Hilbert-Curve Fractal Antenna With Radiation- Pattern Diversity; Single-Camera Panoramic-Imaging Systems; Interface Electronic Circuitry for an Electronic Tongue; Inexpensive Clock for Displaying Planetary or Sidereal Time; Efficient Switching Arrangement for (N + 1)/N Redundancy; Lightweight Reflectarray Antenna for 7.115 and 32 GHz; Opto-Electronic Oscillator Using Suppressed Phase Modulation; Alternative Controller for a Fiber-Optic Switch; Strong, Lightweight, Porous Materials; Nanowicks; Lightweight Thermal Protection System for Atmospheric Entry; Rapid and Quiet Drill; Hydrogen Peroxide Concentrator; MMIC Amplifiers for 90 to 130 GHz; Robot Would Climb Steep Terrain; Measuring Dynamic Transfer Functions of Cavitating Pumps; Advanced Resistive Exercise Device; Rapid Engineering of Three-Dimensional, Multicellular Tissues With Polymeric Scaffolds; Resonant Tunneling Spin Pump; Enhancing Spin Filters by Use of Bulk Inversion Asymmetry; Optical Magnetometer Incorporating Photonic Crystals; WGM-Resonator/Tapered-Waveguide White-Light Sensor Optics; Raman-Suppressing Coupling for Optical Parametric Oscillator; CO2-Reduction Primary Cell for Use on Venus; Cold Atom Source Containing Multiple Magneto- Optical Traps; POD Model Reconstruction for Gray-Box Fault Detection; System for Estimating Horizontal Velocity During Descent; Software Framework for Peer Data-Management Services; Autogen Version 2.0; Tracking-Data-Conversion Tool; NASA Enterprise Visual Analysis; Advanced Reference Counting Pointers for Better Performance; C Namelist Facility; and Efficient Mosaicking of Spitzer Space Telescope Images.

Source record↗
NASA NTRS2015

Design and Demonstration of the Mars Organic Molecule Analyzer (MOMA) on the ExoMars 2018 Rover

The Mars Organic Molecule Analyzer (MOMA) investigation is a key astrobiology experiment scheduled to launch on the joint ESA-Roscosmos ExoMars 2018 rover mission. MOMA will examine the chemical composition of geological samples acquired from depths of up to two meters below the martian surface, where fragile organic molecules may be protected from destructive cosmic radiation and/or oxidative chemical reactions. The heart of the MOMA mass spectrometer subsystem (i.e., MOMA-MS) is a miniaturized linear ion trap (LIT) that supports two distinct modes of operation to detect: i) volatile and semi-volatile, low-to-moderate mass organics (less than or equal to 500 Da) via pyrolysis coupled with gas chromatography mass spectrometry (pyr/GCMS); and, ii) more refractory, moderate-to-high mass compounds (up to 1000 Da) via laser desorption (LDMS) at ambient Mars pressures. Additionally, the LIT mass analyzer enables selective ion trapping via multi-frequency waveform ion excitation (e.g., stored waveform inverse Fourier transform, or SWIFT), and structural characterization of complex molecules using tandem mass spectrometry (MSIMS). A high-fidelity Engineering Test Unit (ETU) of MOMAMS, including the LIT subassembly, dual-gun electron ionization source, micropirani pressure gauge, solenoid-driven aperture valve, redundant detection chains, and control electronics, has been built and tested at NASA GSFC under relevant operational conditions (pressure, temperature, etc.). Spaceflight qualifications of individual hardware components and integrated subassemblies have been validated through vibration, shock, thermal, lifetime, and performance evaluations. The ETU serves as a pathfinder for the flight model buildup, integration and test, as the ETU meets the form, fit and function of the flight unit that will be delivered to MPS in late 2015. To date, the ETU of MOMA-MS has been shown to meet or exceed all functional requirements, including mass range, resolution, accuracy, instrumental drift, and Iimit-of-detection specifications, thereby enabling the primary science objectives of the MOMA investigation and ExoMars 2018 mission.

Mars Organic Molecule Analyzer↗
NASA NTRS1995

High temperature kinetic study of the reactions H + O2 = OH + O and O + H2 = OH + H in H2/O2 system by shock tube-laser absorption spectroscopy

The reactions: (1) H + O2 = OH + O; and (2) O + H2 = OH + H are the most important elementary reactions in gas phase combustion. They are the main chain-branching reaction in the oxidation of H2 and hydrocarbon fuels. In this study, rate coefficients of the reactions and have been measured over a wide range of composition, pressure, density and temperature behind the reflected shock waves. The experiments were performed using the shock tube - laser absorption spectroscopic technique to monitor OH radicals formed in the shock-heated H2/O2/Ar mixtures. The OH radicals were detected using the P(1)(5) line of (0,0) band of the A(exp 2) Sigma(+) from X(exp 2) Pi transition of OH at 310.023 nm (air). The data were analyzed with the aid of computer modeling. In the experiments great care was exercised to obtain high time resolution, linearity and signal-to-noise. The results are well represented by the Arrhenius expressions. The rate coefficient expression for reaction (1) obtained in this study is k(1) = (7.13 +/- 0.31) x 10(exp 13) exp(-6957+/- 30 K/T) cu cm/mol/s (1050 K less than or equal to T less than or equal to 2500 K) and a consensus expression for k(1) from a critical review of the most recent evaluations of k(1) (including our own) is k(1) = 7.82 x 10(exp 13) exp(-7105 K/T) cu cm/mol/s (960 K less than or equal to T less than or equal to 5300 K). The rate coefficient expression of k(2) is given by k(2) = (1.88 +/- 0.07) x 10(exp 14) exp(-6897 +/- 53 K/T) cu cm/mol/s (1424 K less than or equal to T less than or equal to 2427 K). For k(1), the temperature dependent A-factor and the correlation between the values of k(1) and the inverse reactant densities were not found. In the temperature range of this study, non-Arrhenius expression of k(2) which shows the upward curvature was not supported.

Ryu, Si-Ok↗
NASA NTRS2006

M(sub W) = 7.2-7.4 Estimated for A.D. 900 Seattle Fault Earthquake by Modeling the Uplift of a Lidar-Mapped Marine Terrace

Inverse modeling of slip on the Seattle fault system, constrained by elevations of uplifted marine terraces, provides a well-constrained estimate of the magnitude of the largest known upper-crust earthquake in the Puget Sound region within the past 2500 years. The terrace elevations that constrain the slip inversion are extracted from elevation and slope images generated from LIDAR surveys of the Puget Sound collected in 1996-2002. The images reveal a single uplifted terrace, dated to 1000 cal yr B.P. near Restoration Point, which is morphologically continuous along the southern shoreline of Bainbridge Island and is visible at comparable elevations within a 25 km by 12 km region encompassing coastlines of West Seattle, Bremerton, East Bremerton, Port Orchard, and Waterman Point. Considering sea level changes since A.D. 900, the maximum uplift magnitudes of shoreline inner edges approach 9 m and are located at the southernmost coastline of Bainbridge Island and the northern tip of Waterman Point, while tilt magnitudes are modest - approaching 0.1 degrees. For each of several different Seattle fault geometry interpretations, we use a linear inversion code to solve for distributed slip on the fault surfaces. Moment magnitudes of 7.2 to 7.4 are calculated directly from the different slip solutions. In general, the greatest slip of the A.D. 900 event was confined to the frontal thrust of the Seattle fault system and was centered beneath Puget Sound between Restoration Point and Alki Point.

Jordan R. Muller↗
NASA NTRS1993

Receptivity in parallel flows: An adjoint approach

Linear receptivity studies in parallel flows are aimed at understanding how external forcing couples to the natural unstable motions which a flow can support. The vibrating ribbon problem models the original Schubauer and Skramstad boundary layer experiment and represents the classic boundary layer receptivity problem. The process by which disturbances are initiated in convectively-unstable jets and shear layers has also received attention. Gaster was the first to handle the boundary layer analysis with the recognition that spatial modes, rather than temporal modes, were relevant when studying convectively-unstable flows that are driven by a time-harmonic source. The amplitude of the least stable spatial mode, far downstream of the source, is related to the source strength by a coupling coefficient. The determination of this coefficient is at the heart of this type of linear receptivity study. The first objective of the present study was to determine whether the various wave number derivative factors, appearing in the coupling coefficients for linear receptivity problems, could be reexpressed in a simpler form involving adjoint eigensolutions. Secondly, it was hoped that the general nature of this simplification could be shown; indeed, a rather elegant characterization of the receptivity properties of spatial instabilities does emerge. The analysis is quite distinct from the usual Fourier-inversion procedures, although a detailed knowledge of the spectrum of the Orr-Sommerfeld equation is still required. Since the cylinder wake analysis proved very useful in addressing control considerations, the final objective was to provide a foundation upon which boundary layer control theory may be developed.

Hill, D. Christopher↗
NASA NTRSDate not supplied

Feasibility Study of a Multi-Tilt-Rotor Aircraft as the Artemis Lunar Training Vehicle

The Lunar Landing Research Vehicles (LLRVs) and the Lunar Landing Training Vehicles (LLTVs) provided astronauts of the Apollo program with essential experience and confidence required to complete the missions, and contributed to six successful manned landings on the moon. The primary challenge in terrestrial training was being able to replicate the ratio of tilt angle to linear acceleration that a pilot would experience in lunar gravity. Presently, as the Artemis program seeks to return humans to the Moon by 2025, engineers are evaluating suitable platforms to serve as an In-Flight Trainer (IFT) or Artemis Lunar Training Vehicle (ALTV) for astronauts training in the task of manual landing. The program is investigating the viability of current technology in the field of electric vertical takeoff and landing (eVTOL) vehicles and is evaluating using a multi-tilt-rotor aircraft platform as a candidate for a preliminary ALTV. The tilt-rotor capability enables the vehicle attitude to be decoupled from its flight path, which is a crucial requirement in realistically simulating lunar gravity on Earth. Other key considerations include compensating for a lack of aerodynamic forces while flying through the atmosphere of Earth, as well as the ability to simulate the dynamics of multiple different lander designs for the Human Landing System (HLS) program. This paper details the feasibility study and presents a preliminary flight control architecture for an IFT based on a notional multi-tilt-rotor platform. The model-following control law, based on nonlinear dynamic inversion (NDI), removes the need for gain scheduling. The inner-loop dynamic control allocation strategy consists of a static portion that is optimized offline for trim while compensating for the difference in gravity and a dynamic portion that is computed in real time. The reference model consists of the full closed-loop dynamics of a generic HLS design. The modularity of the flight control architecture enables evaluation of multiple HLS concepts with minimal modifications to the control law. Simulation results of the multi-tilt-rotor configuration following the final portion of the Apollo 11 descent trajectory are shown.

Jing Pei↗
NASA NTRSDate not supplied

Feasibility Study of a Multi Tilt-rotor Aircraft as the Artemis Lunar Training Vehicle

The Lunar Landing Research Vehicles (LLRVs) and the Lunar Landing Training Vehicles (LLTVs) provided astronaut candidates for the Apollo program with essential experience and confidence required to complete the missions, and contributed to six successful manned landings on the moon. The primary challenge in terrestrial training was being able to replicate the ratio of bank angle to linear acceleration that a pilot would experience in lunar gravity. Presently, as the Artemis program seeks to return humans to the Moon by 2025, engineers are evaluating suitable platforms to serve as an In-Flight Trainer (IFT) or Artemis Lunar Training Vehicle (ALTV) for astronauts training in the task of manual landing. The program is investigating the viability of current technology in the field of electric vertical takeoff and landing (eVTOL) vehicles and is evaluating using a multi tilt-rotor aircraft platform as a candidate platform for a preliminary ALTV. The tilt-rotor capability enables the vehicle attitude to be decoupled from its flight path, which is a crucial requirement in realistically simulating lunar gravity on Earth. Other key considerations include compensating for a lack of aerodynamic forces while flying through the atmosphere of Earth, as well as the ability to simulate the dynamics of multiple different lander designs for the Human Landing System (HLS) program. This paper details the feasibility study and presents a preliminary flight control architecture for an IFT based on a notional multi tilt-rotor platform. The modeling-following control law, based on nonlinear dynamic inversion (NDI), removes the need for gain scheduling because the vehicle operates across a wide range of flight conditions. The inner-loop dynamic control allocation strategy consists of a static portion that is optimized offline for trim while compensating for the difference in gravity and a dynamic portion that is computed in real time. The reference model consists of the full closed-loop dynamics of a generic HLS design. The modularity of the flight control architecture enables evaluation of multiple HLS concepts with minimal modifications to the control law. Simulation results of the multi tilt-rotor configuration following the final portion of the Apollo 11 descent trajectory are shown.

Jing Pei↗
NASA NTRS2001

Accuracy for detection of simulated lesions: comparison of fluid-attenuated inversion-recovery, proton density--weighted, and T2-weighted synthetic brain MR imaging

OBJECTIVE: The objective of our study was to determine the effects of MR sequence (fluid-attenuated inversion-recovery [FLAIR], proton density--weighted, and T2-weighted) and of lesion location on sensitivity and specificity of lesion detection. MATERIALS AND METHODS: We generated FLAIR, proton density-weighted, and T2-weighted brain images with 3-mm lesions using published parameters for acute multiple sclerosis plaques. Each image contained from zero to five lesions that were distributed among cortical-subcortical, periventricular, and deep white matter regions; on either side; and anterior or posterior in position. We presented images of 540 lesions, distributed among 2592 image regions, to six neuroradiologists. We constructed a contingency table for image regions with lesions and another for image regions without lesions (normal). Each table included the following: the reviewer's number (1--6); the MR sequence; the side, position, and region of the lesion; and the reviewer's response (lesion present or absent [normal]). We performed chi-square and log-linear analyses. RESULTS: The FLAIR sequence yielded the highest true-positive rates (p < 0.001) and the highest true-negative rates (p < 0.001). Regions also differed in reviewers' true-positive rates (p < 0.001) and true-negative rates (p = 0.002). The true-positive rate model generated by log-linear analysis contained an additional sequence-location interaction. The true-negative rate model generated by log-linear analysis confirmed these associations, but no higher order interactions were added. CONCLUSION: We developed software with which we can generate brain images of a wide range of pulse sequences and that allows us to specify the location, size, shape, and intrinsic characteristics of simulated lesions. We found that the use of FLAIR sequences increases detection accuracy for cortical-subcortical and periventricular lesions over that associated with proton density- and T2-weighted sequences.

NASA Discipline Neuroscience↗
NASA NTRS1998

A Conceptual Approach to Assimilating Remote Sensing Data to Improve Soil Moisture Profile Estimates in a Surface Flux/Hydrology Model: Overview - Part 1

Knowledge of the amount of water in the soil is of great importance to many earth science disciplines. Soil moisture is a key variable in controlling the exchange of water and energy between the land surface and the atmosphere. Thus, soil moisture information is valuable in a wide range of applications including weather and climate, runoff potential and flood control, early warning of droughts, irrigation, crop yield forecasting, soil erosion, reservoir management, geotechnical engineering, and water quality. Despite the importance of soil moisture information, widespread and continuous measurements of soil moisture are not possible today. Although many earth surface conditions can be measured from satellites, we still cannot adequately measure soil moisture from space. Research in soil moisture remote sensing began in the mid 1970s shortly after the surge in satellite development. Recent advances in remote sensing have shown that soil moisture can be measured, at least qualitatively, by several methods. Quantitative measurements of moisture in the soil surface layer have been most successful using both passive and active microwave remote sensing, although complications arise from surface roughness and vegetation type and density. Early attempts to measure soil moisture from space-borne microwave instruments were hindered by what is now considered sub-optimal wavelengths (shorter than 5 cm) and the coarse spatial resolution of the measurements. L-band frequencies between 1 and 3 GHz (10-30 cm) have been deemed optimal for detection of soil moisture in the upper few centimeters of soil. The Electronically Steered Thinned Array Radiometer (ESTAR), an aircraft-based instrument operating a 1,4 GHz, has shown great promise for soil moisture determination. Initiatives are underway to develop a similar instrument for space. Existing space-borne synthetic aperture radars (SARS) operating at C- and L-band have also shown some potential to detect surface wetness. The advantage of radar is its much higher resolution than passive microwave systems, but it is currently hampered by surface roughness effects and the lack of a good algorithm based on a single frequency and single polarization. In addition, its repeat frequency is generally low (about 40 days). In the meantime, two new radiometers offer some hope for remote sensing of soil moisture from space. The Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI), launched in November 1997, possesses a 10.65 GHz channel and the Advanced Microwave Scanning Radiometer (AMSR) on both the ADEOS-11 and Earth Observing System AM-1 platforms to be launched in 1999 possesses a 6.9 GHz channel. Aside from issues about interference from vegetation, the coarse resolution of these data will provide considerable challenges pertaining to their application. The resolution of TMI is about 45 km and that of AMSR is about 70 km. These resolutions are grossly inconsistent with the scale of soil moisture processes and the spatial variability of factors that control soil moisture. Scale disparities such as these are forcing us to rethink how we assimilate data of various scales in hydrologic models. Of particular interest is how to assimilate soil moisture data by reconciling the scale disparity between what we can expect from present and future remote sensing measurements of soil moisture and modeling soil moisture processes. It is because of this disparity between the resolution of space-based sensors and the scale of data needed for capturing the spatial variability of soil moisture and related properties that remote sensing of soil moisture has not met with more widespread success. Within a single footprint of current sensors at the wavelengths optimal for this application, in most cases there is enormous heterogeneity in soil moisture created by differences in landcover, soils and topography, as well as variability in antecedent precipitation. It is difficult to interpret the meaning of 'mean' soil moisture under such conditions and even more difficult to apply such a value. Because of the non-linear relationships between near-surface soil moisture and other variables of interest, such as surface energy fluxes and runoff, mean soil moisture has little applicability at such large scales. It is for these reasons that the use of remote sensing in conjunction with a hydrologic model appears to be of benefit in capturing the complete spatial and temporal structure of soil moisture. This paper is Part I of a four-part series describing a method for intermittently assimilating remotely-sensed soil moisture information to improve performance of a distributed land surface hydrology model. The method, summarized in section II, involves the following components, each of which is detailed in the indicated section of the paper or subsequent papers in this series: Forward radiative transfer model methods (section II and Part IV); Use of a Kalman filter to assimilate remotely-sensed soil moisture estimates with the model profile (section II and Part IV); Application of a soil hydrology model to capture the continuous evolution of the soil moisture profile within and below the root zone (section III); Statistical aggregation techniques (section IV and Part II); Disaggregation techniques using a neural network approach (section IV and Part III); and Maximum likelihood and Bayesian algorithms for inversely solving for the soil moisture profile in the upper few cm (Part IV).

Crosson, William L.↗
NASA NTRS2013

Support Routines for In Situ Image Processing

This software consists of a set of application programs that support ground-based image processing for in situ missions. These programs represent a collection of utility routines that perform miscellaneous functions in the context of the ground data system. Each one fulfills some specific need as determined via operational experience. The most unique aspect to these programs is that they are integrated into the large, in situ image processing system via the PIG (Planetary Image Geometry) library. They work directly with space in situ data, understanding the appropriate image meta-data fields and updating them properly. The programs themselves are completely multimission; all mission dependencies are handled by PIG. This suite of programs consists of: (1)marscahv: Generates a linearized, epi-polar aligned image given a stereo pair of images. These images are optimized for 1-D stereo correlations, (2) marscheckcm: Compares the camera model in an image label with one derived via kinematics modeling on the ground, (3) marschkovl: Checks the overlaps between a list of images in order to determine which might be stereo pairs. This is useful for non-traditional stereo images like long-baseline or those from an articulating arm camera, (4) marscoordtrans: Translates mosaic coordinates from one form into another, (5) marsdispcompare: Checks a Left Right stereo disparity image against a Right Left disparity image to ensure they are consistent with each other, (6) marsdispwarp: Takes one image of a stereo pair and warps it through a disparity map to create a synthetic opposite- eye image. For example, a right eye image could be transformed to look like it was taken from the left eye via this program, (7) marsfidfinder: Finds fiducial markers in an image by projecting their approximate location and then using correlation to locate the markers to subpixel accuracy. These fiducial markets are small targets attached to the spacecraft surface. This helps verify, or improve, the pointing of in situ cameras, (8) marsinvrange: Inverse of marsrange . given a range file, re-computes an XYZ file that closely matches the original. . marsproj: Projects an XYZ coordinate through the camera model, and reports the line/sample coordinates of the point in the image, (9) marsprojfid: Given the output of marsfidfinder, projects the XYZ locations and compares them to the found locations, creating a report showing the fiducial errors in each image. marsrad: Radiometrically corrects an image, (10) marsrelabel: Updates coordinate system or camera model labels in an image, (11) marstiexyz: Given a stereo pair, allows the user to interactively pick a point in each image and reports the XYZ value corresponding to that pair of locations. marsunmosaic: Extracts a single frame from a mosaic, which will be created such that it could have been an input to the original mosaic. Useful for creating simulated input frames using different camera models than the original mosaic used, and (12) merinverter: Uses an inverse lookup table to convert 8-bit telemetered data to its 12-bit original form. Can be used in other missions despite the name.

Deen, Robert G.↗
NASA NTRS1986

Aerobic Fitness Does Not Contribute to Prediction of Orthostatic Intolerance

Several investigations have suggested that orthostatic tolerance may be inversely related to aerobic fitness (VO (sub 2max)). To test this hypothesis, 18 males (age 29 to 51 yr) underwent both treadmill VO(sub 2max) determination and graded lower body negative pressures (LBNP) exposure to tolerance. VO(2max) was measured during the last minute of a Bruce treadmill protocol. LBNP was terminated based on pre-syncopal symptoms and LBNP tolerance (peak LBNP) was expressed as the cumulative product of LBNP and time (torr-min). Changes in heart rate, stroke volume cardiac output, blood pressure and impedance rheographic indices of mid-thigh-leg initial accumulation were measured at rest and during the final minute of LBNP. For all 18 subjects, mean (plus or minus SE) fluid accumulation index and leg venous compliance index at peak LBNP were 139 plus or minus 3.9 plus or minus 0.4 ml-torr-min(exp -2) x 10(exp 3), respectively. Pearson product-moment correlations and step-wise linear regression were used to investigate relationships with peak LBNP. Variables associated with endurance training, such as VO(sub 2max) and percent body fat were not found to correlate significantly (P is less than 0.05) with peak LBNP and did not add sufficiently to the prediction of peak LBNP to be included in the step-wise regression model. The step-wise regression model included only fluid accumulation index leg venous compliance index, and blood volume and resulted in a squared multiple correlation coefficient of 0.978. These data do not support the hypothesis that orthostatic tolerance as measured by LBNP is lower in individuals with high aerobic fitness.

Convertino, Victor A.↗
NASA NTRS1994

Gravity waves in Titan's atmosphere

Scintillations (high frequency variations) observed in the radio signal during the occultation of Voyager 1 by Titan (Hinson and Tyler, 1983) provide information concerning neutral atmospheric density fluctuations on scales on hundreds of meters to a few kilometers. Those seen at altitudes higher than 25 km above the surface were interpreted by Hinson and Tyler as being caused by linear, freely propagating (energy-conserving) gravity waves, but this interpretation was found to be inconsistent with the scintillation data below the 25-km altitude level. Here an attempt is made to interpret the entire scintillation profile between the surface and the 90-km altitude level in terms of gravity waves generated at the surface. Numerical calculations of the density fluctuations caused by two-dimensional, nonhydrostatic, finite-amplitude gravity waves propagating vertically through Titan's atmosphere are performed to produce synthetic scintillation profiles for comparison with the observations. The numerical model accurately treats the effects of wave transience, nonlinearity, and breakdown due to convective instability in the overturned part of the wave. The high-altitude scintillation data were accurately recovered with a freely propagating wave solution, confirming the analytic model of Hinson and Tyler. It is found that the low-altitude scintillation data can be fit by a model where a component of the gravity waves becomes convectively unstable and breaks near the 15 km level. The large-scale structure of the observed scintillation profile in the entire altitude range between 5 and 85 km can be simulated by a model where the freely propagating and breaking waves are forced at the surface simultaneously. Further analysis of the Voyager 1 Titan low-altitude scintillation data, using inversion theory appropriate for strong scattering, could potentially remove some of the ambiguities remaining in this analysis and allow a better determination of the strength and source of the waves.

Friedson, A. James↗
NASA NTRS2011

The Impacts of an Observationally-Based Cloud Fraction and Condensate Overlap Parameterization on a GCM's Cloud Radiative Effect

It has been shown that the details of how cloud fraction overlap is treated in GCMs has substantial impact on shortwave and longwave fluxes. Because cloud condensate is also horizontally heterogeneous at GCM grid scales, another aspect of cloud overlap should in principle also be assessed, namely the vertical overlap of hydrometeor distributions. This type of overlap is usually examined in terms of rank correlations, i.e., linear correlations between hydrometeor amount ranks of the overlapping parts of cloud layers at specific separation distances. The cloud fraction overlap parameter and the rank correlation of hydrometeor amounts can be both expressed as inverse exponential functions of separation distance characterized by their respective decorrelation lengths (e-folding distances). Larger decorrelation lengths mean that hydrometeor fractions and probability distribution functions have high levels of vertical alignment. An analysis of CloudSat and CALIPSO data reveals that the two aspects of cloud overlap are related and their respective decorrelation lengths have a distinct dependence on latitude that can be parameterized and included in a GCM. In our presentation we will contrast the Cloud Radiative Effect (CRE) of the GEOS-5 atmospheric GCM (AGCM) when the observationally-based parameterization of decorrelation lengths is used to represent overlap versus the simpler cases of maximum-random overlap and globally constant decorrelation lengths. The effects of specific overlap representations will be examined for both diagnostic and interactive radiation runs in GEOS-5 and comparisons will be made with observed CREs from CERES and CloudSat (2B-FLXHR product). Since the radiative effects of overlap depend on the cloud property distributions of the AGCM, the availability of two different cloud schemes in GEOS-5 will give us the opportunity to assess a wide range of potential cloud overlap consequences on the model's climate.

Oreopoulos, Lazaros↗
NASA NTRS1993

Dust in the Small Magellanic Cloud

Observations of reddened stars in the Small Magellanic Cloud (SMC) indicate that the interstellar grains in that galaxy may show distinct optical properties from those in the Galaxy. In a careful study of three SMC objects, Prevot showed that the UV extinction law in the SMC is almost linear with inverse wavelength and the 2200A feature is generally absent. The first results of a program to determine the wavelength dependence of the interstellar optical polarization in the SMC indicate that highly polarized objects are scarce. Our study has uncovered, however, several objects with optical polarization greater than around 1 percent: AZV 126, AZV 211, AZV 221, AZV 398, and AZV 456. The latter two have already had their UV extinction law determined. Our aim was to obtain International Ultraviolet Explorer (IUE) data and determine the UV extinction law also for AZV 126, AZV 211, and AZV 221. AZV 456, which presents a 'galactic' extinction law, has a 'normal' value for its wavelength of maximum polarization, lmax, while AZV 398, which shows a 'typical' SMC extinction curve, shows a somewhat smaller value for such wavelength. AZV 126, AZV 211, and AZV 221 all present extreme small values of lmax but had not had its extinction curve in the W determined yet. We therefore aimed at ultimately determining the extinction law in the direction of these three objects. Such results, in combination with the optical polarization data, have an important bearing on constraining the composition and size distribution of the interstellar dust in the SMC. In the last report, the images gathered with IUE, their processing, and the extinction curves derived from them were described. Such extinction curves and the theoretical models developed to interpret the SMC extinction and polarization data are discussed. Details are presented in an enclosed preprint. The activities in our ongoing polarimetric program of determining the magnetic field structure of the SMC and the images collected at Cerro Tololo Interamerican Observatory during the period are also briefly described. Other activities are also described.

Magalhaes, A. M.↗
NASA NTRS2013

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗
NASA NTRS2010

Peak Wind Tool for General Forecasting

The expected peak wind speed of the day is an important forecast element in the 45th Weather Squadron's (45 WS) daily 24-Hour and Weekly Planning Forecasts. The forecasts are used for ground and space launch operations at the Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45 WS also issues wind advisories for KSC/CCAFS when they expect wind gusts to meet or exceed 25 kt, 35 kt and 50 kt thresholds at any level from the surface to 300 ft. The 45 WS forecasters have indicated peak wind speeds are challenging to forecast, particularly in the cool season months of October - April. In Phase I of this task, the Applied Meteorology Unit (AMU) developed a tool to help the 45 WS forecast non-convective winds at KSC/CCAFS for the 24-hour period of 0800 to 0800 local time. The tool was delivered as a Microsoft Excel graphical user interface (GUI). The GUI displayed the forecast of peak wind speed, 5-minute average wind speed at the time of the peak wind, timing of the peak wind and probability the peak speed would meet or exceed 25 kt, 35 kt and 50 kt. For the current task (Phase II ), the 45 WS requested additional observations be used for the creation of the forecast equations by expanding the period of record (POR). Additional parameters were evaluated as predictors, including wind speeds between 500 ft and 3000 ft, static stability classification, Bulk Richardson Number, mixing depth, vertical wind shear, temperature inversion strength and depth and wind direction. Using a verification data set, the AMU compared the performance of the Phase I and II prediction methods. Just as in Phase I, the tool was delivered as a Microsoft Excel GUI. The 45 WS requested the tool also be available in the Meteorological Interactive Data Display System (MIDDS). The AMU first expanded the POR by two years by adding tower observations, surface observations and CCAFS (XMR) soundings for the cool season months of March 2007 to April 2009. The POR was expanded again by six years, from October 1996 to April 2002, by interpolating 1000-ft sounding data to 100-ft increments. The Phase II developmental data set included observations for the cool season months of October 1996 to February 2007. The AMU calculated 68 candidate predictors from the XMR soundings, to include 19 stability parameters, 48 wind speed parameters and one wind shear parameter. Each day in the data set was stratified by synoptic weather pattern, low-level wind direction, precipitation and Richardson Number, for a total of 60 stratification methods. Linear regression equations, using the 68 predictors and 60 stratification methods, were created for the tool's three forecast parameters: the highest peak wind speed of the day (PWSD), 5-minute average speed at the same time (A WSD), and timing of the PWSD. For PWSD and A WSD, 30 Phase II methods were selected for evaluation in the verification data set. For timing of the PWSD, 12 Phase\I methods were selected for evaluation. The verification data set contained observations for the cool season months of March 2007 to April 2009. The data set was used to compare the Phase I and II forecast methods to climatology, model forecast winds and wind advisories issued by the 45 WS. The model forecast winds were derived from the 0000 and 1200 UTC runs of the 12-km North American Mesoscale (MesoNAM) model. The forecast methods that performed the best in the verification data set were selected for the Phase II version of the tool. For PWSD and A WSD, linear regression equations based on MesoNAM forecasts performed significantly better than the Phase I and II methods. For timing of the PWSD, none of the methods performed significantly bener than climatology. The AMU then developed the Microsoft Excel and MIDDS GUls. The GUIs display the forecasts for PWSD, AWSD and the probability the PWSD will meet or exceed 25 kt, 35 kt and 50 kt. Since none of the prediction methods for timing of the PWSD performed significantly better thanlimatology, the tool no longer displays this predictand. The Excel and MIDDS GUIs display forecasts for Day-I to Day-3 and Day-I to Day-5, respectively. The Excel GUI uses MesoNAM forecasts as input, while the MIDDS GUI uses input from the MesoNAM and Global Forecast System model. Based on feedback from the 45 WS, the AMU added the daily average wind speed from 30 ft to 60 ft to the tool, which is one of the parameters in the 24-Hour and Weekly Planning Forecasts issued by the 45 WS. In addition, the AMU expanded the MIDDS GUI to include forecasts out to Day-7.

Barrett, Joe H., III↗