SEARCH · Engineering Papers

Results for “generalized linear model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Gain optimization with nonlinear controls

An algorithm has been developed for the analysis and design of controls for nonlinear systems. The technical approach is to use statistical linearization to model the nonlinear dynamics of a system. A covariance analysis is performed to determine the behavior of the dynamical system and a quadratic cost function. Expressions for the cost function and its derivatives are determined so that numerical optimization techniques can be applied to determine optimal feedback laws. The primary application for this report is centered about the design of controls for nominally linear systems but where the controls are saturated or limited by fixed constraints. The analysis is general however and numerical computation requires only that the specific nonlinearity be considered in the analysis.

Slater, G. L.↗

NASA NTRS1983

On the controllability and control law design for an orbiting large flexible antenna system

The controllability and control law synthesis based on a finite element model of the Hoop/Column orbiting antenna system is considered. Graph theoretic techniques are employed to analyze the controllability of the system for posible proposed actuator arrangements which include torquers and point actuators along the mast and a single actuator placed on the hoop assembly. Once controllability is established for a given combination of number and location of actuators, the synthesis of control laws is based on an application of the linear regulator theory and also pole placement techniques. In general, surface torsion and feed mast torsion are among modes having the longest time constants. System transient performance is noticeably degraded when the hoop actuator is not included.

Bainum, P. M.↗

NASA NTRS1984

Finite Element Prediction of Acoustic Scattering and Radiation from Submerged Elastic Structures

A finite element formulation is derived for the scattering and radiation of acoustic waves from submerged elastic structures. The formulation uses as fundamental unknowns the displacement in the structure and a velocity potential in the field. Symmetric coefficient matrices result. The outer boundary of the fluid region is terminated with an approximate local wave-absorbing boundary condition which assumes that outgoing waves are locally planar. The finite element model is capable of predicting only the near-field acoustic pressures. Far-field sound pressure levels may be determined by integrating the surface pressures and velocities over the wet boundary of the structure using the Helmholtz integral. Comparison of finite element results with analytic results show excellent agreement. The coupled fluid-structure problem may be solved with general purpose finite element codes by using an analogy between the equations of elasticity and the wave equation of linear acoustics.

Everstine, G. C.↗

NASA NTRS1985

The causation and sensitivity of the northern winter planetary waves

High-resolution calculations of the wintertime stationary atmospheric response to planetary scale topographic and thermal forcings have been made using a numerical model that solves the spherical primitive equations linearized about observed zonal wind and temperature fields. The model equations are presented and their numerical implementation is discussed. The model inputs and parameters are presented and the model's response to relativistic forcing is analyzed. In middle latitudes the response to topographic forcing strongly dominates the response to thermal forcing. In the midlatitude troposphere, the topographic response is insensitive to changes in the zonal wind. In the stratosphere, both the variance and mean are dominated by the topographic response. In general, the sensitivity of the planetary waves to changes in the zonal wind is found to be much smaller than other recent calculations have indicated.

Jacqmin, D.↗

NASA NTRS1984

The influence of dynamic inflow and torsional flexibility on rotor damping in forward flight from symbolically generated equations

The combined effects of blade torsion and dynamic inflow on the aeroelastic stability of an elastic rotor blade in forward flight are studied. The Helicopter Equations for Stability and Loads (HESL) program is extended to derive the governing equations of motion for the blade, and a Lagrangian formulation is used to obtain the equations in generalized coordinates. The program generates the steady-state and linearized perturbation equations in symbolic form and then codes them into FORTRAN subroutines. The coefficients for each equation and for each mode are identified through a numerical program; the latter can also be used to obtain the harmonic balance equations. The governing multiblade equations are derived explicitly using HESL. These equations can accommodate any number of elastic blade modes. Stability results are presented for several hingeless rotor blade structural models, and the influence of dynamic inflow in forward flight with an elastic hingeless rotor is investigated.

Reddy, T. S. R.↗

NASA NTRS1990

Simplified calculation of wave-coupling between lines in high-speed integrated circuits

The construction of time-domain solutions by the superposition of modes is explained for the general case, which may not have any symmetry in the lines or terminations. This simplified method requires that the logic-gate terminal impedances be approximated by linear effective resistances that can be adequately determined by SPICE analysis of the actual logic gates. Several examples of results of crosstalk calculations made by the simplified method are presented and compared with corresponding results obtained using the SPICE transmission-line methods of Tripathi and Rettig (1985) along with full SPICE models of the logic gates. The simplified approach is shown to give satisfactory agreement for most purposes.

Matthaei, George L.↗

NASA NTRS1991

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulence models. The difficulty lies in the fact that the reaction rate is in general an exponential function of the temperature, and the higher order correlations in the conventional moment closure models of the chemical source term can not be neglected, making the applications of such models impractical. The probability density function (pdf) method offers an attractive alternative: in a pdf model, the chemical source terms are closed and do not require additional models. A grid dependent Monte Carlo scheme was studied, since it is a logical alternative, wherein the number of computer operations increases only linearly with the increase of number of independent variables, as compared to the exponential increase in a conventional finite difference scheme. A new algorithm was devised that satisfies a restriction in the case of pure diffusion or uniform flow problems. Although for nonuniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗

NASA NTRS2004

Use of Radon for Evaluation of Atmospheric Transport Models: Sensitivity to Emissions

This paper presents comparative analyses of atmospheric radon (Rn) distributions simulated using different emission scenarios and the observations. Results indicate that the model generally reproduces observed distributions of Rn but there are some biases in the model related to differences in large-scale and convective transport. Simulations presented here use an off-line three-dimensional chemical transport model driven by assimilated winds and two scenarios of Rn fluxes (atom/cm s) from ice-free land surfaces: (A) globally uniform flux of 1.0, and (B) uniform flux of 1.0 between 60 deg. S and 30 deg. N followed by a sharp linear decrease to 0.2 at 70 deg. N. We considered an additional scenario (C) where Rn emissions for case A were uniformly reduced by 28%. Results show that case A overpredicts observed Rn distributions in both hemispheres. Simulated northern hemispheric (NH) Rn distributions from cases B and C compare better with the observations, but are not discernible from each other. In the southern hemisphere, surface Rn distributions from case C compare better with the observations. We performed a synoptic scale source-receptor analysis for surface Rn to locate regions with ratios B/A and B/C less than 0.5. Considering an uncertainty in regional Rn emissions of a factor of two, our analysis indicates that additional measurements of surface Rn particularly during April-October and north of 50 deg. N over the Pacific as well as Atlantic regions would make it possible to determine if the proposed latitude gradient in Rn emissions is superior to a uniform flux scenario.

Gupta, Mohan L.↗

NASA NTRS2022

A New Machine Learning Based Analysis for Improving Satellite Retrieved Atmospheric Composition Data: OMI SO2 as an Example

Despite recent progress, satellite retrievals of anthropogenic SO2 still suffer from relatively low signal-tonoise ratios. In this study, we demonstrate a new machine learning data analysis method to improve the quality of satellite SO2 products. In the absence of large ground-truth datasets for SO2, we start from SO2 slant column densities (SCDs) retrieved from the Ozone Monitoring Instrument (OMI) using a data-driven, physically based algorithm and calculate the ratio between the SCD and the root mean square (rms) of the fitting residuals for each pixel. To build the training data, we select presumably clean pixels with small SCD / rms ratios (SRRs) and set their target SCDs to zero. For polluted pixels with relatively large SRRs, we set the target to the original retrieved SCDs. We then train neural networks (NNs) to reproduce the target SCDs using predictors including SRRs for individual pixels, solar zenith, viewing zenith and phase angles, scene reflectivity, and O3 column amounts, as well as the monthly mean SRRs. For data analysis, we employ two NNs: (1) one trained daily to produce analyzed SO2 SCDs for polluted pixels each day and (2) the other trained once every month to produce analyzed SCDs for less polluted pixels for the entire month. Test results for 2005 show that our method can significantly reduce noise and artifacts over background regions. Over polluted areas, the monthly mean NN-analyzed and original SCDs generally agree to within ±15 %, indicating that our method can retain SO2 signals in the original retrievals except for large volcanic eruptions. This is further confirmed by running both the NN-analyzed and original SCDs through a topdown emission algorithm to estimate the annual SO2 emissions for ∼ 500 anthropogenic sources, with the two datasets yielding similar results. We also explore two alternative approaches to the NN-based analysis method. In one, we employ a simple linear interpolation model to analyze the original SCD retrievals. In the other, we develop a PCA–NN algorithm that uses OMI measured radiances, transformed and dimension-reduced with a principal component analysis (PCA) technique, as inputs to NNs for SO2 SCD retrievals. While the linear model and the PCA–NN algorithm can reduce retrieval noise, they both underestimate SO2 over polluted areas. Overall, the results presented here demonstrate that our new data analysis method can significantly improve the quality of existing OMI SO2 retrievals. The method can potentially be adapted for other sensors and/or species and enhance the value of satellite data in air quality research and applications.

Can Li↗

NASA NTRS1994

Higher harmonic control analysis for vibration reduction of helicopter rotor systems

An advanced higher harmonic control (HHC) analysis has been developed and applied to investigate its effect on vibration reduction levels, blade and control system fatigue loads, rotor performance, and power requirements of servo-actuators. The analysis is based on a finite element method in space and time. A nonlinear time domain unsteady aerodynamic model, based on the indicial response formulation, is used to calculate the airloads. The rotor induced inflow is computed using a free wake model. The vehicle trim controls and blade steady responses are solved as one coupled solution using a modified Newton method. A linear frequency-domain quasi-steady transfer matrix is used to relate the harmonics of the vibratory hub loads to the harmonics of the HHC inputs. Optimal HHC is calculated from the minimization of the vibratory hub loads expressed in term of a quadratic performance index. Predicted vibratory hub shears are correlated with wind tunnel data. The fixed-gain HHC controller suppresses completely the vibratory hub shears for most of steady or quasi-steady flight conditions. HHC actuator amplitudes and power increase significantly at high forward speeds (above 100 knots). Due to the applied HHC, the blade torsional stresses and control loads are increased substantially. For flight conditions where the blades are stalled considerably, the HHC input-output model is quite nonlinear. For such cases, the adaptive-gain controller is effective in suppressing vibratory hub loads, even though HHC may actually increase stall areas on the rotor disk. The fixed-gain controller performs poorly for such flight conditions. Comparison study of different rotor systems indicates that a soft-inplane hingeless rotor requires less actuator power at high speeds (above 130 knots) than an articulated rotor, and a stiff-inplane hingeless rotor generally requires more actuator power than an articulated or a soft-inplane hingeless rotor. Parametric studies for a hingeless rotor operating in a transition flight regime and for an articulated rotor operating at the level-flight boundary (high speed and high thrust conditions) indicate that blade parameters including flap, lag, torsion stiffness distributions, linear pretwist, chordwise offset of center-of-mass from elastic axis and chordwise offset of elastic axis from aerodynamic center can be selected to minimize the actuator power requirements for HHC.

Nguyen, Khanh Q.↗

NASA NTRS1975

The computer in shell stability analysis

Some examples in which the high-speed computer has been used to improve the static stability analysis capability for general shells are examined. The fundamental concepts of static stability are reviewed with emphasis on the differences between linear bifurcation buckling and nonlinear collapse. The analysis is limited to the stability of conservative systems. Three examples are considered. The problem of cylinders subjected to bending loads is used as an example to illustrate that a simple structure can have a sufficiently complicated nonlinear behavior to require a computer analysis for accurate results. An analysis of the problems involved in the modeling of stiffening elements in plate and shell structures illustrates the necessity that the analyst recognizes all important deformation modes. The stability analysis of the Skylab structure indicates the size of problems that can be solved with current state-of-the-art capability.

Almroth, B. O.↗

NASA NTRS1977

Improved ground truth geoid for the GEOS-3 calibration area

The purpose of this investigation is to develop methods and procedures are reported for computing a detailed geoid to be used as geodetic ground truth for the calibration and verification of GEOS-3 altimeter data. The technique developed is based on rectifying the best available detailed geoid so that the rectified geoid will have correct scale, orientation, shape and position with respect to the geocenter. The approach involved the development of a mathematical model based on a second degree polynomial, in rectangular Cartesian coordinates, describing the geoid undulations at the control stations. A generalized least squares solution was obtained for the polynomial which describes the variation of the undulation differences between the control stations geoid and the gravimetric geoid. Three rectified geoid were determined. These geoids correspond to three sets of tracking station data: (1) WFC/C-band data; (2) GSFC/C-band data; and (3) OSU-275 data. The absolute accuracy of these rectified geoids is linearly correlated with the uncertainties of the tracking station coordinates and, to a certain extent, with those of the detailed geoid being rectified.

Mourad, A. G.↗

NASA NTRS1986

The evaluation of the OSGLR algorithm for restructurable controls

The detection and isolation of commercial aircraft control surface and actuator failures using the orthogonal series generalized likelihood ratio (OSGLR) test was evaluated. The OSGLR algorithm was chosen as the most promising algorithm based on a preliminary evaluation of three failure detection and isolation (FDI) algorithms (the detection filter, the generalized likelihood ratio test, and the OSGLR test) and a survey of the literature. One difficulty of analytic FDI techniques and the OSGLR algorithm in particular is their sensitivity to modeling errors. Therefore, methods of improving the robustness of the algorithm were examined with the incorporation of age-weighting into the algorithm being the most effective approach, significantly reducing the sensitivity of the algorithm to modeling errors. The steady-state implementation of the algorithm based on a single cruise linear model was evaluated using a nonlinear simulation of a C-130 aircraft. A number of off-nominal no-failure flight conditions including maneuvers, nonzero flap deflections, different turbulence levels and steady winds were tested. Based on the no-failure decision functions produced by off-nominal flight conditions, the failure detection performance at the nominal flight condition was determined. The extension of the algorithm to a wider flight envelope by scheduling the linear models used by the algorithm on dynamic pressure and flap deflection was also considered. Since simply scheduling the linear models over the entire flight envelope is unlikely to be adequate, scheduling of the steady-state implentation of the algorithm was briefly investigated.

Bonnice, W. F.↗

NASA NTRS1989

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗

NASA NTRS1996

Algorithms for Automatic Alignment of Arrays

Aggregate data objects (such as arrays) are distributed across the processor memories when compiling a data-parallel language for a distributed-memory machine. The mapping determines the amount of communication needed to bring operands of parallel operations into alignment with each other. A common approach is to break the mapping into two stages: an alignment that maps all the objects to an abstract template, followed by a distribution that maps the template to the processors. This paper describes algorithms for solving the various facets of the alignment problem: axis and stride alignment, static and mobile offset alignment, and replication labeling. We show that optimal axis and stride alignment is NP-complete for general program graphs, and give a heuristic method that can explore the space of possible solutions in a number of ways. We show that some of these strategies can give better solutions than a simple greedy approach proposed earlier. We also show how local graph contractions can reduce the size of the problem significantly without changing the best solution. This allows more complex and effective heuristics to be used. We show how to model the static offset alignment problem using linear programming, and we show that loop-dependent mobile offset alignment is sometimes necessary for optimum performance. We describe an algorithm with for determining mobile alignments for objects within do loops. We also identify situations in which replicated alignment is either required by the program itself or can be used to improve performance. We describe an algorithm based on network flow that replicates objects so as to minimize the total amount of broadcast communication in replication.

Chatterjee, Siddhartha↗

NASA NTRS2013

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗

NASA NTRS1974

A direct approach to the design of linear multivariable systems

Design of multivariable systems is considered and design procedures are formulated in the light of the most recent work on model matching. The word model matching is used exclusively to mean matching the input-output behavior of two systems. The term is used in the frequency domain to indicate the comparison of two transfer matrices containing transfer functions as elements. Design methods where non-interaction is not used as a criteria were studied. Two design methods are considered. The first method of design is based solely upon the specification of generalized error coefficients for each individual transfer function of the overall system transfer matrix. The second design method is called the pole fixing method because all the system poles are fixed at preassigned positions. The zeros of terms either above or below the diagonal are partially fixed via steady state error coefficients. The advantages and disadvantages of each method are discussed and an example is worked to demonstrate their uses. The special cases of triangular decoupling and minimum constraints are discussed.

Agrawal, B. L.↗

NASA NTRS2001

Transient Reliability Analysis Capability Developed for CARES/Life

The CARES/Life software developed at the NASA Glenn Research Center provides a general-purpose design tool that predicts the probability of the failure of a ceramic component as a function of its time in service. This award-winning software has been widely used by U.S. industry to establish the reliability and life of a brittle material (e.g., ceramic, intermetallic, and graphite) structures in a wide variety of 21st century applications.Present capabilities of the NASA CARES/Life code include probabilistic life prediction of ceramic components subjected to fast fracture, slow crack growth (stress corrosion), and cyclic fatigue failure modes. Currently, this code can compute the time-dependent reliability of ceramic structures subjected to simple time-dependent loading. For example, in slow crack growth failure conditions CARES/Life can handle sustained and linearly increasing time-dependent loads, whereas in cyclic fatigue applications various types of repetitive constant-amplitude loads can be accounted for. However, in real applications applied loads are rarely that simple but vary with time in more complex ways such as engine startup, shutdown, and dynamic and vibrational loads. In addition, when a given component is subjected to transient environmental and or thermal conditions, the material properties also vary with time. A methodology has now been developed to allow the CARES/Life computer code to perform reliability analysis of ceramic components undergoing transient thermal and mechanical loading. This means that CARES/Life will be able to analyze finite element models of ceramic components that simulate dynamic engine operating conditions. The methodology developed is generalized to account for material property variation (on strength distribution and fatigue) as a function of temperature. This allows CARES/Life to analyze components undergoing rapid temperature change in other words, components undergoing thermal shock. In addition, the capability has been developed to perform reliability analysis for components that undergo proof testing involving transient loads. This methodology was developed for environmentally assisted crack growth (crack growth as a function of time and loading), but it will be extended to account for cyclic fatigue (crack growth as a function of load cycles) as well.

Nemeth, Noel N.↗