Elasto-plastic analysis of jointed rocks using discrete continuum and equivalent continuum approaches

Results from elasto-plastic numerical simulations of jointed rocks using both the equivalent continuum and discrete continuum approaches are presented, and are compared with experimental measurements. Initially triaxial compression tests on different types of rocks with wide variation in the uniaxial compressive strength are simulated using both the approaches and the results are compared. The applicability and relative merits and limitations of both the approaches for the simulation of jointed rocks are discussed. It is observed that both the approaches are reasonably good in predicting the real response. However, the equivalent continuum approach has predicted somewhat higher stiffness values at low strains. Considering the modelling effort involved in case of discrete continuum approach, for problems with complex geometry, it is suggested that a proper equivalent continuum model can be used, without compromising much on the accuracy of the results. Then the numerical analysis of a tunnel in Japan is taken up using the continuum approach. The deformations predicted are compared well against the field measurements and the predictions from discontinuum analysis. (C) 2012 Elsevier Ltd. All rights reserved.

Bose-Hubbard models in confining potentials: Inhomogeneous mean-field theory

We present an extensive study of Mott insulator (MI) and superfluid (SF) shells in Bose-Hubbard (BH) models for bosons in optical lattices with harmonic traps. For this we apply the inhomogeneous mean-field theory developed by Sheshadri et al. Phys. Rev. Lett. 75, 4075 (1995)]. Our results for the BH model with one type of spinless bosons agree quantitatively with quantum Monte Carlo simulations. Our approach is numerically less intensive than such simulations, so we are able to perform calculations on experimentally realistic, large three-dimensional systems, explore a wide range of parameter values, and make direct contact with a variety of experimental measurements. We also extend our inhomogeneous mean-field theory to study BH models with harmonic traps and (a) two species of bosons or (b) spin-1 bosons. With two species of bosons, we obtain rich phase diagrams with a variety of SF and MI phases and associated shells when we include a quadratic confining potential. For the spin-1 BH model, we show, in a representative case, that the system can display alternating shells of polar SF and MI phases, and we make interesting predictions for experiments in such systems.

Viral Kinetics Suggests a Reconciliation of the Disparate Observations of the Modulation of Claudin-1 Expression on Cells Exposed to Hepatitis C Virus

The tight junction protein claudin-1 (CLDN1) is necessary for hepatitis C virus (HCV) entry into target cells. Recent studies have made disparate observations of the modulation of the expression of CLDN1 on cells following infection by HCV. In one study, the mean CLDN1 expression on cells exposed to HCV declined, whereas in another study HCV infected cells showed increased CLDN1 expression compared to uninfected cells. Consequently, the role of HCV in modulating CLDN1 expression, and hence the frequency of cellular superinfection, remains unclear. Here, we present a possible reconciliation of these disparate observations. We hypothesized that viral kinetics and not necessarily HCV-induced receptor modulation underlies these disparate observations. To test this hypothesis, we constructed a mathematical model of viral kinetics in vitro that mimicked the above experiments. Model predictions provided good fits to the observed evolution of the distribution of CLDN1 expression on cells following exposure to HCV. Cells with higher CLDN1 expression were preferentially infected and outgrown by cells with lower CLDN1 expression, resulting in a decline of the mean CLDN1 expression with time. At the same time, because the susceptibility of cells to infection increased with CLDN1 expression, infected cells tended to have higher CLDN1 expression on average than uninfected cells. Our study thus presents an explanation of the disparate observations of CLDN1 expression following HCV infection and points to the importance of considering viral kinetics in future studies of receptor expression on cells exposed to HCV.

Constraining form factors with the method of unitarity bounds

The availability of a reliable bound on an integral involving the square of the modulus of a form factor on the unitarity cut allows one to constrain the form factor at points inside the analyticity domain and its shape parameters, and also to isolate domains on the real axis and in the complex energy plane where zeros are excluded. In this lecture note, we review the mathematical techniques of this formalism in its standard form, known as the method of unitarity bounds, and recent developments which allow us to include information on the phase and modulus along a part of the unitarity cut. We also provide a brief summary of some results that we have obtained in the recent past, which demonstrate the usefulness of the method for precision predictions on the form factors.

State Dependent Attempt Rate modeling of single cell IEEE 802.11 WLANs with homogeneous nodes and Poisson packet arrivals

We study a State Dependent Attempt Rate (SDAR) approximation to model M queues (one queue per node) served by the Carrier Sense Multiple Access with Collision Avoidance (CSMA/CA) protocol as standardized in the IEEE 802.11 Distributed Coordination Function (DCF). The approximation is that, when n of the M queues are non-empty, the (transmission) attempt probability of each of the n non-empty nodes is given by the long-term (transmission) attempt probability of n saturated nodes. With the arrival of packets into the M queues according to independent Poisson processes, the SDAR approximation reduces a single cell with non-saturated nodes to a Markovian coupled queueing system. We provide a sufficient condition under which the joint queue length Markov chain is positive recurrent. For the symmetric case of equal arrival rates and finite and equal buffers, we develop an iterative method which leads to accurate predictions for important performance measures such as collision probability, throughput and mean packet delay. We replace the MAC layer with the SDAR model of contention by modifying the NS-2 source code pertaining to the MAC layer, keeping all other layers unchanged. By this model-based simulation technique at the MAC layer, we achieve speed-ups (w.r.t. MAC layer operations) up to 5.4. Through extensive model-based simulations and numerical results, we show that the SDAR model is an accurate model for the DCF MAC protocol in single cells. (C) 2012 Elsevier B.V. All rights reserved.

A drop of active matter

We study theoretically the hydrodynamics of a fluid drop containing oriented filaments endowed with active contractile or extensile stresses and placed on a solid surface. The active stresses alter qualitatively the wetting properties of the drop, leading to new spreading laws and novel static drop shapes. Candidate systems for testing our predictions include cytoskeletal extracts with motors and ATP, suspensions of bacteria or pulsatile cells, or fluids laden with artificial self-propelled colloids.

High-pressure synchrotron x-ray diffraction study of RMnO3 (R=Eu, Gd, Tb and Dy) upto 50 GPa

We have carried out synchrotron based high-pressure x-ray diffraction study of orthorhombic EuMnO3, GdMnO3, TbMnO3 and DyMnO3 up to 54.4, 41.6, 47.0 and 50.2 GPa, respectively. The diffraction peaks of all the four manganites shift monotonically to higher diffraction angles and the crystals retain the orthorhombic structure till the highest pressure. We have fitted the observed volume versus pressure data with the Birch-Murnaghan equation of state and determined the bulk modulus to be 185 +/- 6 GPa, 190 +/- 16 GPa, 188 +/- 9 GPa and 192 +/- 8 GPa for EuMnO3, GdMnO3, TbMnO3 and DyMnO3, respectively. The bulk modulus of EuMnO3 is comparable to other manganites, in contrast to theoretical predictions.

Anomalous Stress Profile in a Sheared Granular Column

We present measurements of the stress as a function of vertical position in a column of granular material sheared in a cylindrical Couette device. All three components of the stress tensor on the outer cylinder were measured as a function of distance from the free surface at shear rates low enough that the material was in the dense, slow flow regime. We find that the stress profile differs fundamentally from that of fluids, from the predictions of plasticity theories, and from intuitive expectation. We argue that the anomalous stress profile is due to an anisotropic fabric caused by the combined action of gravity and shear.

A finite population model of mlecular evolution: theory and computation

This article is concerned with the evolution of haploid organisms that reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed the quasispecies model in an attempt to understand such an evolutionary process. Their work has impacted antiviral treatment and vaccine design strategies. Yet, predictions of the quasispecies model are at best viewed as a guideline, primarily because it assumes an infinite population size, whereas realistic population sizes can be quite small. In this paper we consider a population genetics-based model aimed at understanding the evolution of such organisms with finite population sizes and present a rigorous study of the convergence and computational issues that arise therein. Our first result is structural and shows that, at any time during the evolution, as the population size tends to infinity, the distribution of genomes predicted by our model converges to that predicted by the quasispecies model. This justifies the continued use of the quasispecies model to derive guidelines for intervention. While the stationary state in the quasispecies model is readily obtained, due to the explosion of the state space in our model, exact computations are prohibitive. Our second set of results are computational in nature and address this issue. We derive conditions on the parameters of evolution under which our stochastic model mixes rapidly. Further, for a class of widely used fitness landscapes we give a fast deterministic algorithm which computes the stationary distribution of our model. These computational tools are expected to serve as a framework for the modeling of strategies for the deployment of mutagenic drugs.

Premixed flame response to equivalence ratio fluctuations: Comparison between reduced order modeling and detailed computations

Combustion instability events in lean premixed combustion systems can cause spatio-temporal variations in unburnt mixture fuel/air ratio. This provides a driving mechanism for heat-release oscillations when they interact with the flame. Several Reduced Order Modelling (ROM) approaches to predict the characteristics of these oscillations have been developed in the past. The present paper compares results for flame describing function characteristics determined from a ROM approach based on the level-set method, with corresponding results from detailed, fully compressible reacting flow computations for the same two dimensional slot flame configuration. The comparison between these results is seen to be sensitive to small geometric differences in the shape of the nominally steady flame used in the two computations. When the results are corrected to account for these differences, describing function magnitudes are well predicted for frequencies lesser than and greater than a lower and upper cutoff respectively due to amplification of flame surface wrinkling by the convective Darrieus-Landau (DL) instability. However, good agreement in describing function phase predictions is seen as the ROM captures the transit time of wrinkles through the flame correctly. Also, good agreement is seen for both magnitude and phase of the flame response, for large forcing amplitudes, at frequencies where the DL instability has a minimal influence. Thus, the present ROM can predict flame response as long as the DL instability, caused by gas expansion at the flame front, does not significantly alter flame front perturbation amplitudes as they traverse the flame. (C) 2012 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Stability of developing film flow down an inclined surface

Film flows on inclined surfaces are often assumed to be of constant thickness, which ensures that the velocity profile is half-Poiseuille. It is shown here that by shallow water theory, only flows in a portion of Reynolds number-Froude number (Re-Fr) plane can asymptotically attain constant film thickness. In another portion on the plane, the constant thickness solution appears as an unstable fixed point, while in other regions the film thickness seems to asymptote to a positive slope. Our simulations of the Navier-Stokes equations confirm the predictions of shallow water theory at higher Froude numbers, but disagree with them at lower Froude numbers. We show that different regimes of film flow show completely different stability behaviour from that predicted earlier. Supercritical decelerating flows are shown to be always unstable, whereas accelerating flows become unstable below a certain Reynolds number for a given Froude number. Subcritical flows on the other hand are shown to be unstable above a certain Reynolds number. In some range of parameters, two solutions for the base flowexist, and the attached profile is found to be more stable. All flows except those with separation become more stable as they proceed downstream. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4758299]

TURBULENT PUMPING OF MAGNETIC FLUX REDUCES SOLAR CYCLE MEMORY AND THUS IMPACTS PREDICTABILITY OF THE SUN'S ACTIVITY

Prediction of the Sun's magnetic activity is important because of its effect on space environment and climate. However, recent efforts to predict the amplitude of the solar cycle have resulted in diverging forecasts with no consensus. Yeates et al. have shown that the dynamical memory of the solar dynamo mechanism governs predictability, and this memory is different for advection- and diffusion-dominated solar convection zones. By utilizing stochastically forced, kinematic dynamo simulations, we demonstrate that the inclusion of downward turbulent pumping of magnetic flux reduces the memory of both advection- and diffusion-dominated solar dynamos to only one cycle; stronger pumping degrades this memory further. Thus, our results reconcile the diverging dynamo-model-based forecasts for the amplitude of solar cycle 24. We conclude that reliable predictions for the maximum of solar activity can be made only at the preceding minimum-allowing about five years of advance planning for space weather. For more accurate predictions, sequential data assimilation would be necessary in forecasting models to account for the Sun's short memory.

Fluctuations and correlations of conserved charges in the (2+1) Polyakov quark meson model

We consider the (2 + 1) flavor Polyakov quark meson model and study the fluctuations (correlations) of conserved charges up to sixth (fourth) order. A comparison is made with lattice data wherever available and overall good qualitative agreement is found, more so for the case of the normalized susceptibilities. The model predictions for the ratio of susceptibilities go to that of an ideal gas of hadrons as in hadron resonance gas model at low temperatures while at high temperature the values are close to that of an ideal gas of massless quarks. Our study provides a strong basis for the use of the Polyakov quark meson model as an effective model to understand the topology of the QCD phase diagram. DOI: 10.1103/PhysRevD.86.114021 PACS numbers: 12.39.-x, 05.40.-a, 12.38.Aw, 12.38.Mh

Subdiffusion in hair bundle dynamics: The role of protein conformational fluctuations

The detection of sound signals in vertebrates involves a complex network of different mechano-sensory elements in the inner ear. An especially important element in this network is the hair bundle, an antenna-like array of stereocilia containing gated ion channels that operate under the control of one or more adaptation motors. Deflections of the hair bundle by sound vibrations or thermal fluctuations transiently open the ion channels, allowing the flow of ions through them, and producing an electrical signal in the process, eventually causing the sensation of hearing. Recent high frequency (0.1-10 kHz) measurements by Kozlov et al. Proc. Natl. Acad. Sci. U. S. A. 109, 2896 (2012)] of the power spectrum and the mean square displacement of the thermal fluctuations of the hair bundle suggest that in this regime the dynamics of the hair bundle are subdiffusive. This finding has been explained in terms of the simple Brownian motion of a filament connecting neighboring stereocilia (the tip link), which is modeled as a viscoelastic spring. In the present paper, the diffusive anomalies of the hair bundle are ascribed to tip link fluctuations that evolve by fractional Brownian motion, which originates in fractional Gaussian noise and is characterized by a power law memory. The predictions of this model for the power spectrum of the hair bundle and its mean square displacement are consistent with the experimental data and the known properties of the tip link. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4768902]

Microstructural Evolution and Some Unusual Effects During Thermo-Mechanical Cycling of Sn-Ag-Cu Alloys

Sn-Ag-Cu (SAC) solders are susceptible to appreciable microstructural coarsening during storage or service. This results in evolution of joint properties over time, and thereby influences the long-term reliability of microelectronic packages. Accurate prediction of this aging behavior is therefore critical for joint reliability predictions. Here, we study the precipitate coarsening behavior in two Sn-Ag-Cu (SAC) alloys, namely Sn-3.0Ag-0.5Cu and Sn-1.0Cu-0.5Cu, under different thermo-mechanical excursions, including isothermal aging at 150 degrees C for various lengths of time and thermo-mechanical cycling between -25 degrees C and 125 degrees C, with an imposed shear strain of similar to 19.6% per cycle, for different number of cycles. During isothermal aging and the thermo-mechanical cycling up to 200 cycles, Ag3Sn precipitates undergo rapid, monotonous coarsening. However, high number of thermo-mechanical cycling, usually between 200 and 600 cycles, causes dissolution and re-precipitation of precipitates, resulting in a fine and even distribution. Also, recrystallization of Sn-grains near precipitate clusters was observed during severe isothermal aging. Such responses are quite unusual for SAC solder alloys. In the regime of usual precipitate coarsening in these SAC alloys, an explicit parameter, which captures the thermo-mechanical history dependence of Ag3Sn particle size, was defined. Brief mechanistic description for the recrystallization of Sn grains during isothermal aging and reprecipitation of the Ag3Sn due to high number of thermo-mechanical cycles are also presented.