Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional 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 127 records · Page 7

A dynamic analysis of the motion of a low-wing general aviation aircraft about its calculated equilibrium flat spin mode

The properties of the flat spin mode of a general aviation configuration have been studied through analysis of rotary balance data, numerical simulation, and analytical study of the equilibrium state. The equilibrium state is predicted well from rotary balance data. The variations of yawing moment and pitching moment as functions of sideslip have been shown to be of great importance in obtaining accurate modeling. These dependencies are not presently available with sufficient accuracy from previous tests or theories. The stability of the flat spin mode has been examined extensively using numerical linearization, classical perturbation methods, and reduced order modeling. The stability exhibited by the time histories and the eigenvalue analyses is shown to be strongly dependent on certain static cross derivatives and more so on the dynamic derivatives. Explicit stability criteria are obtained from the reduced order models.

Tischler, M. B.↗

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.↗

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗

On the calculation of turbulent heat transport downstream from an abrupt pipe expansion

A numerical study is reported of flow and heat transfer in the separated flow region created by an abrupt pipe expansion. Computations employed an adaptation of the TEACH-2E computer program with the standard model of turbulence. Emphasis is given to the simulation, from both a physical and numerical viewpoint, of the region in the immediate vicinity of the wall where turbulent transport gives way to molecular conduction and diffusion. Wall resistance laws or wall functions used to bridge this near-wall region are based on the idea that, beyond the viscous sublayer, the turbulent length scale is universal, increasing linearly with distance from the wall. Predictions of expermental data for a diameter ratio of 0.54 show generally encouraging agreement with experiment. At a diameter of 0.43 different trends are discernible between measurement and calculation though this appears to be due to effects unconnected with the wall region studied.

Chieng, C. C.↗

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations. Phase 5B: Analysis of gust alleviation methods and rotor dynamic stability

The effects of various gust alleviation methods on the random blade response in flapping were studied analytically, assuming a rigid rotor support. The analytical model assumes rigid flapping blades with elastic root restraints. Linearized equations which are approximately valid at low lift conditions were used. Because of the interblade coupling from the feedback devices, the method of multiblade generalized coordinates was most convenient and was extended to include coning, differential coning, and warping of the rotor. The numerical examples cover the dynamic stability characteristics as affected by feedback gains of three- to six-bladed rotors. The number of blades had large effects on stability limits and modal time functions at these limits. The random flapping response of the blades to atmospheric turbulence was determined at 1.6 rotor advance ratio using feedback gains below the stability limit. The most effective reduction of the flapping response per unit gain was achieved with a rotor coning angle feedback.

Hohenemser, K. H.↗

Probability Bounds Analysis Applied to Multi-Purpose Crew Vehicle Nonlinearity

The Multi-Purpose Crew Vehicle (MPCV) Program Orion vehicle finite element model (FEM) was updated based on a modal test performed by Lockheed Martin. Due to nonlinearity observed in the test results, linear low force level (LL) and high force level (HL) FEMs were developed for use during various Space Launch System (SLS) flight regimes depending on expected forcing levels. Uncertainty models were derived for the combined MPCV and MPCV Stage Adaptor LL and HL Hurty/Craig-Bampton (HCB) components based on the MPCV structural test article Configuration 4 modal test-analysis correlation results. Subsequently, system-level uncertainty quantification analyses were performed using both models for various SLS flight configurations to determine the impact of the nonlinearity on important system metrics. The system metrics included both transfer functions associated with attitude control and dynamic loads associated with aerodynamic buffeting during ascent. In each case, an independent Monte Carlo (MC) analysis was performed, and no attempt was made to combine the results. The Hybrid Parametric Variation (HPV) method was used to develop the LL and HL MPCV HCB uncertainty models. The HPV method provides both parametric and non-parametric components of uncertainty. The non-parametric uncertainty accounts for the difference in model-form between the linearized analytical model and the corresponding linearized component test results in the form of mode shapes and frequencies at that force level. This linear model-form uncertainty is implemented in the HPV method using random matrix theory. However, the HPV uncertainty models developed for the linear LL and HL MPCV components do not account for the nonlinearity in the MPCV. With respect to the linearized models, this nonlinearity is also an uncertainty in model form, but in this case, it must be treated independently as an epistemic uncertainty. It represents a lack of knowledge, in contrast to an aleatory uncertainty due to the randomness of a variable. In the case of an epistemic variable, the true value is unknown, only the interval within which it lies is known. Epistemic uncertainty can be reduced with increased knowledge, while in general, aleatory uncertainty cannot. This work combines the epistemic uncertainty due to the MPCV nonlinearity with the parametric and non-parametric uncertainty within the HPV method using a second order propagation approach. The LL and HL test data is augmented with surrogate test data derived from a nonlinear MPCV representation. The impact of the MPCV nonlinearity on system response statistics is determined using a series of cumulative distribution functions in the form of a horsetail plot, or p-box. This results in an interval of probabilities for a specific response value, or an interval of response values at a specific probability.

Daniel C Kammer↗

Quantitative models of magnetic and electric fields in the magnetosphere

In order to represent the magnetic field B in the magnetosphere various auxiliary functions can be used: the current density, the scalar potential, toroidal and poloidal potentials, and Euler potentials -- or else, the components of B may be expanded directly. The most versatile among the linear representations is the one based on toroidal and poloidal potentials; it has seen relatively little use in the past but appears to be the most promising one for future work. Other classifications of models include simple testbed models vs. comprehensive ones and analytical vs. numerical representations. The electric field E in the magnetosphere is generally assumed to vary only slowly and to be orthogonal to B, allowing the use of a scalar potential which may be deduced from observations in the ionosphere, from the shape of the plasmapause, or from particle observations in synchronous orbits.

Stern, D. P.↗

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

A model of the human observer in failure detection tasks

A model for the human observer in failure detection tasks is proposed which consists of two stages: a linear estimator and a decision mechanism. The estimator is a Kalman filter, and the decision mechanism, which is based on Wald's sequential analysis, leads to a decision function which is the integration of the filter residuals. The final result is a simple detection system which depends on only three parameters, and the sensitivity of the model to these parameters is analyzed. The results of an experiment designed to test the validity of the model are reported. The question of open and closed decision intervals as well as the generalization of the model to more complicated cases is discussed.

Gai, E. G.↗

Jet Surface Interaction-Scrubbing Noise

Generation of sound due to scrubbing of a jet flow past a nearby solid surface is investigated within the framework of the generalized acoustic analogy theory. The analysis applies to the boundary layer noise generated at and near a wall, and excludes the scattered noise component that is produced at the leading or the trailing edge. While compressibility effects are relatively unimportant at very low Mach numbers, frictional heat generation and thermal gradient normal to the surface could play important roles in generation and propagation of sound in high speed jets of practical interest. A general expression is given for the spectral density of the far-field sound as governed by the variable density Pridmore- Brown equation. The propagation Green's function should be solved numerically starting with the boundary conditions on the surface and subject to specified mean velocity and temperature profiles between the surface and the observer. The equivalent sources of aerodynamic sound are associated with non-linear momentum flux and enthalpy flux terms that appear in the linearized Navier-Stokes equations. These multi-pole sources should be modeled and evaluated with input from a Reynolds-Averaged Navier-Stokes (RANS) solver with an appropriate turbulence model.

Khavaran, Abbas↗

A general unified approach to modelling switching-converter power stages

A method for modeling switching-converter power stages is developed, whose starting point is the unified state-space representation of the switched networks and whose end result is either a complete state-space description or its equivalent small-signal low-frequency linear circuit model. A new canonical circuit model is proposed, whose fixed topology contains all the essential input-output and control properties of any dc-to-dc switching converter, regardless of its detailed configuration, and by which different converters can be characterized in the form of a table conveniently stored in a computer data bank to provide a useful tool for computer aided design and optimization. The new canonical circuit model predicts that, in general, switching action introduces both zeros and poles into the duty ratio to output transfer function in addition to those from the effective filter network.

Middlebrook, R. D.↗

Dynamic Reconstruction and Multivariable Control for Force-Actuated, Thin Facesheet Adaptive Optics

The Multiple Mirror Telescope (MMT) under development at the University of Arizona takes a new approach in adaptive optics placing a large (0.65 m) force-actuated, thin facesheet deformable mirror at the secondary of an astronomical telescope, thus reducing the effects of emissivity which are important in IR astronomy. However, The large size of the mirror and low stiffness actuators used drive the natural frequencies of the mirror down into the bandwidth of the atmospheric distortion. Conventional adaptive optics takes a quasi-static approach to controlling the, deformable mirror. However, flexibility within the control bandwidth calls for a new approach to adaptive optics. Dynamic influence functions are used to characterize the influence of each actuator on the surface of the deformable mirror. A linearized model of atmospheric distortion is combined with dynamic influence functions to produce a dynamic reconstructor. This dynamic reconstructor is recognized as an optimal control problem. Solving the optimal control problem for a system with hundreds of actuators and sensors is formidable. Exploiting the circularly symmetric geometry of the mirror, and a suitable model of atmospheric distortion, the control problem is divided into a number of smaller decoupled control problems using circulant matrix theory. A hierarchic control scheme which seeks to emulate the quasi-static control approach that is generally used in adaptive optics is compared to the proposed dynamic reconstruction technique. Although dynamic reconstruction requires somewhat more computational power to implement, it achieves better performance with less power usage, and is less sensitive than the hierarchic technique.

Grocott, Simon C. O.↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Numerical Results - 2

We present numerical results of the diffusion coefficients (DCs) in the coupled diffusion model derived in the preceding paper for a semiconductor quantum well. These include self and mutual DCs in the general two-component case, as well as density- and temperature-related DCs under the single-component approximation. The results are analyzed from the viewpoint of free Fermi gas theory with many-body effects incorporated. We discuss in detail the dependence of these DCs on densities and temperatures in order to identify different roles played by the free carrier contributions including carrier statistics and carrier-LO phonon scattering, and many-body corrections including bandgap renormalization and electron-hole (e-h) scattering. In the general two-component case, it is found that the self- and mutual- diffusion coefficients are determined mainly by the free carrier contributions, but with significant many-body corrections near the critical density. Carrier-LO phonon scattering is dominant at low density, but e-h scattering becomes important in determining their density dependence above the critical electron density. In the single-component case, it is found that many-body effects suppress the density coefficients but enhance the temperature coefficients. The modification is of the order of 10% and reaches a maximum of over 20% for the density coefficients. Overall, temperature elevation enhances the diffusive capability or DCs of carriers linearly, and such an enhancement grows with density. Finally, the complete dataset of various DCs as functions of carrier densities and temperatures provides necessary ingredients for future applications of the model to various spatially inhomogeneous optoelectronic devices.

Li, Jianzhong↗

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.↗

A Self-Stabilizing Synchronization Protocol for Arbitrary Digraphs

This paper presents a self-stabilizing distributed clock synchronization protocol in the absence of faults in the system. It is focused on the distributed clock synchronization of an arbitrary, non-partitioned digraph ranging from fully connected to 1-connected networks of nodes while allowing for differences in the network elements. This protocol does not rely on assumptions about the initial state of the system, other than the presence of at least one node, and no central clock or a centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. There is no theoretical limit on the maximum number of participating nodes. The only constraint on the behavior of the node is that the interactions with other nodes are restricted to defined links and interfaces. This protocol deterministically converges within a time bound that is a linear function of the self-stabilization period. We present an outline of a deductive proof of the correctness of the protocol. A bounded model of the protocol was mechanically verified for a variety of topologies. Results of the mechanical proof of the correctness of the protocol are provided. The model checking results have verified the correctness of the protocol as they apply to the networks with unidirectional and bidirectional links. In addition, the results confirm the claims of determinism and linear convergence. As a result, we conjecture that the protocol solves the general case of this problem. We also present several variations of the protocol and discuss that this synchronization protocol is indeed an emergent system.

Malekpour, Mahyar R.↗

Linear polarization of light by two wheat canopies measured at many view angles

The linear polarization and reflection of visible light by wheat as a function of sun-view directions, crop development stage, and wavelength were examined. Two-hundred spectra were taken continuously in wave-lengths from 0.45 to 0.72 Micron in 33 view directions using an Exotech model 20C spectroradiometer six meters above two wheat canopies in the boot and fully headed maturity stages. The analysis results show that the amount of linearly polarized light from the wheat canopies is greatest in the blue spectral region and decreases gradually with increasing wavelength. The results also show that the linearly polarized light from the canopies is generally greatest in the azimuth direction of the Sun and tends toward zero as the view direction tends toward the direction of the hot spot or anti-solar point. It is demonstrated that the single, angle of incidence of sunlight on the leaf, explains almost all of the variation of the amount of polarized light with Sun-view direction.

Bauer, M. E.↗

Current Limitations on VLBI Accuracy

The contribution of VLBI to geophysics and geodesy arises from its ability to measure distances between stations in a network and to determine the orientation of stations in the network as well as the orientation of the network with respect to the external reference frame of extragalactic radio objects. Integrating nearly two decades of observations provides useful information about station positions and velocities and the orientation of the Earth, but the complications of the real world and the limitations of observing, modeling and analysis prevent recovery of all effects. Of the factors that limit the accuracy of seemingly straightforward geodetic parameters, the neutral propagation medium has been subject to the greatest scrutiny, but the treatment of the mapping function, the wet component and spatial/temporal inhomogeneities is still improving. These affect both the terrestrial scale and consistency over time. The modeling of non-secular site motions (tides and loading) has increased in sophistication, but there are some differences between the models and the observations. VLBI antennas are massive objects, so their behavior is quite unlike GPS monuments, but antenna deformations add some (generally) unmodeled signal. Radio sources used in geodetic VLBI observations are selected for strength and (relative) absence of structure, but apparent changes in position can leak into geodetic parameters. A linear rate of change of baseline or site parameters is the simplest model and its error improves with time span. However, in most cases the VLBI data distribution is insufficient to look for real non-linear behavior that might affect the average rate. A few sites have multiple VLBI antennas, and some show small differences in rate. VLBI intrinsically measures relative positions and velocities, but individual site positions and velocities are generally more useful. The creation of the VLBI terrestrial reference frame, which transforms relative information into individual results, is an empirical process that has intrinsic errors. While UT1 is uniquely measured by VLBI, the geographical distribution and availability of VLBI stations, especially in the southern hemisphere, and the consistency of the VLBI terrestrial reference frame may limit the accuracy of Earth orientation measurements. The effects of particular error sources on geodetic and geophysical parameters derived from VLBI data will be illustrated.

Ma, Chopo↗