Phase-field elasticity model based on mechanical jump conditions

Computational models based on the phase-field method typically operate on a mesoscopic length scale and resolve structural changes of the material and furthermore provide valuable information about microstructure and mechanical property relations. An accurate calculation of the stresses and mechanical energy at the transition region is therefore indispensable. We derive a quantitative phase-field elasticity model based on force balance and Hadamard jump conditions at the interface. Comparing the simulated stress profiles calculated with Voigt/Taylor (Annalen der Physik 274(12):573, 1889), Reuss/Sachs (Z Angew Math Mech 9:49, 1929) and the proposed model with the theoretically predicted stress fields in a plate with a round inclusion under hydrostatic tension, we show the quantitative characteristics of the model. In order to validate the elastic contribution to the driving force for phase transition, we demonstrate the absence of excess energy, calculated by Durga et al. (Model Simul Mater Sci Eng 21(5):055018, 2013), in a one-dimensional equilibrium condition of serial and parallel material chains. To validate the driving force for systems with curved transition regions, we relate simulations to the Gibbs-Thompson equilibrium condition

A Variant of the Hadwiger-Debrunner (p, q)-Problem in the Plane

Let X be a convex curve in the plane (say, the unit circle), and let be a family of planar convex bodies such that every two of them meet at a point of X. Then has a transversal of size at most . Suppose instead that only satisfies the following ``(p, 2)-condition'': Among every p elements of , there are two that meet at a common point of X. Then has a transversal of size . For comparison, the best known bound for the Hadwiger-Debrunner (p, q)-problem in the plane, with , is . Our result generalizes appropriately for if is, for example, the moment curve.

Draft genome of a commonly misdiagnosed multidrug resistant pathogen Candida auris

Background: Candida auris is a multidrug resistant, emerging agent of fungemia in humans. Its actual global distribution remains obscure as the current commercial methods of clinical diagnosis misidentify it as C. haemulonii. Here we report the first draft genome of C. auris to explore the genomic basis of virulence and unique differences that could be employed for differential diagnosis. Results: More than 99.5 % of the C. auris genomic reads did not align to the current whole (or draft) genome sequences of Candida albicans, Candida lusitaniae, Candida glabrata and Saccharomyces cerevisiae; thereby indicating its divergence from the active Candida clade. The genome spans around 12.49 Mb with 8527 predicted genes. Functional annotation revealed that among the sequenced Candida species, it is closest to the hemiascomycete species Clavispora lusitaniae. Comparison with the well-studied species Candida albicans showed that it shares significant virulence attributes with other pathogenic Candida species such as oligopeptide transporters, mannosyl transfersases, secreted proteases and genes involved in biofilm formation. We also identified a plethora of transporters belonging to the ABC and major facilitator superfamily along with known MDR transcription factors which explained its high tolerance to antifungal drugs. Conclusions: Our study emphasizes an urgent need for accurate fungal screening methods such as PCR and electrophoretic karyotyping to ensure proper management of fungemia. Our work highlights the potential genetic mechanisms involved in virulence and pathogenicity of an important emerging human pathogen namely C. auris. Owing to its diversity at the genomic scale; we expect the genome sequence to be a useful resource to map species specific differences that will help develop accurate diagnostic markers and better drug targets.

Nonlinear control of an air-breathing engine including its validation with vehicle guidance

An implementable nonlinear control design approach is presented for a supersonic air-breathing ramjet engine. The primary objective is to ensure that the thrust generated by the engine tracks the commanded thrust without violating the operational constraints. An important constraint is to manage the shock wave location in the intake so that it neither gets detached nor gets too much inside the intake. Both the objectives are achieved by regulating the fuel flow to the combustion chamber and by varying the throat area of the nozzle simultaneously. The design approach accounts for the nonlinear cross-coupling effects and nullifies those. Also, an extended Kalman filter has been used to filter out the sensor and process noises as well as to make the states available for feedback. Furthermore, independent control design has been carried out for the actuators. To test the performance of the engine for a realistic flight trajectory, a representative trajectory is generated through a trajectory optimization process, which is augmented with a newly-developed finite-time state dependent Riccati equation technique for nullifying the perturbations online. Satisfactory overall performance has been obtained during both climb and cruise phases. (C) 2015 Elsevier Masson SAS. All rights reserved.

X-Ha center dot center dot center dot C hydrogen bonds in n-alkane-HX (X = F, OH) complexes are stronger than C-Ha <-X hydrogen bonds

Computational study of X-Ha <-C and C-Ha <-X hydrogen bonds in n-alkane-HX complexes (X =F,OH, alkane =propane, butane, pentane) has been carried out in this work. Ab initio and density functional theories were used for this study. For n-alkane-H2O complexes both Oa <-H-C and O-Ha <-C hydrogen bonded complex have been found, while for n-alkane-HF complexes, our attempt to optimize Fa <-H-C H-bond was not successful. Like most of the hydrogen bonded systems, strong correlation between binding energy and stretching frequency of H-F and O-H stretching mode was observed. The values of electron density and Laplacian of electron density are within the accepted range for hydrogen bonds. In all these cases, X-Ha <-C hydrogen bonds are found to be stronger than C-Ha <-X hydrogen bonds.

Simultaneous acquisition of three NMR spectra in a single experiment for rapid resonance assignments in metabolomics

NMR-based approach to metabolomics typically involves the collection of two-dimensional (2D) heteronuclear correlation spectra for identification and assignment of metabolites. In case of spectral overlap, a 3D spectrum becomes necessary, which is hampered by slow data acquisition for achieving sufficient resolution. We describe here a method to simultaneously acquire three spectra (one 3D and two 2D) in a single data set, which is based on a combination of different fast data acquisition techniques such as G-matrix Fourier transform (GFT) NMR spectroscopy, parallel data acquisition and non-uniform sampling. The following spectra are acquired simultaneously: (1) C-13 multiplicity edited GFT (3,2)D HSQC-TOCSY, (2) 2D H-1- H-1] TOCSY and (3) 2D C-13- H-1] HETCOR. The spectra are obtained at high resolution and provide high-dimensional spectral information for resolving ambiguities. While the GFT spectrum has been shown previously to provide good resolution, the editing of spin systems based on their CH multiplicities further resolves the ambiguities for resonance assignments. The experiment is demonstrated on a mixture of 21 metabolites commonly observed in metabolomics. The spectra were acquired at natural abundance of C-13. This is the first application of a combination of three fast NMR methods for small molecules and opens up new avenues for high-throughput approaches for NMR-based metabolomics.

Intramolecular hydrogen bond: Can it be part of the basis set of valence internal coordinates in normal mode analysis?

It has been shown earlier1] that the relaxed force constants (RFCs) could be used as a measure of bond strength only when the bonds form a part of the complete valence internal coordinates (VIC) basis. However, if the bond is not a part of the complete VIC basis, its RFC is not necessarily a measure of bond strength. Sometimes, it is possible to have a complete VIC basis that does not contain the intramolecular hydrogen bond (IMHB) as part of the basis. This means the RFC of IMHB is not necessarily a measure of bond strength. However, we know that IMHB is a weak bond and hence its RFC has to be a measure of bond strength. We resolve this problem of IMHB not being part of the complete basis by postulating `equivalent' basis sets where IMHB is part of the basis at least in one of the equivalent sets of VIC. As long as a given IMHB appears in one of the equivalent complete VIC basis sets, its RFC could be used as a measure of bond strength parameter.

Does size matter? Comparative population genetics of two butterflies with different wingspans

The dispersal ability of a species is central to its biology, affecting other processes like local adaptation, population and community dynamics, and genetic structure. Among the intrinsic, species-specific factors that affect dispersal ability in butterflies, wingspan was recently shown to explain a high amount of variance in dispersal ability. In this study, a comparative approach was adopted to test whether a difference in wingspan translates into a difference in population genetic structure. Two closely related butterfly species from subfamily Satyrinae, family Nymphalidae, which are similar with respect to all traits that affect dispersal ability except for wingspan, were studied. Melanitis leda (wingspan 60-80 mm) and Ypthima baldus (wingspan 30-40 mm) were collected from the same areas along the Western Ghats of southern India. Amplified fragment length polymorphisms were used to test whether the species with a higher wingspan (M. leda) exhibited a more homogenous population genetic structure, as compared to a species with a shorter wingspan (Y. baldus). In all analyses, Y. baldus exhibited greater degree of population genetic structuring. This study is one of the few adopting a comparative approach to establish the relationship between traits that affect dispersal ability and population genetic structure.

Morphological and functional characterisation of the burrow systems of six earthworm species (Lumbricidae)

Earthworm burrow systems are generally described based on postulated behaviours associated with the three ecological types. In this study, we used X-ray tomography to obtain 3D information on the burrowing behaviour of six very common anecic (Aporrectodea nocturna and Lumbricus terrestris) and endogeic (Aporrectodea rosea, Allolobophora chlorotica, Aporrectodea caliginosa, Aporrectodea icterica) earthworm species, introduced into repacked soil cores for 6 weeks. A simple water infiltration test, the Beerkan method, was also used to assess some functional properties of these burrow systems. Endogeic worms make larger burrow systems, which are more highly branched, less continuous and of smaller diameter, than those of anecic worms. Among the anecic species, L. terrestris burrow systems are shorter (9.2 vs 21.2 m) with a higher number (14.5 vs 23.5) of less branched burrows (12.2 vs 20.2 branches m(-1)), which are also wider (7.78 vs 5.16 mm) than those of A. nocturna. In comparison, the burrow systems made by endogeic species appeared similar to each other. However, A. rosea burrows were short and narrow, whereas A. icterica had a longer burrow system (15.7 m), more intense bioturbation intensity (refilled macropores or soil lateral compaction around them) and thus a greater number of burrows. Regarding water infiltration, anecic burrow systems were far more efficient due to open burrows linking the top and bottom of the cores. For endogeic species, we observed a linear relationship between burrow length and the water infiltration rate (R (2) = 0.49, p < 0.01). Overall, the three main characteristics significantly influencing water infiltration were burrow length, burrow number and bioturbation volume. This last characteristic highlighted the effect of burrow refilling by casts.

Quantitative metabolic profiling of NMR spectral signatures of branched chain amino acids in blood serum

Branched Chain Amino Acids (BCAAs) are related to different aspects of diseases like pathogenesis, diagnosis and even prognosis. While in some diseases, levels of all the BCAAs are perturbed; in some cases, perturbation occurs in one or two while the rest remain unaltered. In case of ischemic heart disease, there is an enhanced level of plasma leucine and isoleucine but valine level remains unaltered. In `Hypervalinemia', valine is elevated in serum and urine, but not leucine and isoleucine. Therefore, identification of these metabolites and profiling of individual BCAA in a quantitative manner in body-fluid like blood plasma/serum have long been in demand. H-1 NMR resonances of the BCAAs overlap with each other which complicates quantification of individual BCAAs. Further, the situation is limited by the overlap of broad resonances of lipoprotein with the resonances of BCAAs. The widely used commercially available kits cannot differentially estimate the BCAAs. Here, we have achieved proper identification and characterization of these BCAAs in serum in a quantitative manner employing a Nuclear Magnetic Resonance-based technique namely T-2-edited Correlation Spectroscopy (COSY). This approach can easily be extended to other body fluids like bile, follicular fluids, saliva, etc.

Coupling geomorphology parameters with lithology for micro-level groundwater resource assessment: a case study from semi-arid hard rock terrain in Tamil Nadu, India

The standard procedure of groundwater resource estimation in India till date is based on the specific yield parameters of each rock type (lithology) derived through pumping test analysis. Using the change in groundwater level, specific yield, and area of influence, groundwater storage change could be estimated. However, terrain conditions in the form of geomorphological variations have an important bearing on the net groundwater recharge. In this study, an attempt was made to use both lithology and geomorphology as input variables to estimate the recharge from different sources in each lithology unit influenced by the geomorphic conditions (lith-geom), season wise separately. The study provided a methodological approach for an evaluation of groundwater in a semi-arid hard rock terrain in Tirunelveli, Tamil Nadu, India. While characterizing the gneissic rock, it was found that the geomorphologic variations in the gneissic rock due to weathering and deposition behaved differently with respect to aquifer recharge. The three different geomorphic units identified in gneissic rock (pediplain shallow weathered (PPS), pediplain moderate weathered (PPM), and buried pediplain moderate (BPM)) showed a significant variation in recharge conditions among themselves. It was found from the study that Peninsular gneiss gives a net recharge value of 0.13 m/year/unit area when considered as a single unit w.r.t. lithology, whereas the same area considered with lith-geom classes gives recharge values between 0.1 and 0.41 m/year presenting a different assessment. It is also found from this study that the stage of development (SOD) for each lith-geom unit in Peninsular gneiss varies from 168 to 230 %, whereas the SOD is 223 % for the lithology as a single unit.

Projected climate change impacts on vegetation distribution over Kashmir Himalayas

Despite high vulnerability, the impact of climate change on Himalayan ecosystem has not been properly investigated, primarily due to the inadequacy of observed data and the complex topography. In this study, we mapped the current vegetation distribution in Kashmir Himalayas from NOAA AVHRR and projected it under A1B SRES, RCP-4.5 and RCP-8.5 climate scenarios using the vegetation dynamics model-IBIS at a spatial resolution of 0.5A degrees. The distribution of vegetation under the changing climate was simulated for the 21st century. Climate change projections from the PRECIS experiment using the HADRM3 model, for the Kashmir region, were validated using the observed climate data from two observatories. Both the observed as well as the projected climate data showed statistically significant trends. IBIS was validated for Kashmir Himalayas by comparing the simulated vegetation distribution with the observed distribution. The baseline simulated scenario of vegetation (1960-1990), showed 87.15 % agreement with the observed vegetation distribution, thereby increasing the credibility of the projected vegetation distribution under the changing climate over the region. According to the model projections, grasslands and tropical deciduous forests in the region would be severely affected while as savannah, shrubland, temperate evergreen broadleaf forest, boreal evergreen forest and mixed forest types would colonize the area currently under the cold desert/rock/ice land cover types. The model predicted that a substantial area of land, presently under the permanent snow and ice cover, would disappear by the end of the century which might severely impact stream flows, agriculture productivity and biodiversity in the region.

Factors affecting distribution of fish within a tidally drained mangrove forest in the Andaman and Nicobar Islands, India

Mangrove forests in meso-tidal areas are completely drained during low tides, forming only temporary habitats for fish. We hypothesised that in such temporary habitats, where stranding risks are high, distance from tidal creeks that provided access to inundated areas during receding tides would be the primary determinant of fish distribution. Factors such as depth, root density and shade were hypothesised to have secondary effects. We tested these hypotheses in a tidally drained mangrove patch in the Andaman Islands, India. Using stake nets, we measured fish abundance and species richness relative to distance from creeks, root density/m(2), shade, water depth and size (total length) of fish. We also predicted that larger fish (including potential predators) would be closer to creeks, as they faced a greater chance of mortality if stranded. Thus we conducted tethering trials to examine if predation would be greater close to the creeks. Generalised linear mixed effects models showed that fish abundance was negatively influenced by increasing creek distance interacting with fish size and positively influenced by depth. Quantile regression analysis showed that species richness was limited by increasing creek distance. Proportion of predation was greatest close to the creeks (0-25 m) and declined with increasing distance. Abundance was also low very close to the creeks, suggesting that close to the creeks predation pressure may be an important determinant of fish abundance. The overall pattern however indicates that access to permanently inundated areas, may be an important determinant of fish distribution in tidally drained mangrove forests.

A hybrid method for stochastic response analysis of a vibrating structure

Response analysis of a linear structure with uncertainties in both structural parameters and external excitation is considered here. When such an analysis is carried out using the spectral stochastic finite element method (SSFEM), often the computational cost tends to be prohibitive due to the rapid growth of the number of spectral bases with the number of random variables and the order of expansion. For instance, if the excitation contains a random frequency, or if it is a general random process, then a good approximation of these excitations using polynomial chaos expansion (PCE) involves a large number of terms, which leads to very high cost. To address this issue of high computational cost, a hybrid method is proposed in this work. In this method, first the random eigenvalue problem is solved using the weak formulation of SSFEM, which involves solving a system of deterministic nonlinear algebraic equations to estimate the PCE coefficients of the random eigenvalues and eigenvectors. Then the response is estimated using a Monte Carlo (MC) simulation, where the modal bases are sampled from the PCE of the random eigenvectors estimated in the previous step, followed by a numerical time integration. It is observed through numerical studies that this proposed method successfully reduces the computational burden compared with either a pure SSFEM of a pure MC simulation and more accurate than a perturbation method. The computational gain improves as the problem size in terms of degrees of freedom grows. It also improves as the timespan of interest reduces.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.