Time Scale Separated Nonlinear Partial Integrated Guidance and Control of Interceptors in the Terminal Phase

The paper proposes a time scale separated partial integrated guidance and control of an interceptor for engaging high speed targets in the terminal phase. In this two loop design, the outer loop is an optimal control formulation based on nonlinear model predictive spread control philosophies. It gives the commanded pitch and yaw rates whereas necessary roll-rate command is generated from a roll-stabilization loop. The inner loop tracks the outer loop commands using the dynamicinversion philosophy. However, unlike conventional designs, in both the loops the Six degree of freedom (Six-DOF) interceptor model is used directly. This intelligent manipulation preserves the inherent time scale separation property between the translational and rotational dynamics, and hence overcomes the deficiency of current IGC designs, while preserving its benefits. Six-DOF simulation studies have been carried out accounting for three dimensional engagement geometry. Different comparison studies were also conducted to measure the performance of the algorithm.

Nonlinear Sliding Mode Control Design for an Air-breathing Engine with State Estimation

A modern system theory based nonlinear control design is discussed in this paper for successful operation of an air-breathing engine operating at supersonic speed. The primary objective of the control design of such an air-breathing engine is to ensure that the engine dynamically produces the thrust that tracks a commanded value of thrust as closely as possible by regulating the fuel flow to the combustion system. However, since the engine operates in the supersonic range, an important secondary objective is to manage the shock wave configuration in the intake section of the engine which is manipulated by varying the throat area of the nozzle. A nonlinear sliding mode control technique has been successfully used to achieve both of the above objectives. In this problem, since the process is faster than the actuators, independent control designs are also carried out for the actuators as well to assure the satisfactory performance of the system. Moreover, to filter out the sensor and process noises and to estimate the states for making the control design operate based on output feedback, an Extended Kalman Filter based state estimation design is also carried out. The promising simulation results suggest that the proposed control design approach is quite successful in obtaining robust performance of the air-breathing engine.

Nonlinear intrinsic instability of solid propellant combustion including gas-phase thermal inertia

The problem of homogeneous solid propellant combustion instability is studied with a one-dimensional flame model, including the effects of gas-phase thermal inertia and nonlinearity. Computational results presented in this paper show nonlinear instabilities inherent in the equations, due to which periodic burning is found even under steady ambient conditions such as pressure. The stability boundary is obtained in terms of Denison-Baum parameters. It is found that inclusion of gas-phase thermal inertia stabilizes the combustion. Also, the effect of a distributed heat release in the gas phase, compared to the flame sheet model, is to destabilize the burning. Direct calculations for finite amplitude pressure disturbances show that two distinct resonant modes exist, the first one near the natural frequency as obtained from intrinsic instability analysis and a second mode occurring at a much higher driving frequency. It is found that er rn in the low frequency region, the response of the propellant is significantly affected by the specific type of gas-phase chemical heat-release model employed. Examination of frequency response function reveals that the role of gas-phase thermal inertia is to stabilize the burning near the first resonant mode. Calculations made for different amplitudes of driving pressure show that the mean burning rate decreases with increasing amplitude. Also, with an increase in the driving amplitude, higher harmonics are generated in the burning rate.

Evolution and decay of cylindrical and spherical nonlinear acoustic waves generated by a sinusoidal source

The present work gives a comprehensive numerical study of the evolution and decay of cylindrical and spherical nonlinear acoustic waves generated by a sinusoidal source. Using pseudospectral and predictor–corrector implicit finite difference methods, we first reproduced the known analytic results of the plane harmonic problem to a high degree of accuracy. The non-planar harmonic problems, for which the amplitude decay is faster than that for the planar case, are then treated. The results are correlated with the known asymptotic results of Scott (1981) and Enflo (1985). The constant in the old-age formula for the cylindrical canonical problem is found to be 1.85 which is rather close to 2, ‘estimated’ analytically by Enflo. The old-age solutions exhibiting strict symmetry about the maximum are recovered; these provide an excellent analytic check on the numerical solutions. The evolution of the waves for different source geometries is depicted graphically.

Trajectory controllability of nonlinear integro-differential system

A stronger concept of complete (exact) controllability which we call Trajectory Controllability is introduced in this paper. We study the Trajectory Controllability of an abstract nonlinear integro-differential system in the finite and infinite dimensional space setting. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

Studies on nonisothermal solid state galvanic cells — effect of gradients on EMF

An expression for the EMF of a nonisothermal galvanic cell, with gradients in both temperature and chemical potential across a solid electrolyte, is derived based on the phenomenological equations of irreversible thermodynamics. The EMF of the nonisothermal cell can be written as a sum of the contributions from the chemical potential gradient and the EMF of a thermocell operating in the same temperature gradient but at unit activity of the neutral form of the migrating species. The validity of the derived equation is confirmed experimentally by imposing nonlinear gradients of temperature and chemical potential across galvanic cells constructed using fully stabilized zirconia as the electrolyte. The nature of the gradient has no effect on the EMF.

Helicopter blade flapping with and without small angle assumption in the presence of dynamic stall

The flapping equation for a rotating rigid helicopter blade is typically derived by considering (1)small flap angle, (2) small induced angle of attack and (3) linear aerodynamics. However, the use of nonlinear aerodynamics such as dynamic stall can make the assumptions of small angles suspect as shown in this paper. A general equation describing helicopter blade flap dynamics for large flap angle and large induced inflow angle of attack is derived. A semi-empirical dynamic stall aerodynamics model (ONERA model) is used. Numerical simulations are performed by solving the nonlinear flapping ordinary differential equation for steady state conditions and the validity of the small angle approximations are examined. It is shown that the small flapping assumption, and to a lesser extent, the small induced angle ofattack assumption, can lead to inaccurate predictions of the blade flap response in certain flight conditions for some rotors when nonlinear aerodynamics is considered. (C) 2010 Elsevier Inc. All rights reserved.

A simple technique for approximate solutions of the Falkner-Skan equation

A simple, sufficiently accurate and efficient method for approximate solutions of the Falkner-Skan equation is proposed here for a wide range of the pressure gradient parameter. The proposed approximate solutions are obtained utilising a known solution of another differential equation.

Nonlinear measures of heart rate time series: influence of posture and controlled breathing

In this study, we investigated measures of nonlinear dynamics and chaos theory in regards to heart rate variability in 27 normal control subjects in supine and standing postures, and 14 subjects in spontaneous and controlled breathing conditions. We examined minimum embedding dimension (MED), largest Lyapunov exponent (LLE) and measures of nonlinearity (NL) of heart rate time series. MED quantifies the system's complexity, LLE predictability and NL, a measure of deviation from linear processes. There was a significant decrease in complexity (P<0.00001), a decrease in predictability (P<0.00001) and an increase in nonlinearity (P=0.00001) during the change from supine to standing posture. Decrease in MED, and increases in NL score and LLE in standing posture appear to be partly due to an increase in sympathetic activity of the autonomous nervous system in standing posture. An improvement in predictability during controlled breathing appears to be due to the introduction of a periodic component. (C) 2000 published by Elsevier Science B.V.

Correlations and the Ring-Kinetic Equation in Dense Sheared Granular Flows

A formal way of deriving fluctuation-correlation relations in dense sheared granular media, starting with the Enskog approximation for the collision integral in the Chapman-Enskog theory, is discussed. The correlation correction to the viscosity is obtained using the ring-kinetic equation, in terms of the correlations in the hydrodynamic modes of the linearised Enskog equation. It is shown that the Green-Kubo formula for the shear viscosity emerges from the two-body correlation function obtained from the ring-kinetic equation.

Nonlinear I-V characteristics of nanocrystalline SnO2

Current versus voltage characteristics (I-V) of nanocrystalline SnO2 materials have been investigated in air at room temperature. The samples were prepared by the inert gas condensation technique (IGCT) as well as by chemical methods. X-ray diffraction studies showed a tetragonal rutile structure for all the samples. Microstructural studies were performed with transmission electron microscopy. All the samples exhibited nonlinear I-V characteristics of the current-controlled negative resistance (CCNR) type. The results show that the threshold field (break down) voltage is higher for the samples prepared by the IGCT method than for those prepared by the chemical method due to the formation of a tin oxide layer over the crystalline tin. It is also found that the threshold field increases with the decrease in grain size.

Fractional quantum Hall edge: Effect of nonlinear dispersion and edge roton

According to Wen's theory, a universal behavior of the fractional quantum Hall edge is expected at sufficiently low energies, where the dispersion of the elementary edge excitation is linear. A microscopic calculation shows that the actual dispersion is indeed linear at low energies, but deviates from linearity beyond certain energy, and also exhibits an "edge roton minimum." We determine the edge exponent from a microscopic approach, and find that the nonlinearity of the dispersion makes a surprisingly small correction to the edge exponent even at energies higher than the roton energy. We explain this insensitivity as arising from the fact that the energy at maximum spectral weight continues to show an almost linear behavior up to fairly high energies. We also study, in an effective-field theory, how interactions modify the exponent for a reconstructed edge with multiple edge modes. Relevance to experiment is discussed.

A note on the equivalence of shock manifold equations

The shock manifold equation is a first order nonlinear partial differential equation, which describes the kinematics of a shockfront in an ideal gas with constant specific heats. However, it was found that there was more than one of these shock manifold equations, and the shock surface could be embedded in a one parameter family of surfaces, obtained as a solution of any of these shock manifold equations. Associated with each shock manifold equation is a set of characteristic curves called lsquoshock raysrsquo. This paper investigates the nature of various associated shock ray equations.

An Application of 3-D Kinematical Conservation Laws: Propagation of a 3-D Wavefront

Three-dimensional (3-D) kinematical conservation laws (KCL) are equations of evolution of a propagating surface Omega(t) in three space dimensions. We start with a brief review of the 3-D KCL system and mention some of its properties relevant to this paper. The 3-D KCL, a system of six conservation laws, is an underdetermined system to which we add an energy transport equation for a small amplitude 3-D nonlinear wavefront propagating in a polytropic gas in a uniform state and at rest. We call the enlarged system of 3-D KCL with the energy transport equation equations of weakly nonlinear ray theory (WNLRT). We highlight some interesting properties of the eigenstructure of the equations of WNLRT, but the main aim of this paper is to test the numerical efficacy of this system of seven conservation laws. We take several initial shapes for a nonlinear wavefront with a suitable amplitude distribution on it and let it evolve according to the 3-D WNLRT. The 3-D WNLRT is a weakly hyperbolic 7 x 7 system that is highly nonlinear. Here we use the staggered Lax-Friedrichs and Nessyahu-Tadmor central schemes and have obtained some very interesting shapes of the wavefronts. We find the 3-D KCL to be suitable for solving many complex problems for which there presently seems to be no other method capable of giving such physically realistic features.

Distributed nonlinear optimization algorithms for general network problems

There are a number of large networks which occur in many problems dealing with the flow of power, communication signals, water, gas, transportable goods, etc. Both design and planning of these networks involve optimization problems. The first part of this paper introduces the common characteristics of a nonlinear network (the network may be linear, the objective function may be non linear, or both may be nonlinear). The second part develops a mathematical model trying to put together some important constraints based on the abstraction for a general network. The third part deals with solution procedures; it converts the network to a matrix based system of equations, gives the characteristics of the matrix and suggests two solution procedures, one of them being a new one. The fourth part handles spatially distributed networks and evolves a number of decomposition techniques so that we can solve the problem with the help of a distributed computer system. Algorithms for parallel processors and spatially distributed systems have been described.There are a number of common features that pertain to networks. A network consists of a set of nodes and arcs. In addition at every node, there is a possibility of an input (like power, water, message, goods etc) or an output or none. Normally, the network equations describe the flows amoungst nodes through the arcs. These network equations couple variables associated with nodes. Invariably, variables pertaining to arcs are constants; the result required will be flows through the arcs. To solve the normal base problem, we are given input flows at nodes, output flows at nodes and certain physical constraints on other variables at nodes and we should find out the flows through the network (variables at nodes will be referred to as across variables).The optimization problem involves in selecting inputs at nodes so as to optimise an objective function; the objective may be a cost function based on the inputs to be minimised or a loss function or an efficiency function. The above mathematical model can be solved using Lagrange Multiplier technique since the equalities are strong compared to inequalities. The Lagrange multiplier technique divides the solution procedure into two stages per iteration. Stage one calculates the problem variables % and stage two the multipliers lambda. It is shown that the Jacobian matrix used in stage one (for solving a nonlinear system of necessary conditions) occurs in the stage two also.A second solution procedure has also been imbedded into the first one. This is called total residue approach. It changes the equality constraints so that we can get faster convergence of the iterations.Both solution procedures are found to coverge in 3 to 7 iterations for a sample network.The availability of distributed computer systems — both LAN and WAN — suggest the need for algorithms to solve the optimization problems. Two types of algorithms have been proposed — one based on the physics of the network and the other on the property of the Jacobian matrix. Three algorithms have been deviced, one of them for the local area case. These algorithms are called as regional distributed algorithm, hierarchical regional distributed algorithm (both using the physics properties of the network), and locally distributed algorithm (a multiprocessor based approach with a local area network configuration). The approach used was to define an algorithm that is faster and uses minimum communications. These algorithms are found to converge at the same rate as the non distributed (unitary) case.