Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulator study”

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 955 records · Page 53

Ground-to-Flight Handling Qualities Comparisons for a High Performance Airplane

A flight test program was conducted in conjunction with a ground-based piloted simulation study to enable a comparison of handling qualities ratings for a variety of maneuvers between flight and simulation of a modern high performance airplane. Specific objectives included an evaluation of pilot-induced oscillation (PIO) tendencies and a determination of maneuver types which result in either good or poor ground-to-flight pilot handling qualities ratings. A General Dynamics F-16XL aircraft was used for the flight evaluations, and the NASA Langley Differential Maneuvering Simulator was employed for the ground based evaluations. Two NASA research pilots evaluated both the airplane and simulator characteristics using tasks developed in the simulator. Simulator and flight tests were all conducted within approximately a one month time frame. Maneuvers included numerous fine tracking evaluations at various angles of attack, load factors and speed ranges, gross acquisitions involving longitudinal and lateral maneuvering, roll angle captures, and an ILS task with a sidestep to landing. Overall results showed generally good correlation between ground and flight for PIO tendencies and general handling qualities comments. Differences in pilot technique used in simulator evaluations and effects of airplane accelerations and motions are illustrated.

Brandon, Jay M.↗

Experiment M-1: Cardiovascular Conditioning

Intermittent venous occlusion of the extremities of man, during weightlessness simulation studies, has been demonstrated to be effective in preventing or mitigating the orthostatic hypotension observed following such simulations. A similar preventive measure was employed on the pilot of Gemini V with a view to determining the efficacy of pulsatile leg cuffs in preventing or lessening the orthostatic hypotension observed following previous space flights. Unfortunately, the cuff device was operative continuously during only the first 4 days of the 8-day mission. Post flight tilt-table responses of the command pilot and the pilot were considerably different, but the data cannot be construed as a conclusive demonstration that the observed differences were the result of the action of the pulsatile cuffs. The differences in the tilt responses of the Gemini V flight crew may be only a reflection of individual variability so commonly observed in biological experimentation. More data shall be required before a judgement can be rendered as to the efficacy of the pulsatile-leg-cuff technique in lessening postflight postural or orthostatic hypotension.

Dietlein, Lawrence F.↗

Intelligent Systems Approach for Automated Identification of Individual Control Behavior of a Human Operator

Results have been obtained using conventional techniques to model the generic human operator?s control behavior, however little research has been done to identify an individual based on control behavior. The hypothesis investigated is that different operators exhibit different control behavior when performing a given control task. Two enhancements to existing human operator models, which allow personalization of the modeled control behavior, are presented. One enhancement accounts for the testing control signals, which are introduced by an operator for more accurate control of the system and/or to adjust the control strategy. This uses the Artificial Neural Network which can be fine-tuned to model the testing control. Another enhancement takes the form of an equiripple filter which conditions the control system power spectrum. A novel automated parameter identification technique was developed to facilitate the identification process of the parameters of the selected models. This utilizes a Genetic Algorithm based optimization engine called the Bit-Climbing Algorithm. Enhancements were validated using experimental data obtained from three different sources: the Manual Control Laboratory software experiments, Unmanned Aerial Vehicle simulation, and NASA Langley Research Center Visual Motion Simulator studies. This manuscript also addresses applying human operator models to evaluate the effectiveness of motion feedback when simulating actual pilot control behavior in a flight simulator.

Zaychik, Kirill B.↗

Analysis of Handling Qualities Design Criteria for Active Inceptor Force-Feel Characteristics

The force-feel system characteristics of the cyclic inceptors of most helicopters are set based on the characteristics of the mechanical components in the control system (mass, springs, friction dampers, etc.). For these helicopters, the force-feel characteristics typically remain constant over the entire flight envelope, with perhaps a trim release to minimize control forces while maneuvering. With the advent of fly-by-wire control systems and active inceptors in helicopters, the force-feel characteristics are now determined by the closed-loop response of the active inceptor itself as defined by the inertia, force/displacement gradient, damping, breakout force and detent shape configuration parameters in the inceptor control laws. These systems give the flexibility to dynamically prescribe different feel characteristics for different control modes or flight conditions, and the ability to provide tactile cueing to the pilot through the actively controlled side-stick or center-stick cyclic inceptor. For rotorcraft, a few studies have been conducted to assess the effects of cyclic force-feel characteristics on handling qualities in flight. An early study provided valuable insight into the static force-deflection characteristics (force gradient) and the number of axes controlled by the side-stick controller for the U.S. Army's Advanced Digital/Optical Control System (ADOCS) demonstrator aircraft [1]. The first of a series of studies providing insight on the inceptor dynamic force-feel characteristics was conducted on the NASA/Army CH-47B variable-stability helicopter [2]. This work led to a proposed requirement that set boundaries based on the cyclic natural frequency and inertia, with the stipulation of a lower damping ratio limit of 0.3 [3]. A second study was conducted by the Canadian Institute for Aerospace Research using their variable-stability Bell 205A helicopter [4]. This research suggested boundaries for stick dynamics based on natural frequency and damping ratio. While these two studies produced boundaries for acceptable/unacceptable stick dynamics for rotorcraft, they were not able to provide guidance on how variations of the stick dynamics in the acceptable region impact handling qualities. More recently, a ground based simulation study [5] suggested little benefit was to be obtained from variations of the damping ratio for a side-stick controller exhibiting high natural frequencies (greater than 17 rad/s) and damping ratios (greater than 2.0). A flight test campaign was conducted concurrently on the RASCAL JUH-60A in-flight simulator and the ACT/FHS EC-135 in flight simulator [6]. Upon detailed analysis of the pilot evaluations the study identified a clear preference for a high damping ratio and natural frequency of the center stick inceptors. Side stick controllers were found to be less sensitive to the damping. While these studies have compiled a substantial amount of data, in the form of qualitative and quantitative pilot opinion, a fundamental analysis of the effect of the inceptor force-feel system on flight control is found to be lacking. The study of Ref. [6] specifically concluded that a systematic analysis was necessary, since discrepancies with the assigned handling qualities showed that proposed analytical design metrics, or criteria, were not suitable. The overall goal of the present study is to develop a clearer fundamental understanding of the underlying mechanisms associated with the inceptor dynamics that govern the handling qualities using a manageable analytical methodology.

Malpica, Carlos A.↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

Analysis of simulated measurements of earth emitted radiation using geometric shape factors, and some results from analysis of ESSA 7 and Nimbus 6 data

In the reported study, simulated measurements of earth radiation with wide field-of-view (WFOV) radiometers are reduced to earth fluxes through geometric shape factors. Estimates of emitted flux for global and zonal scales and for regions smaller than the detector field-of-view are compared to fluxes derived directly from a model. Geometric shape factors are used to estimate global, zonal, and regional averages of earth-emitted flux from measurements with WFOV detectors on the ESSA 7 satellite and the Nimbus 6 satellite. Zonal and global averages of emitted flux for each of the 3-month periods of ESSA 7 data are shown in a graph. Another graph shows zonal and global averages of emitted flux for August 1975 obtained with the aid of Nimbus 6. The data are compared to values reported by Green and Smith (1977) who used a spherical harmonic solution.

Weaver, W. L.↗

Autonomous Aerobraking: Thermal Analysis and Response Surface Development

A high-fidelity thermal model of the Mars Reconnaissance Orbiter was developed for use in an autonomous aerobraking simulation study. Response surface equations were derived from the high-fidelity thermal model and integrated into the autonomous aerobraking simulation software. The high-fidelity thermal model was developed using the Thermal Desktop software and used in all phases of the analysis. The use of Thermal Desktop exclusively, represented a change from previously developed aerobraking thermal analysis methodologies. Comparisons were made between the Thermal Desktop solutions and those developed for the previous aerobraking thermal analyses performed on the Mars Reconnaissance Orbiter during aerobraking operations. A variable sensitivity screening study was performed to reduce the number of variables carried in the response surface equations. Thermal analysis and response surface equation development were performed for autonomous aerobraking missions at Mars and Venus.

Dec, John A.↗

An Evaluation of Automatic Control System Concepts for General Aviation Airplanes

A piloted simulation study of automatic longitudinal control systems for general aviation airplanes has been conducted. These automatic control systems were designed to make the simulated airplane easy to fly for a beginning or infrequent pilot. Different control systems are presented and their characteristics are documented. In a conventional airplane control system each cockpit controller commands combinations of both the airspeed and the vertical speed. The best system in the present study decoupled the airspeed and vertical speed responses to cockpit controller inputs. An important feature of the automatic system was that neither changing flap position nor maneuvering in steeply banked turns affected either the airspeed or the vertical speed. All the pilots who flew the control system simulation were favorably impressed with the very low workload and the excellent handling qualities of the simulated airplane.

Stewart, E. C.↗

Initial Validation of a Simulation System for Studying Interoperability in Future Air Traffic Management Systems

Future air traffic management systems will need to accommodate large numbers of increasingly diverse air vehicles with different operating paradigms. To support this trend, they will digitally share copious amounts of information via a common communication architecture. Operators will deploy programs that create and negotiate flight plans via the architecture’s communication protocols. These programs will autonomously make decisions that must be arbitrated by the architecture and robust to uncertainty. To study interoperability in air traffic management, a new airspace simulation system was composed by integrating a legacy airspace simulation, an air traffic control model, and a new research communication architecture. It was used to evaluate air traffic management concepts by adapting it to handle congested arrival traffic at Newark Liberty International Airport and executing simulations. Results demonstrated the ability of the simulation system to model in detail strategic traffic flow management, predeparture flight planning, and air traffic control working in concert. Subsequent studies can use the simulation system to study interoperability, autonomy, digital communication, and uncertainty in future air traffic management systems.

aircraft scheduling,traffic flow management,autono↗

Initial Validation of a Simulation System for Studying Interoperability in Future Air Traffic Management Systems

Future air traffic management systems will need to accommodate large numbers of increasingly diverse air vehicles with different operating paradigms. To support this trend, they will digitally share copious amounts of information via a common communication architecture. Operators will deploy programs that create and negotiate flight plans via the architecture’s communication protocols. These programs will autonomously make decisions that must be arbitrated by the architecture and robust to uncertainty. To study interoperability in air traffic management, a new airspace simulation system was composed by integrating a legacy airspace simulation, an air traffic control model, and a new research communication architecture. It was used to evaluate air traffic management concepts by adapting it to handle congested arrival traffic at Newark Liberty International Airport and executing simulations. Results demonstrated the ability of the simulation system to model in detail strategic traffic flow management, predeparture flight planning, and air traffic control working in concert. Subsequent studies can use the simulation system to study interoperability, autonomy, digital communication, and uncertainty in future air traffic management systems.

air traffic control↗

Initial Validation of a Simulation System for Studying Interoperability in Future Air Traffic Management Systems

Future air traffic management systems will need to accommodate large numbers of increasingly diverse air vehicles with different operating paradigms. To support this trend, they will digitally share copious amounts of information via a common communication architecture. Operators will deploy programs that create and negotiate flight plans via the architecture’s communication protocols. These programs will autonomously make decisions that must be arbitrated by the architecture and robust to uncertainty. To study interoperability in air traffic management, a new airspace simulation system was composed by integrating a legacy airspace simulation, an air traffic control model, and a new research communication architecture. It was used to evaluate air traffic management concepts by adapting it to handle congested arrival traffic at Newark Liberty International Airport and executing simulations. Results demonstrated the ability of the simulation system to model in detail strategic traffic flow management, predeparture flight planning, and air traffic control working in concert. Subsequent studies can use the simulation system to study interoperability, autonomy, digital communication, and uncertainty in future air traffic management systems.

air traffic control↗

Evaluating Alerting and Guidance Performance of a UAS Detect-And-Avoid System

A key challenge to the routine, safe operation of unmanned aircraft systems (UAS) is the development of detect-and-avoid (DAA) systems to aid the UAS pilot in remaining "well clear" of nearby aircraft. The goal of this study is to investigate the effect of alerting criteria and pilot response delay on the safety and performance of UAS DAA systems in the context of routine civil UAS operations in the National Airspace System (NAS). A NAS-wide fast-time simulation study was conducted to assess UAS DAA system performance with a large number of encounters and a broad set of DAA alerting and guidance system parameters. Three attributes of the DAA system were controlled as independent variables in the study to conduct trade-off analyses: UAS trajectory prediction method (dead-reckoning vs. intent-based), alerting time threshold (related to predicted time to LoWC), and alerting distance threshold (related to predicted Horizontal Miss Distance, or HMD). A set of metrics, such as the percentage of true positive, false positive, and missed alerts, based on signal detection theory and analysis methods utilizing the Receiver Operating Characteristic (ROC) curves were proposed to evaluate the safety and performance of DAA alerting and guidance systems and aid development of DAA system performance standards. The effect of pilot response delay on the performance of DAA systems was evaluated using a DAA alerting and guidance model and a pilot model developed to support this study. A total of 18 fast-time simulations were conducted with nine different DAA alerting threshold settings and two different trajectory prediction methods, using recorded radar traffic from current Visual Flight Rules (VFR) operations, and supplemented with DAA-equipped UAS traffic based on mission profiles modeling future UAS operations. Results indicate DAA alerting distance threshold has a greater effect on DAA system performance than DAA alerting time threshold or ownship trajectory prediction method. Further analysis on the alert lead time (time in advance of predicted loss of well clear at which a DAA alert is first issued) indicated a strong positive correlation between alert lead time and DAA system performance (i.e. the ability of the UAS pilot to maneuver the unmanned aircraft to remain well clear). While bigger distance thresholds had beneficial effects on alert lead time and missed alert rate, it also generated a higher rate of false alerts. In the design and development of DAA alerting and guidance systems, therefore, the positive and negative effects of false alerts and missed alerts should be carefully considered to achieve acceptable alerting system performance by balancing false and missed alerts. The results and methodology presented in this study are expected to help stakeholders, policymakers and standards committees define the appropriate setting of DAA system parameter thresholds for UAS that ensure safety while minimizing operational impacts to the NAS and equipage requirements for its users before DAA operational performance standards can be finalized.

UAS Detect-and-Avoid (DAA) System↗

Solar Dynamic Power System Fault Diagnosis

The objective of this research is to conduct various fault simulation studies for diagnosing the type and location of faults in the power distribution system. Different types of faults are simulated at different locations within the distribution system and the faulted waveforms are monitored at measurable nodes such as at the output of the DDCU's. These fault signatures are processed using feature extractors such as FFT and wavelet transforms. The extracted features are fed to a clustering based neural network for training and subsequent testing using previously unseen data. Different load models consisting of constant impedance and constant power are used for the loads. Open circuit faults and short circuit faults are studied. It is concluded from present studies that using features extracted from wavelet transforms give better success rates during ANN testing. The trained ANN's are capable of diagnosing fault types and approximate locations in the solar dynamic power distribution system.

Momoh, James A.↗

Improving precision and accuracy of genetic mapping with genotyping‐by‐sequencing data in outcrossing species

Abstract Genotyping‐by‐sequencing (GBS) is a widely used strategy for obtaining large numbers of genetic markers in model and non‐model organisms. In crop plants, GBS‐derived marker datasets are frequently used to perform quantitative trait locus (QTL) mapping. In some plant species, however, high heterozygosity and complex genome structure mean that researchers must use care in handling GBS data to conduct QTL mapping most effectively. Such outbred crops include most of the perennial grass and tree species used for bioenergy. To identify strategies for increasing accuracy and precision of QTL mapping using GBS data in outbred crops, we conducted an empirical study of SNP‐calling and genetic map‐building pipeline parameters in a Miscanthus sinensis population, and a complementary simulation study to estimate the relationship between genome‐wide error rate, read depth, and marker number. The bioenergy grass Miscanthus is an obligate outcrossing species with a recent (diploidized) whole‐genome duplication. For the study of empirical M. sinensis data, we compared two SNP‐calling methods (one non‐reference‐based and one reference‐based), a series of depth filters (12×, 20×, 30×, and 40×) and two map‐construction methods (i.e., marker ordering: linkage‐only and order‐corrected based on a reference genome). We found that correcting the order of markers on a linkage map by using a high‐quality reference genome improved QTL precision (shorter confidence intervals). For typical GBS datasets of between 1000 and 5000 markers to build a genetic map for biparental populations, a depth filter set at 30× to 40× applied to outbred populations provided a genome‐wide genotype‐calling error rate of less than 1%, improved accuracy of QTL point estimates and minimized type I errors for identifying QTL. Based on these results, we recommend using a reference genome to correct the marker order of genetic maps and a robust genotype depth filter to improve QTL mapping for outbred crops.

59 BASIC BIOLOGICAL SCIENCES↗

Real-World Driving Features for Identifying Intelligent Driver Model Parameters

Driver behavior models play a significant role in representing different driving styles and the associated relationships with traffic patterns and vehicle energy consumption in simulation studies. The models often serve as a proxy for baseline human driving when assessing energy-saving strategies that alter vehicle velocity. Such models are especially important in connectivity-enabled energy-saving strategy research because they can easily adapt to changing driving conditions like posted speed limits or change in traffic light state. While numerous driver models exist, parametric driver models provide the flexibility required to represent variability in real-world driving through different combinations of model parameters. These model parameters must be informed by a representative set of parameter values for the driver model to adequately represent a real-world driver. It stands to reason that determining the parameter values from real-world driving data would serve the purpose of representing a real-world driver. Although the published literature is replete with techniques and consequences of using real-world driving data to determine parameter values for parametric driver models, none have explored them in the context of using shorter driving features where the parameter values may change over the course of a single trip for the same driver. In this study we consider the “intelligent driver model” (IDM) as our driver behavior model and use real-world driving data from the Transportation Secure Data Center (TSDC) maintained at the National Renewable Energy Laboratory (NREL). The TSDC includes real-world travel data from across the United States, from which NREL has created a wide range of driving routes consisting of road features such as speed limits, stop locations, and turn locations. The real-world driving data are categorized into different driving regimes and extracted into driving features. The driving features are then used to calibrate the parameter values for the IDM. The distribution of the parameters and the relationships among them are reported. The insights obtained from this study enable judicious usage of IDM or similar parametric driver models to represent baseline human driver behavior in simulations.

27 ARPA - Advanced Research Projects Agency-Energy↗

Pilot performance during a simulated standard instrument procedure turn with and without a predictor display

A simulator study was conducted to measure the effectiveness of predictor information incorporated into a CRT display of a computer simulated aircraft's horizontal and vertical situation. Professional pilots served as subjects for the task of executing a standard instrument procedure turn at constant altitude in constant crosswinds with and without their predicted ground track displayed. The results showed that the display with the predicted ground track was markedly and significantly superior to the display without this information and that the subjects were generally satisfied with this type of information. Mean rms lateral path error was independent of the crosswind velocity with the predictor information, and increased without it and with increasing wind velocity. Rms stick activity decreased with the predictor display which also uncoupled aileron and elevator activity.

Kreifeldt, J. G.↗

Pilot Performance During A Simulated Standard Instrument Procedure Turn With and Without A Predictor Display

A simulator study was conducted to measure the effectiveness of predictor information incorporated into a CRT display of a computer simulated aircraft's horizontal and vertical situation. Professional pilots served as subjects for the task of executing a standard instrument procedure turn at constant altitude in constant crosswinds with and without their predicted ground track displayed. The results show that the display with the predicted ground track was markedly and significantly superior to the display without this information and that the subjects were generally satisfied with this type of information. Mean rms lateral path error was independent of the crosswind velocity with the predictor information, and increased without it with increasing wind velocity. Rms stick activity decreased with the predictor display which also uncoupled aileron and elevator activity.

John G Kreifeldt↗

A Mathematical Model for Vertical Attitude Takeoff and Landing (VATOL) Aircraft Simulation: Model Description Application - Volume 1

A mathematical model of a high performance airplane capable of vertical attitude takeoff and landing (VATOL) was developed. An off line digital simulation program incorporating this model was developed to provide trim conditions and dynamic check runs for the piloted simulation studies and support dynamic analyses of proposed VATOL configuration and flight control concepts. Development details for the various simulation component models and the application of the off line simulation program, Vertical Attitude Take-Off and Landing Simulation (VATLAS), to develop a baseline control system for the Vought SF-121 VATOL airplane concept are described.

Fortenbaugh, R. L.↗