Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential mode”

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 541 records · Page 30

Particle Size Distributions Measured in the B757 Engine Plume During EXCAVATE

The Experiment to Characterize Aircraft Volatile Aerosols and Trace Species Emissions (EXCAVATE) took place at NASA Langley Research Center during January 2002. This ground based study was conducted to examine the role of fuel sulfur content on particulate emissions. Size distributions as a function of engine operating conditions were measured in the exhaust plume of a B-757 at four downstream axial locations (1 m, 10 m, 25 m and 35 m). The engine was run on JP-5 with three different sulfur concentrations, 810 ppm, 1050 ppm, 1820 ppm; and was operated over a range of power settings from idle to near-full power. Zalabsky differential-mobility analyzers DMAS), Met One condensation-nuclei counters (CNCs), and a TSI 3022 condensation-particle counter (CPC) were used to measure the size distributions. The total number-count (particle concentration), number-based Emissions Index (EInumber) and mass-based Emissions Index (E1-J increased with fuel sulfur-content and engine pressure ratio (EPR). Count Mean Diameter (Ch4D) also increased with EPR yet remained fairly constant with fuel sulfur-content for a fixed location in the exhaust plume. Also the mode and CMD both increased with distance in the plume.

Sanders, Terry↗

Four-Quadrant Analog Multipliers Using G4-FETs

Theoretical analysis and some experiments have shown that the silicon-on-insulator (SOI) 4-gate transistors known as G4-FETs can be used as building blocks of four-quadrant analog voltage multiplier circuits. Whereas a typical prior analog voltage multiplier contains between six and 10 transistors, it is possible to construct a superior voltage multiplier using only four G4-FETs. A G4-FET is a combination of a junction field-effect transistor (JFET) and a metal oxide/semiconductor field-effect transistor (MOSFET). It can be regarded as a single transistor having four gates, which are parts of a structure that affords high functionality by enabling the utilization of independently biased multiple inputs. The structure of a G4-FET of the type of interest here (see Figure 1) is that of a partially-depleted SOI MOSFET with two independent body contacts, one on each side of the channel. The drain current comprises of majority charge carriers flowing from one body contact to the other that is, what would otherwise be the side body contacts of the SOI MOSFET are used here as the end contacts [the drain (D) and the source (S)] of the G4-FET. What would otherwise be the source and drain of the SOI MOSFET serve, in the G4-FET, as two junction-based extra gates (JG1 and JG2), which are used to squeeze the channel via reverse-biased junctions as in a JFET. The G4-FET also includes a polysilicon top gate (G1), which plays the same role as does the gate in an accumulation-mode MOSFET. The substrate emulates a fourth MOS gate (G2). By making proper choices of G4-FET device parameters in conjunction with bias voltages and currents, one can design a circuit in which two input gate voltages (Vin1,Vin2) control the conduction characteristics of G4-FETs such that the output voltage (Vout) closely approximates a value proportional to the product of the input voltages. Figure 2 depicts two such analog multiplier circuits. In each circuit, there is the following: The input and output voltages are differential, The multiplier core consists of four G4- FETs (M1 through M4) biased by a constant current sink (Ibias), and The G4-FETs in two pairs are loaded by two identical resistors (RL), which convert a differential output current to a differential output voltage. The difference between the two circuits stems from their input and bias configurations. In each case, provided that the input voltages remain within their design ranges as determined by considerations of bias, saturation, and cutoff, then the output voltage is nominally given by Vout = kVin1Vin2, where k is a constant gain factor that depends on the design parameters and is different for the two circuits. In experimental versions of these circuits constructed using discrete G4- FETs and resistors, multiplication of voltages in all four quadrants (that is, in all four combinations of input polarities) was demonstrated, and deviations of the output voltages from linear dependence on the input voltages were found to amount to no more than a few percent. It is anticipated that in fully integrated versions of these circuits, the deviations from linearity will be made considerably smaller through better matching of devices.

Mojarradi, Mohammad↗

Was core formation violent enough to homogenize the early mantle?

The dynamics of iron, its thermal state and its phase in the accreting Earth probably played a major role in the Earth's early thermal evolution. Plausible impact thermal histories make it possible that pure iron was molten in the accreting Earth after it was about 10% grown. Hence, iron eutectic alloys (FeS, FeO) certainly were. Additionally, the initial temperature of the core is an important constraint on the secular cooling of the early Earth and on the strength of the early geodynamo. Whether iron is solid or molten would influence geochemical equilibria in the upper and lower mantle; the mode of core formation, by spherical or near-spherical blobs, stalk-like instabilities, or something more catastrophic would influence the partitioning of siderophiles between silicate and iron phases. Early descent of iron (during accretion) favors partitioning according to low-pressure phase equilibria, whereas late descent favors higher pressure. The later core formation occurs, the greater the heat pulse, due to the strong dependence of gravitational potential energy on planetary radius. The heat may homogenize the mantle if core formation is global; otherwise, heterogeneity of iron differentiation may leave some of the pre-archean mantle unaffected. The larger the chunks of proto-core (and hence smaller surface/volume ratios) the greater the heterogeneity.

Cooperman, S. A.↗

Optimal Force Control of Vibro-Impact Systems for Autonomous Drilling Applications

The need to maintain optimal energy efficiency is critical during the drilling operations performed on future and current planetary rover missions (see figure). Specifically, this innovation seeks to solve the following problem. Given a spring-loaded percussive drill driven by a voice-coil motor, one needs to determine the optimal input voltage waveform (periodic function) and the optimal hammering period that minimizes the dissipated energy, while ensuring that the hammer-to-rock impacts are made with sufficient (user-defined) impact velocity (or impact energy). To solve this problem, it was first observed that when voice-coil-actuated percussive drills are driven at high power, it is of paramount importance to ensure that the electrical current of the device remains in phase with the velocity of the hammer. Otherwise, negative work is performed and the drill experiences a loss of performance (i.e., reduced impact energy) and an increase in Joule heating (i.e., reduction in energy efficiency). This observation has motivated many drilling products to incorporate the standard bang-bang control approach for driving their percussive drills. However, the bang-bang control approach is significantly less efficient than the optimal energy-efficient control approach solved herein. To obtain this solution, the standard tools of classical optimal control theory were applied. It is worth noting that these tools inherently require the solution of a two-point boundary value problem (TPBVP), i.e., a system of differential equations where half the equations have unknown boundary conditions. Typically, the TPBVP is impossible to solve analytically for high-dimensional dynamic systems. However, for the case of the spring-loaded vibro-impactor, this approach yields the exact optimal control solution as the sum of four analytic functions whose coefficients are determined using a simple, easy-to-implement algorithm. Once the optimal control waveform is determined, it can be used optimally in the context of both open-loop and closed-loop control modes (using standard realtime control hardware).

Aldrich, Jack B.↗

The nonlinear evolution of inviscid Goertler vortices in three-dimensional boundary layers

The nonlinear development of inviscid Gortler vortices in a three-dimensional boundary layer is considered. We do not follow the classical approach of weakly nonlinear stability problems and consider a mode which has just become unstable. Instead we extend the method of Blackaby, Dando, and Hall (1992), which considered the closely related nonlinear development of disturbances in stratified shear flows. The Gortler modes we consider are initially fast growing and we assume, following others, that boundary-layer spreading results in them evolving in a linear fashion until they reach a stage where their amplitudes are large enough and their growth rates have diminished sufficiently so that amplitude equations can be derived using weakly nonlinear and non-equilibrium critical-layer theories. From the work of Blackaby, Dando and Hall (1993) is apparent, given the range of parameters for the Gortler problem, that there are three possible nonlinear integro-differential evolution equations for the disturbance amplitude. These are a cubic due to viscous effects, a cubic which corresponds to the novel mechanism investigated in this previous paper, and a quintic. In this paper we shall concentrate on the two cubic integro-differential equations and in particular, on the one due to the novel mechanism as this will be the first to affect a disturbance. It is found that the consideration of a spatial evolution problem as opposed to temporal (as was considered in Blackaby, Dando, and Hall, 1992) causes a number of significant changes to the evolution equations.

Blackaby, Nicholas↗

The role of nonlinear critical layers in boundary layer transition

Asymptotic methods are used to describe the nonlinear self-interaction between pairs of oblique instability modes that eventually develops when initially linear spatially growing instability waves evolve downstream in nominally two-dimensional laminar boundary layers. The first nonlinear reaction takes place locally within a so-called 'critical layer', with the flow outside this layer consisting of a locally parallel mean flow plus a pair of oblique instability waves - which may or may not be accompanied by an associated plane wave. The amplitudes of these waves, which are completely determined by nonlinear effects within the critical layer, satisfy either a single integro-differential equation or a pair of integro-differential equations with quadratic to quartic-type nonlinearities. The physical implications of these equations are discussed.

Goldstein, M.E.↗

Recent advances in capacitance type of blade tip clearance measurements

Two recent electronic advances at NASA-Lewis that meet the blade tip clearance needs of a wide class of fans, compressors, and turbines are described. The first is a frequency modulated (FM) oscillator that requires only a single low cost ultrahigh frequency operational amplifier. Its carrier frequency is 42.8 MHz when used with a 61 cm long hermetically sealed coaxial cable. The oscillator can be calibrated in the static mode and has a negative peak frequency deviation of 400 kHz for a typical rotor blade. High temperature performance tests of the probe and 13 cm of the adjacent cable show good accuracy up to 600 C, the maximum which produces a clearance error of + or - 10 microns at a clearance of 500 microns. In the second advance, a guarded probe configuration allows a longer cable capacitance. The capacitance of the probe is part of a small time constant feedback in a high speed operational amplifier. The solution of the governing differential equation is applied to a ramp type of input. The results show an amplifier output that contains a term which is proportional to the derivative of the feedback capacitance. The capacitance is obtained by subtracting a balancing reference channel followed by an integration stage.

Barranger, John P.↗

Recent advances in capacitance type of blade tip clearance measurements

Two recent electronic advances at NASA-Lewis that meet the blade tip clearance needs of a wide class of fans, compressors, and turbines are described. The first is a frequency modulated (FM) oscillator that requires only a single low cost ultrahigh frequency operational amplifier. Its carrier frequency is 42.8 MHz when used with a 61 cm long hermetically sealed coaxial cable. The oscillator can be calibrated in the static mode and has a negative peak frequency deviation of 400 kHz for a typical rotor blade. High temperature performance tests of the probe and 13 cm of the adjacent cable show good accuracy up to 600 C, the maximum which produces a clearance error of + or - 10 microns at a clearance of 500 microns. In the second advance, a guarded probe configuration allows a longer cable capacitance. The capacitance of the probe is part of a small time constant feedback in a high speed operational amplifier. The solution of the governing differential equation is applied to a ramp type of input. The results show an amplifier output that contains a term which is proportional to the derivative of the feedback capacitance. The capacitance is obtained by subtracting a balancing reference channel followed by an integration stage.

Barranger, John P.↗

The Geophysical Fluid Flow Cell Experiment

The Geophysical Fluid Flow Cell (GFFC) experiment performed visualizations of thermal convection in a rotating differentially heated spherical shell of fluid. In these experiments dielectric polarization forces are used to generate a radially directed buoyancy force. This enables the laboratory simulation of a number of geophysically and astrophysically important situations in which sphericity and rotation both impose strong constraints on global scale fluid motions. During USML-2 a large set of experiments with spherically symmetric heating were carried out. These enabled the determination of critical points for the transition to various forms of nonaxisymmetric convection and, for highly turbulent flows, the transition latitudes separating the different modes of motion. This paper presents a first analysis of these experiments as well as data on the general performance of the instrument during the USML-2 flight.

Hart, J. E.↗

NASA Tech Briefs, July 2013

Dielectrophoresis-Based Particle Sensor Using Nanoelectrode Arrays; Multi-Dimensional Damage Detection for Surfaces and Structures; ULTRA: Underwater Localization for Transit and Reconnaissance Autonomy; Autonomous Cryogenic Leak Detector for Improving Launch Site Operations; Submillimeter Planetary Atmospheric Chemistry Exploration Sounder; Method for Reduction of Silver Biocide Plating on Metal Surfaces; Silicon Micromachined Microlens Array for THz Antennas; Forward-Looking IED Detector Ground Penetrating Radar; Fully Printed, Flexible, Phased Array Antenna for Lunar Surface Communication, Battery Charge Equalizer with Transformer Array; An Efficient, Highly Flexible Multi-Channel Digital Downconverter Architecture; Dimmable Electronic Ballast for a Gas Discharge Lamp; Conductive Carbon Nanotube Inks for Use with Desktop Inkjet Printing Technology; Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates; High-Performance, Low-Temperature-Operating, Long-Lifetime Aerospace Lubricants; Carbon Nanotube Microarrays Grown on Nanoflake Substrates; Differential Muon Tomography to Continuously Monitor Changes in the Composition of Subsurface Fluids; Microgravity Drill and Anchor System; 20 Granular Media-Based Tunable Passive Vibration Suppressor; 21 Miga Aero Actuator and 2D Machined Mechanical Binary Latch; Micro-XRF for In Situ Geological Exploration of Other Planets; Hydrogen-Enhanced Lunar Oxygen Extraction and Storage Using Only Solar Power; Uplift of Ionospheric Oxygen Ions During Extreme Magnetic Storms; Miniaturized, High-Speed, Modulated X-Ray Source; Hollow-Fiber Spacesuit Water Membrane Evaporator 25 High-Power Single-Mode 2.65-micrometers InGaAsSb/AlInGaAsSb Diode Lasers; Optical Device for Converting a Laser Beam Into Two Co-aligned but Oppositely Directed Beams; A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging; X-Ray Diffractive Optics; SynGenics Optimization System (SynOptSys); 29 CFD Script for Rapid TPS Damage Assessment; radEq Add-On Module for CFD Solver Loci-CHEM; Science Opportunity Analyzer (SOA) Version 8; 30 Autonomous Byte Stream Randomizer; Distributed Engine Control Empirical/Analytical Verification Tools; Dynamic Server-Based KML Code Generator Method for Level-of-Detail Traversal of Geospatial Data; Automated Planning of Science Products Based on Nadir Overflights and Alerts for Onboard and Ground Processing; Linked Autonomous Interplanetary Satellite Orbit Navigation; Risk-Constrained Dynamic Programming for Optimal Mars Entry, Descent, and Landing; Scheduling Operations for Massive Heterogeneous Clusters; Deepak Condenser Model (DeCoM); Flight Software Math Library; Recirculating 1-K-Pot for Pulse-Tube Cryostats; 35 Method for Processing Lunar Regolith Using Microwaves; Wells for In Situ Extraction of Volatiles from Regolith (WIEVR); and Estimating the Backup Reaction Wheel Orientation Using Reaction Wheel Spin Rates Flight Telemetry from a Spacecraft.

Source record↗

The effect of an unsteady three-dimensional wake on elastic blade-flapping eigenvalues in hover

This paper describes the formulation of a finite-state inflow model based on an acceleration potential and a helical wake geometry. The states of the model are coefficients of an inflow expansion in terms of a Fourier series (azimuthally) and of special polynomials (radially). The integrals over the wake are done in closed-form to obtain a set of ordinary differential equations for the inflow coefficients. The forcing functions for these equations are generalized forces which are integrals of the blade loading exactly as in structural dynamics. This model implicitly includes (for the hover case) Prandtl-Goldstein tip losses, dynamic inflow, and Theodorsen/Loewy lift deficiency. Thus, it is a fully three-dimensional unsteady wake model. Here, this model is coupled with elastic-blade equations in hover and eigenvalues are found. The results show that the three-dimensional wake has a large effect on the flap damping of all modes.

Peters, David A.↗

Effect of chordwise forces and deformations and deformations due to steady lift on wing flutter

This investigation explores the effects of chordwise forces and deformations and steady-state deformation due to lift on the static and dynamic aeroelastic stability of a uniform cantilever wing. Results of this analysis are believed to have practical applications for high-performance sailplanes and certain RPV's. The airfoil cross section is assumed to be symmetric and camber bending is neglected. Motions in vertical bending, fore-and-aft bending, and torsion are considered. A differential equation model is developed, which included the nonlinear elastic bending-torsion coupling that accompanies even moderate deflections. A linearized expansion in small time-dependent deflections is made about a steady flight condition. The stability determinant of the linearized system then contains coefficients that depend on steady displacements. Loads derived from two-dimensional incompressible aerodynamic theory are used to obtain the majority of the results, but cases using three-dimensional subsonic compressible theory are also studied. The stability analysis is carried out in terms of the dynamically uncoupled natural modes of vibration of the uniform cantilever.

Boyd, W. N.↗

MAP Attitude Control System Design and Analysis

The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The MAP spacecraft will perform its mission in a Lissajous orbit around the Earth-Sun L(sub 2) Lagrange point to suppress potential instrument disturbances. To make a full-sky map of cosmic microwave background fluctuations, a combination fast spin and slow precession motion will be used. MAP requires a propulsion system to reach L(sub 2), to unload system momentum, and to perform stationkeeping maneuvers once at L(sub 2). A minimum hardware, power and thermal safe control mode must also be provided. Sufficient attitude knowledge must be provided to yield instrument pointing to a standard deviation of 1.8 arc-minutes. The short development time and tight budgets require a new way of designing, simulating, and analyzing the Attitude Control System (ACS). This paper presents the design and analysis of the control system to meet these requirements.

Andrews, S. F.↗

Experiment 7: The Geophysical Fluid Flow Cell Experiment on USML-2

The Geophysical Fluid Flow Cell (GFFC) experiment performed visualizations of thermal convection in a rotating differentially heated spherical shell of fluid. In these experiments dielectric polarization forces are used to generate a radially directed buoyancy force. This enables the laboratory simulation of a number of geophysically and astrophysically important situations in which sphericity and rotation both impose strong constraints on global scale fluid motions. During USML-2 a large set of experiments with spherically symmetric heating were carried out. These enabled the determination of critical points for the transition to various forms of non-axisymmetric convection and, for highly turbulent flows, the transition latitudes separating the different modes of motion. This paper presents a first analysis of these experiments as well as data on the general performance of the instrument during the USML-2 flight.

Hart, John E.↗

Broadband Noise of Fans - With Unsteady Coupling Theory to Account for Rotor and Stator Reflection/Transmission Effects

This report examines the effects on broadband noise generation of unsteady coupling between a rotor and stator in the fan stage of a turbofan engine. Whereas previous acoustic analyses treated the blade rows as isolated cascades, the present work accounts for reflection and transmission effects at both blade rows by tracking the mode and frequency scattering of pressure and vortical waves. The fan stage is modeled in rectilinear geometry to take advantage of a previously existing unsteady cascade theory for 3D perturbation waves and thereby use a realistic 3D turbulence spectrum. In the analysis, it was found that the set of participating modes divides itself naturally into "independent mode subsets" that couple only among themselves and not to the other such subsets. This principle is the basis for the analysis and considerably reduces computational effort. It also provides a simple, accurate scheme for modal averaging for further efficiency. Computed results for a coupled fan stage are compared with calculations for isolated blade rows. It is found that coupling increases downstream noise by 2 to 4 dB. Upstream noise is lower for isolated cascades and is further reduced by including coupling effects. In comparison with test data, the increase in the upstream/downstream differential indicates that broadband noise from turbulent inflow at the stator dominates downstream noise but is not a significant contributor to upstream noise.

Hanson, Donald B.↗

Separation and Transition on the ROTEX-T Cone-Flare

As part of NATO STO AVT-346 “Predicting Hypersonic Boundary-Layer Transition on Complex Geometries,” coordinated experimental and computational studies were conducted on the ROTEX-T, a cone-flare geometry used in a successful flight-test experiment. At the as-flown conditions, a separation bubble existed at the compression corner. Separation, reattachment, and the multifaceted linear instability paths leading this bubble to transition to turbulence are challenging to predict, but have significant impact on surface pressure and heat flux. High-resolution background-oriented schlieren and infrared thermography measurements were made in the AFOSR–Notre Dame Large Mach-6 Quiet Tunnel at freestream unit Reynolds numbers from 5.8×10(exp 6) to 12.2×10(exp 6) m-1and nominally zero angle of attack. High-speed self-aligned focusing schlieren, infrared thermography, and focused laser differential interferometry measurements were made in the AFRL Mach-6 Ludwieg Tube from 2.2×10(exp 6) to 24.7×10(exp 6) m-1and nominally zero degrees angle of attack. The surface heat-flux and Stanton number distributions were computed. Separation and reattachment locations, as well as the flow state at each, were determined from the combination of surface and off-wall measurements. The convective and global boundary-layer instabilities of the axisymmetric laminar flow at the experimental conditions were investigated computationally. Amplification of Mack’s first and second modes were observed to have logarithmic amplification factors between 5 to 7.5 at the separation location, depending on conditions. The flow was found to be globally unstable to stationary three-dimensional disturbances concentrated in the reattachment region. Previous analysis of the ROTEX-T flight data had not assessed reattachment location or the flow state upon reattachment. Thanks to the insights gained from the coordinated, on- and off-wall ground-test measurements, these evaluations have now been made. The separation location indicated by laminar simulations is consistently numerically predicted to be upstream of the experimentally observed location for a transitional separation bubble. The cause of this difference is understood to lie in the steady-state and axisymmetric assumptions made by both solvers employed to compute the basic states analyzed, as flow topology considerations assert that unsteady two-dimensional or axisymmetric separation bubbles are structurally unstable and will become three-dimensional. Computed laminar heating rates prior to separation agreed well with experiment; transitional heating rates after reattachment were between laminar and turbulent computations.

CFD Validation↗

Separation and Transition on the ROTEX-T Cone-Flare

As part of NATO STO AVT-346 “Predicting Hypersonic Boundary-Layer Transition on Complex Geometries,” coordinated experimental and computational studies were conducted on the ROTEX-T, a cone-flare geometry used in a successful flight-test experiment. At the as-flown conditions, a separation bubble existed at the compression corner. Separation, reattachment, and the multifaceted linear instability paths leading this bubble to transition to turbulence are challenging to predict, but have significant impact on surface pressure and heat flux. High-resolution background-oriented schlieren and infrared thermography measurements were made in the AFOSR–Notre Dame Large Mach-6 Quiet Tunnel at freestream unit Reynolds numbers from 5.8×10(exp 6) to 12.2×10(exp 6) m-1and nominally zero angle of attack. High-speed self-aligned focusing schlieren, infrared thermography, and focused laser differential interferometry measurements were made in the AFRL Mach-6 Ludwieg Tube from 2.2×10(exp 6) to 24.7×10(exp 6) m-1and nominally zero degrees angle of attack. The surface heat-flux and Stanton number distributions were computed. Separation and reattachment locations, as well as the flow state at each, were determined from the combination of surface and off-wall measurements. The convective and global boundary-layer instabilities of the axisymmetric laminar flow at the experimental conditions were investigated computationally. Amplification of Mack’s first and second modes were observed to have logarithmic amplification factors between 5 to 7.5 at the separation location, depending on conditions. The flow was found to be globally unstable to stationary three-dimensional disturbances concentrated in the reattachment region. Previous analysis of the ROTEX-T flight data had not assessed reattachment location or the flow state upon reattachment. Thanks to the insights gained from the coordinated, on- and off-wall ground-test measurements, these evaluations have now been made. The separation location indicated by laminar simulations is consistently numerically predicted to be upstream of the experimentally observed location for a transitional separation bubble. The cause of this difference is understood to lie in the steady-state and axisymmetric assumptions made by both solvers employed to compute the basic states analyzed, as flow topology considerations assert that unsteady two-dimensional or axisymmetric separation bubbles are structurally unstable and will become three-dimensional. Computed laminar heating rates prior to separation agreed well with experiment; transitional heating rates after reattachment were between laminar and turbulent computations.

CFD Validation↗

Determination of eigenvalues of dynamical systems by symbolic computation

A symbolic computation technique for determining the eigenvalues of dynamical systems is described wherein algebraic operations, symbolic differentiation, matrix formulation and inversion, etc., can be performed on a digital computer equipped with a formula-manipulation compiler. An example is included that demonstrates the facility with which the system dynamics matrix and the control distribution matrix from the state space formulation of the equations of motion can be processed to obtain eigenvalue loci as a function of a system parameter. The example chosen to demonstrate the technique is a fourth-order system representing the longitudinal response of a DC 8 aircraft to elevator inputs. This simplified system has two dominant modes, one of which is lightly damped and the other well damped. The loci may be used to determine the value of the controlling parameter that satisfied design requirements. The results were obtained using the MACSYMA symbolic manipulation system.

Howard, J. C.↗