Taming Control Divergence in GPUs through Control Flow Linearization

Branch divergence is a very commonly occurring performance problem in GPGPU in which the execution of diverging branches is serialized to execute only one control flow path at a time. Existing hardware mechanism to reconverge threads using a stack causes duplicate execution of code for unstructured control flow graphs. Also the stack mechanism cannot effectively utilize the available parallelism among diverging branches. Further, the amount of nested divergence allowed is also limited by depth of the branch divergence stack. In this paper we propose a simple and elegant transformation to handle all of the above mentioned problems. The transformation converts an unstructured CFG to a structured CFG without duplicating user code. It incurs only a linear increase in the number of basic blocks and also the number of instructions. Our solution linearizes the CFG using a predicate variable. This mechanism reconverges the divergent threads as early as possible. It also reduces the depth of the reconvergence stack. The available parallelism in nested branches can be effectively extracted by scheduling the basic blocks to reduce the effect of stalls due to memory accesses. It can also increase execution efficiency of nested loops with different trip counts for different threads. We implemented the proposed transformation at PTX level using the Ocelot compiler infrastructure. We evaluated the technique using various benchmarks to show that it can be effective in handling the performance problem due to divergence in unstructured CFGs.

Stability of Poiseuille flow in a fluid overlying an anisotropic and inhomogeneous porous layer

We present the linear stability analysis of horizontal Poiseuille flow in a fluid overlying a porous medium with anisotropic and inhomogeneous permeability. The generalized Darcy model is used to describe the flow in the porous medium with the Beavers-Joseph condition at the interface of the two layers and the eigenvalue problem is solved numerically. The effect of major system parameters on the stability characteristics is addressed in detail. It is shown that the anisotropic and inhomogeneous modulation of the permeability of the underlying porous layer provides an effective means for passive control of the flow stability.

Field-Portable Microfluidics-Based Imaging Flow Cytometer

Clinical microscopy is a versatile and robust tool used for the diagnosis of a plethora of diseases. However, due to various reasons, it remains inaccessible in resource limited settings. In this paper, we present an automated and cost-effective alternative to microscopy for use in clinical diagnostics. With the use of custom optics and microfluidics, we demonstrate a field-portable imaging flow cytometry system. Using the presented system, we have been able to image 586 cells per second. We demonstrate the clinical relevance of the proposed system by differentiating between suspensions of healthy and sphered RBCs based on high-throughput morphometric analysis. The instrument presented here is a major advancement in the domain of field portable diagnostics as it enables fast and robust quantitative diagnostic testing at the point-of-care.

Deformability Measurement of Single-Cells at High-Throughput With Imaging Flow Cytometry

Disease conditions like malaria, sickle cell anemia, diabetes mellitus, cancer, etc., are known to significantly alter the deformability of certain types of cells (red blood cells, white blood cells, circulating tumor cells, etc.). To determine the cellular deformability, techniques like micropipette aspiration, atomic force microscopy, optical tweezers, quantitative phase imaging have been developed. Many of these techniques have an advantage of determining the single cell deformability with ultrahigh precision. However, the suitability of these techniques for the realization of a deformability based diagnostic tool is questionable as they are expensive and extremely slow to operate on a huge population of cells. In this paper, we propose a technique for high-throughput (800 cells/s) determination of cellular deformability on a single cell basis. This technique involves capturing the image(s) of cells in flow that have undergone deformation under the influence of shear gradient generated by the fluid flowing through the microfluidic channels. Deformability indices of these cells can be computed by performing morphological operations on these images. We demonstrate the applicability of this technique for examining the deformability index on healthy, diabetic, and sphered red blood cells. We believe that this technique has a strong role to play in the realization of a potential tool that uses deformability as one of the important criteria in disease diagnosis.

Investigation of schemes for incorporating generator Q limits in the fast decoupled load flow method

Fast Decoupled Load Flow (FDLF) is a very popular and widely used power flow analysis method because of its simplicity and efficiency. Even though the basic FDLF algorithm is well investigated, the same is not true in the case of additional schemes/modifications required to obtain adjusted load flow solutions using the FDLF method. Handling generator Q limits is one such important feature needed in any practical load flow method. This paper presents a comprehensive investigation of two classes of schemes intended to handle this aspect i.e. the bus type switching scheme and the sensitivity scheme. We propose two new sensitivity based schemes and assess their performance in comparison with the existing schemes. In addition, a new scheme to avoid the possibility of anomalous solutions encountered while using the conventional schemes is also proposed and evaluated. Results from extensive simulation studies are provided to highlight the strengths and weaknesses of these existing and proposed schemes, especially from the point of view of reliability.

Atmospheric pressure plasma chemical vapor deposition reactor for 100 mm wafers, optimized for minimum contamination at low gas flow rates

Gas discharge plasmas used for thinfilm deposition by plasma-enhanced chemical vapor deposition (PECVD) must be devoid of contaminants, like dust or active species which disturb the intended chemical reaction. In atmospheric pressure plasma systems employing an inert gas, the main source of such contamination is the residual air inside the system. To enable the construction of an atmospheric pressure plasma (APP) system with minimal contamination, we have carried out fluid dynamic simulation of the APP chamber into which an inert gas is injected at different mass flow rates. On the basis of the simulation results, we have designed and built a simple, scaled APP system, which is capable of holding a 100 mm substrate wafer, so that the presence of air (contamination) in the APP chamber is minimized with as low a flow rate of argon as possible. This is examined systematically by examining optical emission from the plasma as a function of inert gas flow rate. It is found that optical emission from the plasma shows the presence of atmospheric air, if the inlet argon flow rate is lowered below 300 sccm. That there is minimal contamination of the APP reactor built here, was verified by conducting an atmospheric pressure PECVD process under acetylene flow, combined with argon flow at 100 sccm and 500 sccm. The deposition of a polymer coating is confirmed by infrared spectroscopy. X-ray photoelectron spectroscopy shows that the polymer coating contains only 5% of oxygen, which is comparable to the oxygen content in polymer deposits obtained in low-pressure PECVD systems. (C) 2015 AIP Publishing LLC.

After transition in a soft-walled microchannel

In comparison to the flow in a rigid channel, there is a multifold reduction in the transition Reynolds number for the flow in a microchannel when one of the walls is made sufficiently soft, due to a dynamical instability induced by the fluid-wall coupling, as shown by Verma & Kumaran (J. Fluid Mech., vol. 727, 2013, pp. 407-455). The flow after transition is characterised using particle image velocimetry in the x-y plane, where x is the streamwise direction and y is the cross-stream coordinate along the small dimension of the channel of height 0.2-0.3 mm. The flow after transition is characterised by a mean velocity profile that is flatter at the centre and steeper at the walls in comparison to that for a laminar flow. The root mean square of the streamwise fluctuating velocity shows a characteristic sharp increase away from the wall and a maximum close to the wall, as observed in turbulent flows in rigid-walled channels. However, the profile is asymmetric, with a significantly higher maximum close to the soft wall in comparison to that close to the hard wall, and the Reynolds stress is found to be non-zero at the soft wall, indicating that there is a stress exerted by fluid velocity fluctuations on the wall. The maximum of the root mean square of the velocity fluctuations and the Reynolds stress (divided by the fluid density) in the soft-walled microchannel for Reynolds numbers in the range 250-400, when scaled by suitable powers of the maximum velocity, are comparable to those in a rigid channel at Reynolds numbers in the range 5000-20 000. The near-wall velocity profile shows no evidence of a viscous sublayer for (y upsilon(*)/nu) as low as two, but there is a logarithmic layer for (y upsilon(*)/nu) up to approximately 30, where the von Karman constants are very different from those for a rigid-walled channel. Here, upsilon(*) is the friction velocity, nu is the kinematic viscosity and y is the distance from the soft surface. The surface of the soft wall in contact with the fluid is marked with dye spots to monitor the deformation and motion along the fluid-wall interface. Low-frequency oscillations in the displacement of the surface are observed after transition in both the streamwise and spanwise directions, indicating that the velocity fluctuations are dynamically coupled to motion in the solid.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.

DTM AT 25: ESSAYS ON THEMES AND FUTURE DIRECTIONS

This paper describes the development and evolution of research themes in the Design Theory and Methodology (DTM) conference. Essays containing reflections on the history of DTM, supported by an analysis of session titles and papers winning the ``best paper award'', describe the development of the research themes. A second set of essays describes the evolution of several key research themes. Two broad trends in research themes are evident, with a third one emerging. The topics of the papers in the first decade or so reflect an underlying aim to apply artificial intelligence toward developing systems that could `design'. To do so required understanding how human designers behave, formalizing design processes so that they could be computed, and formalizing representations of design knowledge. The themes in the first DTM conference and the recollections of the DTM founders reflect this underlying aim. The second decade of DTM saw the emergence of product development as an underlying concern and included a growth in a systems view of design. More recently, there appears to be a trend toward design-led innovation, which entails both executing the design process more efficiently and understanding the characteristics of market-leading designs so as to produce engineered products and systems of exceptional levels of quality and customer satisfaction.

Flow Cytometry Analysis of Cytotoxicity In Vitro and Long-Term Toxicity of HA-40wt% BaTiO3 Nanoparticles In Vivo

The development of new implantable biomaterials requires bone-mimicking physical properties together with desired biocompatible property. In continuation to our earlier published research to establish compositional dependent multifunctional bone-like properties and cytocompatibility response of hydroxyapatite (HA)-BaTiO3 composites, the toxicological property evaluation, both invitro and invivo, were conducted on HA-40wt% BaTiO3 and reported in this work. In particular, this work reports invitro cytotoxicity of mouse myoblast cells as well as invivo long-term tissue and nanoparticles interaction of intra-articularly injected HA-40wt% BaTiO3 and BaTiO3 up to the concentration of 25mg/mL in physiological saline over 12weeks in mouse model. The careful analysis of flow cytometry results could not reveal any statistically significant difference in terms of early/late apoptotic cells or necrotic cells over 8d in culture. Extensive histological analysis could not record any signature of cellular level toxicity or pronounced inflammatory response in vital organs as well as at knee joints of Balb/c mice after 12weeks. Taken together, this study establishes nontoxic nature of HA-40wt% BaTiO3 and therefore, HA-40wt% BaTiO3 can be used safely for various biomedical applications.

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.

Effect of Channel Confinement on Mixed Convective Flow Past an Equilateral Triangular Cylinder

The present work investigates the mixed convective flow and heat transfer characteristics past a triangular cylinder placed symmetrically in a vertical channel. At a representative Reynolds number, Re = 100, simulations are carried out for the blockage ratios beta = 1/3; 1/4; and 1/6. Effect of aiding and opposing buoyancy is brought about by varying the Richardson number in the range -1.0 <= Ri <= 1.0. At a blockage ratio of 1/3, suppression of vortex shedding is found at Ri = 1, whereas von Karman vortex street is seen both at beta = 1/4 and 1/6, respectively. This is the first time that such behavior of blockage ratio past a triangular cylinder in the present flow configuration is reported. Drag coefficient increases progressively with increasing Ri and a slightly higher value is noticed at beta = 1/3. For all b, heat transfer increases with increasing Ri. Flattening of Nu(avg)-Ri curve beyond Ri > 0: 75 is observed at beta = 1/3.

Shock tunnel measurements of surface pressures in shock induced separated flow field using MEMS sensor array

Characterized not just by high Mach numbers, but also high flow total enthalpies-often accompanied by dissociation and ionization of flowing gas itself-the experimental simulation of hypersonic flows requires impulse facilities like shock tunnels. However, shock tunnel simulation imposes challenges and restrictions on the flow diagnostics, not just because of the possible extreme flow conditions, but also the short run times-typically around 1 ms. The development, calibration and application of fast response MEMS sensors for surface pressure measurements in IISc hypersonic shock tunnel HST-2, with a typical test time of 600 mu s, for the complex flow field of strong (impinging) shock boundary layer interaction with separation close to the leading edge, is delineated in this paper. For Mach numbers 5.96 (total enthalpy 1.3 MJ kg(-1)) and 8.67 (total enthalpy 1.6 MJ kg(-1)), surface pressures ranging from around 200 Pa to 50 000 Pa, in various regions of the flow field, are measured using the MEMS sensors. The measurements are found to compare well with the measurements using commercial sensors. It was possible to resolve important regions of the flow field involving significant spatial gradients of pressure, with a resolution of 5 data points within 12 mm in each MEMS array, which cannot be achieved with the other commercial sensors. In particular, MEMS sensors enabled the measurement of separation pressure (at Mach 8.67) near the leading edge and the sharply varying pressure in the reattachment zone.

Effect of absolute pressure on flow through a textured hydrophobic microchannel

The potential of textured hydrophobic surfaces to provide substantial drag reduction has been attributed to the presence of air bubbles trapped on the surface cavities. In this paper, we present results on water flow past a textured hydrophobic surface, while systematically varying the absolute pressure close to the surface. Trapped air bubbles on the surface are directly visualized, along with simultaneous pressure drop measurements across the surface in a microchannel configuration. We find that varying the absolute pressure within the channel greatly influences the trapped air bubble behavior, causing a consequent effect on the pressure drop (drag). When the absolute pressure within the channel is maintained below atmospheric pressure, we find that the air bubbles grow in size, merge and eventually detach from the surface. This growth and subsequent merging of the air bubbles leads to a substantial increase in the pressure drop. On the other hand, a pressure above the atmospheric pressure within the channel leads to gradual shrinkage and eventual disappearance of trapped air bubbles. We find that in this case, air bubbles do cause reduction in the pressure drop with the minimum pressure drop (or maximum drag reduction) occurring when the bubbles are flush with the surface. These results show that the trapped air bubble dynamics and the pressure drop across a textured hydrophobic microchannel are very significantly dependent on the absolute pressure within the channel. The results obtained hold important implications toward achieving sustained drag reduction in microfluidic applications.

An alternative approach to study nonlinear inviscid flow over arbitrary bottom topography

This paper deals with a new approach to study the nonlinear inviscid flow over arbitrary bottom topography. The problem is formulated as a nonlinear boundary value problem which is reduced to a Dirichlet problem using certain transformations. The Dirichlet problem is solved by applying Plemelj-Sokhotski formulae and it is noticed that the solution of the Dirichlet problem depends on the solution of a coupled Fredholm integral equation of the second kind. These integral equations are solved numerically by using a modified method. The free-surface profile which is unknown at the outset is determined. Different kinds of bottom topographies are considered here to study the influence of bottom topography on the free-surface profile. The effects of the Froude number and the arbitrary bottom topography on the free-surface profile are demonstrated in graphical forms for the subcritical flow. Further, the nonlinear results are validated with the results available in the literature and compared with the results obtained by using linear theory. (C) 2015 Elsevier Inc. All rights reserved.