Functional significance of stiffness in adaptation of multijoint arm movements to stable and unstable dynamics.

This study compared the mechanisms of adaptation to stable and unstable dynamics from the perspective of changes in joint mechanics. Subjects were instructed to make point to point movements in force fields generated by a robotic manipulandum which interacted with the arm in either a stable or an unstable manner. After subjects adjusted to the initial disturbing effects of the force fields they were able to produce normal straight movements to the target. In the case of the stable interaction, subjects modified the joint torques in order to appropriately compensate for the force field. No change in joint torque or endpoint force was required or observed in the case of the unstable interaction. After adaptation, the endpoint stiffness of the arm was measured by applying displacements to the hand in eight different directions midway through the movements. This was compared to the stiffness measured similarly during movements in a null force field. After adaptation, the endpoint stiffness under both the stable and unstable dynamics was modified relative to the null field. Adaptation to unstable dynamics was achieved by selective modification of endpoint stiffness in the direction of the instability. To investigate whether the change in endpoint stiffness could be accounted for by change in joint torque or endpoint force, we estimated the change in stiffness on each trial based on the change in joint torque relative to the null field. For stable dynamics the change in endpoint stiffness was accurately predicted. However, for unstable dynamics the change in endpoint stiffness could not be reproduced. In fact, the predicted endpoint stiffness was similar to that in the null force field. Thus, the change in endpoint stiffness seen after adaptation to stable dynamics was directly related to changes in net joint torque necessary to compensate for the dynamics in contrast to adaptation to unstable dynamics, where a selective change in endpoint stiffness occurred without any modification of net joint torque.

Different mechanisms involved in adaptation to stable and unstable dynamics.

Recently, we demonstrated that humans can learn to make accurate movements in an unstable environment by controlling magnitude, shape, and orientation of the endpoint impedance. Although previous studies of human motor learning suggest that the brain acquires an inverse dynamics model of the novel environment, it is not known whether this control mechanism is operative in unstable environments. We compared learning of multijoint arm movements in a "velocity-dependent force field" (VF), which interacted with the arm in a stable manner, and learning in a "divergent force field" (DF), where the interaction was unstable. The characteristics of error evolution were markedly different in the 2 fields. The direction of trajectory error in the DF alternated to the left and right during the early stage of learning; that is, signed error was inconsistent from movement to movement and could not have guided learning of an inverse dynamics model. This contrasted sharply with trajectory error in the VF, which was initially biased and decayed in a manner that was consistent with rapid feedback error learning. EMG recorded before and after learning in the DF and VF are also consistent with different learning and control mechanisms for adapting to stable and unstable dynamics, that is, inverse dynamics model formation and impedance control. We also investigated adaptation to a rotated DF to examine the interplay between inverse dynamics model formation and impedance control. Our results suggest that an inverse dynamics model can function in parallel with an impedance controller to compensate for consistent perturbing force in unstable environments.

Adaptation to stable and unstable dynamics achieved by combined impedance control and inverse dynamics model.

This study compared adaptation in novel force fields where trajectories were initially either stable or unstable to elucidate the processes of learning novel skills and adapting to new environments. Subjects learned to move in a null force field (NF), which was unexpectedly changed either to a velocity-dependent force field (VF), which resulted in perturbed but stable hand trajectories, or a position-dependent divergent force field (DF), which resulted in unstable trajectories. With practice, subjects learned to compensate for the perturbations produced by both force fields. Adaptation was characterized by an initial increase in the activation of all muscles followed by a gradual reduction. The time course of the increase in activation was correlated with a reduction in hand-path error for the DF but not for the VF. Adaptation to the VF could have been achieved solely by formation of an inverse dynamics model and adaptation to the DF solely by impedance control. However, indices of learning, such as hand-path error, joint torque, and electromyographic activation and deactivation suggest that the CNS combined these processes during adaptation to both force fields. Our results suggest that during the early phase of learning there is an increase in endpoint stiffness that serves to reduce hand-path error and provides additional stability, regardless of whether the dynamics are stable or unstable. We suggest that the motor control system utilizes an inverse dynamics model to learn the mean dynamics and an impedance controller to assist in the formation of the inverse dynamics model and to generate needed stability.

Generalization in Adaptation to Stable and Unstable Dynamics

Humans skillfully manipulate objects and tools despite the inherent instability. In order to succeed at these tasks, the sensorimotor control system must build an internal representation of both the force and mechanical impedance. As it is not practical to either learn or store motor commands for every possible future action, the sensorimotor control system generalizes a control strategy for a range of movements based on learning performed over a set of movements. Here, we introduce a computational model for this learning and generalization, which specifies how to learn feedforward muscle activity in a function of the state space. Specifically, by incorporating co-activation as a function of error into the feedback command, we are able to derive an algorithm from a gradient descent minimization of motion error and effort, subject to maintaining a stability margin. This algorithm can be used to learn to coordinate any of a variety of motor primitives such as force fields, muscle synergies, physical models or artificial neural networks. This model for human learning and generalization is able to adapt to both stable and unstable dynamics, and provides a controller for generating efficient adaptive motor behavior in robots. Simulation results exhibit predictions consistent with all experiments on learning of novel dynamics requiring adaptation of force and impedance, and enable us to re-examine some of the previous interpretations of experiments on generalization. © 2012 Kadiallah et al.

How are internal models of unstable tasks formed?

The results of recent studies suggest that humans can form internal models that they use in a feedforward manner to compensate for both stable and unstable dynamics. To examine how internal models are formed, we performed adaptation experiments in novel dynamics, and measured the endpoint force, trajectory and EMG during learning. Analysis of reflex feedback and change of feedforward commands between consecutive trials suggested a unified model of motor learning, which can coherently unify the learning processes observed in stable and unstable dynamics and reproduce available data on motor learning. To our knowledge, this algorithm, based on the concurrent minimization of (reflex) feedback and muscle activation, is also the first nonlinear adaptive controller able to stabilize unstable dynamics.

Fluid-structure interaction with mean flow: Over-scattering and unstable resonance growth

It is known theoretically [1-3] that infinitely long fluid loaded plates in mean flow exhibit a range of unusual phenomena in the 'long time' limit. These include convective instability, absolute instability and negative energy waves which are destabilized by dissipation. However, structures are necessarily of finite length and may have discontinuities. Moreover, linear instability waves can only grow over a limited number of cycles before non-linear effects become dominant. We have undertaken an analytical and computational study to investigate the response of finite, discontinuous plates to ascertain if these unusual effects might be realized in practice. Analytically, we take a "wave scattering" [2,4] - as opposed to a "modal superposition" [5] - view of the fluttering plate problem. First, we solve for the scattering coefficients of localized plate discontinuities and identify a range of parameter space, well outside the convective instability regime, where over-scattering or amplified reflection/transmission occurs. These are scattering processes that draw energy from the mean flow into the plate. Next, we use the Wiener-Hopf technique to solve for the scattering coefficients from the leading and trailing edges of a baffled plate. Finally, we construct the response of a finite, baffled plate by a superposition of infinite plate propagating waves continuously scattering off the plate ends and solve for the unstable resonance frequencies and temporal growth rates for long plates. We present a comparison between our computational results and the infinite plate theory. In particular, the resonance response of a moderately sized plate is shown to be in excellent agreement with our long plate analytical predictions. Copyright © 2010 by ASME.

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

Lyapunov's stability and operational margin of grid-connected photovoltaic modules

Analyses of photovoltaic power generation based on Lyapunov's theorems are presented. The characteristics of the photovoltaic module and the power conditioning unit are analyzed in order to establish energy functions that assess the stability of solutions and define safe regions of operation. Furthermore, it is shown that grid-connected photovoltaic modules driven at maximum power may become unstable under normal grid transients. In such cases, stability can be maintained by allowing an operational margin defined as the energy difference between the stable and the unstable solutions of the system. Simulations show that modules cope well with grid transients when a sufficiently large margin is used.

The stochastic implicit Euler method - A stable coupling scheme for Monte Carlo burnup calculations

Existing Monte Carlo burnup codes use various schemes to solve the coupled criticality and burnup equations. Previous studies have shown that the coupling schemes of the existing Monte Carlo burnup codes can be numerically unstable. Here we develop the Stochastic Implicit Euler method - a stable and efficient new coupling scheme. The implicit solution is obtained by the stochastic approximation at each time step. Our test calculations demonstrate that the Stochastic Implicit Euler method can provide an accurate solution to problems where the methods in the existing Monte Carlo burnup codes fail. © 2013 Elsevier Ltd. All rights reserved.

Numerically stable Monte Carlo-burnup-thermal hydraulic coupling schemes

This paper presents stochastic implicit coupling method intended for use in Monte-Carlo (MC) based reactor analysis systems that include burnup and thermal hydraulic (TH) feedbacks. Both feedbacks are essential for accurate modeling of advanced reactor designs and analyses of associated fuel cycles. In particular, we investigate the effect of different burnup-TH coupling schemes on the numerical stability and accuracy of coupled MC calculations. First, we present the beginning of time step method which is the most commonly used. The accuracy of this method depends on the time step length and it is only conditionally stable. This work demonstrates that even for relatively short time steps, this method can be numerically unstable. Namely, the spatial distribution of neutronic and thermal hydraulic parameters, such as nuclide densities and temperatures, exhibit oscillatory behavior. To address the numerical stability issue, new implicit stochastic methods are proposed. The methods solve the depletion and TH problems simultaneously and use under-relaxation to speed up convergence. These methods are numerically stable and accurate even for relatively large time steps and require less computation time than the existing methods. © 2013 Elsevier Ltd. All rights reserved.

Calculations and predications of combustion oscillations

The combustion oscillations are the phenomena which we may meet in developing the clean, safe and efficient power and propulsion systems. This paper summarizes authors' systematic work on fuel spray combustion oscillations in the recent years. Combining CFD calculations and stability analysis, a new approach of predicting combustion stabilities was developed. With this approach, detailed flow information and unstable modes can be obtained by CFD and solving perturbation equations, respectively. The results provide the guidance on understanding combustion instability mechanisms and developing the control strategies.

Interfacial immobilization of monoclonal antibody and detection of human prostate-specific antigen.

Antibody orientation and its antigen binding efficiency at interface are of particular interest in many immunoassays and biosensor applications. In this paper, spectroscopic ellipsometry (SE), neutron reflection (NR), and dual polarization interferometry (DPI) have been used to investigate interfacial assembly of the antibody [mouse monoclonal anti-human prostate-specific antigen (anti-hPSA)] at the silicon oxide/water interface and subsequent antigen binding. It was found that the mass density of antibody adsorbed at the interface increased with solution concentration and adsorption time while the antigen binding efficiency showed a steady decline with increasing antibody amount at the interface over the concentration range studied. The amount of antigen bound to the interfacial immobilized antibody reached a maximum when the surface-adsorbed amount of antibody was around 1.5 mg/m(2). This phenomenon is well interpreted by the interfacial structural packing or crowding. NR revealed that the Y-shaped antibody laid flat on the interface at low surface mass density with a thickness around 40 Å, equivalent to the short axial length of the antibody molecule. The loose packing of the antibody within this range resulted in better antigen binding efficiency, while the subsequent increase of surface-adsorbed amount led to the crowding or overlapping of antibody fragments, hence reducing the antigen binding due to the steric hindrance. In situ studies of antigen binding by both NR and DPI demonstrated that the antigen inserted into the antibody layer rather than forming an additional layer on the top. Stability assaying revealed that the antibody immobilized at the silica surface remained stable and active over the monitoring period of 4 months. These results are useful in forming a general understanding of antibody interfacial behavior and particularly relevant to the control of their activity and stability in biosensor development.

Central representation of dynamics when manipulating handheld objects.

To explore the neural mechanisms related to representation of the manipulation dynamics of objects, we performed whole-brain fMRI while subjects balanced an object in stable and highly unstable states and while they balanced a rigid object and a flexible object in the same unstable state, in all cases without vision. In this way, we varied the extent to which an internal model of the manipulation dynamics was required in the moment-to-moment control of the object's orientation. We hypothesized that activity in primary motor cortex would reflect the amount of muscle activation under each condition. In contrast, we hypothesized that cerebellar activity would be more strongly related to the stability and complexity of the manipulation dynamics because the cerebellum has been implicated in internal model-based control. As hypothesized, the dynamics-related activation of the cerebellum was quite different from that of the primary motor cortex. Changes in cerebellar activity were much greater than would have been predicted from differences in muscle activation when the stability and complexity of the manipulation dynamics were contrasted. On the other hand, the activity of the primary motor cortex more closely resembled the mean motor output necessary to execute the task. We also discovered a small region near the anterior edge of the ipsilateral (right) inferior parietal lobule where activity was modulated with the complexity of the manipulation dynamics. We suggest that this is related to imagining the location and motion of an object with complex manipulation dynamics.

Stability enhancement by casing grooves: The importance of stall inception mechanism and solidity

This paper concerns the optimisation of casing grooves and the important influence of stall inception mechanism on groove performance. Installing casing grooves is a well known technique for improving the stable operating range of a compressor, but the wide-spread use of grooves is restricted by the loss of efficiency and flow capacity. In this paper, laboratory tests are used to examine the conditions under which casing treatment can be used to greatest effect. The use of a single casing groove was investigated in a recently published companion paper. The current work extends this to multiple-groove treatments and considers their performance in relation to stall inception mechanisms. Here it is shown that the stall margin gain from multiple grooves is less than the sum of the gains if the grooves were used individually. By contrast, the loss of efficiency is additive as the number of grooves increases. It is then shown that casing grooves give the greatest stall margin improvement when used in a compressor which exhibits spike-type stall inception, while modal activity before stall can dramatically reduce the effectiveness of the grooves. This finding highlights the importance of being able to predict the stall inception mechanism which might occur in a given compressor before and after grooves are added. Some published prediction techniques are therefore examined, but found wanting. Lastly, it is shown that casing grooves can, in some cases, be used to remove rotor blades and produce a more efficient, stable and light-weight rotor. © 2010 by ASME.

Obtaining bifurcation diagrams with a thermoacoustic network model

Linear techniques can predict whether the non-oscillating (steady) state of a thermoacoustic system is stable or unstable. With a sufficiently large impulse, however, a thermoacoustic system can reach a stable oscillating state even when the steady state is also stable. A nonlinear analysis is required to predict the existence of this oscillating state. Continuation methods are often used for this but they are computationally expensive. In this paper, an acoustic network code called LOTAN is used to obtain the steady and the oscillating solutions for a horizontal Rijke tube. The heat release is modelled as a nonlinear function of the mass flow rate. Several test cases from the literature are analysed in order to investigate the effect of various nonlinear terms in the flame model. The results agree well with the literature, showing that LOTAN can be used to map the steady and oscillating solutions as a function of the control parameters. Furthermore, the nature of the bifurcation between steady and oscillating states can be predicted directly from the nonlinear terms inside the flame model. Copyright © 2012 by ASME.