Viscous flow computations using a meshless solver, LSFD-U

The Upwind-Least Squares Finite Difference (LSFD-U) scheme has been successfully applied for inviscid flow computations. In the present work, we extend the procedure for computing viscous flows. Different ways of discretizing the viscous fluxes are analysed for the positivity, which determines the robustness of the solution procedure. The scheme which is found to be more positive is employed for viscous flux computation. The numerical results for validating the procedure are presented.

Novel composite coded pattern for small angle measurement using imaging method -art.no.62891D

A novel approach for measurement of small rotation angles using imaging method is proposed and demonstrated. A plane mirror placed on a precision rotating table is used for imaging the newly designed composite coded pattern. The imaged patterns are captured with the help of a CCD camera. The angular rotation of the plane mirror is determined from a pair of the images of the pattern, captured once before and once after affecting the tilt of the mirror. Both simulation and experimental results suggest that the proposed approach not only retains the advantages of the original imaging method but also contributes significantly to the enhancement of its measuring range (+/- 4.13 degrees with accuracy of the order of 1 arcsec).

Elastic property estimation using Ultrasound Assisted Optical Elastography through Remote Palpation - A simulation study - art. no. 61640Q

We propose an effective elastography technique in which an acoustic radiation force is used for remote palpation to generate localized tissue displacements, which are directly correlated to localized variations of tissue stiffness and are measured using a light probe in the same direction of ultrasound propagation. The experimental geometry has provision to input light beam along the ultrasound propagation direction, and hence it can be prealigned to ensure proper interception of the focal region by the light beam. Tissue-mimicking phantoms with homogeneous and isotropic mechanical properties of normal and malignant breast tissue are considered for the study. Each phantom is insonified by a focusing ultrasound transducer (1 MHz). The focal volume of the transducer and the ultrasound radiation force in the region are estimated through solving acoustic wave propagation through medium assuming average acoustic properties. The forward elastography problem is solved for the region of insonification assuming the Lame's parameters and Poisson's ratio, under Dirichlet boundary conditions which gives a distribution of displacement vectors. The direction of displacement, though presented spatial variation, is predominantly towards the ultrasound propagation direction. Using Monte Carlo (MC) simulation we have traced the photons through the phantom and collected the photons arriving at the detector on the boundary of the object in the direction of ultrasound. The intensity correlations are then computed from detected photons. The intensity correlation function computed through MC simulation showed a modulation whose strength is found to be proportional to the amplitude of displacement and inversely related to the storage (elastic) modulus. It is observed that when the storage modulus in the focal region is increased the computed displacement magnitude, as indicated by the depth of modulation in the intensity autocorrelation, decreased and the trend is approximately exponential.

Computationally efficient optical tomographic reconstructions through waveletizing the normalized quadratic perturbation equation - art. no.61640R

In this paper, we present a wavelet - based approach to solve the non-linear perturbation equation encountered in optical tomography. A particularly suitable data gathering geometry is used to gather a data set consisting of differential changes in intensity owing to the presence of the inhomogeneous regions. With this scheme, the unknown image, the data, as well as the weight matrix are all represented by wavelet expansions, thus yielding the representation of the original non - linear perturbation equation in the wavelet domain. The advantage in use of the non-linear perturbation equation is that there is no need to recompute the derivatives during the entire reconstruction process. Once the derivatives are computed, they are transformed into the wavelet domain. The purpose of going to the wavelet domain, is that, it has an inherent localization and de-noising property. The use of approximation coefficients, without the detail coefficients, is ideally suited for diffuse optical tomographic reconstructions, as the diffusion equation removes most of the high frequency information and the reconstruction appears low-pass filtered. We demonstrate through numerical simulations, that through solving merely the approximation coefficients one can reconstruct an image which has the same information content as the reconstruction from a non-waveletized procedure. In addition we demonstrate a better noise tolerance and much reduced computation time for reconstructions from this approach.

Elimination of cross-talk in diffuse optical tomographic reconstructions through a weighted update procedure: results of simulation study - art. no. 61640S

Reconstructions in optical tomography involve obtaining the images of absorption and reduced scattering coefficients. The integrated intensity data has greater sensitivity to absorption coefficient variations than scattering coefficient. However, the sensitivity of intensity data to scattering coefficient is not zero. We considered an object with two inhomogeneities (one in absorption and the other in scattering coefficient). The standard iterative reconstruction techniques produced results, which were plagued by cross talk, i.e., the absorption coefficient reconstruction has a false positive corresponding to the location of scattering inhomogeneity, and vice-versa. We present a method to remove cross talk in the reconstruction, by generating a weight matrix and weighting the update vector during the iteration. The weight matrix is created by the following method: we first perform a simple backprojection of the difference between the experimental and corresponding homogeneous intensity data. The built up image has greater weightage towards absorption inhomogeneity than the scattering inhomogeneity and its appropriate inverse is weighted towards the scattering inhomogeneity. These two weight matrices are used as multiplication factors in the update vectors, normalized backprojected image of difference intensity for absorption inhomogeneity and the inverse of the above for the scattering inhomogeneity, during the image reconstruction procedure. We demonstrate through numerical simulations, that cross-talk is fully eliminated through this modified reconstruction procedure.

Sensitive high-resolution ion microprobe U-Pb dating of prograde and retrograde ultrahigh-temperature metamorphism as exemplified by Sri Lankan granulites

Ultrahigh-temperature (UHT) granulites of the central Highland Complex, Sri Lanka, underwent some of the highest known peak temperatures of crustal metamorphism. Zircon and monazite U-Pb systems in granulites near Kandy, the highest grade region (similar to 1050 degrees C; 0.9 GPa), preserve both a record of the timing of prograde and retrograde phases of UHT metamorphism and evidence for the ages of older protolith components. Zircon grains from a quartz-saturated granulite containing relics of the peak UHT assemblage have remnant detrital cores with dates of ca. 2.5-0.83 Ga. Date clusters of ca. 1.7 and 1.04-0.83 Ga record episodes of zircon growth in the source region of the protolith sediment. Two generations of overgrowths with contrasting Th/U record metamorphic zircon growth at 569 +/- 5 and 551 +/- 7 Ma, probably in the absence and presence of monazite, respectively. The age of coexisting metamorphic monazite (547 +/- 7 Ma) is indistinguishable from that of the younger, low-Th/U zircon overgrowths. Zircon from a quartz-undersaturated monazite-absent UHT granulite with a mainly retrograde assemblage is mostly metamorphic (551 +/- 5 Ma). The ca. 570 Ma zircon overgrowths in the quartz-saturated granulite probably record partial melting just before or at the metamorphic peak. The ca. 550 Ma zircon in both rocks, and the ca. 550 Ma monazite in the quartz-saturated sample, record post-peak isothermal decompression. A possible model for this pressure-temperature-time evolution is ultrahot collisional orogeny during the assembly of Gondwana, locally superheated by basaltic underplating, followed by fast extensional exhumation.

Formation and preservation of pedogenic carbonates in South India, links with paleo-monsoon and pedological conditions: Clues from Sr isotopes, U-Th series and REEs

The influence of the pedogenic and climatic contexts on the formation and preservation of pedogenic carbonates in a climosequence in the Western Ghats (Karnataka Plateau, South West India) has been studied. Along the climosequence, the current mean annual rainfall (MAR) varies within a 80 km transect from 6000 mm at the edge of the Plateau to 500 mm inland. Pedogenic carbonates occur in the MAR range of 500-1200 mm. In the semi-arid zone (MAR: 500-900 mm), carbonates occur (i) as rhick hardpan calcretes on pediment slopes and (ii) as nodular horizons in polygenic black soils (i.e. vertisols). In the sub-humid zone (MAR: 900-1500 mm), pedogenic carbonates are disseminated in the black soil matrices either as loose, irregular and friable nodules of millimetric size or as indurated botryoidal nodules of centimetric to pluricentimetric size. They also occur at the top layers of the saprolite either as disseminated pluricentimetric indurated nodules or carbonate-cemented lumps of centimetric to decimetric size. Chemical and isotopic (Sr-87/Sr-86) compositions of the carbonate fraction were determined after leaching with 0.25 N HCl. The corresponding residual fractions containing both primary minerals and authigenic clays were digested separately and analyzed. The trend defined by the Sr-87/Sr-86 signatures of both labile carbonate fractions and corresponding residual fractions indicates that a part of the labile carbonate fraction is genetically linked to the local soil composition. Considering the residual fraction of each sample as the most likely lithogenic source of Ca in carbonates, it is estimated that from 24% to 82% (55% on average) of Ca is derived from local bedrock weathering, leading to a consumption of an equivalent proportion of atmospheric CO2. These values indicate that climatic conditions were humid enough to allow silicate weathering: MAR at the time of carbonate formation likely ranged from 400 to 700 mm, which is 2- to 3-fold less than the current MAR at these locations. The Sr, U and Mg contents and the (U-234/U-238) activity ratio in the labile carbonate fraction help to understand the conditions of carbonate formation. The relatively high concentrations of Sr, U and Mg in black soil carbonates may indicate fast growth and accumulation compared to carbonates in saprolite, possibly due to a better confinement of the pore waters which is supported by their high (U-234/U-238) signatures, and/or to higher content of dissolved carbonates in the pore waters. The occurrence of Ce, Mn and Fe oxides in the cracks of carbonate reflects the existence of relatively humid periods after carbonate formation. The carbonate ages determined by the U-Th method range from 1.33 +/- 0.84 kyr to 7.5 +/- 2.7 kyr and to a cluster of five ages around 20 kyr, i.e. the Last Glacial Maximum period. The young occurrences are only located in the black soils, which therefore constitute sensitive environments for trapping and retaining atmospheric CO2 even on short time scales. The maximum age of carbonates depends on their location in the climatic gradient: from about 20 kyr for centimetric nodules at Mule Hole (MAR = 1100 mm/yr) to 200 kyr for the calcrete at Gundlupet (MAR = 700 mm/yr, Durand et al., 2007). The intensity of rainfall during wet periods would indeed control the lifetime of pedogenic carbonates and thus the duration of inorganic carbon storage in soils. (C) 2010 Elsevier Ltd. All rights reserved.

Turbulent flow computations on a hybrid cartesian point distribution using meshless solver LSFD-U

This paper may be considered as a sequel to one of our earlier works pertaining to the development of an upwind algorithm for meshless solvers. While the earlier work dealt with the development of an inviscid solution procedure, the present work focuses on its extension to viscous flows. A robust viscous discretization strategy is chosen based on positivity of a discrete Laplacian. This work projects meshless solver as a viable cartesian grid methodology. The point distribution required for the meshless solver is obtained from a hybrid cartesian gridding strategy. Particularly considering the importance of an hybrid cartesian mesh for RANS computations, the difficulties encountered in a conventional least squares based discretization strategy are highlighted. In this context, importance of discretization strategies which exploit the local structure in the grid is presented, along with a suitable point sorting strategy. Of particular interest is the proposed discretization strategies (both inviscid and viscous) within the structured grid block; a rotated update for the inviscid part and a Green-Gauss procedure based positive update for the viscous part. Both these procedures conveniently avoid the ill-conditioning associated with a conventional least squares procedure in the critical region of structured grid block. The robustness and accuracy of such a strategy is demonstrated on a number of standard test cases including a case of a multi-element airfoil. The computational efficiency of the proposed meshless solver is also demonstrated. (C) 2010 Elsevier Ltd. All rights reserved.

Exotic Physics in the Negative-U, Extended-Hubbard Model for Barium Bismuthates?

We point out possibilities for exotic physics in barium bismuthates, from a detailed study of the negative-U, extended-Hubbard model proposed for these systems. We emphasize the different consequences of electronic and phononic mechanisms for negative U. We show that, for an electronic mechanism, the semiconducting phases must be unique, with their transport properties dominated by charge ± 2e Cooperon bound states. This can explain the observed difference between the optical and transport gaps. We propose other experimental tests for this novel mechanism of charge transport.

Squeezed states, photon-number distributions, and U(l) invariance

We study the photon-number distribution in squeezed states of a single-mode radiation field. A U(l)-invariant squeezing criterion is compared and contrasted with a more restrictive criterion, with the help of suggestive geometric representations. The U(l) invariance of the photon-number distribution in a squeezed coherent state, with arbitrary complex squeeze and displacement parameters, is explicitly demonstrated. The behavior of the photon-number distribution for a representative value of the displacement and various values of the squeeze parameter is numerically investigated. A new kind of giant oscillation riding as an envelope over more rapid oscillations in this distribution is demonstrated.

Prediction of the Indian summer monsoon rainfall using a state-of-the-art coupled ocean-atmosphere model

A state-of-the-art model of the coupled ocean-atmosphere system, the climate forecast system (CFS), from the National Centres for Environmental Prediction (NCEP), USA, has been ported onto the PARAM Padma parallel computing system at the Centre for Development of Advanced Computing (CDAC), Bangalore and retrospective predictions for the summer monsoon (June-September) season of 2009 have been generated, using five initial conditions for the atmosphere and one initial condition for the ocean for May 2009. Whereas a large deficit in the Indian summer monsoon rainfall (ISMR; June-September) was experienced over the Indian region (with the all-India rainfall deficit by 22% of the average), the ensemble average prediction was for above-average rainfall during the summer monsoon. The retrospective predictions of ISMR with CFS from NCEP for 1981-2008 have been analysed. The retrospective predictions from NCEP for the summer monsoon of 1994 and that from CDAC for 2009 have been compared with the simulations for each of the seasons with the stand-alone atmospheric component of the model, the global forecast system (GFS), and observations. It has been shown that the simulation with GFS for 2009 showed deficit rainfall as observed. The large error in the prediction for the monsoon of 2009 can be attributed to a positive Indian Ocean Dipole event seen in the prediction from July onwards, which was not present in the observations. This suggests that the error could be reduced with improvement of the ocean model over the equatorial Indian Ocean.

Observation of asymmetry in U�I characteristics of contacts of the form metal 1�superconductor�metal 2

During experiments carried out to find out a suitable contact metal for electronic components based on high-T(c) superconductor films (Y-Ba-Cu-O), it is observed that there is an asymmetry in the U-I characteristics if the two contacts are made of different metals. The asymmetry is more pronounced if one of the contact metals is aluminium. The asymmetry is lowest if one of the contact metals is silver and the other gold.

On the solution of the equation ut+unux+H(x, t, u)=0

We consider the equation u(t) + u(n)u(x) + H(x, t, u) = 0 and derive a transformation relating it to u(t) + u(n)u(x) = 0. Special cases of the equation appearing in applications are discussed. Initial value problems and asymptotic behaviour of the solution are studied.

Role of zirconium addition on stability of UAl3 phase in Al-U system

The microstructural changes of Al-22 wt%U and Al-46 wt%U alloys containing 3 wt% Zr were investigated after heat treatment at 620 degrees C for 1 to 45 days, Though it is reported that addition of similar to 3 wt% Zr stabilizes the (U,Zr)Al-3 phase at room temperature, the present investigation shows that the (U,Zr)Al-3 phase is not stable but slowly transforms to the U0.9Al4 phase, The high temperature creep curves generated for these ternary alloys showed a wavy pattern which also suggests that the (U,Zr)Al-3 phase is not stable.

Microstructure evolution during solidification of DRMMCs (discontinuously reinforced metal matrix composites): State of art

Distribution of particle reinforcements in cast composites is determined by the morphology of the solidification front. Interestingly, during solidification, the morphology of the interface is intrinsically affected by the presence of dispersed reinforcements. Thus the dispersoid distribution and length scale of matrix microstructure is a result of the interplay between these two. A proper combination of material and process parameters can be used to obtain composites with tailored microstructures. This requires the generation of a broad data base and optimization of the complete solidification process. The length scale of soldification microtructure has a large influence on the mechanical properties of the composites. This presentation addresses the concept of a particle distribution map which can help in predicting particle distribution under different solidification conditions Future research directions have also been indicated.