Phase-plane techniques for the solution of transient-stability problems

The paper presents a graphical-numerical method for determining the transient stability limits of a two-machine system under the usual assumptions of constant input, no damping and constant voltage behind transient reactance. The method presented is based on the phase-plane criterion,1, 2 in contrast to the usual step-by-step and equal-area methods. For the transient stability limit of a two-machine system, under the assumptions stated, the sum of the kinetic energy and the potential energy, at the instant of fault clearing, should just be equal to the maximum value of the potential energy which the machines can accommodate with the fault cleared. The assumption of constant voltage behind transient reactance is then discarded in favour of the more accurate assumption of constant field flux linkages. Finally, the method is extended to include the effect of field decrement and damping. A number of examples corresponding to each case are worked out, and the results obtained by the proposed method are compared with those obtained by the usual methods.

Techniques for Analysis of Certain Nonlinear Systems

This paper suggests the use of simple transformations like ÿ=kx, kx2 for second-order nonlinear differential equations to effect rapid plotting of the phase-plane trajectories. The method is particularly helpful in determining quickly the trajectory slopes along simple curves in any desired region of the phase plane. New planes such as the tÿ-x, tÿ2-x are considered for the study of some groups of nonlinear time-varying systems. Suggestions for solving certain higher-order nonlinear systems are also made.

The x¿n-x Plane for Analysis of Certain Second-Order Nonlinear Systems

Analysis of certain second-order nonlinear systems, not easily amenable to the phase-plane methods, and described by either of the following differential equations xÿn-2ÿ+ f(x)xÿ2n+g(x)xÿn+h(x)=0 ÿ+f(x)xÿn+h(x)=0 n≫0 can be effected easily by drawing the entire portrait of trajectories on a new plane; that is, on one of the xÿnÿx planes. Simple equations are given to evaluate time from a trajectory on any of these n planes. Poincaré's fundamental phase plane xÿÿx is conceived of as the simplest case of the general xÿnÿx plane.

Rocking response of rectangular rigid blocks under random noise base excitations

The response of a rigid rectangular block resting on a rigid foundation and acted upon simultaneously by a horizontal and a vertical random white-noise excitation is considered. In the equation of motion, the energy dissipation is modeled through a viscous damping term. Under the assumption that the body does not topple, the steady-state joint probability density function of the rotation and the rotational velocity is obtained using the Fokker-Planck equation approach. Closed form solution is obtained for a specific combination of system parameters. A more general but approximate solution to the joint probability density function based on the method of equivalent non-linearization is also presented. Further, the problem of overturning of the block is approached in the framework of the diffusion methods for first passage failure studies. The overturning of the block is deemed incipient when the response trajectories in the phase plane cross the separatrix of the conservative unforced system. Expressions for the moments of first passage time are obtained via a series solution to the governing generalized Pontriagin-Vitt equations. Numerical results illustra- tive of the theoretical solutions are presented and their validity is examined through limited amount of digital simulations.

The R-x, R2-x Planes; Some New Ideas for the Study of Nonlinear Systems

Use of some new planes such as the R-x, R2-x (where R represents in the n-dimensional phase space, the radius vector from the origin to any point on the trajectory described by the system) is suggested for analysis of nonlinear systems of any kind. The stability conditions in these planes are given. For easy understanding of the method, the transformation from the phase plane to the R-x, R2-x planes is brought out for second-order systems. In general, while these planes serve as useful as the phase plane, they have proved to be simpler in determining quickly the general behavior of certain classes of second-order nonlinear systems. A chart and a simple formula are suggested to evaluate time easily from the R-x and R2-x trajectories, respectively. A means of solving higher-order nonlinear systems is also illustrated. Finally, a comparative study of the trajectories near singular points on the phase plane and on the new planes is made.

A class of one-dimensional steady-state flows in radiation-gas-dynamics

The study of steady-state flows in radiation-gas-dynamics, when radiation pressure is negligible in comparison with gas pressure, can be reduced to the study of a single first-order ordinary differential equation in particle velocity and radiation pressure. The class of steady flows, determined by the fact that the velocities in two uniform states are real, i.e. the Rankine-Hugoniot points are real, has been discussed in detail in a previous paper by one of us, when the Mach number M of the flow in one of the uniform states (at x=+∞) is greater than one and the flow direction is in the negative direction of the x-axis. In this paper we have discussed the case when M is less than or equal to one and the flow direction is still in the negative direction of the x-axis. We have drawn the various phase planes and the integral curves in each phase plane give various steady flows. We have also discussed the appearance of discontinuities in these flows.

Exact nonlinear travelling hydromagnetic wave solutions

Exact travelling wave solutions for hydromagnetic waves in an exponentially stratified incompressible medium are obtained. With the help of two integrals it becomes possible to reduce the system of seven nonlinear PDE's to a second order nonlinear ODE which describes an one dimensional harmonic oscillator with a nonlinear friction term. This equation is studied in detail in the phase plane. The travelling waves are periodic only when they propagate either horizontally or vertically. The reduced second order nonlinear differential equation describing the travelling waves in inhomogeneous conducting media has rather ubiquitous nature in that it also appears in other geophysical systems such as internal waves, Rossby waves and topographic Rossby waves in the ocean.

Multiphase models for simulating hot and slightly miscible DNAPL (dense non-aqueous phase fluids) in a saturated rock fracture under deformation

In the present study a two dimensional model is first developed to show the behaviour of dense non-aqueous phase liquids (DNAPL) within a rough fracture. To consider the rough fracture, the fracture is imposed with variable apertures along its plane. It is found that DNAPL follows preferential pathways. In next part of the study the above model is further extended for non-isothermal DNAPL flow and DNAPL-water interphase mass transfer phenomenon. These two models are then coupled with joint deformation due to normal stresses. The primary focus of these models is specifically to elucidate the influence of joint alteration due to external stress and fluid pressures on flow driven energy transport and interphase mass transfer. For this, it is assumed that the critical value for joint alteration is associated with external stress and average of water and DNAPL pressures in multiphase system and the temporal and spatial evolution of joint alteration are determined for its further influence on energy transport and miscible phase transfer. The developed model has been studied to show the influence of deformation on DNAPL flow. Further this preliminary study demonstrates the influence of joint deformation on heat transport and phase miscibility via multiphase flow velocities. It is seen that the temperature profile changes and shows higher diffusivity due to deformation and although the interphase miscibility value decreases but the lateral dispersion increases to a considerably higher extent.

Stress-induced martensitic phase transformation in Cu-Zr nanowires

A novel stress-induced martensitic phase transformation in an initial < 100 >/{100} B2-CuZr nanowire is reported for the first time in this letter. Such behavior is observed in a nanowire with cross-sectional dimensions of 19.44 x 19.44 angstrom(2) over a temperature range of 100-400 K and at a strain rate of 1 x 10(9) s(-1) using atomistic simulations. Phase transformation from an initial B2 phase to a BCT (Body-Centered-Tetragonal) phase is observed via nucleation and propagation of {100} twinning plane under high strain rate tensile deformation. (C) 2009 Elsevier B.V. All rights reserved.

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.

Optimal Nonlinear Control and Estimation for a Reusable Launch Vehicle During Reentry Phase

In this paper a nonlinear optimal controller has been designed for aerodynamic control during the reentry phase of the Reusable Launch Vehicle (RLV). The controller has been designed based on a recently developed technique Optimal Dynamic Inversion (ODI). For full state feedback the controller has required full information about the system states. In this work an Extended Kalman filter (EKF) is developed to estimate the states. The vehicle (RLV) has been has been consider as a nonlinear Six-Degree-Of-Freedom (6-DOF) model. The simulation results shows that EKF gives a very good estimation of the states and it is working well with ODI. The resultant trajectories are very similar to those obtained by perfect state feedback using ODI only.

Single and multi-step phase transformation in CuZr nanowire under compressive/tensile loading

A novel stress induced martenistic phase transformation is reported in an initial B2-CuZr nanowire of cross-sectional dimensions in the range of 19.44 x 19.44-38.88 x 38.88 angstrom(2) and temperature in the range of 10-400 K under both tensile and compressive loading. Extensive Molecular Dynamic simulations are performed using an inter-atomic potential of type Finnis and Sinclair. The nanowire shows a phase transformation from an initial B2 phase to BCT (body-centered-tetragonal) phase with failure strain of similar to 40% in tension, whereas in compression, comparatively a small B2 -> BCT phase transformation is observed with failure strain of similar to 25%. Size and temperature dependent deformation mechanisms which control ultimately the B2 -> BCT phase transformation are found to be completely different for tensile and compressive loadings. Under tensile loading, small cross-sectional nanowire shows a single step phase transformation, i.e. B2 -> BCT via twinning along {100} plane, whereas nanowires with larger cross-sectional area show a two step phase transformation, i.e. B2 -> R phase -> BCT along with intermediate hardening. In the first step, nanowire shows phase transformation from B2 -> R phase via twinning along {100} plane, afterwards the nanowire deforms via twinning along {110} plane which cause further transformation from R phase -> BCT phase. Under compressive loading, the nanowire shows crushing along {100} plane after a single step phase transformation from B2 -> BCT. Proper tailoring of such size and temperature dependent phase transformation can be useful in designing nanowire for high strength applications with corrosion and fatigue resistance. (C) 2009 Elsevier Ltd. All rights reserved.

Phase diagram of a hard-sphere system in a quenched random potential: A numerical study

We report numerical results for the phase diagram in the density-disorder plane of a hard-sphere system in the presence of quenched, random, pinning disorder. Local minima of a discretized version of the Ramakrishnan-Yussouff free energy functional are located numerically and their relative stability is studied as a function of the density and the strength of disorder. Regions in the phase diagram corresponding to liquid, glassy, and nearly crystalline states are mapped out, and the nature of the transitions is determined. The liquid to glass transition changes from first to second order as the strength of the disorder is increased. For weak disorder, the system undergoes a first-order crystallization transition as the density is increased. Beyond a critical value of the disorder strength, this transition is replaced by a continuous glass transition. Our numerical results are compared with those of analytical work on the same system. Implications of our results for the field-temperature phase diagram of type-II superconductors are discussed.

Suboptimal reentry guidance of a reusable launch vehicle using pitch plane maneuver

Using the recently developed model predictive static programming (MPSP) technique, a nonlinear suboptimal reentry guidance scheme is presented in this paper for a reusable launch vehicle (RLV). Unlike traditional RLV guidance, the problem considered over here is restricted only to pitch plane maneuver of the vehicle, which allows simpler mission planning and vehicle load management. The computationally efficient MPSP technique brings in the philosophy of trajectory optimization into the framework of guidance design, which in turn results in very effective guidance schemes in general. In the problem addressed in this paper, it successfully guides the RLV through the critical reentry phase both by constraining it to the allowable narrow flight corridor as well as by meeting the terminal constraints at the end of the reentry segment. The guidance design is validated by considering possible aerodynamic uncertainties as well as dispersions in the initial conditions. (C) 2010 Elsevier Masson SAS. All rights reserved.

Observations on transition in plane bubble plumes

Flow visualization studies of plane laminar bubble plumes have been conducted to yield quantitative data on transition height, wavelength and wave velocity of the most unstable disturbance leading to transition. These are believed to be the first results of this kind. Most earlier studies are restricted to turbulent bubble plumes. In the present study, the bubble plumes were generated by electrolysis of water and hence very fine control over bubble size distribution and gas flow rate was possible to enable studies with laminar bubble plumes. Present observations show that (a) the dominant mode of instability in plane bubble plumes is the sinuous mode, (b) transition height and wavelength are related linearly with the proportionality constant being about 4, (c) wave velocity is about 40 % of the mean plume velocity, and (d) normalized transition height data correlate very well with a source Grashof number. Some agreement and some differences in transition characteristics of bubble plumes have been observed compared to those for similar single-phase flows.