Inverse Sensitivity Analysis of Singular Solutions of FRF matrix in Structural System Identification

The problem of structural damage detection based on measured frequency response functions of the structure in its damaged and undamaged states is considered. A novel procedure that is based on inverse sensitivity of the singular solutions of the system FRF matrix is proposed. The treatment of possibly ill-conditioned set of equations via regularization scheme and questions on spatial incompleteness of measurements are considered. The application of the method in dealing with systems with repeated natural frequencies and (or) packets of closely spaced modes is demonstrated. The relationship between the proposed method and the methods based on inverse sensitivity of eigensolutions and frequency response functions is noted. The numerical examples on a 5-degree of freedom system, a one span free-free beam and a spatially periodic multi-span beam demonstrate the efficacy of the proposed method and its superior performance vis-a-vis methods based on inverse eigensensitivity.

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.

Parametrized actor-critic algorithms for finite-horizon MDPs

Due to their non-stationarity, finite-horizon Markov decision processes (FH-MDPs) have one probability transition matrix per stage. Thus the curse of dimensionality affects FH-MDPs more severely than infinite-horizon MDPs. We propose two parametrized 'actor-critic' algorithms to compute optimal policies for FH-MDPs. Both algorithms use the two-timescale stochastic approximation technique, thus simultaneously performing gradient search in the parametrized policy space (the 'actor') on a slower timescale and learning the policy gradient (the 'critic') via a faster recursion. This is in contrast to methods where critic recursions learn the cost-to-go proper. We show w.p 1 convergence to a set with the necessary condition for constrained optima. The proposed parameterization is for FHMDPs with compact action sets, although certain exceptions can be handled. Further, a third algorithm for stochastic control of stopping time processes is presented. We explain why current policy evaluation methods do not work as critic to the proposed actor recursion. Simulation results from flow-control in communication networks attest to the performance advantages of all three algorithms.

A study of symbolic processing and computational aspects in helicopter dynamics

Even research models of helicopter dynamics often lead to a large number of equations of motion with periodic coefficients; and Floquet theory is a widely used mathematical tool for dynamic analysis. Presently, three approaches are used in generating the equations of motion. These are (1) general-purpose symbolic processors such as REDUCE and MACSYMA, (2) a special-purpose symbolic processor, DEHIM (Dynamic Equations for Helicopter Interpretive Models), and (3) completely numerical approaches. In this paper, comparative aspects of the first two purely algebraic approaches are studied by applying REDUCE and DEHIM to the same set of problems. These problems range from a linear model with one degree of freedom to a mildly non-linear multi-bladed rotor model with several degrees of freedom. Further, computational issues in applying Floquet theory are also studied, which refer to (1) the equilibrium solution for periodic forced response together with the transition matrix for perturbations about that response and (2) a small number of eigenvalues and eigenvectors of the unsymmetric transition matrix. The study showed the following: (1) compared to REDUCE, DEHIM is far more portable and economical, but it is also less user-friendly, particularly during learning phases; (2) the problems of finding the periodic response and eigenvalues are well conditioned.

Automatic identification of modal damping from Floquet analysis

A practical method is proposed to identify the mode associated with the frequency part of the eigenvalue of the Floquet transition matrix (FTM). From the FTM eigenvector, which contains the states and their derivatives, the ratio of the derivative and the state corresponding to the largest component is computed. The method exploits the fact that the imaginary part of this (complex) ratio closely approximates the frequency of the mode. It also lends itself well to automation and has been tested over a large number of FTMs of order as high as 250.

Magnetization properties of some quantum spin ladders

The experimental realization of various spin ladder systems has prompted their detailed theoretical investigations. Hen we study the evolution of ground-state magnetization with an external magnetic field for two different antiferromagnetic systems: a three-legged spin-1/2 ladder, and a two-legged spin-1/2 ladder with an additional diagonal interaction. The finite system density-matrix renormalization-group method is employed for numerical studies of the three-chain system, and an effective low-energy Hamiltonian is used in the limit of strong interchain coupling to study the two- and three-chain systems. The three-chain system has a magnetization plateau at one-third of the saturation magnetization. The two-chain system has a plateau at zero magnetization due to a gap above the singlet ground state. It also has a plateau at half of the saturation magnetization for a certain range of values of the couplings. We study the regions of transitions between plateaus numerically and analytically, and find that they are described, at first order in a strong-coupling expansion, by an XXZ spin-1/2 chain in a magnetic field; the second-order terms give corrections to the XXZ model, We also study numerically some low-temperature properties of the three-chain system, such as the magnetization, magnetic susceptibility and specific heat. [S0163-1829(99)303001-5].

Markov chains, R-trivial monoids and representation theory

We develop a general theory of Markov chains realizable as random walks on R-trivial monoids. It provides explicit and simple formulas for the eigenvalues of the transition matrix, for multiplicities of the eigenvalues via Mobius inversion along a lattice, a condition for diagonalizability of the transition matrix and some techniques for bounding the mixing time. In addition, we discuss several examples, such as Toom-Tsetlin models, an exchange walk for finite Coxeter groups, as well as examples previously studied by the authors, such as nonabelian sandpile models and the promotion Markov chain on posets. Many of these examples can be viewed as random walks on quotients of free tree monoids, a new class of monoids whose combinatorics we develop.

A finite variant of the Toom Model

We present results for a finite variant of the one-dimensional Toom model with closed boundaries. We show that the steady state distribution is not of product form, but is nonetheless simple. In particular, we give explicit formulas for the densities and some nearest neighbour correlation functions. We also give exact results for eigenvalues and multiplicities of the transition matrix using the theory of R-trivial monoids in joint work with A. Schilling, B. Steinberg and N. M. Thiery.

NMR relaxation study of disorder in the mixed system BP x GPI(1-x) and the low temperature transition from classical to quantum rotation

1H NMR spin-lattice relaxation time (T1) studies have been carried out in the temperature range 100 K to 4 K, at two Larmor frequencies 11.4 and 23.3 MHz, in the mixed system of betaine phosphate and glycine phosphite (BPxGPI(1-x)), to study the effects of disorder on the proton group dynamics. Analysis of T1 data indicates the presence of a number of inequivalent methyl groups and a gradual transition from classical reorientations to quantum tunneling rotations. At lower temperatures, microstructural disorder in the local environments of the methyl groups, result in a distribution in the activation energy (Ea) and the torsional energy gap (E01). For certain values of x, the magnetisation recovery shows biexponential behaviour at lower temperatures.

Sol → gel transition with rings in a finite polycondensing system in solution

We report a theoretical formulation for the mean cluster size distribution in a finite polycondensing system. Expressions for the mean number of n-mers with j bonds ( nj) are developed. Numerical calculations show that while the non-cyclic molecules make the dominant contribution to the small clusters, the large clusters are dominated by cyclic structures. The number of particles in ringless chains, n n,n-1, decays monotonically with n at all extents of reaction, but n n becomes bimodal near the gel point. We also find that the solvent plays an important role in the cluster size distribution.

Study of linear and nonlinear optical properties of dendrimers using density matrix renormalization group method

We have used the density matrix renormalization group (DMRG) method to study the linear and nonlinear optical responses of first generation nitrogen based dendrimers with donor acceptor groups. We have employed Pariser–Parr–Pople Hamiltonian to model the interacting pi electrons in these systems. Within the DMRG method we have used an innovative scheme to target excited states with large transition dipole to the ground state. This method reproduces exact optical gaps and polarization in systems where exact diagonalization of the Hamiltonian is possible. We have used a correction vector method which tacitly takes into account the contribution of all excited states, to obtain the ground state polarizibility, first hyperpolarizibility, and two photon absorption cross sections. We find that the lowest optical excitations as well as the lowest excited triplet states are localized. It is interesting to note that the first hyperpolarizibility saturates more rapidly with system size compared to linear polarizibility unlike that of linear polyenes.

Structure and phase transition of the (1-x)PbTiO3-(x)BiAlO3 system

The system (1 - x)PbTiO3 - (x)BiAlO3 has been studied with regard to its structure and phase transition behaviour using x-ray, neutron and dielectric measurements. The structure is tetragonal within the solid solubility limit (x < 0.25). Interatomic distance analysis revealed that the Pb-O bond lengths remain unaffected and the only the Ti-O lengths are altered by BiAlO3 substitution. The results are suggestive of a crossover from and A-site & B site driven ferroelectric system for x = 0 to a dominant A-site driven ferroelectric system for higher x. This cross-over is brought about by (i) a reduction in the contribution to the ferroelectric stability from B-site cations due to dilution of the Ti-sublattice by Al and (ii) a reinforcement to the stability of the ferroelectric state by the A-site cations by the Bi+3 cations..

S-matrix for magnons in the D1-D5 system

We show that integrability and symmetries of the near horizon geometry of the D1-D5 system determine the S-matrix for the scattering of magnons with polarizations in AdS(3) x S-3 completely up to a phase. Using semi-classical methods we evaluate the phase to the leading and to the one-loop approximation in the strong coupling expansion. We then show that the phase obeys the unitarity constraint implied by the crossing relations to the one-loop order. We also verify that the dispersion relation obeyed by these magnons is one-loop exact at strong coupling which is consistent with their BPS nature.

Large negative magnetoresistance and metal-insulator transition observed in MnO/C composite coatings grown on ceramic alumina by metalorganic chemical vapor deposition

MnO/C composite coatings were grown by the metalorganic chemical vapor deposition process on ceramic alumina in argon ambient. Characterization by various techniques confirms that these coatings are homogeneous composites comprising nanometer-sized MnO particles embedded in a matrix of nanometer-sized graphite. Components of the MnO/C composite coating crystalline disordered, but are electrically quite conductive. Resistance vs. temperature measurements show that coating resistance increases exponentially from a few hundred ohms at room temperature to a few megaohms at 30 K. Logarithmic plots of reduced activation energy vs. temperature show that the coating material undergoes a metal-insulator transition. The reduced activation energy exponent for the film under zero magnetic field was 2.1, which is unusually high, implying that conduction is suppressed at much faster rate than the Mott or the Efros-Shklovskii hopping mechanism. Magnetoconductance us. magnetic field plots obtained at various temperatures show a high magnetoconductance (similar to 28.8%) at 100 K, which is unusually large for a disordered system, wherein magnetoresistance is attributed typically to weak localization. A plausible explanation for the unusual behavior observed in the carbonaceous disordered composite material is proposed. (C) 2010 Elsevier Ltd. All rights reserved.

Molecular Dynamics Investigation of Structure and Transport in the K2O-2SiO2 System Using a Partial Charge Based Model Potential

Potassium disilicate glass and melt have been investigated by using a new partial charge based potential model in which nonbridging oxygens are differentiated from bridging oxygens by their charges. The model reproduces the structural data pertaining to the coordination polyhedra around potassium and the various bond angle distributions excellently. The dynamics of the glass has been studied by using space and time correlation functions. It is found that K ions migrate by a diffusive mechanism in the melt and by hops below the glass transition temperature. They are also found to migrate largely through nonbridging oxygenrich sites in the silicate matrix, thus providing support to the predictions of the modified random network model.