Molecular dynamics Investigation of structure and transport in the K2O-2SiO2 system using a partial charge based model potential

Potassium disilicate glass and melt have been investigated by using anew partial charge based potential model in which nonbridging oxygens are differentiated from bridging oxygens by their charges. The model reproduces the structural data pertaining to the coordination polyhedra around potassium and the various bond angle distributions excellently. The dynamics of the glass has been studied by using space and time correlation functions. It is found that K ions migrate by a diffusive mechanism in the melt and by hops below the glass transition temperature. They are also found to migrate largely through nonbridging oxygen-rich sites in the silicate matrix, thus providing support to the predictions of the modified random network model.

Precipitation in small systems-I. stochastic analysis

Precipitation in small droplets involving emulsions, microemulsions or vesicles is important for Producing multicomponent ceramics and nanoparticles. Because of the random nature of nucleation and the small number of particles in a droplet, the use of a deterministic population balance equation for predicting the number density of particles may lead to erroneous results even for evaluating the mean behavior of such systems. A comparison between the predictions made through stochastic simulation and deterministic population balance involving small droplets has been made for two simple systems, one involving crystallization and the other a single-component precipitation. The two approaches have been found to yield quite different results under a variety of conditions. Contrary to expectation, the smallness of the population alone does not cause these deviations. Thus, if fluctuation in supersaturation is negligible, the population balance and simulation predictions concur. However, for large fluctuations in supersaturation, the predictions differ significantly, indicating the need to take the stochastic nature of the phenomenon into account. This paper describes the stochastic treatment of populations, which involves a sequence of so-called product density equations and forms an appropriate framework for handling small systems.

Rank dependence of orientational relaxation in dipolar systems

An important yet unsolved problem in the field of orientational relaxation in dipolar liquids is the dependence of the correlation functions C(l)(t), C(l)(t) = [4pi/(2l + 1)SIGMA(m = -l)l [Y(lm)(OMEGA(0)Y(lm)(OMEGA(t))] on the rank l (where Y(lm)(OMEGA) are the usual spherical harmonics). The existing theories on this effect differ in their predictions. To investigate this, we have carried out extensive computer simulations of a Brownian dipolar lattice. The dielectric friction was found to decrease rapidly with increasing l, in qualitative agreement with the predictions of Hubbard-Wolynes. However, the observed effect is much stronger than the predictions of the existing theories.

The gasification of wood-char spheres in CO2---N2 mixtures: analysis and experiments

The gasification of charcoal spheres in an atmosphere of carbon-dioxide-nitrogen mixture involving diffusion and reactions in the pores is modelled and the results are compared with experiments of Standish and Tanjung and those performed in the laboratory on wood-char spheres to determine the effects of diameter, density, gas composition and flow. The results indicate that the conversion time, t(c) approximately d1.03 for large particles (> 5 mm), departing substantially from the t(c) approximately d2 law valid for diffusion limited conditions. The computational studies indicate that the kinetic limit for the particle is below 100 mum. The conversion time varies inversely as the initial char density as expected in the model. Predictions from the model show that there is no significant change in conversion time up to 60% N2 consistent with the CO2-N2 experiments. The variation of diameter and density with time are predicted. The peculiar dependence of conversion time on flow velocity in the experiments is sought to be explained by opposing free and forced convection heat transfer and the attempt is only partly successful. The studies also indicate that the dependence on the CO concentration with low CO2 is significant, indicating the need for multistep reaction mechanism against the generally accepted single-step reaction.

Vertical uplift capacity of horizontal anchors

The method of characteristics coupled with a log-spiral failure surface was used to develop a theory for vertical uplift capacity of shallow horizontal strip anchors in a general c-phi soil. Uplift-capacity factors F(c), F(q) and F(gamma), for the effects of cohesion, surcharge, and density, respectively, have been established as functions of embedment ratio lambda and angle of friction phi. The extent of the failure surface at the ground has also been determined. Comparisons made with existing test results support the predictive capability of the theory, and comparisons with the analysis proposed by Meyerhof and Adams show the proposed analysis provides slightly more conservative predictions of pullout capacity.

Acoustic performance of hoses - A parametric study

Modeling of wave propagation in hoses, unlike in rigid pipes or waveguides, introduces a coupling between the inside medium, the hose wall, and the outside medium, This alters the axial wave number and thence the corresponding effective speed of sound inside the hose resulting in sound radiation into the outside medium, also called the breakout or shell noise, The existing literature on the subject is such that a hose cannot be integrated into the,whole piping system made up of sections of hoses, pipes, and mufflers to predict the acoustical performance in terms of transmission loss (TL), The present paper seeks to fill this gap, Three one-dimensional coupled wave equations are written to account for the presence of a yielding wall with a finite lumped transverse impedance of the hose material, The resulting wave equation can readily be reduced to a transfer matrix form using an effective wave number for a moving medium in a hose section, Incorporating the effect of fluid loading due to the outside medium also allows prediction of the transverse TL and the breakout noise, Axial TL and transverse TL have been combined into net TL needed by designers, Predictions of the axial as well as transverse TL are shown to compare well with those of a rigorous 3-D analysis using only one-hundredth of the computation time, Finally, results of some parametric studies are reported for engineers involved in the acoustical design of hoses. (C) 1996 Institute of Noise Control Engineering.

A new model for drainage of static foams

Existing theories of foam drainage assume bubbles as pentagonal dodecahedrons, though a close-packed structure built with cells of this shape is not space-filling. The present work develops a theory for calculating drainage rates based on the more realistic beta-tetrakaidecahedral shape for the bubbles. In contrast with the earlier works, three types of films, and Plateau borders had to be considered in view of the more complex shape used in the present work. The exchange of liquid between Plateau borders was treated in a way different From earlier theories, using the idea that the volume of junctions of Plateau borders is negligible. For foams made of large bubble sizes, the present model performs as well as the previous models, but when bubble size is small, its predictions of drainage rates from static foams are in better agreement with the experimental observations.

A numerical study of division of flow in open channels

A two-dimensional numerical model which employs the depth-averaged forms of continuity and momentum equations along with k-e turbulence closure scheme is used to simulate the flow at the open channel divisions. The model is generalised to flows of arbitrary geometries and MacCormack finite volume method is used for solving governing equations. Application of cartesian version of the model to analyse the flow at right-angled junction is presented. The numerical predictions are compared with experimental data of earlier investigators and measurements made as part of the present study. Performance of the model in predicting discharge distribution, surface profiles, separation zone parameters and energy losses is evaluated and discussed in detail. To illustrate the application of the numerical model to analyse the flow in acute angled offtakes and streamlined branch entries, a few computational results are presented.

Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Field evaluation of unsaturated hydraulic conductivity models and parameter estimation from retention data

Predictions of two popular closed-form models for unsaturated hydraulic conductivity (K) are compared with in situ measurements made in a sandy loam field soil. Whereas the Van Genuchten model estimates were very close to field measured values, the Brooks-Corey model predictions were higher by about one order of magnitude in the wetter range. Estimation of parameters of the Van Genuchten soil moisture characteristic (SMC) equation, however, involves the use of non-linear regression techniques. The Brooks-Corey SMC equation has the advantage of being amenable to application of linear regression techniques for estimation of its parameters from retention data. A conversion technique, whereby known Brooks-Corey model parameters may be converted into Van Genuchten model parameters, is formulated. The proposed conversion algorithm may be used to obtain the parameters of the preferred Van Genuchten model from in situ retention data, without the use of non-linear regression techniques.

Modeling Of Hydrolytic Polymerization In A Semibatch Nylon 6 Reactor

Hydrolytic polymerization of caprolactam to Nylon 6 in a semibatch reactor is carried out by heating a mixture of water and caprolactam. Evaporation of volatiles caused by heating results in a pressure build-up. After the pressure reaches a predetermined value, vapors are vented to keep the pressure constant for some time, and thereafter, to lower the pressure to a value slightly above atmospheric in a preprogrammed manner. The characteristics of the polymer are determined by the chemical reactions and the vaporization of water and caprolactam. The semibatch operation has been simulated and the predictions have been compared with industria data. The observed temperature and pressure histories were predicted with a fair degree of accuracy. It was found that the predictions of the degree of polymerization however are sensitive to the vapor-liquid equilibrium relations. A comparison with an earlier model, which neglected mass transfer resistance, indicates that simulation using the VLE data of Giori and Hayes and accounting for mass transfer resistance is more reliable.

Climate change and agricultrue-An Indian perspective

In this paper, we examine the major predictions made so far regarding the nature of climate change and its impacts on our region in the light of the known errors of the set of models and the observations over this century. The major predictions of the climate models about the impact of increased concentration of greenhouse gases ave at variance with the observations over the Indian region during the last century characterized by such increases and global warming. It is important to note that as far as the Indian region is concerned, the impact of year-to-year variation of the monsoon will continue to be dominant over longer period changes even in the presence of global warming. Recent studies have also brought out the uncertainties in the yields simulated by crop models. It is suggested that a deeper understanding of the links between climate and agricultural productivity is essential for generating reliable predictions of impact of climate change. Such an insight is also required for identifying cropping patterns and management practices which are tailored for sustained maximum yield in the face of the vagaries of the monsoon.

Industrial validation of processing maps of 316L stainless steel using hot forging, rolling, and extrusion

The development of microstructure in 316L stainless steel during industrial hot forming operations including press forging (strain rate of 0 . 15 s(-1)), rolling/extrusion (strain rate of 2-8 . 8 s(-1)), and hammer forging (strain rate of 100 s(-1)) at different temperatures in the range 600-1200 degrees C was studied with a view to validating the predictions of the processing map. The results showed that good col relation existed between the regimes indicated in the map and the product microstructures. The 316L stainless steel exhibited unstable flow in the form of flow localisation when hammer forged at temperatures above 900 degrees C, rolled below 1000 degrees C, or press forged below 900 degrees C. All these conditions must therefore be avoided in mechanical processing of the material. Conversely, in order to obtain defect free microstructures, ideally the material should be rolled at temperatures above 1100 degrees C, press forged at temperatures above 1000 degrees C, or hammer forged in the temperature range 600-900 degrees C. (C) 1996 The Institute of Materials.

Magnetic and spin evolution of pulsars

We explore the consequences of the model of spin-down-induced flux expulsion for the magnetic field evolution in solitary as well as in binary neutron stars. The spin evolution of pulsars, allowing for their field evolution according to this model, is shown to be consistent with the existing observational constraints in both low- and high-mass X-ray binary systems. The contribution from pulsars recycled in massive binaries to the observed excess in the number of low-field (10(11)-10(12) G) solitary pulsars is argued to be negligible in comparison with that of normal pulsars undergoing a 'restricted' field decay predicted by the adopted field decay model. Magnetic fields of neutron stars born in close binaries with intermediate- or high-mass main-sequence companions are predicted to decay down to values as low as similar to 10(6) G, which would leave them unobservable as pulsars during most of their lifetimes. The post-recycling evolution of some of these systems can, however, account for the observed binary pulsars having neutron star or massive white dwarf companions. Pulsars recycled in the disc population low-mass binaries are expected to have residual fields greater than or similar to 10(8) G, while for those processed in globular clusters larger residual fields are predicted because of the lower field strength of the neutron star at the epoch of binary formation. A value of tau similar to 1-2 x 10(7) yr for the mean value of the Ohmic decay time-scale in the crusts of neutron stars is suggested, based on the consistency of the model predictions with the observed distribution of periods and magnetic fields in the single and binary pulsars.