A molecular dynamics study of ambient and high pressure phases of silica: Structure and enthalpy variation with molar volume

Extensive molecular dynamics studies of 13 different silica polymorphs are reported in the isothermal-isobaric ensemble with the Parrinello-Rahman variable shape simulation cell. The van Beest-Kramer-van Santen (BKS) potential is shown to predict lattice parameters for most phases within 2%-3% accuracy, as well as the relative stabilities of different polymorphs in agreement with experiment. Enthalpies of high-density polymorphs - CaCl2-type, alpha-PbO2-type, and pyrite-type for which no experimental data are available as yet, are predicted here. Further, the calculated enthalpies exhibit two distinct regimes as a function of molar volume-for low and medium-density polymorphs, it is almost independent of volume, while for high-pressure phases a steep dependence is seen. A detailed analysis indicates that the increased short-range contributions to enthalpy in the high-density phases arise not only from an increased coordination number of silicon but also shorter Si-O bond lengths. Our results indicate that amorphous phases of silica exhibit better optimization of short-range interactions than crystalline phases at the same density while the magnitude of Coulombic contributions is lower in the amorphous phase. (C) 2014 AIP Publishing LLC.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Effect of hydrogen on the yielding behavior and shear transformation zone volume in metallic glass ribbons

The effect of hydrogen (H) charging on the shear yield strength (tau(max)) and shear transformation zone volume (Omega) of Ni-Nb-Zr metallic glass ribbons, with varying Zr content, were studied through the first pop-in loads during nanoindentation. Weight gain measurements after H charging and desorption studies were utilized to identify how the total H absorbed during charging is partitioned into mobile and immobile (or trapped) parts. These, in turn, indicate the significant role of H mobility in the amorphous structure on the yielding behavior. In high-Zr alloys, tau(max) increases significantly whereas Omega decreases. In low-Zr alloys, a slight decrease in tau(max) and increase in Omega were noted. These experimental observations are rationalized in terms of the mobility of the absorbed H in the amorphous structure and the possible role of it in the shear transformation zone dynamics during deformation of the metallic glass. (C) 2014 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Breakup processes of pressure swirl spray in gaseous cross-flow

Injection of liquid fuel in cross flowing air has been a strategy for future aircraft engines in order to control the emissions. In this context, breakup of a pressure swirl spray in gaseous cross-flow is investigated experimentally. The atomizer discharges a conical swirling sheet of liquid that interacts with cross-flowing air. This complex interaction and the resulting spray structures at various flow conditions are studied through flow visualization using still as well as high speed photography. Experiments are performed over a wide range of aerodynamic Weber number (2-300) and liquid-to-air momentum flux ratio (5-150). Various breakup regimes exhibiting different breakup processes are mapped on a parameter space based on flow conditions. This map shows significant variations from breakup regime map for a plain liquid jet in cross-flow. It is observed that the breakup of leeward side of the sheet is dominated by bag breakup and the windward side of the sheet undergoes breakup through surface waves. Similarities and differences between bag breakup present in plain liquid jet in cross-flow and swirl spray in cross-flow are explained. Multimodal drop size distribution from bag breakup, frequency of bag breakup, wavelength of surface waves and trajectory of spray in cross-flow are measured by analyzing the spray images and parametric study of their variations is also presented. (C) 2014 Elsevier Ltd. All rights reserved.

Tin-Incorporation Induced Changes in the Microstructural, Optical, and Electrical Behavior of Tungsten Oxide Nanocrystalline Thin Films Grown Via Spray Pyrolysis

Undoped and Sn-doped WO3 thin films were grown on cleaned glass substrates by chemical spray pyrolysis, using ammonium tungstate (NH4)(2)WO4 as the host precursor and tin chloride (SnCl4 center dot 5H(2)O) as the source of dopant. The XRD spectra confirm the monoclinic structure with a sharp narrow peak along (200) direction along with other peaks of low relative intensities for all the samples. On Sn doping, the films exhibit reduced crystallinity relative to the undoped film. The standard deviation for relative peak intensity with dopant concentration shows enhancement in heterogeneous nucleation growth. As evident from SEM images, on Sn doping, appearance of island-like structure (i.e., cluster of primary crystallites at few places) takes place. The transmittance has been found to decrease in all the Sn-doped films. The optical band gap has been calculated for both direct and indirect transitions. On Sn doping, the direct band gap shows a red shift and becomes 2.89 eV at 2 at.% doping. Two distinct peaks, one blue emission at 408 nm and other green emission at 533 nm, have been found in the PL spectra. Electrical conductivity has been found to increase with Sn doping.

Estimating the volume of glaciers in the Himalayan-Karakoram region using different methods

Ice volume estimates are crucial for assessing water reserves stored in glaciers. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan-Karakoram (HK) region. In this study, different existing methodologies are used to estimate the ice reserves: three area-volume relations, one slope-dependent volume estimation method, and two ice-thickness distribution models are applied to a recent, detailed, and complete glacier inventory of the HK region, spanning over the period 2000-2010 and revealing an ice coverage of 40 775 km(2). An uncertainty and sensitivity assessment is performed to investigate the influence of the observed glacier area and important model parameters on the resulting total ice volume. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 to 4737 km(3), depending on the approach. This range is lower than most previous estimates. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses. However, total volume estimates from area-related relations are larger than those from other approaches. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation.

NUMERICAL SIMULATIONS OF LIQUID JET BREAK UP IN A CROSSFLOW

Atomization is the process of disintegration of a liquid jet into ligaments and subsequently into smaller droplets. A liquid jet injected from a circular orifice into cross flow of air undergoes atomization primarily due to the interaction of the two phases rather than an intrinsic break up. Direct numerical simulation of this process resolving the finest droplets is computationally very expensive and impractical. In the present study, we resort to multiscale modelling to reduce the computational cost. The primary break up of the liquid jet is simulated using Gerris, an open source code, which employs Volume-of-Fluid (VOF) algorithm. The smallest droplets formed during primary atomization are modeled as Lagrangian particles. This one-way coupling approach is validated with the help of the simple test case of tracking a particle in a Taylor-Green vortex. The temporal evolution of the liquid jet forming the spray is captured and the flattening of the cylindrical liquid column prior to breakup is observed. The size distribution of the resultant droplets is presented at different distances downstream from the location of injection and their spatial evolution is analyzed.

SWEETENING OF BIOGAS IN SPRAY TOWERS FOR POWER GENERATION

This paper presents the experience of the new design of using impinging jet spray columns for scrubbing hydrogen sulfide from biogas that has been developed by Indian Institute of Science and patented. The process uses a chelated polyvalent metal ion which oxidizes the hydrogen sulfide to sulfur as a precipitate. The sulfur generated is filtered and the scrubbing liquid recycled after oxidation. The process involves in bringing contact the sour gas with chelated liquid in the spray columns where H2S reacts with chelated Fe3+ and precipitates as sulfur, whereas Fe3+ gets reduced to Fe2+. Fe2+ is regenerated to Fe3+ by reaction of oxygen in air in a separate packed column. The regenerated liquid is recirculated. Sulfur is filtered and separated as a byproduct. The paper presents the experience in using the spray towers for hydrogen sulfide removal and further use of the clean gas for generating power using gas engines. The maximum allowable limit of H2S for the gas engine is 200 ppm (v/v) in order to prevent any corrosion of engine parts and fouling of the lubricating oil. With the current ISET process, the hydrogen sulfide from the biogas is cleaned to less than 100 ppm (v/v) and the sweet gas is used for power generation. The system is designed for 550 NM3/hr of biogas and inlet H2S concentration of 2.5 %. The inlet concentration of the H2S is about 1 - 1.5 % and average measured outlet concentration is about 30 ppm, with an average gas flow of about 300 - 350 NM3/hr, which is the current gas production rate. The sweet gas is used for power generation in a 1.2 MWe V 12 engine. The average power generation is about 650 - 750 kWe, which is the captive load of the industry. The plant is a CHP (combined heat power) unit with heat from the cylinder cooling and flue being recovered for hot water and steam generation respectively. The specific fuel consumption is 2.29 kWh/m(3) of gas. The system has been in operation for more than 13,000 hours in last one year in the industry. About 8.4 million units of electricity has been generated scrubbing about 2.1 million m3 of gas. Performance of the scrubber and the engine is discussed at daily performance level and also the overall performance with an environment sustenance by precipitating over 27 tons of sulfur.

Airblast spray in crossflow - Structure, trajectory and droplet sizing

This study reports results of an experimental investigation of airblast spray of water and ethanol in crossflow. Laser shadowgraphy and Particle/Droplet Imaging Analysis (PDIA) are used to derive spray trajectory and drop size information while Particle Tracking Velocimetry (PTV) is used to measure droplet velocities. A new phenomenon of spray bifurcation is observed for low Gas to Liquid Ratio (GLR) cases. The reasons for the spatial bifurcation can be attributed to a combination of reasons. These are (a) presence of large ligaments and droplets in the near-nozzle region for low GLRs (b) secondary breakup experienced by ligaments/droplets leading to formation of a large number of small droplets, and (c) the crossflow causing differential dispersion of the small and large droplets. A novel correlation for spray trajectory is proposed incorporating the momentum ratio and liquid surface tension. This correlation is shown to be effective in predicting the non-linear spray trajectory over a large range of conditions for not only water but ethanol and Jet-A also. It is observed that the larger droplets penetrate further into the crossflow, in the direction of injection. Thus, with increase in height of the measurement location from the injection plane, the droplet Sauter Mean Diameter (SMD) is found to increase. Moreover, as the droplets travel downstream in the crossflow direction, the droplet SMD is observed to decrease. The effect of drag is assessed by comparing velocity of different sizes of droplets at various locations. Smaller droplets are entrained into the crossflow at much lower elevations, whereas larger droplets tend to penetrate further into the crossflow. (C) 2015 Elsevier Ltd. All rights reserved.

Spray: A Multi-Modal Localization System for Stationary Sensor Network Deployment

We present a localization system that targets rapid deployment of stationary wireless sensor networks (WSN). The system uses a particle filter to fuse measurements from multiple localization modalities, such as RF ranging, neighbor information or maps, to obtain position estimations with higher accuracy than that of the individual modalities. The system isolates different modalities into separate components which can be included or excluded independently to tailor the system to a specific scenario. We show that position estimations can be improved with our system by combining multiple modalities. We evaluate the performance of the system in both an indoor and outdoor environment using combinations of five different modalities. Using two anchor nodes as reference points and combining all five modalities, we obtain RMS (Root Mean Square) estimation errors of approximately 2.5m in both cases, while using the components individually results in errors within the range of 3.5 and 9 m.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

Metal-insulator transition characteristics of vanadium dioxide thin films synthesized by ultrasonic nebulized spray pyrolysis of an aqueous combustion mixture

We report the synthesis of high quality vanadium dioxide (VO2) thin films by a novel spray pyrolysis technique, namely ultrasonic nebulized spray pyrolysis of aqueous combustion mixture (UNSPACM). This simple and cost effective two step process involves synthesis of a V2O5 film on an LaAlO3 substrate followed by a controlled reduction to form single phase VO2. The formation of M1 phase (p21/c) is confirmed by Raman spectroscopic studies. A thermally activated metal-insulator transition (MIT) was observed at 61 degrees C, where the resistivity changes by four orders of magnitude. Activation energies for the low conduction phase and the high conduction phase were obtained from temperature variable resistance measurements. The infrared spectra also show a dramatic change in reflectance from 13% to over 90% in the wavelength range of 7-15 mu m. This indicates the suitability of the films for optical switching applications at infrared frequencies.

Role of the Molar Volume on Estimated Diffusion Coefficients

The role of the molar volume on the estimated diffusion parameters has been speculated for decades. The Matano-Boltzmann method was the first to be developed for the estimation of the variation of the interdiffusion coefficients with composition. However, this could be used only when the molar volume varies ideally or remains constant. Although there are no such systems, this method is still being used to consider the ideal variation. More efficient methods were developed by Sauer-Freise, Den Broeder, and Wagner to tackle this problem. However, there is a lack of research indicating the most efficient method. We have shown that Wagner's method is the most suitable one when the molar volume deviates from the ideal value. Similarly, there are two methods for the estimation of the ratio of intrinsic diffusion coefficients at the Kirkendall marker plane proposed by Heumann and van Loo. The Heumann method, like the Matano-Boltzmann method, is suitable to use only when the molar volume varies more or less ideally or remains constant. In most of the real systems, where molar volume deviates from the ideality, it is safe to use the van Loo method. We have shown that the Heumann method introduces large errors even for a very small deviation of the molar volume from the ideal value. On the other hand, the van Loo method is relatively less sensitive to it. Overall, the estimation of the intrinsic diffusion coefficient is more sensitive than the interdiffusion coefficient.

THE EFFECTS OF SURFACTANT ON SIMPLEX NOZZLE SPRAY BEHAVIOR AND ITS COMPARISON TO LIQUID FUELS

Pressure-swirl nozzles (simplex nozzles) are used in various field applications such as aero-engines, power generation, spray painting and agricultural irrigation. For this particular nozzle, research in the past decade has dealt with the development of numerical models for predicting droplet distribution profiles. Although these results have been valuable, the experimental results have been contradictory, therefore fundamental understanding of the influence of properties in nozzle is important. This paper experimentally investigates the effect of surfactants on breakup and coalescence. Since most of the fuels and biofuels have low surface tension compared to water, a comparative analysis between a surfactant solution and a liquid fuel is imperative. For this experimental study, a simplex nozzle characterized as flow number 0.4 will be utilized. The injection pressures will range from 0.3 - 4Mpa while altering the surface tension from 72 to 28mN/m. By applying Phase Doppler Particle Anemometry (PDPA) which is a non-intrusive laser diagnostic technique, the differences in spray characteristics due to spray surface tension can be highlighted. The average droplet diameter decreases for a low surface tension fluid in the axial direction in comparison to pure water. The average velocity of droplets is surprisingly lower in the same spray zone. Measurements made in the radial direction show no significant changes, but at the locations close to the nozzle, water droplets have larger diameter and velocity. The results indicate the breakup and coalescence regimes have been altered when surface tension is lowered. A decrease in surface tension alters the breakup length while increasing the spray angle. Moreover, higher injection pressure shortens the breakup length and decrease in overall diameter of the droplets. By performing this experimental study the fundamentals of spray dynamics, such as spray formation, liquid breakup length, and droplet breakup regimes can be observed as a function of surface tension and how a surrogate fuel compares with a real fuel for experimental purposes. This knowledge potentially will lead to designing a better atomizer or new biofuels.

Experimental studies on evaporation of fuel droplets under forced convection using spray in crossflow methodology

Experimental data on evaporation of droplets of decane, Jet-A1, and Jet-A1 surrogate are generated using a spray in crossflow configuration. The advantage of a crossflow configuration is that it enables us to study droplet evaporation under forced convective conditions involving droplet diameters of size relevant in practical combustors. Specifically, spray from an airblast atomizer is injected into a preheated crossflow of air and the resulting spray is characterized in terms of spray structure along with droplet size and velocity. An existing correlation for the spray trajectory is modified to incorporate the effect of elevated temperature, and is found to be in good agreement with the experimental data. Droplet sizes and velocities are measured at different locations along the crossflow direction to assess droplet evaporation. Specifically, droplets having size less than 25-mu m are selected for further analysis since these droplets are observed to exhibit velocities which are aligned with the crossflow. By comparing the droplet diameter profiles at upstream and downstream locations, the evaporation constant k for the d(2)-law is obtained iteratively. To assess the efficacy of the values of k obtained, the calculated droplet size distribution using the proposed k values at the downstream location is compared with the measured droplet size distribution at that location. A reasonably good match is found for all the three liquids confirming the validity of the analysis. (C) 2015 Elsevier Ltd. All rights reserved.