Analytical solutions for transient temperature distribution in a geothermal reservoir due to cold water injection

An analytical solution to describe the transient temperature distribution in a geothermal reservoir in response to injection of cold water is presented. The reservoir is composed of a confined aquifer, sandwiched between rocks of different thermo-geological properties. The heat transport processes considered are advection, longitudinal conduction in the geothermal aquifer, and the conductive heat transfer to the underlying and overlying rocks of different geological properties. The one-dimensional heat transfer equation has been solved using the Laplace transform with the assumption of constant density and thermal properties of both rock and fluid. Two simple solutions are derived afterwards, first neglecting the longitudinal conductive heat transport and then heat transport to confining rocks. Results show that heat loss to the confining rock layers plays a vital role in slowing down the cooling of the reservoir. The influence of some parameters, e.g. the volumetric injection rate, the longitudinal thermal conductivity and the porosity of the porous media, on the transient heat transport phenomenon is judged by observing the variation of the transient temperature distribution with different values of the parameters. The effects of injection rate and thermal conductivity have been found to be profound on the results.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

Modeling of Bubble-Size Distribution in Water and Freon Co-Blown Free Rise Polyurethane Foams

A model has been developed to simulate the foam characteristics obtained, when chemical (water) and physical (Freon) blowing agents are used together for the formation of polyurethane foams. The model considers the rate of reaction, the consequent rise in temperature of the reaction mixture, nucleation of bubbles, and mass transfer of CO2 and Freon to them till the time of gelation. The model is able to explain the experimental results available in literature. It further predicts that the nucleation period gets reduced with increase in water (at constant Freon content), whereas with increase in Freon (at constant water) concentration nucleation period decreases marginally leading to narrower bubble-size distribution. By the use of uniform sized nuclei added initially, the model predicts that the bubble-size distribution can be made independent of the rate of homogeneous nucleation and can, thus, offer an extra parameter for its control. (C) 2014 Wiley Periodicals, Inc.

On the possible variation of the lightning striking distance as assumed in the IEC lightning protection standard as a function of structure height

The effect of structure height on the lightning striking distance is estimated using a lightning strike model that takes into account the effect of connecting leaders. According to the results, the lightning striking distance may differ significantly from the values assumed in the IEC standard for structure heights beyond 30m. However, for structure heights smaller than about 30m, the results show that the values assumed by IEC do not differ significantly from the predictions based on a lightning attachment model taking into account the effect of connecting leaders. However, since IEC assumes a smaller striking distance than the ones predicted by the adopted model one can conclude that the safety is not compromised in adhering to the IEC standard. Results obtained from the model are also compared with Collection Volume Method (CVM) and other commonly used lightning attachment models available in the literature. The results show that in the case of CVM the calculated attractive distances are much larger than the ones obtained using the physically based lightning attachment models. This indicates the possibility of compromising the lightning protection procedures when using CVM. (C) 2014 Elsevier B.V. All rights reserved.

Fish distribution dynamics in the Aghanashini estuary of Uttara Kannada, west coast of India

Fish diversity (77 species) in the Aghanashini River estuary of the Indian west coast is linked to variable salinity conditions and zones I, II and III for high, medium and low salinity respectively. Zone I, the junction between Arabian Sea and the estuary, had all species in yearly succession due to freshwater conditions in monsoon to high salinity in pre-monsoon. The medium (zone II) and low (zone III) salinity mid and upstream portions had maximum of 67 and 39 fish species respectively. Maintenance of natural salinity regimes in estuary, among other ecological factors, is critical for its fish diversity.

Lantana camara L. (Verbenaceae) invasion along streams in a heterogeneous landscape

Streams are periodically disturbed due to flooding, act as edges between habitats and also facilitate the dispersal of propagules, thus being potentially more vulnerable to invasions than adjoining regions. We used a landscape-wide transect-based sampling strategy and a mixed effects modelling approach to understand the effects of distance from stream, a rainfall gradient, light availability and fire history on the distribution of the invasive shrub Lantana camara L.(lantana) in the tropical dry forests of Mudumalai in southern India. The area occupied by lantana thickets and lantana stem abundance were both found to be highest closest to streams across this landscape with a rainfall gradient. There was no advantage in terms of increased abundance or area occupied by lantana when it grew closer to streams in drier areas as compared to moister areas. On an average, the area covered by lantana increased with increasing annual rainfall. Areas that experienced greater number of fires during 1989-2010 had lower lantana stem abundance irrespective of distance from streams. In this landscape, total light availability did not affect lantana abundance. Understanding the spatially variable environmental factors in a heterogeneous landscape influencing the distribution of lantana would aid in making informed management decisions at this scale.

Give-and-take based peer-to-peer content distribution networks

Content Distribution Networks (CDNs) are widely used to distribute data to large number of users. Traditionally, content is being replicated among a number of surrogate servers, leading to high operational costs. In this context, Peer-to-Peer (P2P) CDNs have emerged as a viable alternative. An issue of concern in P2P networks is that of free riders, i.e., selfish peers who download files and leave without uploading anything in return. Free riding must be discouraged. In this paper, we propose a criterion, the Give-and-Take (G&T) criterion, that disallows free riders. Incorporating the G&T criterion in our model, we study a problem that arises naturally when a new peer enters the system: viz., the problem of downloading a `universe' of segments, scattered among other peers, at low cost. We analyse this hard problem, and characterize the optimal download cost under the G&T criterion. We propose an optimal algorithm, and provide a sub-optimal algorithm that is nearly optimal, but runs much more quickly; this provides an attractive balance between running time and performance. Finally, we compare the performance of our algorithms with that of a few existing P2P downloading strategies in use. We also study the computation time for prescribing the strategy for initial segment and peer selection for the newly arrived peer for various existing and proposed algorithms, and quantify cost-computation time trade-offs.

Dense shallow granular flows

Simplified equations are derived for a granular flow in the `dense' limit where the volume fraction is close to that for dynamical arrest, and the `shallow' limit where the stream-wise length for flow development (L) is large compared with the cross-stream height (h). The mass and diameter of the particles are set equal to 1 in the analysis without loss of generality. In the dense limit, the equations are simplified by taking advantage of the power-law divergence of the pair distribution function chi proportional to (phi(ad) - phi)(-alpha), and a faster divergence of the derivativ rho(d chi/d rho) similar to (d chi/d phi), where rho and phi are the density and volume fraction, and phi(ad) is the volume fraction for arrested dynamics. When the height h is much larger than the conduction length, the energy equation reduces to an algebraic balance between the rates of production and dissipation of energy, and the stress is proportional to the square of the strain rate (Bagnold law). In the shallow limit, the stress reduces to a simplified Bagnold stress, where all components of the stress are proportional to (partial derivative u(x)/partial derivative y)(2), which is the cross-stream (y) derivative of the stream-wise (x) velocity. In the simplified equations for dense shallow flows, the inertial terms are neglected in the y momentum equation in the shallow limit because the are O(h/L) smaller than the divergence of the stress. The resulting model contains two equations, a mass conservation equations which reduces to a solenoidal condition on the velocity in the incompressible limit, and a stream-wise momentum equation which contains just one parameter B which is a combination of the Bagnold coefficients and their derivatives with respect to volume fraction. The leading-order dense shallow flow equations, as well as the first correction due to density variations, are analysed for two representative flows. The first is the development from a plug flow to a fully developed Bagnold profile for the flow down an inclined plane. The analysis shows that the flow development length is ((rho) over barh(3)/B) , where (rho) over bar is the mean density, and this length is numerically estimated from previous simulation results. The second example is the development of the boundary layer at the base of the flow when a plug flow (with a slip condition at the base) encounters a rough base, in the limit where the momentum boundary layer thickness is small compared with the flow height. Analytical solutions can be found only when the stream-wise velocity far from the surface varies as x(F), where x is the stream-wise distance from the start of the rough base and F is an exponent. The boundary layer thickness increases as (l(2)x)(1/3) for all values of F, where the length scale l = root 2B/(rho) over bar. The analysis reveals important differences between granular flows and the flows of Newtonian fluids. The Reynolds number (ratio of inertial and viscous terms) turns out to depend only on the layer height and Bagnold coefficients, and is independent of the flow velocity, because both the inertial terms in the conservation equations and the divergence of the stress depend on the square of the velocity/velocity gradients. The compressibility number (ratio of the variation in volume fraction and mean volume fraction) is independent of the flow velocity and layer height, and depends only on the volume fraction and Bagnold coefficients.

Sensitivity of simulated climate to latitudinal distribution of solar insolation reduction in solar radiation management

Solar radiation management (SRM) geoengineering has been proposed as a potential option to counteract climate change. We perform a set of idealized geoengineering simulations using Community Atmosphere Model version 3.1 developed at the National Center for Atmospheric Research to investigate the global hydrological implications of varying the latitudinal distribution of solar insolation reduction in SRM methods. To reduce the solar insolation we have prescribed sulfate aerosols in the stratosphere. The radiative forcing in the geoengineering simulations is the net forcing from a doubling of CO2 and the prescribed stratospheric aerosols. We find that for a fixed total mass of sulfate aerosols (12.6 Mt of SO4), relative to a uniform distribution which nearly offsets changes in global mean temperature from a doubling of CO2, global mean radiative forcing is larger when aerosol concentration is maximum at the poles leading to a warmer global mean climate and consequently an intensified hydrological cycle. Opposite changes are simulated when aerosol concentration is maximized in the tropics. We obtain a range of 1 K in global mean temperature and 3% in precipitation changes by varying the distribution pattern in our simulations: this range is about 50% of the climate change from a doubling of CO2. Hence, our study demonstrates that a range of global mean climate states, determined by the global mean radiative forcing, are possible for a fixed total amount of aerosols but with differing latitudinal distribution. However, it is important to note that this is an idealized study and thus not all important realistic climate processes are modeled.

Flocking at a distance in active granular matter

The self-organized motion of vast numbers of creatures in a single direction is a spectacular example of emergent order. Here, we recreate this phenomenon using actuated nonliving components. We report here that millimetre-sized tapered rods, rendered motile by contact with an underlying vibrated surface and interacting through a medium of spherical beads, undergo a phase transition to a state of spontaneous alignment of velocities and orientations above a threshold bead area fraction. Guided by a detailed simulation model, we construct an analytical theory of this flocking transition, with two ingredients: a moving rod drags beads; neighbouring rods reorient in the resulting flow like a weathercock in the wind. Theory and experiment agree on the structure of our phase diagram in the plane of rod and bead concentrations and power-law spatial correlations near the phase boundary. Our discovery suggests possible new mechanisms for the collective transport of particulate or cellular matter.

Rapid divergence of ecotypes of an invasive plant

Invasive species demonstrate rapid evolution within a very short period of time allowing one to understand the underlying mechanism(s). Lantana camara, a highly invasive plant of the tropics and subtropics, has expanded its range and successfully established itself almost throughout India. In order to uncover the processes governing the invasion dynamics, 218 individuals from various locations across India were characterized with six microsatellites. By integrating genetic data with niche modelling, we examined the effect of drift and environmental selection on genetic divergence. We found multiple genetic clusters that were non-randomly distributed across space. Spatial autocorrelation revealed a strong fine-scale structure, i.e. isolation by distance. In addition, we obtained evidence of inhibitory effects of selection on gene flow, i.e. isolation by environmental distance. Perhaps, local adaptation in response to selection is offsetting gene flow and causing the populations to diverge. Niche models suggested that temperature and precipitation play a major role in the observed spatial distribution of this plant. Based on a non-random distribution of clusters, unequal gene flow among them and different bioclimatic niche requirements, we concluded that the emergence of ecotypes represented by two genetic clusters is underway. They may be locally adapted to specific climatic conditions, and perhaps at the very early stages of ecological divergence.

Distribution, Demography, and Conservation of Lion-tailed Macaques (Macaca silenus) in the Anamalai Hills Landscape, Western Ghats, India

The status of the endemic and endangered lion-tailed macaque (Macaca silenus) has not been properly assessed in several regions of the Western Ghats of southern India. We conducted a study in Parambikulam Forest Reserve in the state of Kerala to determine the distribution, demography, and status of lion-tailed macaques. We laid 5km(2) grid cells on the map of the study area (644km(2)) and made four replicated walks in each grid cell using GPS. We gathered data on lion-tailed macaque group locations, demography, and site covariates including trail length, duration of walk, proportion of evergreen forest, height of tallest trees, and human disturbance index. We also performed occupancy modeling using PRESENCE ver. 3.0. We estimated a minimum of 17 groups of macaques in these hills. Low detection and occupancy probabilities indicated a low density of lion-tailed macaques in the study area. Height of the tallest trees correlated positively whereas human disturbance and proportion of evergreen forest correlated negatively with occupancy in grid cells. We also used data from earlier studies carried out in the surrounding Anamalai Tiger Reserve and Nelliyampathy Hills to discuss the conservation status in the large Anamalai Hills Landscape. This landscape harbors an estimated population of 1108 individuals of lion-tailed macaques, which is about one third of the entire estimated wild population of this species. A conservation plan for this landscape could be used as a model for conservation in other regions of the Western Ghats.

Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

Role of interface curvature on stress distribution under indentation for ZrN/Zr multilayer coating

Contact damage in curved interface nano-layeredmetal/nitride (150 (ZrN)/10 (Zr) nm) multilayer is investigated in order to understand the role of interface morphology on contact damage under indentation. A finite element method (FEM) model was formulated with different wavelengths of 1000 nm, 500 nm, 250 nm and common height of 50 nm, which gives insight on the effect of different curvature on stress field generated under indentation. Elastic-plastic properties were assigned to the metal layer and substrate while the nitride layer was assigned perfectly elastic properties. Curved interface multilayers show delamination along the metal/nitride interface and vertical cracks emanating from the ends of the delamination. FEM revealed the presence of tensile stress normal to the interface even under the contact, along with tensile radial stresses, both present at the valley part of the curve, which leads to vertical cracks associated with interfacial delamination. Stress enhancement was seen to be relatively insensitive to curvature. (C) 2014 Elsevier B.V. All rights reserved.

Double Gaussian distribution of barrier height observed in densely packed GaN nanorods over Si (111) heterostructures

GaN nanorods were grown by plasma assisted molecular beam epitaxy on intrinsic Si (111) substrates which were characterized by powder X-ray diffraction, field emission scanning electron microscopy, and photoluminescence. The current-voltage characteristics of the GaN nanorods on Si (111) heterojunction were obtained from 138 to 493K which showed the inverted rectification behavior. The I-V characteristics were analyzed in terms of thermionic emission model. The temperature variation of the apparent barrier height and ideality factor along with the non-linearity of the activation energy plot indicated the presence of lateral inhomogeneities in the barrier height. The observed two temperature regimes in Richardson's plot could be well explained by assuming two separate Gaussian distribution of the barrier heights. (C) 2014 AIP Publishing LLC.