Two-dimensional moist stratified turbulence and the emergence of vertically sheared horizontal flows

Moist stratified turbulence is studied in a two-dimensional Boussinesq system influenced by condensation and evaporation. The problem is set in a periodic domain and employs simple evaporation and condensation schemes, wherein both the processes push parcels towards saturation. Numerical simulations demonstrate the emergence of a moist turbulent state consisting of ordered structures with a clear power-law type spectral scaling from initially spatially uncorrelated conditions. An asymptotic analysis in the limit of rapid condensation and strong stratification shows that, for initial conditions with enough water substance to saturate the domain, the equations support a straightforward state of moist balance characterized by a hydrostatic, saturated, vertically sheared horizontal flow (VSHF). For such initial conditions, by means of long time numerical simulations, the emergence of moist balance is verified. Specifically, starting from uncorrelated data, subsequent to the development of a moist turbulent state, the system experiences a rather abrupt transition to a regime which is close to saturation and dominated by a strong VSHF. On the other hand, initial conditions which do not have enough water substance to saturate the domain, do not attain moist balance. Rather, the system is observed to remain in a turbulent state and oscillates about moist balance. Even though balance is not achieved with these general initial conditions, the time scale of oscillation about moist balance is much larger than the imposed time scale of condensation and evaporation, thus indicating a distinct dominant slow component in the moist stratified two-dimensional turbulent system. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3694805]

Fatigue crack growth due to overloads in plain concrete using scaling laws

Scaling laws are represented in power law form and can be utilized to extract the characteristic properties of a new phenomenon with the help of self-similar solutions. In this work, an attempt has been made to propose a scaling law analytically, for plain concrete when subjected to variable amplitude loading. Due to the application of overload on concrete structures, acceleration in the crack growth process takes place. A closed form expression has been developed to capture the acceleration in crack growth rate in conjunction with the principles of dimensional analysis and self-similarity. The proposed model accounts for parameters such as, the tensile strength, fracture toughness, overload effect and the structural size. Knowing the governed and the governing parameters of the physical problem and by using the concepts of self-similarity, a relationship is obtained between the different parameters involved. The predicted results are compared with experimental crack growth data for variable amplitude loading and are found to capture the overload effect with sufficient accuracy. Through a sensitivity analysis, fracture toughness is found to be the most dominant parameter in accelerating the crack length due to application of overload.

Practical fully three-dimensional reconstruction algorithms for diffuse optical tomography

We have developed an efficient fully three-dimensional (3D) reconstruction algorithm for diffuse optical tomography (DOT). The 3D DOT, a severely ill-posed problem, is tackled through a pseudodynamic (PD) approach wherein an ordinary differential equation representing the evolution of the solution on pseudotime is integrated that bypasses an explicit inversion of the associated, ill-conditioned system matrix. One of the most computationally expensive parts of the iterative DOT algorithm, the reevaluation of the Jacobian in each of the iterations, is avoided by using the adjoint-Broyden update formula to provide low rank updates to the Jacobian. In addition, wherever feasible, we have also made the algorithm efficient by integrating along the quadratic path provided by the perturbation equation containing the Hessian. These algorithms are then proven by reconstruction, using simulated and experimental data and verifying the PD results with those from the popular Gauss-Newton scheme. The major findings of this work are as follows: (i) the PD reconstructions are comparatively artifact free, providing superior absorption coefficient maps in terms of quantitative accuracy and contrast recovery; (ii) the scaling of computation time with the dimension of the measurement set is much less steep with the Jacobian update formula in place than without it; and (iii) an increase in the data dimension, even though it renders the reconstruction problem less ill conditioned and thus provides relatively artifact-free reconstructions, does not necessarily provide better contrast property recovery. For the latter, one should also take care to uniformly distribute the measurement points, avoiding regions close to the source so that the relative strength of the derivatives for measurements away from the source does not become insignificant. (c) 2012 Optical Society of America

Projecting low-dimensional chaos from spatiotemporal dynamics in a model for plastic instability

We investigate the possibility of projecting low-dimensional chaos from spatiotemporal dynamics of a model for a kind of plastic instability observed under constant strain rate deformation conditions. We first discuss the relationship between the spatiotemporal patterns of the model reflected in the nature of dislocation bands and the nature of stress serrations. We show that at low applied strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of mobile dislocation density and the stress drops. We then show that the model equations are spatiotemporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov dimension scale with the system size at low and high strain rates. Using a modified algorithm for calculating correlation dimension density, we show that the stress-strain signals at low applied strain rates corresponding to spatially uncorrelated dislocation bands exhibit features of low-dimensional chaos. This is made quantitative by demonstrating that the model equations can be approximately reduced to space-independent model equations for the average dislocation densities, which is known to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks with increasing applied strain rate due to increasing propensity for propagation of the dislocation bands.

Fidelity susceptibility of one-dimensional models with twisted boundary conditions

Recently it has been shown that the fidelity of the ground state of a quantum many-body system can be used todetect its quantum critical points (QCPs). If g denotes the parameter in the Hamiltonian with respect to which the fidelity is computed, we find that for one-dimensional models with large but finite size, the fidelity susceptibility chi(F) can detect a QCP provided that the correlation length exponent satisfies nu < 2. We then show that chi(F) can be used to locate a QCP even if nu >= 2 if we introduce boundary conditions labeled by a twist angle N theta, where N is the system size. If the QCP lies at g = 0, we find that if N is kept constant, chi(F) has a scaling form given by chi(F) similar to theta(-2/nu) f (g/theta(1/nu)) if theta << 2 pi/N. We illustrate this both in a tight-binding model of fermions with a spatially varying chemical potential with amplitude h and period 2q in which nu = q, and in a XY spin-1/2 chain in which nu = 2. Finally we show that when q is very large, the model has two additional QCPs at h = +/- 2 which cannot be detected by studying the energy spectrum but are clearly detected by chi(F). The peak value and width of chi(F) seem to scale as nontrivial powers of q at these QCPs. We argue that these QCPs mark a transition between extended and localized states at the Fermi energy. DOI: 10.1103/PhysRevB.86.245424

Generalized bidirectional multi-pair multi-antennawireless network coding

In this paper, we focus on increasing the throughput and diversity of network coded MIMO transmissions in bidirectional multi-pair wireless relay networks. All nodes have multi-antenna capability. Pairs of nodes want to exchange messages via a relay having multi-antenna and encoding/decoding capability. Nodes transmit their messages to the relay in the first (MAC) phase. The relay decodes all the messages and XORs them and broadcasts the XORed message in the second (BC) phase. We develop a generalized framework for bidirectional multi-pair multi-antenna wireless network coding, which models different MIMO transmission schemes including spatial multiplexing (V-BLAST), orthogonal STBC (OSTBC), and non-orthogonal STBC (NO-STBC) in a unified way. Enhanced throughputs are achieved by allowing all nodes to simultaneously transmit at their full rate. High diversity orders are achieved through the use of NO-STBCs, characterized by full rate and full transmit diversity. We evaluate and compare the performance of VBLAST, OSTBC, and NO-STBC schemes in one-dimensional 1-pair linear network (one pair of nodes and a relay) and two-dimensional 2-pair `cross' network (two pairs of nodes and a relay).

Simulation of length-preserving motions of flexible one dimensional oObjects using optimization

During the motion of one dimensional flexible objects such as ropes, chains, etc., the assumption of constant length is realistic. Moreover,their motion appears to be naturally minimizing some abstract distance measure, wherein the disturbance at one end gradually dies down along the curve defining the object. This paper presents purely kinematic strategies for deriving length-preserving transformations of flexible objects that minimize appropriate ‘motion’. The strategies involve sequential and overall optimization of the motion derived using variational calculus. Numerical simulations are performed for the motion of a planar curve and results show stable converging behavior for single-step infinitesimal and finite perturbations 1 as well as multi-step perturbations. Additionally, our generalized approach provides different intuitive motions for various problem-specific measures of motion, one of which is shown to converge to the conventional tractrix-based solution. Simulation results for arbitrary shapes and excitations are also included.

Projecting low and extensive dimensional chaos from spatio-temporal dynamics

We review the spatio-temporal dynamical features of the Ananthakrishna model for the Portevin-Le Chatelier effect, a kind of plastic instability observed under constant strain rate deformation conditions. We then establish a qualitative correspondence between the spatio-temporal structures that evolve continuously in the instability domain and the nature of the irregularity of the scalar stress signal. Rest of the study is on quantifying the dynamical information contained in the stress signals about the spatio-temporal dynamics of the model. We show that at low applied strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of mobile dislocation density and the stress drops. We then show that the model equations are spatio-temporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov dimension scale with the system size at low and high strain rates. Using a modified algorithm for calculating correlation dimension density, we show that the stress-strain signals at low applied strain rates corresponding to spatially uncorrelated dislocation bands exhibit features of low dimensional chaos. This is made quantitative by demonstrating that the model equations can be approximately reduced to space independent model equations for the average dislocation densities, which is known to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks with increasing applied strain rate due to increasing propensity for propagation of the dislocation bands. The stress signals in the partially propagating to fully propagating bands turn to have features of extensive chaos.

Joint Self-Iterating Equalization and Detection for Two-Dimensional Intersymbol-Interference Channels

We develop several novel signal detection algorithms for two-dimensional intersymbol-interference channels. The contribution of the paper is two-fold: (1) We extend the one-dimensional maximum a-posteriori (MAP) detection algorithm to operate over multiple rows and columns in an iterative manner. We study the performance vs. complexity trade-offs for various algorithmic options ranging from single row/column non-iterative detection to a multi-row/column iterative scheme and analyze the performance of the algorithm. (2) We develop a self-iterating 2-D linear minimum mean-squared based equalizer by extending the 1-D linear equalizer framework, and present an analysis of the algorithm. The iterative multi-row/column detector and the self-iterating equalizer are further connected together within a turbo framework. We analyze the combined 2-D iterative equalization and detection engine through analysis and simulations. The performance of the overall equalizer and detector is near MAP estimate with tractable complexity, and beats the Marrow Wolf detector by about at least 0.8 dB over certain 2-D ISI channels. The coded performance indicates about 8 dB of significant SNR gain over the uncoded 2-D equalizer-detector system.

Turbulence in the two-dimensional Fourier-truncated Gross-Pitaevskii equation

We undertake a systematic, direct numerical simulation of the twodimensional, Fourier-truncated, Gross-Pitaevskii equation to study the turbulent evolutions of its solutions for a variety of initial conditions and a wide range of parameters. We find that the time evolution of this system can be classified into four regimes with qualitatively different statistical properties. Firstly, there are transients that depend on the initial conditions. In the second regime, powerlaw scaling regions, in the energy and the occupation-number spectra, appear and start to develop; the exponents of these power laws and the extents of the scaling regions change with time and depend on the initial condition. In the third regime, the spectra drop rapidly for modes with wave numbers k > kc and partial thermalization takes place for modes with k < kc; the self-truncation wave number kc(t) depends on the initial conditions and it grows either as a power of t or as log t. Finally, in the fourth regime, complete thermalization is achieved and, if we account for finite-size effects carefully, correlation functions and spectra are consistent with their nontrivial Berezinskii-Kosterlitz-Thouless forms. Our work is a natural generalization of recent studies of thermalization in the Euler and other hydrodynamical equations; it combines ideas from fluid dynamics and turbulence, on the one hand, and equilibrium and nonequilibrium statistical mechanics on the other.

On the k(1)(-1) scaling in sink-flow turbulent boundary layers

Scaling of the streamwise velocity spectrum phi(11)(k(1)) in the so-called sink-flow turbulent boundary layer is investigated in this work. The present experiments show strong evidence for the k(1)(-1) scaling i.e. phi(11)(k(1)) = Lambda(1)U(tau)(2)k(1)(-1), where k(1)(-1) is the streamwise wavenumber and U-tau is the friction velocity. Interestingly, this k(1)(-1) scaling is observed much farther from the wall and at much lower flow Reynolds number (both differing by almost an order of magnitude) than what the expectations from experiments on a zero-pressure-gradient turbulent boundary layer flow would suggest. Furthermore, the coefficient A(1) in the present sink-flow data is seen to be non-universal, i.e. A(1) varies with height from the wall; the scaling exponent -1 remains universal. Logarithmic variation of the so-called longitudinal structure function, which is the physical-space counterpart of spectral k(1)(-1) scaling, is also seen to be non-universal, consistent with the non-universality of A(1). These observations are to be contrasted with the universal value of A(1) (along with the universal scaling exponent of 1) reported in the literature on zero-pressure-gradient turbulent boundary layers. Theoretical arguments based on dimensional analysis indicate that the presence of a streamwise pressure gradient in sink-flow turbulent boundary layers makes the coefficient A(1) non-universal while leaving the scaling exponent -1 unaffected. This effect of the pressure gradient on the streamwise spectra, as discussed in the present study (experiments as well as theory), is consistent with other recent studies in the literature that are focused on the structural aspects of turbulent boundary layer flows in pressure gradients (Harun etal., J. Flui(d) Mech., vol. 715, 2013, pp. 477-498); the present paper establishes the link between these two. The variability of A(1) accommodated in the present framework serves to clarify the ideas of universality of the k(1)(-1) scaling.

Characterization, charge transport and magnetic properties of multi-walled carbon nanotube-polyvinyl chloride nanocomposites

Multi-walled carbon nanotube (MWCNT)-polyvinyl chloride (PVC) nanocomposites, with MWCNT loading up to 44.4 weight percent (wt%), were prepared by the solvent mixing and casting method. Electron microscopy indicates high degree of dispersion of MWCNT in PVC matrix, achieved by ultrasonication without using any surfactants. Thermogravimetric analysis showed a significant monotonic enhancement in the thermal stability of nanocomposites by increasing the wt% of MWCNT. Electrical conductivity of nanocomposites followed the classical percolation theory and the conductivity prominently improved from 10(-7) to 9 S/cm as the MWCNT loading increased from 0.1 to 44.4 wt%. Low value of electrical percolation threshold similar to 0.2 wt% is achieved which is attributed to high aspect ratio and homogeneous dispersion of MWCNT in PVC. The analysis of the low temperature electrical resistivity data shows that sample of 1.9 wt% follows three dimensional variable range hopping model whereas higher wt% nanocomposite samples follow power law behavior. The magnetization versus applied field data for both bulk MWCNTs and nanocomposite of 44.4 wt% display ferromagnetic behavior with enhanced coercivities of 1.82 and 1.27 kOe at 10 K, respectively. The enhancement in coercivity is due to strong dipolar interaction and shape anisotropy of rod-shaped iron nanoparticles. (C) 2013 Elsevier B.V. All rights reserved.

Electrostatic modulation of periodic potentials in a two-dimensional electron gas: from antidot lattice to quantum dot lattice

We use a dual gated device structure to introduce a gate-tuneable periodic potential in a GaAs/AlGaAs two dimensional electron gas (2DEG). Using only a suitable choice of gate voltages we can controllably alter the potential landscape of the bare 2DEG, inducing either a periodic array of antidots or quantum dots. Antidots are artificial scattering centers, and therefore allow for a study of electron dynamics. In particular, we show that the thermovoltage of an antidot lattice is particularly sensitive to the relative positions of the Fermi level and the antidot potential. A quantum dot lattice, on the other hand, provides the opportunity to study correlated electron physics. We find that its current-voltage characteristics display a voltage threshold, as well as a power law scaling, indicative of collective Coulomb blockade in a disordered background.

Automatic Data Allocation and Buffer Management for Multi-GPU Machines

Multi-GPU machines are being increasingly used in high-performance computing. Each GPU in such a machine has its own memory and does not share the address space either with the host CPU or other GPUs. Hence, applications utilizing multiple GPUs have to manually allocate and manage data on each GPU. Existing works that propose to automate data allocations for GPUs have limitations and inefficiencies in terms of allocation sizes, exploiting reuse, transfer costs, and scalability. We propose a scalable and fully automatic data allocation and buffer management scheme for affine loop nests on multi-GPU machines. We call it the Bounding-Box-based Memory Manager (BBMM). BBMM can perform at runtime, during standard set operations like union, intersection, and difference, finding subset and superset relations on hyperrectangular regions of array data (bounding boxes). It uses these operations along with some compiler assistance to identify, allocate, and manage data required by applications in terms of disjoint bounding boxes. This allows it to (1) allocate exactly or nearly as much data as is required by computations running on each GPU, (2) efficiently track buffer allocations and hence maximize data reuse across tiles and minimize data transfer overhead, and (3) and as a result, maximize utilization of the combined memory on multi-GPU machines. BBMM can work with any choice of parallelizing transformations, computation placement, and scheduling schemes, whether static or dynamic. Experiments run on a four-GPU machine with various scientific programs showed that BBMM reduces data allocations on each GPU by up to 75% compared to current allocation schemes, yields performance of at least 88% of manually written code, and allows excellent weak scaling.

Selection and thermodynamic analysis of a turbocharger for a producer gas-fuelled multi-cylinder engine

The current work addresses the use of producer gas, a bio-derived gaseous alternative fuel, in engines designed for natural gas, derived from diesel engine frames. Impact of the use of producer gas on the general engine performance with specific focus on turbo-charging is addressed. The operation of a particular engine frame with diesel, natural gas and producer gas indicates that the peak load achieved is highest with diesel fuel (in compression ignition mode) followed by natural gas and producer gas (both in spark ignite mode). Detailed analysis of the engine power de-rating on fuelling with natural gas and producer gas indicates that the change in compression ratio (migration from compression to spark ignited mode), difference in mixture calorific value and turbocharger mismatch are the primary contributing factors. The largest de-rating occurs due to turbocharger mismatch. Turbocharger selection and optimization is identified as the strategy to recover the non-thermodynamic power loss, identified as the recovery potential (the loss due to mixture calorific value and turbocharger mismatch) on operating the engine with a fuel different from the base fuel. A turbocharged after-cooled six cylinder, 5.9 l, 90 kWe (diesel rating) engine (12.2 bar BMEP) is available commercially as a naturally aspirated natural gas engine delivering a peak load of 44.0 kWe (6.0 bar BMEP). The engine delivers a load of 27.3 kWe with producer gas under naturally aspirated mode. On charge boosting the engine with a turbocharger similar in configuration to the diesel engine turbocharger, the peak load delivered with producer gas is 36 kWe (4.8 bar BMEP) indicating a de-rating of about 60% over the baseline diesel mode. Estimation of knock limited peak load for producer gas-fuelled operation on the engine frame using a Wiebe function-based zero-dimensional code indicates a knock limited peak load of 76 kWe, indicating the potential to recover about 40 kWe. As a part of the recovery strategy, optimizing the ignition timing for maximum brake torque based on both spark sweep tests and established combustion descriptors and engine-turbocharger matching for producer gas-fuelled operation resulted in a knock limited peak load of 72.8 kWe (9.9 bar BMEP) at a compressor pressure ratio of 2.30. The de-rating of about 17.0 kWe compared to diesel rating is attributed to the reduction in compression ratio. With load recovery, the specific biomass consumption reduces from 1.2 kg/kWh to 1.0 kg/kWh, an improvement of over 16% while the engine thermal efficiency increases from 28% to 32%. The thermodynamic analysis of the compressor and the turbine indicates an isentropic efficiency of 74.5% and 73%, respectively.