Planform structure of turbulent free convection on horizontal surfaces

The planform structure of turbulent free convection over a heated horizontal surface has been visualized and analyzed for different boundary conditions at the top and for different aspect ratios, for flux Rayleigh numbers ranging from 10 exp 8 - 10 exp 10. The different boundary conditions correspond to Rayleigh-Benard convection, open convection with evaporation at the top and with an imposed external flow on the heated boundary. Without the external flow the planform is one randomly oriented line plume. At large Ra, these line plumes seem to align along the diagonal, persumably due to a large-scale flow along as visualized in the side view. When the external flow is imposed, the line plumes clearly align in the direction of external flow. Flow visualization reveals that at these Ra, the shear tends to break the plumes which otherwise would reach the opposite boundary. (Author)

Numerical simulation of nano-indentation on rough surfaces

Nano-indentation is a technique used to measure various mechanical properties like hardness, Young's modulus and the adherence of thin films and surface layers. It can be used as a quality control tool for various surface modification techniques like ion-implantation, film deposition processes etc. It is important to characterise the increasing scatter in the data measured at lower penetration depths observed in the nano-indentation, for the technique to be effectively applied. Surface roughness is one of the parameters contributing for the scatter. This paper is aimed at quantifying the nature and the amount of scatter that will be introduced in the measurement due to the roughness of the surface on which the indentation is carried out. For this the surface is simulated using the Weierstrass-Mandelbrot function which gives a self-affine fractal. The contact area of this surface with a conical indenter with a spherical cap at the tip is measured numerically. The indentation process is simulated using the spherical cavity model. This eliminates the indentation size effect observed at the micron and sub-micron scales. It has been observed that there exists a definite penetration depth in relation to the surface roughness beyond which the scatter is reduced such that reliable data could be obtained.

Plumes Dynamics and Heat Transfer over Horizontal Grooved Surfaces

This study investigates the free convection and plumes dynamics over horizontal surfaces with parallel V-grooves. The convection is studied in a tank of water with the bottom surface being a smooth or grooved surface and the top of the water surface exposed to ambient. Two groove heights were used-10 mm and 3 mm-and the experiment was done with two values of aspect ratio-2.9 and 1.8 (aspect ratio is the width of the fluid layer/height of fluid layer). Heat flux at the bottom surface was from electrical heating. Beyond a certain critical temperature difference, enhanced heat transfer is obtained on the grooved surface compared to a smooth surface. Nusselt numbers are evaluated for both smooth and grooved surfaces and correlated using modified Rayleigh numbers. Visualization shows that the enhanced heat transport in the rough cavities cannot be ascribed to the increase in the contact area; rather, it must be the local dynamics of the thermal boundary layer.

Flow visualization in turbulent free convection over horizontal smooth and grooved surfaces

The present paper discusses the flow visualization for turbulent free convection in a tank of water with the bottom surface being a smooth or a grooved surface and the top of the water surface exposed to ambient. The grooved surface is of parallel 90 degrees V-grooves with groove height of 10 mm and groove width of 20 mm. The experiment is carried out with aspect ratio (AR) of 2.9 and Rayleigh number (Ra) in the range, 1.3 x 10(7) - 4 x 10(7). Here AR is the aspect ratio (= width of fluid layer/height of fluid layer). Heat flux at the bottom surface is from electrical heating. From the pH-dye visualization, interesting flow structures are observed and these structures are analyzed with the help of plumes dynamics and temperature variations with time. (C) 2011 Elsevier Ltd. All rights reserved.

Junction between surfaces of two topological insulators

We study the properties of a line junction which separates the surfaces of two three-dimensional topological insulators. The velocities of the Dirac electrons on the two surfaces may be unequal and may even have opposite signs. For a time-reversal invariant system, we show that the line junction is characterized by an arbitrary parameter alpha which determines the scattering from the junction. If the surface velocities have the same sign, we show that there can be edge states which propagate along the line junction with a velocity and spin orientation which depend on alpha and the ratio of the velocities. Next, we study what happens if the two surfaces are at an angle phi with respect to each other. We study the scattering and differential conductance through the line junction as functions of phi and alpha. We also find that there are edge states which propagate along the line junction with a velocity and spin orientation which depend on phi. Finally, if the surface velocities have opposite signs, we find that the electrons must transmit into the two-dimensional interface separating the two topological insulators.

Mixed convection adjacent to non-isothermal vertical surfaces

Study of laminar boundary layer in mixed convection from vertical plates is carried out. The surface temperature along the vertical plate is assumed to vary arbitrarily with vertical distance. Perturbation technique is used to solve the governing boundary layer equations. The differentials of the wall temperature are used as perturbation elements, which are functions of vertical distance, to obtain universal functions. The universal functions are valid for any type of vertical wall temperature variation. Heat transfer rates and fluid velocity inside the boundary layer can be expressed and calculated using these universal functions. Heat transfer rates are obtained for the special cases of power-law variation of the wall temperature. The effect of the governing parameter (Gr(y)/Re-y(2)) and the power index of the power-law wall temperature variation on heat transfer rates is studied. For the purpose of validation, the mixed convection results obtained by the present technique pertaining to the special cases of isothermal vertical wall are compared with those obtained by similarity analysis reported in literature, and the agreement is found to be good. (C) 2012 Elsevier Ltd. All rights reserved.

Directional motion of impacting drops on dual-textured surfaces

In this work, we analyze the directional movement of impacting liquid drops on dual-textured solid surfaces comprising two different surface morphologies: a textured surface and a smooth surface. The dynamics of liquid drops impacting onto the junction line between the two parts of the dual-textured surfaces is studied experimentally for varying drop impact velocity. The dual-textured surfaces used here featured a variation in their textures' geometrical parameters as well as their surface chemistry. Two types of liquid drop differing in their surface tension were used. The impact process develops a net horizontal drop velocity towards the higher-wettability surface portion and results in a bulk movement of the impacting drop liquid. The final distance moved by the impacting drop from the junction line decreases with increasing impacting drop Weber number We. A fully theoretical model, employing a balance of forces acting at the drop contact line as well as energy conservation, is formulated to determine the variation, with We, of net horizontal drop velocity and subsequent movement of the impacting drop on the dual-textured surfaces.

Combining heterogeneous classifiers for relational databases

Practical usage of machine learning is gaining strategic importance in enterprises looking for business intelligence. However, most enterprise data is distributed in multiple relational databases with expert-designed schema. Using traditional single-table machine learning techniques over such data not only incur a computational penalty for converting to a flat form (mega-join), even the human-specified semantic information present in the relations is lost. In this paper, we present a practical, two-phase hierarchical meta-classification algorithm for relational databases with a semantic divide and conquer approach. We propose a recursive, prediction aggregation technique over heterogeneous classifiers applied on individual database tables. The proposed algorithm was evaluated on three diverse datasets. namely TPCH, PKDD and UCI benchmarks and showed considerable reduction in classification time without any loss of prediction accuracy. (C) 2012 Elsevier Ltd. All rights reserved.

3-MANIFOLD GROUPS, KÄHLER GROUPS AND COMPLEX SURFACES

Let G be a Kahler group admitting a short exact sequence 1 -> N -> G -> Q -> 1 where N is finitely generated. (i) Then Q cannot be non-nilpotent solvable. (ii) Suppose in addition that Q satisfies one of the following: (a) Q admits a discrete faithful non-elementary action on H-n for some n >= 2. (b) Q admits a discrete faithful non-elementary minimal action on a simplicial tree with more than two ends. (c) Q admits a (strong-stable) cut R such that the intersection of all conjugates of R is trivial. Then G is virtually a surface group. It follows that if Q is infinite, not virtually cyclic, and is the fundamental group of some closed 3-manifold, then Q contains as a finite index subgroup either a finite index subgroup of the three-dimensional Heisenberg group or the fundamental group of the Cartesian product of a closed oriented surface of positive genus and the circle. As a corollary, we obtain a new proof of a theorem of Dimca and Suciu in Which 3-manifold groups are Kahler groups? J. Eur. Math. Soc. 11 (2009) 521-528] by taking N to be the trivial group. If instead, G is the fundamental group of a compact complex surface, and N is finitely presented, then we show that Q must contain the fundamental group of a Seifert-fibered 3-manifold as a finite index subgroup, and G contains as a finite index subgroup the fundamental group of an elliptic fibration. We also give an example showing that the relation of quasi-isometry does not preserve Kahler groups. This gives a negative answer to a question of Gromov which asks whether Kahler groups can be characterized by their asymptotic geometry.

Automatic compilation of MATLAB programs for synergistic execution on heterogeneous processors

MATLAB is an array language, initially popular for rapid prototyping, but is now being increasingly used to develop production code for numerical and scientific applications. Typical MATLAB programs have abundant data parallelism. These programs also have control flow dominated scalar regions that have an impact on the program's execution time. Today's computer systems have tremendous computing power in the form of traditional CPU cores and throughput oriented accelerators such as graphics processing units(GPUs). Thus, an approach that maps the control flow dominated regions to the CPU and the data parallel regions to the GPU can significantly improve program performance. In this paper, we present the design and implementation of MEGHA, a compiler that automatically compiles MATLAB programs to enable synergistic execution on heterogeneous processors. Our solution is fully automated and does not require programmer input for identifying data parallel regions. We propose a set of compiler optimizations tailored for MATLAB. Our compiler identifies data parallel regions of the program and composes them into kernels. The problem of combining statements into kernels is formulated as a constrained graph clustering problem. Heuristics are presented to map identified kernels to either the CPU or GPU so that kernel execution on the CPU and the GPU happens synergistically and the amount of data transfer needed is minimized. In order to ensure required data movement for dependencies across basic blocks, we propose a data flow analysis and edge splitting strategy. Thus our compiler automatically handles composition of kernels, mapping of kernels to CPU and GPU, scheduling and insertion of required data transfer. The proposed compiler was implemented and experimental evaluation using a set of MATLAB benchmarks shows that our approach achieves a geometric mean speedup of 19.8X for data parallel benchmarks over native execution of MATLAB.

Time dependent superhydrophobicity of drag reducing surfaces

Air can be trapped on the crevices of specially textured hydrophobic surfaces immersed in water. This heterogenous state of wetting in which the water is in contact with both the solid surface and the entrapped air is not stable. Diffusion of air into the surrounding water leads to gradual reduction in the size and numbers of the air bubbles. The sustainability of the entrapped air on such surfaces is important for many underwater applications in which the surfaces have to remain submersed for longer time periods. In this paper we explore the suitability of different classes of surface textures towards the drag reduction application by evaluating the time required for the disappearance of the air bubbles under hydrostatic conditions. Different repetitive textures consisting of holes, pillars and ridges of different sizes have been generated in silicon, aluminium and brass by isotropic etching, wire EDM and chemical etching respectively. These surfaces were rendered hydrophobic with self-assembled layer of fluorooctyl trichlorosilane for silicon and aluminium surfaces and 1-dodecanethiol for brass surfaces. Using total internal reflection the air bubbles are visualized with the help of a microscope and time lapse photography. Irrespective of the texture, both the size and the number of air pockets were found to decrease with time gradually and eventually disappear. In an attempt to reverse the diffusion we explore the possibility of using electrolysis to generate gases at the textured surfaces. The gas bubbles are nucleated everywhere on the surface and as they grow they coalesce with each other and get pinned at the texture edges.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Runtime dependence computation and execution of loops on heterogeneous systems

GPUs have been used for parallel execution of DOALL loops. However, loops with indirect array references can potentially cause cross iteration dependences which are hard to detect using existing compilation techniques. Applications with such loops cannot easily use the GPU and hence do not benefit from the tremendous compute capabilities of GPUs. In this paper, we present an algorithm to compute at runtime the cross iteration dependences in such loops. The algorithm uses both the CPU and the GPU to compute the dependences. Specifically, it effectively uses the compute capabilities of the GPU to quickly collect the memory accesses performed by the iterations by executing the slice functions generated for the indirect array accesses. Using the dependence information, the loop iterations are levelized such that each level contains independent iterations which can be executed in parallel. Another interesting aspect of the proposed solution is that it pipelines the dependence computation of the future level with the actual computation of the current level to effectively utilize the resources available in the GPU. We use NVIDIA Tesla C2070 to evaluate our implementation using benchmarks from Polybench suite and some synthetic benchmarks. Our experiments show that the proposed technique can achieve an average speedup of 6.4x on loops with a reasonable number of cross iteration dependences.

Heterogeneous locational optimisation using a generalised Voronoi partition

In this paper a generalisation of the Voronoi partition is used for locational optimisation of facilities having different service capabilities and limited range or reach. The facilities can be stationary, such as base stations in a cellular network, hospitals, schools, etc., or mobile units, such as multiple unmanned aerial vehicles, automated guided vehicles, etc., carrying sensors, or mobile units carrying relief personnel and materials. An objective function for optimal deployment of the facilities is formulated, and its critical points are determined. The locally optimal deployment is shown to be a generalised centroidal Voronoi configuration in which the facilities are located at the centroids of the corresponding generalised Voronoi cells. The problem is formulated for more general mobile facilities, and formal results on the stability, convergence and spatial distribution of the proposed control laws responsible for the motion of the agents carrying facilities, under some constraints on the agents' speed and limit on the sensor range, are provided. The theoretical results are supported with illustrative simulation results.

Selenium adsorption on Au(111) and Ag(111) surfaces: adsorbed selenium and selenide films

Results of a high resolution photoemission and electrochemistry study of Se adsorption Au(111) and Ag(111) surfaces performed by immersion of pristine samples into an aqeuous solution of Na2Se are presented. Cyclic voltammetry on Au shows formation of selenium adsorbed species and the structures observed in reductive desorption are to the atomic and polymeric species observed in XPS. In the case of Au(111) XPS spectra in the Se(3d) region indeed show two main features attributed to Se chemisorbed atomically and polymeric Se-8 features.' Smaller structures due to other types of Se conformations were also observed. The Au(4f) peak line, shape does not show core level, shifts: indicative of Au selenide formation the case of silver, XPS spectra for the Ag(3d) show a broadening of the peak and a deconvolution into Ag-B bulk like Ag-Se components shows that the Ag-Se is located at a lower binding energy, an effect similar to oxidation and sulfidation of Ag. The Se(3d) XPS spectrum is found to be substantially different from the Au case and dominated by atomic type Se due to the selenide, though a smaller intensity Se structure at an energy similar to the Se-8 structure for Au is also observed. Changes in the valence band region. related to Se adsorption are reported.