Multi-time scale analysis applied to long time simulation of power systems

The results are presented of applying multi-time scale analysis using the singular perturbation technique for long time simulation of power system problems. A linear system represented in state-space form can be decoupled into slow and fast subsystems. These subsystems can be simulated with different time steps and then recombined to obtain the system response. Simulation results with a two-time scale analysis of a power system show a large saving in computational costs.

Weak discontinuities in relativisitic MHD

By using singular surface theory and ray theory the speeds of propagation of fast and slow waves, propagating into a medium in arbitrary motion, have been obtained in relativistic magnetohydrodynamics. The differential equation governing the growth of these waves along the rays has been derived and the solution has been presented in integral form.

Biaxial-load effects on isopachics for a crack at the interface of dissimilar media

The plane problem of two dissimilar materials, bonded together and containing a crack along their common interface, which were subjected to a biaxial load at infinity, is examined by giving a closed-form expression for the first stress invariant of the normal stresses, which is equally valid everywhere, near to, and far from, the crack-tip region. This exact expression for the first-stress invariant is compared by constructing the respective isopachic-fringe patterns, to the approximate expression with non-singular terms, due to the biaxiality factor, for the same quantity. Significant differences between respective isopachic-patterns were found and their dependence on the elastic properties of both materials and the applied loads was demonstrated. The relative errors between the computedK I - andK II -components by using the approximate expression for the first stress-invariant and the accurate one, derived from closed-form solution along either isopachic-fringes or along circles and radii from the crack-tip have been given, indicating in some cases large discrepancies between exact and approximate solutions.

On a finite element model for the analysis of through cracks in laminated anisotropic cylindrical shells

A finite element model for the analysis of laminated composite cylindrical shells with through cracks is presented. The analysis takes into account anisotropic elastic behaviour, bending-extensional coupling and transverse shear deformation effects. The proposed finite element model is based on the approach of dividing a cracked configuration into triangular shaped singular elements around the crack tip with adjoining quadrilateral shaped regular elements. The parabolic isoparametric cylindrical shell elements (both singular and regular) used in this model employ independent displacement and rotation interpolation in the shell middle surface. The numerical comparisons show the evidence to the conclusion that the proposed model will yield accurate stress intensity factors from a relatively coarse mesh. Through the analysis of a pressurised fibre composite cylindrical shell with an axial crack, the effect of material orthotropy on the crack tip stress intensity factors is shown to be quite significant.

Regularized numerical integration of multibody dynamics with the generalized alpha method*

This paper discusses the consistent regularization property of the generalized α method when applied as an integrator to an initial value high index and singular differential-algebraic equation model of a multibody system. The regularization comes from within the discretization itself and the discretization remains consistent over the range of values the regularization parameter may take. The regularization involves increase of the smallest singular values of the ill-conditioned Jacobian of the discretization and is different from Baumgarte and similar techniques which tend to be inconsistent for poor choice of regularization parameter. This regularization also helps where pre-conditioning the Jacobian by scaling is of limited effect, for example, when the scleronomic constraints contain multiple closed loops or singular configuration or when high index path constraints are present. The feed-forward control in Kane's equation models is additionally considered in the numerical examples to illustrate the effect of regularization. The discretization presented in this work is adopted to the first order DAE system (unlike the original method which is intended for second order systems) for its A-stability and same order of accuracy for positions and velocities.

Application of Artificial Viscosity in Establishing Supercritical Solutions to the Transonic Integra

The nonlinear singular integral equation of transonic flow is examined in the free-stream Mach number range where only solutions with shocks are known to exist. It is shown that, by the addition of an artificial viscosity term to the integral equation, even the direct iterative scheme, with the linear solution as the initial iterate, leads to convergence. Detailed tables indicating how the solution varies with changes in the parameters of the artificial viscosity term are also given. In the best cases (when the artificial viscosity is smallest), the solutions compare well with known results, their characteristic feature being the representation of the shock by steep gradients rather than by abrupt discontinuities. However, 'sharp-shock solutions' have also been obtained by the implementation of a quadratic iterative scheme with the 'artificial viscosity solution' as the initial iterate; the converged solution with a sharp shock is obtained with only a few more iterates. Finally, a review is given of various shock-capturing and shock-fitting schemes for the transonic flow equations in general, and for the transonic integral equation in particular, frequent comparisons being made with the approach of this paper.

Constrained linear minimax filter for systems with large plant uncertainties

This paper presents a method of designing a minimax filter in the presence of large plant uncertainties and constraints on the mean squared values of the estimates. The minimax filtering problem is reformulated in the framework of a deterministic optimal control problem and the method of solution employed, invokes the matrix Minimum Principle. The constrained linear filter and its relation to singular control problems has been illustrated. For the class of problems considered here it is shown that the filter can he constrained separately after carrying out the mini maximization. Numorieal examples are presented to illustrate the results.

Planar pursuit-evasion with variable speeds, part 1, extremal trajectory maps

Pursuit evasion in a plane is formulated with both players allowed to vary their speeds between fixed limits. A suitable choice of real-space coordinates confers open-loop optimality on the game. The solution in the small is described in terms of the individual players'' extremal trajectory maps (ETM). Each map is independent of role, adversary, and capture radius. An ETM depicts the actual real-space trajectories. A template method of generating constant control arcs is described. Examples of ETM for an aircraft flying at a constant altitude with fixed and varying speeds are presented.

Resistance Minimum in Dilute Magnetic Alloys in the Absence of Impurity Moments

A possible mechanism for the resistance minimum in dilute alloys in which the localized impurity states are non-magnetic is suggested. The fact is considered that what is essential to the Kondo-like behaviour is the interaction of the conduction electron spin s with the internal dynamical degrees of freedom of the impurity centre. The necessary internal dynamical degrees of freedom are provided by the dynamical Jahn-Teller effect associated with the degenerate 3d-orbitals of the transition-metal impurities interacting with the surrounding (octahedral) complex of the nearest-neighbour atoms. The fictitious spin I characterizing certain low-lying vibronic states of the system is shown to couple with the conduction electron spin s via s-d mixing and spin-orbit coupling, giving rise to a singular temperature-dependent exchange-like interaction. The resistivity so calculated is in fair agreement with the experimental results of Cape and Hake for Ti containing 0.2 at% of Fe.

The jump in vorticity across a shock wave in relativistic hydrodynamics

Using the singular surface theory, an expression for the jump in vorticity across a shock wave of arbitrary shape propagating in a uniform, perfect fluid occupying the space-time of special relativity, has been derived. It has been shown that the jump in vorticity across a shock of given strength and curvature depends only on the velocity of the medium ahead of the shock.

Symmetries and constraints in generalized Hamiltonian dynamics

A general analysis of symmetries and constraints for singular Lagrangian systems is given. It is shown that symmetry transformations can be expressed as canonical transformations in phase space, even for such systems. The relation of symmetries to generators, constraints, commutators, and Dirac brackets is clarified.

Discontinuities in an arbitrarily moving gas in special relativity

Singular surface theory and ray theory are used to study the propagation of a weak discontinuity in an arbitrarily moving gas within the framework of special relativity. A differential equation is obtained describing the variation of the strength of the discontinuity along rays.

Korteweg-De Vries Equation for Non-Linear Ion-Acoustic Waves in a Plasma with Negative Ions

The theoretical analysis, based on the perturbation technique, of ion-acoustic waves in the vicinity of a Korteweg-de Vries (K-dV) equation derived in a plasma with some negative ions has been made. The investigation shows that the negative ions in plasma with isothermal electrons introduced a critical concentration at which the ion-acoustic wave plays an important role of wave-breaking and forming a precursor while the plasma with non-isothermal electrons has no such singular behaviour of the wave. These two distinct features of ion waves lead to an overall different approach of present study of ion-waves. A distinct feature of non-uniform transition from the nonisothermal case to isothermal case has been shown. Few particular plasma models have been chosen to show the characteristics behaviour of the ion-waves existing in different cases

Semi-Classical Mechanics in Phase Space: A Study of Wigner's Function

We explore the semi-classical structure of the Wigner functions ($\Psi $(q, p)) representing bound energy eigenstates $|\psi \rangle $ for systems with f degrees of freedom. If the classical motion is integrable, the classical limit of $\Psi $ is a delta function on the f-dimensional torus to which classical trajectories corresponding to ($|\psi \rangle $) are confined in the 2f-dimensional phase space. In the semi-classical limit of ($\Psi $ ($\hslash $) small but not zero) the delta function softens to a peak of order ($\hslash ^{-\frac{2}{3}f}$) and the torus develops fringes of a characteristic 'Airy' form. Away from the torus, $\Psi $ can have semi-classical singularities that are not delta functions; these are discussed (in full detail when f = 1) using Thom's theory of catastrophes. Brief consideration is given to problems raised when ($\Psi $) is calculated in a representation based on operators derived from angle coordinates and their conjugate momenta. When the classical motion is non-integrable, the phase space is not filled with tori and existing semi-classical methods fail. We conjecture that (a) For a given value of non-integrability parameter ($\epsilon $), the system passes through three semi-classical regimes as ($\hslash $) diminishes. (b) For states ($|\psi \rangle $) associated with regions in phase space filled with irregular trajectories, ($\Psi $) will be a random function confined near that region of the 'energy shell' explored by these trajectories (this region has more than f dimensions). (c) For ($\epsilon \neq $0, $\hslash $) blurs the infinitely fine classical path structure, in contrast to the integrable case ($\epsilon $ = 0, where $\hslash $ )imposes oscillatory quantum detail on a smooth classical path structure.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.