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.

Pillaring of smectites using an aluminium oligomer: A study of pillar density and thermal stability

Two smectite samples having different layer charges were pillared using hydroxy aluminium oligomers at a OH/Al ratio of 2.5 and at pH 4.3 to 4.6. Pillaring was carried out at different conditions such as ageing, temperature and base addition time of the pillaring solution, and also in the presence of nonionic surfactant polyoxyethylene sorbitanmonooleate (Tween-80). The primary objective of preparing at different conditions was to introduce varied quantities of aluminium oligomer between the layers and to study its effect on the properties of the pillared products. A simple method has been followed to estimate the amount of interlayer aluminium. A quantity called pillar density number (PDN) based on the ratio of interlayer Al adsorbed to CEC of the parent clay has been effectively used to evaluate the nature of the resulting pillared product. PDN, for a given clay, was found to correlate well with the sharpness of the d(001) peaks for the air dried samples. The calculated number of pillars, varied from 3.00 x 10(18) to 5.32 x 10(18) per meq charge. The present study shows that a higher value of PDN is indicative of better thermal stability. Pillar density number may be conveniently used as a measure of the thermal stability of pillared samples.

Response, loads and stability of rotors with interconnected blades

An isolated rotor with blades interconnected through viscoelastic elements is analyzed for response, loads and stability in moment trim under forward flight conditions. A conceptual model of a multibladed rotor with rigid flap and lag motions is considered, Although the interconnecting elements are placed in the In-plane direction, considerable coupling between the flap-lag motions of the blades can occur in certain ranges of interblade element stiffness. Interblade coupling can yield significant changes in the response, loads and stability which are dependent on the interblade element and rotor parameters.

Role of zirconium addition on stability of UAl3 phase in Al-U system

The microstructural changes of Al-22 wt%U and Al-46 wt%U alloys containing 3 wt% Zr were investigated after heat treatment at 620 degrees C for 1 to 45 days, Though it is reported that addition of similar to 3 wt% Zr stabilizes the (U,Zr)Al-3 phase at room temperature, the present investigation shows that the (U,Zr)Al-3 phase is not stable but slowly transforms to the U0.9Al4 phase, The high temperature creep curves generated for these ternary alloys showed a wavy pattern which also suggests that the (U,Zr)Al-3 phase is not stable.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Stability analysis of stochastically parametered nonconservative columns

Nonconservatively loaded columns. which have stochastically distributed material property values and stochastic loadings in space are considered. Young's modulus and mass density are treated to constitute random fields. The support stiffness coefficient and tip follower load are considered to be random variables. The fluctuations of external and distributed loadings are considered to constitute a random field. The variational formulation is adopted to get the differential equation and boundary conditions. The non self-adjoint operators are used at the boundary of the regularity domain. The statistics of vibration frequencies and modes are obtained using the standard perturbation method, by treating the fluctuations to be stochastic perturbations. Linear dependence of vibration and stability parameters over property value fluctuations and loading fluctuations are assumed. Bounds for the statistics of vibration frequencies are obtained. The critical load is first evaluated for the averaged problem and the corresponding eigenvalue statistics are sought. Then, the frequency equation is employed to transform the eigenvalue statistics to critical load statistics. Specialization of the general procedure to Beck, Leipholz and Pfluger columns is carried out. For Pfluger column, nonlinear transformations are avoided by directly expressing the critical load statistics in terms of input variable statistics.

Stability of laminated composite plates with cut-outs

Critical buckling loads of laminated fibre-reinforced plastic square panels have been obtained using the finite element method. Various boundary conditions, lay-up details, fibre orientations, cut-out sizes are considered. A 36 degrees of freedom triangular element, based on the classical lamination theory (CLT) has been used for the analysis. The performance of this element is validated by comparing results with some of those available in literature. New results have been given for several cases of boundary conditions for [0°/ ± 45°/90°]s laminates. The effect of fibre-orientation in the ply on the buckling loads has been investigated by considering [±?]6s laminates.

Stability of stochastic Leipholz column with stochastic loading

The Leipholz column which is having the Young modulus and mass per unit length as stochastic processes and also the distributed tangential follower load behaving stochastically is considered. The non self-adjoint differential equation and boundary conditions are considered to have random field coefficients. The standard perturbation method is employed. The non self-adjoint operators are used within the regularity domain. Full covariance structure of the free vibration eigenvalues and critical loads is derived in terms of second order properties of input random fields characterizing the system parameter fluctuations. The mean value of critical load is calculated using the averaged problem and the corresponding eigenvalue statistics are sought. Through the frequency equation a transformation is done to yield load parameter statistics. A numerical study incorporating commonly observed correlation models is reported which illustrates the full potentials of the derived expressions.

Engineered dimer interface mutants of triosephosphate isomerase: the role of inter-subunit interactions in enzyme function and stability

The role of inter-subunit interactions in maintaining optimal catalytic activity in triosephosphate isomerase (TIM) has been probed, using the Plasmodium falciparum enzyme as a model. Examination of subunit interface contacts in the crystal structures suggests that residue 75 (Thr, conserved) and residue 13 (Cys, variable) make the largest number of inter-subunit contacts. The mutants Cys13Asp (C13D) and Cys13Glu (C13E) have been constructed and display significant reduction in catalytic activity when compared with wild-type (WT) enzyme (similar to 7.4-fold decrease in k(cat) for the C13D and similar to 3.3-fold for the C13E mutants). Analytical gel filtration demonstrates that the C13D mutant dissociates at concentrations < 1.25 mu M, whereas the WT and the C13E enzymes retain the dimeric structure. The order of stability of the mutants in the presence of chemical denaturants, like urea and guanidium chloride, is WT > Cys13Glu > Cys13Asp. Irreversible thermal precipitation temperatures follow the same order as well. Modeling studies establish that the Cys13Asp mutation is likely to cause a significantly greater structural perturbation than Cys13Glu. Analysis of sequence and structural data for TIMs from diverse sources suggests that residues 13 and 82 form a pair of proximal sites, in which a limited number of residue pairs may be accommodated.

Microstructural stability and nanoindentation of an electrodeposited nanocrystalline Ni-1.5wt.%P alloy

Microstructural stability of nanocrystalline Ni-1.5wt.%P alloy with an initial grain size of 3 nm processed by pulsed electrodeposition was studied using differential scanning calorimetry (DSC) and annealing. Microstructural characterization suggests that the observed exothermic peak during heating in DSC is related to both concurrent grain growth and Ni3P formation. Nanoindentation on samples with grain sizes from 3 to 50 nm revealed a breakdown in Hall-Petch strengthening in nano Ni-P alloy at grain sizes <= 10 nm, consistent with some previous observations. It is concluded that there is a grain boundary weakening regime for grain sizes < 10 nm, based on analysis which show that the data cannot be rationalized in terms of microstrain relaxation, variation in elastic modulus, texture evolution and duplex structure formation.

Thermal Stability of Spherical Nanoporous Aggregates and Formation of Hollow Structures by Sintering-A Phase-Field Study

Nanoporous structures are widely used for many applications and hence it Is important to investigate their thermal stability. We study the stability of spherical nanoporous aggregates using phase-field simulations that explore systematically the effect of grain boundary diffusion, surface diffusion, and grain boundary mobility on the pathways for microstructural evolution. Our simulations for different combinations of surface and GB diffusivity and GB mobility show four distinct microstructural pathways en route to 100% density: multiple dosed pores, hollow shells, hollow shells with a core, and multiple interconnected pores. The microstructures from our simulations are consistent with experimental observations in several different systems. Our results have important implications for rational synthesis of hollow nanostructures or aggregates with open pores, and for controlling the stability of nanoporous aggregates that are widely used for many applications.

Optical and electrical properties of hydrogen-passivated gallium antimonide

The effect of hydrogen-plasma passivation on the optical and electrical properties of gallium antimonide bulk single crystals is presented. Fundamental changes of the radiative recombination after hydrogenation in undoped, zinc-doped, tellurium-doped, and codoped (with Zn and Te) GaSb are reported. The results of optical measurements indicate that passivation of acceptors is more efficient than that of the donors and, in general, the passivation efficiency depends on the doping level. Passivation of deep nonradiative centers is reflected by the gain of photoluminescence intensity and decrease in deep-level transient spectroscopy peak height. Extended defects like grain boundaries and dislocations have also been found to be passivated. The thermal stability of the passivated deep level and extended defects is higher than that of the shallow level. The kinetics of thermally released hydrogen in the bulk has been studied by reverse-bias annealing experiments.

Stability of spatially developing boundary layers in pressure gradients

A new formulation of the stability of boundary-layer flows in pressure gradients is presented, taking into account the spatial development of the flow and utilizing a special coordinate transformation. The formulation assumes that disturbance wavelength and eigenfunction vary downstream no more rapidly than the boundary-layer thickness, and includes all terms nominally of order R(-1) in the boundary-layer Reynolds number R. In Blasius flow, the present approach is consistent with that of Bertolotti et al. (1992) to O(R(-1)) but simpler (i.e. has fewer terms), and may best be seen as providing a parametric differential equation which can be solved without having to march in space. The computed neutral boundaries depend strongly on distance from the surface, but the one corresponding to the inner maximum of the streamwise velocity perturbation happens to be close to the parallel flow (Orr-Sommerfeld) boundary. For this quantity, solutions for the Falkner-Skan flows show the effects of spatial growth to be striking only in the presence of strong adverse pressure gradients. As a rational analysis to O(R(-1)) demands inclusion of higher-order corrections on the mean flow, an illustrative calculation of one such correction, due to the displacement effect of the boundary layer, is made, and shown to have a significant destabilizing influence on the stability boundary in strong adverse pressure gradients. The effect of non-parallelism on the growth of relatively high frequencies can be significant at low Reynolds numbers, but is marginal in other cases. As an extension of the present approach, a method of dealing with non-similar flows is also presented and illustrated. However, inherent in the transformation underlying the present approach is a lower-order non-parallel theory, which is obtained by dropping all terms of nominal order R(-1) except those required for obtaining the lowest-order solution in the critical and wall layers. It is shown that a reduced Orr-Sommerfeld equation (in transformed coordinates) already contains the major effects of non-parallelism.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.