Improving GPGPU Concurrency with Elastic Kernels

Each new generation of GPUs vastly increases the resources available to GPGPU programs. GPU programming models (like CUDA) were designed to scale to use these resources. However, we find that CUDA programs actually do not scale to utilize all available resources, with over 30% of resources going unused on average for programs of the Parboil2 suite that we used in our work. Current GPUs therefore allow concurrent execution of kernels to improve utilization. In this work, we study concurrent execution of GPU kernels using multiprogram workloads on current NVIDIA Fermi GPUs. On two-program workloads from the Parboil2 benchmark suite we find concurrent execution is often no better than serialized execution. We identify that the lack of control over resource allocation to kernels is a major serialization bottleneck. We propose transformations that convert CUDA kernels into elastic kernels which permit fine-grained control over their resource usage. We then propose several elastic-kernel aware concurrency policies that offer significantly better performance and concurrency compared to the current CUDA policy. We evaluate our proposals on real hardware using multiprogrammed workloads constructed from benchmarks in the Parboil 2 suite. On average, our proposals increase system throughput (STP) by 1.21x and improve the average normalized turnaround time (ANTT) by 3.73x for two-program workloads when compared to the current CUDA concurrency implementation.

Parametrisation-free determination of the shape parameters for the pion electromagnetic form factor

Recent data from high-statistics experiments that have measured the modulus of the pion electromagnetic form factor from threshold to relatively high energies are used as input in a suitable mathematical framework of analytic continuation to find stringent constraints on the shape parameters of the form factor at t = 0. The method uses also as input a precise description of the phase of the form factor in the elastic region based on Fermi-Watson theorem and the analysis of the pi pi scattering amplitude with dispersive Roy equations, and some information on the spacelike region coming from recent high precision experiments. Our analysis confirms the inconsistencies of several data on the modulus, especially from low energies, with analyticity and the input phase, noted in our earlier work. Using the data on the modulus from energies above 0.65 GeV, we obtain, with no specific parametrisation, the prediction < r(pi)(2)> is an element of (0.42, 0.44) fm(2) for the charge radius. The same formalism leads also to very narrow allowed ranges for the higher-order shape parameters at t = 0, with a strong correlation among them.

MODELING AND SIMULATION OF HYDRODYNAMIC INTERACTION OF DNA IN A MICRO-FLUIDIC CHANNEL

The motion of DNA (in the bulk solution) and the non-Newtonian effective fluid behavior are considered separately and self-consistently with the fluid motion satisfying the no-slip boundary condition on the surface of the confining geometry in the presence of channel pressure gradients. A different approach has been developed to model DNA in the micro-channel. In this study the DNA is assumed as an elastic chain with its characteristic Young's modulus, Poisson's ratio and density. The force which results from the fluid dynamic pressure, viscous forces and electromotive forces is applied to the elastic chain in a coupled manner. The velocity fields in the micro-channel are influenced by the transport properties. Simulations are carried out for the DNAs attached to the micro-fluidic wall. Numerical solutions based on a coupled multiphysics finite element scheme are presented. The modeling scheme is derived based on mass conservation including biomolecular mass, momentum balance including stress due to Coulomb force field and DNA-fluid interaction, and charge transport associated to DNA and other ionic complexes in the fluid. Variation in the velocity field for the non-Newtonian flow and the deformation of the DNA strand which results from the fluid-structure interaction are first studied considering a single DNA strand. Motion of the effective center of mass is analyzed considering various straight and coil geometries. Effects of DNA statistical parameters (geometry and spatial distribution of DNAs along the channel) on the effective flow behavior are analyzed. In particular, the dynamics of different DNA physical properties such as radius of gyration, end-to-end length etc. which are obtained from various different models (Kratky-Porod, Gaussian bead-spring etc.) are correlated to the nature of interaction and physical properties under the same background fluid environment.

Elastic-plastic contact conditions for frictionally constrained bodies under cyclic tangential loading

A Frictionally constrained condition implies dependence of friction force on tangential displacement amplitude. The condition may occur due to chemical, physical, and/or mechanical interaction between the contacting surfaces. The condition, sometimes also referred to as the presliding condition or partial slip condition, is characterized under fretting. Under such conditions, various experimental studies indicate the existence of two distinguishable regions, that is, stick region and slip region. In the present study, frictionally constrained conditions are identified and the evolutions of stick-slip regions are investigated in detail. Investigations have been performed on self-mated stainless steel and chromium carbide coated surfaces mated against stainless steel, under both vacuum and ambient conditions. Contact conditions prevailing at the contact interface were identified based on the mechanical responses and were correlated with the surface damage observed. Surface degradation has been observed in the form of microcracks and material transfer. Detailed numerical analysis has also been performed in order to understand the energy dissipation and the damage mode involved in the surface or subsurface damage. It has been observed that under frictionally constrained conditions, the occurrence of annular slip features are mainly due to the junction growth, resulting from elastic-plastic deformation at the contact interface. Ratcheting has been observed as the governing damage mode under cyclic tangential loading condition.

Elastic thickness structure of the Andaman subduction zone: Implications for convergence of the Ninetyeast Ridge

We use the Bouguer coherence (Morlet isostatic response function) technique to compute the spatial variation of effective elastic thickness (T-e) of the Andaman subduction zone. The recovered T-e map resolves regional-scale features that correlate well with known surface structures of the subducting Indian plate and the overriding Burma plate. The major structure on the India plate, the Ninetyeast Ridge (NER), exhibits a weak mechanical strength, which is consistent with the expected signature of an oceanic ridge of hotspot origin. However, a markedly low strength (0< T-e <3 km) in that region, where the NER is close to the Andaman trench (north of 10 N), receives our main attention in this study. The subduction geometry derived from the Bouguer gravity forward modeling suggests that the NER has indented beneath the Andaman arc. We infer that the bending stresses of the viscous plate, which were reinforced within the subducting oceanic plate as a result of the partial subduction of the NER buoyant load, have reduced the lithospheric strength. The correlation, T-e < T-s (seismogenic thickness) reveals that the upper crust is actively deforming beneath the frontal arc Andaman region. The occurrence of normal-fault earthquakes in the frontal arc, low Te zone, is indicative of structural heterogeneities within the subducting plate. The fact that the NER along with its buoyant root is subducting under the Andaman region is inhibiting the subduction processes, as suggested by the changes in trench line, interrupted back-arc volcanism, variation in seismicity mechanism, slow subduction, etc. The low T-e and thinned crustal structure of the Andaman back-arc basin are attributed to a thermomechanically weakened lithosphere. The present study reveals that the ongoing back-arc spreading and strike-slip motion along the West Andaman Fault coupled with the ridge subduction exerts an important control on the frequency and magnitude of seismicity in the Andaman region. (C) 2013 Elsevier Ltd. All rights reserved.

Transverse elastic response of bundled coated cylinders

The overall elastic response of a bundle of coated cylinders is a major aspect of thermal, nuclear and automotive engineering designs. This paper extends the previous work on tubular bundles to assess the effect of coating material and thickness. A major contribution from this paper is determining the overall transverse elastic response of coated thick cylinders by extending the Michell stress function approach in conjunction with contact mechanics. Finite element results using contact elements pave the way for applying the contact stress boundary conditions for Michell analysis. Theoretical and finite element analyses overall give results consistent with the previous work, and the results also fall within the well-established Voigt-Reuss bounds. (C) 2013 Elsevier Ltd. All rights reserved.

Energy Loss in the Impact of Elastic Spheres on a Rigid Half-Space in Presence of Adhesion

Adhesion can cause energy losses in asperities or particles coming into dynamic contact resulting in frictional dissipation, even if the deformation occurring is purely elastic. Such losses are of special significance in impact of nanoparticles and friction between surfaces under low contact pressure to hardness ratio. The objective of this work is to study the effect of adhesion during the normal impact of elastic spheres on a rigid half-space, with an emphasis on understanding the mechanism of energy loss. We use finite element method for modeling the impact phenomenon, with the adhesion due to van der Waals force and the short-range repulsion included as body forces distributed over the volume of the sphere. This approach, in contrast with commonly used surface force approximation, helps to model the interactions in a more precise way. We find that the energy loss in impact of elastic spheres is negligible unless there are adhesion-induced instabilities. Significant energy loss through elastic stress waves occurs due to jump-to-contact and jump-out-of-contact instabilities and can even result in capture of the elastic sphere on the half-space.

Spatial variations of effective elastic thickness over the Ninetyeast Ridge and implications for its structure and tectonic evolution

We present new data on the strength of oceanic lithosphere along the Ninetyeast Ridge (NER) from two independent methods: spectral analysis (Bouguer coherence) using the fan wavelet transform technique, and spatial analysis (flexure inversion) with the convolution method. The two methods provide effective elastic thickness (T-e) patterns that broadly complement each other, and correlate well with known surface structures and regional-scale features. Furthermore, our study presents a new high resolution database on the Moho configuration, which obeys flexural isostasy, and exhibit regional correlations with the T-e variations. A continuous ridge structure with a much lower T-e value than that of normal oceanic lithosphere provides strong support for the hotspot theory. The derived T-e values vary over the northern (higher T-e similar to 10-20 km), central (anomalously low T-e similar to 0-5 km), and southern (low T-e similar to 5 km) segments of the NER. The lack of correlation of the T-e value with the progressive aging of the lithosphere implies differences in thermo-mechanical setting of the crust and underlying mantle in different parts of the NER, again indicating diversity in their evolution. The anomalously low T-e and deeper Moho (similar to 22 km) estimates of the central NER (between 0.5 degrees N and 17 degrees S) are attributed to the interaction of a hotspot with the Wharton spreading ridge that caused significant thermal rejuvenation and hence weakening of the lithosphere. The higher mechanical strength values in the northern NER (north of 0.5 degrees N) may support the idea of off-ridge emplacement and a relatively large plate motion at the time of volcanism. The low T-e and deeper Moho (similar to 22 km) estimates in the southern part (south of 17 degrees S) suggest that the lithosphere was weak and therefore younger at the time of volcanism, and this supports the idea that the southern NER was emplaced on the edge of the Indian plate. (C) 2013 Elsevier B.V. All rights reserved.

Two-pion low-energy contribution to the muon g-2 with improved precision from analyticity and unitarity

The two-pion contribution from low energies to the muon magnetic moment anomaly, although small, has a large relative uncertainty since in this region the experimental data on the cross sections are neither sufficient nor precise enough. It is therefore of interest to see whether the precision can be improved by means of additional theoretical information on the pion electromagnetic form factor, which controls the leading-order contribution. In the present paper, we address this problem by exploiting analyticity and unitarity of the form factor in a parametrization-free approach that uses the phase in the elastic region, known with high precision from the Fermi-Watson theorem and Roy equations for pi pi elastic scattering as input. The formalism also includes experimental measurements on the modulus in the region 0.65-0.70 GeV, taken from the most recent e(+)e(-) ->pi(+)pi(-) experiments, and recent measurements of the form factor on the spacelike axis. By combining the results obtained with inputs from CMD2, SND, BABAR, and KLOE, we make the predictions a(mu)(pi pi,LO)2m(pi), 0.30 GeV] = (0.553 +/- 0.004) x 10(-10) and a(mu)(pi pi,LO)0.30 GeV; 0.63 GeV] = (133.083 +/- 0.837) x 10(-10). These are consistent with the other recent determinations and have slightly smaller errors.

Energy loss due to adhesion in longitudinal impact of elastic cylinders

Adhesive interaction between impacting bodies can cause energy loss, even in an otherwise elastic impact. Adhesion force induces tensile stress in the bodies, which modifies the stress wave profile and influences the restitution behavior. We investigate this effect by developing a finite element framework, which incorporates a Lennard-Jones-type potential for modeling the adhesive interaction between volume elements. With this framework, the classical problems in contact mechanics can be revisited without the restrictive surface-force approximation. In this paper, we study the longitudinal impact of an elastic cylinder on a rigid half-space with adhesion. In the absence of adhesion, this problem reduces to the impact between two identical cylinders in which there is no energy loss. Adhesion causes a fraction of energy in the stress waves to remain in the cylinder as residual stress waves. This apparent loss in kinetic energy is shown to be a unique function of maximum tensile strain energy. We have developed a 1-D model in terms of interaction force parameters, velocity and material properties to estimate the tensile stain energy. We show that this model can be used to predict practically important phenomena like capture wherein the impacting bodies stick together. (C) 2013 Elsevier Masson SAS. All rights reserved.

Tectonic inheritance of the Indian Shield: New insights from its elastic thickness structure

A new evaluation of the elastic thickness (Te) structure of the Indian Shield, derived from isotropic fan wavelet methodology, documents spatial variations of lithospheric deformation in different tectonic provinces correlated with episodic tectono-thermal events. The Te variations corroborated by shear velocity, crustal thickness, and seismogenic thickness reveal the heterogeneous rheology of the Indian lithosphere. The thinned, attenuated lithosphere beneath Peninsular India is considered to be the reason for its mechanically weak strength (<30 km), where a decoupled crust-mantle rheology under different surface/subsurface loading structures may explain the prominent low Te patterns. The arcuate Te structure of the Western Dharwar province and a NNE-trending band of low Te anomaly in the Southern Granulite Terrane are intriguing patterns. The average Te values (40-50 km) of the Central Indian Tectonic Zone, the Bastar Craton, and the northern Eastern Ghats Mobile Belt are suggestive of old, stable, Indian lithosphere, which was not affected by any major tectono-thermal events after cratonic stabilization. We propose that the anomalously high Te (60-85 km) and high S-wave velocity zone to the north of the Narmada-Son Lineament, mainly in NW Himalaya, and the northern Aravalli and Bundelkhand Cratons, suggest that Archean lithosphere characterized by a high velocity mantle keel supports the orogenic topographic loads in/near the Himalaya. The Te map clearly segments the volcanic provinces of the Indian Shield, where the signatures of the Reunion, Marion, and Kerguelen hotspots are indicated by significantly low Te patterns that correlate with plume- and rift-related thermal and mechanical rejuvenation, magmatic underplating, and crustal necking. The correlations between Te variations and the occurrence of seismicity over seismically active zones reveal different causal relationships, which led to the current seismogenic zonation of the Indian Shield. (C) 2013 Elsevier B.V. All rights reserved.

Large-Scale Elastic Net Regularized Linear Classification SVMs and Logistic Regression

Elastic Net Regularizers have shown much promise in designing sparse classifiers for linear classification. In this work, we propose an alternating optimization approach to solve the dual problems of elastic net regularized linear classification Support Vector Machines (SVMs) and logistic regression (LR). One of the sub-problems turns out to be a simple projection. The other sub-problem can be solved using dual coordinate descent methods developed for non-sparse L2-regularized linear SVMs and LR, without altering their iteration complexity and convergence properties. Experiments on very large datasets indicate that the proposed dual coordinate descent - projection (DCD-P) methods are fast and achieve comparable generalization performance after the first pass through the data, with extremely sparse models.

Elastic Resources Framework in IaaS, preserving performance SLAs

Elasticity in cloud systems provides the flexibility to acquire and relinquish computing resources on demand. However, in current virtualized systems resource allocation is mostly static. Resources are allocated during VM instantiation and any change in workload leading to significant increase or decrease in resources is handled by VM migration. Hence, cloud users tend to characterize their workloads at a coarse grained level which potentially leads to under-utilized VM resources or under performing application. A more flexible and adaptive resource allocation mechanism would benefit variable workloads, such as those characterized by web servers. In this paper, we present an elastic resources framework for IaaS cloud layer that addresses this need. The framework provisions for application workload forecasting engine, that predicts at run-time the expected demand, which is input to the resource manager to modulate resource allocation based on the predicted demand. Based on the prediction errors, resources can be over-allocated or under-allocated as compared to the actual demand made by the application. Over-allocation leads to unused resources and under allocation could cause under performance. To strike a good trade-off between over-allocation and under-performance we derive an excess cost model. In this model excess resources allocated are captured as over-allocation cost and under-allocation is captured as a penalty cost for violating application service level agreement (SLA). Confidence interval for predicted workload is used to minimize this excess cost with minimal effect on SLA violations. An example case-study for an academic institute web server workload is presented. Using the confidence interval to minimize excess cost, we achieve significant reduction in resource allocation requirement while restricting application SLA violations to below 2-3%.

Diffusion in an elastic medium: A model for macromolecule transport across the nuclear pore complex

Nuclear pore complexes (NPCs) are very selective filters that sit on the membrane of the nucleus and monitor the transport between the cytoplasm and the nucleoplasm. For the central plug of NPC two models have been suggested in the literature. The first suggests that the plug is a reversible hydrogel while the other suggests that it is a polymer brush. Here we propose a model for the transport of a protein through the plug, which is general enough to cover both the models. The protein stretches the plug and creates a local deformation, which together with the protein, we refer to as the bubble. We start with the free energy for creation of the bubble and consider its motion within the plug. The relevant coordinate is the center of the bubble which executes random walk. We find that for faster relaxation of the gel, the diffusion of the bubble is greater. (C) 2014 Elsevier-B.V. All rights reserved.

First-Principles Study of Structure, Vibrational, and Elastic Properties of Stoichiometric and Calcium-Deficient Hydroxyapatite

Hydroxyapatite (HAp), a primary constituent of human bone, is usually nonstoichiometric with varying Ca/P molar ratios, with the well-known fact that Ca deficiency can cause marked reductions in its mechanical properties. To gain insights into the mechanism of this degradation, we employ first-principles calculations based on density functional theory and determine the effects of Ca deficiency on structure, vibrational, and elastic properties of HAp. Our simulation results confirm a considerable reduction in the elastic constants of HAp due to Ca deficiency, which was experimentally reported earlier. Stress-induced transformation of the Ca-deficient defected structure into a metastable state upon the application of stress could be a reason for this. Local structural stability of HAp and Ca-deficient HAp structures is assessed with full phonon dispersion studies. Further, specific signatures in the computed vibrational spectra for Ca deficiency in HAp can be utilized in experimental characterization of different types of defected HAp.