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.

A Numerical Scheme for Three-Dimensional Front Propagation and Control of Jordan Mode

As an example of a front propagation, we study the propagation of a three-dimensional nonlinear wavefront into a polytropic gas in a uniform state and at rest. The successive positions and geometry of the wavefront are obtained by solving the conservation form of equations of a weakly nonlinear ray theory. The proposed set of equations forms a weakly hyperbolic system of seven conservation laws with an additional vector constraint, each of whose components is a divergence-free condition. This constraint is an involution for the system of conservation laws, and it is termed a geometric solenoidal constraint. The analysis of a Cauchy problem for the linearized system shows that when this constraint is satisfied initially, the solution does not exhibit any Jordan mode. For the numerical simulation of the conservation laws we employ a high resolution central scheme. The second order accuracy of the scheme is achieved by using MUSCL-type reconstructions and Runge-Kutta time discretizations. A constrained transport-type technique is used to enforce the geometric solenoidal constraint. The results of several numerical experiments are presented, which confirm the efficiency and robustness of the proposed numerical method and the control of the Jordan mode.

On methods of calculating successive positions of a shock front

In this paper we have discussed limits of the validity of Whitham's characteristic rule for finding successive positions of a shock in one space dimension. We start with an example for which the exact solution is known and show that the characteristic rule gives correct result only if the state behind the shock is uniform. Then we take the gas dynamic equations in two cases: one of a shock propagating through a stratified layer and other down a nonuniform tube and derive exact equations for the evolution of the shock amplitude along a shock path. These exact results are then compared with the results obtained by the characteristic rule. The characteristic rule not only incorrectly accounts for the deviation of the state behind the shock from a uniform state but also gives a coefficient in the equation which differ significantly from the exact coefficients for a wide range of values of the shock strength.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear flow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear how of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

The effect of boundaries on the plane Couette flow of granular materials: a bifurcation analysis

The tendency of granular materials in rapid shear ow to form non-uniform structures is well documented in the literature. Through a linear stability analysis of the solution of continuum equations for rapid shear flow of a uniform granular material, performed by Savage (1992) and others subsequently, it has been shown that an infinite plane shearing motion may be unstable in the Lyapunov sense, provided the mean volume fraction of particles is above a critical value. This instability leads to the formation of alternating layers of high and low particle concentrations oriented parallel to the plane of shear. Computer simulations, on the other hand, reveal that non-uniform structures are possible even when the mean volume fraction of particles is small. In the present study, we have examined the structure of fully developed layered solutions, by making use of numerical continuation techniques and bifurcation theory. It is shown that the continuum equations do predict the existence of layered solutions of high amplitude even when the uniform state is linearly stable. An analysis of the effect of bounding walls on the bifurcation structure reveals that the nature of the wall boundary conditions plays a pivotal role in selecting that branch of non-uniform solutions which emerges as the primary branch. This demonstrates unequivocally that the results on the stability of bounded shear flow of granular materials presented previously by Wang et al. (1996) are, in general, based on erroneous base states.

Optimal infinite-horizon feedback laws for a general class of constrained discrete-time systems: stability and moving-horizon approximations

Stability results are given for a class of feedback systems arising from the regulation of time-varying discrete-time systems using optimal infinite-horizon and moving-horizon feedback laws. The class is characterized by joint constraints on the state and the control, a general nonlinear cost function and nonlinear equations of motion possessing two special properties. It is shown that weak conditions on the cost function and the constraints are sufficient to guarantee uniform asymptotic stability of both the optimal infinite-horizon and movinghorizon feedback systems. The infinite-horizon cost associated with the moving-horizon feedback law approaches the optimal infinite-horizon cost as the moving horizon is extended.

Effect of inter-edge Coulomb interactions on transport through a point contact in nu=5/2 quantum Hall state

We study transport across a point contact separating two line junctions in a nu = 5/2 quantum Hall system. We analyze the effect of inter-edge Coulomb interactions between the chiral bosonic edge modes of the half-filled Landau level (assuming a Pfaffian wave function for the half-filled state) and of the two fully filled Landau levels. In the presence of inter-edge Coulomb interactions between all the six edges participating in the line junction, we show that the stable fixed point corresponds to a point contact that is neither fully opaque nor fully transparent. Remarkably, this fixed point represents a situation where the half-filled level is fully transmitting, while the two filled levels are completely backscattered; hence the fixed point Hall conductance is given by G(H) = 1/2e(2)/h. We predict the non-universal temperature power laws by which the system approaches the stable fixed point from the two unstable fixed points corresponding to the fully connected case (G(H) = 5/2e(2)/h) and the fully disconnected case (G(H) = 0).

On plane-wave propagation in a uniform pipe in the presence of a mean flow and a temperature gradient

A theoretical study on the propagation of plane waves in the presence of a hot mean flow in a uniform pipe is presented. The temperature variation in the pipe is taken to be a linear temperature gradient along the axis. The theoretical studies include the formulation of a wave equation based on continuity, momentum, and state equation, and derivation of a general four-pole matrix, which is shown to yield the well-known transfer matrices for several other simpler cases.

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.

Elastic behaviour of matter under very high pressures - Uniform compression

In order to study the elastic behaviour of matter when subjected to very large pressures, such as occur for example in the interior of the earth, and to provide an explanation for phenomena like earthquakes, it is essential to be able to calculate the values of the elastic constants of a substance under a state of large initial stress in terms of the elastic constants of a natural or stress-free state. An attempt has been made in this paper to derive expressions for these quantities for a substance of cubic symmetry on the basis of non-linear theory of elasticity and including up to cubic powers of the strain components in the strain energy function. A simple method of deriving them directly from the energy function itself has been indicated for any general case and the same has been applied to the case of hydrostatic compression. The notion of an effective elastic energy-the energy require to effect an infinitesimal deformation over a state of finite strain-has been introduced, the coefficients in this expression being the effective elastic constants. A separation of this effective energy function into normal co-ordinates has been given for the particular case of cubic symmetry and it has been pointed out, that when any of such coefficients in this normal form becomes negative, elastic instability will set in, with associated release of energy.

Impedance parameters and the state-of-charge. II. Lead-acid battery

The determination of the state-of-charge of the lead-acid battery has been examined from the viewpoint of internal impedance. It is shown that the impedance is controlled by charge transfer and to a smaller extent by diffusion processes in the frequency range 15–100 Hz. The equivalent series/parallel capacitance as well as the a.c. phase-shift show a parabolic dependence upon the state-of-charge, with a maximum or minimum at 50% charge. These results are explained on the basis of a uniform transmission-line analog equivalent circuit for the battery electrodes.

Impedance parameters and the state-of charge. I. Nickel-cadmium battery

The problem of nondestructive determination of the state-of-charge of nickel-cadmium batteries has been examined experimentally as well as theoretically from the viewpoint of internal impedance. It is shown that the modulus of the impedance is mainly controlled by diffusion at all states of charge. Even so, a prediction of the state of charge is possible if the equivalent series/parallel capacitance or the alternating current phase shift is measured at a sufficiently low a.c. test frequency (5–30 Hz) which also avoids inductive effects. These results are explained on the basis of a uniform transmission-line analog equivalent circuit for the battery electrodes.

New uniform solid catalyst for the low-temperature oxidation of carbon monoxide: A triple-layered rare earth perovskite containing Co and Cu ions

Oxygen reactivity and catalytic activity of the cobalt-containing layered defect perovskites, YBa2Cu2CoO7+delta and LaBa2Cu2CoO7+delta, in comparison with LaBa2Cu3O7-delta have been investigated employing temperature-programmed desorption (TPD) and temperature-programmed surface reactions (TPSR) in the stoichiometric and catalytic mode using carbon monoxide as a probe molecule. TPD studies showed evidence for the presence of two distinct labile oxygen species, one at (0 0 1/2) sites and the other at (0 1/2 0) sites in LaBa2Cu2CoO7+delta against a single labile species at (0 1/2 0) in the case of two other oxides. The activation energies for the catalytic oxidation of carbon monoxide by oxygen over LaBa2Cu3O7-delta, YBa2Cu2CoO7+delta, and LaBa2Cu2CoO7+delta have been estimated to be 24.2, 15.9, and 13.6 kcal/mol, respectively. The reactivity and catalytic activity of the oxide systems have been interpreted in terms of the structural changes brought about by substituents, guided by a directing effect of the larger rare earth cation. TPSR profiles, structural analysis, and infrared spectroscopic investigations suggest that the oxygen present at (0 0 1/2) sites in the case of LaBa2Cu2CoO7+delta is accessible to catalytic oxidation of CO through a Mars-Van Krevelen pathway. Catalytic conversion of CO to CO2 over LaBa2Cu2CoO7+delta occurs at 200 degrees C. The enhanced reactivity is explained in terms of changes brought about in the coordination polyhedra around transition metals, enhanced basal plane oxygen diffusivity, and redox potentials of the different transition metal cations.

Friction Stir-Welded Dissimilar Aluminum Alloys: Microstructure, Mechanical Properties, and Physical State

A356 and 6061 aluminum alloys were joined by friction stir welding at constant tool rotational rate with different tool-traversing speeds. Thermomechanical data of welding showed that increment in tool speed reduced the pseudo heat index and temperature at weld nugget (WN). On the other hand, volume of material within extrusion zone, strain rate, and Zenner Hollomon parameter were reduced with decrease in tool speed. Optical microstructure of WN exhibited nearly uniform dispersion of Si-rich particles, fine grain size of 6061 Al alloy, and disappearance of second phase within 6061 Al alloy. With enhancement in welding speed, matrix grain size became finer, yet size of Si-rich particles did not reduce incessantly. Size of Si-rich particles was governed by interaction time between tool and substrate. Mechanical property of WN was evaluated. It has been found that the maximum joint efficiency of 116% with respect to that of 6061 alloy was obtained at an intermediate tool-traversing speed, where matrix grain size was significantly fine and those of Si-rich particles were substantially small.

Zero-sum risk-sensitive stochastic games on a countable state space

Infinite horizon discounted-cost and ergodic-cost risk-sensitive zero-sum stochastic games for controlled Markov chains with countably many states are analyzed. Upper and lower values for these games are established. The existence of value and saddle-point equilibria in the class of Markov strategies is proved for the discounted-cost game. The existence of value and saddle-point equilibria in the class of stationary strategies is proved under the uniform ergodicity condition for the ergodic-cost game. The value of the ergodic-cost game happens to be the product of the inverse of the risk-sensitivity factor and the logarithm of the common Perron-Frobenius eigenvalue of the associated controlled nonlinear kernels. (C) 2013 Elsevier B.V. All rights reserved.