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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗
NASA NTRS2012

Bending of Light in Modified Gravity at Large Distances

We discuss the bending of light in a recent model for gravity at large distances containing a Rindler type acceleration proposed by Grumiller. We consider the static, spherically symmetric metric with cosmological constant and Rindler-like term 2ar presented in this model, and we use the procedure by Rindler and Ishak. to obtain the bending angle of light in this metric. Earlier work on light bending in this model by Carloni, Grumiller, and Preis, using the method normally employed for asymptotically flat space-times, led to a conflicting result (caused by the Rindler-like term in the metric) of a bending angle that increases with the distance of closest approach r(sub 0) of the light ray from the centrally concentrated spherically symmetric matter distribution. However, when using the alternative approach for light bending in nonasymptotically flat space-times, we show that the linear Rindler-like term produces a small correction to the general relativistic result that is inversely proportional to r(sub 0). This will in turn affect the bounds on Rindler acceleration obtained earlier from light bending and casts doubts on the nature of the linear term 2ar in the metric

Grumiller↗
NASA NTRS2004

Eruption of a Multiple-Turn Helical Magnetic Flux Tube in a Large Flare: Evidence for External and Internal Reconnection that Fits the Breakout Model of Solar Magnetic Eruptions

We present observations and an interpretation of a unique multiple-turn spiral flux tube eruption from active region 10030 on 2002 July 15. The TRACE C IV observations clearly show a flux tube that is helical and erupting from within a sheared magnetic field. These observations are interpreted in the context of the breakout model for magnetic field explosions. The initiation of the helix eruption. as determined by a linear backward extrapolation, starts 25 s after the peak of the flare's strongest impulsive spike of microwave gyrosynchrotron radiation early in the flare s explosive phase, implying that the sheared core field is not the site of the initial reconnection. Within the quadrupolar configuration of the active region, the external and internal reconnection sites are identified in each of two consecutive eruptive flares that produce a double coronal mass ejection (CME). The first external breakout reconnection apparently releases an underlying sheared core field and allows it to erupt, leading to internal reconnection in the wake of the erupting helix. This internal reconnection releases the helix and heats the two-ribbon flare. These events lead to the first CME and are followed by a second breakout that initiates a second and larger halo CME. The strong magnetic shear in the region is compatible with the observed rapid proper motion and evolution of the active region. The multiple-turn helix originates from above a sheared-field magnetic inversion line within a filament channel. and starts to erupt only after fast breakout reconnection has started. These observations are counter to the standard flare model and support the breakout model for eruptive flare initiation.

Gary, G. Allen↗
NASA NTRS1977

Unsteady aerodynamic modeling and active aeroelastic control

Unsteady aerodynamic modeling techniques are developed and applied to the study of active control of elastic vehicles. The problem of active control of a supercritical flutter mode poses a definite design goal stability, and is treated in detail. The transfer functions relating the arbitrary airfoil motions to the airloads are derived from the Laplace transforms of the linearized airload expressions for incompressible two dimensional flow. The transfer function relating the motions to the circulatory part of these loads is recognized as the Theodorsen function extended to complex values of reduced frequency, and is termed the generalized Theodorsen function. Inversion of the Laplace transforms yields exact transient airloads and airfoil motions. Exact root loci of aeroelastic modes are calculated, providing quantitative information regarding subcritical and supercritical flutter conditions.

Edwards, J. W.↗
NASA NTRS1988

Mass and energy flow near sunspots

Sunspots block the flow of energy to the solar surface. The blocked energy heats the volume beneath the spot, producing a pressure excess which drives an outflow of mass. Linear numerical models of the mass and energy flow around spots were constructed to estimate the predictions of this physical picture against the observed properties of sunspot bright rings and moat flows. The width of the bright ring and moat are predicted to be proportional to the depth of the spot penumbra, in conflict with the observed proportionality of the moat width to the spot diameter. Postulating that spot depths are proportional to spot diameters would bury the moat flow too deeply to be observed, because the radial velocity at the surface is found to be inversely proportional to the depth of the spot penumbra. The radial velocity at the surface is of order a few hundred meters per second after 1 day, in agreement with the observed excess of moat velocities over supergranule velocities.

Nye, Alan↗
NASA NTRS2012

Application of Artificial Neural Networks to the Development of Improved Multi-Sensor Retrievals of Near-Surface Air Temperature and Humidity Over Ocean

Improved estimates of near-surface air temperature and air humidity are critical to the development of more accurate turbulent surface heat fluxes over the ocean. Recent progress in retrieving these parameters has been made through the application of artificial neural networks (ANN) and the use of multi-sensor passive microwave observations. Details are provided on the development of an improved retrieval algorithm that applies the nonlinear statistical ANN methodology to a set of observations from the Advanced Microwave Scanning Radiometer (AMSR-E) and the Advanced Microwave Sounding Unit (AMSU-A) that are currently available from the NASA AQUA satellite platform. Statistical inversion techniques require an adequate training dataset to properly capture embedded physical relationships. The development of multiple training datasets containing only in-situ observations, only synthetic observations produced using the Community Radiative Transfer Model (CRTM), or a mixture of each is discussed. An intercomparison of results using each training dataset is provided to highlight the relative advantages and disadvantages of each methodology. Particular emphasis will be placed on the development of retrievals in cloudy versus clear-sky conditions. Near-surface air temperature and humidity retrievals using the multi-sensor ANN algorithms are compared to previous linear and non-linear retrieval schemes.

Roberts, J. Brent↗
NASA NTRS1978

Evidence for a non-random magnetization of the moon

An approximate inversion, using a specific model, of Apollo subsatellite vector magnetometer data has provided estimates for the locations, scale sizes, and magnetization characteristics of strongly magnetized lunar crustal materials across a portion of the lunar far side. The nearly model-independent local bulk directions of magnetization, while otherwise random, exhibit a tendency to lie in or near planes perpendicular to the present lunar spin axis. The apparent existence of such a preferred plane of magnetization must be regarded as unfavorable for proposed local, presumably random sources of large-scale, coherent magnetization and implies that any assumed internal magnetizing dipole or external linear magnetizing field must also have been oriented mainly perpendicular to the present spin axis. The otherwise random orientations of these vectors when projected onto the equatorial plane suggests a simple explanation for the negligible lunar global dipole moment.

Hood, L.↗
NASA NTRS2005

Adjoint Sensitivity Analysis of Orbital Mechanics: Application to Computations of Observables' Partials with Respect to Harmonics of the Planetary Gravity Fields

An approach is presented to the inversion of gravity fields based on evaluation of partials of observables with respect to gravity harmonics using the solution of adjoint problem of orbital dynamics of the spacecraft. Corresponding adjoint operator is derived directly from the linear operator of the linearized forward problem of orbital dynamics. The resulting adjoint problem is similar to the forward problem and can be solved by the same methods. For given highest degree N of gravity harmonics desired, this method involves integration of N adjoint solutions as compared to integration of N2 partials of the forward solution with respect to gravity harmonics in the conventional approach. Thus, for higher resolution gravity models, this approach becomes increasingly more effective in terms of computer resources as compared to the approach based on the solution of the forward problem of orbital dynamics.

gravity fields↗
NASA NTRS1972

The general linear inverse problem - Implication of surface waves and free oscillations for earth structure.

The discrete general linear inverse problem reduces to a set of m equations in n unknowns. There is generally no unique solution, but we can find k linear combinations of parameters for which restraints are determined. The parameter combinations are given by the eigenvectors of the coefficient matrix. The number k is determined by the ratio of the standard deviations of the observations to the allowable standard deviations in the resulting solution. Various linear combinations of the eigenvectors can be used to determine parameter resolution and information distribution among the observations. Thus we can determine where information comes from among the observations and exactly how it constraints the set of possible models. The application of such analyses to surface-wave and free-oscillation observations indicates that (1) phase, group, and amplitude observations for any particular mode provide basically the same type of information about the model; (2) observations of overtones can enhance the resolution considerably; and (3) the degree of resolution has generally been overestimated for many model determinations made from surface waves.

Wiggins, R. A.↗
NASA NTRS1996

Possible effects of protracted exposure on the additivity of risks from space radiations

Conventional radiation risk assessments are presently based on the additivity assumption. This assumption states that risks from individual components of a complex radiation field involving many different types of radiation can be added to yield the total risk of the complex radiation field. If the assumption is not correct, the summations and integrations performed to obtain the presently quoted risk estimates are not appropriate. This problem is particularly important in the area of space radiation risk evaluation because of the many different types of high- and low-LET radiation present in the galactic cosmic ray environment. For both low- and high-LET radiations at low enough dose rates, the present convention is that the addivity assumption holds. Mathematically, the total risk, Rtot is assumed to be Rtot = summation (i) Ri where the summation runs over the different types of radiation present. If the total dose (or fluence) from each component is such that the interaction between biological lesions caused by separate single track traversals is negligible within a given cell, it is presently considered to be reasonable to accept the additivity assumption. However, when the exposure is protracted over many cell doubling times (as will be the case for extended missions to the moon or Mars), the possibility exists that radiation effects that depend on multiple cellular events over a long time period, such as is probably the case in radiation-induced carcinogenesis, may not be additive in the above sense and the exposure interval may have to be included in the evaluation procedure. It is shown, however, that "inverse" dose-rate effects are not expected from intermediate LET radiations arising from the galactic cosmic ray environment due to the "sensitive-window-in-the-cell-cycle" hypothesis.

NASA Discipline Number 45-10↗
NASA NTRS2010

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station

The expected peak wind speed for the day is an important element in the daily morning forecast for ground and space launch operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45th Weather Squadron (45 WS) must issue forecast advisories for KSC/CCAFS when they expect peak gusts for >= 25, >= 35, and >= 50 kt thresholds at any level from the surface to 300 ft. In Phase I of this task, the 45 WS tasked the Applied Meteorology Unit (AMU) to develop a cool-season (October - April) tool to help forecast the non-convective peak wind from the surface to 300 ft at KSC/CCAFS. During the warm season, these wind speeds are rarely exceeded except during convective winds or under the influence of tropical cyclones, for which other techniques are already in use. The tool used single and multiple linear regression equations to predict the peak wind from the morning sounding. The forecaster manually entered several observed sounding parameters into a Microsoft Excel graphical user interface (GUI), and then the tool displayed the forecast peak wind speed, average wind speed at the time of the peak wind, the timing of the peak wind and the probability the peak wind will meet or exceed 35, 50 and 60 kt. The 45 WS customers later dropped the requirement for >= 60 kt wind warnings. During Phase II of this task, the AMU expanded the period of record (POR) by six years to increase the number of observations used to create the forecast equations. A large number of possible predictors were evaluated from archived soundings, including inversion depth and strength, low-level wind shear, mixing height, temperature lapse rate and winds from the surface to 3000 ft. Each day in the POR was stratified in a number of ways, such as by low-level wind direction, synoptic weather pattern, precipitation and Bulk Richardson number. The most accurate Phase II equations were then selected for an independent verification. The Phase I and II forecast methods were compared using an independent verification data set. The two methods were compared to climatology, wind warnings and advisories issued by the 45 WS, and North American Mesoscale (NAM) model (MesoNAM) forecast winds. The performance of the Phase I and II methods were similar with respect to mean absolute error. Since the Phase I data were not stratified by precipitation, this method's peak wind forecasts had a large negative bias on days with precipitation and a small positive bias on days with no precipitation. Overall, the climatology methods performed the worst while the MesoNAM performed the best. Since the MesoNAM winds were the most accurate in the comparison, the final version of the tool was based on the MesoNAM winds. The probability the peak wind will meet or exceed the warning thresholds were based on the one standard deviation error bars from the linear regression. For example, the linear regression might forecast the most likely peak speed to be 35 kt and the error bars used to calculate that the probability of >= 25 kt = 76%, the probability of >= 35 kt = 50%, and the probability of >= 50 kt = 19%. The authors have not seen this application of linear regression error bars in any other meteorological applications. Although probability forecast tools should usually be developed with logistic regression, this technique could be easily generalized to any linear regression forecast tool to estimate the probability of exceeding any desired threshold . This could be useful for previously developed linear regression forecast tools or new forecast applications where statistical analysis software to perform logistic regression is not available. The tool was delivered in two formats - a Microsoft Excel GUI and a Tool Command Language/Tool Kit (Tcl/Tk) GUI in the Meteorological Interactive Data Display System (MIDDS). The Microsoft Excel GUI reads a MesoNAM text file containing hourly forecasts from 0 to 84 hours, from one model run (00 or 12 UTC). The GUI then displays e peak wind speed, average wind speed, and the probability the peak wind will meet or exceed the 25-, 35- and 50-kt thresholds. The user can display the Day-1 through Day-3 peak wind forecasts, and separate forecasts are made for precipitation and non-precipitation days. The MIDDS GUI uses data from the NAM and Global Forecast System (GFS), instead of the MesoNAM. It can display Day-1 and Day-2 forecasts using NAM data, and Day-1 through Day-5 forecasts using GFS data. The timing of the peak wind is not displayed, since the independent verification showed that none of the forecast methods performed significantly better than climatology. The forecaster should use the climatological timing of the peak wind (2248 UTC) as a first guess and then adjust it based on the movement of weather features.

Barrett, Joe H., III↗
NASA NTRS1991

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗
NASA NTRS2016

GOCI Yonsei Aerosol Retrieval (YAER) Algorithm and Validation During the DRAGON-NE Asia 2012 Campaign

The Geostationary Ocean Color Imager (GOCI) onboard the Communication, Ocean, and Meteorological Satellite (COMS) is the first multi-channel ocean color imager in geostationary orbit. Hourly GOCI top-of-atmosphere radiance has been available for the retrieval of aerosol optical properties over East Asia since March 2011. This study presents improvements made to the GOCI Yonsei Aerosol Retrieval (YAER) algorithm together with validation results during the Distributed Regional Aerosol Gridded Observation Networks - Northeast Asia 2012 campaign (DRAGONNE Asia 2012 campaign). The evaluation during the spring season over East Asia is important because of high aerosol concentrations and diverse types of Asian dust and haze. Optical properties of aerosol are retrieved from the GOCI YAER algorithm including aerosol optical depth (AOD) at 550 nm, fine-mode fraction (FMF) at 550 nm, single-scattering albedo (SSA) at 440 nm, Angstrom exponent (AE) between 440 and 860 nm, and aerosol type. The aerosol models are created based on a global analysis of the Aerosol Robotic Networks (AERONET) inversion data, and covers a broad range of size distribution and absorptivity, including nonspherical dust properties. The Cox-Munk ocean bidirectional reflectance distribution function (BRDF) model is used over ocean, and an improved minimum reflectance technique is used over land. Because turbid water is persistent over the Yellow Sea, the land algorithm is used for such cases. The aerosol products are evaluated against AERONET observations and MODIS Collection 6 aerosol products retrieved from Dark Target (DT) and Deep Blue (DB) algorithms during the DRAGON-NE Asia 2012 campaign conducted from March to May 2012. Comparison of AOD from GOCI and AERONET resulted in a Pearson correlation coefficient of 0.881 and a linear regression equation with GOCI AOD = 1.083 x AERONET AOD - 0.042. The correlation between GOCI and MODIS AODs is higher over ocean than land. GOCI AOD shows better agreement with MODIS DB than MODIS DT. The other GOCI YAER products (AE, FMF, and SSA) show lower correlation with AERONET than AOD, but still show some skills for qualitative use.

Choi, Myungje↗
NASA NTRS1989

Characterization of robotics parallel algorithms and mapping onto a reconfigurable SIMD machine

The kinematics, dynamics, Jacobian, and their corresponding inverse computations are six essential problems in the control of robot manipulators. Efficient parallel algorithms for these computations are discussed and analyzed. Their characteristics are identified and a scheme on the mapping of these algorithms to a reconfigurable parallel architecture is presented. Based on the characteristics including type of parallelism, degree of parallelism, uniformity of the operations, fundamental operations, data dependencies, and communication requirement, it is shown that most of the algorithms for robotic computations possess highly regular properties and some common structures, especially the linear recursive structure. Moreover, they are well-suited to be implemented on a single-instruction-stream multiple-data-stream (SIMD) computer with reconfigurable interconnection network. The model of a reconfigurable dual network SIMD machine with internal direct feedback is introduced. A systematic procedure internal direct feedback is introduced. A systematic procedure to map these computations to the proposed machine is presented. A new scheduling problem for SIMD machines is investigated and a heuristic algorithm, called neighborhood scheduling, that reorders the processing sequence of subtasks to reduce the communication time is described. Mapping results of a benchmark algorithm are illustrated and discussed.

Lee, C. S. G.↗
NASA NTRS1993

Are the 1986-1988 changes in solar free-oscillation frequency splitting significant?

The solar normal-mode splitting coefficients deduced from Big Bear Solar Observatory (BBSO) data differ between 1986 and 1988; inversions for equatorial rotation are slower at depth and faster near the surface in 1988 than in 1986. The significance of the change has been disputed. The data sets overlap for five splitting coefficients (a(j))super 5 sub j = 1 associated with 710 multiplets. On the assumption that rotation rate varies smoothly with radius, both data sets are satisfied by the same rotation model at all colatitudes except near 30-40 deg and near 70 deg (and at their southern hemisphere reflections 140-150 deg and 110 deg). The evidence for equatorial change is weak. Nonparametric tests show a significant offset in the magnitudes of a(1), a(2), and a(4), and of linear combinations sensitive to rotation at colatitudes of 60-80 deg (and 120 deg). Nonparametric tests show significant radial trends in the changes to a(2), a(4), and (less significantly) a(5). There is strong anticorrelation between a(2) and a(4), a(1) and a(3), and a(3) and a(5), suggesting that the estimates are not independent. Individual coefficients a(j) show more evidence for change than do 'physical' linear combinations, adding weight to this hypothesis. Some of the changes in splitting might be related to solar activity, which changed most near colatitude 70 deg from 1986 to 1988.

Gough, Douglas↗
NASA NTRS1995

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗
NASA NTRS1990

Interpretation of ground-based radiometric observations in terms of a gravity wave model

An analysis is presented of 2-hour and 4-hour segments of data taken at Denver, Colorado, on February 3, 1984, by a ground-based radiometer designed and operated by the Wave Propagation Laboratory of the National Oceanic and Atmospheric Administration (NOAA). The zenith-viewing instrument has two moisture-sensing and four temperature-sensing channels. It is demonstrated that a peak at a period of 10 min, present in the spectra of the measured brightness temperature and of the derived geopotential heights, thicknesses, and vertically integrated water vapor content, is due to an internal gravity wave generated by wind shear in the jet aloft. This analysis shows that the radiometer has the sensitivity to detect such disturbances and that the mathematical inversion technique used to retrieve the geopotential field and other integrated quantities retains the derived information as well. Finally, a linear expression is derived which relates the brightness temperature to the atmospheric temperature, density, humidity, and cloud liquid perturbation fields.

Canavero, F. G.↗