Dependence of diffusivity on density and solute diameter in liquid phase: A molecular dynamics study of Lennard-Jones system

Investigations into the variation of self-diffusivity with solute radius, density, and degree of disorder of the host medium is explored. The system consists of a binary mixture of a relatively smaller sized solute, whose size is varied and a larger sized solvent interacting via Lennard-Jones potential. Calculations have been performed at three different reduced densities of 0.7, 0.8, and 0.933. These simulations show that diffusivity exhibits a maximum for some intermediate size of the solute when the solute diameter is varied. The maximum is found at the same size of the solute at all densities which is at variance with the prediction of the levitation effect. In order to understand this anomaly, additional simulations were carried out in which the degree of disorder has been varied while keeping the density constant. The results show that the diffusivity maximum gradually disappears with increase in disorder. Disorder has been characterized by means of the minimal spanning tree. Simulations have also been carried out in which the degree of disorder is constant and only the density is altered. The results from these simulations show that the maximum in diffusivity now shifts to larger distances with decrease in density. This is in agreement with the changes in void and neck distribution with density of the host medium. These results are in excellent agreement with the predictions of the levitation effect. They suggest that the effect of disorder is to shift the maximum in diffusivity towards smaller solute radius while that of the decrease in density is to shift it towards larger solute radius. Thus, in real systems where the degree of disorder is lower at higher density and vice versa, the effect due to density and disorder have opposing influences. These are confirmed by the changes seen in the velocity autocorrelation function, self part of the intermediate scattering function and activation energy. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.3701619]

Diffusion ordered spectroscopy for resolution of double bonded cis, trans-isomers

NMR spectroscopic separation of double bonded cis- and trans-isomers, that have different molecular shapes but identical mass have been carried out using Diffusion Ordered Spectroscopy (DOSY). The mixtures of fumaric acid and maleic acid, that have similar hydrodynamic radii, have resolved been on the basis of their diffusion coefficients arising due to their different tendencies to associate with micelles or reverse micelles. Sodium dodecyl sulfate (SDS) and Dioctyl sulfosuccinate sodium salt (AOT) have been used as the media to mimic the chromatographic conditions, modify the average mobility and to achieve differential diffusion rates. The best separation of the components has been achieved by Dioctyl sulfosuccinate sodium salt (AOT) in D2O solution. (C) 2012 Elsevier B.V. All rights reserved.

Robust solver based on modified particle swarm optimization for improved solution of diffusion transport through containment facilities

Accurate estimation of mass transport parameters is necessary for overall design and evaluation processes of the waste disposal facilities. The mass transport parameters, such as effective diffusion coefficient, retardation factor and diffusion accessible porosity, are estimated from observed diffusion data by inverse analysis. Recently, particle swarm optimization (PSO) algorithm has been used to develop inverse model for estimating these parameters that alleviated existing limitations in the inverse analysis. However, PSO solver yields different solutions in successive runs because of the stochastic nature of the algorithm and also because of the presence of multiple optimum solutions. Thus the estimated mean solution from independent runs is significantly different from the best solution. In this paper, two variants of the PSO algorithms are proposed to improve the performance of the inverse analysis. The proposed algorithms use perturbation equation for the gbest particle to gain information around gbest region on the search space and catfish particles in alternative iterations to improve exploration capabilities. Performance comparison of developed solvers on synthetic test data for two different diffusion problems reveals that one of the proposed solvers, CPPSO, significantly improves overall performance with improved best, worst and mean fitness values. The developed solver is further used to estimate transport parameters from 12 sets of experimentally observed diffusion data obtained from three diffusion problems and compared with published values from the literature. The proposed solver is quick, simple and robust on different diffusion problems. (C) 2012 Elsevier Ltd. All rights reserved.

Diffusion mechanism in XSi2 and X5Si3 (X= Nb, Mo, V) phases

In view of the importance of the suicides in the high temperature applications, the diffusion behaviour is compared in different systems for two types of silicides, XSi2 and X5Si3 (X=Nb, Mo, V). Atomic mechanism of diffusion and defects present in the structure are discussed. In both the phases, Si has faster diffusion rate than the metal species. This is expected from the nearest neighbour (NN) bonds present in the XSi2 phase but rather unusual in the X5Si3 phase. Relative mobilities of the species calculated indicate the presence of high concentration of Si antisites. Moreover, the concentration of the defects is different in different systems to find different diffusion rates.

Diffusion mechanism in the mu phase in Nb-X (X = Ni, Co, Fe) systems

Nb is one of the common refractory elements added in Ni, Co and Fe based superalloys. This lead to the formation of brittle topological close packed (tcp) mu phase, which is deleterious to the structure. It mainly grows by interdiffusion and in the present article, the interdiffusion process in different Nb-X (X=Ni, Co, Fe) systems is discussed. The activation energy for interdiffusion is lower in the Co-Nb system (173 kJ/mol) than Fe-Nb system (233 kJ/mol), which is again lower than the value found in the Ni-Nb system (319.7 kJ/mol). The mole fraction of Nb in this phase is less than Fe or Co at stoichiometric compositions in the Nb-Fe (that is Fe7Nb6) and Nb-Co (that is Co7Nb6) systems. On the other hand, the mole fraction of Nb is higher than Ni in the same phase (Ni6Nb2) in Ni-Nb system. However, in all the phases, Nb has lower diffusion rate. Possible diffusion mechanism in this phase is discussed with respect to the crystal structure.

A versatile resolving agent for diffusion edited separation of enantiomers, complex mixtures and constitutional isomers

The study is the first report of the utilization of a crown ether as a new and versatile resolving agent for the diffusion edited separation of enantiomers, complex mixtures and constitutional isomers. As a consequence of different binding affinities of enantiomers of a chiral molecule and individual components of the complex mixtures with the crown ether, the molecules diffuse at different rates. The enhanced separation achieved due to matrix assisted diffusion permitted their separation in the diffusion dimension. The generality and wide utility of the new resolving agent and the methodology are demonstrated on diverse examples, such as an organic chiral molecule, constitutional isomers and complex mixture of molecules possessing different functional groups that possess nearly identical molecular weights.

Use of Fick's law and Maxwell–Stefan equations in computation of multicomponent diffusion

This article does not have an abstract.

Diffusion mechanism in the gold-copper system

Interdiffusion study is conducted in the Au-Cu system, which has complete solid solution in the higher temperature range and ordered phases in the lower temperature range. First experiments are conducted at higher temperatures, where atoms can diffuse randomly. Higher values of interdiffusion coefficients are found in the range of 40-50 at.% Cu. This trend is explained with the help of thermodynamic factor and possible concentration of vacancies. Following an experiment is conducted at 623 K (350 degrees C), where the ordered phases are grown. The interdiffusion coefficients at this temperature are compared after extrapolating the data calculated at higher temperatures.

Diffusion and growth mechanism of phases in the Pd-Sn system

Growth mechanism of phases and atomic mechanism of diffusion are discussed in the Pd-Sn system. The Kirkendall marker plane location indicates that the PdSn4 phase grows because of diffusion of Sn. Atomic arrangement in the crystal indicates that Sn can diffuse through its own sublattice but Pd cannot diffuse unless antisites are present. The negligible diffusion of Pd indicates the absence of Pd antisites. The activation energy value indicates that the contribution from grain boundary diffusion cannot be neglected although experiments were conducted in the homologous temperature range of 0.7-0.79.

Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.

Temporal Moments for Reactive Transport through Fractured Impermeable/Permeable Formations

The transport of reactive solutes through fractured porous formations has been analyzed. The transport through the porous block is represented by a general multiprocess nonequilibrium equation (MPNE), which, for the fracture, is represented by an advection-dispersion equation with linear equilibrium sorption and first-order transformation. An implicit finite-difference technique has been used to solve the two coupled equations. The transport characteristics have been analyzed in terms of zeroth, first, and second temporal moments of the solute in the fracture. The solute behavior for fractured impermeable and fractured permeable formations are first compared and the effects of various fracture and matrix transport parameters are analyzed. Subsequently, the transport through a fractured permeable formation is analyzed to ascertain the effect of equilibrium sorption, rate-limited sorption, and the multiprocess nonequilibrium transport process. It was found that the temporal moments were nearly identical for the fractured impermeable and permeable formations when both the diffusion coefficient and the first-order transformation coefficient were relatively large. The multiprocess nonequilibrium model resulted in a smaller mass recovery in the fracture and higher dispersion than the equilibrium and rate-limited sorption models. DOI: 10.1061/(ASCE)HE.19435584.0000586. (C) 2012 American Society of Civil Engineers.

Optical properties change with the addition and diffusion of Bi to As2S3 in the Bi/As2S3 bilayer thin film

The role of Bi layer (thickness similar to 7 nm) on As2S3 film was extensively studied for different optical applications in which Bi (top layer) as active and diffusing layer and As2S3 as barrier (matrix) layer. Bilayer thin films of Bi/As2S3 were prepared from Bi and As2S3 by thermal evaporation technique under high vacuum. The decrease of optical band gap with the addition of Bi to As2S3 has been explained on the basis of density of states and the increase in disorder in the system. It was found that the efficient changes of optical parameters (transmission, optical band gap, refraction) could be realized due to the photo induced diffusion activated by the focused 532 nm laser irradiation and formation of different bonds. The diffusion of Bi into As2S3 matrix increases the optical band gap producing photo bleaching effect. The changes were characterised by different experimental techniques. (C) 2012 Elsevier B.V. All rights reserved.

Doping by diffusion of dopants from the substrate: synthesis of doped ZnO nanowires

We report on the substrate assisted doping of ZnO nanowires grown by a vapor transport technique. The nanowires were grown non-catalytically on multiwalled carbon nanotubes (MWCNTs) and soda lime glass (SLG). Carbon from MWCNTs and sodium from SLG diffuse into ZnO during the growth and are distributed uniformly and provide doping. An advantage associated with the technique is that no conventional external dopant source is required to obtain doped ZnO nanowires. The diameter, length and hence the aspect ratio can easily be varied by changing the growth conditions. The transport studies on both carbon and sodium doped ZnO support the p-type nature of ZnO. The p-type nature of carbon doped ZnO is stable for at least eight months.

A computational study of heat transfer, soot production and transport in an acetylene-air laminar diffusion flame in the field of a line vortex

The present work involves a computational study of soot formation and transport in case of a laminar acetylene diffusion flame perturbed by a co nvecting line vortex. The topology of the soot contours (as in an earlier experimental work [4]) have been investigated. More soot was produced when vortex was introduced from the air si de in comparison to a fuel side vortex. Also the soot topography was more diffused in case of the air side vortex. The computational model was found to be in good agreement with the ex perimental work [4]. The computational simulation enabled a study of the various parameters affecting soot transport. Temperatures were found to be higher in case of air side vortex as compared to a fuel side vortex. In case of the fuel side vortex, abundance of fuel in the vort ex core resulted in stoichiometrically rich combustion in the vortex core, and more discrete so ot topography. Overall soot production too was low. In case of the air side vortex abundan ce of air in the core resulted in higher temperatures and more soot yield. Statistical techniques like probability density fun ction, correlation coefficient and conditional probability function were introduced to explain the transient dependence of soot yield and transport on various parameters like temperature, a cetylene concentration.

A computational study of flame dynamics on soot production and transport for an acetylene-air laminar diffusion flame in the field of a convecting line vortex

The present work involves a computational study of soot formation and transport in case of a laminar acetylene diffusion flame perturbed by a convecting line vortex. The topology of the soot contours (as in an earlier experimental work [4]) have been investigated. More soot was produced when vortex was introduced from the air side in comparison to a fuel side vortex. Also the soot topography was more diffused in case of the air side vortex. The computational model was found to be in good agreement with the experimental work [4]. The computational simulation enabled a study of the various parameters affecting soot transport. Temperatures were found to be higher in case of air side vortex as compared to a fuel side vortex. In case of the fuel side vortex, abundance of fuel in the vort ex core resulted in stoichiometrically rich combustion in the vortex core, and more discrete soot topography. Overall soot production too was low. In case of the air side vortex abundance of air in the core resulted in higher temperatures and more soot yield. Statistical techniques like probability density function, correlation coefficient and conditional probability function were introduced to explain the transient dependence of soot yield and transport on various parameters like temperature, a cetylene concentration.