Comparison of BR and QR Eigenvalue Algorithms for Power System Small Signal Stability Analysis

The BR algorithm is a novel and efficient method to find all eigenvalues of upper Hessenberg matrices and has never been applied to eigenanalysis for power system small signal stability. This paper analyzes differences between the BR and the QR algorithms with performance comparison in terms of CPU time based on stopping criteria and storage requirement. The BR algorithm utilizes accelerating strategies to improve its performance when computing eigenvalues of narrowly banded, nearly tridiagonal upper Hessenberg matrices. These strategies significantly reduce the computation time at a reasonable level of precision. Compared with the QR algorithm, the BR algorithm requires fewer iteration steps and less storage space without depriving of appropriate precision in solving eigenvalue problems of large-scale power systems. Numerical examples demonstrate the efficiency of the BR algorithm in pursuing eigenanalysis tasks of 39-, 68-, 115-, 300-, and 600-bus systems. Experiment results suggest that the BR algorithm is a more efficient algorithm for large-scale power system small signal stability eigenanalysis.

A study of the stability of subcycling algorithms in structural dynamics

Algorithms for explicit integration of structural dynamics problems with multiple time steps (subcycling) are investigated. Only one such algorithm, due to Smolinski and Sleith has proved to be stable in a classical sense. A simplified version of this algorithm that retains its stability is presented. However, as with the original version, it can be shown to sacrifice accuracy to achieve stability. Another algorithm in use is shown to be only statistically stable, in that a probability of stability can be assigned if appropriate time step limits are observed. This probability improves rapidly with the number of degrees of freedom in a finite element model. The stability problems are shown to be a property of the central difference method itself, which is modified to give the subcycling algorithm. A related problem is shown to arise when a constraint equation in time is introduced into a time-continuous space-time finite element model. (C) 1998 Elsevier Science S.A.

Structure determination of the three disulfide bond isomers of alpha-conotoxin GI: A model for the role of disulfide bonds in structural stability

The three possible disulfide bonded isomers of alpha-conotoxin GI have been selectively synthesised and their structures determined by H-1 NMR spectroscopy. alpha-Conotoxin GI derives from the venom of Conus geographus and is a useful neuropharmacological tool as it selectively binds to the nicotinic acetylcholine receptor (nAChR), a ligand-gated ion channel involved in nerve signal transmission. The peptide has the sequence ECCNPACGRHYSC-NH2, and the three disulfide bonded isomers are referred to as GI(2-7;3-13), GI(2-13;3-7) and GI(2-3;7-13). The NMR structure for the native isomer GI(2-7;3-13) is of excellent quality, with a backbone pairwise RMSD of 0.16 Angstrom for a family of 35 structures, and comprises primarily a distorted 3(10),, helix between residues 5 to 11. The two non-native isomers exhibit multiple conformers in solution, with the major populated forms being different in structure both from each other and from the native form. Structure-activity relationships for the native GI(2-7;3-13) as well as the role of the disulfide bonds on folding and stability of the three isomers are examined. It is concluded that the disulfide bonds in alpha-conotoxin GI play a crucial part in determining both the structure and stability of the peptide. A trend for increased conformational heterogeneity was observed in the order of GI(2-7;3-13) < GI(2-13;3-7) < GI(2-3;7-13). It was found that the peptide bond joining Cys2 to Cys3 in GI(2-3;7-13) is predominantly trans, rather than cis as theoretically predicted. These structural data are used to interpret the varying nAChR binding of the non-native forms. A model for the binding of native GI(2-7;3-13) to the mammalian nAChR is proposed, with an alpha-subunit binding face made up of Cys2, Asn4, Pro5, Ala6 and Cys7 and a selectivity face, comprised of Arg9 and His10. These two faces orient the molecule between the alpha and delta subunits of the receptor. The structure of the CCNPAC sequence of the native GI(2-7;3-13) is compared to the structure of the identical sequence from the toxic domain of heat-stable enterotoxins, which forms part of the receptor binding region of the enterotoxins, but which has a different disulfide connectivity. (C) 1998 Academic Press Limited.

Thickness and stability of adsorbed film in cylindrical mesopores

The adsorbed film in small cylindrical mesopores is studied by using MCM-41 samples of uniform cylindrical channels as model systems. It is found that at a given relative pressure, the smaller the pore radius, the thicker the adsorbed film is, as postulated by Broekhoff and De Beer. Thermodynamics analysis established that the stability of the adsorbed film is determined by interface curvature and the potential of interaction between adsorbate and adsorbent. A semiempirical equation is proposed to describe the state of stable adsorbed films in cylindrical mesopores. It is also shown to be useful in calculations of pore size distributions of mesoporous solids.

Analysis of small signal stability margins using genetic optimization

Power system small signal stability analysis aims to explore different small signal stability conditions and controls, namely: (1) exploring the power system security domains and boundaries in the space of power system parameters of interest, including load flow feasibility, saddle node and Hopf bifurcation ones; (2) finding the maximum and minimum damping conditions; and (3) determining control actions to provide and increase small signal stability. These problems are presented in this paper as different modifications of a general optimization to a minimum/maximum, depending on the initial guesses of variables and numerical methods used. In the considered problems, all the extreme points are of interest. Additionally, there are difficulties with finding the derivatives of the objective functions with respect to parameters. Numerical computations of derivatives in traditional optimization procedures are time consuming. In this paper, we propose a new black-box genetic optimization technique for comprehensive small signal stability analysis, which can effectively cope with highly nonlinear objective functions with multiple minima and maxima, and derivatives that can not be expressed analytically. The optimization result can then be used to provide such important information such as system optimal control decision making, assessment of the maximum network's transmission capacity, etc. (C) 1998 Elsevier Science S.A. All rights reserved.

The Modified Card Sorting Test: Test-retest stability and relationships with demographic variables in a healthy older adult sample

Objectives. To investigate the test-retest stability of a standardized version of Nelson's (1976) Modified Card Sorting Test (MCST) and its relationships with demographic variables in a sample of healthy older adults. Design. A standard card order and administration were devised for the MCST and administered to participants at an initial assessment, and again at a second session conducted a minimum of six months later in order to examine its test-retest stability. Participants were also administered the WAIS-R at initial assessment in order to provide a measure of psychometric intelligence. Methods. Thirty-six (24 female, 12 male) healthy older adults aged 52 to 77 years with mean education 12.42 years (SD = 3.53) completed the MCST on two occasions approximately 7.5 months (SD = 1.61) apart. Stability coefficients and test-retest differences were calculated for the range of scores. The effect of gender on MCST performance was examined. Correlations between MCST scores and age, education and WAIS-R IQs were also determined. Results. Stability coefficients ranged from .26 for the percent perseverative errors measure to .49 for the failure to maintain set measure. Several measures were significantly correlated with age, education and WAIS-R IQs, although no effect of gender on MCST performance was found. Conclusions. None of the stability coefficients reached the level required for clinical decision making. The results indicate that participants' age, education, and intelligence need to be considered when interpreting MCST performance. Normative studies of MCST performance as well as further studies with patients with executive dysfunction are needed.

Preparatory trunk motion accompanies rapid upper limb movement

Evaluation of trunk movements, trunk muscle activation, intra-abdominal pressure and displacement of centres of pressure and mass was undertaken to determine whether trunk orientation is a controlled variable prior to and during rapid bilateral movement of the upper limbs. Standing subjects performed rapid bilateral symmetrical upper limb movements in three directions (flexion, abduction and extension). The results indicated a small (0.4-3.3 degrees) but consistent initial angular displacement between the segments of the trunk in a direction opposite to that produced by the reactive moments resulting from limb movement. Phasic activation of superficial trunk muscles was consistent with this pattern of preparatory motion and with the direction of motion of the centre of mass. In contrast, activation of the deep abdominal muscles was independent of the direction of limb motion, suggesting a non-direction specific contribution to spinal stability. The results support the opinion that feedforward postural responses result in trunk movements, and that orientation of the trunk and centre of mass are both controlled variables in relation to rapid limb movements.

Complexes of 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) - Synthesis, structure and stability constants

The 12-membered macrocyclic ligand 1-thia-4,7, 10-triazacyclododecane ([12]aneN(3)S) has been synthesised, although upon crystallization from acetonitrile a product in which carbon dioxide had added to one secondary amine in the macrocyclic ring (H[12]aneN(3)SCO(2). H2O) was isolated and subsequently characterised by X-ray crystallography. The protonation constants for [12]aneN(3)S and stability constants with Zn(II), Pb(II), Cd(II) and Cu(II) have been determined either potentiometrically or spectrophotometrically in aqueous solution, and compared with those measured or reported for the ligands 1-oxa-4,7,10-triazacyclododecane ([12]aneN(3)O) and 1,4,7,10-tetraazacyclododecane ([12]aneN(4)). The magnitudes of the stability constants are consistent with trends observed previously for macrocyclic ligands as secondary amine donors are replaced with oxygen and thioether donors although the stability constant for the [Hg([12]aneN(4))](2+) complex has been estimated from an NMR experiment to be at least three orders of magnitude larger than reported previously. Zinc(II), mercury(II), lead(II), copper(II) and nickel(II) complexes of [12]aneN(3)S have been isolated and characterised by X-ray crystallography. In the case of copper(II), two complexes [Cu([12]aneN(3)S)(H2O)](ClO4)(2) and [Cu-2([12]aneN(3)S)(2)(OH)(2)](ClO4)(2) were isolated, depending on the conditions employed. Molecular mechanics calculations have been employed to investigate the relative metal ion size preferences of the [3333], asym-[2424] and sym-[2424] conformation isomers. The calculations predict that the asym-[2424] conformer is most stable for M-N bond lengths in the range 2.00-2.25 Angstrom whilst for the larger metal ions the [3333] conformer is dominant. The disorder seen in the structure of the [Zn([12]aneN(3)S)(NO3)](+) complex is also explained by the calculations. (C) 1999 Elsevier Science Ltd. All rights reserved.

Activation of the human diaphragm during a repetitive postural task

1. The co-ordination between respiratory and postural functions of the diaphragm was investigated during repetitive upper Limb movement. It was hypothesised that diaphragm activity would occur either tonically or phasically in association with the forces from each movement and that this activity would combine with phasic respiratory activity. 2. Movements of the upper limb and ribcage were measured while standing subjects performed repetitive upper limb movements 'as fast as possible'. Electromyographic (EMG) recordings of the costal diaphragm were made using intramuscular electrodes in four subjects. Surface electrodes were placed over the deltoid and erector spinae muscles. 3. In contrast to standing at rest, diaphragm activity was present throughout expiration at 78 +/- 17% (mean +/- S.D.) of its peak inspiratory magnitude during repeated upper limb movement. 4. Bursts of deltoid and erector spinae EMG activity occurred at the Limb movement frequency (similar to 2.9 Hz). Although the majority of diaphragm EMG power was at the respiratory frequency (similar to 0.4 Hz), a peak was also present at the movement frequency. This finding was corroborated by averaged EMG activity triggered from upper limb movement. In addition, diaphragm EMG activity was coherent with ribcage motion at the respiratory frequency and with upper limb movement at the movement frequency. 5. The diaphragm response was similar when movement was performed while sitting. In addition, when subjects moved with increasing frequency the peak upper limb acceleration correlated with diaphragm EMG amplitude. These findings support the argument that diaphragm contraction is related to trunk control. 6. The results indicate that activity of human phrenic motoneurones is organised such that it contributes to both posture and respiration during a task which repetitively challenges trunk posture.

In vivo measurement of the effect of intra-abdominal pressure on the human spine

In humans, intra-abdominal pressure (IAP) is elevated during many everyday activities. This experiment aimed to investigate the extent to which increased IAP-without concurrent activity of the abdominal or back extensor muscles-produces an extensor torque. With subjects positioned in side lying on a swivel table with its axis at L3, moments about this vertebral level were measured when IAP was transiently increased by electrical stimulation of the diaphragm via the phrenic nerve. There was no electromyographic activity in abdominal and back extensor muscles. When IAP was increased artificially to similar to 15% of the maximum IAP amplitude that could be generated voluntarily with the trunk positioned in flexion, a trunk extensor moment (similar to6 Nm) was recorded. The size of the effect was proportional to the increase in pressure. The extensor moment was consistent with that predicted from a model based on measurements of abdominal cross-sectional area and IAP moment arm. When IAP was momentarily increased while the trunk was flexed passively at a constant velocity, the external torque required to maintain the velocity was increased. These results provide the first in vivo data of the amplitude of extensor moment that is produced by increased IAP. Although the net effect of this extensor torque in functional tasks would be dependent on the muscles used to increase the IAP and their associated flexion torque, the data do provide evidence that IAP contributes, at least in part, to spinal stability. (C) 2001 Elsevier Science Ltd. All rights reserved.

Determination of the stability constants of mercury(II) complexes of mixed donor macrobicyclic encapsulating ligands

The reactions of mercury(II) with the mixed donor encapsulating ligands 3,6,16-trithia-6,11,19-triazabicyclo[6.6.6]icosane (AMN(3)S(3)sar) and 1-amino-8-methyl-6,19-dithia-3,10,13,16-tetraazabicyclo[6.6.6]icosane (AMN(4)S(2)sar) have been studied. NMR ligand-ligand competition experiments with the ligands 1,4,8,11-tetraazaeyclotetradecane ([14]aneN(4)), 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) and ethylenediaminetetraacetic acid (EDTA) with AMN(3)S(3)sar and Hg(II) indicated that [14]aneN(4) would be an appropriate competing ligand for the, determination of the Hg(II) stability constant. Calculations indicated the ratio of concentrations of AMN3S3sar, [14]aneN(4) and Hg(II) required for the determination of the stability constant ranged from 1:1:1 to 1:5:1. Refinement of the titration curves yielded log(10)K[Hg(AMN(3)S(3)sar)](2+) = 17.7. A similar competition titration resulted in the determination of the stability constant for the AMN(4)S(2)sar system as log(10)K[Hg(AMN(4)S(2)sar)](2+) = 19.5. The observed binding constants for the mixed N/S donor systems and the hexaaza analogues sar (3,6,10,13,16,19-hexaazabicyclo [6.6.6]icosane) and diamsar (1,8-diamino-3,6,10,13,16,19 -hexazabicyclo [6.6.6] icosane (log(10)K-[Hg(diamsar)](2+) = 26.4; log(10)K[Hg(sar)](2+) = 28.1) differ by approximately ten orders of magnitude. The difference is ascribed not to a cryptate effect but to a mismatch in the Hg-N and Hg-S bond lengths in the N/S systems.

Thermal stability of mixed-cation alpha-sialon ceramics

A series of alpha-sialon (alpha') compositions containing mixed stabilising cations were prepared, by introducing additional CaO to a basic Sm alpha-sialon compositions. The thermal stability of these Sm-Ca-containing alpha-sialon phases was investigated using XRD, SEM and EDXS techniques. It was found that the addition of calcium into the Sm alpha-sialon systems greatly improved the stability of the alpha-sialon phases. Calcium was found to be incorporated into the alpha-sialon structure, coexistent with the samarium, and partitioning of the calcium and samarium was observed between the alpha' phase and grain boundary phases. This indicates a technique which may be used to improve the thermal stability of the alpha' phase while maintaining good refractory phases at the gialon grain boundaries. (C) 2003 Elsevier Science B.V. All rights reserved.

α′ phase stability in Nd-Li-sialon systems

The formability and stability of the alpha-sialon (alpha') phase was investigated in multi-cation Nd-Li-sialon systems. Four samples were prepared, ranging from a pure Nd-sialon to a pure Li-sialon, with two intermediate samples being prepared with either lithium or neodymium replacing the other alpha'-stabilising additive by 20 eq.%, as to maintain an equivalent design composition in all samples. After sintering, all samples were subsequently heat treated up to 192 h at 1450 and 1300 degreesC. While significant quantities of the beta'-sialon (beta' phase were found in most samples, the high-lithium Li-Nd-sialon sample was found to be almost pure a' phase after sintering. Furthermore, the long-term stability of the a' phase on heat treatment was also found to be superior in both multi-cation samples than in either of the single-alpha'-stabilising-cation samples. This is thought to be related to improved retention of the lithium in the multi-cation systems, as much of the lithium was found to volatilise during sintering in the neodymium-free sample. (C) 2002 Elsevier Science Ltd. All rights reserved.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.