Bispectra of A Tropical Coupled Ocean-Atmosphere System

It has recently been proposed that the broad spectrum of interannual variability in the tropics with a peak around four years results from an interaction between the linear low-frequency oscillatory mode of the coupled system and the nonlinear higher-frequency modes of the system. In this study we determine the bispectrum of the conceptual model consisting of a nonlinear low-order model coupled to a linear oscillator for various values of the coupling constants.

Linear dynamics and active control of an elastically supported travelling string

In this paper, the linear dynamics and active control of a string travelling with uniform velocity is presented. Discrete elastic supports are introduced along the length of the string. Finite element formulation is adopted to obtain the governing equations of motion. The velocity of translation introduces gyroscopic terms in the system equations. The effect of translation and the discrete elastic supports on the free vibration solution is studied. The solution is utilized in actively controlling the string vibrations due to an initial disturbance. The control, affected in modal space, is optimal with respect to a quadratic performance index. Numerical results are presented to demonstrate the effectiveness of the control strategy in regulating the travelling string vibrations.

Lyapunov exponents and predictability of the tropical coupled ocean-atmosphere system

It has recently been proposed that the broad spectrum of interannual variability in the tropics with a peak around four years results from an interaction between the linear low-frequency oscillatory mode of the coupled system and the nonlinear higher-frequency modes of the system. In this study we determine the Lyapunov exponents of the conceptual model consisting of a nonlinear low-order model coupled to a linear oscillator for various values of the coupling constants.

On evaluating erosion by sand particles in polythene system without and with ceramic particles

The effect of the inclusion of ceramic particles in polythene material on the response to erosion due to impingement by sand particles at three angles is investigated. It is seen that erosion resistance varies with ceramic inclusions. The work also considers the limitations posed by the system in adopting weight change measurements as a measure to follow erosive wear owing to the softer nature of the matrix material. Consequently, the investigation looks at two other experimental parameter, that can readily be measured to quantify erosion. Of the two approaches. the advantages of following wear through measuring linear dimension of the resulting crater is stressed in this work. The study also highlights the problems associated in assessing the depth of the crater as a parameter to express the extent of erosion owing to the phenomenon of material flow suggested and schematically illustrated in the work. Corroborative evidence for this flow behaviour through scanning electron microscopic studies is presented. (C) 2002 Elsevier Science Ltd. All rights reserved.

Phase diagram of a two-species lattice model with a linear instability

We discuss the properties of a one-dimensional lattice model of a driven system with two species of particles in which the mobility of one species depends on the density of the other. This model was introduced by Lahiri and Ramaswamy (Phys. Rev. Lett., 79, 1150 (1997)) in the context of sedimenting colloidal crystals, and its continuum version was shown to exhibit an instability arising from linear gradient couplings. In this paper we review recent progress in understanding the full phase diagram of the model. There are three phases. In the first, the steady state can be determined exactly along a representative locus using the condition of detailed balance. The system shows phase separation of an exceptionally robust sort, termed strong phase separation, which survives at all temperatures. The second phase arises in the threshold case where the first species evolves independently of the second, but the fluctuations of the first influence the evolution of the second, as in the passive scalar problem. The second species then shows phase separation of a delicate sort, in which long-range order coexists with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase ordering is associated with power law decays in cluster size distributions and a breakdown of the Porod law. The third phase is one with a uniform overall density, and along a representative locus the steady state is shown to have product measure form. Density fluctuations are transported by two kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation properties are governed by the symmetries of these couplings, which depend on the overall densities. In the most interesting case,, the dissipation of the two modes is characterized by different critical exponents, despite the nonlinear coupling.

Decentralized control of linear multivariable systems

A strongly connected decentralized control system may be made single channel controllable and observable with respect to any channel by decentralized feedbacks. It is noted here that the system example considered by Corfmat and Morse to illustrate this fact is already single channel controllable and observable, with respect to one of the channels. An alternate example which fits into the situation is presented in this item.

Fault-Tolerant Average Execution Time Optimization for General Purpose Multi-Processor System-on Chips

Fault-tolerance is due to the semiconductor technology development important, not only for safety-critical systems but also for general-purpose (non-safety critical) systems. However, instead of guaranteeing that deadlines always are met, it is for general-purpose systems important to minimize the average execution time (AET) while ensuring fault-tolerance. For a given job and a soft (transient) error probability, we define mathematical formulas for AET that includes bus communication overhead for both voting (active replication) and rollback-recovery with checkpointing (RRC). And, for a given multi-processor system-on-chip (MPSoC), we define integer linear programming (ILP) models that minimize AET including bus communication overhead when: (1) selecting the number of checkpoints when using RRC, (2) finding the number of processors and job-to-processor assignment when using voting, and (3) defining fault-tolerance scheme (voting or RRC) per job and defining its usage for each job. Experiments demonstrate significant savings in AET.

Analysis of an LMS linear equalizer for fading channels in decision directed mode

We consider a time varying wireless fading channel, equalized by an LMS linear equalizer in decision directed mode (DD-LMS-LE). We study how well this equalizer tracks the optimal Wiener equalizer. Initially we study a fixed channel.For a fixed channel, we obtain the existence of DD attractors near the Wiener filter at high SNRs using an ODE (Ordinary Differential Equation) approximating the DD-LMS-LE. We also show, via examples, that the DD attractors may not be close to the Wiener filters at low SNRs. Next we study a time varying fading channel modeled by an Auto-regressive (AR) process of order 2. The DD-LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs. We show via examples that the LMS equalizer ODE show tracks the ODE corresponding to the instantaneous Wiener filter when the SNR is high. This may not happen at low SNRs.

Tracking performance of an LMS-Linear Equalizer for fading channels

We consider a time varying wireless fading channel, equalized by an LMS linear equalizer. We study how well this equalizer tracks the optimal Wiener equalizer. We model the channel by an Auto-regressive (AR) process. Then the LMS equalizer and the AR process are jointly approximated by the solution of a system of ODEs (ordinary differential equations). Using these ODEs, the error between the LMS equalizer and the instantaneous Wiener filter is shown to decay exponentially/polynomially to zero unless the channel is marginally stable in which case the convergence may not hold.Using the same ODEs, we also show that the corresponding Mean Square Error (MSE) converges towards minimum MSE(MMSE) at the same rate for a stable channel. We further show that the difference between the MSE and the MMSE does not explode with time even when the channel is unstable. Finally we obtain an optimum step size for the linear equalizer in terms of the AR parameters, whenever the error decay is exponential.

Solute-solute and solute-solvent interactions in transition metal alloys: Pt-Ti system

An isothermal section of the phase diagram for (silver + rhodium + oxygen) at T = 1173 K has been established by equilibration of samples representing twelve different compositions, and phase identification after quenching by optical and scanning electron microscopy (s.e.m.), X-ray diffraction (x.r.d.), and energy dispersive analysis of X-rays (e.d.x.), Only one ternary oxide, AgRhO2, was found to be stable and a three phase region involving Ag, AgRhO2 and Rh2O3 was identified. The thermodynamic properties of AgRhO2 were measured using a galvanic cell in the temperature range 980 K to 1320 K. Yttria-stabilized zirconia was used as the solid electrolyte and pure oxygen gas at a pressure of 0.1 MPa was used as the reference electrode. The Gibbs free energy of formation of the ternary oxide from the elements, ΔfGo (AgRhO2), can be represented by two linear equations that join at the melting temperature of silver. In the temperature range 980 K to 1235 K, ΔfGo(AgRhO2)/(J . mol-1) = -249080 + 179.08 T/K (±120). Above the melting temperature of silver, in the temperature range 1235 K to 1320 K, ΔfGo(AgRhO2)/(J . mol-1) = -260400 + 188.24 T/K (±95). The thermodynamic properties of AgRhO2 at T = 298.15 K were evaluated from the high temperature data. The chemical potential diagram for (silver + rhodium + oxygen) at T = 1200 K was also computed on the basis of the results of this study.

Accelerating Numerical Linear Algebra Kernels on a Scalable Run Time Reconfigurable Platform

Numerical Linear Algebra (NLA) kernels are at the heart of all computational problems. These kernels require hardware acceleration for increased throughput. NLA Solvers for dense and sparse matrices differ in the way the matrices are stored and operated upon although they exhibit similar computational properties. While ASIC solutions for NLA Solvers can deliver high performance, they are not scalable, and hence are not commercially viable. In this paper, we show how NLA kernels can be accelerated on REDEFINE, a scalable runtime reconfigurable hardware platform. Compared to a software implementation, Direct Solver (Modified Faddeev's algorithm) on REDEFINE shows a 29X improvement on an average and Iterative Solver (Conjugate Gradient algorithm) shows a 15-20% improvement. We further show that solution on REDEFINE is scalable over larger problem sizes without any notable degradation in performance.

Matrix Processors Using p-adic Arithmetic for Exact Linear Computations

A unique code (called Hensel's code) is derived for a rational number by truncating its infinite p-adic expansion. The four basic arithmetic algorithms for these codes are described and their application to rational matrix computations is demonstrated by solving a system of linear equations exactly, using the Gaussian elimination procedure.

2.5PN linear momentum flux from inspiralling compact binaries in quasicircular orbits and associated recoil: Nonspinning case

Anisotropic emission of gravitational waves (GWs) from inspiralling compact binaries leads to the loss of linear momentum and hence gravitational recoil of the system. The loss rate of linear momentum in the far-zone of the source (a nonspinning binary system of black holes in quasicircular orbit) is investigated at the 2.5 post-Newtonian (PN) order and used to provide an analytical expression in harmonic coordinates for the 2.5PN accurate recoil velocity of the binary accumulated in the inspiral phase. The maximum recoil velocity of the binary system at the end of its inspiral phase (i.e at the innermost stable circular orbit (ISCO)) estimated by the 2.5PN formula is of the order of 4 km s(-1) which is smaller than the 2PN estimate of 22 km s(-1). Going beyond inspiral, we also provide an estimate of the more important contribution to the recoil velocity from the plunge phase. The maximum recoil velocity at the end of the plunge, involving contributions both from inspiral and plunge phase, for a binary with symmetric mass ratio nu = 0.2 is of the order of 182 km s(-1).

Radial infall of two compact objects: 2.5-post-Newtonian linear momentum flux and associated recoil

The loss rate of linear momentum from a binary system composed of compact objects (radially falling towards each other under mutual gravitational influence) has been investigated using the multipolar post-Minkowskian approach. The 2.5PN accurate analytical formula for the linear momentum flux is provided, in terms of the separation of the two objects, in harmonic coordinates, both for a finite and an infinite initial separation. The 2.5PN formulas for the linear momentum flux are finally used to estimate the recoil velocity accumulated during a premerger phase of the binary evolution.

Linearization and Reconstruction of Non-linear Diffuse Optical Tomographic Image

Diffuse optical tomography (DOT) is one of the ways to probe highly scattering media such as tissue using low-energy near infra-red light (NIR) to reconstruct a map of the optical property distribution. The interaction of the photons in biological tissue is a non-linear process and the phton transport through the tissue is modelled using diffusion theory. The inversion problem is often solved through iterative methods based on nonlinear optimization for the minimization of a data-model misfit function. The solution of the non-linear problem can be improved by modeling and optimizing the cost functional. The cost functional is f(x) = x(T)Ax - b(T)x + c and after minimization, the cost functional reduces to Ax = b. The spatial distribution of optical parameter can be obtained by solving the above equation iteratively for x. As the problem is non-linear, ill-posed and ill-conditioned, there will be an error or correction term for x at each iteration. A linearization strategy is proposed for the solution of the nonlinear ill-posed inverse problem by linear combination of system matrix and error in solution. By propagating the error (e) information (obtained from previous iteration) to the minimization function f(x), we can rewrite the minimization function as f(x; e) = (x + e)(T) A(x + e) - b(T)(x + e) + c. The revised cost functional is f(x; e) = f(x) + e(T)Ae. The self guided spatial weighted prior (e(T)Ae) error (e, error in estimating x) information along the principal nodes facilitates a well resolved dominant solution over the region of interest. The local minimization reduces the spreading of inclusion and removes the side lobes, thereby improving the contrast, localization and resolution of reconstructed image which has not been possible with conventional linear and regularization algorithm.