Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inversion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

Interannual Variation in Phytoplankton Concentration and Community in the Pacific Ocean

Climate events such as El Nino have been shown to have an effect on the biology of our ocean. Because of the lack of data, we still have very little knowledge about the spatial and temporal effect these climate events may have on biological marine systems. In this study, we used the NASA Ocean Biogeochemical Model (NOBM) to assess the interannual variability in phytoplankton community in the Pacific Ocean between 1998 and 2005. In the North Central and Equatorial Pacific Ocean, changes in the Multivariate El Nino Index were associated with changes in phytoplankton composition. The model identified an increase in diatoms of approx.33 % in the equatorial Pacific in 1999 during a La Nina event. This increase in diatoms coincided with a decrease of approx.11 % in cyanobacteria concentration. The inverse relationship between cyanobacteria and diatoms concentration was significant (p<0.05) throughout the period of study. The use of a numerical model allows us to assess the impact climate variability has on key phytoplankton groups known to lead to contrasting food chain at a spatial and temporal resolution unachievable when relying solely on in-situ observations.

Rousseaux, C. S.↗

Far-infrared spectral observations of M42 and M17

Results are reported for coarse-resolution far-IR spectral observations of the H II regions M17 and M41 in the wavelength range from 42 to 115 microns. The observations were made at an altitude of 13.7 km with the 30-cm telescope of the NASA Lear jet. The spectrum of M42 is shown to be quite smooth and to fit a diluted 100-K blackbody very well. It is found that the emission from M17 falls steeply from 40 to 63 microns, then rises again to 80 microns before falling off once more toward longer wavelengths. The best fit of this spectrum for an inverse-square-wavelength emissivity law is found to require a two-temperature model with temperatures below 30 K for the cold dust component. It is concluded that the spectrum of M17 seems to require the presence of a material with a better far-IR emissivity than those for which detailed laboratory measurements have been made.

Ward, D. B.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

An analysis of the shortening and lengthening of the phase of satellite-to-satellite (SST) data that passed within 40-700 km above the earth surface during its ATS-6 to GEOS-3 to ATS-6 path resulted in refractivity vs height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbations from the moon, sun, and a 15th-order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occulation period and a model ionosphere was estimated by a least-squares adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the two methods were caused by orbital errors, which propagated into the solution.

Liu, A. S.↗

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

When the radio link between two satellites (GEOS-3/ATS-6) is intercepted by the earth's ionosphere and neutral atmosphere, a change in the Doppler frequency results. Travel through the atmosphere causes the Doppler phase to be advanced in the ionosphere's portion and retarded in the neutral portion of the atmosphere. Analysis of the shortening and lengthening of the phase of the Satellite-to-Satellite Tracking (SST) data that passed within 40-700 km above the earth's surface during its ATS-6 to GEOS-3 to ATS-6 path, caused by the atmosphere, results in refractivity versus height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbation from the Moon, Sun and a 15th order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occultation period and a model ionosphere was estimated by a least-square adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the 2 methods were caused by orbital errors, which propagated into the solution. The SST data yielded refractive index profiles in a novel economical manner because no additional or special on-board equipment were required.

Liu, A. S.↗

Near-infrared observations of IRAS sources in and near dense cores

We report 0.4 to 20 micron photometry of 27 Infrared Astronomy Satellite (IRAS) point sources associated with dense cores in nearby dark clouds. It is found that these objects have a bimodal distribution of spectral slopes. The objects in the group with steep spectral slopes are typically within one half power radius of the core peak and are not visible on the POSS. The remainder of the objects are typically further from their cores, are optically visible and have shallow spectral slopes (s is approximately 0.6). Most of the sources in this group which have been previously identified are T Tauri stars. Both the groups of objects have essentially the same luminosity function with median luminosity of 1 to 2 L(solar), which is similiar to the luminosity function for T Tauri stars. The near-infrared (J,H,K,L) colors of the objects have been used to estimate the visual extinction to the stars. This indicates that the typical steep spectrum source has extinction A(upsilon) of approximately 30 magnitudes, which is larger by a factor of 3 to 4 than a uniform core can provide, thus the density must rise steeply in the vicinity of the star. For those objects where optical estimates of the extinction are available the optical estimates are typically smaller than our IR estimates by a factor of 1.6. For multiple scattering by grains that scatter primarily in the forward direction, this implies a grain albedo of 0.4 at V, consistent with theoretical predictions for bare grains. Using this derived extinction, the FIR emission from a spherical shell with an inverse power law density which contributes the correct amount of extinction has been modelled. It is shown that to be consistent with both the observations of the density of the core at a scale of 1.5 x 10(17) cm and the extinction derived from the near-IR colors there must be a circumstellar hole of radius 10-100 A.U. for all density laws with exponents greater than 1.2. It is also shown that a component in addition to the star and circumstellar shell is necessary to match the observed spectrum of steep spectrum sources.

Fuller, G. A.↗

PNS calculations for 3-D hypersonic corner flow with two turbulence models

A three-dimensional parabolized Navier-Stokes code has been used as a testbed to investigate two turbulence models, the McDonald Camarata and Bushnell Beckwith model, in the hypersonic regime. The Bushnell Beckwith form factor correction to the McDonald Camarata mixing length model has been extended to three-dimensional flow by use of an inverse averaging of the resultant length scale contributions from each wall. Two-dimensional calculations are compared with experiment for Mach 18 helium flow over a 4-deg wedge. Corner flow calculations have been performed at Mach 11.8 for a Reynolds number of .67 x 10 to the 6th, based on the duct half-width, and a freestream stagnation temperature of 1750-deg Rankine.

Smith, Gregory E.↗

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗

An adaptive inverse kinematics algorithm for robot manipulators

An adaptive algorithm for solving the inverse kinematics problem for robot manipulators is presented. The algorithm is derived using model reference adaptive control (MRAC) theory and is computationally efficient for online applications. The scheme requires no a priori knowledge of the kinematics of the robot if Cartesian end-effector sensing is available, and it requires knowledge of only the forward kinematics if joint position sensing is used. Computer simulation results are given for the redundant seven-DOF robotics research arm, demonstrating that the proposed algorithm yields accurate joint angle trajectories for a given end-effector position/orientation trajectory.

Colbaugh, R.↗

Calculation of unsteady transonic flows with mild separation by viscous-inviscid interaction

This paper presents a method for calculating viscous effects in two- and three-dimensional unsteady transonic flow fields. An integral boundary-layer method for turbulent viscous flow is coupled with the transonic small-disturbance potential equation in a quasi-steady manner. The viscous effects are modeled with Green's lag-entrainment equations for attached flow and an inverse boundary-layer method for flows that involve mild separation. The boundary-layer method is used stripwise to approximate three-dimensional effects. Applications are given for two-dimensional airfoils, aileron buzz, and a wing planform. Comparisons with inviscid calculations, other viscous calculation methods, and experimental data are presented. The results demonstrate that the present technique can economically and accurately calculate unsteady transonic flow fields that have viscous-inviscid interactions with mild flow separation.

Howlett, James T.↗

Singular Isothermal Disks: Nonaxiymmetric Bifurcations and Equilibria - Paper 2

We review the difficulties of the classical fission and fragmentation hypotheses for the formation of binary and multiple stars. A crucial missing ingredient in previous theoretical studies is the inclusion of dynamically important levels of magnetic fields. As a minimal model for a candidate presursor to the formation of binary and multiple stars, we therefore formulate and solve the problem of the equilibria of isopedically magnetized, singular isothermal disks, without the assumption of axial symmetry. Considerable analytical progress can be made if we restrict our attention to models that are scale-free, i.e., that have surface densities that vary inversely with distance omega from the rotation axis of the system. In agreement with earlier analysis by Syer and Tremaine, we find that lopsided (M = 1) configurations exist at any dimensionless rotation rate, including zero. Multiple-lobed (M = 2, 3, 4, ...) configurations bifurcate from an underlying axisymmetric sequence at progressively higher dimensionless rates of rotation, but such nonaxisymmetric sequences always terminate in shockwaves before they have a chance to fission into M = 2, 3, 4, ... separate bodies. On the basis of our experience in this paper, we advance the hypothesis that binary and multiple star-formation from smooth (i.e., not highly turbulent) starting states that are supercritical but in unstable mechanical balance requires the rapid (i.e., dynamical) loss of magnetic flux at some stage of the ensuing gravitational collapse.

Galli, Danielle↗

Tidal Implications for the GRACE Satellite Mission

A fundamental part of the processing of the GRACE satellite data will be correction for the gravitational effect of tides. This includes solid earth, ocean, and atmospheric tides. The most problematic are ocean tides in shallow seas and in high latitudes (above the 66$^{\circ}$ inclination of Topex). Somewhat smaller, but still important, errors are associated with the S$_1$ air tide and the ocean's dynamic response to it. In this paper we explore several facets of the tidal error budget. The ocean tide error spectrum is estimated (a) by inverse methods and (b) by simple differencing of several of the best tidal models. The tidal error spectrum is flatter than the tidal signal spectrum, and it exceeds 10\% of the signal at degrees 15 and above. The tide errors also exceed the anticipated sensitivity of GRACE for all degrees below 40, and possibly below 50. The assumption of constant seawater density (of about 1035 kg\,M$^{-3}$) is a minor error source. Owing to GRACE's high inclination, solar tides are aliased to long periods, which is potentially problematic for geophysical interpretation of the gravity signal. Simulations of the effect of high-latitude solar-tide errors on gravity estimates around Antarctica will be discussed.

Ray, Richard D.↗

Flame Chemiluminescence Rate Constants for Quantitative Microgravity Combustion Diagnostics

Absolute excited state concentrations of OH(A), CH(A), and C2(d) were determined in three low pressure premixed methane-air flames. Two dimensional images of chemiluminescence from these states were recorded by a filtered CCD camera, processed by Abel inversion, and calibrated against Rayleigh scattering, Using a previously validated 1-D flame model with known chemistry and excited state quenching rate constants, rate constants are extracted for the reactions CH + O2 (goes to) OH(A) + CO and C2H + O (goes to) CH(A) + CO at flame temperatures. Variations of flame emission intensities with stoichiometry agree well with model predictions.

Luque, Jorge↗

Device Concepts Based on Spin-dependent Transmission in Semiconductor Heterostructures

We examine zero-magnetic-field spin-dependent transmission in nonmagnetic semiconductor heterostructures with structural inversion asymmetry (SIA) and bulk inversion asymmetry (BIA), and report spin devices concepts that exploit their properties. Our modeling results show that several design strategies could be used to achieve high spin filtering efficiencies. The current spin polarization of these devices is electrically controllable, and potentially amenable to highspeed spin modulation, and could be integrated in optoelectronic devices for added functionality.

resonant tunneling↗

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↗

Shape and Stress Sensing of Multilayered Composite and Sandwich Structures Using an Inverse Finite Element Method

The marked increase in the use of composite and sandwich material systems in aerospace, civil, and marine structures leads to the need for integrated Structural Health Management systems. A key capability to enable such systems is the real-time reconstruction of structural deformations, stresses, and failure criteria that are inferred from in-situ, discrete-location strain measurements. This technology is commonly referred to as shape- and stress-sensing. Presented herein is a computationally efficient shape- and stress-sensing methodology that is ideally suited for applications to laminated composite and sandwich structures. The new approach employs the inverse Finite Element Method (iFEM) as a general framework and the Refined Zigzag Theory (RZT) as the underlying plate theory. A three-node inverse plate finite element is formulated. The element formulation enables robust and efficient modeling of plate structures instrumented with strain sensors that have arbitrary positions. The methodology leads to a set of linear algebraic equations that are solved efficiently for the unknown nodal displacements. These displacements are then used at the finite element level to compute full-field strains, stresses, and failure criteria that are in turn used to assess structural integrity. Numerical results for multilayered, highly heterogeneous laminates demonstrate the unique capability of this new formulation for shape- and stress-sensing.

Cerracchio, Priscilla↗

Saturn Ring Seismology: Interpreting the Seismogram

Marley (1990) and Marley and Porco (1993) proposed that f-mode oscillations of Saturn could excite resonant density and bending waves in the inner C-ring. They hypothesized that certain wave features discovered by Rosen et al. (1991) that were not associated with known satellite resonances could be the result of such resonant interactions with the planetary oscillation modes. They also predicted that if this was the case the waves would be found to be density (and not bending) waves by Cassini and predicted the azimuthal wave number of the C-ring waves m. Employing Cassini VIMS stellar occultation data Hedman and Nicholson (2013) have now confirmed the predictions and demonstrated that at least some of the C-ring features identified by Rosen et al. are indeed likely caused by resonant oscillation modes of Saturn. Given this context we have taken a fresh look at the Saturn ring seismology. First we propose that an apparent bending wave denoted 'j' by Rosen may be a second order outer vertical resonance with the l=3, m=2 f-mode of Saturn and discuss the locations of other plausible second order resonances in the rings. Since only a handful of ring resonances have been identified, measuring even one or two additional planetary mode frequencies would substantially assist the process of inverting mode frequencies to constrain Saturn interior's structure. Using the available mode frequencies, modern inversion technique employed in stellar seismology, and a recent set of Saturn interior models we provide an initial estimation of what available mode frequencies are telling us about the interior structure of the planet. Since the f-modes are confined relatively closely to the planetary surface, most of the observed modes probe only the outermost layers of the planet that are already comparatively well understood. However the l = 2 mode does probe relatively deeply into the planet and we will discuss the potential the measurement of this mode frequency has for placing new constraints on the interior structure.

Seismogram↗

Uncertainty Quantification of GEOS-5 L-band Radiative Transfer Model Parameters Using Bayesian Inference and SMOS Observations

Uncertainties in L-band (1.4 GHz) radiative transfer modeling (RTM) affect the simulation of brightness temperatures (Tb) over land and the inversion of satellite-observed Tb into soil moisture retrievals. In particular, accurate estimates of the microwave soil roughness, vegetation opacity and scattering albedo for large-scale applications are difficult to obtain from field studies and often lack an uncertainty estimate. Here, a Markov Chain Monte Carlo (MCMC) simulation method is used to determine satellite-scale estimates of RTM parameters and their posterior uncertainty by minimizing the misfit between long-term averages and standard deviations of simulated and observed Tb at a range of incidence angles, at horizontal and vertical polarization, and for morning and evening overpasses. Tb simulations are generated with the Goddard Earth Observing System (GEOS-5) and confronted with Tb observations from the Soil Moisture Ocean Salinity (SMOS) mission. The MCMC algorithm suggests that the relative uncertainty of the RTM parameter estimates is typically less than 25 of the maximum a posteriori density (MAP) parameter value. Furthermore, the actual root-mean-square-differences in long-term Tb averages and standard deviations are found consistent with the respective estimated total simulation and observation error standard deviations of m3.1K and s2.4K. It is also shown that the MAP parameter values estimated through MCMC simulation are in close agreement with those obtained with Particle Swarm Optimization (PSO).

MCMC↗