A simplified kinetic model for oxidative dehydrogenation of ethylbenzene over Pd-NaBr/Al2O3 catalyst

The oxidative dehydrogenation of ethylbenzene is gaining considerable importance in recent years as a promising alternative for styrene production. This vapour phase reaction has been studied over Pd-NaBr/Al2O3 catalyst in the temperature range 623-793 K in a fixed bed reactor. Kinetic analysis of this reaction has been done using a recursion procedure developed in this work from first principles. The advantage of this method is the absence of any restriction on the conversion level as it uses an integrated rate equation. The rate of styrene formation was found to follow a linear relationship with concentration of ethylbenzene and shows a Langmuir type dependence on the concentration of oxygen.

Matrix method for eigenstructure assignment - The multi-input case with application

The eigenvalue and eigenstructure assignment procedure has found application in a wide variety of control problems. In this paper a method for assigning eigenstructure to a linear time invariant multi-input system is proposed. The algorithm determines a matrix that has eigenvalues and eigenvectors at the desired locations. It is obtained from the knowledge of the open-loop system and the desired eigenstructure. Solution of the matrix equation, involving unknown controller gams, open-loop system matrices, and desired eigenvalues and eigenvectors, results hi the state feedback controller. The proposed algorithm requires the closed-loop eigenvalues to be different from those of the open-loop case. This apparent constraint can easily be overcome by a negligible shift in the values. Application of the procedure is illustrated through the offset control of a satellite supported, from an orbiting platform, by a flexible tether.

Matrix method for eigenstructure assignment - The multi-input case with application

The eigenvalue and eigenstructure assignment procedure has found application in a wide variety of control problems. In this paper a method for assigning eigenstructure to a Linear time invariant multi-input system is proposed. The algorithm determines a matrix that has eigenvalues and eigenvectors at the desired locations. It is obtained from the knowledge of the open-loop system and the desired eigenstructure. solution of the matrix equation, involving unknown controller gains, open-loop system matrices, and desired eigenvalues and eigenvectors, results in the state feedback controller. The proposed algorithm requires the closed-loop eigenvalues to be different from those of the open-loop case. This apparent constraint can easily be overcome by a negligible shift in the values. Application of the procedure is illustrated through the offset control of a satellite supported, from an orbiting platform, by a flexible tether,

Low-Complexity Detection and Performance in Multi-Gigabit High Spectral Efficiency Wireless Systems

Recently, we reported a low-complexity likelihood ascent search (LAS) detection algorithm for large MIMO systems with several tens of antennas that can achieve high spectral efficiencies of the order of tens to hundreds of bps/Hz. Through simulations, we showed that this algorithm achieves increasingly near SISO AWGN performance for increasing number of antennas in Lid. Rayleigh fading. However, no bit error performance analysis of the algorithm was reported. In this paper, we extend our work on this low-complexity large MIMO detector in two directions: i) We report an asymptotic bit error probability analysis of the LAS algorithm in the large system limit, where N-t, N-r -> infinity keeping N-t = N-r, where N-t and N-r are the number of transmit and receive antennas, respectively. Specifically, we prove that the error performance of the LAS detector for V-BLAST with 4-QAM in i.i.d. Rayleigh fading converges to that of the maximum-likelihood (ML) detector as N-t, N-r -> infinity keeping N-t = N-r ii) We present simulated BER and nearness to capacity results for V-BLAST as well as high-rate non-orthogonal STBC from Division Algebras (DA), in a more realistic spatially correlated MIMO channel model. Our simulation results show that a) at an uncoded BER of 10(-3), the performance of the LAS detector in decoding 16 x 16 STBC from DA with N-t = = 16 and 16-QAM degrades in spatially correlated fading by about 7 dB compared to that in i.i.d. fading, and 19) with a rate-3/4 outer turbo code and 48 bps/Hz spectral efficiency, the performance degrades by about 6 dB at a coded BER of 10(-4). Our results further show that providing asymmetry in number of antennas such that N-r > N-t keeping the total receiver array length same as that for N-r = N-t, the detector is able to pick up the extra receive diversity thereby significantly improving the BER performance.

Conjectures about phase turbulence in the complex Ginzburg-Landau equation

In the complex Ginzburg-Landau equation, we consider possible ''phase turbulent'' regimes, where asymptotic correlations are controlled by phase fluctuations rather than by topological defects. Conjecturing that the decay of such correlations is governed by the Kardar-Parisi-Zhang (KPZ) model of growing interfaces, we derive the following results: (1) A scaling ansatz implies that equal-time spatial correlations in 1d, 2d, and 3d decay like e(-Ax2 zeta), where A is a nonuniversal constant, and zeta=1/2 in 1d. (2) Temporal correlations decay as exp(-t(2 beta)h(t/L(z))), with the scaling law <(beta)over bar> = <(zeta)over bar>/z, where z = 3/2, 1.58..., and 1.66..., for d = 1,2, and 3 respectively. The scaling function h(y) approaches a constant as y --> 0, and behaves like y(2(beta-<(beta)over bar>)), for large y. If in 3d the associated KPZ model turns out to be in its weak-coupling (''smooth'') phase, then, instead of the above behavior, the CGLE exhibits rotating long-range order whose connected correlations decay like 1/x in space or 1/t(1/2) in time. (3) For system sizes, L, and times t respectively less than a crossover length, L(c), and time, t(c), correlations are governed by the free-field or Edwards-Wilkinson (EW) equation, rather than the KPZ model. In 1d, we find that L(c) is large: L(c) similar to 35,000; for L < L(c) we show numerical evidence for stretched exponential decay of temporal correlations with an exponent consistent with the EW value beta(EW)= 1/4.

The mean and turbulence structure simulation of the monsoon trough boundary layer using a one-dimensional model with e-l and e-epsilon closures

An attempt has been made here to study the sensitivity of the mean and the turbulence structure of the monsoon trough boundary layer to the choice of the constants in the dissipation equation for two stations Delhi and Calcutta, using one-dimensional atmospheric boundary layer model with e-epsilon turbulence closure. An analytical discussion of the problems associated with the constants of the dissipation equation is presented. It is shown here that the choice of the constants in the dissipation equation is quite crucial and the turbulence structure is very sensitive to these constants. The modification of the dissipation equation adopted by earlier studies, that is, approximating the Tke generation (due to shear and buoyancy production) in the epsilon-equation by max (shear production, shear + buoyancy production), can be avoided by a suitable choice of the constants suggested here. The observed turbulence structure is better simulated with these constants. The turbulence structure simulation with the constants recommended by Aupoix et al (1989) (which are interactive in time) for the monsoon region is shown to be qualitatively similar to the simulation obtained with the constants suggested here, thus implying that no universal constants exist to regulate dissipation rate. Simulations of the mean structure show little sensitivity to the type of the closure parameterization between e-l and e-epsilon closures. However the turbulence structure simulation with e-epsilon closure is far better compared to the e-l model simulations. The model simulations of temperature profiles compare quite well with the observations whenever the boundary layer is well mixed (neutral) or unstable. However the models are not able to simulate the nocturnal boundary layer (stable) temperature profiles. Moisture profiles are simulated reasonably better. With one-dimensional models, capturing observed wind variations is not up to the mark.

Nonequilibrium phase transition in the kinetic Ising model: Critical slowing down and the specific-heat singularity

The nonequilibrium dynamic phase transition, in the kinetic Ising model in the presence of an oscillating magnetic field has been studied both by Monte Carlo simulation and by solving numerically the mean-field dynamic equation of motion for the average magnetization. In both cases, the Debye ''relaxation'' behavior of the dynamic order parameter has been observed and the ''relaxation time'' is found to diverge near the dynamic transition point. The Debye relaxation of the dynamic order parameter and the power law divergence of the relaxation time have been obtained from a very approximate solution of the mean-field dynamic equation. The temperature variation of appropriately defined ''specific heat'' is studied by the Monte Carlo simulation near the transition point. The specific heat has been observed to diverge near the dynamic transition point.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

Accelerating multi-core simulators

Simulation is an important means of evaluating new microarchitectures. With the invention of multi-core (CMP) platforms, simulators are becoming larger and more complex. However, with the availability of CMPs with larger caches and higher operating frequency, the wall clock time required for simulating an application has become comparatively shorter. Reducing this simulation time further is a great challenge, especially in the case of multi-threaded workload due to indeterminacy introduced due to simultaneously executing various threads. In this paper, we propose a technique for speeding multi-core simulation. The model of the processor core and cache are replaced with functional models, to achieve speedup. A timed Petri net model is used to estimate the execution time of the processor and the memory access latencies are estimated using hit/miss information obtained from the functional model of the cache. This model can be used to predict performance of data parallel applications or multiprogramming workload on CMP platform with various cache hierarchies and shared bus interconnect. The error in estimation of the execution time of an application is within 6%. The speedup achieved ranges between an average of 2x--4x over the cycle accurate simulator.

Indian Ocean dipole mode events in an ocean general circulation model

The evolution of the dipole mode (DM) events in the Indian Ocean is examined using an ocean model that is driven by the NCEP fluxes for the period 1975-1998. The positive DM events during 1997, 1994 and 1982 and negative DM events during 1996 and 1984-1985 are captured by the model and it reproduces both the surface and subsurface features associated with these events. In its positive phase, the DM is characterized by warmer than normal SST in the western Indian Ocean and cooler than normal SST in the eastern Indian Ocean. The DM events are accompanied by easterly wind anomalies along the equatorial Indian Ocean and upwelling-favorable alongshore wind anomalies along the coast of Sumatra. The Wyrtki jets are weak during positive DM events, and the thermocline is shallower than normal in the eastern Indian Ocean and deeper in the west. This anomaly pattern reverses during negative DM events. During the positive phase of the DM easterly wind anomalies excite an upwelling equatorial Kelvin wave. This Kelvin wave reflects from the eastern boundary as an upwelling Rossby wave which propagates westward across the equatorial Indian Ocean. The anomalies in the eastern Indian Ocean weaken after the Rossby wave passes. A similar process excites a downwelling Rossby wave during the negative phase. This Rossby wave is much weaker but wind forcing in the central equatorial Indian Ocean amplifies the downwelling and increases its westward phase speed. This Rossby wave initiates the deepening of the thermocline in the western Indian Ocean during the following positive phase of the DM. Rossby wave generated in the southern tropical Indian Ocean by Ekman pumping contributes to this warming. Concurrently, the temperature equation of the model shows upwelling and downwelling to be the most important mechanism during both positive events of 1994 and 1997. (C) 2002 Elsevier Science Ltd. All rights reserved.

Model exact low-lying states and spin dynamics in ferric wheels: Fe-6 to Fe-12

Using an efficient numerical scheme that exploits spatial symmetries and spin parity, we have obtained the exact low-lying eigenstates of exchange Hamiltonians for ferric wheels up to Fe-12. The largest calculation involves the Fe-12 ring which spans a Hilbert space dimension of about 145x10(6) for the M-S=0 subspace. Our calculated gaps from the singlet ground state to the excited triplet state agree well with the experimentally measured values. Study of the static structure factor shows that the ground state is spontaneously dimerized for ferric wheels. The spin states of ferric wheels can be viewed as quantized states of a rigid rotor with the gap between the ground and first excited states defining the inverse of the moment of inertia. We have studied the quantum dynamics of Fe-10 as a representative of ferric wheels. We use the low-lying states of Fe-10 to solve exactly the time-dependent Schrodinger equation and find the magnetization of the molecule in the presence of an alternating magnetic field at zero temperature. We observe a nontrivial oscillation of the magnetization which is dependent on the amplitude of the ac field. We have also studied the torque response of Fe-12 as a function of a magnetic field, which clearly shows spin-state crossover.

An extended Cahn-Hilliard model for interfaces with cubic anisotropy

For studying systems with a cubic anisotropy in interfacial energy sigma, we extend the Cahn-Hilliard model by including in it a fourth-rank term, namely, gamma (ijlm) [partial derivative (2) c/(partial derivativex(i) partial derivativex(j))] [partial derivative (2) c/(partial derivativex(l) partial derivativex(m))]. This term leads to an additional linear term in the evolution equation for the composition parameter field. It also leads to an orientation-dependent effective fourth-rank coefficient gamma ([hkl]) in the governing equation for the one-dimensional composition profile across a planar interface. The main effect of a non-negative gamma ([hkl]) is to increase both sigma and interfacial width w, each of which, upon suitable scaling, is related to gamma ([hkl]) through a universal scaling function. In this model, sigma is a differentiable function of interface orientation (n) over cap, and does not exhibit cusps; therefore, the equilibrium particle shapes (Wulff shapes) do not contain planar facets. However, the anisotropy in the interfacial energy can be large enough to give rise to corners in the Wulff shapes in two dimensions. In particles of finite sizes, the corners become rounded, and their shapes tend towards the Wulff shape with increasing particle size.

A nonlinear spectral finite element model for analysis of wave propagation in solid with internal friction and dissipation

A geometrically non-linear Spectral Finite Flement Model (SFEM) including hysteresis, internal friction and viscous dissipation in the material is developed and is used to study non-linear dissipative wave propagation in elementary rod under high amplitude pulse loading. The solution to non-linear dispersive dissipative equation constitutes one of the most difficult problems in contemporary mathematical physics. Although intensive research towards analytical developments are on, a general purpose cumputational discretization technique for complex applications, such as finite element, but with all the features of travelling wave (TW) solutions is not available. The present effort is aimed towards development of such computational framework. Fast Fourier Transform (FFT) is used for transformation between temporal and frequency domain. SFEM for the associated linear system is used as initial state for vector iteration. General purpose procedure involving matrix computation and frequency domain convolution operators are used and implemented in a finite element code. Convergnence of the spectral residual force vector ensures the solution accuracy. Important conclusions are drawn from the numerical simulations. Future course of developments are highlighted.

A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

Particle dynamics in a turbulent particle–gas suspension at high Stokes number. Part 2. The fluctuating-force model

A fluctuating-force model is developed for representing the effect of the turbulent fluid velocity fluctuations on the particle phase in a turbulent gas–solid suspension in the limit of high Stokes number, where the particle relaxation time is large compared with the correlation time for the fluid velocity fluctuations. In the model, a fluctuating force is incorporated in the equation of motion for the particles, and the force distribution is assumed to be an anisotropic Gaussian white noise. It is shown that this is equivalent to incorporating a diffusion term in the Boltzmann equation for the particle velocity distribution functions. The variance of the force distribution, or equivalently the diffusion coefficient in the Boltzmann equation, is related to the time correlation functions for the fluid velocity fluctuations. The fluctuating-force model is applied to the specific case of a Couette flow of a turbulent particle–gas suspension, for which both the fluid and particle velocity distributions were evaluated using direct numerical simulations by Goswami & Kumaran (2010). It is found that the fluctuating-force simulation is able to quantitatively predict the concentration, mean velocity profiles and the mean square velocities, both at relatively low volume fractions, where the viscous relaxation time is small compared with the time between collisions, and at higher volume fractions, where the time between collisions is small compared with the viscous relaxation time. The simulations are also able to predict the velocity distributions in the centre of the Couette, even in cases in which the velocity distribution is very different from a Gaussian distribution.