Boundary perturbations and the Helmholtz equation in three dimensions

We propose an analytic perturbative scheme in the spirit of Lord Rayleigh's work for determining the eigenvalues of the Helmholtz equation in three dimensions inside an arbitrary boundary where the eigenfunction satisfies either the Dirichlet boundary condition or the Neumann boundary condition. Although numerous works are available in the literature for arbitrary boundaries in two dimensions, to the best of our knowledge the formulation in three dimensions is proposed for the first time. In this novel prescription, we have expanded the arbitrary boundary in terms of spherical harmonics about an equivalent sphere and obtained perturbative closed-form solutions at each order for the problem in terms of corrections to the equivalent spherical boundary for both the boundary conditions. This formulation is in parallel with the standard time-independent Rayleigh-Schrodinger perturbation theory. The efficacy of the method is tested by comparing the perturbative values against the numerically calculated eigenvalues for spheroidal, superegg and superquadric shaped boundaries. It is shown that this perturbation works quite well even for wide departure from spherical shape and for higher excited states too. We believe this formulation would find applications in the field of quantum dots and acoustical cavities.

Suppression of instabilities in burning droplets using preferential acoustic perturbations

Self-induced internal boiling in burning functional droplets has been observed to induce severe bulk shape oscillations in droplets with characteristic bubble ejection events that corrugate the droplet surface. Such bubble-droplet interactions are characterized by a distinct regime of a single bubble growing inside the droplet where evaporative Darrieus-Landau instability occurs at the bubble-droplet interface. In this regime the bubble-droplet system behaves as a self-excited coupled oscillator. In this study, we report the external flame-acoustic interaction with bubbles inside the droplet resulting in controlled droplet deformation. In particular, by exciting the droplet flame in a critical, responsive frequency range (80 Hz <= f(p) <= 120 Hz) the droplet deformation cycle could be altered through suppression of these self-excited instabilities and intensity/frequency of bubble ejection events. This selective acoustic tuning also enabled the control of bubble dynamics, bulk droplet-shape distortion and the final precipitate morphology even in burning nanoparticle laden droplets. (C) 2014 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

NESTED SPARSE BAYESIAN LEARNING FOR BLOCK-SPARSE SIGNALS WITH INTRA-BLOCK CORRELATION

In this work, we address the recovery of block sparse vectors with intra-block correlation, i.e., the recovery of vectors in which the correlated nonzero entries are constrained to lie in a few clusters, from noisy underdetermined linear measurements. Among Bayesian sparse recovery techniques, the cluster Sparse Bayesian Learning (SBL) is an efficient tool for block-sparse vector recovery, with intra-block correlation. However, this technique uses a heuristic method to estimate the intra-block correlation. In this paper, we propose the Nested SBL (NSBL) algorithm, which we derive using a novel Bayesian formulation that facilitates the use of the monotonically convergent nested Expectation Maximization (EM) and a Kalman filtering based learning framework. Unlike the cluster-SBL algorithm, this formulation leads to closed-form EMupdates for estimating the correlation coefficient. We demonstrate the efficacy of the proposed NSBL algorithm using Monte Carlo simulations.

The cold mode: a phenomenological model for the evolution of density perturbations in the intracluster medium

Cool cluster cores are in global thermal equilibrium but are locally thermally unstable. We study a non-linear phenomenological model for the evolution of density perturbations in the intracluster medium (ICM) due to local thermal instability and gravity. We have analysed and extended a model for the evolution of an overdense blob in the ICM. We find two regimes in which the overdense blobs can cool to thermally stable low temperatures. One for large t(cool)/t(ff) (t(cool) is the cooling time and t(ff) is the free-fall time), where a large initial overdensity is required for thermal runaway to occur; this is the regime which was previously analysed in detail. We discover a second regime for t(cool)/t(ff) less than or similar to 1 (in agreement with Cartesian simulations of local thermal instability in an external gravitational field), where runaway cooling happens for arbitrarily small amplitudes. Numerical simulations have shown that cold gas condenses out more easily in a spherical geometry. We extend the analysis to include geometrical compression in weakly stratified atmospheres such as the ICM. With a single parameter, analogous to the mixing length, we are able to reproduce the results from numerical simulations; namely, small density perturbations lead to the condensation of extended cold filaments only if t(cool)/t(ff) less than or similar to 10.

Layered Exact-Repair Regenerating Codes via Embedded Error Correction and Block Designs

A new class of exact-repair regenerating codes is constructed by stitching together shorter erasure correction codes, where the stitching pattern can be viewed as block designs. The proposed codes have the help-by-transfer property where the helper nodes simply transfer part of the stored data directly, without performing any computation. This embedded error correction structure makes the decoding process straightforward, and in some cases the complexity is very low. We show that this construction is able to achieve performance better than space-sharing between the minimum storage regenerating codes and the minimum repair-bandwidth regenerating codes, and it is the first class of codes to achieve this performance. In fact, it is shown that the proposed construction can achieve a nontrivial point on the optimal functional-repair tradeoff, and it is asymptotically optimal at high rate, i.e., it asymptotically approaches the minimum storage and the minimum repair-bandwidth simultaneously.

Adaptive executions of hyperbolic block-structured AMR applications on GPU systems

A block-structured adaptive mesh refinement (AMR) technique has been used to obtain numerical solutions for many scientific applications. Some block-structured AMR approaches have focused on forming patches of non-uniform sizes where the size of a patch can be tuned to the geometry of a region of interest. In this paper, we develop strategies for adaptive execution of block-structured AMR applications on GPUs, for hyperbolic directionally split solvers. While effective hybrid execution strategies exist for applications with uniform patches, our work considers efficient execution of non-uniform patches with different workloads. Our techniques include bin-packing work units to load balance GPU computations, adaptive asynchronism between CPU and GPU executions using a knapsack formulation, and scheduling communications for multi-GPU executions. Our experiments with synthetic and real data, for single-GPU and multi-GPU executions, on Tesla S1070 and Fermi C2070 clusters, show that our strategies result in up to a 3.23 speedup in performance over existing strategies.

Seismic response of reduced-scale modular block and rigid faced reinforced walls through shaking table tests

This paper focuses on understanding the seismic response of geosynthetic reinforced retaining walls through shaking table tests on models of modular block and rigid faced reinforced retaining walls. Reduced-scale models of retaining walls reinforced with geogrid layers were constructed in a laminar box mounted on a uniaxial shaking table and subjected to various levels of sinusoidal base shaking. Models were instrumented with ultrasonic displacement sensors, earth pressure sensors and accelerometers. Effects of backfill density, number of reinforcement layers and reinforcement type on the performance of rigid faced and modular block walls were studied through different series of model tests. Performances of the walls were assessed in terms of face deformations, crest settlement and acceleration amplification at different elevations and compared. Modular block walls performed better than the rigid faced walls for the same level of base shaking because of the additional support derived by stacking the blocks with an offset. Type and quantity of reinforcement has significant effect on the seismic performance of both the types of walls. Displacements are more sensitive to relative density of the backfill and decrease with increasing relative density, the effect being more pronounced in case of unreinforced walls compared to the reinforced ones. Acceleration amplifications are not affected by the wall facing and inclusion of reinforcement. (C) 2015 Elsevier Ltd. All rights reserved.

Neoarchean arc magmatism followed by high-temperature, high-pressure metamorphism in the Nilgiri Block, southern India

The Nilgiri Block, southern India is an exhumed lower crust formed through arc magmatic processes in the Neoarchean. The main lithologies in this terrane include charnockites, gneisses, volcanic tuff, metasediments, banded iron formation and mafic-ultramafic bodies. Mafic-ultramafic rocks are present towards the northern and central part of the Nilgiri Block. We examine the evolution of these mafic granulites/metagabbros by phase diagram modeling and U-Pb sensitive high resolution ion microprobe (SHRIMP) dating. They consist of a garnet-clinopyroxene-plagioclase-hornblende-ilmenite +/- orthopyroxene +/- rutile assemblage. Garnet and clinopyroxene form major constituents with labradorite and orthopyroxene as the main mineral inclusions. Labradorite, identified using Raman analysis, shows typical peaks at 508 cm(-1), 479 cm(-1), 287 cm(-1) and 177 cm(-1). It is stable along with orthopyroxene towards the low-pressure high-temperature region of the granulite fades (M1 stage). Subsequently, orthopyroxene reacted with plagioclase to form the peak garnet + clinopyroxene + rutile assemblage (M2 stage). The final stage is represented by amphibolite facies-hornblende and plagioclase-rim around the garnet-clinopyroxene assemblage (M3 stage). Phase diagram modeling shows that these mafic granulites followed an anticlockwise P-T-t path during their evolution. The initial high-temperature metamorphism (M1 stage) was at 850-900 degrees C and similar to 9 kbar followed by high-pressure granulite fades metamorphism (M2 stage) at 850-900 degrees C and 14-15 kbar. U-Pb isotope studies of zircons using SHRIMP revealed late Neoarchean to early paleoproterozoic ages of crystallization and metamorphism respectively. The age data shows that these mafic granulites have undergone arc magmatism at ca. 25392 +/- 3 Ma and high-temperature, high-pressure metamorphism at ca. 2458.9 +/- 8.6 Ma. Thus our results suggests a late Neoarchean arc magmatism followed by early paleoproterozoic high-temperature, high-pressure granulite facies metamorphism due to the crustal thickening and suturing of the Nilgiri Block onto the Dharwar Craton. (C) 2015 Elsevier B.V. All rights reserved.

Archean tectonics and crustal evolution of the Biligiri Rangan Block, southern India

The Southern Granulite Terrain in India is a collage of crustal blocks ranging in age from Archean to Neoproterozoic. This study investigate the tectonic evolution of one of the northernmost block- the Biligiri Block (BRB) through a multidisciplinary approach involving field investigation, petrographic studies, LA-ICPMS zircon U-Pb geochronology, Hf isotopic analyses, metamorphic P-T phase diagram computations, and crustal thickness modeling. The garnet bearing quartzofeldspathic gneiss from the central BRB preserve Mesoarchean magmatic zircons with ages between 3207 and 2806 Ma and positive epsilon Hf value (+2.7) which possibly indicates vestiges of a Mesoarchean primitive continental crust. The occurrence of quartzite-iron formation intercalation as well as ultramafic lenses along the western boundary of the BRB is interpreted to indicate that the Kollegal structural lineament is a possible paleo-suture. Phase diagram computation of a metagabbro from the southwestern periphery of the Kollegal suture zone reveals high-pressure (similar to 18.5 kbar) and medium-temperature (similar to 840 degrees C) metamorphism, likely during eastward subduction of the Western Dharwar oceanic crust beneath the Mesoarchean BRB. In the model presented here, slab subduction, melting and underplating processes generated arc magmatism and subsequent charnockitization within the BRB between ca. 2650 Ma and ca. 2498 Ma. These results thus reveal Meso- to Neoarchean tectonic evolution of the BRB. The spatial variation of crustal thickness, derived from flexure inversion technique, provides additional constraints on the tectonic linkage of the BRB with its surrounding terrains. In conjunction with published data, the Moyar and the Kollegal suture zones are considered to mark the trace of ocean closure along which the Nilgiri and Biligiri Rangan Blocks accreted on to the Western Dharwar Craton. (C) 2016 Elsevier B.V. All rights reserved.

Stalagmite growth perturbations from the Kumaun Himalaya as potential earthquake recorders

The central part of the Himalaya (Kumaun and Garhwal Provinces of India) is noted for its prolonged seismic quiescence, and therefore, developing a longer-term time series of past earthquakes to understand their recurrence pattern in this segment assumes importance. In addition to direct observations of offsets in stratigraphic exposures or other proxies like paleoliquefaction, deformation preserved within stalagmites (speleothems) in karst system can be analyzed to obtain continuous millennial scale time series of earthquakes. The Central Indian Himalaya hosts natural caves between major active thrusts forming potential storehouses for paleoseismological records. Here, we present results from the limestone caves in the Kumaun Himalaya and discuss the implications of growth perturbations identified in the stalagmites as possible earthquake recorders. This article focuses on three stalagmites from the Dharamjali Cave located in the eastern Kumaun Himalaya, although two other caves, one of them located in the foothills, were also examined for their suitability. The growth anomalies in stalagmites include abrupt tilting or rotation of growth axes, growth termination, and breakage followed by regrowth. The U-Th age data from three specimens allow us to constrain the intervals of growth anomalies, and these were dated at 4273 +/- 410 years BP (2673-1853 BC), 2782 +/- 79 years BP (851-693 BC), 2498 +/- 117 years BP (605-371 BC), 1503 +/- 245 years BP (262-752 AD), 1346 +/- 101 years BP (563-765 AD), and 687 +/- 147 years BP (1176-1470 AD). The dates may correspond to the timings of major/great earthquakes in the region and the youngest event (1176-1470 AD) shows chronological correspondence with either one of the great medieval earthquakes (1050-1250 and 1259-1433 AD) evident from trench excavations across the Himalayan Frontal Thrust.

Effect of Rock Mass Structure and Block Size on the Slope Stability-Physical Modeling and Discrete Element Simulation

This paper studies the stability of jointed rock slopes by using our improved three-dimensional discrete element methods (DEM) and physical modeling. Results show that the DEM can simulate all failure modes of rock slopes with different joint configurations. The stress in each rock block is not homogeneous and blocks rotate in failure development. Failure modes depend on the configuration of joints. Toppling failure is observed for the slope with straight joints and sliding failure is observed for the slope with staged joints. The DEM results are also compared with those of limit equilibrium method (LEM). Without considering the joints in rock masses, the LEM predicts much higher factor of safety than physical modeling and DEM. The failure mode and factor of safety predicted by the DEM are in good agreement with laboratory tests for any jointed rock slope.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.