A fully discrete C-0 interior penalty Galerkin approximation of the extended Fisher-Kolmogorov equation

A fully discrete C-0 interior penalty finite element method is proposed and analyzed for the Extended Fisher-Kolmogorov (EFK) equation u(t) + gamma Delta(2)u - Delta u + u(3) - u = 0 with appropriate initial and boundary conditions, where gamma is a positive constant. We derive a regularity estimate for the solution u of the EFK equation that is explicit in gamma and as a consequence we derive a priori error estimates that are robust in gamma. (C) 2013 Elsevier B.V. All rights reserved.

Reconstructing the properties of dark energy using standard sirens

Future space-based gravity wave (GW) experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of GW sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source in only 68% of the cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini et al. T.D. Saini, S.K. Sethi, and V. Sahni, Phys. Rev. D 81, 103009 (2010)] argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as type Ia supernovae, baryon acoustic oscillations, cosmic microwave background, etc. DOI:10.1103/PhysRevD.87.083001

Chapman-Enskog analysis for finite-volume formulation of lattice Boltzmann equation

The classical Chapman-Enskog expansion is performed for the recently proposed finite-volume formulation of lattice Boltzmann equation (LBE) method D.V. Patil, K.N. Lakshmisha, Finite volume TVD formulation of lattice Boltzmann simulation on unstructured mesh, J. Comput. Phys. 228 (2009) 5262-5279]. First, a modified partial differential equation is derived from a numerical approximation of the discrete Boltzmann equation. Then, the multi-scale, small parameter expansion is followed to recover the continuity and the Navier-Stokes (NS) equations with additional error terms. The expression for apparent value of the kinematic viscosity is derived for finite-volume formulation under certain assumptions. The attenuation of a shear wave, Taylor-Green vortex flow and driven channel flow are studied to analyze the apparent viscosity relation.

Rheometry of dense granular materials: the crucial effects of gravity and confining walls

We describe here the rheological response of dense, slowly deforming granular materials to shear in a cylindrical Couette cell. All components of the stress on the outer cylinder are measured pointwise as a function of the depth, for different methods of construction of the bed that presumably lead to distinct fabrics. The static stress profiles for the different construction protocols are different, but a stress profile that is independent of construction history emerges when the granular column is sheared for sufficient time, in accord with the predictions of plasticity theories. However the qualitative features of the the stress profile under shear differs radically from the predictions of plasticity theories and data reported in earlier studies. We discuss a hypothesis for the anomalous stress profiles that was proposed recently by us, and the ways in which further experiments may to conducted to verify it.

Rheometry of granular materials in cylindrical Couette cells: Anomalous stress caused by gravity and shear

The cylindrical Couette device is commonly employed to study the rheology of fluids, but seldom used for dense granular materials. Plasticity theories used for granular flows predict a stress field that is independent of the shear rate, but otherwise similar to that in fluids. In this paper we report detailed measurements of the stress as a function of depth, and show that the stress profile differs fundamentally from that of fluids, from the predictions of plasticity theories, and from intuitive expectation. In the static state, a part of the weight of the material is transferred to the walls by a downward vertical shear stress, bringing about the well-known Janssen saturation of the stress in vertical columns. When the material is sheared, the vertical shear stress changes sign, and the magnitudes of all components of the stress rise rapidly with depth. These qualitative features are preserved over a range of the Couette gap and shear rate, for smooth and rough walls and two model granular materials. To explain the anomalous rheological response, we consider some hypotheses that seem plausibleapriori, but showthat none survive after careful analysis of the experimental observations. We argue that the anomalous stress is due to an anisotropic fabric caused by the combined actions of gravity, shear, and frictional walls, for which we present indirect evidence from our experiments. A general theoretical framework for anisotropic plasticity is then presented. The detailed mechanics of how an anisotropic fabric is brought about by the above-mentioned factors is not clear, and promises to be a challenging problem for future investigations. (C) 2013 AIP Publishing LLC.

Solution of Time Dependent Joule Heat Equation for a Graphene Sheet Under Thomson Effect

We address a physics-based solution of joule heating phenomenon in a single-layer graphene (SLG) sheet under the presence of Thomson effect. We demonstrate that the temperature in an isotopically pure (containing only C-12) SLG sheet attains its saturation level quicker than when doped with its isotopes (C-13). From the solution of the joule heating equation, we find that the thermal time constant of the SLG sheet is in the order of tenths of a nanosecond for SLG dimensions of a few micrometers. These results have been formulated using the electron interactions with the inplane and flexural phonons to demonstrate a field-dependent Landauer transmission coefficient. We further develop an analytical model of the SLG specific heat using the quadratic (out of plane) phonon band structure over the room temperature. Additionally, we show that a cooling effect in the SLG sheet can be substantially enhanced with the addition of C-13. The methodologies as discussed in this paper can be put forward to analyze the graphene heat spreader theory.

Ground motion prediction equation considering combined dataset of recorded and simulated ground motions

Himalayan region is one of the most active seismic regions in the world and many researchers have highlighted the possibility of great seismic event in the near future due to seismic gap. Seismic hazard analysis and microzonation of highly populated places in the region are mandatory in a regional scale. Region specific Ground Motion Predictive Equation (GMPE) is an important input in the seismic hazard analysis for macro- and micro-zonation studies. Few GMPEs developed in India are based on the recorded data and are applicable for a particular range of magnitudes and distances. This paper focuses on the development of a new GMPE for the Himalayan region considering both the recorded and simulated earthquakes of moment magnitude 5.3-8.7. The Finite Fault simulation model has been used for the ground motion simulation considering region specific seismotectonic parameters from the past earthquakes and source models. Simulated acceleration time histories and response spectra are compared with available records. In the absence of a large number of recorded data, simulations have been performed at unavailable locations by adopting Apparent Stations concept. Earthquakes recorded up to 2007 have been used for the development of new GMPE and earthquakes records after 2007 are used to validate new GMPE. Proposed GMPE matched very well with recorded data and also with other highly ranked GMPEs developed elsewhere and applicable for the region. Comparison of response spectra also have shown good agreement with recorded earthquake data. Quantitative analysis of residuals for the proposed GMPE and region specific GMPEs to predict Nepal-India 2011 earthquake of Mw of 5.7 records values shows that the proposed GMPE predicts Peak ground acceleration and spectral acceleration for entire distance and period range with lower percent residual when compared to exiting region specific GMPEs. Crown Copyright (C) 2013 Published by Elsevier Ltd. All rights reserved.

Entanglement entropy from the holographic stress tensor

We consider entanglement entropy in the context of gauge/gravity duality for conformal field theories in even dimensions. The holographic prescription due to Ryu and Takayanagi (RT) leads to an equation describing how the entangling surface extends into the bulk geometry. We show that setting to zero, the timetime component of the Brown-York stress tensor evaluated on the co-dimension 1 entangling surface, leads to the same equation. By considering a spherical entangling surface as an example, we observe that the Euclidean actionmethods in AdS/CFT will lead to the RT area functional arising as a counterterm needed to regularize the stress tensor. We present arguments leading to a justification for the minimal area prescription.

Turbulence in the two-dimensional Fourier-truncated Gross-Pitaevskii equation

We undertake a systematic, direct numerical simulation of the twodimensional, Fourier-truncated, Gross-Pitaevskii equation to study the turbulent evolutions of its solutions for a variety of initial conditions and a wide range of parameters. We find that the time evolution of this system can be classified into four regimes with qualitatively different statistical properties. Firstly, there are transients that depend on the initial conditions. In the second regime, powerlaw scaling regions, in the energy and the occupation-number spectra, appear and start to develop; the exponents of these power laws and the extents of the scaling regions change with time and depend on the initial condition. In the third regime, the spectra drop rapidly for modes with wave numbers k > kc and partial thermalization takes place for modes with k < kc; the self-truncation wave number kc(t) depends on the initial conditions and it grows either as a power of t or as log t. Finally, in the fourth regime, complete thermalization is achieved and, if we account for finite-size effects carefully, correlation functions and spectra are consistent with their nontrivial Berezinskii-Kosterlitz-Thouless forms. Our work is a natural generalization of recent studies of thermalization in the Euler and other hydrodynamical equations; it combines ideas from fluid dynamics and turbulence, on the one hand, and equilibrium and nonequilibrium statistical mechanics on the other.

Solutions of a system of forced Burgers equation

In this article, we obtain explicit solutions of a system of forced Burgers equation subject to some classes of bounded and compactly supported initial data and also subject to certain unbounded initial data. In a series of papers, Rao and Yadav (2010) 1-3] obtained explicit solutions of a nonhomogeneous Burgers equation in one dimension subject to certain classes of bounded and unbounded initial data. Earlier Kloosterziel (1990) 4] represented the solution of an initial value problem for the heat equation, with initial data in L-2 (R-n, e(vertical bar x vertical bar 2/2)), as a series of self-similar solutions of the heat equation in R-n. Here we express the solutions of certain classes of Cauchy problems for a system of forced Burgers equation in terms of self-similar solutions of some linear partial differential equations. (C) 2013 Elsevier Inc. All rights reserved.

On generalized gravitational entropy, squashed cones and holography

We consider generalized gravitational entropy in various higher derivative theories of gravity dual to four dimensional CFTs using the recently proposed regularization of squashed cones. We derive the universal terms in the entanglement entropy for spherical and cylindrical surfaces. This is achieved by constructing the Fefferman-Graham expansion for the leading order metrics for the bulk geometry and evaluating the generalized gravitational entropy. We further show that the Wald entropy evaluated in the bulk geometry constructed for the regularized squashed cones leads to the correct universal parts of the entanglement entropy for both spherical and cylindrical entangling surfaces. We comment on the relation with the Iyer-Wald formula for dynamical horizons relating entropy to a Noether charge. Finally we show how to derive the entangling surface equation in Gauss-Bonnet holography.

A Constitutive Equation for Grain Boundary Sliding: An Experimental Approach

Although grain boundary sliding (GBS) has been recognized as an important process during high-temperature deformation in crystalline materials, there is paucity in experimental data for characterizing a constitutive equation for GBS. High-temperature tensile creep experiments were conducted, together with measurements of GBS at different strains, stresses, grain sizes, and temperatures. Experimental data obtained on a Mg AZ31 alloy demonstrate that, for the first time, dynamic recrystallization during creep does not alter the contribution of GBS to creep during high-temperature deformation. The experimentally observed invariance of the sliding contribution with strain was used together with the creep data for developing a constitutive equation for GBS in a manner similar to the standard creep equation. Using this new approach, it is demonstrated that the stress, grain size, and temperature dependence for creep and GBS are identical. This is rationalized by a model based on GBS controlled by dislocations, within grains or near-grain boundaries. (C) The Minerals, Metals & Materials Society and ASM International 2013

Nonlinear gravity-wave interactions in stratified turbulence

To investigate the dynamics of gravity waves in stratified Boussinesq flows, a model is derived that consists of all three-gravity-wave-mode interactions (the GGG model), excluding interactions involving the vortical mode. The GGG model is a natural extension of weak turbulence theory that accounts for exact three-gravity-wave resonances. The model is examined numerically by means of random, large-scale, high-frequency forcing. An immediate observation is a robust growth of the so-called vertically sheared horizontal flow (VSHF). In addition, there is a forward transfer of energy and equilibration of the nonzero-frequency (sometimes called ``fast'') gravity-wave modes. These results show that gravity-wave-mode interactions by themselves are capable of systematic interscale energy transfer in a stratified fluid. Comparing numerical simulations of the GGG model and the full Boussinesq system, for the range of Froude numbers (Fr) considered (0.05 a parts per thousand currency sign Fr a parts per thousand currency sign 1), in both systems the VSHF is hardest to resolve. When adequately resolved, VSHF growth is more vigorous in the GGG model. Furthermore, a VSHF is observed to form in milder stratification scenarios in the GGG model than the full Boussinesq system. Finally, fully three-dimensional nonzero-frequency gravity-wave modes equilibrate in both systems and their scaling with vertical wavenumber follows similar power-laws. The slopes of the power-laws obtained depend on Fr and approach -2 (from above) at Fr = 0.05, which is the strongest stratification that can be properly resolved with our computational resources.

Constraining gravity using entanglement in AdS/CFT

We investigate constraints imposed by entanglement on gravity in the context of holography. First, by demanding that relative entropy is positive and using the Ryu-Takayanagi entropy functional, we find certain constraints at a nonlinear level for the dual gravity. Second, by considering Gauss-Bonnet gravity, we show that for a class of small perturbations around the vacuum state, the positivity of the two point function of the field theory stress tensor guarantees the positivity of the relative entropy. Further, if we impose that the entangling surface closes off smoothly in the bulk interior, we find restrictions on the coupling constant in Gauss-Bonnet gravity. We also give an example of an anisotropic excited state in an unstable phase with broken conformal invariance which leads to a negative relative entropy.

SODAR observations of inertia- gravity waves in the atmospheric boundary layer during the passage of tropical cyclone

This study reports characteristics of inertia-gravity waves (IGWs) in the atmospheric boundary layer during the passage of Tropical Cylone-03B, using the Doppler Sound Detection and Ranging (SODAR) observations at the Indian tropical station of Gadanki (13.45 degrees N, 79.2 degrees E; near the east coast of India). Wavelet analysis of horizontal winds indicates significant wave motion (60h) near the characteristic inertial period. The hodograph analysis of the filtered winds shows an anti-cyclonic turning of horizontal wind with height and time, indicating the presence of IGW. This study finds important implications in boundary layer dynamics during the passage of tropical cyclones.