Parallel computing concepts and methods for floquet analysis of helicopter trim and stability

Floquet analysis is widely used for small-order systems (say, order M < 100) to find trim results of control inputs and periodic responses, and stability results of damping levels and frequencies, Presently, however, it is practical neither for design applications nor for comprehensive analysis models that lead to large systems (M > 100); the run time on a sequential computer is simply prohibitive, Accordingly, a massively parallel Floquet analysis is developed with emphasis on large systems, and it is implemented on two SIMD or single-instruction, multiple-data computers with 4096 and 8192 processors, The focus of this development is a parallel shooting method with damped Newton iteration to generate trim results; the Floquet transition matrix (FTM) comes out as a byproduct, The eigenvalues and eigenvectors of the FTM are computed by a parallel QR method, and thereby stability results are generated, For illustration, flap and flap-lag stability of isolated rotors are treated by the parallel analysis and by a corresponding sequential analysis with the conventional shooting and QR methods; linear quasisteady airfoil aerodynamics and a finite-state three-dimensional wake model are used, Computational reliability is quantified by the condition numbers of the Jacobian matrices in Newton iteration, the condition numbers of the eigenvalues and the residual errors of the eigenpairs, and reliability figures are comparable in both the parallel and sequential analyses, Compared to the sequential analysis, the parallel analysis reduces the run time of large systems dramatically, and the reduction increases with increasing system order; this finding offers considerable promise for design and comprehensive-analysis applications.

Surface-ray tracing solution of ellipsoid of revolution for UTD mutual coupling applications

An analytical surface-ray tracing has been carried out for the prolate ellipsoid of revolution using a novel geodesic constant method. This method yields closed form expressions for all the ray-geometric parameters required for the UTD mutual coupling calculations for the antennas located arbitrarily in three dimensions, on the ellipsoid of revolution.

Singularity structure and chaotic behavior of the homopolar disk dynamo

We study in great detail a system of three first-order ordinary differential equations describing a homopolar disk dynamo (HDD). This system displays a large variety of behaviors, both regular and chaotic. Existence of periodic solutions is proved for certain ranges of parameters. Stability criteria for periodic solutions are given. The nonintegrability aspects of the HDD system are studied by investigating analytically the singularity structure of the system in the complex domain. Coexisting attractors (including period-doubling sequence) and coexisting strange attractors appear in some parametric regimes. The gluing of strange attractors and the ungluing of a strange attractor are also shown to occur. A period of bifurcation leading to chaos, not observed for other chaotic systems, is shown to characterize the chaotic behavior in some parametric ranges. The limiting case of the Lorenz system is also studied and is related to HDD.

Time-Resolved Photoluminescence and Microwave Conductivity at Semiconductor Electrodes: Depletion Layer Effects

Analytical expressions which include depletion layer effects on low-injection carrier relaxation are being presented for the first time here. Starting from the continuity equation for the minority carriers, we derive expressions for the output signal pertinent to time-resolved microwave and luminescence experiments. These are valid for the time domain that usually overlaps with the time scales of surface processes, such as charge transfer and trapping. Apart from the usual pulse form of illumination, theoretical expressions pertaining to other forms of illumination such as switch-on and switch-off transient modes, a periodic mode, and a steady state and their various inter-relationships are derived here. The expressions obtained are seen to be generalizations of existing flat-band low-injection results in the Limit of early or initial band bendings. The importance of the depletion layer as an experimental parameter is clearly seen in the limit of larger band bendings wherein it is shown, unlike the flat-band case, to exhibit pure exponential forms of carrier relaxation. Our results are consistent with the main conclusions of the numerical and experimental work published recently. Furthermore, this work provides the actual functional relationships between the applied potential and observed carrier decay. This should enable one to extract the surface kinetic parameters, after deciding on the dominant mode of carrier relaxation at the interface, whether charge transfer or trapping, by studying the potential dependence of the fate of relaxation.

Plate-shaped transformation products in zirconium-base alloys

Plate-shaped products resulting from martensitic, diffusional, and mixed mode transformations in zirconium-base alloys are compared. in the present study. These alloys are particularly suitable for the comparison in view of the fact that the lattice correspondence between the parent beta (bcc) and the product alpha (hcp) or gamma-hydride (fct) phases are remarkably similar for different types of transformations. Crystallographic features such as orientation relations, habit planes, and interface structures associated with these transformations have been compared:, with a view toward examining whether the transformation mechanisms have characteristic imprints on these experimental observables. Martensites exhibiting dislocated lath, internally twinned plate, and self-accommodating three-plate cluster morphologies have been encountered in Zr-2.5Nb alloy. Habit planes corresponding to all these morphologies have been found to be consistent with the predictions based on the invariant plane strain (IFS) criterion. Different morphologies have been found to reflect the manner in which the neighboring martensite variants are assembled. Lattice-invariant shears (LISs) for all these cases have been identified to be either {10 (1) over bar 1}(alpha) ((1) over bar 123)(alpha) slip or twinning on (10 (1) over bar 1)(alpha) planes. Widmanstatten alpha precipitates, forming in a step-quenching treatment, have been shown to have a lath morphology, the alpha/beta interface being decorated with a periodic array of (c + a) dislocations at a spacing of 8 to 10 nm. The line vectors of these dislocations are nearly parallel to the invariant lines. The alpha precipitates, forming in the retained beta phase on aging, exhibit an internally twinned structure with a zigzag habit plane. Average habit planes for the morphologies have been found to lie near the {103}(beta) - {113}(beta) poles, which are close to the specific variant of the {112}(beta) plane, which transforms into a prismatic plane of the type {1 (1) over bar 00}(alpha). The crystallography of the formation of the gamma-hydride phase (fct) from both the alpha and beta phases is seen to match the IFS predictions. While the beta-gamma transformation can be treated approximately as a simple shear on the basal plane involving a change in the stacking sequence, the alpha-gamma transformation call be conceptually broken into a alpha --> beta transformation following the Burgers correspondence and the simple beta-gamma shear process. The active eutectoid decomposition in the Zr-Cu system, beta --> alpha + beta', has been described in terms of cooperative growth of the alpha phase from the beta phase through the Burgers correspondence and of the partially ordered beta' (structurally similar to the equilibrium Zr2Cu phase) through an ordering process. Similarities and differences in crystallographic features of these transformations have been discussed. and the importance of the invariant line vector in deciding the geometry of the corresponding habit planes has been pointed out.

Halogen Bonding in 2,5-Dichloro-1,4-benzoquinone: Insights from Experimental and Theoretical Charge Density Analysis

Experimental charge density distribution in 2, 5-dichloro-1, 4-benzoquinone has been carried out using high resolution X-ray diffraction data at 90 K to quantitatively evaluate the nature of C-Cl center dot center dot center dot O=C halogen bond in molecular crystals. Additionally, the halogen bond is studied from geometrical point of view and the same has been visualized using Hirshfeld surface analysis. The obtained results from experimental charge density analysis are compared with periodic quantum calculations using B3LYP 6-31G(d,p) level of theory. The topological values at bond critical point, three-dimensional static deformation density features and electrostatic potential isosurfaces unequivocally establish the attractive nature of C-Cl center dot center dot center dot O=C halogen bond in crystalline lattice.

Spectral asymptotics of the Helmholtz model in fluid-solid structures

A model representing the vibrations of a coupled fluid-solid structure is considered. This structure consists of a tube bundle immersed in a slightly compressible fluid. Assuming periodic distribution of tubes, this article describes the asymptotic nature of the vibration frequencies when the number of tubes is large. Our investigation shows that classical homogenization of the problem is not sufficient for this purpose. Indeed, our end result proves that the limit spectrum consists of three parts: the macro-part which comes from homogenization, the micro-part and the boundary layer part. The last two components are new. We describe in detail both macro- and micro-parts using the so-called Bloch wave homogenization method. Copyright (C) 1999 John Wiley & Sons, Ltd.

The 2D Coulomb gas on a 1D lattice

The statistical mechanics of a two-dimensional Coulomb gas confined to one dimension is studied, wherein hard core particles move on a ring. Exact self-duality is shown for a version of the sine-Gordon model arising in this context, thereby locating the transition temperature exactly. We present asymptotically exact results for the correlations in the model and characterize the low- and high-temperature phases. Numerical simulations provide support to these renormalization group calculations. Connections with other interesting problems, such as the quantum Brownian motion of a panicle in a periodic potential and impurity problems, are pointed out.

Zero-point entropy in 'spin ice'

Common water ice (ice I-h) is an unusual solid-the oxygen atoms form a periodic structure but the hydrogen atoms are highly disordered due to there being two inequivalent O-H bond lengths'. Pauling showed that the presence of these two bond lengths leads to a macroscopic degeneracy of possible ground states(2,3), such that the system has finite entropy as the temperature tends towards zero. The dynamics associated with this degeneracy are experimentally inaccessible, however, as ice melts and the hydrogen dynamics cannot be studied independently of oxygen motion(4). An analogous system(5) in which this degeneracy can be studied is a magnet with the pyrochlore structure-termed 'spin ice'-where spin orientation plays a similar role to that of the hydrogen position in ice I-h. Here we present specific-heat data for one such system, Dy2Ti2O7, from which we infer a total spin entropy of 0.67Rln2. This is similar to the value, 0.71Rln2, determined for ice I-h, SO confirming the validity of the correspondence. We also find, through application of a magnetic field, behaviour not accessible in water ice-restoration of much of the ground-state entropy and new transitions involving transverse spin degrees of freedom.

A differential geometric method for kinematic analysis of two- and three-degree-of-freedom rigid body motions

In this paper, we present a novel differential geometric characterization of two- and three-degree-of-freedom rigid body kinematics, using a metric defined on dual vectors. The instantaneous angular and linear velocities of a rigid body are expressed as a dual velocity vector, and dual inner product is defined on this dual vector, resulting in a positive semi-definite and symmetric dual matrix. We show that the maximum and minimum magnitude of the dual velocity vector, for a unit speed motion, can be obtained as eigenvalues of this dual matrix. Furthermore, we show that the tip of the dual velocity vector lies on a dual ellipse for a two-degree-of-freedom motion and on a dual ellipsoid for a three-degree-of-freedom motion. In this manner, the velocity distribution of a rigid body can be studied algebraically in terms of the eigenvalues of a dual matrix or geometrically with the dual ellipse and ellipsoid. The second-order properties of the two- and three-degree-of-freedom motions of a rigid body are also obtained from the derivatives of the elements of the dual matrix. This results in a definition of the geodesic motion of a rigid body. The theoretical results are illustrated with the help of a spatial 2R and a parallel three-degree-of-freedom manipulator.

Driven Heisenberg magnets: Nonequilibrium criticality, spatiotemporal chaos and control

We drive a d-dimensional Heisenberg magnet using an anisotropic current. The continuum Langevin equation is analysed using a dynamical renormalization group and numerical simulations. We discover a rich steady-state phase diagram, including a critical point in a new nonequilibrium universality class, and a spatiotemporally chaotic phase. The latter may be controlled in a robust manner to target spatially periodic steady states with helical order.

Minimizing buffer requirements under rate-optimal schedule in regular dataflow networks

Large-grain synchronous dataflow graphs or multi-rate graphs have the distinct feature that the nodes of the dataflow graph fire at different rates. Such multi-rate large-grain dataflow graphs have been widely regarded as a powerful programming model for DSP applications. In this paper we propose a method to minimize buffer storage requirement in constructing rate-optimal compile-time (MBRO) schedules for multi-rate dataflow graphs. We demonstrate that the constraints to minimize buffer storage while executing at the optimal computation rate (i.e. the maximum possible computation rate without storage constraints) can be formulated as a unified linear programming problem in our framework. A novel feature of our method is that in constructing the rate-optimal schedule, it directly minimizes the memory requirement by choosing the schedule time of nodes appropriately. Lastly, a new circular-arc interval graph coloring algorithm has been proposed to further reduce the memory requirement by allowing buffer sharing among the arcs of the multi-rate dataflow graph. We have constructed an experimental testbed which implements our MBRO scheduling algorithm as well as (i) the widely used periodic admissible parallel schedules (also known as block schedules) proposed by Lee and Messerschmitt (IEEE Transactions on Computers, vol. 36, no. 1, 1987, pp. 24-35), (ii) the optimal scheduling buffer allocation (OSBA) algorithm of Ning and Gao (Conference Record of the Twentieth Annual ACM SIGPLAN-SIGACT Symposium on Principles of Programming Languages, Charleston, SC, Jan. 10-13, 1993, pp. 29-42), and (iii) the multi-rate software pipelining (MRSP) algorithm (Govindarajan and Gao, in Proceedings of the 1993 International Conference on Application Specific Array Processors, Venice, Italy, Oct. 25-27, 1993, pp. 77-88). Schedules generated for a number of random dataflow graphs and for a set of DSP application programs using the different scheduling methods are compared. The experimental results have demonstrated a significant improvement (10-20%) in buffer requirements for the MBRO schedules compared to the schedules generated by the other three methods, without sacrificing the computation rate. The MBRO method also gives a 20% average improvement in computation rate compared to Lee's Block scheduling method.

Persistent passive hopping and juggling is possible even with plastic collisions

We describe simple one-dimensional models of passive (no energy input, no control), generally dissipative, vertical hopping and one-ball juggling. The central observation is that internal passive system motions can conspire to eliminate collisions in these systems. For hopping, two point masses are connected by a spring and the lower mass has inelastic collisions with the ground. For juggling, a lower point-mass hand is connected by a spring to the ground and an upper point-mass ball is caught with an inelastic collision and then re-thrown into gravitational free flight. The two systems have identical dynamics. Despite inelastic collisions between non-zero masses, these systems have special symmetric energy-conserving periodic motions where the collision is at zero relative velocity. Additionally, these special periodic motions have a non-zero sized, one-sided region of attraction on the higher-energy side. For either very large or very small mass ratios, the one-sided region of attraction is large. These results persist for mildly non-linear springs and non-constant gravity. Although non-collisional damping destroys the periodic motions, small energy injection makes the periodic motions stable, with a two-sided region of attraction. The existence of such special energy conserving solutions for hopping and juggling points to possibly useful strategies for both animals and robots. The lossless motions are demonstrated with a table-top experiment.

Quantum spin models with exact dimer ground states

Inspired by the exact solution of the Majumdar-Ghosh model, a family of one-dimensional, translationally invariant spin Hamiltonians is constructed. The exchange coupling in these models is antiferromagnetic, and decreases linearly with the separation between the spins. The coupling becomes identically zero beyond a certain distance. It is rigorously proved that the dimer configuration is an exact, superstable ground-state configuration of all the members of the family on a periodic chain. The ground state is twofold degenerate, and there exists an energy gap above the ground state. The Majumdar-Ghosh Hamiltonian with a twofold degenerate dimer ground state is just the first member of the family. The scheme of construction is generalized to two and three dimensions, and illustrated with the help of some concrete examples. The first member in two dimensions is the Shastry-Sutherland model. Many of these models have exponentially degenerate, exact dimer ground states.

Ratchet for energy transport between identical reservoirs

A one-dimensional periodic array of elastically colliding hard points, with a noncentrosymmetric unit cell, connected at its two ends to identical but nonthermal energy reservoirs, is shown to carry a sustained unidirectional energy current.