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 649 records · Page 36

Intraplate deformation due to continental collisions: A numerical study of deformation in a thin viscous sheet

A model of crustal deformation from continental collision that involves the penetration of a rigid punch into a deformable sheet is investigated. A linear viscous flow law is used to compute the magnitude and rate of change of crustal thickness, the velocity of mass points, strain rates and their principal axes, modes of deformation, areal changes, and stress. In general, a free lateral boundary reduces the magnitude of changes in crustal thickening by allowing material to more readily escape the advancing punch. The shearing that occurs diagonally in front of the punch terminates in compression or extension depending on whether the lateral boundary is fixed or free. When the ratio of the diameter of the punch to that of the sheet exceeds one-third, the deformation is insenstive to the choice of lateral boundary conditions. When the punch is rigid with sharply defined edges, deformation is concentrated near the punch corners. With non-rigid punches, shearing results in deformation being concentrated near the center of the punch. Variations with respect to linearity and nonlinearity of flow are discussed.

Cohen, S. C.↗

Large liquid rocket engine transient performance simulation system

A simulation system, ROCETS, was designed and developed to allow cost-effective computer predictions of liquid rocket engine transient performance. The system allows a user to generate a simulation of any rocket engine configuration using component modules stored in a library through high-level input commands. The system library currently contains 24 component modules, 57 sub-modules and maps, and 33 system routines and utilities. FORTRAN models from other sources can be operated in the system upon inclusion of interface information on comment cards. Operation of the simulation is simplified for the user by run, execution, and output processors. The simulation system makes available steady-state trim balance, transient operation, and linear partial generation. The system utilizes a modern equation solver for efficient operation of the simulations. Transient integration methods include integral and differential forms for the trapezoidal, first order Gear, and second order Gear corrector equations. A detailed technology test bed engine (TTBE) model was generated to be used as the acceptance test of the simulation system. The general level of model detail was that reflected in the Space Shuttle Main Engine DTM. The model successfully obtained steady-state balance in main stage operation and simulated throttle transients, including engine starts and shutdown. A NASA FORTRAN control model was obtained, ROCETS interface installed in comment cards, and operated with the TTBE model in closed-loop transient mode.

Mason, J. R.↗

Software For Integer Programming

Improved Exploratory Search Technique for Pure Integer Linear Programming Problems (IESIP) program optimizes objective function of variables subject to confining functions or constraints, using discrete optimization or integer programming. Enables rapid solution of problems up to 10 variables in size. Integer programming required for accuracy in modeling systems containing small number of components, distribution of goods, scheduling operations on machine tools, and scheduling production in general. Written in Borland's TURBO Pascal.

Fogle, F. R.↗

Model reduction of nonsquare linear MIMO systems using multipoint matrix continued-fraction expansions

This paper deals with the multipoint Cauer matrix continued-fraction expansion (MCFE) for model reduction of linear multi-input multi-output (MIMO) systems with various numbers of inputs and outputs. A salient feature of the proposed MCFE approach to model reduction of MIMO systems with square transfer matrices is its equivalence to the matrix Pade approximation approach. The Cauer second form of the ordinary MCFE for a square transfer function matrix is generalized in this paper to a multipoint and nonsquare-matrix version. An interesting connection of the multipoint Cauer MCFE method to the multipoint matrix Pade approximation method is established. Also, algorithms for obtaining the reduced-degree matrix-fraction descriptions and reduced-dimensional state-space models from a transfer function matrix via the multipoint Cauer MCFE algorithm are presented. Practical advantages of using the multipoint Cauer MCFE are discussed and a numerical example is provided to illustrate the algorithms.

Guo, Tong-Yi↗

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗

Use of North American and European Air Quality Networks to Evaluate Global Chemistry-Climate Modeling of Surface Ozone

We test the current generation of global chemistry-climate models in their ability to simulate observed, present-day surface ozone. Models are evaluated against hourly surface ozone from 4217 stations in North America and Europe that are averaged over 1 degree by 1 degree grid cells, allowing commensurate model-measurement comparison. Models are generally biased high during all hours of the day and in all regions. Most models simulate the shape of regional summertime diurnal and annual cycles well, correctly matching the timing of hourly (approximately 15:00 local time (LT)) and monthly (mid-June) peak surface ozone abundance. The amplitude of these cycles is less successfully matched. The observed summertime diurnal range (25 ppb) is underestimated in all regions by about 7 parts per billion, and the observed seasonal range (approximately 21 parts per billion) is underestimated by about 5 parts per billion except in the most polluted regions, where it is overestimated by about 5 parts per billion. The models generally match the pattern of the observed summertime ozone enhancement, but they overestimate its magnitude in most regions. Most models capture the observed distribution of extreme episode sizes, correctly showing that about 80 percent of individual extreme events occur in large-scale, multi-day episodes of more than 100 grid cells. The models also match the observed linear relationship between episode size and a measure of episode intensity, which shows increases in ozone abundance by up to 6 parts per billion for larger-sized episodes. We conclude that the skill of the models evaluated here provides confidence in their projections of future surface ozone.

diurnal variations↗

Novel Estimation of Pilot Performance Characteristics

Two mechanisms internal to the pilot that affect performance during a tracking task are: 1) Pilot equalization (i.e. lead/lag); and 2) Pilot gain (i.e. sensitivity to the error signal). For some applications McRuer's Crossover Model can be used to anticipate what equalization will be employed to control a vehicle's dynamics. McRuer also established approximate time delays associated with different types of equalization - the more cognitive processing that is required due to equalization difficulty, the larger the time delay. However, the Crossover Model does not predict what the pilot gain will be. A nonlinear pilot control technique, observed and coined by the authors as 'amplitude clipping', is shown to improve stability, performance, and reduce workload when employed with vehicle dynamics that require high lead compensation by the pilot. Combining linear and nonlinear methods a novel approach is used to measure the pilot control parameters when amplitude clipping is present, allowing precise measurement in real time of key pilot control parameters. Based on the results of an experiment which was designed to probe workload primary drivers, a method is developed that estimates pilot spare capacity from readily observable measures and is tested for generality using multi-axis flight data. This paper documents the initial steps to developing a novel, simple objective metric for assessing pilot workload and its variation over time across a wide variety of tasks. Additionally, it offers a tangible, easily implementable methodology for anticipating a pilot's operating parameters and workload, and an effective design tool. The model shows promise in being able to precisely predict the actual pilot settings and workload, and observed tolerance of pilot parameter variation over the course of operation. Finally, an approach is proposed for generating Cooper-Harper ratings based on the workload and parameter estimation methodology.

Human pilot models↗

Synchrotron sources. I - Extension of theory for small pitch angles.

Some of the features of synchrotron radiation are examined which occur when relativistic particles move at small pitch angles in a magnetic field. The emission from individual particles is computed. The general properties of the emission from distributions of particles with small pitch angles is discussed, and the emissivities and degrees of polarization for several specific source models are presented. For sources which mainly contain particles with small pitch angles, the degree of circular polarization will normally be high, whereas, if the magnetic field is somewhat irregular, the degree of linear polarization will be small. If the particle distribution also has a low-energy cutoff, the emissivity will vary linearly with frequency at low frequencies.

Epstein, R. I.↗

ERTS-1 imagery use in reconnaissance prospecting: Evaluation of commercial utility of ERTS-1 imagery in structural reconnaissance for minerals and petroleum

The author has identified the following significant results. This study was performed to investigate applications of ERTS-1 imagery in commercial reconnaissance for mineral and hydrocarbon resources. ERTS-1 imagery collected over five areas in North America (Montana; Colorado; New Mexico-West Texas; Superior Province, Canada; and North Slope, Alaska) has been analyzed for data content including linears, lineaments, and curvilinear anomalies. Locations of these features were mapped and compared with known locations of mineral and hydrocarbon accumulations. Results were analyzed in the context of a simple-shear, block-coupling model. Data analyses have resulted in detection of new lineaments, some of which may be continental in extent, detection of many curvilinear patterns not generally seen on aerial photos, strong evidence of continental regmatic fracture patterns, and realization that geological features can be explained in terms of a simple-shear, block-coupling model. The conculsions are that ERTS-1 imagery is of great value in photogeologic/geomorphic interpretations of regional features, and the simple-shear, block-coupling model provides a means of relating data from ERTS imagery to structures that have controlled emplacement of ore deposits and hydrocarbon accumulations, thus providing a basis for a new approach for reconnaissance for mineral, uranium, gas, and oil deposits and structures.

Saunders, D. F.↗

The three dimensional motion and stability of a rotating space station: Cable-counterweight configuration

The three dimensional equations of motion for a cable connected space station--counterweight system are developed using a Lagrangian formulation. The system model employed allows for cable and end body damping and restoring effects. The equations are then linearized about the equilibrium motion and nondimensionalized. To first degree, the out-of-plane equations uncouple from the inplane equations. Therefore, the characteristic polynomials for the in-plane and out-of-plane equations are developed and treated separately. From the general in-plane characteristic equation, necessary conditions for stability are obtained. The Routh-Hurwitz necessary and sufficient conditions for stability are derived for the general out-of-plane characteristic equation. Special cases of the in-plane and out-of-plane equations (such as identical end masses, and when the cable is attached to the centers of mass of the two end bodies) are then examined for stability criteria.

Bainum, P. M.↗

Rotorcraft derivative identification from analytical models and flight test data

A general procedure is presented for systematic development of rotorcraft models for use in systems identification, which includes fuselage and rotor degrees of freedom (DOF). Formulations for rigid blade flap and lag as well as the normal mode representation of an elastic blade are developed for hingeless and articulated rotor systems. The method of multiblade coordinates is used to obtain linear constant coefficient state variable models of various levels of approximation. Two of the approximate models, a 6 DOF, are identified from a nonlinear articulated helicopter computer simulation. The results demonstrate the accuracy attainable for each model. Advanced results outline the status of rotorcraft modeling and systems identification and indicate areas that require further investigation.

Molusis, J. A.↗

Overview of magnetic suspension research at Langley Research Center

An overview of research in small- and large-gap magnetic suspension systems at LaRC is presented. The overview is limited to systems which have been built as laboratory models or engineering models. Small-gap systems applications include the Annular Momentum Control Device (AMCD), which is a momentum storage device for the stabilization and control of spacecraft, and the Annular Suspension and Pointing System (ASPS), which is a general purpose pointing mount designed to provide orientation, mechanical isolation, and fine pointing space experiments. These devices are described and control and linearization approaches for the magnetic suspension systems for these devices are discussed. Large-gap systems applications at LaRC have been almost exclusively wind tunnel magnetic suspension systems. A brief description of these efforts is also presented.

Groom, Nelson J.↗

Galactic cosmic ray abundances and spectra behind defined shielding

Linear Energy Transfer (LET) spectra have been measured for lunar missions and for several near Earth orbits ranging from 28 deg to 83 deg inclination. In some of the experiments the flux of Galactic Cosmic Rays (GCR) was determined separately from contributions caused by interactions in the detector material. Results of these experiments are compared to model calculations. The general agreement justifies the use of the model to calculate GCR fluxes. The magnitude of variations caused by solar modulation, geomagnetic shielding, and shielding by matter determined from calculated LET spectra is generally in agreement with experimental data. However, more detailed investigations show that there are some weak points in modeling solar modulation and shielding by material. These points are discussed in more detail.

Heinrich, W.↗

A Thermodynamic Theory of Solid Viscoelasticity: Nonlinear Thermo-viscoelasticity - Part II:

This paper, second in the series of three papers, develops a general, nonlinear, non-isothermal, compressible theory for finite rubber viscoelasticity and specifies it in a form convenient for solving problems important to the rubber, tire, automobile, and air-space industries, among others. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory of differential type has been developed for arbitrary non-isothermal deformations of viscoelastic solids. In this theory, the constitutive equations were presented as the sum of a rubber elastic (equilibrium) and a liquid type viscoelastic (non-equilibrium) terms. These equations have then been simplified using several modeling and simplicity arguments.

Freed, Alan D.↗

Dissipative rendering and neural network control system design

Model-based control system designs are limited by the accuracy of the models of the plant, plant uncertainty, and exogenous signals. Although better models can be obtained with system identification, the models and control designs still have limitations. One approach to reduce the dependency on particular models is to design a set of compensators that will guarantee robust stability to a set of plants. Optimization over the compensator parameters can then be used to get the desired performance. Conservativeness of this approach can be reduced by integrating fundamental properties of the plant models. This is the approach of dissipative control design. Dissipative control designs are based on several variations of the Passivity Theorem, which have been proven for nonlinear/linear and continuous-time/discrete-time systems. These theorems depend not on a specific model of a plant, but on its general dissipative properties. Dissipative control design has found wide applicability in flexible space structures and robotic systems that can be configured to be dissipative. Currently, there is ongoing research to improve the performance of dissipative control designs. For aircraft systems that are not dissipative active control may be used to make them dissipative and then a dissipative control design technique can be used. It is also possible that rendering a system dissipative and dissipative control design may be combined into one step. Furthermore, the transformation of a non-dissipative system to dissipative can be done robustly. One sequential design procedure for finite dimensional linear time-invariant systems has been developed. For nonlinear plants that cannot be controlled adequately with a single linear controller, model-based techniques have additional problems. Nonlinear system identification is still a research topic. Lacking analytical models for model-based design, artificial neural network algorithms have recently received considerable attention. Using their universal approximation property, neural networks have been introduced into nonlinear control designs in several ways. Unfortunately, little work has appeared that analyzes neural network control systems and establishes margins for stability and performance. One approach for this analysis is to set up neural network control systems in the framework presented above. For example, one neural network could be used to render a system to be dissipative, a second strictly dissipative neural network controller could be used to guarantee robust stability.

Gonzalez, Oscar R.↗

DYSCO: A software system for modeling general dynamic systems

The DYSCO program has been under development since 1979. It has been funded by Army and Air Force laboratories and by the Kaman Aerospace Corporation. It is presently available at a number of government and nongovernment installations. It has been used to analyze a very broad range of dynamics problems. A principle feature of the software design of DYSCO is the separation of the executive from the technology. The executive, which controls all the operations, is intelligent in the sense that it knows that its function is to assemble differential equations and to prepare them for solution. The technology library contains FORTRAN routines which perform standard functions, such as, computing the equation coefficients of an element (or component) given the local state at any time. The technology library also contains algorithms and procedures for solving the coupled system equations. The system was designed to allow easy additional of technology to the library. Any linear or nonlinear structural entity, control system, or set of ordinary differential equations may be simply coded and added to the library, as well as algorithms for time or frequency domain solution. The program is described with emphasis on its usefulness in easily modeling unusual concepts and configurations, performing analysis of damage, evaluating new algorithms, and simulating dynamic tests.

Berman, Alex↗

Simultaneous observations of energetic protons close to the bow shock and far upstream

Four upstream energetic proton events (30 keV-75 keV), which were simultaneously observed by the ISEE-1 and ISEE-3 satellites, are investigated. A comparison was made between the absolute flux values about 200 earth radii from the bow shock and the flux values several earth radii in front of the bow shock at two different energies. Close to the bow shock the particle distribution is more or less isotropic and indicates relatively strong scattering of these particles in the upstream wave field. At ISEE-3 between 100 and 200 earth radii upstream from the earth's bow shock, the particles move essentially scatter-free from the general bow shock direction. The proton differential intensity at ISEE-3 is a factor of about 4-15 less than at ISEE-1 at 30 keV. The spectra at ISEE-1 are steeper than the spectra at ISEE-3. Flux ratios and spectra are discussed in terms of a first order Fermi acceleration model with a diffusion coefficient increasing approximately linearly with energy and a free escape boundary at some distance upstream.

Scholer, M.↗

The surface and through crack problems in layered orthotropic plates

An analytical method is developed for a relatively accurate calculation of Stress Intensity Factors in a laminated orthotropic plate containing a through or part-through crack. The laminated plate is assumed to be under bending or membrane loading and the mode 1 problem is considered. First three transverse shear deformation plate theories (Mindlin's displacement based first-order theory, Reissner's stress-based first-order theory, and a simple-higher order theory due to Reddy) are reviewed and examined for homogeneous, laminated and heterogeneous orthotropic plates. Based on a general linear laminated plate theory, a method by which the stress intensity factors can be obtained in orthotropic laminated and heterogeneous plates with a through crack is developed. Examples are given for both symmetrically and unsymmetrically laminated plates and the effects of various material properties on the stress intensity factors are studied. In order to implement the line-spring model which is used later to study the surface crack problem, the corresponding plane elasticity problem of a two-bonded orthotropic plated containing a crack perpendicular to the interface is also considered. Three different crack profiles: an internal crack, an edge crack, and a crack terminating at the interface are considered. The effect of the different material combinations, geometries, and material orthotropy on the stress intensity factors and on the power of stress singularity for a crack terminating at the interface is fully examined. The Line Spring model of Rice and Levy is used for the part-through crack problem. The surface crack is assumed to lie in one of the two-layered laminated orthotropic plates due to the limitation of the available plane strain results. All problems considered are of the mixed boundary value type and are reduced to Cauchy type of singular integral equations which are then solved numerically.

Erdogan, Fazil↗