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.

Impact of physical processes on the chlorophyll distribution in the Bay of Bengal

Satellite-derived chlorophyll a concentration (chl a) maps show three regions with high chl a in the Bay of Bengal. First among these is close to the coast, particularly off river mouths, with high values coinciding with the season of peak discharge; second is in the southwestern bay during the northeast monsoon, which is forced by local Ekman pumping; and the third is to the east of Sri Lanka in response to the summer monsoon winds. Chlorophyll-rich water from the mouths of rivers flows either along the coast or in an offshore direction, up to several hundred kilometers, depending on the prevailing ocean current pattern. The Irrawady River plume flows toward offshore and then turns northwestward during October–December, but it flows along the coast into the Andaman Sea for the rest of the year. From the Ganga-Brahmaputra river mouth, chl a–rich water flows directly southward into the open bay during spring but along the Indian coast during summer and winter. Along the Indian coast, the flow of chl a–rich water is determined by the East India Coastal Current (EICC). Whenever the EICC meanders off the Indian coast, it leads to an offshore outbreak of chl a–rich water from the coastal region into open ocean. The EICC as well as open ocean circulation in the bay is made up of several eddies, and these eddies show relatively higher chl a. Eddies near the coast, however, can often have higher chl a because of advection from the coastal region rather than generation within the eddy itself. The bay experiences several cyclones in a year, most of them occurring during October–November. These cyclones cause a drop in the sea surface temperature, a dip in the sea level, and a local increase in chl a. The impact of a cyclone is weaker in the northern part of the bay because of stronger stratification compared to the southern parts.

Surface acoustic wave distribution and acousto-optic interaction in proton exchanged LiNbO3 waveguides

The efficiency of acoustooptic (AO) interaction in YZ-cut proton exchanged (PE) LiNbO3 waveguides is theoretically analysed by determining the overlap between the optical and acoustic field distributions. The present analysis takes into account the perturbed SAW field distribution due to the presence of the PE layer on the LiNbO3 substrate determined by the rigorous layered medium approach. The overlap is found to be significant upto very high acoustic frequencies of the order of 5 GHz, whereas in the earlier analysis by vonHelmolt and Schaffer [6] for diffused waveguides, it was shown that the overlap integral rolls down to nearly zero at this high frequency range.

Numerical Electromagnetic Analysis of Current and Voltage Distribution along Down Conductors and Towers Hit by Direct Lightning

In order to answer the practically important question of whether the down conductors of lightning protection systems to tall towers and buildings can be electrically isolated from the structure itself, this work is conducted. As a first step in this regard, it is presumed that the down conductor placed on metallic tower will be a pessimistic representation of the actual problem. This opinion was based on the fact that the proximity of heavy metallic structure will have a large damping effect. The post-stroke current distributions along the down conductors and towers, which can be quite different from that in the lightning channel, govern the post-stroke near field and the resulting gradient in the soil. Also, for a reliable estimation of the actual stroke current from the measured down conductor currents, it is essential to know the current distribution characteristics along the down conductors. In view of these, the present work attempts to deduce the post-stroke current and voltage distribution along typical down conductors and towers. A solution of the governing field equations on an electromagnetic model of the system is sought for the investigation. Simulation results providing the spatio-temporal distribution of the post-stroke current and voltage has provided very interesting results. It is concluded that it is almost impossible to achieve electrical isolation between the structure and the down conductor. Furthermore, there will be significant induction into the steel matrix of the supporting structure.

Texture evolution in commercially pure titanium after warm equal channel angular extrusion

Texture development in commercially pure titanium during equal channel angular extrusion (ECAE) through Routes A, Be and C has been studied up to three passes at 400 C. Textures were measured using X-ray diffraction, while the microstructural analyses were performed using electron back-scattered diffraction as well as transmission electron microscopy. Occurrences of dynamic restoration processes (recovery and recrystallization) were clearly noticed at all levels of deformations. Finally, the textures were simulated using a viscoplastic polycrystal self-consistent (VPSC) model. Simulations were performed incorporating basal, prismatic and pyramidal slip systems as well as tensile and compressive twinning. The simulated textures corroborate well with experimental textures in spite of the occurrence of dynamic restoration processes. (C) 2010 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Microstructure evolution and hardness variation during annealing of equal channel angular pressed ultra-fine grained nickel subjected to 12 passes

The microstructure, thermal stability and hardness of ultra-fine grained (UFG) Ni produced by 12 passes of equal channel angular pressing (ECAP) through the route Bc were studied. Comparing the microstructure and hardness of the as-ECAPed samples with the published data on UFG Ni obtained after 8 passes of ECAP through the route Bc reveals a smaller average grain size (230 nm in the present case compared with 270 nm in 8-pass Ni), significantly lower dislocation density (1.08 x 10(14) m(-2) compared with 9 x 10(14) m(-2) in 8-pass Ni) and lower hardness (2 GPa compared with 2.45 GPa for 8-pass Ni). Study of the thermal stability of the 12-pass UFG Ni revealed that recovery is dominant in the temperature range 150-250A degrees C and recrystallisation occurred at temperatures > 250 A degrees C. The UFG microstructure is relatively stable up to about 400 A degrees C. Due to the lower dislocation density and consequently a lower stored energy, the recrystallisation of 12-pass ECAP Ni occurred at a higher temperature (similar to 250 A degrees C) compared with the 8-pass Ni (similar to 200 A degrees C). In the 12-pass Nickel, hardness variation shows that its dependence on grain size is inversely linear rather than the common grain size(-0.5) dependence.

Geographical Distribution Patterns of Amphibians in the Western Ghats, India

The Western Ghats of India are very rich in amphibian species with 117 species of frogs, toads and caecilians. Eighty-nine species are endemic to this biogeographical region. Analysis of ranges and patterns of geographical distribution of amphibians on the Western Ghats suggest that the southern half of the Western Ghats and the low-medium elevation hills are more diverse in species than the northern half and higher hills. This is attributed to the more widespread rainfall and the less variable climatic conditions in the south. About half the species are apparently localized. Of those species with wider ranges, a majority show patchy distribution. Species preferring the moist evergreen forests as habitats tend to have patchy distributions. This appears to be a result of habitat destruction and fragmentation. The overall patterns of species richness and local endemism are rather different from those of the angiosperms and birds. In birds and angiosperms, a significant proportion of endemics are found on the higher hills. On the contrary, endemic amphibian species are found in the lower altitudinal range of 0-1000 m, with a majority between 800 and 1000 m.

An unusual distribution of glucose-6-phosphate dehydrogenase deficiency of south Indian newborn population.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is seen at a higher frequency in many national and ethnic groups in areas of current or former malaria endemicity. A screening programme undertaken to evaluate the gene frequencies for this deficiency in the highly inbred South Indian population of Karnataka revealed that of the 5140 neonates screened, 7.8% were G6PD deficient with no correlation between the reported level of inbreeding and enzyme deficiency. An interesting finding was the equal number of male (198) and female (207) individuals, with G6PD activity of less than 3 IU. The possible implications of this finding with regard to the expression of G6PD gene is discussed.

Distribution of ionic charge carriers and migration barriers in binary alkali silicate glasses

Likely spatial distributions of network-modifying (and mobile) cations in (oxide) glasses are discussed here. At very low modifier concentrations, the ions form dipoles with non-bridging oxygen centres while, at higher levels of modification, the cations tend to order as a result of Coulombic interactions. Activation energies for cation migration are calculated, assuming that the ions occupy (face-sharing) octahedral sites. It is found that conductivity activation energy decreases markedly with increasing modifier content, in agreement with experiment.

Xenon in sodium Y zeolite. 2. Arrhenius relation, mechanism, and barrier height distribution for cage-to-cage diffusion

Various geometrical and energetic distribution functions and other properties connected with the cage-to-cage diffusion of xenon in sodium Y zeolite have been obtained from long molecular dynamics calculations. Analysis of diffusion pathways reveals two interesting mechanisms-surface-mediated and centralized modes for cage-to-cage diffusion. The surface-mediated mode of diffusion exhibits a small positive barrier, while the centralized diffusion exhibits a negative barrier for the sorbate to diffuse across the 12-ring window. In both modes, however, the sorbate has to be activated from the adsorption site to enable it to gain mobility. The centralized diffusion additionally requires the sorbate to be free of the influence of the surface of the cage as well. The overall rate for cage-to-cage diffusion shows an Arrhenius temperature dependence with E(a) = 3 kJ/mol. It is found that the decay in the dynamical correction factor occurs on a time scale comparable to the cage residence time. The distributions of barrier heights have been calculated. Functions reflecting the distribution of the sorbate-zeolite interaction at the window and the variations of the distance between the sorbate and the centers of the parent and daughter cages are presented.

Parameter estimation in water-distribution systems by least squares

The weighted-least-squares method using sensitivity-analysis technique is proposed for the estimation of parameters in water-distribution systems. The parameters considered are the Hazen-Williams coefficients for the pipes. The objective function used is the sum of the weighted squares of the differences between the computed and the observed values of the variables. The weighted-least-squares method can elegantly handle multiple loading conditions with mixed types of measurements such as heads and consumptions, different sets and number of measurements for each loading condition, and modifications in the network configuration due to inclusion or exclusion of some pipes affected by valve operations in each loading condition. Uncertainty in parameter estimates can also be obtained. The method is applied for the estimation of parameters in a metropolitan urban water-distribution system in India.

Solvation dynamics, energy distribution and trapping of a light solute ion

In the theoretical treatments of the dynamics of solvation of a newly created ion in a dipolar solvent, the self-motion of the solute is usually ignored. Recently, it has been shown that for a light ion the translational motion of the ion can significantly enhance its own rate of solvation. Therefore, solvation itself may not be the rate determining step in the equilibration. Instead, the rate determining step is the search of the low energy configuration which serves to localize the light ion. In this article a microscopic calculation of the probability distribution of the interaction energy of the nascent charge with the dipolar solvent molecules is presented in order to address this problem of solute trapping. It is found that to a good approximation, this distribution is Gaussian and the second moment of this distribution is exactly equal to the half of its own solvation energy. It is shown that this is in excellent agreement with the simulation results that are available for the model Brownian dipolar lattice and for liquid acetonitrile. If the distortion of the solvent by the ion is negligible then the same relation gives the energy distribution for the solvated ion, with the average centered at the final equilibrium solvation energy. These results are expected to be useful in understanding various chemical processes in dipolar liquids. Another interesting outcome of the present study is a simple dynamic argument that supports Onsager's ''inverse snow-ball'' conjecture of solvation of a light ion. A simple derivation of the semi-phenomenological relation between the solvation time correlation function and the single particle orientation, reported recently by Maroncelli et al. (J. Phys. Chem. 97 (1993) 13), is also presented.

Synthesis, Crystal Structure, and Magnetic Properties of a Ferromagnetically Coupled Angular Trinuclear Copper(11) Complex [Cu3(02CMe)4(bpy)s(HzO)](PF)62

The synthesis, X-ray crystal structure, and magnetic properties of an angular trinuclear copper(II) complex [Cu3(O2CMC)4(bpy)3(H2O)](PF6)2 (1), obtained from a reaction of Cu2(O2CMe)4(H2O)2 With 2,2'-bipyridine (bpy) and NH4PF6 in ethanol, are reported. Complex 1 crystallizes in triclinic space group P1BAR with a = 11.529(1) angstrom, b = 12.121(2) angstrom, c = 17.153(2) angstrom, alpha = 82.01(1)-degrees, beta = 79.42(1)-degrees, gamma = 89.62(1)-degrees, and Z = 2. A total of 6928 data with I > 2.5sigma(I) were refined to R = 0.0441 and R(w) = 0.0557. The structure consists of a trinuclear core bridged by four acetate ligands showing different bonding modes. The coordination geometry at each copper is distorted square-pyramidal with a CuN2O2...O chromophore. The Cu...Cu distances are 3.198(1) angstrom, 4.568(1) angstrom, and 6.277(1) angstrom. There are two monoatomic acetate bridges showing Cu-O-Cu angles of 93.1(1) and 97.5(1)-degrees. Magnetic studies in the temperature range 39-297 K show the presence of a strong ferromagnetically coupled dicopper(II) unit (2J = +158 cm-1) and an essentially isolated copper(II) center (2J' = -0.4 cm-1) in 1. The EPR spectra display an axial spectrum giving g(parallel-to) = 2.28 (A(parallel-to) = 160 X 10(-4) cm-1) and g(perpendicular-to) = 2.06 (A(perpendicular-to) = 12 X 10(-4) cm-1) for the normal copper and two intense isotropic signals with g values 2.70 and 1.74 for the strongly coupled copper pair. The structural features of 1 compare well with the first generation models for ascorbate oxidase.

Single-molecule thermodynamics: the heat distribution function of a charged particle in a static magnetic field

Interest in the applicability of fluctuation theorems to the thermodynamics of single molecules in external potentials has recently led to calculations of the work and total entropy distributions of Brownian oscillators in static and time-dependent electromagnetic fields. These calculations, which are based on solutions to a Smoluchowski equation, are not easily extended to a consideration of the other thermodynamic quantity of interest in such systems-the heat exchanges of the particle alone-because of the nonlinear dependence of the heat on a particle's stochastic trajectory. In this paper, we show that a path integral approach provides an exact expression for the distribution of the heat fluctuations of a charged Brownian oscillator in a static magnetic field. This approach is an extension of a similar path integral approach applied earlier by our group to the calculation of the heat distribution function of a trapped Brownian particle, which was found, in the limit of long times, to be consistent with experimental data on the thermal interactions of single micron-sized colloids in a viscous solvent.