A demonstration of artificial dissipation effects in some finite volume solution methods for the Navier-Stokes equations

The artificial dissipation effects in some solutions obtained with a Navier-Stokes flow solver are demonstrated. The solvers were used to calculate the flow of an artificially dissipative fluid, which is a fluid having dissipative properties which arise entirely from the solution method itself. This was done by setting the viscosity and heat conduction coefficients in the Navier-Stokes solvers to zero everywhere inside the flow, while at the same time applying the usual no-slip and thermal conducting boundary conditions at solid boundaries. An artificially dissipative flow solution is found where the dissipation depends entirely on the solver itself. If the difference between the solutions obtained with the viscosity and thermal conductivity set to zero and their correct values is small, it is clear that the artificial dissipation is dominating and the solutions are unreliable.

The acquisition of movement skills: Practice enhances the dynamic stability of bimanual coordination

During bimanual movements, two relatively stable inherent patterns of coordination (in-phase and anti-phase) are displayed (e.g., Kelso, Am. J. Physiol. 246 (1984) R1000). Recent research has shown that new patterns of coordination can be learned. For example, following practice a 90 degrees out-of-phase pattern can emerge as an additional, relatively stable, state (e.g., Zanone & Kelso, J. Exp. Psychol.: Human Performance and Perception 18 (1992) 403). On this basis, it has been concluded that practice leads to the evolution and stabilisation of the newly learned pattern and that this process of learning changes the entire attractor layout of the dynamic system. A general feature of such research has been to observe the changes of the targeted pattern's stability characteristics during training at a single movement frequency. The present study was designed to examine how practice affects the maintenance of a coordinated pattern as the movement frequency is scaled. Eleven volunteers were asked to perform a bimanual forearm pronation-supination task. Time to transition onset was used as an index of the subjects' ability to maintain two symmetrically opposite coordinated patterns (target task - 90 degrees out-of-phase - transfer task - 270 degrees out-of-phase). Their ability to maintain the target task and the transfer task were examined again after five practice sessions each consisting of 15 trials of only the 90 degrees out-of-phase pattern. Concurrent performance feedback (a Lissajous figure) was available to the participants during each practice trial. A comparison of the time to transition onset showed that the target task was more stable after practice (p = 0.025). These changes were still observed one week (p = 0.05) and two months (p = 0.075) after the practice period. Changes in the stability of the transfer task were not observed until two months after practice (p = 0.025). Notably, following practice, transitions from the 90 degrees pattern were generally to the anti-phase (180 degrees) pattern, whereas, transitions from the 270 degrees pattern were to the 90 degrees pattern. These results suggest that practice does improve the stability of a 90 degrees pattern, and that such improvements are transferable to the performance of the unpractised 270 degrees pattern. In addition, the anti-phase pattern remained more stable than the practised 90 degrees pattern throughout. (C) 2001 Elsevier Science B.V. All rights reserved.

Controlling the complex Lorenz equations by modulation

We demonstrate that a system obeying the complex Lorenz equations in the deep chaotic regime can be controlled to periodic behavior by applying a modulation to the pump parameter. For arbitrary modulation frequency and amplitude there is no obvious simplification of the dynamics. However, we find that there are numerous windows where the chaotic system has been controlled to different periodic behaviors. The widths of these windows in parameter space are narrow, and the positions are related to the ratio of the modulation frequency of the pump to the average pulsation frequency of the output variable. These results are in good agreement with observations previously made in a far-infrared laser system.

Subcycling first- and second-order generalizations of the trapezoidal rule

Subcycling algorithms which employ multiple timesteps have been previously proposed for explicit direct integration of first- and second-order systems of equations arising in finite element analysis, as well as for integration using explicit/implicit partitions of a model. The author has recently extended this work to implicit/implicit multi-timestep partitions of both first- and second-order systems. In this paper, improved algorithms for multi-timestep implicit integration are introduced, that overcome some weaknesses of those proposed previously. In particular, in the second-order case, improved stability is obtained. Some of the energy conservation properties of the Newmark family of algorithms are shown to be preserved in the new multi-timestep extensions of the Newmark method. In the first-order case, the generalized trapezoidal rule is extended to multiple timesteps, in a simple way that permits an implicit/implicit partition. Explicit special cases of the present algorithms exist. These are compared to algorithms proposed previously. (C) 1998 John Wiley & Sons, Ltd.

Heat flow for the Yang-Mills-Higgs field and the Hermitian Yang-Mills-Higgs metric

For a parameter lambda > 0, we study a type of vortex equations, which generalize the well-known Hermitian-Einstein equation, for a connection A and a section phi of a holomorphic vector bundle E over a Kahler manifold X. We establish a global existence of smooth solutions to heat flow for a self-dual Yang-Mills-Higgs field on E. Assuming the lambda -stability of (E, phi), we prove the existence of the Hermitian Yang-Mills-Higgs metric on the holomorphic bundle E by studying the limiting behaviour of the gauge flow.

Interaction of directional, neuromuscular and egocentric constraints on the stability of preferred bimanual coordination patterns

We investigated how the relative direction of limb movements in external space (iso- and non-isodirectionality), muscular constraints (the relative timing of homologous muscle activation) and the egocentric frame of reference (moving simultaneously toward/away the longitudinal axis of the body) contribute to the stability of coordinated movements. In the first experiment, we attempted to determine the respective stability of isodirectional and non-isodirectional movements in between-persons coordination. In a second experiment, we determined the effect of the relative direction in external space, and of muscular constraints, on pattern stability during a within-person bimanual coordination task. In the third experiment we dissociated the effects on pattern stability of the muscular constraints, relative direction and egocentric frame of reference. The results showed that (1) simultaneous activation of homologous muscles resulted in more stable performance than simultaneous activation of non-homologous muscles during within-subject coordination, and that (2) isodirectional movements were more stable than non-isodirectional movements during between-persons coordination, confirming the role of the relative direction of the moving limbs in the stability of bimanual coordination. Moreover, the egocentric constraint was to some extent found distinguishable from the effect of the relative direction of the moving limbs in external space, and from the effect of the relative timing of muscle activation. In summary, the present study showed that relative direction of the moving limbs in external space and muscular constraints may interact either to stabilize or destabilize coordination patterns. (C) 2003 Published by Elsevier B.V.

Dynamic stability of laminated FGM plates based on higher-order shear deformation theory

This paper conducts a dynamic stability analysis of symmetrically laminated FGM rectangular plates with general out-of-plane supporting conditions, subjected to a uniaxial periodic in-plane load and undergoing uniform temperature change. Theoretical formulations are based on Reddy's third-order shear deformation plate theory, and account for the temperature dependence of material properties. A semi-analytical Galerkin-differential quadrature approach is employed to convert the governing equations into a linear system of Mathieu-Hill equations from which the boundary points on the unstable regions are determined by Bolotin's method. Free vibration and bifurcation buckling are also discussed as subset problems. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with FGM layers made of silicon nitride and stainless steel. The influences of various parameters such as material composition, layer thickness ratio, temperature change, static load level, boundary constraints on the dynamic stability, buckling and vibration frequencies are examined in detail through parametric studies.

Thermal, chemical, and enzymatic stability of the cyclotide kalata B1: The importance of the cyclic cystine knot

The cyclotides constitute a recently discovered family of plant-derived peptides that have the unusual features of a head-to-tail cyclized backbone and a cystine knot core. These features are thought to contribute to their exceptional stability, as qualitatively observed during experiments aimed at sequencing and characterizing early members of the family. However, to date there has been no quantitative study of the thermal, chemical, or enzymatic stability of the cyclotides. In this study, we demonstrate the stability of the prototypic cyclotide kalata B1 to the chaotropic agents 6 M guanidine hydrochloride (GdHCl) and 8 M urea, to temperatures approaching boiling, to acid, and following incubation with a range of proteases, conditions under which most proteins readily unfold. NMR spectroscopy was used to demonstrate the thermal stability, while fluorescence and circular dichroism were used to monitor the chemical stability. Several variants of kalata B1 were also examined, including kalata 132, which has five amino acid substitutions from B1, two acyclic permutants in which the backbone was broken but the cystine knot was retained, and a two-disulfide bond mutant. Together, these allowed determinations of the relative roles of the cystine knot and the circular backbone on the stability of the cyclotides. Addition of a denaturant to kalata B1 or an acyclic permutant did not cause unfolding, but the two-disulfide derivative was less stable, despite having a similar three-dimensional structure. It appears that the cystine knot is more important than the circular backbone in the chemical stability of the cyclotides. Furthermore, the cystine knot of the cyclotides is more stable than those in similar-sized molecules, judging by a comparison with the conotoxin PVIIA. There was no evidence for enzymatic digestion of native kalata B1 as monitored by LC-MS, but the reduced form was susceptible to proteolysis by trypsin, endoproteinase Glu-C, and thermolysin. Fluorescence spectra of kalata B1 in the presence of dithiothreitol, a reducing agent, showed a marked increase in intensity thought to be due to removal of the quenching effect on the Trp residue by the neighboring Cys5-Cys17 disulfide bond. In general, the reduced peptides were significantly more susceptible to chemical or enzymatic breakdown than the oxidized species.

Robust stability of sequential multi-input multi-output quantitative feedback theory designs

This paper re-examines the stability of multi-input multi-output (MIMO) control systems designed using sequential MIMO quantitative feedback theory (QFT). In order to establish the results, recursive design equations for the SISO equivalent plants employed in a sequential MIMO QFT design are established. The equations apply to sequential MIMO QFT designs in both the direct plant domain, which employs the elements of plant in the design, and the inverse plant domain, which employs the elements of the plant inverse in the design. Stability theorems that employ necessary and sufficient conditions for robust closed-loop internal stability are developed for sequential MIMO QFT designs in both domains. The theorems and design equations facilitate less conservative designs and improved design transparency.

Scaling laws for the critical rupture thickness or common thin films

Despite decades of experimental and theoretical investigation on thin films, considerable uncertainty exists in the prediction of their critical rupture thickness. According to the spontaneous rupture mechanism, common thin films become unstable when capillary waves. at the interfaces begin to grow. In a horizontal film with symmetry at the midplane. unstable waves from adjacent interfaces grow towards the center of the film. As the film drains and becomes thinner, unstable waves osculate and cause the film to rupture, Uncertainty sterns from a number of sources including the theories used to predict film drainage and corrugation growth dynamics. In the early studies, (lie linear stability of small amplitude waves was investigated in the Context of the quasi-static approximation in which the dynamics of wave growth and film thinning are separated. The zeroth order wave growth equation of Vrij predicts faster wave growth rates than the first order equation derived by Sharma and Ruckenstein. It has been demonstrated in an accompanying paper that film drainage rates and times measured by numerous investigations are bounded by the predictions of the Reynolds equation and the more recent theory of Manev, Tsekov, and Radoev. Solutions to combinations of these equations yield simple scaling laws which should bound the critical rupture thickness of foam and emulsion films, In this paper, critical thickness measurements reported in the literature are compared to predictions from the bounding scaling equations and it is shown that the retarded Hamaker constants derived from approximate Lifshitz theory underestimate the critical thickness of foam and emulsion films, The non-retarded Hamaker constant more adequately bounds the critical thickness measurements over the entire range of film radii reported in the literature. This result reinforces observations made by other independent researchers that interfacial interactions in flexible liquid films are not adequately represented by the retarded Hamaker constant obtained from Lifshitz theory and that the interactions become significant at much greater separations than previously thought. (c) 2005 Elsevier B.V. All rights reserved.

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

Factors influencing the stability of cyclotides: Proteins with a circular backbone and cystine knot motif

Cyclotides are a large family of mini-proteins that have the distinguishing features of a head-to-tail cyclised backbone and a cystine knot formed by six conserved cysteine residues. They are present in plants from the Rubiaceae, Violaceae and Cucurbitaceae families. The unique structural features of the cyclotides make them extremely resistant to chemical, thermal and proteolytic degradation. In this article we review recent Studies from our laboratory that dissect the role of the individual structural elements in defining the stability of cyclotides. The resistance of cyclotides to chemical and proteolytic degradation is in large part due to the cystine knot, whereas the thermal stability is I composite of several features including the cystine knot, the cyclic backbone and the hydrogen bonding network. A range of biological activities of cyclotides is critically dependent oil the presence of the cyclic backbone.