Efficient asynchronous executions of AMR computations and visualization on a GPU system

Adaptive Mesh Refinement is a method which dynamically varies the spatio-temporal resolution of localized mesh regions in numerical simulations, based on the strength of the solution features. In-situ visualization plays an important role for analyzing the time evolving characteristics of the domain structures. Continuous visualization of the output data for various timesteps results in a better study of the underlying domain and the model used for simulating the domain. In this paper, we develop strategies for continuous online visualization of time evolving data for AMR applications executed on GPUs. We reorder the meshes for computations on the GPU based on the users input related to the subdomain that he wants to visualize. This makes the data available for visualization at a faster rate. We then perform asynchronous executions of the visualization steps and fix-up operations on the CPUs while the GPU advances the solution. By performing experiments on Tesla S1070 and Fermi C2070 clusters, we found that our strategies result in 60% improvement in response time and 16% improvement in the rate of visualization of frames over the existing strategy of performing fix-ups and visualization at the end of the timesteps.

Natural motion of one-dimensional flexible objects using minimization approaches

For one-dimensional flexible objects such as ropes, chains, hair, the assumption of constant length is realistic for large-scale 3D motion. Moreover, when the motion or disturbance at one end gradually dies down along the curve defining the one-dimensional flexible objects, the motion appears ``natural''. This paper presents a purely geometric and kinematic approach for deriving more natural and length-preserving transformations of planar and spatial curves. Techniques from variational calculus are used to determine analytical conditions and it is shown that the velocity at any point on the curve must be along the tangent at that point for preserving the length and to yield the feature of diminishing motion. It is shown that for the special case of a straight line, the analytical conditions lead to the classical tractrix curve solution. Since analytical solutions exist for a tractrix curve, the motion of a piecewise linear curve can be solved in closed-form and thus can be applied for the resolution of redundancy in hyper-redundant robots. Simulation results for several planar and spatial curves and various input motions of one end are used to illustrate the features of motion damping and eventual alignment with the perturbation vector.

Out of equilibrium plasticity dynamics and the annealing of supersolidity in solid He-4

We present a numerical study of a continuum plasticity field coupled to a Ginzburg-Landau model for superfluidity. The results suggest that a supersolid fraction may appear as a long-lived transient during the time evolution of the plasticity field at higher temperatures where both dislocation climb and glide are allowed. Supersolidity, however, vanishes with annealing. As the temperature is decreased, dislocation climb is arrested and any residual supersolidity due to incomplete annealing remains frozen. Our results may provide a resolution of many perplexing issues concerning a variety of experiments on bulk solid He-4.

Regionalization of precipitation in India-a review

Regionalization of precipitation refers to delineation of rain gauges in an area into homogeneous groups (clusters or regions). Various regionalization procedures are employed by researchers in hydrometeorology for addressing a wide spectrum of problems. This paper provides an overview of underlying concepts as well as advantages and limitations of procedures that have been developed over the past six decades. Emphasis is given to studies that have been carried out in India. Following this, gaps where more research needs to be focussed are highlighted, and challenges for regionalization in a climate change scenario are discussed.

Investigations of Ramachandran disallowed conformations in protein domain families

In peptide and protein structures, occurrence of (phi,psi.) angles in the disallowed region of the Ramachandran map almost always suggests local regions of error or poor accuracy. However, very rarely genuine disallowed conformations occur as noted in the current study in proteins of known structure available at ultra-high resolution (<= 1.2 (A) over circle). In the current work, extent of conservation of genuine disallowed conformations in homologous proteins of known structures has been analyzed. From a dataset of 124 protein domain families, with structure of at least one constituent member in each family available at a resolution of 1.2 (A) over circle or better, we have analyzed the conservation of 221 disallowed conformations. It is observed that the disallowed conformation is only moderately conservedin protein domain families. In the gross dataset no particular residue type adopting disallowed conformation elicit high conservation of residue type though there are alignment positions in the dataset with complete conservation of both the residue type and the disallowed conformation. Conserved disallowed conformation in protein domain families play biologically significant role in roughly 50% of the cases. The residues with the disallowed conformation or its flanking residues are often located within or around the functional site of the protein. (C) 2013 Elsevier B.V. All rights reserved.

A simple model for spatial disaggregation of evaporative fraction: Comparative study with thermal sharpened land surface temperature data over India

[1] Evaporative fraction (EF) is a measure of the amount of available energy at the earth surface that is partitioned into latent heat flux. The currently operational thermal sensors like the Moderate Resolution Imaging Spectroradiometer (MODIS) on satellite platforms provide data only at 1000 m, which constraints the spatial resolution of EF estimates. A simple model (disaggregation of evaporative fraction (DEFrac)) based on the observed relationship between EF and the normalized difference vegetation index is proposed to spatially disaggregate EF. The DEFrac model was tested with EF estimated from the triangle method using 113 clear sky data sets from the MODIS sensor aboard Terra and Aqua satellites. Validation was done using the data at four micrometeorological tower sites across varied agro-climatic zones possessing different land cover conditions in India using Bowen ratio energy balance method. The root-mean-square error (RMSE) of EF estimated at 1000 m resolution using the triangle method was 0.09 for all the four sites put together. The RMSE of DEFrac disaggregated EF was 0.09 for 250 m resolution. Two models of input disaggregation were also tried with thermal data sharpened using two thermal sharpening models DisTrad and TsHARP. The RMSE of disaggregated EF was 0.14 for both the input disaggregation models for 250 m resolution. Moreover, spatial analysis of disaggregation was performed using Landsat-7 (Enhanced Thematic Mapper) ETM+ data over four grids in India for contrasted seasons. It was observed that the DEFrac model performed better than the input disaggregation models under cropped conditions while they were marginally similar under non-cropped conditions.

Structural and kinetic studies on adenylosuccinate lyase from Mycobacterium smegmatis and Mycobacterium tuberculosis provide new insights on the catalytic residues of the enzyme

Adenylosuccinate lyase (ASL), an enzyme involved in purine biosynthesis, has been recognized as a drug target against microbial infections. In the present study, ASL from Mycobacteriumsmegmatis (MsASL) and Mycobacteriumtuberculosis (MtbASL) were cloned, purified and crystallized. The X-ray crystal structure of MsASL was determined at a resolution of 2.16 angstrom. It is the first report of an apo-ASL structure with a partially ordered active site C3 loop. Diffracting crystals of MtbASL could not be obtained and a model for its structure was derived using MsASL as a template. These structures suggest that His149 and either Lys285 or Ser279 of MsASL are the residues most likely to function as the catalytic acid and base, respectively. Most of the active site residues were found to be conserved, with the exception of Ser148 and Gly319 of MsASL. Ser148 is structurally equivalent to a threonine in most other ASLs. Gly319 is replaced by an arginine residue in most ASLs. The two enzymes were catalytically much less active compared to ASLs from other organisms. Arg319Gly substitution and reduced flexibility of the C3 loop might account for the low catalytic activity of mycobacterial ASLs. The low activity is consistent with the slow growth rate of Mycobacteria and their high GC containing genomes, as well as their dependence on other salvage pathways for the supply of purine nucleotides. Structured digital abstract andby()

Freeform Skeletal Shape Optimization of Compliant Mechanisms

Compliant mechanisms are elastic continua used to transmit or transform force and motion mechanically. The topology optimization methods developed for compliant mechanisms also give the shape for a chosen parameterization of the design domain with a fixed mesh. However, in these methods, the shapes of the flexible segments in the resulting optimal solutions are restricted either by the type or the resolution of the design parameterization. This limitation is overcome in this paper by focusing on optimizing the skeletal shape of the compliant segments in a given topology. It is accomplished by identifying such segments in the topology and representing them using Bezier curves. The vertices of the Bezier control polygon are used to parameterize the shape-design space. Uniform parameter steps of the Bezier curves naturally enable adaptive finite element discretization of the segments as their shapes change. Practical constraints such as avoiding intersections with other segments, self-intersections, and restrictions on the available space and material, are incorporated into the formulation. A multi-criteria function from our prior work is used as the objective. Analytical sensitivity analysis for the objective and constraints is presented and is used in the numerical optimization. Examples are included to illustrate the shape optimization method.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

Desingularization of the Milne Universe

Resolution of cosmological singularities is an important problem in any full theory of quantum gravity. The Milne orbifold is a cosmology with a big-bang/big-crunch singularity, but being a quotient of flat space it holds potential for resolution in string theory. It is known, however, that some perturbative string amplitudes diverge in the Milne geometry. Here we show that flat space higher spin theories can effect a simple resolution of the Milne singularity when one embeds the latter in 2 + 1 dimensions. We explain how to reconcile this with the expectation that non-perturbative string effects are required for resolving Milne. Along the way, we introduce a Grassmann realization of the inonfi-Wigner contraction to export much of the AdS technology to -our flat space computation. (C) 2014 The Authors. Published by Elsevier BAT.

A non-resonant mass sensor to eliminate the ``missing mass'' effect during mass measurement of biological materials

Resonant sensors and crystal oscillators for mass detection need to be excited at very high natural frequencies (MHz). Use of such systems to measure mass of biological materials affects the accuracy of mass measurement due to their viscous and/or viscoelastic properties. The measurement limitation of such sensor system is the difficulty in accounting for the ``missing mass'' of the biological specimen in question. A sensor system has been developed in this work, to be operated in the stiffness controlled region at very low frequencies as compared to its fundamental natural frequency. The resulting reduction in the sensitivity due to non-resonant mode of operation of this sensor is compensated by the high resolution of the sensor. The mass of different aged drosophila melanogaster (fruit fly) is measured. The difference in its mass measurement during resonant mode of operation is also presented. That, viscosity effects do not affect the working of this non-resonant mass sensor is clearly established by direct comparison. (C) 2014 AIP Publishing LLC.

Super Resolution Using a Single Image Dictionary

To perform super resolution of low resolution images, state-of-the-art methods are based on learning a pair of lowresolution and high-resolution dictionaries from multiple images. These trained dictionaries are used to replace patches in lowresolution image with appropriate matching patches from the high-resolution dictionary. In this paper we propose using a single common image as dictionary, in conjunction with approximate nearest neighbour fields (ANNF) to perform super resolution (SR). By using a common source image, we are able to bypass the learning phase and also able to reduce the dictionary from a collection of hundreds of images to a single image. By adapting recent developments in ANNF computation, to suit super-resolution, we are able to perform much faster and accurate SR than existing techniques. To establish this claim, we compare the proposed algorithm against various state-of-the-art algorithms, and show that we are able to achieve b etter and faster reconstruction without any training.

Design of a Portable Compliant Device for Estimating the Failure-Load of Mesoscale Cemented Sand Specimens

In this paper, we present the design and development of a portable, hand-operated composite compliant mechanism for estimating the failure-load of cm-sized stiff objects whose stiffness is of the order of 10 s of kN/m. The motivation for the design comes from the need to estimate the failure-load of mesoscale cemented sand specimens in situ, which is not possible with traditional devices used for large specimens or very small specimens. The composite compliant device, developed in this work, consists of two compliant mechanisms: a force-amplifying compliant mechanism (FaCM) to amplify sufficiently the force exerted by hand in order to break the specimen and a displacement-amplifying compliant mechanism (DaCM) to enable measurement of the force using a proximity sensor. The two mechanisms are designed using the selection-maps technique to amplify the force up to 100N by about a factor of 3 and measure the force with a resolution of 15 mN. The composite device, made using a FaCM, a DaCM, and a Hall effect-based proximity sensor, was tested on mesoscale cemented sand specimens that were 10mm in diameter and 20mm in length. The results are compared with those of a large commercial instrument. Through the experiments, it was observed that the failure-load of the cemented sand specimens varied from 0.95N to 24.33 N, depending on the percentage of cementation and curing period. The estimation of the failure-load using the compliant device was found to be within 1.7% of the measurements obtained using the commercial instrument and thus validating the design. The details of the design, prototyping, specimen preparation, testing, and the results comprise the paper.

Micro-scale composite compliant mechanisms for evaluating the bulk stiffness of MCF-7 cells

Biomechanical assays offer a good alternative to biochemical assays in diagnosing disease states and assessing the efficacy of drugs. In view of this, we have developed a miniature compliant tool to estimate the bulk stiffness of cells, particularly MCF-7 (Michigan Cancer Foundation) cells whose diameter is 12-15 mu m in suspension. The compliant tool comprises a gripper and a displacement-amplifying compliant mechanism (DaCM), where the former helps in grasping the cell and the latter enables vision-based force-sensing. A DaCM is necessary because the microscope's field of view at the required magnification is not sufficient to simultaneously observe the cell and the movement of a point on the gripper, in order to estimate the force. Therefore, a DaCMis strategically embedded within an existing gripper design leading to a composite compliant mechanism. The DaCM is designed using the kinetoelastostatic map technique to achieve a 42 nN resolution of the force. The gripper, microfabricated with SU-8 using photolithography, is within the footprint of about 10 mm by 10 mm with the smallest feature size of about 5 mu m. The experiments with MCF-7 cells suggest that the bulk stiffness of these is in the range of 8090 mN/m. The details of design, prototyping and testing comprise the paper. (C) 2015 Elsevier Ltd. All rights reserved.

Dynamic range tuning of graphene nanoresonators

From sensing perspective, smaller electromechanical devices, in general, are expected to be more responsive to the stimuli. This enhanced performance, however, is contingent upon the noise sources remaining unchanged and the onset of nonlinear behavior not being precipitated by miniaturization. In this paper, we study the effect of strain on the nonlinearities and dynamic range in graphene nanoresonators. The dynamic response and the onset of nonlinearity in these devices are sensitive both to the electrostatic field used to actuate the device and the strain. By tuning the strain of the device by two orders of magnitude, we observe an enhancement of 25 dB in the dynamic range leading to a mass resolution of 100 yoctogram. The increase in dynamic range in our devices is modeled as a combined effect of strain and partial cancellation of elastic and electrostatic nonlinearities. (C) 2015 AIP Publishing LLC.