Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “non linear dynamics”

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.

537 records · Page 30

System Identification of X-33 Neural Network

Modern flight control research has improved spacecraft survivability as its goal. To this end we need to have a failure detection system on board. In case the spacecraft is performing imperfectly, reconfiguration of control is needed. For that purpose we need to have parameter identification of spacecraft dynamics. Parameter identification of a system is called system identification. We treat the system as a black box which receives some inputs that lead to some outputs. The question is: what kind of parameters for a particular black box can correlate the observed inputs and outputs? Can these parameters help us to predict the outputs for a new given set of inputs? This is the basic problem of system identification. The X33 was supposed to have the onboard capability of evaluating the current performance and if needed to take the corrective measures to adapt to desired performance. The X33 is comprised of both rocket and aircraft vehicle design characteristics and requires, in general, analytical methods for evaluating its flight performance. Its flight consists of four phases: ascent, transition, entry and TAEM (Terminal Area Energy Management). It spends about 200 seconds in ascent phase, reaching an altitude of about 180,000 feet and a speed of about 10 to 15 Mach. During the transition phase which lasts only about 30 seconds, its altitude may increase to about 190,000 feet but its speed is reduced to about 9 Mach. At the beginning of this phase, the Main Engine is Cut Off (MECO) and the control is reconfigured with the help of aerosurfaces (four elevons, two flaps and two rudders) and reaction control system (RCS). The entry phase brings down the altitude of X33 to about 90,000 feet and its speed to about Mach 3. It spends about 250 seconds in this phase. Main engine is still cut off and the vehicle is controlled by complex maneuvers of aerosurfaces. The last phase TAEM lasts for about 450 seconds and the altitude and speed, both are reduced to zero. The present attempt, as a start, focuses only on the entry phase. Since the main engine remains cut off in this phase, there is no thrust acting on the system. This considerably simplifies the equations of motion. We introduce another simplification by assuming the system to be linear after some non-linearities are removed analytically from our consideration. Under these assumptions, the problem could be solved by Classical Statistics by employing the least sum of squares approach. Instead we chose to use the Neural Network method. This method has many advantages. It is modern, more efficient, can be adapted to work even when the assumptions are diluted. In fact, Neural Networks try to model the human brain and are capable of pattern recognition.

Aggarwal, Shiv↗

Effect of acute exposure to hypergravity (GX vs. GZ) on dynamic cerebral autoregulation

We examined the effects of 30 min of exposure to either +3GX (front-to-back) or +GZ (head-to-foot) centrifugation on cerebrovascular responses to 80 degrees head-up tilt (HUT) in 14 healthy individuals. Both before and after +3 GX or +3 GZ centrifugation, eye-level blood pressure (BP(eye)), end tidal PCO2 (PET(CO2)), mean cerebral flow velocity (CFV) in the middle cerebral artery (transcranial Doppler ultrasound), cerebral vascular resistance (CVR), and dynamic cerebral autoregulatory gain (GAIN) were measured with subjects in the supine position and during subsequent 80 degrees HUT for 30 min. Mean BP(eye) decreased with HUT in both the GX (n = 7) and GZ (n = 7) groups (P < 0.001), with the decrease being greater after centrifugation only in the GZ group (P < 0.05). PET(CO2) also decreased with HUT in both groups (P < 0.01), but the absolute level of decrease was unaffected by centrifugation. CFV decreased during HUT more significantly after centrifugation than before centrifugation in both groups (P < 0.02). However, these greater decreases were not associated with greater increases in CVR. In the supine position after centrifugation compared with before centrifugation, GAIN increased in both groups (P < 0.05, suggesting an autoregulatory deficit), with the change being correlated to a measure of otolith function (the linear vestibulo-ocular reflex) in the GX group (r = 0.76, P < 0.05) but not in the GZ group (r = 0.24, P = 0.60). However, GAIN was subsequently restored to precentrifugation levels during postcentrifugation HUT (i.e., as BP(eye) decreased), suggesting that both types of centrifugation resulted in a leftward shift of the cerebral autoregulation curve. We speculate that this leftward shift may have been due to vestibular activation (especially during +GX) or potentially to an adaptation to reduced cerebral perfusion pressure during +GZ.

NASA Discipline Neuroscience↗

Ventilatory responses to dynamic exercise elicited by intramuscular sensors

PURPOSE: Eight subjects, aged 27.0+/-1.6 yr, performed incremental workload cycling to investigate the contribution of skeletal muscle mechano- and metaboreceptors to ventilatory control during dynamic exercise. METHODS: Each subject performed four bouts of exercise: exercise with no intervention (CON); exercise with bilateral thigh cuffs inflated to 90 mm Hg (CUFF); exercise with application of lower-body positive pressure (LBPP) to 45 torr (PP); and exercise with 90 mm Hg thigh cuff inflation and 45 torr LBPP (CUFF+PP). Ventilatory responses and pulmonary gas exchange variables were collected breath-by-breath with concomitant measurement of leg intramuscular pressure. RESULTS: Ventilation (VE) was significantly elevated from CON during PP and CUFF+PP at workloads corresponding to > or = 60% CON peak oxygen uptake (VO2peak) and during CUFF at workloads > or = 80% CON VO2peak, P < 0.05. The VO2 at which ventilatory threshold occurred was significantly reduced from CON (2.17+/-0.28 L x min(-1)) to 1.60+/-0.19 L x min(-1), 1.45+/-0.15 L x min(-1), and 1.15+/-0.11 L x min(-1) during CUFF, PP, and CUFF+PP, respectively. The slope of the linear regression describing the VE/CO2 output relationship was increased from CON by approximately 22% during CUFF, 40% during PP, and 41% during CUFF+PP. CONCLUSIONS: As intramuscular pressure was significantly elevated immediately upon application of LBPP during PP and CUFF+PP without a concomitant increase in VE, it seems unlikely that LBPP-induced increases in VE can be attributed to activation of the mechanoreflex. These findings suggest that LBPP-induced reductions in perfusion pressure and decreases in venous outflow resulting from inflation of bilateral thigh cuffs may generate a metabolite sensitive intramuscular ventilatory stimulus.

Non-NASA Center↗

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗

Buoyancy Effects on Flow Structure and Instability of Low-Density Gas Jets

A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.

Pasumarthi, Kasyap Sriramachandra↗

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗

Experimental tests of a superposition hypothesis to explain the relationship between the vestibuloocular reflex and smooth pursuit during horizontal combined eye-head tracking in humans

1. We used a modeling approach to test the hypothesis that, in humans, the smooth pursuit (SP) system provides the primary signal for cancelling the vestibuloocular reflex (VOR) during combined eye-head tracking (CEHT) of a target moving smoothly in the horizontal plane. Separate models for SP and the VOR were developed. The optimal values of parameters of the two models were calculated using measured responses of four subjects to trials of SP and the visually enhanced VOR. After optimal parameter values were specified, each model generated waveforms that accurately reflected the subjects' responses to SP and vestibular stimuli. The models were then combined into a CEHT model wherein the final eye movement command signal was generated as the linear summation of the signals from the SP and VOR pathways. 2. The SP-VOR superposition hypothesis was tested using two types of CEHT stimuli, both of which involved passive rotation of subjects in a vestibular chair. The first stimulus consisted of a "chair brake" or sudden stop of the subject's head during CEHT; the visual target continued to move. The second stimulus consisted of a sudden change from the visually enhanced VOR to CEHT ("delayed target onset" paradigm); as the vestibular chair rotated past the angular position of the stationary visual stimulus, the latter started to move in synchrony with the chair. Data collected during experiments that employed these stimuli were compared quantitatively with predictions made by the CEHT model. 3. During CEHT, when the chair was suddenly and unexpectedly stopped, the eye promptly began to move in the orbit to track the moving target. Initially, gaze velocity did not completely match target velocity, however; this finally occurred approximately 100 ms after the brake onset. The model did predict the prompt onset of eye-in-orbit motion after the brake, but it did not predict that gaze velocity would initially be only approximately 70% of target velocity. One possible explanation for this discrepancy is that VOR gain can be dynamically modulated and, during sustained CEHT, it may assume a lower value. Consequently, during CEHT, a smaller-amplitude SP signal would be needed to cancel the lower-gain VOR. This reduction of the SP signal could account for the attenuated tracking response observed immediately after the brake. We found evidence for the dynamic modulation of VOR gain by noting differences in responses to the onset and offset of head rotation in trials of the visually enhanced VOR.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

A dual-input nonlinear system analysis of autonomic modulation of heart rate

Linear analyses of fluctuations in heart rate and other hemodynamic variables have been used to elucidate cardiovascular regulatory mechanisms. The role of nonlinear contributions to fluctuations in hemodynamic variables has not been fully explored. This paper presents a nonlinear system analysis of the effect of fluctuations in instantaneous lung volume (ILV) and arterial blood pressure (ABP) on heart rate (HR) fluctuations. To successfully employ a nonlinear analysis based on the Laguerre expansion technique (LET), we introduce an efficient procedure for broadening the spectral content of the ILV and ABP inputs to the model by adding white noise. Results from computer simulations demonstrate the effectiveness of broadening the spectral band of input signals to obtain consistent and stable kernel estimates with the use of the LET. Without broadening the band of the ILV and ABP inputs, the LET did not provide stable kernel estimates. Moreover, we extend the LET to the case of multiple inputs in order to accommodate the analysis of the combined effect of ILV and ABP effect on heart rate. Analyzes of data based on the second-order Volterra-Wiener model reveal an important contribution of the second-order kernels to the description of the effect of lung volume and arterial blood pressure on heart rate. Furthermore, physiological effects of the autonomic blocking agents propranolol and atropine on changes in the first- and second-order kernels are also discussed.

NASA Discipline Regulatory Physiology↗

Comparative transduction mechanisms of hair cells in the bullfrog utriculus. II. Sensitivity and response dynamics to hair bundle displacement

1. Hair cells in whole-mount in vitro preparations of the utricular macula of the bullfrog (Rana catesbeiana) were selected according to their macular location and hair bundle morphology. The sensitivity and response dynamics of selected hair cells to natural stimulation were examined by recording their voltage responses to step and sinusoidal hair bundle displacements applied to their longest stereocilia. 2. The voltage responses of 31 hair cells to sinusoidal hair bundle displacements were characterized by their gains and phases, taken with respect to peak hair bundle displacement. The gains of Type B and Type C cells at both 0.5 and 5.0 Hz were markedly lower than those of Type F and Type E cells. Phases, with the exception of Type C cells, lagged hair bundle displacement at 0.5 Hz. Type C cells had phase leads of 25-40 degrees. At 5.0 Hz, response phases in all cells were phase lagged with respect to those at 0.5 Hz. Type C cells had larger gains and smaller phase leads at 5.0 Hz than at 0.5 Hz, suggesting the presence of low-frequency adaptation. 3. Displacement-response curves, derived from the voltage responses to 5.0-Hz sinusoids, were sigmoidal in shape and asymmetrical, with the depolarizing response having a greater magnitude and saturating less abruptly than the hyperpolarizing response. When normalized to their largest displacement the linear ranges of these curves varied from < 0.5 to 1.25 microns and were largest in Type B and smallest in Type F and Type E cells. Sensitivity, defined as the slope of the normalized displacement-response curve, was inversely correlated with linear range. 4. The contribution of geometric factors associated with the hair bundle to linear range and sensitivity were predicted from realistic models of utricular hair bundles created using morphological data obtained from light and electron microscopy. Three factors, including 1) the inverse ratio of the lengths of the kinocilium and longest stereocilia, representing the lever arm between kinociliary and stereociliary displacement; 2) tip link extension/linear displacement, largely a function of stereociliary height and separation; and 3) stereociliary number, an estimate of the number of transduction channels, were considered in this analysis. The first of these factors was quantitatively more important than the latter two factors and their total contribution was largest in Type B and Type C cells. Theoretical models were also used to calculate the relation between rotary and linear displacement.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

Role of vestibular information in initiation of rapid postural responses

Patients with bilateral vestibular loss have difficulty maintaining balance without stepping when standing in tandem, on compliant surfaces, across narrow beams, or on one foot, especially with eyes closed. Normal individuals (with no sensory impairment) maintain balance in these tasks by employing quick, active hip rotation (a "hip strategy"). The absence of a hip strategy in vestibular patients responding to translations of a short support surface has previously been taken as evidence that the use of hip strategy requires an intact vestibular system. However, many tasks requiring hip strategy alter one or a combination of important system characteristics, such as initial state of the body (tandem stance), dynamics (compliant surfaces), or biomechanical limits of stability (narrow beams). Therefore, the balance deficit in these tasks may result from a failure to account for these support surface alterations when planning and executing sensorimotor responses. In this study, we tested the hypothesis that vestibular information is critical to trigger a hip strategy even on an unaltered support surface, which imposes no changes on the system characteristics. We recorded the postural responses of vestibular patients and control subjects with eyes closed to rearward support surface translations of varying velocity, in erect stance on a firm, flat surface. Subjects were instructed to maintain balance without stepping, if possible. Faster translation velocities (25 cm/s or more) produced a consistent pattern of early hip torque (first 400 ms) in control subjects (i.e., a hip strategy). Most of the patients with bilateral vestibular loss responded to the same translation velocities with similar torques. Contrary to our hypothesis, we conclude that vestibular function is not necessary to trigger a hip strategy. We postulate, therefore, that the balance deficit previously observed in vestibular patients during postural tasks that elicit a hip strategy may have been due to the sensorimotor consequences of the system alterations imposed by the postural tasks used in those studies. Preliminary results from two younger patients who lost vestibular function as infants indicate that age, duration of vestibular loss, and/or the timing of the loss may also be factors that can influence the use of hip strategy as a rapid postural response.

NASA Discipline Neuroscience↗

Motor adaptation to Coriolis force perturbations of reaching movements: endpoint but not trajectory adaptation transfers to the nonexposed arm

1. Reaching movements made in a rotating room generate Coriolis forces that are directly proportional to the cross product of the room's angular velocity and the arm's linear velocity. Such Coriolis forces are inertial forces not involving mechanical contact with the arm. 2. We measured the trajectories of arm movements made in darkness to a visual target that was extinguished at the onset of each reach. Prerotation subjects pointed with both the right and left arms in alternating sets of eight movements. During rotation at 10 rpm, the subjects reached only with the right arm. Postrotation, the subjects pointed with the left and right arms, starting with the left, in alternating sets of eight movements. 3. The initial perrotary reaching movements of the right arm were highly deviated both in movement path and endpoint relative to the prerotation reaches of the right arm. With additional movements, subjects rapidly regained straight movement paths and accurate endpoints despite the absence of visual or tactile feedback about reaching accuracy. The initial postrotation reaches of the left arm followed straight paths to the wrong endpoint. The initial postrotation reaches of the right arm had paths with mirror image curvature to the initial perrotation reaches of the right arm but went to the correct endpoint. 4. These observations are inconsistent with current equilibrium point models of movement control. Such theories predict accurate reaches under our experimental conditions. Our observations further show independent implementation of movement and posture, as evidenced by transfer of endpoint adaptation to the nonexposed arm without transfer of path adaptation. Endpoint control may occur at a relatively central stage that represents general constraints such as gravitoinertial force background or egocentric direction relative to both arms, and control of path may occur at a more peripheral stage that represents moments of inertia and muscle dynamics unique to each limb. 5. Endpoint and path adaptation occur despite the absence both of mechanical contact cues about the perturbing force and visual or tactile cues about movement accuracy. These findings point to the importance of muscle spindle signals, monitoring of motor commands, and possibly joint and tendon receptors in a detailed trajectory monitoring process. Muscle spindle primary and secondary afferent signals may differentially influence adaptation of movement shape and endpoint, respectively.

Non-NASA Center↗

A management information system to study space diets

A management information system (MIS), including a database management system (DBMS) and a decision support system (DSS), was developed to dynamically analyze the variable nutritional content of foods grown and prepared in an Advanced Life Support System (ALSS) such as required for long-duration space missions. The DBMS was designed around the known nutritional content of a list of candidate crops and their prepared foods. The DSS was designed to determine the composition of the daily crew diet based on crop and nutritional information stored in the DBMS. Each of the selected food items was assumed to be harvested from a yet-to-be designed ALSS biomass production subsystem and further prepared in accompanying food preparation subsystems. The developed DBMS allows for the analysis of the nutrient composition of a sample 20-day diet for future Advanced Life Support missions and is able to determine the required quantities of food needed to satisfy the crew's daily consumption. In addition, based on published crop growth rates, the DBMS was able to calculate the required size of the biomass production area needed to satisfy the daily food requirements for the crew. Results from this study can be used to help design future ALSS for which the integration of various subsystems (e.g., biomass production, food preparation and consumption, and waste processing) is paramount for the success of the mission.

NASA Discipline Life Support Systems↗

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

Spatial orientation of the vestibular system

1. A simplified three-dimensional state space model of visual vestibular interaction was formulated. Matrix and dynamical system operators representing coupling from the semicircular canals and the visual system to the velocity storage integrator were incorporated into the model. 2. It was postulated that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix while another rotates the pitch and roll eigenvectors with the head, while maintaining the yaw axis eigenvector approximately spatially invariant. Using this representation, the response characteristics of the pitch, roll, and yaw eye velocity were obtained in terms of the eigenvalues and associated eigenvectors. 3. Using OKAN data obtained from monkeys and comparing to the model predictions, the eigenvalues and eigenvectors of the system matrix were identified as a function of tilt to the side or of tilt to the prone positions, using a modification of the Marquardt algorithm. The yaw eigenvector for right-side-down tilt and for downward pitch cross-coupling was approximately 30 degrees from the spatial vertical. For the prone position, the eigenvector was computed to be approximately 20 degrees relative to the spatial vertical. For both side-down and prone positions, oblique OKN induced along eigenvector directions generated OKAN which decayed to zero along a straight line with approximately a single time constant. This was verified by a spectral analysis of the residual sequence about the straight line fit to the decaying data. The residual sequence was associated with a narrow autocorrelation function and a wide power spectrum. 4. Parameters found using the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced to represent the direct pathway and the coupling of the visual system to the integrator. Model simulations predicted the behavior of yaw and pitch OKN and OKAN when the animal was upright, as well as the cross-coupling in the tilted position. The trajectories in velocity space were also accurately simulated. 5. There were similarities between the monkey eigenvectors and human perception of the spatial vertical. For side-down tilts and downward eye velocity cross-coupling, there was only an Aubert (A) effect. For upward eye velocity cross-coupling there were both Muller (E) and Aubert (A) effects. The mean of the eigenvectors for upward and downward eye velocities overlay human 1 x g perceptual data.(ABSTRACT TRUNCATED AT 400 WORDS).

NASA Discipline Neuroscience↗

Three-dimensional organization of otolith-ocular reflexes in rhesus monkeys. I. Linear acceleration responses during off-vertical axis rotation

1. The dynamic properties of otolith-ocular reflexes elicited by sinusoidal linear acceleration along the three cardinal head axes were studied during off-vertical axis rotations in rhesus monkeys. As the head rotates in space at constant velocity about an off-vertical axis, otolith-ocular reflexes are elicited in response to the sinusoidally varying linear acceleration (gravity) components along the interaural, nasooccipital, or vertical head axis. Because the frequency of these sinusoidal stimuli is proportional to the velocity of rotation, rotation at low and moderately fast speeds allows the study of the mid-and low-frequency dynamics of these otolith-ocular reflexes. 2. Animals were rotated in complete darkness in the yaw, pitch, and roll planes at velocities ranging between 7.4 and 184 degrees/s. Accordingly, otolith-ocular reflexes (manifested as sinusoidal modulations in eye position and/or slow-phase eye velocity) were quantitatively studied for stimulus frequencies ranging between 0.02 and 0.51 Hz. During yaw and roll rotation, torsional, vertical, and horizontal slow-phase eye velocity was sinusoidally modulated as a function of head position. The amplitudes of these responses were symmetric for rotations in opposite directions. In contrast, mainly vertical slow-phase eye velocity was modulated during pitch rotation. This modulation was asymmetric for rotations in opposite direction. 3. Each of these response components in a given rotation plane could be associated with an otolith-ocular response vector whose sensitivity, temporal phase, and spatial orientation were estimated on the basis of the amplitude and phase of sinusoidal modulations during both directions of rotation. Based on this analysis, which was performed either for slow-phase eye velocity alone or for total eye excursion (including both slow and fast eye movements), two distinct response patterns were observed: 1) response vectors with pronounced dynamics and spatial/temporal properties that could be characterized as the low-frequency range of "translational" otolith-ocular reflexes; and 2) response vectors associated with an eye position modulation in phase with head position ("tilt" otolith-ocular reflexes). 4. The responses associated with two otolith-ocular vectors with pronounced dynamics consisted of horizontal eye movements evoked as a function of gravity along the interaural axis and vertical eye movements elicited as a function of gravity along the vertical head axis. Both responses were characterized by a slow-phase eye velocity sensitivity that increased three- to five-fold and large phase changes of approximately 100-180 degrees between 0.02 and 0.51 Hz. These dynamic properties could suggest nontraditional temporal processing in utriculoocular and sacculoocular pathways, possibly involving spatiotemporal otolith-ocular interactions. 5. The two otolith-ocular vectors associated with eye position responses in phase with head position (tilt otolith-ocular reflexes) consisted of torsional eye movements in response to gravity along the interaural axis, and vertical eye movements in response to gravity along the nasooccipital head axis. These otolith-ocular responses did not result from an otolithic effect on slow eye movements alone. Particularly at high frequencies (i.e., high speed rotations), saccades were responsible for most of the modulation of torsional and vertical eye position, which was relatively large (on average +/- 8-10 degrees/g) and remained independent of frequency. Such reflex dynamics can be simulated by a direct coupling of primary otolith afferent inputs to the oculomotor plant. (ABSTRACT TRUNCATED).

NASA Discipline Neuroscience↗