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 919 records · Page 51

Modeling Progressive Failure of Bonded Joints Using a Single Joint Finite Element

Enhanced finite elements are elements with an embedded analytical solution which can capture detailed local fields, enabling more efficient, mesh-independent finite element analysis. In the present study, an enhanced finite element is applied to generate a general framework capable of modeling an array of joint types. The joint field equations are derived using the principle of minimum potential energy, and the resulting solutions for the displacement fields are used to generate shape functions and a stiffness matrix for a single joint finite element. This single finite element thus captures the detailed stress and strain fields within the bonded joint, but it can function within a broader structural finite element model. The costs associated with a fine mesh of the joint can thus be avoided while still obtaining a detailed solution for the joint. Additionally, the capability to model non-linear adhesive constitutive behavior has been included within the method, and progressive failure of the adhesive can be modeled by using a strain-based failure criteria and re-sizing the joint as the adhesive fails. Results of the model compare favorably with experimental and finite element results.

Stapleton, Scott E.↗

Integrated control-system design via generalized LQG (GLQG) theory

Thirty years of control systems research has produced an enormous body of theoretical results in feedback synthesis. Yet such results see relatively little practical application, and there remains an unsettling gap between classical single-loop techniques (Nyquist, Bode, root locus, pole placement) and modern multivariable approaches (LQG and H infinity theory). Large scale, complex systems, such as high performance aircraft and flexible space structures, now demand efficient, reliable design of multivariable feedback controllers which optimally tradeoff performance against modeling accuracy, bandwidth, sensor noise, actuator power, and control law complexity. A methodology is described which encompasses numerous practical design constraints within a single unified formulation. The approach, which is based upon coupled systems or modified Riccati and Lyapunov equations, encompasses time-domain linear-quadratic-Gaussian theory and frequency-domain H theory, as well as classical objectives such as gain and phase margin via the Nyquist circle criterion. In addition, this approach encompasses the optimal projection approach to reduced-order controller design. The current status of the overall theory will be reviewed including both continuous-time and discrete-time (sampled-data) formulations.

Bernstein, Dennis S.↗

Generation and Calibration of Linear Models of Aircraft with Highly Coupled Aeroelastic and Flight Dynamics

The lightweight structures and unconventional configurations being considered for the next generation of aircraft mean that any effort to predict or control the flight dynamics is impacted by the structural dynamics. One of the most severe forms of coupling between aeroelasticity and flight dynamics is an instability called body freedom flutter. The existing tools often assume a relatively weak effect of structural dynamics on the flight dynamics, and are therefore incapable of modeling strong interactions like body freedom flutter. A method of combining different sources of data traditionally used for aeroelasticity and flight dynamics is described by reconciling many of the differences between these models. By building upon past modeling efforts, a level of familiarity in the approach is achieved. Generally the differences from the traditional approaches are subtle but significant. The traditional frequency domain flutter model in a modal coordinate system is converted to a form consistent with a time domain flight dynamics model. The time domain rational function approximation about a non-inertial coordinate system and the unique constraints for the conversion between the inertial and non-inertial coordinate systems are discussed. A consistent transformation of the states of aeroelastic models to flight dynamics models is derived, which enables the integration of data from higher fidelity computational fluid dynamics models or wind-tunnel testing. The present method of integrating multidisciplinary data was used to create models that compare well with X-56A flight-test data, including conditions past the flutter speed.

Ouellette, Jeffrey↗

Effect of wing planform and canard location and geometry on the longitudinal aerodynamic characteristics of a close-coupled canard wing model at subsonic speeds

A generalized wind-tunnel model with canard and wing planforms typical of highly maneuverable aircraft was tested in the Langley 7- by 10-foot high-speed tunnel at a Mach number of 0.30 to determine the effect of canard location, canard size, wing sweep, and canard strake on canard-wing interference to high angles of attack. The major results of this investigation may be summarized as follows: the high-canard configuration (excluding the canard strake and canard flap), for both the 60 deg and 44 deg swept leading-edge wings, produced the highest maximum lift coefficient and the most linear pitching-moment curves; substantially larger gains in the canard lift and total lift were obtained by adding a strake to the canard located below the wing chord plane rather than by adding a strake to the canard located above the wing chord plane.

Gloss, B. B.↗

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗

Stability and control characteristics of a monoplanar missile configuration with triform-tail-fin arrangements at Mach numbers from 1.70 to 2.86

A wind-tunnel missile model with either a lower vertical tail fin with a pair of horizontal fins having 0 deg, 22.5 deg, or 30 deg dihedral or an upper vertical tail fin with horizontal fins having 0 deg, -22.5 deg, or -30 deg dihedral was investigated. The results indicated that those configurations with horizontal fins at or below the horizontal plane had nearly linear pitching-moment characteristics, while those with the horizontal fins above the horizontal plane experienced pitch-up which increased with increasing horizontal-fin-dihedral angle. At zero angle of attack, the configurations were directionally stable at most test Mach numbers. Generally, those configurations with the upper vertical fin had positive effective dihedral at zero angle of attack, while those with he lower vertical fin had negative effective dihedral. For roll control, three deflected tail fins produced more total roll control than two horizontal fins. For yaw control, three tail fins deflected equally or differentially produced more total yaw control than the single vertical fin.

Lamb, M.↗

Stability of Detached Solidification

Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.

Mazuruk, K.↗

Development and validation of a piloted simulation of a helicopter and external sling load

A generalized, real time, piloted, visual simulation of a single rotor helicopter, suspension system, and external load is described and validated for the full flight envelope of the U.S. Army CH-54 helicopter and cargo container as an example. The mathematical model described uses modified nonlinear classical rotor theory for both the main rotor and tail rotor, nonlinear fuselage aerodynamics, an elastic suspension system, nonlinear load aerodynamics, and a loadground contact model. The implementation of the mathematical model on a large digital computing system is described, and validation of the simulation is discussed. The mathematical model is validated by comparing measured flight data with simulated data, by comparing linearized system matrices, eigenvalues, and eigenvectors with manufacturers' data, and by the subjective comparison of handling characteristics by experienced pilots. A visual landing display system for use in simulation which generates the pilot's forward looking real world display was examined and a special head up, down looking load/landing zone display is described.

Shaughnessy, J. D.↗

Modeling the effects of latitudinal gradients in stellar winds, with application to the solar wind

A steady, axisymmetric, quasi-radial, global model previously developed for stellar winds with embedded magnetic fields has been extended to include latitudinal gradient effects on the azimuthal velocity and magnetic field. The linear results at large radii are presented for large-amplitude latitudinal variations in the radial magnetic field, mass loss rate, and radial velocity of the wind. The magnetohydrodynamic (MHD) equations predict meridional flows that develop naturally from internal magnetic stresses. The flows open flux tubes in the star's equatorial plane, redistributing mass and magnetic flux as a function of stellar latitude. The plasma spins up to conserve angular momentum in fields and plasma. The results are generally applicable to stellar winds (including radiatively driven winds), provided that the internal structure is not dominated by rotation. The asymptotic solutions do not explicitly depend on the form of the energy equation, although the assumed O(1) state which drives these solutions depends on the deposition of energy and momentum throughout the wind.

Nerney, S.↗

Exact finite element method analysis of viscoelastic tapered structures to transient loads

A general method is presented for determining the dynamic torsional/axial response of linear structures composed of either tapered bars or shafts to transient excitations. The method consists of formulating and solving the dynamic problem in the Laplace transform domain by the finite element method and obtaining the response by a numerical inversion of the transformed solution. The derivation of the torsional and axial stiffness matrices is based on the exact solution of the transformed governing equation of motion, and it consequently leads to the exact solution of the problem. The solution permits treatment of the most practical cases of linear tapered bars and shafts, and employs modeling of structures with only one element per member which reduces the number of degrees of freedom involved. The effects of external viscous or internal viscoelastic damping are also taken into account.

Spyrakos, Constantine Chris↗

Linear and nonlinear interpretation of CV-580 lightning data

Numerical models developed for the study of lightning strike data acquired by in-flight aircraft are applied to the data measured on the CV-580. The basic technique used is the three dimensional time domain finite difference solution of Maxwell's equations. Both linear and nonlinear models are used in the analysis. In the linear model, the lightning channel and the aircraft are assumed to form a linear time invariant system. A transfer function technique can then be used to study the response of the aircraft to a given lightning strike current. Conversely, the lightning current can be inferred from the measured response. In the nonlinear model, the conductivity of air in the vicinity of the aircraft is calculated and incorporated into the solution of the Maxwell's equations. The nonlinear model thus simulates corona formation and air breakdown. Results obtained from the models are in reasonable agreement with the measured data. This study provides another validation of the models and increases confidence that the models may be used to predict aircraft response to any general lightning strike.

Ng, Poh H.↗

Performance of a Bounce-Averaged Global Model of Super-Thermal Electron Transport in the Earth's Magnetic Field

In this paper, we report the results of our recent research on the application of a multiprocessor Cray T916 supercomputer in modeling super-thermal electron transport in the earth's magnetic field. In general, this mathematical model requires numerical solution of a system of partial differential equations. The code we use for this model is moderately vectorized. By using Amdahl's Law for vector processors, it can be verified that the code is about 60% vectorized on a Cray computer. Speedup factors on the order of 2.5 were obtained compared to the unvectorized code. In the following sections, we discuss the methodology of improving the code. In addition to our goal of optimizing the code for solution on the Cray computer, we had the goal of scalability in mind. Scalability combines the concepts of portabilty with near-linear speedup. Specifically, a scalable program is one whose performance is portable across many different architectures with differing numbers of processors for many different problem sizes. Though we have access to a Cray at this time, the goal was to also have code which would run well on a variety of architectures.

McGuire, Tim↗

Water Isotopes in Precipitation: Data/Model Comparison for Present-Day and Past Climates

Variations of HDO and H2O-18 concentrations are observed in precipitation both on a geographical and on a temporal basis. These variations, resulting from successive isotopic fractionation processes at each phase change of water during its atmospheric cycle, are well documented through the IAEA/WMO network and other sources. Isotope concentrations are, in middle and high latitudes, linearly related to the annual mean temperature at the precipitation site. Paleoclimatologists have used this relationship to infer paleotemperatures from isotope paleodata extractable from ice cores, deep groundwater and other such sources. For this application to be valid, however, the spatial relationship must also hold in time at a given location as the location undergoes a series of climatic changes. Progress in water isotope modeling aimed at examining and evaluating this assumption has been recently reviewed with a focus on polar regions and, more specifically, on Greenland. This article was largely based on the results obtained using the isotopic version of the NASA/GISS Atmospheric General Circulation Model (AGCM) fitted with isotope tracer diagnostics. We extend this review in comparing the results of two different isotopic AGCMs (NASA/GISS and ECHAM) and in examining, with a more global perspective, the validity of the above assumption, i.e. the equivalence of the spatial and temporal isotope-temperature relationship. We also examine recent progress made in modeling the relationship between the conditions prevailing in moisture source regions for precipitation and the deuterium-excess of that precipitation.

Jouzel, J.↗

Integrated Modeling Tools for Thermal Analysis and Applications

Integrated modeling of spacecraft systems is a rapidly evolving area in which multidisciplinary models are developed to design and analyze spacecraft configurations. These models are especially important in the early design stages where rapid trades between subsystems can substantially impact design decisions. Integrated modeling is one of the cornerstones of two of NASA's planned missions in the Origins Program -- the Next Generation Space Telescope (NGST) and the Space Interferometry Mission (SIM). Common modeling tools for control design and opto-mechanical analysis have recently emerged and are becoming increasingly widely used. A discipline that has been somewhat less integrated, but is nevertheless of critical concern for high precision optical instruments, is thermal analysis and design. A major factor contributing to this mild estrangement is that the modeling philosophies and objectives for structural and thermal systems typically do not coincide. Consequently the tools that are used in these discplines suffer a degree of incompatibility, each having developed along their own evolutionary path. Although standard thermal tools have worked relatively well in the past. integration with other disciplines requires revisiting modeling assumptions and solution methods. Over the past several years we have been developing a MATLAB based integrated modeling tool called IMOS (Integrated Modeling of Optical Systems) which integrates many aspects of structural, optical, control and dynamical analysis disciplines. Recent efforts have included developing a thermal modeling and analysis capability, which is the subject of this article. Currently, the IMOS thermal suite contains steady state and transient heat equation solvers, and the ability to set up the linear conduction network from an IMOS finite element model. The IMOS code generates linear conduction elements associated with plates and beams/rods of the thermal network directly from the finite element structural model. Conductances for temperature varying materials are accommodated. This capability both streamlines the process of developing the thermal model from the finite element model, and also makes the structural and thermal models compatible in the sense that each structural node is associated with a thermal node. This is particularly useful when the purpose of the analysis is to predict structural deformations due to thermal loads. The steady state solver uses a restricted step size Newton method, and the transient solver is an adaptive step size implicit method applicable to general differential algebraic systems. Temperature dependent conductances and capacitances are accommodated by the solvers. In addition to discussing the modeling and solution methods. applications where the thermal modeling is "in the loop" with sensitivity analysis, optimization and optical performance drawn from our experiences with the Space Interferometry Mission (SIM), and the Next Generation Space Telescope (NGST) are presented.

Milman, Mark H.↗

Stability Assessment and Tuning of an Adaptively Augmented Classical Controller for Launch Vehicle Flight Control

Recently, a robust and practical adaptive control scheme for launch vehicles [ [1] has been introduced. It augments a classical controller with a real-time loop-gain adaptation, and it is therefore called Adaptive Augmentation Control (AAC). The loop-gain will be increased from the nominal design when the tracking error between the (filtered) output and the (filtered) command trajectory is large; whereas it will be decreased when excitation of flex or sloshing modes are detected. There is a need to determine the range and rate of the loop-gain adaptation in order to retain (exponential) stability, which is critical in vehicle operation, and to develop some theoretically based heuristic tuning methods for the adaptive law gain parameters. The classical launch vehicle flight controller design technics are based on gain-scheduling, whereby the launch vehicle dynamics model is linearized at selected operating points along the nominal tracking command trajectory, and Linear Time-Invariant (LTI) controller design techniques are employed to ensure asymptotic stability of the tracking error dynamics, typically by meeting some prescribed Gain Margin (GM) and Phase Margin (PM) specifications. The controller gains at the design points are then scheduled, tuned and sometimes interpolated to achieve good performance and stability robustness under external disturbances (e.g. winds) and structural perturbations (e.g. vehicle modeling errors). While the GM does give a bound for loop-gain variation without losing stability, it is for constant dispersions of the loop-gain because the GM is based on frequency-domain analysis, which is applicable only for LTI systems. The real-time adaptive loop-gain variation of the AAC effectively renders the closed-loop system a time-varying system, for which it is well-known that the LTI system stability criterion is neither necessary nor sufficient when applying to a Linear Time-Varying (LTV) system in a frozen-time fashion. Therefore, a generalized stability metric for time-varying loop=gain perturbations is needed for the AAC.

VanZwieten, Tannen↗

Expanding the application of soil moisture monitoring systems through regression-based transformation

Relative to other geophysical variables, soil moisture (SM) estimates derived from land surface models (LSMs) and land data assimilation systems (LDAS) are difficult to transfer between platforms and applications. This difficulty stems from the highly model-dependent nature of LSM SM estimates and differences in the vertical support of discretized SM values. As a result, operational SM estimates generated by one LSM (or LDAS) cannot generally be directly applied to a hydrologic monitoring or forecast system designed around a second LSM. This lack of transferability is particularly problematic for LDAS applications, where the time, expertise, and computational resources required to generate an operational LDAS analysis cannot be practically duplicated for every LSM-specific application. Here, we develop a set of simple regression tools for translating SM estimates between LSMs and multiple LDAS analyses. Results demonstrate that simple multivariate linear regression - utilizing independent variables based on multi-layer and temporally lagged SM estimates - can significantly improve upon baseline transformation approaches using direct percentile matching. The proposed regression approaches are effective for both the LSM-to-LSM and LDAS-to-LDAS transformation of multi-layer SM percentiles. Application of this approach will expand the utility of existing, high-quality (but LSM-specific) operational sources of SM information like the NASA Soil Moisture Active Passive Level-4 Soil Moisture product.

Soil Moisture↗

Automatic voice recognition using traditional and artificial neural network approaches

The main objective of this research is to develop an algorithm for isolated-word recognition. This research is focused on digital signal analysis rather than linguistic analysis of speech. Features extraction is carried out by applying a Linear Predictive Coding (LPC) algorithm with order of 10. Continuous-word and speaker independent recognition will be considered in future study after accomplishing this isolated word research. To examine the similarity between the reference and the training sets, two approaches are explored. The first is implementing traditional pattern recognition techniques where a dynamic time warping algorithm is applied to align the two sets and calculate the probability of matching by measuring the Euclidean distance between the two sets. The second is implementing a backpropagation artificial neural net model with three layers as the pattern classifier. The adaptation rule implemented in this network is the generalized least mean square (LMS) rule. The first approach has been accomplished. A vocabulary of 50 words was selected and tested. The accuracy of the algorithm was found to be around 85 percent. The second approach is in progress at the present time.

Botros, Nazeih M.↗

Limits of adaptation, residual interferences

Methods of determining linear residual wall interference appear to be well established theoretically; however they need to be validated, for example by comparative studies of test data on the same model in different adaptive-wall wind tunnels as well as in passive, ventilated-wall tunnels. The GARTEur CAST 7 and the CAST 10/DOA 2 investigations are excellent examples of such comparative studies. Results to date in both one-variable and two-variable methods for nonlinear wall interference indicate that a great deal more research and validation are required. The status in 2D flow is advanced over that in 3D flow as is the case generally with adaptive-wall development. Nevertheless, it is now well established that for transonic testing with extensive supercritical flow present, significant wall interference is likely to exist in conventional ventilated test sections. Consequently, residual correction procedures require further development hand-in-hand with further adaptive-wall development.

Mokry, Miroslav↗