Parallelization of an unstructured data based cell centre finite volume code, HIFUN-3D

This work describes the parallelization of High Resolution flow solver on unstructured meshes, HIFUN-3D, an unstructured data based finite volume solver for 3-D Euler equations. For mesh partitioning, we use METIS, a software based on multilevel graph partitioning. The unstructured graph used for partitioning is associated with weights both on its vertices and edges. The data residing on every processor is split into four layers. Such a novel procedure of handling data helps in maintaining the effectiveness of the serial code. The communication of data across the processors is achieved by explicit message passing using the standard blocking mode feature of Message Passing Interface (MPI). The parallel code is tested on PACE++128 available in CFD Center

Influence of viscoelastic nature on the intermittent peel-front dynamics of adhesive tape

We investigate the influence of viscoelastic nature of the adhesive on the intermittent peel front dynamics by extending a recently introduced model for peeling of an adhesive tape. As time and rate-dependent deformation of the adhesives are measured in stationary conditions, a crucial step in incorporating the viscoelastic effects applicable to unstable intermittent peel dynamics is the introduction of a dynamization scheme that eliminates the explicit time dependence in terms of dynamical variables. We find contrasting influences of viscoelastic contribution in different regions of tape mass, roller inertia, and pull velocity. As the model acoustic energy dissipated depends on the nature of the peel front and its dynamical evolution, the combined effect of the roller inertia and pull velocity makes the acoustic energy noisier for small tape mass and low-pull velocity while it is burstlike for low-tape mass, intermediate values of the roller inertia and high-pull velocity. The changes are quantified by calculating the largest Lyapunov exponent and analyzing the statistical distributions of the amplitudes and durations of the model acoustic energy signals. Both single and two stage power-law distributions are observed. Scaling relations between the exponents are derived which show that the exponents corresponding to large values of event sizes and durations are completely determined by those for small values. Th scaling relations are found to be satisfied in all cases studied. Interestingly, we find only five types of model acoustic emission signals among multitude of possibilities of the peel front configurations.

Chronopotentiometry with power-law perturbation functions at an expanding plane electrode with and without a preceding blank period for systems with a coupled first-order homogeneous chemical reaction

The systems formalism is used to obtain the interfacial concentration transients for power-law current input at an expanding plane electrode. The explicit results for the concentration transients obtained here pertain to arbitrary homogeneous reaction schemes coupled to the oxidant and reductant of a single charge-transfer step and the power-law form without and with a preceding blank period (for two types of power-law current profile, say, (i) I(t) = I0(t−t0)q for t greater-or-equal, slanted t0, I(t) = 0 for t < t0; and (ii) I(t) = I0tq for t greater-or-equal, slanted t0, I(t) = 0 for t < t0). Finally the potential transients are obtained using Padé approximants. The results of Galvez et al. (for E, CE, EC, aC) (J. Electroanal. Chem., 132 (1982) 15; 146 (1983) 221, 233, 243), Molina et al. (for E) (J. Electroanal. Chem., 227 (1987) 1 and Kies (for E) (J. Electroanal. Chem., 45 (1973) 71) are obtained as special cases.

Quantum Error Correction via Codes over GF(2)

It is well known that n-length stabilizer quantum error correcting codes (QECCs) can be obtained via n-length classical error correction codes (CECCs) over GF(4), that are additive and self-orthogonal with respect to the trace Hermitian inner product. But, most of the CECCs have been studied with respect to the Euclidean inner product. In this paper, it is shown that n-length stabilizer QECCs can be constructed via 371 length linear CECCs over GF(2) that are self-orthogonal with respect to the Euclidean inner product. This facilitates usage of the widely studied self-orthogonal CECCs to construct stabilizer QECCs. Moreover, classical, binary, self-orthogonal cyclic codes have been used to obtain stabilizer QECCs with guaranteed quantum error correcting capability. This is facilitated by the fact that (i) self-orthogonal, binary cyclic codes are easily identified using transform approach and (ii) for such codes lower bounds on the minimum Hamming distance are known. Several explicit codes are constructed including two pure MDS QECCs.

Multiple contour extraction from graylevel images using an artificial neural network

For active contour modeling (ACM), we propose a novel self-organizing map (SOM)-based approach, called the batch-SOM (BSOM), that attempts to integrate the advantages of SOM- and snake-based ACMs in order to extract the desired contours from images. We employ feature points, in the form of ail edge-map (as obtained from a standard edge-detection operation), to guide the contour (as in the case of SOM-based ACMs) along with the gradient and intensity variations in a local region to ensure that the contour does not "leak" into the object boundary in case of faulty feature points (weak or broken edges). In contrast with the snake-based ACMs, however, we do not use an explicit energy functional (based on gradient or intensity) for controlling the contour movement. We extend the BSOM to handle extraction of contours of multiple objects, by splitting a single contour into as many subcontours as the objects in the image. The BSOM and its extended version are tested on synthetic binary and gray-level images with both single and multiple objects. We also demonstrate the efficacy of the BSOM on images of objects having both convex and nonconvex boundaries. The results demonstrate the superiority of the BSOM over others. Finally, we analyze the limitations of the BSOM.

STBCs with Reduced Sphere Decoding Complexity for Two-User MIMO-MAC

In this paper, Space-Time Block Codes (STBCs) with reduced Sphere Decoding Complexity (SDC) are constructed for two-user Multiple-Input Multiple-Output (MIMO) fading multiple access channels. In this set-up, both the users employ identical STBCs and the destination performs sphere decoding for the symbols of the two users. First, we identify the positions of the zeros in the R matrix arising out of the Q-R decomposition of the lattice generator such that (i) the worst case SDC (WSDC) and (ii) the average SDC (ASDC) are reduced. Then, a set of necessary and sufficient conditions on the lattice generator is provided such that the R matrix has zeros at the identified positions. Subsequently, explicit constructions of STBCs which results in the reduced ASDC are presented. The rate (in complex symbols per channel use) of the proposed designs is at most 2/N-t where N-t denotes the number of transmit antennas for each user. We also show that the class of STBCs from complex orthogonal designs (other than the Alamouti design) reduce the WSDC but not the ASDC.

A numerical study of laminar axisymmetric plumes due to the combined buoyancy of heat and mass diffusion

This paper presents the results of a computational study of laminar axisymmetric plumes generated by the simultaneous diffusion of thermal energy and chemical species. Species concentrations are assumed small. The plume is treated as a boundary layer. Boussinesq approximations are incorporated and the governing conservation equations of mass, momentum, energy and species are suitably non-dimensionalised. These equations are solved using one time-step-forward explicit finite-difference method. Upwind differencing is employed for convective terms. The results thus obtained are explained in terms of the basic physical mechanisms that govern these flows. They show many interesting aspects of the complex interaction of the two buoyant mechanisms.

Optimal Weighted Antenna Selection For Imperfect Channel Knowledge From Training

Receive antenna selection (AS) reduces the hardware complexity of multi-antenna receivers by dynamically connecting an instantaneously best antenna element to the available radio frequency (RF) chain. Due to the hardware constraints, the channels at various antenna elements have to be sounded sequentially to obtain estimates that are required for selecting the ``best'' antenna and for coherently demodulating data. Consequently, the channel state information at different antennas is outdated by different amounts. We show that, for this reason, simply selecting the antenna with the highest estimated channel gain is not optimum. Rather, the channel estimates of different antennas should be weighted differently, depending on the training scheme. We derive closed-form expressions for the symbol error probability (SEP) of AS for MPSK and MQAM in time-varying Rayleigh fading channels for arbitrary selection weights, and validate them with simulations. We then derive an explicit formula for the optimal selection weights that minimize the SEP. We find that when selection weights are not used, the SEP need not improve as the number of antenna elements increases, which is in contrast to the ideal channel estimation case. However, the optimal selection weights remedy this situation and significantly improve performance.

On the admittance of lipid bilayer membranes: Part I. Membrane-permeable ions

The modular formalism of Rangarajan [J. Electroanal. Chem., 55 (1974) 297] has been applied to the admittance of lipid bilayer membranes. The method leads to equations which clearly show the interrelations between the various partial processes involved in ion transport, and which allow examination of model assumptions without the need for a complete rederivation of the membrane admittance. Explicit expressions are given for both the continuum and single jump models. The former includes the ionic displacement component, important mostly at high frequencies.

Completeness of the Coulomb wave functions in quantum mechanics

We give an explicit, direct, and fairly elementary proof that the radial energy eigenfunctions for the hydrogen atom in quantum mechanics, bound and scattering states included, form a complete set. The proof uses only some properties of the confluent hypergeometric functions and the Cauchy residue theorem from analytic function theory; therefore it would form useful supplementary reading for a graduate course on quantum mechanics.

On the admittance of lipid bilayer membranes: Part I. Membrane-permeable ions

The modular formalism of Rangarajan [J. Electroanal. Chem., 55 (1974) 297] has been applied to the admittance of lipid bilayer membranes. The method leads to equations which clearly show the interrelations between the various partial processes involved in ion transport, and which allow examination of model assumptions without the need for a complete rederivation of the membrane admittance. Explicit expressions are given for both the continuum and single jump models. The former includes the ionic displacement component, important mostly at high frequencies.

Investigation of the effect of nonlocal scale on ultrasonic wave dispersion characteristics of a monolayer graphene

In this article, an ultrasonic wave propagation in graphene sheet is studied using nonlocal elasticity theory incorporating small scale effects. The graphene sheet is modeled as an isotropic plate of one-atom thick. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of wave propagation model is also derived for the graphene sheet. The nonlocal scale parameter introduces certain band gap region in in-plane and flexural wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. The explicit expressions for cutoff frequencies and escape frequencies are derived. The escape frequencies are mainly introduced because of the nonlocal elasticity. Obviously these frequencies are function of nonlocal scaling parameter. It has also been obtained that these frequencies are independent of y-directional wavenumber. It means that for any type of nanostructure, the escape frequencies are purely a function of nonlocal scaling parameter only. It is also independent of the geometry of the structure. It has been found that the cutoff frequencies are function of nonlocal scaling parameter (e(0)a) and the y-directional wavenumber (k(y)). For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(o)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

Sequential design of decentralized control systems

Our main result is a new sequential method for the design of decentralized control systems. Controller synthesis is conducted on a loop-by-loop basis, and at each step the designer obtains an explicit characterization of the class C of all compensators for the loop being closed that results in closed-loop system poles being in a specified closed region D of the s-plane, instead of merely stabilizing the closed-loop system. Since one of the primary goals of control system design is to satisfy basic performance requirements that are often directly related to closed-loop pole location (bandwidth, percentage overshoot, rise time, settling time), this approach immediately allows the designer to focus on other concerns such as robustness and sensitivity. By considering only compensators from class C and seeking the optimum member of that set with respect to sensitivity or robustness, the designer has a clearly-defined limited optimization problem to solve without concern for loss of performance. A solution to the decentralized tracking problem is also provided. This design approach has the attractive features of expandability, the use of only 'local models' for controller synthesis, and fault tolerance with respect to certain types of failure.

Time Dependent Operator-split and Unsplit Schemes for One Dimensional Premixed Flames

This paper is concerned with a study of an operator split scheme and unsplit scheme for the computation of adiabatic freely propagating one-dimensional premixed flames. The study uses unsteady method for both split and unsplit schemes employing implicit chemistry and explicit diffusion, a combination which is stable and convergent. Solution scheme is not sensitive to the initial starting estimate and provides steady state even with straight line profiles (far from steady state) in small number of time steps. Two systems H2-Air and H2-NO (involving complex nitrogen chemistry) are considered in presentinvestigation. Careful comparison shows that the operator split approach is slightly superior than the unsplit when chemistry becomes complex. Comparison of computational times with those of existing steady and unsteady methods seems to suggest that the method employing implicit-explicit algorithm is very efficient and robust.

Volume effects and associations in liquid alloys

The regular associated solution model for binary systems has been modified by incorporating the size of the complex as an explicit variable. The thermodynamic properties of the liquid alloy and the interactions between theA ?B type of complex and the unassociated atoms in anA-B binary have been evaluated as a function of relative size of the complex using the activity coefficients at infinite dilution and activity data at one other composition in the binary. The computational procedure adopted for determining the concentration of clusters and interaction energies in the associated liquid is similar to that proposed by Lele and Rao. The analysis has been applied to the thermodynamic mixing functions of liquid Al-Ca alloys believed to contain Al2Ca associates. It is found that the size of the cluster significantly affects the interaction energies between the complex and the unassociated atoms, while the equilibrium constant and enthalpy change for the association reaction exhibit only minor variation, when the equations are fitted to experimental data. The interaction energy between unassociated free atoms remains virtually unaltered as the size of the complex is varied between extreme values. Accurate data on free energy, enthalpy, and volume of mixing at the same temperature on alloy systems with compound forming tendency would permit a rigorous test of the proposed model.