Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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 865 records · Page 48

Forced Patterns of Sea Level Rise in the Community Earth System Model Large Ensemble From 1920 to 2100

To provide context for observed sea level rise, the forced response (FR) in dynamic sea level (DSL) during the 20th and 21st centuries is examined in the Community Earth System Model Large Ensemble (LE). After accounting for simulation drift, which in the LE is sizable, the DSL FR is found to be complex, both in space and time. Its evolving character is suggested to arise from both the diversity and transient evolution of climate forcing agents and the slow adjustment timescales of the intermediate and deep oceans. Nonetheless, various intervals of spatially coherent change simulated for the recent past and near future are identified, and their characteristics and associated driving mechanisms are identified and discussed. The mid- to late-20th century DSL FR is characterized by an hemispherically asymmetric pattern of change, with depressed rates of rise in the northern oceans. There is also a dipole of change in the Southern Ocean caused by changes in near-surface zonal winds. Through the late 20th and early 21st centuries, a different pattern of rise emerges with elevated rates in the tropics and depressed rates at high latitudes. Zonal and interbasin variations characterize both intervals and involve the pattern and depth of anomalous ocean heat content storage and spatial contrasts in the expansion coefficient tied mainly to base state temperature, with greater rates of rise in warm regions per unit warming. Finally, the relative roles of surface flux and ocean convergence anomalies are examined.

54 ENVIRONMENTAL SCIENCES↗

Analysis of hourglass instabilities and control in underintegrated finite element methods

Belytschko et al. (1981, 1984) has developed stabilization methods for the treatment of underintegrated FEM problems; these methods involve the computation of an underintegrated stiffness matrix, which is rank-deficient, and the addition of a stabilization matrix which effectively eliminates the spurious modes. An attempt is presently made to give this a priori stabilization method a mathematical means of support. Attention is also given to an a posteriority stabilization method for hourglass control, in which an approximate solution of the underintegrated system is obtained and then subjected to a special projection in order to eliminate the hourglass modes. A proof is obtained for the convergence of this stabilized underintegrated approximation to the exact solution of a model problem at almost the same rate (as the mesh is refined) as the fully integrated solutions.

Jacquotte, O.-P.↗

The solution of linear systems of equations with a structural analysis code on the NAS CRAY-2

Two methods for solving linear systems of equations on the NAS Cray-2 are described. One is a direct method; the other is an iterative method. Both methods exploit the architecture of the Cray-2, particularly the vectorization, and are aimed at structural analysis applications. To demonstrate and evaluate the methods, they were installed in a finite element structural analysis code denoted the Computational Structural Mechanics (CSM) Testbed. A description of the techniques used to integrate the two solvers into the Testbed is given. Storage schemes, memory requirements, operation counts, and reformatting procedures are discussed. Finally, results from the new methods are compared with results from the initial Testbed sparse Choleski equation solver for three structural analysis problems. The new direct solvers described achieve the highest computational rates of the methods compared. The new iterative methods are not able to achieve as high computation rates as the vectorized direct solvers but are best for well conditioned problems which require fewer iterations to converge to the solution.

Poole, Eugene L.↗

Improved collisional excitation rates for interstellar water

Theoretical rate constants among the lowest 45 para and 45 ortho rotational levels of water in collisions with He atoms have been calculated for temperatures between 20 and 2000 K using a recently improved theoretical interaction potential. These values are about 30-40 percent larger than those reported previously but relative sizes of different state-to-state rates have not changed significantly. Successive improvements to the theoretical description of this system now appear to have converged.

Green, Sheldon↗

Effect of metamorphism on isolated olivine grains in CO3 chondrites

The presence of a metamorphic sequence in the CO3 chondrite group has been shown previously to result in changes in properties of chondrule silicates. However, the role of isolated olivine grains during metamorphism of these chondrites has not been addressed. Isolated olivine grains in two metamorphosed CO3 chondrites, Lance and Isna, have been investigated in this study in order to assess the compositional properties of isolated olivine grains that may be attributable to metamorphism. Compositional changes in isolated olivines with increasing petrologic subtype are very similar to changes in chondrule olivines in the same chondrites. Olivine compositions from all occurrences (chondrules, isolated grains, and matrix) converge with increasing petrologic subtype. The degree of equilibration of minor elements is qualitatively related to the diffusion rate of each element in olivine, suggesting that diffusion-controlled processes are the most important processes responsible for compositional changes within the metamorphic sequence. The data are consistent with metamorphism taking place in a closed system on the CO3 chondrite parent body. Fe-poor olivine grains in metamorphosed chondrites are characterized by an Fe-rich rim, which is the result of diffusion of Fe into the grains from Fe-rich matrix. In some instances, 'complex', Fe-rich rims have been identified, which appear to have originated as igneous overgrowths and subsequently to have been overprinted by diffusion processes during metamorphism. Processes experienced by CO3 chondrites are more similar to those experienced by the ordinary chondrites than to those encountered by other carbonaceous chondrites, such as the CV3 group.

Jones, Rhian H.↗

Coherent Detection of High-Rate Optical PPM Signals

A method of coherent detection of high-rate pulse-position modulation (PPM) on a received laser beam has been conceived as a means of reducing the deleterious effects of noise and atmospheric turbulence in free-space optical communication using focal-plane detector array technologies. In comparison with a receiver based on direct detection of the intensity modulation of a PPM signal, a receiver based on the present method of coherent detection performs well at much higher background levels. In principle, the coherent-detection receiver can exhibit quantum-limited performance despite atmospheric turbulence. The key components of such a receiver include standard receiver optics, a laser that serves as a local oscillator, a focal-plane array of photodetectors, and a signal-processing and data-acquisition assembly needed to sample the focal-plane fields and reconstruct the pulsed signal prior to detection. The received PPM-modulated laser beam and the local-oscillator beam are focused onto the photodetector array, where they are mixed in the detection process. The two lasers are of the same or nearly the same frequency. If the two lasers are of different frequencies, then the coherent detection process is characterized as heterodyne and, using traditional heterodyne-detection terminology, the difference between the two laser frequencies is denoted the intermediate frequency (IF). If the two laser beams are of the same frequency and remain aligned in phase, then the coherent detection process is characterized as homodyne (essentially, heterodyne detection at zero IF). As a result of the inherent squaring operation of each photodetector, the output current includes an IF component that contains the signal modulation. The amplitude of the IF component is proportional to the product of the local-oscillator signal amplitude and the PPM signal amplitude. Hence, by using a sufficiently strong local-oscillator signal, one can make the PPM-modulated IF signal strong enough to overcome thermal noise in the receiver circuits: this is what makes it possible to achieve near-quantum-limited detection in the presence of strong background. Following quantum-limited coherent detection, the outputs of the individual photodetectors are automatically aligned in phase by use of one or more adaptive array compensation algorithms [e.g., the least-mean-square (LMS) algorithm]. Then the outputs are combined and the resulting signal is processed to extract the high-rate information, as though the PPM signal were received by a single photodetector. In a continuing series of experiments to test this method (see Fig. 1), the local oscillator has a wavelength of 1,064 nm, and another laser is used as a signal transmitter at a slightly different wavelength to establish an IF of about 6 MHz. There are 16 photodetectors in a 4 4 focal-plane array; the detector outputs are digitized at a sampling rate of 25 MHz, and the signals in digital form are combined by use of the LMS algorithm. Convergence of the adaptive combining algorithm in the presence of simulated atmospheric turbulence for optical PPM signals has already been demonstrated in the laboratory; the combined output is shown in Fig. 2(a), and Fig. 2(b) shows the behavior of the phase of the combining weights as a function of time (or samples). We observe that the phase of the weights has a sawtooth shape due to the continuously changing phase in the down-converted output, which is not exactly at zero frequency. Detailed performance analysis of this coherent free-space optical communication system in the presence of simulated atmospheric turbulence is currently under way.

Vilnrotter, Victor↗

Alaska Crustal Deformation: Finite Element Modeling Constrained by Geologic and Very Long Baseline Interferometry Data

We compute crustal motions in Alaska by calculating the finite element solution for an elastic spherical shell problem. The method we use allows the finite element mesh to include faults and very long baseline interferometry (VLBI) baseline rates of change. Boundary conditions include Pacific-North American (PA-NA) plate motions. The solution is constrained by the oblique orientation of the Fairweather-Queen Charlotte strike-slip faults relative to the PA-NA relative motion direction and the oblique orientation from normal convergence of the eastern Aleutian trench fault systems, as well as strike-shp motion along the Denali and Totschunda fault systems. We explore the effects that a range of fault slip constraints and weighting of VLBI rates of change has on the solution. This allows us to test the motion on faults, such as the Denali fault, where there are conflicting reports on its present-day slip rate. We find a pattern of displacements which produce fault motions generally consistent with geologic observations. The motion of the continuum has the general pattern of radial movement of crust to the NE away from the Fairweather-Queen Charlotte fault systems in SE Alaska and Canada. This pattern of crustal motion is absorbed across the Mackenzie Mountains in NW Canada, with strike-slip motion constrained along the Denali and Tintina fault systems. In south central Alaska and the Alaska forearc oblique convergence at the eastern Aleutian trench and the strike-shp motion of the Denali fault system produce a counterclockwise pattern of motion which is partially absorbed along the Contact and related fault systems in southern Alaska and is partially extruded into the Bering Sea and into the forearc parallel the Aleutian trench from the Alaska Peninsula westward. Rates of motion and fault slip are small in western and northern Alaska, but the motions we compute are consistent with the senses of strike-slip motion inferred geologically along the Kaltag, Kobuk Trench, and Thompson Creek faults and with the normal faulting observed in NW Alaska near Nome. The nonrigid behavior of our finite element solution produces patterns of motion that would not have been expected from rigid block models: strike-slip faults can exist in a continuum that has motion mostly perpendicular to their strikes, and faults can exhibit along-strike differences in magnitudes and directions.

Lundgren, Paul↗

Marshall Convergent Spray Formulation Improvement for High Temperatures

The Marshall Convergent Coating-1 (MCC-1) formulation was produced in the 1990s, and uses a standard bisphenol A epoxy resin system with a triamine accelerator. With the increasing heat rates forecast for the next generation of vehicles, higher-temperature sprayable coatings are needed. This work substitutes the low-temperature epoxy resins used in the MCC-1 coating with epoxy phenolic, epoxy novalac, or resorcinolinic resins (higher carbon content), which will produce a higher char yield upon exposure to high heat and increased glass transition temperature. High-temperature filler materials, such as granular cork and glass ecospheres, are also incorporated as part of the convergent spray process, but other sacrificial (ablative) materials are possible. In addition, the use of polyhedral oligomeric silsesquioxanes (POSS) nanoparticle hybrids will increase both reinforcement aspects and contribute to creating a tougher silacious char, which will reduce recession at higher heat rates. Use of expanding epoxy resin (lightweight MCC) systems are also useful in that they reduce system weight, have greater insulative properties, and a decrease in application times can be realized.

Scarpa, Jack↗

Distributed Inference with Sparse and Quantized Communication

Here, we consider the problem of distributed inference where agents in a network observe a stream of private signals generated by an unknown state, and aim to uniquely identify this state from a finite set of hypotheses. We focus on scenarios where communication between agents is costly, and takes place over channels with finite bandwidth. To reduce the frequency of communication, we develop a novel event-triggered distributed learning rule that is based on the principle of diffusing low beliefs on each false hypothesis. Building on this principle, we design a trigger condition under which an agent broadcasts only those components of its belief vector that have adequate innovation, to only those neighbors that require such information. We prove that our rule guarantees convergence to the true state exponentially fast almost surely despite sparse communication, and that it has the potential to significantly reduce information flow from uninformative agents to informative agents. Next, to deal with finite-precision communication channels, we propose a distributed learning rule that leverages the idea of adaptive quantization. We show that by sequentially refining the range of the quantizers, every agent can learn the truth exponentially fast almost surely, while using just 1 bit to encode its belief on each hypothesis. For both our proposed algorithms, we rigorously characterize the trade-offs between communication-efficiency and the learning rate.

42 ENGINEERING↗

Application of axisymmetric analogue for calculating heating in three-dimensional flows

A rapid, approximate method has been developed for calculating the heating rates on three-dimensional vehicles such as the Space Shuttle Orbiter and other advanced reentry configurations. The method is based on the axisymmetric analogue for three-dimensional boundary layers. It uses information obtained from a three-dimensional inviscid flowfield solution, such as HALIS, to calculate inviscid surface streamlines along which approximate heating rates are calculated independent of what happens along other streamlines. Three-dimensional effects are included through the metric coefficient that describes the divergence or convergence of streamlines. Boundary-layer edge properties are obtained from the inviscid flowfield solution by interpolating in the inviscid flowfield at a distance equal to the boundary-layer thickness away from the wall. This accounts, approximately, for the variable boundary-layer edge entropy. Using this method, heating calculations can be made along a typical streamline in a few seconds. This method has been used to accurately predict heating rates for simple shapes such as a spherically blunted cone and more complex shapes such as the Shuttle Orbiter for a variety of wind-tunnel and flight conditions. A unique feature of the method is its ability to accurately predict heating rates on the Shuttle Orbiter wing.

Hamilton, H. H.↗

COBE ground segment gyro calibration

Discussed here is the calibration of the scale factors and rate biases for the Cosmic Background Explorer (COBE) spacecraft gyroscopes, with the emphasis on the adaptation for COBE of an algorithm previously developed for the Solar Maximum Mission. Detailed choice of parameters, convergence, verification, and use of the algorithm in an environment where the reference attitudes are determined form the Sun, Earth, and star observations (via the Diffuse Infrared Background Experiment (DIRBE) are considered. Results of some recent experiments are given. These include tests where the gyro rate data are corrected for the effect of the gyro baseplate temperature on the spacecraft electronics.

Freedman, I.↗

Convergence in simulating global soil organic carbon by structurally different models after data assimilation

Abstract Current biogeochemical models produce carbon–climate feedback projections with large uncertainties, often attributed to their structural differences when simulating soil organic carbon (SOC) dynamics worldwide. However, choices of model parameter values that quantify the strength and represent properties of different soil carbon cycle processes could also contribute to model simulation uncertainties. Here, we demonstrate the critical role of using common observational data in reducing model uncertainty in estimates of global SOC storage. Two structurally different models featuring distinctive carbon pools, decomposition kinetics, and carbon transfer pathways simulate opposite global SOC distributions with their customary parameter values yet converge to similar results after being informed by the same global SOC database using a data assimilation approach. The converged spatial SOC simulations result from similar simulations in key model components such as carbon transfer efficiency, baseline decomposition rate, and environmental effects on carbon fluxes by these two models after data assimilation. Moreover, data assimilation results suggest equally effective simulations of SOC using models following either first‐order or Michaelis–Menten kinetics at the global scale. Nevertheless, a wider range of data with high‐quality control and assurance are needed to further constrain SOC dynamics simulations and reduce unconstrained parameters. New sets of data, such as microbial genomics‐function relationships, may also suggest novel structures to account for in future model development. Overall, our results highlight the importance of observational data in informing model development and constraining model predictions.

54 ENVIRONMENTAL SCIENCES↗

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗

Investigating On-Orbit Attitude Determination Anomalies for the Solar Dynamics Observatory Mission

The Solar Dynamics Observatory (SDO) was launched on February 11, 2010 from Kennedy Space Center on an Atlas V launch vehicle into a geosynchronous transfer orbit. SDO carries a suite of three scientific instruments, whose observations are intended to promote a more complete understanding of the Sun and its effects on the Earth's environment. After a successful launch, separation, and initial Sun acquisition, the launch and flight operations teams dove into a commissioning campaign that included, among other things, checkout and calibration of the fine attitude sensors and checkout of the Kalman filter (KF) and the spacecraft s inertial pointing and science control modes. In addition, initial calibration of the science instruments was also accomplished. During that process of KF and controller checkout, several interesting observations were noticed and investigated. The SDO fine attitude sensors consist of one Adcole Digital Sun Sensor (DSS), two Galileo Avionica (GA) quaternion-output Star Trackers (STs), and three Kearfott Two-Axis Rate Assemblies (hereafter called inertial reference units, or IRUs). Initial checkout of the fine attitude sensors indicated that all sensors appeared to be functioning properly. Initial calibration maneuvers were planned and executed to update scale factors, drift rate biases, and alignments of the IRUs. After updating the IRU parameters, the KF was initialized and quickly reached convergence. Over the next few hours, it became apparent that there was an oscillation in the sensor residuals and the KF estimation of the IRU bias. A concentrated investigation ensued to determine the cause of the oscillations, their effect on mission requirements, and how to mitigate them. The ensuing analysis determined that the oscillations seen were, in fact, due to an oscillation in the IRU biases. The low frequencies of the oscillations passed through the KF, were well within the controller bandwidth, and therefore the spacecraft was actually following the oscillating biases, resulting in movement of the spacecraft on the order of plus or minus 20 arcsec. Though this level of error met the ACS attitude knowledge requirement of [35, 70, 70] arcsec, 3 sigma, the desire of the ACS and instrument teams was to remove as much of the oscillation as possible. The Kearfott IRUs have an internal temperature controller, designed to maintain the IRU temperature at a constant temperature of approximately 70 C, thus minimizing the change in the bias drift and scale factors of the mechanical gyros. During ground testing of the observatory, it was discovered that the 83-Hz control cycle of the IRU heaters put a tremendous amount of stress on the spacecraft battery. Analysis by the power systems team indicated that the constant charge/discharge on the battery due to the IRU thermal control cycle could potentially limit the life of the battery. After much analysis, the decision was made not to run the internal IRU heaters. Analysis of on orbit data revealed that the oscillations in the IRU bias had a connection to the temperature of the IRU; changes in IRU temperature resulted in changes in the amplitude and period of the IRU biases. Several mitigating solutions were investigated, the result of which was to tune the KF with larger IRU noise assumptions which allows the KF to follow and correct for the time-varying IRU biases.

Vess, Melissa F.↗

Reliability enhancement of Navier-Stokes codes through convergence acceleration

Methods for enhancing the reliability of Navier-Stokes computer codes through improving convergence characteristics are presented. The improving of these characteristics decreases the likelihood of code unreliability and user interventions in a design environment. The problem referred to as a 'stiffness' in the governing equations for propulsion-related flowfields is investigated, particularly in regard to common sources of equation stiffness that lead to convergence degradation of CFD algorithms. Von Neumann stability theory is employed as a tool to study the convergence difficulties involved. Based on the stability results, improved algorithms are devised to ensure efficient convergence in different situations. A number of test cases are considered to confirm a correlation between stability theory and numerical convergence. The examples of turbulent and reacting flow are presented, and a generalized form of the preconditioning matrix is derived to handle these problems, i.e., the problems involving additional differential equations for describing the transport of turbulent kinetic energy, dissipation rate and chemical species. Algorithms for unsteady computations are considered. The extension of the preconditioning techniques and algorithms derived for Navier-Stokes computations to three-dimensional flow problems is discussed. New methods to accelerate the convergence of iterative schemes for the numerical integration of systems of partial differential equtions are developed, with a special emphasis on the acceleration of convergence on highly clustered grids.

Merkle, Charles L.↗

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vestibular convergence patterns in vestibular nuclei neurons of alert primates

Sensory signal convergence is a fundamental and important aspect of brain function. Such convergence may often involve complex multidimensional interactions as those proposed for the processing of otolith and semicircular canal (SCC) information for the detection of translational head movements and the effective discrimination from physically congruent gravity signals. In the present study, we have examined the responses of primate rostral vestibular nuclei (VN) neurons that do not exhibit any eye movement-related activity using 0.5-Hz translational and three-dimensional (3D) rotational motion. Three distinct neural populations were identified. Approximately one-fourth of the cells exclusively encoded rotational movements (canal-only neurons) and were unresponsive to translation. The canal-only central neurons encoded head rotation in SCC coordinates, exhibited little orthogonal canal convergence, and were characterized with significantly higher sensitivities to rotation as compared to primary SCC afferents. Another fourth of the neurons modulated their firing rates during translation (otolith-only cells). During rotations, these neurons only responded when the axis of rotation was earth-horizontal and the head was changing orientation relative to gravity. The remaining one-half of VN neurons were sensitive to both rotations and translations (otolith + canal neurons). Unlike primary otolith afferents, however, central neurons often exhibited significant spatiotemporal (noncosine) tuning properties and a wide variety of response dynamics to translation. To characterize the pattern of SCC inputs to otolith + canal neurons, their rotational maximum sensitivity vectors were computed using exclusively responses during earth-vertical axis rotations (EVA). Maximum sensitivity vectors were distributed throughout the 3D space, suggesting strong convergence from multiple SCCs. These neurons were also tested with earth-horizontal axis rotations (EHA), which would activate both vertical canals and otolith organs. However, the recorded responses could not be predicted from a linear combination of EVA rotational and translational responses. In contrast, one-third of the neurons responded similarly during EVA and EHA rotations, although a significant response modulation was present during translation. Thus this subpopulation of otolith + canal cells, which included neurons with either high- or low-pass dynamics to translation, appear to selectively ignore the component of otolith-selective activation that is due to changes in the orientation of the head relative to gravity. Thus contrary to primary otolith afferents and otolith-only central neurons that respond equivalently to tilts relative to gravity and translational movements, approximately one-third of the otolith + canal cells seem to encode a true estimate of the translational component of the imposed passive head and body movement.

Non-NASA Center↗