Validation of a Modified Eddy Dissipation Concept Model for Stationary and Nonstationary Turbulent Diffusion Flames

A transient flame simulation tool based on unsteady Reynolds average Navier Stokes (RANS) is characterized for stationary and nonstationary flame applications with a motivation of performing computationally affordable flame stability studies. Specifically, the KIVA-3V code is utilized with incorporation of a recently proposed modified eddy dissipation concept for simulating turbulence-chemistry interaction along with a model for radiation loss. Detailed comparison of velocities, turbulent kinetic energies, temperature, and species are made with the experimental data of the turbulent, non-premixed DLR_A CH4/H-2/N-2 jet flame. The comparison shows that the model is able to predict flame structure very well. The effect of some of the modeling assumptions is assessed, and strategies to model a stationary diffusion flame are recommended. Unsteady flame simulation capabilities of the numerical model are assessed by simulating an acoustically excited, experimental, oscillatory H-2-air diffusion flame. Comparisons are made with oscillatory velocity field and OH plots, and the numerical code is observed to predict transient flame structure well.

Growth mechanism of the Nb(X)Si-2 and [Nb(X)](5)Si-3 phases by reactive diffusion in Nb (X = Ti, Mo, or Zr)-Si systems

The effects of Mo, Ti, and Zr on the diffusion and growth of the Nb(X)Si-2 and Nb(X)(5)Si-3 phases in an Nb(X)-Si system are analyzed. The integrated diffusion coefficients are determined from diffusion couple experiments and compared with the data previously calculated in a binary Nb-Si system. The growth rates of both phases are affected by the addition of Mo and Zr, whereas the addition of Ti has no effect. The atomic mechanism of diffusion is also discussed based on the crystal structure and the possible changes in the defect concentrations due to alloying. Finally, the growth mechanism of the phases is discussed on the basis of a physico-chemical approach. (C) 2011 Elsevier Ltd. All rights reserved.

Decoupling of diffusion from viscosity: Difference scenario for translational and rotational motions

Recent experiments have indicated a dramatically different viscosity dependence of the translational and the rotational diffusion coefficients in a supercooled liquid as the glass transition temperature is approached from above. While the translational motion seems to be decoupled from the rising viscosity (eta), the rotational motion seems to remain firmly coupled to eta. In order to understand the microscopic origin of this behavior, we have carried nut detailed theoretical calculations of both the quantities by using a self-consistent mode-coupling theory (MCT). it is found that when the size of the solute is same as that of the solvent molecules, the conventional MCT fails to predict the observed decoupling. The solvent inhomogeneity is found to play a decisive role in determining the decoupling. The difference in the viscosity dependence between rotation and translational diffusion coefficient is discussed.

Self-diffusion of colloids at freezing

This is an introduction to the theory of interacting Brownian particles, as applied to charge-stabilised colloidal suspensions near their equilibrium liquid-solid transition. The density functional approach to the statics of the transition is reviewed briefly, and the generalised Langevin equation method for the dynamics presented in detail. Work with A.V. Indrani [1] on a self-consistent approach for calculating the excess single-particle friction is presented, which explains the observed [2] ''universal'' suppression of self-diffusion at freezing as a consequence of the universal structure-factor height at this transition. Criticisms, open questions, and challenges to theory are discussed.

Shear-enhanced diffusion in charge-stabilised colloids

The generalised Langevin equation method for the dynamics of interacting colloids presented in my previous lecture is extended here to the case of a sheared suspension. A calculation of shear-dependent diffusivities using these methods is found to account for puzzling observations in experiments and simulations. The limitations of the method are discussed, and important unresolved questions presented. This lecture summarises work done in collaboration with A.V. Indrani [1].

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.