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 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Collisional excitation of interstellar molecules - Linear molecules CO, CS, OCS, and HC3N

Rates for excitation of CO, CS, OCS, and HC3N by collisions with He atoms and H2 molecules are presented. These have been obtained from extensive theoretical calculations using methods of known reliability. In general, rates necessary for modeling interstellar molecular clouds with kinetic temperatures to 100 K have been included. Methods for fitting these rates in more compact form are discussed.

Green, S.↗

NASA NTRS1981

The Total Synthesis Problem of linear multivariable control. II - Unity feedback and the design morphism

Zames (1981) has observed that there is, in general, no 'separation principle' to guarantee optimality of a division between control law design and filtering of plant uncertainty. Peczkowski and Sain (1978) have solved a model matching problem using transfer functions. Taking into consideration this investigation, Peczkowski et al. (1979) proposed the Total Synthesis Problem (TSP), wherein both the command/output-response and command/control-response are to be synthesized, subject to the plant constraint. The TSP concept can be subdivided into a Nominal Design Problem (NDP), which is not dependent upon specific controller structures, and a Feedback Synthesis Problem (FSP), which is. Gejji (1980) found that NDP was characterized in terms of the plant structural matrices and a single, 'good' transfer function matrix. Sain et al. (1981) have extended this NDP work. The present investigation is concerned with a study of FSP for the unity feedback case. NDP, together with feedback synthesis, is understood as a Total Synthesis Problem.

Sain, M. K.↗

NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗

NASA NTRS2010

A Coupled Aeroelastic Model for Launch Vehicle Stability Analysis

A technique for incorporating distributed aerodynamic normal forces and aeroelastic coupling effects into a stability analysis model of a launch vehicle is presented. The formulation augments the linear state-space launch vehicle plant dynamics that are compactly derived as a system of coupled linear differential equations representing small angular and translational perturbations of the rigid body, nozzle, and sloshing propellant coupled with normal vibration of a set of orthogonal modes. The interaction of generalized forces due to aeroelastic coupling and thrust can be expressed as a set of augmenting non-diagonal stiffness and damping matrices in modal coordinates with no penalty on system order. While the eigenvalues of the structural response in the presence of thrust and aeroelastic forcing can be predicted at a given flight condition independent of the remaining degrees of freedom, the coupled model provides confidence in closed-loop stability in the presence of rigid-body, slosh, and actuator dynamics. Simulation results are presented that characterize the coupled dynamic response of the Ares I launch vehicle and the impact of aeroelasticity on control system stability margins.

Orr, Jeb S.↗

NASA NTRS1972

On the synthesis of multivariable systems.

Description of a general synthesis procedure for the compensation of linear multivariable systems through the combined use of dynamic feed-forward compensation and linear state variable feedback. Applications of the synthesis algorithm presented to problems of decoupling and exact model matching illustrate its use.

Wolovich, W. A.↗

NASA NTRS1985

A sequential linear optimization approach for controller design

A linear optimization approach with a simple real arithmetic algorithm is presented for reliable controller design and vibration suppression of flexible structures. Using first order sensitivity of the system eigenvalues with respect to the design parameters in conjunction with a continuation procedure, the method converts a nonlinear optimization problem into a maximization problem with linear inequality constraints. The method of linear programming is then applied to solve the converted linear optimization problem. The general efficiency of the linear programming approach allows the method to handle structural optimization problems with a large number of inequality constraints on the design vector. The method is demonstrated using a truss beam finite element model for the optimal sizing and placement of active/passive-structural members for damping augmentation. Results using both the sequential linear optimization approach and nonlinear optimization are presented and compared. The insensitivity to initial conditions of the linear optimization approach is also demonstrated.

Horta, L. G.↗

NASA NTRS1987

Passive stabilization for large space systems

The optimal tuning of multiple tuned-mass dampers for the transient vibration damping of large space structures is investigated. A multidisciplinary approach is used. Structural dynamic techniques are applied to gain physical insight into absorber/structure interaction and to optimize specific cases. Modern control theory and parameter optimization techniques are applied to the general optimization problem. A design procedure for multi-absorber multi-DOF vibration damping problems is presented. Classical dynamic models are extended to investigate the effects of absorber placement, existing structural damping, and absorber cross-coupling on the optimal design synthesis. The control design process for the general optimization problem is formulated as a linear output feedback control problem via the development of a feedback control canonical form. The techniques are applied to sample micro-g and pointing problems on the NASA dual keel space station.

Sesak, J. R.↗

NASA NTRS1989

Probabilistic boundary element method

The purpose of the Probabilistic Structural Analysis Method (PSAM) project is to develop structural analysis capabilities for the design analysis of advanced space propulsion system hardware. The boundary element method (BEM) is used as the basis of the Probabilistic Advanced Analysis Methods (PADAM) which is discussed. The probabilistic BEM code (PBEM) is used to obtain the structural response and sensitivity results to a set of random variables. As such, PBEM performs analogous to other structural analysis codes such as finite elements in the PSAM system. For linear problems, unlike the finite element method (FEM), the BEM governing equations are written at the boundary of the body only, thus, the method eliminates the need to model the volume of the body. However, for general body force problems, a direct condensation of the governing equations to the boundary of the body is not possible and therefore volume modeling is generally required.

Cruse, T. A.↗

NASA NTRS1992

Progress in integrated analysis with adaptive unstructured meshing

Design of lightweight structures and thermal protection systems for hypersonic vehicles depend on accurate prediction of aerothermal loads, structural temperatures and their gradients, and structural deformations and stresses. Concentration is on an alternative meshing technique which generates an entirely new adaptive unstructured mesh based on the solution obtained from the earlier mesh. The technique combined with the finite element method has been shown to significantly improve the efficiency and accuracy of the fluid, thermal, and structural analyses. Current capability of the adaptive unstructured meshing technique for the integrated fluid-thermal-structural analysis is described first. The technique was extended to transient thermal analysis of structures with time-dependent adaptive meshing to capture the detailed temperature response with a minimum number of unknowns and computational cost. Both linear and higher-order finite elements are implemented to demonstrate the generality of the technique and to investigate their solution accuracy. Currently, the adaptive meshing technique is being developed for plane structures that can be modeled with membrane elements and built-up structures modeled with membrane and bending elements. The capability of the technique to these different disciplinary problems is demonstrated by several examples.

Dechaumphai, Pramote↗

NASA NTRS1961

Pressure Distributions and Wave Drag Due to Two-Dimensional Fabrication-Type Surface Roughness on an Ogive Cylinder at Mach Numbers of 1.61 and 2.01

An investigation has been made at Mach numbers of 1.61 and 2.01 and over a range of free-stream Reynolds number per foot from about 1.2 x 10(exp 6) to 8.3 x 10(exp 6) to determine the pressure distributions and wave drags due to two-dimensional fabrication-type surface roughness. Ten types of surface roughness, including step, wave, crease, and swept configurations were investigated. The tests were made on an ogive cylinder of fineness ratio 12.2, the roughness elements covering the cylindrical portion of the model. The results indicate that wave drag is the major component of the drag due to roughness at supersonic speeds. The pressure distributions over the roughness elements were generally found to be in good agreement with linearized two-dimensional theory except for regions of the elements affected by boundary-layer separation and shock detachment. There was little or no effect of Reynolds number except on the pressures within the regions influenced by separation or shock detachment. Inasmuch as most of the roughness configurations were affected by flow separation and shock detachment, there was generally an effect of Reynolds number on the roughness wave drag. This wave drag decreased as the free-stream Reynolds number was decreased.

Czarnecki, K. R.↗

NASA NTRS1994

Computer modeling the fatigue crack growth rate behavior of metals in corrosive environments

The objective of this task was to develop a method to digitize FCP (fatigue crack propagation) kinetics data, generally presented in terms of extensive da/dN-Delta K pairs, to produce a file for subsequent linear superposition or curve-fitting analysis. The method that was developed is specific to the Numonics 2400 Digitablet and is comparable to commercially available software products as Digimatic(sup TM 4). Experiments demonstrated that the errors introduced by the photocopying of literature data, and digitization, are small compared to those inherent in laboratory methods to characterize FCP in benign and aggressive environments. The digitizing procedure was employed to obtain fifteen crack growth rate data sets for several aerospace alloys in aggressive environments.

Richey, Edward, III↗

NASA NTRS1991

Mixed convective/dynamic roll vortices and their effects on initial wind and temperature profiles

The onset and development of both dynamically and convectively forced boundary layer rolls are studied with linear and nonlinear analyses of a truncated spectral model of shallow Boussinesq flow. Emphasis is given here on the energetics of the dominant roll modes, on the magnitudes of the roll-induced modifications of the initial basic state wind and temperature profiles, and on the sensitivity of the linear stability results to the use of modified profiles as basic states. It is demonstrated that the roll circulations can produce substantial changes to the cross-roll component of the initial wind profile and that significant changes in orientation angle estimates can result from use of a roll-modified profile in the stability analysis. These results demonstrate that roll contributions must be removed from observed background wind profiles before using them to investigate the mechanisms underlying actual secondary flows in the boundary layer. The model is developed quite generally to accept arbitrary basic state wind profiles as dynamic forcing. An Ekman profile is chosen here merely to provide a means for easy comparison with other theoretical boundary layer studies; the ultimate application of the model is to study observed boundary layer profiles. Results of the analytic stability analysis are validated by comparing them with results from a larger linear model. For an appropriate Ekman depth, a complete set of transition curves is given in forcing parameter space for roll modes driven both thermally and dynamically. Preferred orientation angles, horizontal wavelengths and propagation frequencies, as well as energetics and wind profile modifications, are all shown to agree rather well with results from studies on Ekman layers as well as with studies on near-neutral and convective atmospheric boundary layers.

Haack, Tracy↗

NASA NTRS1998

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

NASA NTRS1989

Computational fluid dynamic control

A general technique is presented for modeling fluid, or gas, dynamic systems specifically for the development of control systems. The numerical methods which are generally used in computational fluid dynamics are borrowed to create either continuous-time or discrete-time models of the particular fluid system. The resulting equations can be either left in a nonlinear form, or easily linearized about an operating point. As there are typically very many states in these systems, the usual linear model reduction methods can be used on them to allow a low-order controller to be designed. A simple example is given which typifies many internal flow control problems. The resulting control is termed computational fluid dynamic control.

Hartley, Tom T.↗

NASA NTRS1990

State-variable models of structures having rigid-body modes

In cases where the equations of motion of a structure having rigid-body freedom are cast in state-variable form, generalized state rigid-body modes may be needed. It is possible to find a linearly-independent set of generalized vectors which transform an n x n matrix into the almost-diagonal Jordan form. Attention is presently given to equations governing these generalized eigenvectors, together with illustrative examples of the damped and undamped structure cases.

Craig, Roy R., Jr.↗

NASA NTRS2000

Turbulence and Laminar Structures: Can They Co-Exist?

Schwarzschild first suggested that the laminar structures observed in the high Reynolds number Re = UL/nu approx. = (10(exp 12)) solar photosphere are the result of turbulence rather than a proof of its absence. He reasoned that since turbulence generates large turbulent viscosities nu(sub t) much greater than nu, the "effective" Reynolds number Re = UL/nu(sub t) approx. = O(1). Schwarzschild's argument is, however, incomplete for it assumes that the entire role of the non-linear interactions is to "enhance" viscosity. While this is not true in general, we present a proof of how and why it may occur, thus completing Schwarzschild's argument. We further discuss the fact that the same model non-local turbulence models have been shown to reproduce LES data for a variety of flows pertaining to astrophysics, geophysics and laboratory situations (at a fraction of the time).

Canuto, V. M.↗

NASA NTRS2011

Geochemical Mapping of 4 Vesta Begins

By December, the NASA Dawn spacecraft will have descended to a low altitude mapping orbit (LAMO), where the Gamma Ray and Neutron Detector (GRaND) will acquire global mapping data for up to four months. Measurements by GRaND will help answer elusive questions about how Vesta differentiated and the nature of processes that shaped Vesta s surface. The data will be analyzed to determine the abundances of Mg, Si, Fe, K, Th, and H at a spatial resolution of roughly 300 km full-width-at-half-maximum from a 465 km radius orbit. Thermal and fast neutron counting data will be analyzed to determine the neutron macroscopic absorption cross section and average atomic mass, providing constraints on additional elements, such as Ca and Al. GRaND will quantify the elemental composition of coarse spatial units identified by Dawn s Framing Camera (FC) and the Visible & Infrared Spectrometer (VIR). In addition, GRaND will map the mixing ratio of compositional end members selected from the howardite, eucrite and diogenite (HED) meteorites, determine the relative proportions of plagioclase and mafic minerals, and search for compositions that are absent or under-represented in the meteorite collection. While it is generally thought that Vesta s crust on a regional scale should be well-represented by linear mixing of HED whole-rock compositions, there are hints that Vesta may be more diverse than implied by this model. For example, the discovery of K-rich impact glasses in howardites suggests that K-rich rocks may be present on a portion of Vesta s surface, and the analysis of diogenites indicates considerable variability in the magmatic processes that formed them. The chemical composition of materials within Vesta s south polar structure may provide further clues to how it formed. An impact might have exposed mantle and lower crustal materials, which should have a distinctive compositional signature. We present the analysis of data acquired by GRaND from cruise through the descent to LAMO, including GRaND s sensitivity to different elements and geochemical processes.

Prettyman, Thomas H.↗

NASA NTRS1950

Free-Flight Investigation of the Static and Dynamic Longitudinal Stability Characteristics of 1/3.7-Scale Rocket-Powered Models of the Bell MX-776A

An investigation of the static and dynamic longitudinal stability characteristics of 1/3.7 scale rocket-powered model of the Bell MX-776A has been made for a Mach number range from 0.8 to 1.6. Two models were tested with all control surfaces at 0 degree deflection and centers of gravity located 1/4 and 1/2 body diameters, respectively, ahead of the equivalent design location. Both models were stable about the trim conditions but did not trim at 0 degree angle of attack because of slight constructional asymmetries. The results indicated that the variation of lift and pitching moment was not linear with angle of attack. Both lift-curve slope and pitching-moment-curve slope were of the smallest magnitude near 0 degree angle of attack. In general, an increase in angle of attack was accompanied by a rearward movement of the aerodynamic center as the rear wing moved out of the downwash from the forward surfaces. This characteristic was more pronounced in the transonic region. The dynamic stability in the form of total damping factor varied with normal-force coefficient but was greatest for both models at a Mach number of approximately 1.25. The damping factor was greater at the lower trim normal-force coefficients except at a Mach number of 1.0. At that speed the damping factor was of about the same magnitude for both models. The drag coefficient increased with trim normal-force coefficient and was largest in the transonic region.

Michal, David H.↗