Engineering Papers⌕ Search

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.

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

At least 721 records · Page 40

Experimental Investigation of a Hypersonic Glider Configuration at a Mach Number of 6 and at Full-Scale Reynolds Numbers

The aerodynamic characteristics of a hypersonic glider configuration, consisting of a slender ogive cylinder with three highly swept wings, spaced 120 apart, with the wing chord equal to the body length, were investigated experimentally at a Mach number of 6 and at Reynolds numbers from 6 to 16 million. The objectives were to evaluate the theoretical procedures which had been used to estimate the performance of the glider, and also to evaluate the characteristics of the glider itself. A principal question concerned the viscous drag at full-scale Reynolds number, there being a large difference between the total drags for laminar and turbulent boundary layers. It was found that the procedures which had been applied for estimating minimum drag, drag due to lift, lift curve slope, and center of pressure were generally accurate within 10 percent. An important exception was the non-linear contribution to the lift coefficient which had been represented by a Newtonian term. Experimentally, the lift curve was nearly linear within the angle-of-attack range up to 10 deg. This error affected the estimated lift-drag ratio. The minimum drag measurements indicated that substantial amounts of turbulent boundary layer were present on all models tested, over a range of surface roughness from 5 microinches maximum to 200 microinches maximum. In fact, the minimum drag coefficients were nearly independent of the surface smoothness and fell between the estimated values for turbulent and laminar boundary layers, but closer to the turbulent value. At the highest test Reynolds numbers and at large angles of attack, there was some indication that the skin friction of the rough models was being increased by the surface roughness. At full-scale Reynolds number, the maximum lift-drag ratio with a leading edge of practical diameter (from the standpoint of leading-edge heating) was 4.0. The configuration was statically and dynamically stable in pitch and yaw, and the center of pressure was less than 2-percent length ahead of the centroid of plan-form area.

Seiff, Alvin↗

Long-Wave Infrared Dyson Spectrometer

Preliminary results are presented for an ultra compact long-wave infrared slit spectrometer based on the dyson concentric design. The dyson spectrometer has been integrated in a dewar environment with a quantum well infrared photodetecor (QWIP), concave electron beam fabricated diffraction grating and ultra precision slit. The entire system is cooled to cryogenic temperatures to maximize signal to noise ratio performance, hence eliminating thermal signal from transmissive elements and internal stray light. All of this is done while maintaining QWIP thermal control. A general description is given of the spectrometer, alignment technique and predicated performance. The spectrometer has been designed for optimal performance with respect to smile and keystone distortion. A spectral calibration is performed with NIST traceable targets. A 2-point non-uniformity correction is performed with a precision blackbody source to provide radiometric accuracy. Preliminary laboratory results show excellent agreement with modeled noise equivalent delta temperature and detector linearity over a broad temperature range.

imaging↗

Understanding Climate-Vegetation Interactions in Global Rainforests Through a GP-Tree Analysis

The tropical rainforests are the largest reserves of terrestrial carbon and, therefore, the future of these rainforests is a question that is of immense importance in the geoscience research community. With the recent severe Amazonian droughts in 2005 and 2010 and on-going drought in the Congo region for more than two decades, there is growing concern that these forests could succumb to precipitation reduction, causing extensive carbon release and feedback to the carbon cycle. However, there is no single ecosystem model that quantifies the relationship between vegetation health in these rainforests and climatic factors. Small scale studies have used statistical correlation measure and simple linear regression to model climate-vegetation interactions, but suffer from the lack of comprehensive data representation as well as simplistic assumptions about dependency of the target on the covariates. In this paper we use genetic programming (GP) based symbolic regression for discovering equations that govern the vegetation climate dynamics in the rainforests. Expecting micro-regions within the rainforests to have unique characteristics compared to the overall general characteristics, we use a modified regression-tree based hierarchical partitioning of the space to build individual models for each partition. The discovery of these equations reveal very interesting characteristics about the Amazon and the Congo rainforests. Our method GP-tree shows that the rainforests exhibit tremendous resiliency in the face of extreme climatic events by adapting to changing conditions.

genetic programming↗

Instabilities in line-driven stellar winds. I - Dependence on perturbation wavelength

An analytical study is presented of the stability of absorption line-driven flows, such as found in stellar winds, in the presence of small-amplitude disturbances. A generalized calculation of the perturbed direct extinction force is performed and the evolution of the perturbation into the nonlinear regime and the dynamical results of perturbation in the scattered radiation field are examined. An expression is derived for the wavenumber variation of the perturbed line force from the milieu of nonoverlapping lines which have a power-law distribution in opacity. A linear dispersion analysis is carried out to model the growth and propagation of radiative-acoustic waves in absorption line-driven flows, which are found unstable to perturbations that may be smaller than the Sobolev length. No damping mechanism was found that would eliminate the absorption line-driven flows.

Owocki, S. P.↗

PIFCGT: A PIF autopilot design program for general aviation aircraft

This report documents the PIFCGT computer program. In FORTRAN, PIFCGT is a computer design aid for determing Proportional-Integral-Filter (PIF) control laws for aircraft autopilots implemented with a Command Generator Tracker (CGT). The program uses Linear-Quadratic-Regulator synthesis algorithms to determine feedback gains, and includes software to solve the feedforward matrix equation which is useful in determining the command generator tracker feedforward gains. The program accepts aerodynamic stability derivatives and computes the corresponding aerodynamic linear model. The nine autopilot modes that can be designed include four maneuver modes (ROLL SEL, PITCH SEL, HDG SEL, ALT SEL), four final approach models (APR GS, APR LOCI, APR LOCR, APR LOCP), and a BETA HOLD mode. The program has been compiled and executed on a CDC computer.

Broussard, J. R.↗

Review and evaluation of recent developments in melic inlet dynamic flow distortion prediction and computer program documentation and user's manual estimating maximum instantaneous inlet flow distortion from steady-state total pressure measurements with full, limited, or no dynamic data

A review of the Melick method of inlet flow dynamic distortion prediction by statistical means is provided. These developments include the general Melick approach with full dynamic measurements, a limited dynamic measurement approach, and a turbulence modelling approach which requires no dynamic rms pressure fluctuation measurements. These modifications are evaluated by comparing predicted and measured peak instantaneous distortion levels from provisional inlet data sets. A nonlinear mean-line following vortex model is proposed and evaluated as a potential criterion for improving the peak instantaneous distortion map generated from the conventional linear vortex of the Melick method. The model is simplified to a series of linear vortex segments which lay along the mean line. Maps generated with this new approach are compared with conventionally generated maps, as well as measured peak instantaneous maps. Inlet data sets include subsonic, transonic, and supersonic inlets under various flight conditions.

Schweikhard, W. G.↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, J. P.↗

A methodology for computing uncertainty bounds of multivariable systems based on sector stability theory concepts

The application of a sector-based stability theory approach to the formulation of useful uncertainty descriptions for linear, time-invariant, multivariable systems is explored. A review of basic sector properties and sector-based approach are presented first. The sector-based approach is then applied to several general forms of parameter uncertainty to investigate its advantages and limitations. The results indicate that the sector uncertainty bound can be used effectively to evaluate the impact of parameter uncertainties on the frequency response of the design model. Inherent conservatism is a potential limitation of the sector-based approach, especially for highly dependent uncertain parameters. In addition, the representation of the system dynamics can affect the amount of conservatism reflected in the sector bound. Careful application of the model can help to reduce this conservatism, however, and the solution approach has some degrees of freedom that may be further exploited to reduce the conservatism.

Waszak, Martin R.↗

Feedback Control of Plasma Instabilities with Charged Particle Beams and Study of Plasma Turbulence

A new non-perturbing technique for feedback control of plasma instabilities has been developed in the Columbia Linear Machine (CLM). The feedback control scheme involves the injection of a feedback modulated ion beam as a remote suppressor. The ion beam was obtained from a compact ion beam source which was developed for this purpose. A Langmuir probe was used as the feedback sensor. The feedback controller consisted of a phase-shifter and amplifiers. This technique was demonstrated by stabilizing various plasma instabilities to the background noise level, like the trapped particle instability, the ExB instability and the ion-temperature-gradient (ITG) driven instability. An important feature of this scheme is that the injected ion beam is non-perturbing to the plasma equilibrium parameters. The robustness of this feedback stabilization scheme was also investigated. The principal result is that the scheme is fairly robust, tolerating about 100% variation about the nominal parameter values. Next, this scheme is extended to the unsolved general problem of controlling multimode plasma instabilities simultaneously with a single sensor-suppressor pair. A single sensor-suppressor pair of feedback probes is desirable to reduce the perturbation caused by the probes. Two plasma instabilities the ExB and the ITG modes, were simultaneously stabilized. A simple 'state' feedback type method was used where more state information was generated from the single sensor Langmuir probe by appropriate signal processing, in this case, by differentiation. This proof-of-principle experiment demonstrated for the first time that by designing a more sophisticated electronic feedback controller, many plasma instabilities may be simultaneously controlled. Simple theoretical models showed generally good agreement with the feedback experimental results. On a parallel research front, a better understanding of the saturated state of a plasma instability was sought partly with the help of feedback. A plasma instability is usually observed in its saturated state and appears as a single feature in the frequency spectrum with a single azimuthal and parallel wavenumbers. The physics of the non-zero spectral width was investigated in detail because the finite spectral width can cause "turbulent" transport. One aspect of the "turbulence" was investigated by obtaining the scaling of the linear growth rate of the instabilities with the fluctuation levels. The linear growth rates were measured with the established gated feedback technique. The research showed that the ExB instability evolves into a quasi-coherent state when the fluctuation level is high. The coherent aspects were studied with a bispectral analysis. Moreover, the single spectral feature was discovered to be actually composed of a few radial harmonics. The radial harmonics play a role in the nonlinear saturation of the instability via three-wave coupling.

Tham, Philip Kin-Wah↗

On the sputtering of binary compounds

A simple physical model is presented to describe some aspects of the sputtering of compound targets. In particular, expressions are developed for the partial sputtering yields for binary systems in terms of the elemental sputtering rates, the stoichiometric concentrations and surface binding energy. The partial yields depend non-linearly on the bulk target concentrations. Comparison of the theoretical predictions with the data on sputtering of PtSi, NiSi and Cu3Au indicates that the general features are well described.

Haff, P. K.↗

Dynamic response and stability analysis of flexible, multibody systems

A general version of Lagrange's equations, including auxiliary nonholonomic, rheonomic conditions of constraint, is used in the dynamic simulation and stability analysis of interconnected flexible bodies. Modeling of the nonlinear flexible/rigid dynamic coupling effects, the interaction forces/torques, and the elastic deformation effects is discussed. A digital computer program is developed to obtain time-domain solution for the nonlinear response of systems represented as a collection of individual bodies, numerical linearization of system-governing equations, time-domain solution for the perturbation response about a nominal state, and a frequency-domain stability analysis corresponding to the linearization. The digital simulation code is employed to study the dynamic behavior of a typical satellite and a spacecraft with deployable experiment booms.

Bodley, C. S.↗

The 3-D inelastic analysis methods for hot section components: Brief description

Advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate yet more cost-effective analysis of components subjected to severe thermal gradients and loads in the presence of mechanical loads, with steep stress and strain gradients are being developed. Anisotropy, time and temperature dependent plasticity and creep effects are also addressed. The approach is to develop four different theories, one linear and three higher order theories (polynomial function, special function, general function). The theories are progressively more complex from linear to general function in order to provide streamlined analysis capability with increasing accuracy for each hot section component and for different parts of the same component according to the severity of the local stress, strain and temperature gradients associated with hot spots, cooling holes and surface coating cracks. To further enhance the computational effectiveness, the higher order theories will have embedded singularities (cooling passages, for example) in the generic modeling region. Each of the four theories consists of three formulation models derivable from independent theoretical formulations. These formulation models are based on: (1) mechanics of materials; (2) special finite elements; and (3) an advanced formulation to be recommended by the contractor.

Chamis, C. C.↗

Frequency Domain Control Design of Large Space Structures: A Practical Approach

Requirements indicate the need for much larger, more accurate, and in some cases, very dynamic satellites. Large control system bandwidths are needed to meet accuracy and response requirements while maintaining tight control over appendage oscillations. Studies in recent years have shown that linear quadratic Gaussian (LQG) controllers can achieve the desired performance if the system is linearized and if the system model is accurate. Results of an LQG controller applied to a single axis satellite with large solar arrays, are given. A reduced order model (ROM) comprises rigid body motion with dominant structural modes. Optimal control and estimation gains are calculated based on an extremely conservative 0.0005 critical damping ratio. In order to examine stability characteristics, single-input single-output (SISO) frequency response concepts are generalized to develop a method of displaying open loop frequency response of multi-input multi-output (MIMO) control system.

Harding, R.↗

Power flow for spacecraft power systems

A method for constructing the generalized system-level admittance matrix for use with a Newton-Raphson power flow is presented. The network modeling technique presented does not use the standard pi-equivalent models, which assume a lossless return path, for the transmission line and transformer. If the return path cannot be assumed lossless, then the standard algorithms for constructing the system admittance matrix cannot be used. The method presented here uses concepts from linear graph theory to combine network modules to form the system-level admittance matrix. The modeling technique is presented, and the resulting matrix is used with a standard Newton-Raphson power flow to calculate all system voltages and current (power) flows.

Halpin, S. M.↗

Systems engineering studies of on-orbit assembly operation

While the practice of construction has a long history, the underlying theory of construction is relatively young. Very little has been documented as to techniques of logistic support, construction planning, construction scheduling, construction testing, and inspection. The lack of 'systems approaches' to construction processes is certainly one of the most serious roadblocks to the construction of space structures. System engineering research efforts at CSC are aimed at developing concepts and tools which contribute to a systems theory of space construction. The research is also aimed at providing means for trade-offs of design parameters for other research areas in CSC. Systems engineering activity at CSC has divided space construction into the areas of orbital assembly, lunar base construction, interplanetary transport vehicle construction, and Mars base construction. A brief summary of recent results is given. Several models for 'launch-on-time' were developed. Launch-on-time is a critical concept to the assembly of such Earth-orbiting structures as the Space Station Freedom, and to planetary orbiters such as the Mars transfer vehicle. CSC has developed a launch vehicle selection model which uses linear programming to find optimal combinations of launch vehicles of various sizes (Atlas, Titan, Shuttles, HLLV's) to support SEI missions. Recently, the Center developed a cost trade-off model for studying on orbit assembly logistics. With this model it was determined that the most effective size of the HLLV would be in the range of 120 to 200 metric tons to LEO, which is consistent with the choices of General Stafford's Synthesis Group Report. A second-generation Dynamic Construction Activities Model ('DYCAM') process model has been under development, based on our past results in interruptability and our initial DYCAM model. This second-generation model is built on the paradigm of knowledge-based expert systems. It is aimed at providing answers to two questions: (1) what are some necessary or sufficient conditions for judging conceptual designs of spacecraft?, and (2) can a methodology be formulated such that these conditions may be used to provide computer-aided tools for evaluating conceptual designs and planning for space assembly sequences? Early simulation results indicate that the DYCAM model has a clear ability to emulate and simulate human orbital construction processes.

Morgenthaler, George W.↗

A fixed-lag Kalman smoother for retrospective data assimilation

Analyses incorporating future data, as well as current and past data, are termed retrospective analyses. In this paper, the fixed-lag Kalman smoother (FLKS) is proposed as a means of providing retrospective analysis capability in data assimilation. The FLKS is a direct generalization of the Kalman filter. It incorporates all data observed up to and including some fixed amount of time past each analysis time. A computationally efficient form of the FLKS is derived. A simple scalar examination of the FLKS demonstrates that incorporating future data improves analyses the most in the presence of dynamical instabilities, for accurate models and for accurate observations. An implementation of the FLKS for two-dimensional linear shallow-water model corroborates the scalar analysis. The numerical experiments also demonstrate the ability of the FLKS to propagate information upstream as well as downstream, thus improving analysis quality substantially in data voids.

Cohn, Stephen E.↗

Effects of Shocks on Emission from Central Engines of Active Galactic Nuclei. I

In this paper we show that perturbations of the accretion flow within the central engines of some active galactic nuclei (AGNS) are likely to form shock waves in the accreting plasma. Such shocks, which may be either collisional or collisionless, can contribute to the observed high-energy temporal and spectral variability. Our rationale is the following: Observations show that the continuum emission probably originates in an optically thin, hot plasma in the AGN central engine. The flux and spectrum from this hot plasma varies significantly over light crossing timescales. Several authors have suggested that macroscopic perturbations contained within this plasma are the sources of this variability. In order to produce the observed emission the perturbations must be radiatively coupled with the optically thin hot matter and must also move with high velocities. We suggest that shocks, which can be very effective in randomizing the bulk motion of the perturbations, are responsible for this coupling. Shocks should form in the central engine, because the temperatures and magnetic fields are probably reduced below their virial values by radiative dissipation. Perturbations moving at Keplerian speeds, or strong non-linear excitations, result in supersonic and super-Alfvenic velocities leading to shock waves within the hot plasma. We show that even a perturbation smaller than the emitting region can form a shock that significantly modifies the continuum emission in an AGN, and that the spectral and temporal variability from such a shock generally resembles those of radio-quiet AGNS. As an example, the shock inducing perturbation in our model is a small main-sequence star, the capturing and eventual accretion of which are known to be a plausible process. We argue that shocks in the central engine may also provide a natural triggering mechanism for the "cold" component of Guilbert & Rees two-phase medium and an efficient mecha- nism for angular momentum transfer. Current and future missions, such as ASCA, XTE, XMM, AXAF, and ASTRO-E may determine the importance of shock-related emission from the central engines of AGNS.

Sivron, R.↗