Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear programming”

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,495 records · Page 83

Numerical solution of differential equations by artificial neural networks

Conventionally programmed digital computers can process numbers with great speed and precision, but do not easily recognize patterns or imprecise or contradictory data. Instead of being programmed in the conventional sense, artificial neural networks (ANN's) are capable of self-learning through exposure to repeated examples. However, the training of an ANN can be a time consuming and unpredictable process. A general method is being developed by the author to mate the adaptability of the ANN with the speed and precision of the digital computer. This method has been successful in building feedforward networks that can approximate functions and their partial derivatives from examples in a single iteration. The general method also allows the formation of feedforward networks that can approximate the solution to nonlinear ordinary and partial differential equations to desired accuracy without the need of examples. It is believed that continued research will produce artificial neural networks that can be used with confidence in practical scientific computing and engineering applications.

Meade, Andrew J., Jr.↗

SONEX Airborne Mission and Coordinated POLINAT-2 Activity: Overview and Accomplishments

The SASS (Subsonic Assessment) Ozone and NO(x) Experiment (SONEX) was an airborne field campaign conducted in October-November 1997 in the vicinity of the North Atlantic Flight Corridor Lo study the impact of aircraft emissions on NOx and ozone (03). A fully instrumented NASA DC-8 aircraft was used as the primary SONEX platform. SONEX activities were closely coordinated with the European POLINAT-2 (Pollution from Aircraft Emissions in the North Atlantic Flight Corridor) program, which used a Falcon-20 aircraft and an instrumented in-service Swissair B-747. Both campaigns focused on the upper troposphere/"lowermost" stratosphere (UT/LS) as the region of greatest interest. Specific sampling goals were achieved with the aid of a state-of-the art modeling and meteorological support system, which allowed targeted sampling of air parcels with desired characteristics. A substantial impact of aircraft emissions on NO(x) and O3 in the UT/LS of the study region is shown to be present. It is further shown that the NO(x)- HO(x)-O3 relationships are highly nonlinear and must be accurately simulated to make meaningful future predictions with global models. SONEXIPOLINAT-2 results are being published in Special Sections of GRL and JGR. Here we provide a brief overview of SONEX design, implementation, and expected results to provide a context within which these publications can be understood.

Singh, Hanwant B.↗

Dynamic analysis of flexible multibody structures

A system composed of several interconnected elastic components that may experience large angular motion relative to each other during operation is referred to as a flexible multibody structure. Several formulations were proposed for the determination of the dynamic response of controlled flexible multibody structures. In general, these formulations consist of superposing elastic deformations of the component body onto the large rigid body motion of the component. It was shown that this particular methodology for combining linear structural deformations with nonlinear kinematics can lead to erroneous response predictions when either the beam member is very flexible or the rotational speed is high. In addition, previous formulations introduce constraint equations to define the interrelations among system components. This approach increases the number of equations that must be solved, and may result in contraint violation when numerical error accumulates during the integration process. In order to overcome the difficulties, a new approach was suggested. The approach is essentially a finite element formulation which takes advantage of the fact that many multibody structures are joint dominated. The Large Angle Transient Dynamic Analysis (LATDYN) program for clarity of documentation, ease of use, user friendliness, modeling generality, and accuracy of results was evaluated. This required gaining a working familiarity with the code and performing several case studies.

Hernried, Alan G.↗

Correlation of Structural Analysis and Test Results for the McDonnell Douglas Stitched/RFI All-Composite Wing Stub Box

The analytical and experimental results of an all-composite wing stub box are presented in this report. The wing stub box, which is representative of an inboard portion of a commercial transport high-aspect-ratio wing, was fabricated from stitched graphite-epoxy material with a Resin Film Infusion manufacturing process. The wing stub box was designed and constructed by the McDonnell Douglas Aerospace Company as part of the NASA Advanced Composites Technology program. The test article contained metallic load-introduction structures on the inboard and outboard ends of the graphite-epoxy wing stub box. The root end of the inboard load introduction structure was attached to a vertical reaction structure, and an upward load was applied to the outermost tip of the outboard load introduction structure to induce bending of the wing stub box. A finite element model was created in which the center portion of the wing-stub-box upper cover panel was modeled with a refined mesh. The refined mesh was required to represent properly the geometrically nonlinear structural behavior of the upper cover panel and to predict accurately the strains in the stringer webs of the stiffened upper cover panel. The analytical and experimental results for deflections and strains are in good agreement.

Wang, John T.↗

Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments

The new Bayesian Monte Carlo (MC) evaluation framework described in this abstract has been conceived as an attempt to improve nuclear data evaluations of differential crosssection data by removing the following two approximations conventionally employed for nuclear data evaluations: all probability density functions (PDFs) of all data and model parameters, both prior and posterior, are assumed to be normal (i.e., Gaussian) PDFs, and all uncertainties and covariances are propagated using a linear approximation. With these approximations removed, the Bayesian MC (BMC) framework could be used to account for nonlinear effects and would enable improved evaluations of differential cross sections and IBE data that are presently performed based on the assumptions itemized above. The BMC would also improve upon the uniform sampling of IBE parameters from within ranges defined by their evaluated uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Metal colloids and semiconductor quantum dots: Linear and nonlinear optical properties

One aspect of this project involves a collaborative effort with the Solid State Division of ORNL. The thrust behind this research is to develop ion implantion for synthesizing novel materials (quantum dots wires and wells, and metal colloids) for applications in all optical switching devices, up conversion, and the synthesis of novel refractory materials. In general the host material is typically a glass such as optical grade silica. The ions of interest are Au, Ag, Cd, Se, In, P, Sb, Ga and As. An emphasis is placed on host guest interactions between the matrix and the implanted ion and how the matrix effects and implantation parameters can be used to obtain designer level optical devices tailored for specific applications. The specific materials of interest are: CdSe, CdTe, InAs, GaAs, InP, GaP, InSb, GaSb and InGaAs. A second aspect of this research program involves using porous glass (25-200 A) for fabricating materials of finite size. In this part of the program, we are particularly interested in characterizing the thermodynamic and optical properties of these non-composite materials. We also address how phase diagram of the confined material is altered by the interfacial properties between the confined material and the pore wall.

Henderson, D. O.↗

Stacking-sequence optimization for buckling of laminated plates by integer programming

Integer-programming formulations for the design of symmetric and balanced laminated plates under biaxial compression are presented. Both maximization of buckling load for a given total thickness and the minimization of total thickness subject to a buckling constraint are formulated. The design variables that define the stacking sequence of the laminate are zero-one integers. It is shown that the formulation results in a linear optimization problem that can be solved on readily available software. This is in contrast to the continuous case, where the design variables are the thicknesses of layers with specified ply orientations, and the optimization problem is nonlinear. Constraints on the stacking sequence such as a limit on the number of contiguous plies of the same orientation and limits on in-plane stiffnesses are easily accommodated. Examples are presented for graphite-epoxy plates under uniaxial and biaxial compression using a commercial software package based on the branch-and-bound algorithm.

Haftka, Raphael T.↗

Full-Scale Crash Test and Finite Element Simulation of a Composite Prototype Helicopter

A full-scale crash test of a prototype composite helicopter was performed at the Impact Dynamics Research Facility at NASA Langley Research Center in 1999 to obtain data for validation of a finite element crash simulation. The helicopter was the flight test article built by Sikorsky Aircraft during the Advanced Composite Airframe Program (ACAP). The composite helicopter was designed to meet the stringent Military Standard (MIL-STD-1290A) crashworthiness criteria and was outfitted with two crew and two troop seats and four anthropomorphic dummies. The test was performed at 38-ft/s vertical and 32.5-ft/s horizontal velocity onto a rigid surface. An existing modal-vibration model of the Sikorsky ACAP helicopter was converted into a model suitable for crash simulation. A two-stage modeling approach was implemented and an external user-defined subroutine was developed to represent the complex landing gear response. The crash simulation was executed with a nonlinear, explicit transient dynamic finite element code. Predictions of structural deformation and failure, the sequence of events, and the dynamic response of the airframe structure were generated and the numerical results were correlated with the experimental data to validate the simulation. The test results, the model development, and the test-analysis correlation are described.

Jackson, Karen E.↗

Preliminary Sensorimotor and Cardiovascular Results from the Joint Russian and U.S. Pilot Field Test with Planning for the Full Field Test Beginning with the Year Long Intenational Space Station

Ongoing collaborative research efforts between NASA's Neuroscience and Cardiovascular Laboratories, and the Institute of Biomedical Problems' (IBMP) Sensory-Motor and Countermeasures Laboratories have been measuring functional sensorimotor, cardiovascular and strength responses following bed rest, dry immersion, short duration (Space Shuttle) and long duration (Mir and International Space Station) space flights. While the unloading paradigms associated with dry immersion and bed rest does serve as acceptable flight analogs, testing of crew responses following the long duration flights previously has not been possible until a minimum of 24 hours after landing. As a result, it is not possible to estimate the nonlinear trend of the early (<24 hr) recovery process, nor is it possible to accurately assess the full impact of the decrements associated with long duration flight. To overcome these limitations, both the Russian and U.S. sides have implemented testing at landing site. By joint agreement, this research effort has been identified as the functional Field Test (FT). For practical reasons the FT has been divided into two phases: the full FT and a preliminary pilot version (PFT) of the FT that is reduced in both length and scope. The primary goal of this research is to determine functional abilities in long duration space flight crews beginning as soon after landing as possible (< 2 hr) with one to three immediate follow-up measurements on the day of landing. This goal has both sensorimotor and cardiovascular elements, including evaluations of NASA's new anti-orthostatic compression garment and the Russian Kentavr garment. Functional sensorimotor measurements will include, but are not limited to, assessment of hand/eye coordination, ability to egress from a seated position, walk normally without falling, measurement of dynamic visual acuity, ability to discriminate different forces generated with both the hands and legs, recovery from a fall, a coordinated walk involving tandem heel-to-toe placement, and determination of postural ataxia while standing. The cardiovascular portion of the investigation includes blood pressure and heart rate measurements during a timed stand test in conjunction with postural ataxia testing (quiet stance sway) as well as cardiovascular responses during other functional tasks. In addition to the immediate post-landing collection of data for the full FT, postflight data will be acquired at a minimum of one to three more other times within the 24 hr following landing and continue until functional sensorimotor and cardiovascular responses have returned to preflight normative values. The PFT represents a single trial run comprised of jointly agreed upon subset of tests from the full FT and relies heavily on IBMP's Sensory-Motor and Countermeasures Laboratories for content and implementation. The PFT was first conducted following the September 2013 landing of the Soyuz spacecraft (34S) and again following the landing of Soyuz 35S in November. Testing included: (1) a sit-tostand test, (2) recovery from a fall where the crewmember began in the prone position on the ground and then stood for 3 min while cardiovascular stability was determined and postural ataxia data were acquired, and (3) a tandem heel-to-toe walk test to determine changes in the central locomotor program. Video, cardiovascular parameters (heart rate and blood pressure), data from body-worn inertial sensors and severity of postflight motion sickness were collected for analysis. In summary, the level of functional deficit is expected to be most profound during the acquisition of gravity loads immediately after landing when the demands for crew intervention in response to emergency operations will be greatest. Clearly measureable performance parameters such as ability to perform a seat egress, recover from a fall or the ability to see clearly when walking, and related physiological data (orthostatic responses) are required to provide an evidence base for characterizing programmatic risks and the degree of variability among crewmembers for exploration missions where the crew will be unassisted after landing. Overall, these early functional and related physiological measurements will allow estimation of nonlinear sensorimotor and cardiovascular recovery trends to an accuracy that has not been previously captured in over 50 years of space flight.

Reschke, M. F.↗

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗

Investigation of Nonlinear Difference-Frequency Wave Excitation on a Semisubmersible Offshore-Wind Platform with Bichromatic-Wave CFD Simulations: Preprint

The natural surge and pitch frequencies of semisubmersible offshore wind platforms are typically designed to be below the wave frequencies to avoid direct excitation. However, surge or pitch resonance can be excited by the nonlinear low-frequency loads generated by irregular incident waves. Second-order potential-flow models with added Morison drag have been found to underpredict this low-frequency excitation and response. As part of the OC6 project, Computational Fluid Dynamics (CFD) simulations were performed to enable a better understanding of the low-frequency loads and the limitations of lower-fidelity models. The focus of the current paper is to set up a computationally cost-effective CFD simulation of a fixed semisubmersible platform to investigate nonlinear difference-frequency loads and establish the corresponding uncertainty in the results. Due to the high computing cost, CFD simulations of irregular waves can be challenging. Instead, simulations were performed with bichromatic waves having a shorter repeat period. A preliminary comparison with Quadratic Transfer Functions from second-order potential-flow theory shows that CFD models consistently predict higher nonlinear wave loads at the difference frequency, likely due to flow separation and viscous drag not accounted for in potential-flow theory.

2nd order↗

Single-Atom Doping at the Molecule–Metal Interface: How Rh Affects Surface Explosion Kinetics

Controlling chemical reactions at interfaces is central to many fields, including atmospheric, environmental, and biological chemistry as well as catalysis and corrosion. We performed a fundamental study of how catalytically active Rh atoms placed at and above a molecule-metal interface influence the surface reaction kinetics of tartaric acid (TA) decomposition. Specifically, TA decomposition on Cu(100) exhibits autocatalytic decomposition kinetics involving a slow initiation step followed by fast decomposition of the molecular layer. These so-called “surface explosions” are extremely sensitive to initial conditions, making it fundamentally interesting to study how the placement of reactive atoms affects the reaction rate. Temperature-programmed reaction (TPR) experiments reveal that Rh atoms embedded in the Cu(100) surface beneath TA (TA/Rh/Cu) or atop the TA layer (Rh/TA/Cu) reduce the decomposition temperature and modify the reaction rate constants. Isothermal TPR analysis reveals that Rh enhances both the initiation rate constant, k i , and the explosion rate constant, k e , facilitating explosive decomposition at lower temperatures. This result provides evidence that both the initiation and explosion steps occur at the metal-molecule interface and are accelerated by the presence of Rh at this interface. This study illustrates how small amounts of reactive elements may be used to control nonlinear kinetic processes at interfaces.

desorption↗

A New Hybrid Quantum-Classical Algorithm for Solving the Unit Commitment Problem

Solving problems related to planning and operations of large-scale power systems is challenging on classical computers due to their inherent nature as mixed-integer and nonlinear problems. Quantum computing provides new avenues to approach these problems. We develop a hybrid quantum-classical algorithm for the Unit Commitment (UC) problem in power systems which aims at minimizing the total cost while optimally allocating generating units to meet the hourly demand of the power loads. The hybrid algorithm combines a variational quantum algorithm (VQA) with a classical Benders-type heuristic. The resulting algorithm computes approximate solutions to UC in three stages: i) a collection of UC vectors capable meeting the power demand with lowest possible operating costs is generated based on VQA; ii) a classical sequential least squares programming (SLSQP) routine is leveraged to find the optimal power level corresponding to a predetermined number of candidate vectors; iii) in the last stage, the approximate solution of UC along with generating units power level combination is given. To demonstrate the effectiveness of the presented method, three different systems with 3 generating units, 10 generating units, and 26 generating units were tested for different time periods. In addition, convergence of the hybrid quantum-classical algorithm for select time periods is proven out on IonQ's Forte system.

Aboumrad, Willie [IonQ, Inc]↗

Detection and Perception of Sound by Eagles and Surrogate Raptors

One overarching objective of this program of study was the accumulation of objective, scientifically valid information relating to auditory performance of bald and golden eagles that may be used to guide the development of acoustic alerting/deterrence technologies intended to discourage encroachment into wind energy air spaces. To that end, analyses aimed at the characterization of sensitivity to sound in bald and golden eagles, along with findings in the supra-threshold, dynamic frequency spaces related to response latencies and amplitudes, leads us to conclude that bald, and golden eagles navigate the same basic working auditory space, as in other known and thus far characterized members of the diurnal raptor family. Specifically, bald and golden eagles, along with other raptor species within the group, operate in an auditory space characterized by a frequency band at least four octaves wide and centered on 2 kHz, with an upper frequency limit between 6 and 10 kHz at 80 dB SPL and a lower frequency limit that almost certainly extends below 0.2 kHz. Consequently, we recommend that signal designers use these data as a guideline in efforts to design effective and efficient acoustic alerting/deterrent systems. It is important to note that signal energy broadcast outside of this frequency band at moderate levels will not contribute to the efficacy of a deterrent but will add an unnecessary fraction to the overall acoustic pollution budget. The importance of this consideration is heightened by contemporaneous concerns related to the transmission of noise broadcast by wind energy farms. In addition, based on analyses of data acquired from red-tailed hawks using the same experimental paradigm and data acquisition system, we conclude that auditory function in the red-tailed hawk is sufficiently like that observed in bald and golden eagles to permit its use as a surrogate species. Response waveforms, threshold-frequency curves, and input-output characteristics match those of eagles closely. It should be noted however, that differences in sensitivity and slightly extended high-frequency limits of hearing should be taken into account when extrapolating findings from one species to the others. Although the inclusion of behavioral tests of red-tailed hawks to acoustic stimuli was beyond the scope of this investigation, future efforts to assess response parameters like signal-type preference and habituation rate will further elucidate their suitability to serve as eagle surrogates in behavioral studies; nonetheless, the species in question are well matched with respect to basic auditory performance. A second essential objective of this program of study was the acoustic characterization of a subset of calls comprising the vocal repertoires of bald and golden eagles that may be used to supplement auditory performance findings in the effort to guide the development of acoustic alerting signals. With regard to that objective, the vocal repertoires of both bald and golden eagle species are rich and varied. While similar in spectrographic structure, distinctive differences are also clear. Generally, golden eagles produce some calls with shorter durations, and similar “sounding” calls exhibit distinctively different spectrographic patterns than those of bald eagles. Both species produce calls that contain a wide variety of nonlinear elements that operate to enhance the rich and varied nature of commonly observed vocal products. Comparison of the average power spectra of commonly observed bald and golden eagle calls with threshold-frequency curves leads to the conclusion that call energies fall within the frequency bounds of hearing. Further, the acoustic energy of calls considered in this report tend to fall into overlapping, but different frequency ranges of the acoustic sensitivity curve. This condition may encourage signal designers to vary the frequency content of acoustic deterrence signals in the field. Finally, preliminary observations relating to the tendencies and proclivities of bald eagles to attend to the acoustic landscape lead to the conclusion that eagles monitor their immediate sound environment assiduously. Individuals respond to a variety of natural and synthetic sound signals reliably and, perhaps most relevant in the context of the engineering of acoustic alerting/deterrence technologies, habituation to most sounds considered in this effort was minimal. These preliminary results, while calling for extended behavioral testing, are promising and set the stage for the exportation of behavioral studies into real world scenarios.

17 WIND ENERGY↗

Trajectory Design Employing Convex Optimization for Landing on Irregularly Shaped Asteroids

Mission proposals that land on asteroids are becoming popular. However, in order to have a successful mission the spacecraft must reliably and softly land at the intended landing site. The problem under investigation is how to design a fuel-optimal powered descent trajectory that can be quickly computed on- board the spacecraft, without interaction from ground control. An optimal trajectory designed immediately prior to the descent burn has many advantages. These advantages include the ability to use the actual vehicle starting state as the initial condition in the trajectory design and the ease of updating the landing target site if the original landing site is no longer viable. For long trajectories, the trajectory can be updated periodically by a redesign of the optimal trajectory based on current vehicle conditions to improve the guidance performance. One of the key drivers for being completely autonomous is the infrequent and delayed communication between ground control and the vehicle. Challenges that arise from designing an asteroid powered descent trajectory include complicated nonlinear gravity fields, small rotating bodies and low thrust vehicles. There are two previous studies that form the background to the current investigation. The first set looked in-depth at applying convex optimization to a powered descent trajectory on Mars with promising results.1, 2 This showed that the powered descent equations of motion can be relaxed and formed into a convex optimization problem and that the optimal solution of the relaxed problem is indeed a feasible solution to the original problem. This analysis used a constant gravity field. The second area applied a successive solution process to formulate a second order cone program that designs rendezvous and proximity operations trajectories.3, 4 These trajectories included a Newtonian gravity model. The equivalence of the solutions between the relaxed and the original problem is theoretically established. The proposed solution for designing the asteroid powered descent trajectory is to use convex optimization, a gravity model with higher fidelity than Newtonian, and an iterative solution process to design the fuel optimal trajectory. The solution to the convex optimization problem is the thrust profile, magnitude and direction, that will yield the minimum fuel trajectory for a soft landing at the target site, subject to various mission and operational constraints. The equations of motion are formulated in a rotating coordinate system and includes a high fidelity gravity model. The vehicle's thrust magnitude can vary between maximum and minimum bounds during the burn. Also, constraints are included to ensure that the vehicle does not run out of propellant, or go below the asteroid's surface, and any vehicle pointing requirements. The equations of motion are discretized and propagated with the trapezoidal rule in order to produce equality constraints for the optimization problem. These equality constraints allow the optimization algorithm to solve the entire problem, without including a propagator inside the optimization algorithm.

Pinson, Robin M.↗

Battery asset management with cycle life prognosis

We report Battery Asset Management problem determines the minimum cost replacement schedules for each individual asset in a group of battery assets that operate in parallel. Battery cycle life varies under different operating conditions including temperature, depth of discharge (DOD), charge rate, etc., and a battery deteriorates due to usage, which cannot be handled by current asset management models. This paper presents a new battery asset management methodology where battery cycle life prognosis is integrated with parallel asset management to reduce lifecycle cost of the Battery Energy Storage Systems (BESS). For the battery failure time prognosis, a nonlinear physics-based battery capacity fade model is developed and incorporated in parallel asset management model to update battery capacity over time. Experiment results have shown that the developed battery asset management methodology can be conveniently used to facilitate BESS asset management decision making thereby decreasing asset lifecycle costs.

25 ENERGY STORAGE↗

Inversion of Time-Lapse Seismic Reservoir Monitoring Data Using CycleGAN: A Deep Learning-Based Approach for Estimating Dynamic Reservoir Property Changes

Carbon capture and storage is being pursued globally as a geoengineering measure for reducing the emission of anthropogenic CO 2 the atmosphere. Comprehensive monitoring, verification, and accounting programs must be established for demonstrating the safe storage of injected CO 2 . One of the most commonly deployed monitoring techniques is time-lapse seismic reservoir monitoring (also known as 4-D seismic), which involves comparing 3-D seismic survey data taken at the same study site but over different times. Analyses of 4-D seismic data volumes can help improve the quality of storage reservoir characterization, track the movement of injected CO 2 plume, and identify potential CO 2 spillover/leakage from the storage reservoirblue. However, the derivation of high-resolution CO 2 saturation maps from 4-D seismic data is a highly nonlinear and ill-posed inverse problem, often requiring significant computational effort. In this research, we apply a physics-based deep learning method to facilitate the solution of both the forward and inverse problems in seismic inversion while honoring physical constraints. A cycle generative adversarial neural network (CycleGAN) model is trained to learn the bidirectional functional mappings between the reservoir dynamic property changes and seismic attribute changes, such that both forward and inverse solutions can be obtained efficiently from the trained model. We show that our CycleGAN-based approach not only improves the reliability of 4-D seismic inversion but also expedites the quantitative interpretation. Our deep learning-based workflow is generic and can be readily used for reservoir characterization and reservoir model updates involving the use of 4-D seismic data.

58 GEOSCIENCES↗

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. I. Normal responses

The horizontal angular vestibuloocular reflex (VOR) evoked by high-frequency, high-acceleration rotations was studied in five squirrel monkeys with intact vestibular function. The VOR evoked by steps of acceleration in darkness (3,000 degrees /s(2) reaching a velocity of 150 degrees /s) began after a latency of 7.3 +/- 1.5 ms (mean +/- SD). Gain of the reflex during the acceleration was 14.2 +/- 5.2% greater than that measured once the plateau head velocity had been reached. A polynomial regression was used to analyze the trajectory of the responses to steps of acceleration. A better representation of the data was obtained from a polynomial that included a cubic term in contrast to an exclusively linear fit. For sinusoidal rotations of 0.5-15 Hz with a peak velocity of 20 degrees /s, the VOR gain measured 0.83 +/- 0.06 and did not vary across frequencies or animals. The phase of these responses was close to compensatory except at 15 Hz where a lag of 5.0 +/- 0.9 degrees was noted. The VOR gain did not vary with head velocity at 0.5 Hz but increased with velocity for rotations at frequencies of >/=4 Hz (0. 85 +/- 0.04 at 4 Hz, 20 degrees /s; 1.01 +/- 0.05 at 100 degrees /s, P < 0.0001). No responses to these rotations were noted in two animals that had undergone bilateral labyrinthectomy indicating that inertia of the eye had a negligible effect for these stimuli. We developed a mathematical model of VOR dynamics to account for these findings. The inputs to the reflex come from linear and nonlinear pathways. The linear pathway is responsible for the constant gain across frequencies at peak head velocity of 20 degrees /s and also for the phase lag at higher frequencies being less than that expected based on the reflex delay. The frequency- and velocity-dependent nonlinearity in VOR gain is accounted for by the dynamics of the nonlinear pathway. A transfer function that increases the gain of this pathway with frequency and a term related to the third power of head velocity are used to represent the dynamics of this pathway. This model accounts for the experimental findings and provides a method for interpreting responses to these stimuli after vestibular lesions.

NASA Program Biomedical Research and Countermeasur↗