Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear systems”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Aeroelasticity and structural optimization of composite helicopter rotor blades with swept tips

This report describes the development of an aeroelastic analysis capability for composite helicopter rotor blades with straight and swept tips, and its application to the simulation of helicopter vibration reduction through structural optimization. A new aeroelastic model is developed in this study which is suitable for composite rotor blades with swept tips in hover and in forward flight. The hingeless blade is modeled by beam type finite elements. A single finite element is used to model the swept tip. Arbitrary cross-sectional shape, generally anisotropic material behavior, transverse shears and out-of-plane warping are included in the blade model. The nonlinear equations of motion, derived using Hamilton's principle, are based on a moderate deflection theory. Composite blade cross-sectbnal properties are calculated by a separate linear, two-dimensional cross section analysis. The aerodynamic loads are obtained from quasi-steady, incompressible aerodynamics, based on an implicit formulation. The trim and steady state blade aeroelastic response are solved in a fully coupled manner. In forward flight, where the blade equations of motion are periodic, the coupled trim-aeroelastic response solution is obtained from the harmonic balance method. Subsequently, the periodic system is linearized about the steady state response, and its stability is determined from Floquet theory.

Yuan, K. A.↗
NASA NTRS1994

Certification of windshear performance with RTCA class D radomes

Superposition testing of detection range performance forms a digital signal for input into a simulation of signal and data processing equipment and algorithms to be employed in a sensor system for advanced warning of hazardous windshear. For suitable pulse-Doppler radar, recording of the digital data at the input to the digital signal processor furnishes a realistic operational scenario and environmentally responsive clutter signal including all sidelobe clutter, ground moving target indications (GMTI), and large signal spurious due to mainbeam clutter and/or RFI respective of the urban airport clutter and aircraft scenarios (approach and landing antenna pointing). For linear radar system processes, a signal at the same point in the process from a hazard phenomena may be calculated from models of the scattering phenomena, for example, as represented in fine 3 dimensional reflectivity and velocity grid structures. Superposition testing furnishes a competing signal environment for detection and warning time performance confirmation of phenomena uncontrollable in a natural environment.

Mathews, Bruce D.↗
NASA NTRS1992

Autonomous space processor for orbital debris advanced design project in support of solar system exploration

This paper is regarding a project in the Advanced Design Program at the University of Arizona. The project is named the Autonomous Space Processor for Orbital Debris (ASPOD) and is a NASA/Universities Space Research Association (USRA) sponsored design project. The development of ASPOD and the students' abilities in designing and building a prototype spacecraft are the ultimate goals of this project. This year's focus entailed the development of a secondary robotic arm and end-effector to work in tandem with an existent arm in the removal of orbital debris. The new arm features the introduction of composite materials and a linear drive system, thus producing a light-weight and more accurate prototype. The main characteristic of the end-effector design is that it incorporates all of the motors and gearing internally, thus not subjecting them to the harsh space environment. Furthermore, the arm and the end-effector are automated by a control system with positional feedback. This system is composed of magnetic and optical encoders connected to a 486 PC via two servo-motor controller cards. Programming a series of basic routines and sub-routines has allowed the ASPOD prototype to become more autonomous. The new system is expected to perform specified tasks with a positional accuracy of 0.5 cm.

Ramohalli, Kumar↗
NASA NTRS1992

Aerodynamic shape optimization via sensitivity analysis on decomposed computational domains

Direct and iterative method considered to be most applicable to large systems of linear equations arising in discrete sensitivity analysis are assessed. Based on a single-domain grid, computations are performed using a banded matrix solver and an iterative solver, the generalized minimum residual (GMRES) method. The banded matrix solver is found to be generally the most economical method for those applications where the number of right-hand sides is large (i.e., a large number of design variables or a large number of adjoint vectors). For systems of equations that are too large to be solved by direct methods, an approach is proposed whereby the computational domain is divided into small subdomains, and each subdomain is solved separately.

Eleshaky, Mohamed E.↗
NASA NTRS1996

Robust Stabilization of Uncertain Systems Based on Energy Dissipation Concepts

Robust stability conditions obtained through generalization of the notion of energy dissipation in physical systems are discussed in this report. Linear time-invariant (LTI) systems which dissipate energy corresponding to quadratic power functions are characterized in the time-domain and the frequency-domain, in terms of linear matrix inequalities (LMls) and algebraic Riccati equations (ARE's). A novel characterization of strictly dissipative LTI systems is introduced in this report. Sufficient conditions in terms of dissipativity and strict dissipativity are presented for (1) stability of the feedback interconnection of dissipative LTI systems, (2) stability of dissipative LTI systems with memoryless feedback nonlinearities, and (3) quadratic stability of uncertain linear systems. It is demonstrated that the framework of dissipative LTI systems investigated in this report unifies and extends small gain, passivity, and sector conditions for stability. Techniques for selecting power functions for characterization of uncertain plants and robust controller synthesis based on these stability results are introduced. A spring-mass-damper example is used to illustrate the application of these methods for robust controller synthesis.

Gupta, Sandeep↗
NASA NTRS1982

On conforming mixed finite element methods for incompressible viscous flow problems

The application of conforming mixed finite element methods to obtain approximate solutions of linearized Navier-Stokes equations is examined. Attention is given to the convergence rates of various finite element approximations of the pressure and the velocity field. The optimality of the convergence rates are addressed in terms of comparisons of the approximation convergence to a smooth solution in relation to the best approximation available for the finite element space used. Consideration is also devoted to techniques for efficient use of a Gaussian elimination algorithm to obtain a solution to a system of linear algebraic equations derived by finite element discretizations of linear partial differential equations.

Gunzburger, M. D↗
NASA NTRS1977

Development of a parameter optimization technique for the design of automatic control systems

Parameter optimization techniques for the design of linear automatic control systems that are applicable to both continuous and digital systems are described. The model performance index is used as the optimization criterion because of the physical insight that can be attached to it. The design emphasis is to start with the simplest system configuration that experience indicates would be practical. Design parameters are specified, and a digital computer program is used to select that set of parameter values which minimizes the performance index. The resulting design is examined, and complexity, through the use of more complex information processing or more feedback paths, is added only if performance fails to meet operational specifications. System performance specifications are assumed to be such that the desired step function time response of the system can be inferred.

Whitaker, P. H.↗
NASA NTRS1990

Linear and non-linear dynamic models of a geared rotor-bearing system

A three degree of freedom non-linear model of a geared rotor-bearing system with gear backlash and radial clearances in rolling element bearings is proposed here. This reduced order model can be used to describe the transverse-torsional motion of the system. It is justified by comparing the eigen solutions yielded by corresponding linear model with the finite element method results. Nature of nonlinearities in bearings is examined and two approximate nonlinear stiffness functions are proposed. These approximate bearing models are verified by comparing their frequency responses with the results given by the exact form of nonlinearity. The proposed nonlinear dynamic model of the geared rotor-bearing system can be used to investigate the dynamic behavior and chaos.

Kahraman, Ahmet↗
NASA NTRS2000

Experiments on the Richtmyer-Meshkov Instability of Incompressible Fluids

Richtmyer-Meshkov (R-M) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and possibly become turbulent given the proper initial conditions. R-M instability is similar to the Rayleigh-Taylor (R-T) instability, which is generated when the two fluids undergo a constant acceleration. R-M instability is a fundamental fluid instability that is important to fields ranging from astrophysics to high-speed combustion. For example, R-M instability is currently the limiting factor in achieving a net positive yield with inertial confinement fusion. The experiments described here utilize a novel technique that circumvents many of the experimental difficulties previously limiting the study of the R-M instability. A Plexiglas tank contains two unequal density liquids and is gently oscillated horizontally to produce a controlled initial fluid interface shape. The tank is mounted to a sled on a high speed, low friction linear rail system, constraining the main motion to the vertical direction. The sled is released from an initial height and falls vertically until it bounces off of a movable spring, imparting an impulsive acceleration in the upward direction. As the sled travels up and down the rails, the spring retracts out of the way, allowing the instability to evolve in free-fall until impacting a shock absorber at the end of the rails. The impulsive acceleration provided to the system is measured by a piezoelectric accelerometer mounted on the tank, and a capacitive accelerometer measures the low-level drag of the bearings. Planar Laser-Induced Fluorescence is used for flow visualization, which uses an Argon ion laser to illuminate the flow and a CCD camera, mounted to the sled, to capture images of the interface. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The amplitude of the disturbance during the experiment is measured and compared to theory. The results show good agreement (within 10%) with linear stability theory up to nondimensional amplitude ka = 0.7 (wavenumber x amplitude). These results hold true for an initial ka (before acceleration) of -0.7 less than ka less than -0.06, while the linear theory was developed for absolute value of ka much less than 1. In addition, a third order weakly nonlinear perturbation theory is shown to be accurate for amplitudes as large as ka = 1.3, even though the interface becomes double-valued at ka = 1.1. As time progresses, the vorticity on the interface concentrates, and the interface spirals around the alternating sign vortex centers to form a mushroom pattern. At higher Reynolds Number (based on circulation), an instability of the vortex cores has been observed. While time limitations of the apparatus prevent determination of a critical Reynolds Number, the lowest Reynolds Number this vortex instability has been observed at is 5000.

Jacobs, J.↗
NASA NTRS1990

On the continuous dependence with respect to sampling of the linear quadratic regulator problem for distributed parameter systems

The convergence of solutions to the discrete or sampled time linear quadratic regulator problem and associated Riccati equation for infinite dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero (infinity) is established. Both the finite and infinite time horizon problems are studied. In the finite time horizon case, strong continuity of the operators which define the control system and performance index together with a stability and consistency condition on the sampling scheme are required. For the infinite time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary of delay system, and a flexible beam are presented and discussed.

Rosen, I. G.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the 'loop unrolling' technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large-scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the loop unrolling technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
NASA NTRS2019

Indoor Ground Testing of a Small UAS Sense and Avoid Airborne Doppler Radar

The National Aeronautics and Space Administration’s Unmanned Aircraft System (UAS) Traffic Management (UTM) project is researching prototype technologies needed to ensure safe integration of UAS operations into the National Airspace System (NAS). Within the UTM Concept of Operations, UAS would be equipped with on-board Sense and Avoid (SAA) technology to continually monitor for manned and unmanned aircraft in its vicinity while operating beyond visual line of sight in uncontrolled airspace. To support this effort, a candidate commercially available 24.5 GHz Doppler radar was selected and evaluated to determine if the technology could reliably support minimum requirements for SAA applications of small UAS (sUAS). Indoor ground tests were conducted inside the NASA Langley Research Center’s Experimental Test Range (ETR) from a stationary platform to evaluate the Doppler radar performance characteristics and gain operational proficiency before the radar was authorized to transmit outdoors. A high speed linear rail system was developed for the radar evaluation and was shown to be an effective method to generate Doppler radar targets of known radar cross section. The accuracy of the range and velocity reported by the radar was shown to be dependent on the Kalman filter state variance parameter settings. Antenna measurements were collected with the radar installed both on and off a sUAS to quantify the relative antenna gain, beam width and side lobe levels of the radar’s Metamaterial Electronically Scanning Array (MESA) antennas at boresight and extreme field of view pointing vectors. The relative antenna gain measured 2.6 dB lower at extreme field view angles compared to the boresight radiation pattern.

Szatkowski, George N.↗
NASA NTRS1998

The Control System for the X-33 Linear Aerospike Engine

The linear aerospike engine is being developed for single-stage -to-orbit (SSTO) applications. The primary advantages of a linear aerospike engine over a conventional bell nozzle engine include altitude compensation, which provides enhanced performance, and lower vehicle weight resulting from the integration of the engine into the vehicle structure. A feature of this integration is the ability to provide thrust vector control (TVC) by differential throttling of the engine combustion elements, rather than the more conventional approach of gimballing the entire engine. An analysis of the X-33 flight trajectories has shown that it is necessary to provide +/- 15% roll, pitch and yaw TVC authority with an optional capability of +/- 30% pitch at select times during the mission. The TVC performance requirements for X-33 engine became a major driver in the design of the engine control system. The thrust level of the X-33 engine as well as the amount of TVC are managed by a control system which consists of electronic, instrumentation, propellant valves, electro-mechanical actuators, spark igniters, and harnesses. The engine control system is responsible for the thrust control, mixture ratio control, thrust vector control, engine health monitoring, and communication to the vehicle during all operational modes of the engine (checkout, pre-start, start, main-stage, shutdown and post shutdown). The methodology for thrust vector control, the health monitoring approach which includes failure detection, isolation, and response, and the basic control system design are the topic of this paper. As an additional point of interest a brief description of the X-33 engine system will be included in this paper.

Jackson, Jerry E.↗
NASA NTRS2017

Background and Lunar Neutron Populations Detected by LEND and Average Concentration of Near-Surface Hydrogen near the Moon's Poles

Neutron flux measurements by the Lunar Exploration Neutron Detector (LEND) on the Lunar Reconnaissance Orbiter (LRO) enable quantifying hydrogen-bearing volatiles in the lunar surface from orbit. Accurately determining hydrogen abundance requires discriminating between the instrument background detection rate and the population of lunar-sourced neutrons that are sensitive to surficial hydrogen. We have investigated the detection rate for lunar and non-lunar (spacecraft-sourced) neutrons in LEND by modeling maps of measured count rate in three LEND detector systems using linear combinations of maps compiled from LEND detectors and from the Lunar Prospector Neutron Spectrometer. We find that 30% of the global-average 24.926 +/- 0.020 neutron counts per second (cps) detected by the LEND STN3 thermal-energy neutron sensor are lunar-sourced neutrons in the thermal energy range (E < 0.4 eV), 65% are lunar-sourced neutrons in the epithermal and fast energy range (E > 0.4 eV), and 5% are from spacecraft-sourced background signal. In the SETN epithermal neutron detector, 90% of the 10.622 +/- 0.002 neutron detections per second are consistent with a lunar source of epithermal and fast neutrons combined (E > 0.4 eV), with 3% due to lunar-sourced thermal neutron leakage into the detector (E < 0.4 eV), and background signal accounting for 7% of total detections. Background signal due to spacecraft-derived neutrons is substantial in the CSETN collimated detector system, accounting for 57% of the global average detection rate of 5.082 +/- 0.001 cps, greater than the 48% estimated from cruise-phase data. Lunar-sourced epithermal and fast neutrons account for 43% of detected neutrons, including neutrons in collimation as well as neutrons that penetrate the collimator wall to reach the detector. We estimate a lower limit of 17% of lunar-sourced neutrons detected by CSETN are epithermal neutrons in collimation (0.37 cps), with an upper limit estimate of 54 +/- 11% of lunar-sourced neutrons received in collimation, or 1.2 +/- 0.2 cps global average. The pole-to-equator contrast ratio in epithermal and high-energy epithermal neutron flux indicates that the average concentration of hydrogen in the polar regolith above 80deg north or south latitude is ∼105 ppmw (parts per million by weight), or 0.095 +/- 0.01 wt% water-equivalent hydrogen. Above 88deg north or south, the concentration increases to ∼140 ppmw, or 0.13 +/- 0.02 wt% water-equivalent hydrogen. The similar pattern of neutron flux suppression at both poles suggests that hydrogen concentration generally increases nearer the pole and is not closely associated with a specific feature such as Shackleton Crater at the lunar south pole that has no northern counterpart. Epithermal neutron flux decreases with increasing latitude outside the polar regions, consistent with surface hydration that increases with latitude if that hydration extends to ∼13-40 cm into the surface.

T A Livengood↗
NASA NTRS2018

Background and Lunar Neutron Populations Detected by LEND and Average Concentration of Near-Surface Hydrogen near the Moon's Poles

Neutron flux measurements by the Lunar Exploration Neutron Detector (LEND) on the Lunar Reconnaissance Orbiter (LRO) enable quantifying hydrogen-bearing volatiles in the lunar surface from orbit. Accurately determining hydrogen abundance requires discriminating between the instrument background detection rate and the population of lunar-sourced neutrons that are sensitive to surficial hydrogen. We have investigated the detection rate for lunar and non-lunar (spacecraft-sourced) neutrons in LEND by modeling maps of measured count rate in three LEND detector systems using linear combinations of maps compiled from LEND detectors and from the Lunar Prospector Neutron Spectrometer. We find that 30% of the global-average 24.926 0.020 neutron counts per second (cps) detected by the LEND STN3 thermal-energy neutron sensor are lunar-sourced neutrons in the thermal energy range (E < 0.4 eV), 65% are lunar-sourced neutrons in the epithermal and fast energy range (E > 0.4 eV), and 5% are from spacecraft-sourced background signal. In the SETN epithermal neutron detector, 90% of the 10.622 0.002 neutron detections per second are consistent with a lunar source of epithermal and fast neutrons combined (E > 0.4 eV), with 3% due to lunar-sourced thermal neutron leakage into the detector (E < 0.4 eV), and background signal accounting for 7% of total detections. Background signal due to spacecraft-derived neutrons is substantial in the CSETN collimated detector system, accounting for 57% of the global average detection rate of 5.082 0.001 cps, greater than the 48% estimated from cruise-phase data. Lunar-sourced epithermal and fast neutrons account for 43% of detected neutrons, including neutrons in collimation as well as neutrons that penetrate the collimator wall to reach the detector. We estimate a lower limit of 17% of lunar-sourced neutrons detected by CSETN are epithermal neutrons in collimation (0.37 cps), with an upper limit estimate of 54 11% of lunar-sourced neutrons received in collimation, or 1.2 0.2 cps global average. The pole-to-equator contrast ratio inepithermal and high-energy epithermal neutron flux indicates that the average concentration of hydrogen in the polar regolith above 80 north or south latitude is ~105 ppmw (parts per million by weight), or 0.095 0.01 wt% water-equivalent hydrogen. Above 88 north or south, the concentration increases to ~140 ppmw, or 0.13 0.02 wt% water-equivalent hydrogen. The similar pattern of neutron flux suppression at both poles suggests that hydrogen concentration generally increases nearer the pole and is not closely associated with a specific feature such as Shackleton Crater at the lunar south pole that has no northern counterpart. Epithermal neutron flux decreases with increasing latitude outside the polar regions, consistent with surface hydration that increases with latitude if that hydration extends to ~13-40 cm into the surface.

Livengood, Timothy A.↗
NASA NTRSDate not supplied

Verifying PLC Programs via Monitors: Extending the Integration of FRET and PLCverif

Verification of Programmable Logic Controller (PLC) programs requires reasoning about propositions qualified in terms of time. CERN’s PLCverif, an open-source tool for the analysis of safety-critical PLC systems, uses Linear Temporal Logic (LTL) for the specification of properties. Until now, PLCverif depended on third-party tools that accept LTL specifications to perform verification. However, our experience with industrial PLC programs shows that, to overcome analysis limitations, a wide range of techniques are needed to successfully verify complex properties. In this paper, we extend PLCverif to enable PLC program verification of pure-past LTL (PLTL) safety properties with assertion-based verification tools. To this end, we take an algorithm from the runtime-monitoring domain, apply it to bounded model checking of PLC programs, and implement it in PLCverif. We extend the integration of NASA’s Formal Requirements Elicitation Tool (FRET) into PLCverif to use PLTL properties generated with FRET. In addition, we leverage the program structure induced by the PLC scan-cycle for a state-space reduction. Finally, we expose the algorithm to a real-world case study of critical systems at CERN.

Formal verification↗