THREE-LAYER FLUID FLOW OVER A SMALL OBSTRUCTION ON THE BOTTOM OF A CHANNEL

Many boundary value problems occur in a natural way while studying fluid flow problems in a channel. The solutions of two such boundary value problems are obtained and analysed in the context of flow problems involving three layers of fluids of different constant densities in a channel, associated with an impermeable bottom that has a small undulation. The top surface of the channel is either bounded by a rigid lid or free to the atmosphere. The fluid in each layer is assumed to be inviscid and incompressible, and the flow is irrotational and two-dimensional. Only waves that are stationary with respect to the bottom profile are considered in this paper. The effect of surface tension is neglected. In the process of obtaining solutions for both the problems, regular perturbation analysis along with a Fourier transform technique is employed to derive the first-order corrections of some important physical quantities. Two types of bottom topography, such as concave and convex, are considered to derive the profiles of the interfaces. We observe that the profiles are oscillatory in nature, representing waves of variable amplitude with distinct wave numbers propagating downstream and with no wave upstream. The observations are presented in tabular and graphical forms.

CHEXVIS: a tool for molecular channel extraction and visualization

Background: Understanding channel structures that lead to active sites or traverse the molecule is important in the study of molecular functions such as ion, ligand, and small molecule transport. Efficient methods for extracting, storing, and analyzing protein channels are required to support such studies. Further, there is a need for an integrated framework that supports computation of the channels, interactive exploration of their structure, and detailed visual analysis of their properties. Results: We describe a method for molecular channel extraction based on the alpha complex representation. The method computes geometrically feasible channels, stores both the volume occupied by the channel and its centerline in a unified representation, and reports significant channels. The representation also supports efficient computation of channel profiles that help understand channel properties. We describe methods for effective visualization of the channels and their profiles. These methods and the visual analysis framework are implemented in a software tool, CHEXVIS. We apply the method on a number of known channel containing proteins to extract pore features. Results from these experiments on several proteins show that CHEXVIS performance is comparable to, and in some cases, better than existing channel extraction techniques. Using several case studies, we demonstrate how CHEXVIS can be used to study channels, extract their properties and gain insights into molecular function. Conclusion: CHEXVIS supports the visual exploration of multiple channels together with their geometric and physico-chemical properties thereby enabling the understanding of the basic biology of transport through protein channels. The CHEXVIS web-server is freely available at http://vgl.serc.iisc.ernet.in/chexvis/. The web-server is supported on all modern browsers with latest Java plug-in.

Behaviour of turbulent Prandtl/Schmidt number in compressible mixing layer

The behaviour of turbulent Prandtl/Schmidt number is explored through the model-free simulation results. It has been observed that compressibility affects the Reynolds scalar flux vectors. Reduced peak values are also observed for compressible convective Mach number mixing layer as compared with the incompressible convective Mach number counterpart, indicating a reduction in the mixing of enthalpy and species. Pr-t and Sc-t variations also indicate a reduction in mixing. It is observed that unlike the incompressible case, it is difficult to assign a constant value to these numbers due to their continuous variation in space. Modelling of Pr-t and Sc-t would be necessary to cater for this continuous spatial variation. However, the turbulent Lewis number is evaluated to be near unity for the compressible case, making it necessary to model only one of the Pr-t and Sc-t..

Role of deformation temperature on the evolution and heterogeneity of texture during equal channel angular pressing of magnesium

Investigations on texture evolution and through-thickness texture heterogeneity during equal channel angular pressing (ECAP) of pure magnesium at 200 degrees C, 150 degrees C and room temperature (RT) was carried out by neutron, high energy synchrotron X-ray and electron back-scatter diffraction. Irrespective of the ECAP temperature, a distinctive basal (B) and pyramidal (C-2) II type of fibers forms. The texture differs in the bottom 1 mm portion, where the B-fiber is shifted similar to 55 degrees due to negative shear attributed to friction. (C) 2015 Elsevier Inc. All rights reserved.

Outer Bounds on the Secrecy Rate of the 2-User Symmetric Deterministic Interference Channel with Transmitter Cooperation

This paper derives outer bounds for the 2-user symmetric linear deterministic interference channel (SLDIC) with limited-rate transmitter cooperation and perfect secrecy constraints at the receivers. Five outer bounds are derived, under different assumptions of providing side information to receivers and partitioning the encoded message/output depending on the relative strength of the signal and the interference. The usefulness of these outer bounds is shown by comparing the bounds with the inner bound on the achievable secrecy rate derived by the authors in a previous work. Also, the outer bounds help to establish that sharing random bits through the cooperative link can achieve the optimal rate in the very high interference regime.

Relay Selection with Channel Probing in Sleep-Wake Cycling Wireless Sensor Networks

In geographical forwarding of packets in a large wireless sensor network (WSN) with sleep-wake cycling nodes, we are interested in the local decision problem faced by a node that has ``custody'' of a packet and has to choose one among a set of next-hop relay nodes to forward the packet toward the sink. Each relay is associated with a ``reward'' that summarizes the benefit of forwarding the packet through that relay. We seek a solution to this local problem, the idea being that such a solution, if adopted by every node, could provide a reasonable heuristic for the end-to-end forwarding problem. Toward this end, we propose a local relay selection problem consisting of a forwarding node and a collection of relay nodes, with the relays waking up sequentially at random times. At each relay wake-up instant, the forwarder can choose to probe a relay to learn its reward value, based on which the forwarder can then decide whether to stop (and forward its packet to the chosen relay) or to continue to wait for further relays to wake up. The forwarder's objective is to select a relay so as to minimize a combination of waiting delay, reward, and probing cost. The local decision problem can be considered as a variant of the asset selling problem studied in the operations research literature. We formulate the local problem as a Markov decision process (MDP) and characterize the solution in terms of stopping sets and probing sets. We provide results illustrating the structure of the stopping sets, namely, the (lower bound) threshold and the stage independence properties. Regarding the probing sets, we make an interesting conjecture that these sets are characterized by upper bounds. Through simulation experiments, we provide valuable insights into the performance of the optimal local forwarding and its use as an end-to-end forwarding heuristic.

NUMERICAL STUDIES OF DYNAMO ACTION IN A TURBULENT SHEAR FLOW. I.

We perform numerical experiments to study the shear dynamo problem where we look for the growth of a large-scale magnetic field due to non-helical stirring at small scales in a background linear shear flow in previously unexplored parameter regimes. We demonstrate the large-scale dynamo action in the limit where the fluid Reynolds number (Re) is below unity while the magnetic Reynolds number (Rm) is above unity; the exponential growth rate scales linearly with shear, which is consistent with earlier numerical works. The limit of low Re is particularly interesting, as seeing the dynamo action in this limit would provide enough motivation for further theoretical investigations, which may focus attention on this analytically more tractable limit of Re < 1 compared to the more formidable limit of Re > 1. We also perform simulations in the regimes where (i) both (Re, Rm) < 1, and (ii) Re > 1 and Rm < 1, and compute all of the components of the turbulent transport coefficients (alpha(ij) and alpha(ij)) using the test-field method. A reasonably good agreement is observed between our results and the results of earlier analytical works in similar parameter regimes.

Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics

The interaction of a single bubble with a single vortex ring in water has been studied experimentally. Measurements of both the bubble dynamics and vorticity dynamics have been done to help understand the two-way coupled problem. The circulation strength of the vortex ring (Gamma) has been systematically varied, while keeping the bubble diameter (D-b) constant, with the bubble volume to vortex core volume ratio (V-R) also kept fixed at about 0.1. The other important parameter in the problem is a Weber number based on the vortex ring strength. (We = 0.87 rho(Gamma/2 pi a)(2)/(sigma/D-b); a = vortex core radius, sigma = surface tension), which is varied over a large range, We = 3-406. The interaction between the bubble and ring for each of the We cases broadly falls into four stages. Stage I is before capture of the bubble by the ring where the bubble is drawn into the low-pressure vortex core, while in stage II the bubble is stretched in the azimuthal direction within the ring and gradually broken up into a number of smaller bubbles. Following this, in stage III the bubble break-up is complete and the resulting smaller bubbles slowly move around the core, and finally in stage IV the bubbles escape. Apart from the effect of the ring on the bubble, the bubble is also shown to significantly affect the vortex ring, especially at low We (We similar to 3). In these low-We cases, the convection speed drops significantly compared to the base case without a bubble, while the core appears to fragment with a resultant large decrease in enstrophy by about 50 %. In the higher-We cases (We > 100), there are some differences in convection speed and enstrophy, but the effects are relatively small. The most dramatic effects of the bubble on the ring are found for thicker core rings at low We (We similar to 3) with the vortex ring almost stopping after interacting with the bubble, and the core fragmenting into two parts. The present idealized experiments exhibit many phenomena also seen in bubbly turbulent flows such as reduction in enstrophy, suppression of structures, enhancement of energy at small scales and reduction in energy at large scales. These similarities suggest that results from the present experiments can be helpful in better understanding interactions of bubbles with eddies in turbulent flows.

Inner Bound on the GDOF of the K-User MIMO Gaussian Symmetric Interference Channel

The K-user multiple input multiple output (MIMO) Gaussian symmetric interference channel where each transmitter has M antennas and each receiver has N antennas is studied from a generalized degrees of freedom (GDOF) perspective. An inner bound on the GDOF is derived using a combination of techniques such as treating interference as noise, zero forcing (ZF) at the receivers, interference alignment (IA), and extending the Han-Kobayashi (HK) scheme to K users, as a function of the number of antennas and the log INR/log SNR level. Several interesting conclusions are drawn from the derived bounds. It is shown that when K > N/M + 1, a combination of the HK and IA schemes performs the best among the schemes considered. When N/M < K <= N/M + 1, the HK-scheme outperforms other schemes and is found to be GDOF optimal in many cases. In addition, when the SNR and INR are at the same level, ZF-receiving and the HK-scheme have the same GDOF performance.

Outer Bounds on the Sum Rate of the K-User MIMO Gaussian Interference Channel

This paper derives outer bounds on the sum rate of the K-user MIMO Gaussian interference channel (GIC). Three outer bounds are derived, under different assumptions of cooperation and providing side information to receivers. The novelty in the derivation lies in the careful selection of side information, which results in the cancellation of the negative differential entropy terms containing signal components, leading to a tractable outer bound. The overall outer bound is obtained by taking the minimum of the three outer bounds. The derived bounds are simplified for the MIMO Gaussian symmetric IC to obtain outer bounds on the generalized degrees of freedom (GDOF). The relative performance of the bounds yields insight into the performance limits of multiuser MIMO GICs and the relative merits of different schemes for interference management. These insights are confirmed by establishing the optimality of the bounds in specific cases using an inner bound on the GDOF derived by the authors in a previous work. It is also shown that many of the existing results on the GDOF of the GIC can be obtained as special cases of the bounds, e. g., by setting K = 2 or the number of antennas at each user to 1.

HYDROMAGNETICS OF ADVECTIVE ACCRETION FLOWS AROUND BLACK HOLES: REMOVAL OF ANGULAR MOMENTUM BY LARGE-SCALE MAGNETIC STRESSES

We show that the removal of angular momentum is possible in the presence of large-scale magnetic stresses in geometrically thick, advective, sub-Keplerian accretion flows around black holes in steady state, in the complete absence of alpha-viscosity. The efficiency of such an angular momentum transfer could be equivalent to that of alpha-viscosity with alpha = 0.01-0.08. Nevertheless, the required field is well below its equipartition value, leading to a magnetically stable disk flow. This is essentially important in order to describe the hard spectral state of the sources when the flow is non/sub-Keplerian. We show in our simpler 1.5 dimensional, vertically averaged disk model that the larger the vertical-gradient of the azimuthal component of the magnetic field is, the stronger the rate of angular momentum transfer becomes, which in turn may lead to a faster rate of outflowing matter. Finding efficient angular momentum transfer in black hole disks via magnetic stresses alone, is very interesting when the generic origin of alpha-viscosity is still being explored.

Origin of nonlinearity and plausible turbulence by hydromagnetic transient growth in accretion disks: Faster growth rate than magnetorotational instability

We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that a mechanism, known as magnetorotational instability (MRI), is responsible for transporting matter in the presence of a weak magnetic field. However, there are some shortcomings, which question the effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g., transient growth (TG), can play an important role in bringing nonlinearity into the system, even at weak magnetic fields. In addition, it should be determined whether MRI or TG is primarily responsible for revealing nonlinearity in order to make the flow turbulent. Our results prove explicitly that the flows with a high Reynolds number (Re), which is the case for realistic astrophysical accretion disks, exhibit nonlinearity via TG of perturbation modes faster than that by modes producing MRI. For a fixed wave vector, MRI dominates over transient effects only at low Re, lower than the value expected to be in astrophysical accretion disks, and low magnetic fields. This calls into serious question the (overall) persuasiveness of MRI in astrophysical accretion disks.

Joint Channel Estimation and Data Detection in MIMO-OFDM Systems: A Sparse Bayesian Learning Approach

The impulse response of wireless channels between the N-t transmit and N-r receive antennas of a MIMO-OFDM system are group approximately sparse (ga-sparse), i.e., NtNt the channels have a small number of significant paths relative to the channel delay spread and the time-lags of the significant paths between transmit and receive antenna pairs coincide. Often, wireless channels are also group approximately cluster-sparse (gac-sparse), i.e., every ga-sparse channel consists of clusters, where a few clusters have all strong components while most clusters have all weak components. In this paper, we cast the problem of estimating the ga-sparse and gac-sparse block-fading and time-varying channels in the sparse Bayesian learning (SBL) framework and propose a bouquet of novel algorithms for pilot-based channel estimation, and joint channel estimation and data detection, in MIMO-OFDM systems. The proposed algorithms are capable of estimating the sparse wireless channels even when the measurement matrix is only partially known. Further, we employ a first-order autoregressive modeling of the temporal variation of the ga-sparse and gac-sparse channels and propose a recursive Kalman filtering and smoothing (KFS) technique for joint channel estimation, tracking, and data detection. We also propose novel, parallel-implementation based, low-complexity techniques for estimating gac-sparse channels. Monte Carlo simulations illustrate the benefit of exploiting the gac-sparse structure in the wireless channel in terms of the mean square error (MSE) and coded bit error rate (BER) performance.

Turbulent states and their transitions in mathematical models for ventricular tissue: The effects of random interstitial fibroblasts

We study the dynamical behaviors of two types of spiral-and scroll-wave turbulence states, respectively, in two-dimensional (2D) and three-dimensional (3D) mathematical models, of human, ventricular, myocyte cells that are attached to randomly distributed interstitial fibroblasts; these turbulence states are promoted by (a) the steep slope of the action-potential-duration-restitution (APDR) plot or (b) early afterdepolarizations (EADs). Our single-cell study shows that (1) the myocyte-fibroblast (MF) coupling G(j) and (2) the number N-f of fibroblasts in an MF unit lower the steepness of the APDR slope and eliminate the EAD behaviors of myocytes; we explore the pacing dependence of such EAD suppression. In our 2D simulations, we observe that a spiral-turbulence (ST) state evolves into a state with a single, rotating spiral (RS) if either (a) G(j) is large or (b) the maximum possible number of fibroblasts per myocyte N-f(max) is large. We also observe that the minimum value of G(j), for the transition from the ST to the RS state, decreases as N-f(max) increases. We find that, for the steep-APDR-induced ST state, once the MF coupling suppresses ST, the rotation period of a spiral in the RS state increases as (1) G(j) increases, with fixed N-f(max), and (2) N-f(max) increases, with fixed G(j). We obtain the boundary between ST and RS stability regions in the N-f(max)-G(j) plane. In particular, for low values of N-f(max), the value of G(j), at the ST-RS boundary, depends on the realization of the randomly distributed fibroblasts; this dependence decreases as N-f(max) increases. Our 3D studies show a similar transition from scroll-wave turbulence to a single, rotating, scroll-wave state because of the MF coupling. We examine the experimental implications of our study and propose that the suppression (a) of the steep slope of the APDR or (b) EADs can eliminate spiral-and scroll-wave turbulence in heterogeneous cardiac tissue, which has randomly distributed fibroblasts.

Two-Phase Flow With Phase Change in Porous Channels

Phase change heat transfer in porous media finds applications in various geological flows and modern heat pipes. We present a study to show the effect of phase change on heat transfer in a porous channel. We show that the ratio of Jakob numbers based on wall superheat and inlet fluid subcooling governs the liquid-vapor interface location in the porous channel and below a critical value of the ratio, the liquid penetrates all the way to the extent of the channel in the flow direction. In such cases, the Nusselt number is higher due to the proximity of the liquid-vapor interface to the heat loads. For higher heat loads or lower subcooling of the liquid, the liquid-vapor interface is pushed toward the inlet, and heat transfer occurs through a wider vapor region thus resulting in a lower Nusselt number. This study is relevant in the designing of efficient two-phase heat exchangers such as capillary suction based heat pipes where a prior estimation of the interface location for the maximum heat load is required to ensure that the liquid-vapor interface is always inside the porous block for its operation.