Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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,531 records · Page 85

Performance of the NASA Airborne Radar with the Windshear Database for Forward-Looking Systems

This document describes the simulation approach used to test the performance of the NASA airborne windshear radar. An explanation of the actual radar hardware and processing algorithms provides an understanding of the parameters used in the simulation program. This report also contains a brief overview of the NASA airborne windshear radar experimental flight test results. A description of the radar simulation program shows the capabilities of the program and the techniques used for certification evaluation. Simulation of the NASA radar is comprised of three steps. First, the choice of the ground clutter data must be made. The ground clutter is the return from objects in or nearby an airport facility. The choice of the ground clutter also dictates the aircraft flight path since ground clutter is gathered while in flight. The second step is the choice of the radar parameters and the running of the simulation program which properly combines the ground clutter data with simulated windshear weather data. The simulated windshear weather data is comprised of a number of Terminal Area Simulation System (TASS) model results. The final step is the comparison of the radar simulation results to the known windshear data base. The final evaluation of the radar simulation is based on the ability to detect hazardous windshear with the aircraft at a safe distance while at the same time not displaying false alerts.

Switzer, George F.↗

3D Structured Grid Adaptation

Grid adaptation for structured meshes is the art of using information from an existing, but poorly resolved, solution to automatically redistribute the grid points in such a way as to improve the resolution in regions of high error, and thus the quality of the solution. This involves: (1) generate a grid vis some standard algorithm, (2) calculate a solution on this grid, (3) adapt the grid to this solution, (4) recalculate the solution on this adapted grid, and (5) repeat steps 3 and 4 to satisfaction. Steps 3 and 4 can be repeated until some 'optimal' grid is converged to but typically this is not worth the effort and just two or three repeat calculations are necessary. They also may be repeated every 5-10 time steps for unsteady calculations.

Banks, D. W.↗

Dynamics of Superfluid Helium in Low-Gravity

This report summarizes the work performed under a contract entitled 'Dynamics of Superfluid Helium in Low Gravity'. This project performed verification tests, over a wide range of accelerations of two Computational Fluid Dynamics (CFD) codes of which one incorporates the two-fluid model of superfluid helium (SFHe). Helium was first liquefied in 1908 and not until the 1930s were the properties of helium below 2.2 K observed sufficiently to realize that it did not obey the ordinary physical laws of physics as applied to ordinary liquids. The term superfluidity became associated with these unique observations. The low temperature of SFHe and it's temperature unifonrmity have made it a significant cryogenic coolant for use in space applications in astronomical observations with infrared sensors and in low temperature physics. Superfluid helium has been used in instruments such as the Shuttle Infrared Astronomy Telescope (IRT), the Infrared Astronomy Satellite (IRAS), the Cosmic Background Observatory (COBE), and the Infrared Satellite Observatory (ISO). It is also used in the Space Infrared Telescope (SIRTF), Relativity Mission Satellite formally called Gravity Probe-B (GP-B), and the Test of the Equivalence Principle (STEP) presently under development. For GP-B and STEP, the use of SFHE is used to cool Superconducting Quantum Interference Detectors (SQUIDS) among other parts of the instruments. The Superfluid Helium On-Orbit Transfer (SHOOT) experiment flown in the Shuttle studied the behavior of SFHE. This experiment attempted to get low-gravity slosh data, however, the main emphasis was to study the low-gravity transfer of SFHE from tank to tank. These instruments carried tanks of SFHE of a few hundred liters to 2500 liters. The capability of modeling the behavior of SFHE is important to spacecraft control engineers who must design systems that can overcome disturbances created by the movement of the fluid. In addition instruments such as GP-B and STEP are very sensitive to quasi-steady changes in the mass distribution of the liquid. The CFD codes were used to model the fluid's dynamic motion. Tests in one-g were performed with the main emphasis on being able to compute the actual damping of the fluid. A series of flights on the NASA Lewis reduced gravity DC-9 aircraft were performed with the Jet Propulsion Laboratory (JPL) Low Temperature Flight Facility and a superfluid Test Cell. The data at approximately 0.04g, lg and 2g were used to determine if correct fundamental frequencies can be predicted based on the acceleration field. Tests in zero gravity were performed to evaluate zero gravity motion.

Frank, David J.↗

Insoluble Coatings for Stirling Engine Heat Pipe Condenser Surfaces

The principal objective of this Phase 2 SBIR program was to develop and demonstrate a practically insoluble coating for nickel-based superalloys for Stirling engine heat pipe applications. Specific technical objectives of the program were: (1) Determine the solubility corrosion rates for Nickel 200, Inconel 718, and Udimet 72OLI in a simulated Stirling engine heat pipe environment, (2) Develop coating processes and techniques for capillary groove and screen wick structures, (3) Evaluate the durability and solubility corrosion rates for capillary groove and screen wick structures coated with an insoluble coating in cylindrical heat pipes operating under Stirling engine conditions, and (4) Design and fabricate a coated full-scale, partial segment of the current Stirling engine heat pipe for the Stirling Space Power Convertor program. The work effort successfully demonstrated a two-step nickel aluminide coating process for groove wick structures and interior wall surfaces in contact with liquid metals; demonstrated a one-step nickel aluminide coating process for nickel screen wick structures; and developed and demonstrated a two-step aluminum-to-nickel aluminide coating process for nickel screen wick structures. In addition, the full-scale, partial segment was fabricated and the interior surfaces and wick structures were coated. The heat pipe was charged with sodium, processed, and scheduled to be life tested for up to ten years as a Phase 3 effort.

Dussinger, Peter M.↗

Aerobic Fitness Does Not Contribute to Prediction of Orthostatic Intolerance

Several investigations have suggested that orthostatic tolerance may be inversely related to aerobic fitness (VO (sub 2max)). To test this hypothesis, 18 males (age 29 to 51 yr) underwent both treadmill VO(sub 2max) determination and graded lower body negative pressures (LBNP) exposure to tolerance. VO(2max) was measured during the last minute of a Bruce treadmill protocol. LBNP was terminated based on pre-syncopal symptoms and LBNP tolerance (peak LBNP) was expressed as the cumulative product of LBNP and time (torr-min). Changes in heart rate, stroke volume cardiac output, blood pressure and impedance rheographic indices of mid-thigh-leg initial accumulation were measured at rest and during the final minute of LBNP. For all 18 subjects, mean (plus or minus SE) fluid accumulation index and leg venous compliance index at peak LBNP were 139 plus or minus 3.9 plus or minus 0.4 ml-torr-min(exp -2) x 10(exp 3), respectively. Pearson product-moment correlations and step-wise linear regression were used to investigate relationships with peak LBNP. Variables associated with endurance training, such as VO(sub 2max) and percent body fat were not found to correlate significantly (P is less than 0.05) with peak LBNP and did not add sufficiently to the prediction of peak LBNP to be included in the step-wise regression model. The step-wise regression model included only fluid accumulation index leg venous compliance index, and blood volume and resulted in a squared multiple correlation coefficient of 0.978. These data do not support the hypothesis that orthostatic tolerance as measured by LBNP is lower in individuals with high aerobic fitness.

Convertino, Victor A.↗

Hybrid Analytical Technique for Nonlinear Vibration Analysis of Thin-Walled Beams

A two-step hybrid analytical technique is presented for the nonlinear vibration analysis of thin-walled beams. The first step involves the generation of various-order perturbation functions using the Linstedt-Poincare perturbation technique. The second step consists of using the perturbation functions as coordinate (or approximation) functions and then computing both the amplitudes of these functions and the nonlinear frequency of vibration via a direct variational procedure. The analytical formulation is based on a form of the geometrically nonlinear beam theory with the effects of in-plane inertia, rotatory inertia, and transverse shear deformation included. The effectiveness of the proposed technique is demonstrated by means of a numerical example of thin-walled beam with a doubly symmetric I-section. The solutions obtained using a single-spatial mode were compared with those obtained using multiple-spatial modes. The standard of comparison was taken to be the frequencies obtained by the direct integration/fast Fourier transform (FFT) technique. The nonlinear frequencies obtained by the hybrid technique were shown to converge to the corresponding ones obtained by the direct integration/fast Fourier transform (FFT) technique well beyond the range of applicability of the perturbation technique. The frequencies and total strain energy of the beam were overestimated by using a single-spatial mode.

Noor, Ahmed K.↗

A Geometry Based Infra-structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis (CFD, structural analysis, and etc.) are: Surface Generation - usually by employing a CAD system; Grid Generation - preparing the volume for the simulation; Flow Solver - producing the results at the specified operational point; and Post-processing Visualization - interactively attempting to understand the results For structural analysis, integrated systems can be obtained from a number of commercial vendors. For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. Specifically the problems with this procedure are: (1) File based. Information flows from one step to the next via data files with formats specified for that procedure. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. (3) One-Way communication. All information travels on from one phase to the next. Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive.

Haimes, Robert↗

Physiological Response of Plants Grown on Porous Ceramic Tubes

This research involves the manipulation of the root-zone water potential for the purposes of discriminating the rate limiting step in the inorganic nutrient uptake mechanism utilized by higher plants. This reaction sequence includes the pathways controlled by the root-zone conditions such as water tension and gradient concentrations. Furthermore, plant based control mechanisms dictated by various protein productions are differentiated as well. For the nutrients limited by the environmental availability, the kinetics were modeled using convection and diffusion equations. Alternatively, for the nutrients dependent upon enzyme manipulations, the uptakes are modeled using Michaelis-Menten kinetics. In order to differentiate between these various mechanistic steps, an experimental apparatus known as the Porous Ceramic Tube - Nutrient Delivery System (PCT-NDS) was used. Manipulation of the applied suction pressure circulating a nutrient solution through this system imposes a change in the matric component of the water potential. This compensates for the different osmotic components of water potential dictated by nutrient concentration. By maintaining this control over the root-zone conditions, the rate limiting steps in the uptake of the essential nutrients into tomato plants (Lycopersicon esculentum cv. Cherry Elite) were differentiated. Results showed that the uptake of some nutrients were mass transfer limited while others were limited by the enzyme kinetics. Each of these were adequately modeled with calculations and discussions of the parameter estimations provided.

Tsao, David↗

Multi Objective Controller Design for Linear System via Optimal Interpolation

We propose a methodology for the design of a controller which satisfies a set of closed-loop objectives simultaneously. The set of objectives consists of: (1) pole placement, (2) decoupled command tracking of step inputs at steady-state, and (3) minimization of step response transients with respect to envelope specifications. We first obtain a characterization of all controllers placing the closed-loop poles in a prescribed region of the complex plane. In this characterization, the free parameter matrix Q(s) is to be determined to attain objectives (2) and (3). Objective (2) is expressed as determining a Pareto optimal solution to a vector valued optimization problem. The solution of this problem is obtained by transforming it to a scalar convex optimization problem. This solution determines Q(O) and the remaining freedom in choosing Q(s) is used to satisfy objective (3). We write Q(s) = (l/v(s))bar-Q(s) for a prescribed polynomial v(s). Bar-Q(s) is a polynomial matrix which is arbitrary except that Q(O) and the order of bar-Q(s) are fixed. Obeying these constraints bar-Q(s) is now to be 'shaped' to minimize the step response characteristics of specific input/output pairs according to the maximum envelope violations. This problem is expressed as a vector valued optimization problem using the concept of Pareto optimality. We then investigate a scalar optimization problem associated with this vector valued problem and show that it is convex. The organization of the report is as follows. The next section includes some definitions and preliminary lemmas. We then give the problem statement which is followed by a section including a detailed development of the design procedure. We then consider an aircraft control example. The last section gives some concluding remarks. The Appendix includes the proofs of technical lemmas, printouts of computer programs, and figures.

Ozbay, Hitay↗

SeaWiFS Science Algorithm Flow Chart

This flow chart describes the baseline science algorithms for the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) Data Processing System (SDPS). As such, it includes only processing steps used in the generation of the operational products that are archived by NASA's Goddard Space Flight Center (GSFC) Distributed Active Archive Center (DAAC). It is meant to provide the reader with a basic understanding of the scientific algorithm steps applied to SeaWiFS data. It does not include non-science steps, such as format conversions, and places the greatest emphasis on the geophysical calculations of the level-2 processing. Finally, the flow chart reflects the logic sequences and the conditional tests of the software so that it may be used to evaluate the fidelity of the implementation of the scientific algorithm. In many cases however, the chart may deviate from the details of the software implementation so as to simplify the presentation.

Darzi, Michael↗

Fully Threaded Tree for Adaptive Refinement Fluid Dynamics Simulations

A fully threaded tree (FTT) for adaptive refinement of regular meshes is described. By using a tree threaded at all levels, tree traversals for finding nearest neighbors are avoided. All operations on a tree including tree modifications are O(N), where N is a number of cells, and are performed in parallel. An efficient implementation of the tree is described that requires 2N words of memory. A filtering algorithm for removing high frequency noise during mesh refinement is described. A FTT can be used in various numerical applications. In this paper, it is applied to the integration of the Euler equations of fluid dynamics. An adaptive mesh time stepping algorithm is described in which different time steps are used at different l evels of the tree. Time stepping and mesh refinement are interleaved to avoid extensive buffer layers of fine mesh which were otherwise required ahead of moving shocks. Test examples are presented, and the FTT performance is evaluated. The three dimensional simulation of the interaction of a shock wave and a spherical bubble is carried out that shows the development of azimuthal perturbations on the bubble surface.

FINITE ELEMENT ANALYSIS↗

Photographic Evidence of Streamwise Arrays of Vortices in Boundary-Layer Flow

Photographs are presented of various models coated with fluorescent oil to show evidence of surface vortices at a Mach number of 3.03. Vortex formation was evidently present on models with forward-facing steps, rearward-facing steps, wires, discrete surface particles, or unswept flat surfaces with sharp leading edges. Some photographs are also presented for the models coated with a sublimation material which clearly indicates the location of boundary-layer transition; however, it does not show the vortices as clearly as the fluorescent oil. The study was made on the models at an angle of attack of 0 deg at unit Reynolds numbers of 7.7 and 10.7 million per foot. The spacing of the vortices as indicated by the flow studies on the unswept model was smaller at the higher Reynolds number in accordance with Gortler's theory. The flow studies also indicated that stable surface vortices produced by either steps or surface roughness persisted over model areas known to have turbulent boundary layers.

Hopkins, Edward J.↗

Determining the Molecular Growth Mechanisms of Tetragonal Lysozyme Crystals

Studies of the growth of tetragonal lysozyme crystals employing atomic force microscopy (AFM) have shown the advantages of this technique in investigating the growth mechanisms of protein crystals [1]. The resolution of these studies was in the micron range, which revealed surface features such as the occurrence of dislocations and 2D nucleation islands, similar to those found in inorganic systems. They clearly showed that the crystals grew by these surface growth mechanisms. However, the studies also revealed some surprising features, such as bimolecular growth step heights and pronounced growth anisotropies on the (110) face, which could not be explained. In previous studies we employed Periodic Bond Chain (PBC) theory to tetragonal lysozyme crystal growth and found that the crystals were constructed by strongly bonded molecular chains forming helices about the 43 axes [2,3]. The helices were connected to each other with weaker bonds. The growth process was shown to proceed by the formation of these 43 helices, resulting in bimolecular growth steps on the (110) face. It was also shown to explain many other observations on tetragonal lysozyme crystal growth. Although PBC analysis is not a new technique [4], it has not been widely used as the mechanisms predicted from it could not be experimentally verified. In this study the growth process of these crystals was investigated, particularly for the (110) face, employing some newly developed high resolution AFM techniques. These techniques allowed individual lysozyme molecules on the crystal faces to be resolved and predictions from PBC analyses to be tested. The analyses had shown that of the two possible packing arrangements on (110) faces, only one would actually occur. Employing the first of the newly developed techniques, these faces were scanned by high resolution AFM. The resulting images were then compared with the theoretically constructed images for the two possible packing arrangements on the (110) face. The theoretical images were constructed by convolution of the crystal surface shape obtained from crystallographic data with the AFM tip shape. The comparison confirmed the prediction that the molecular packing arrangement of these faces corresponded to that for complete 43 helices. The second AFM technique that was developed was used to follow the growth process by measuring the dimensions of individual growth units on the (110) face. Linescans across a growth step, performed near the saturation limit of the crystals, allowed the growth unit dimensions to be measured. These measurements showed that growth on the (110) face proceeded by the formation of new 43 helices from the addition of at least tetramer units in the [110] direction. In the [001] direction growth proceeded by the addition of various aggregate units corresponding to the 4(sub 3) helices.

Li, Huayu↗

Parallel 3D Multi-Stage Simulation of a Turbofan Engine

A 3D multistage simulation of each component of a modern GE Turbofan engine has been made. An axisymmetric view of this engine is presented in the document. This includes a fan, booster rig, high pressure compressor rig, high pressure turbine rig and a low pressure turbine rig. In the near future, all components will be run in a single calculation for a solution of 49 blade rows. The simulation exploits the use of parallel computations by using two levels of parallelism. Each blade row is run in parallel and each blade row grid is decomposed into several domains and run in parallel. 20 processors are used for the 4 blade row analysis. The average passage approach developed by John Adamczyk at NASA Lewis Research Center has been further developed and parallelized. This is APNASA Version A. It is a Navier-Stokes solver using a 4-stage explicit Runge-Kutta time marching scheme with variable time steps and residual smoothing for convergence acceleration. It has an implicit K-E turbulence model which uses an ADI solver to factor the matrix. Between 50 and 100 explicit time steps are solved before a blade row body force is calculated and exchanged with the other blade rows. This outer iteration has been coined a "flip." Efforts have been made to make the solver linearly scaleable with the number of blade rows. Enough flips are run (between 50 and 200) so the solution in the entire machine is not changing. The K-E equations are generally solved every other explicit time step. One of the key requirements in the development of the parallel code was to make the parallel solution exactly (bit for bit) match the serial solution. This has helped isolate many small parallel bugs and guarantee the parallelization was done correctly. The domain decomposition is done only in the axial direction since the number of points axially is much larger than the other two directions. This code uses MPI for message passing. The parallel speed up of the solver portion (no 1/0 or body force calculation) for a grid which has 227 points axially.

Turner, Mark G.↗

A Hierarchical Approach to Buckling Load Calculations

The advantages of using a hierarchical analysis approach to calculate the buckling load of an axially compressed composite cylindrical shell is demonstrated using an example taken from a recent experimental program. The Delft Interactive Shell DEsign COde (DISDECO) shell design code is used for this hierarchical analysis approach to provide an accurate prediction of the critical buckling load of the actual shell structure. DISDECO includes the influence of the boundary conditions, initial geometric imperfections, the effects of stiffener and load eccentricities, and the effects of prebuckling deformations caused by edge constraints in the analysis. It is shown that the use of DISDECO makes it relatively simple to proceed step by step from simple to more complex models and solution procedures. As a final step in the hierarchical analysis approach, the critical buckling load and the estimated imperfection sensitivity of the shell are verified by conducting an analysis of a large finite element model with one of the current generation two-dimensional shell analysis codes with advanced capabilities needed to represent both geometric and material nonlinearities.

Arbocz, Johann↗

Precision Lower Stratospheric Temperature Monitoring with the MSU: Technique, Validation, and Results, 1979-1991

Microwave Sounding Unit channel 4 data from the TIROS-N (Television and Infrared Operational Satellite)-N series of NOAA (National Oceanic and Atmospheric Administration) satellites are intercalibrated to provide a continuous global record of deep-layer averaged lower stratospheric temperatures during 1979-1991. A 13-year record of temperature anomalies is time averaged into pentads and months on a 2.5 deg. grid. The monthly gridpoint anomalies are validated with ten years of radiosonde data during 1979-88. The calibration stability of each satellite's measurements is evaluated during satellite overlap periods, the longest of which reveal no measurable instrumental drift at the level of 0.01 C yr (exp -1). Intercomparisons between NOAA-6 and NOAA- 7 anomalies indicate monthly gridpoint precision of 0.05 C in the tropics to around 0.10 C in the extratropics, and signal-to-noise ratios precision and stability statistics are much better than have been previously reported by other investigators for MSU channel 4. Pentad precision is about 0.10 C in the tropics to around 0.25 C at high latitudes and signal-to- noise ratios generally over 250 in the tropics and high latitudes, but 100-200 in the middle latitudes. Radiosonde comparisons to the monthly gridpoint anomalies have correlations ranging from 0.90 in the tropics (where the interannual variability is smallest) to as high as 0.99 at high-latitude stations. The corresponding standard error of estimate is generally around 0.3 C. A significant difference in decadal trends is found between the satellite and radiosonde systems, with a step change of 0.21 C (sondes cooler) compared to the satellite measurements. Investigations of the possible sources of the discrepancy lead us to suspect that the gradual transition from on-site calibration of sondes with thermometers to factory calibration of sondes around 1982 might have caused a change in the calibration, although this conclusion must be viewed as tentative. The largest globally averaged temperature variations during 1979-91 occur after the El Chichon (1982) and Pinatubo (1991) volcanic eruptions. These warm events are superimposed upon a net downward trend in temperatures during the period. This cooling trend has more of a step function than linear character, with the step occurring during the El Chichon warm event. It is strongest in polar regions and the Northern Hemisphere middle latitudes. These characteristics are qualitatively consistent with radiative adjustments expected to occur with observed ozone depletions.

Spencer, Roy W.↗

Impact Crises, Mass Extinctions, and Galactic Dynamics: A Unified Theory

A general hypothesis linking mass extinctions of life with impacts of large asteroids and comets is based on astronomical data, impact dynamics, and geological information. The waiting times of large-body impacts on the Earth, derived from the flux of Earth-crossing asteroids and comets, and the estimated size of impacts capable of causing large-scale environmental disasters predict that impacts of objects (sup 3)5 km in diameter ((sup 3)10(exp 7) Mt TNT equivalent) could be sufficient to explain the record of about 25 extinction pulses in the last 540 m.y., with the five recorded major mass extinctions related to the impacts of the largest objects of (sup 3)10 km in diameter ( (sup 3)10(exp 8) Mt events). Smaller impacts (about 10(exp 6)-10(exp 7) Mt), with significant regional and even global environmental effects, could be responsible for the lesser boundaries in the geologic record. Tests of the "kill curve" relationship for impact-induced extinctions based on new data on extinction intensities and several well-dated large impact craters suggest that major mass extinctions require large impacts, and that a step in the kill curve may exist at impacts that produce craters of -100 km diameter, with smaller impacts capable of only relatively weak extinction pulses. Single impact craters < about 60 km in diameter should not be associated with global extinction pulses detectable in the Sepkoski database (although they may explain stage and zone boundaries marked by lesser faunal turnover), but multiple impacts in that size range may produce significant stepped extinction pulses. Statistical tests of the last occurrences of species at mass-extinction boundaries are generally consistent with predictions for abrupt or stepped extinctions, and several boundaries are known to show "catastrophic" signatures of environmental disasters and biomass crash, impoverished postextinction fauna and flora dominated by stress-tolerant and opportunistic species, and gradual ecological recovery and radiation of new taxa. Isotopic and other geochemical signatures are also generally consistent with the expected after-effects of catastrophic impacts. Seven of the recognized extinction pulses are associated with concurrent (in some cases multiple) stratigraphic impact markers (e.g., layers with high Ir, shocked minerals, microtektites), and/or large, dated impact craters. Other less-well-studied crisis intervals show elevated Ir, still well below that of the K/T spike, which might be explained by low-Ir impactors, ejecta blowoff, or the sedimentary reworking and dilution of impact signatures. The best explanation for a possible periodic component of about 30 m.y. in mass extinctions and clusters of impacts is the modulation of the comet flux associated with the solar system's periodic passage through the plane of the Milky Way Galaxy. The quantitative agreement among paleontological, geological, and astronomical data suggests an important underlying unification of the processes involved.

Rampino, M.R.↗

Maximum Normalized Rate as a Flying Qualities Parameter

Discrete attitude commands have become a standard task for flying qualities evaluation and control system testing. Much pilot opinion data is now available for ground-based and in-flight simulations, but adequate performance measures and prediction methods have not been established. The Step Target Tracking Prediction method, introduced in 1978, correlated time-on-target and rms tracking data with NT-33 in-flight longitudinal simulations, but did not employ parameters easily measured in manned flight and simulation. Recent application of the Step Target Tracking Prediction method to lateral flying qualities analysis has led to a new measure of performance. This quantity, called Maximum Normalized Rate (MNR), reflects the greatest attitude rate a pilot can employ during a discrete maneuver without excessive overshoot and oscillation. MNR correlates NT-33 lateral pilot opinion ratings well, and is easily measured during night test or simulation. Furthermore, the Step Target MNR method can be used to analyze large amplitude problems concerning rate limiting and nonlinear aerodynamics.

Onstott, E. D.↗