A generalized equation for diffusion in liquids

The association parameter in the diffuswn equaiior, dye fo Wiike one Chong has been interpreted in deferminable properties, thus permitting easily the calculation of the same for unknown systems. The proposed eqyotion a!se holds goods for water as soiute in organic solvenfs. The over-all percentage error remains the sarrse as that of the original equation.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Catalyst effectiveness factor for redox model kinetic equation

A generalized isothermal effectiveness factor correlation has been proposed for catalytic reactions whose intrinsic kinetics are based on the redox model. In this correlation which is exact for asymptotic values of the Thiele parameter the effect of the parameters appearing in the model, the order of the reaction and particle geometry are incorporated in a modified form of Thiele parameter. The relationship takes the usual form: Image and predicts effectiveness factor with an error of less than 2% in a range of Thiele parameter that accommodates both the kinetic and diffusion control regimes.

Sintering of beta titanium: Role of dislocation pipe diffusion

Sintering of titanium in its high temperature beta phase was studied by isothermal dilatometry. The sintering shrinkage y did not follow the normal time exponent type of behaviour, instead being described by the equation y = Kt(m)/[1-(A+Bt)(2)], where m = 1.93 +/- 0.07, with an activation energy of 62-90 kJ mol(-1). A detailed analysis of these results, based on the 'anomalous' diffusion behaviour reported for beta titanium, is carried out. It is shown that the generation of a high density of dislocations during the alpha --> beta phase transformation, coupled with sluggish recovery at the sintering necks, enables sintering mass transport by pipe diffusion through dislocation cores from sources of matter within the particles to become dominant.

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].

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.

Optimal error estimate using mesh equidistribution technique for singularly perturbed system of reaction-diffusion boundary-value problems

In this article, we study the problem of determining an appropriate grading of meshes for a system of coupled singularly perturbed reaction-diffusion problems having diffusion parameters with different magnitudes. The central difference scheme is used to discretize the problem on adaptively generated mesh where the mesh equation is derived using an equidistribution principle. An a priori monitor function is obtained from the error estimate. A suitable a posteriori analogue of this monitor function is also derived for the mesh construction which will lead to an optimal second-order parameter uniform convergence. We present the results of numerical experiments for linear and semilinear reaction-diffusion systems to support the effectiveness of our preferred monitor function obtained from theoretical analysis. (C) 2014 Elsevier Inc. All rights reserved.

Determination of diffusion parameters and activation energy of diffusion in V3Si phase with A15 crystal structure

Diffusion such is the integrated diffusion coefficient of the phase, the tracer diffusion coefficient of species at different temperatures and the activation energy for diffusion, are determined in V3Si phase with A15 crystal structure. The tracer diffusion coefficient of Si Was found to be negligible compared to the tracer diffusion coefficient of V. The calculated diffusion parameters will help to validate the theoretical analysis of defect structure of the phase, which plays an important role in the superconductivity.

Molecular theory for the effects of specific solute-solvent interaction on the diffusion of a solute particle in a molecular liquid

An understanding of the effect of specific solute-solvent interactions on the diffusion of a solute probe is a long standing problem of physical chemistry. In this paper a microscopic treatment of this effect is presented. The theory takes into account the modification of the solvent structure around the solute due to this specific interaction between them. It is found that for strong, attractive interaction, there is an enhanced coupling between the solute and the solvent dynamic modes (in particular, the density mode), which leads to a significant increase in the friction on the solute. The diffusion coefficient of the solute is found to depend strongly and nonlinearly on the magnitude of the attractive interaction. An interesting observation is that specific solute-solvent interaction can induce a crossover from a sliplike to a sticklike diffusion. In the limit of strong attractive interaction, we recover a dynamic version of the solvent-berg picture. On the other hand, for repulsive interaction, the diffusion coefficient of the solute increases. These results are in qualitative agreement with recent experimental observations.

Stability of a Random Diffusion with Linear Drift

We consider a Linear system with Markovian switching which is perturbed by Gaussian type noise, If the linear system is mean square stable then we show that under certain conditions the perturbed system is also stable, We also shaw that under certain conditions the linear system with Markovian switching can be stabilized by such noisy perturbation.

Interdiffusion and Growth of the Superconductor Nb3Sn in Nb/Cu(Sn) Diffusion Couples

The confusion over the growth rate of the Nb3Sn superconductor compound following the bronze technique is addressed. Furthermore, a possible explanation for the corrugated structure of the product phase in the multifilamentary structure is discussed. Kirkendall marker experiments are conducted to study the relative mobilities of the species, which also explains the reason for finding pores in the product phase layer. The movement of the markers after interdiffusion reflects that Sn is the faster diffusing species. Furthermore, different concentrations of Sn in the bronze alloy are considered to study the effect of Sn content on the growth rate. Based on the parabolic growth constant at different temperatures, the activation energy for the growth is determined.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

A comment on the consistency of truncated nonlinear integral equation based theories of freezing

We report the results of two studies of aspects of the consistency of truncated nonlinear integral equation based theories of freezing: (i) We show that the self-consistent solutions to these nonlinear equations are unfortunately sensitive to the level of truncation. For the hard sphere system, if the Wertheim–Thiele representation of the pair direct correlation function is used, the inclusion of part but not all of the triplet direct correlation function contribution, as has been common, worsens the predictions considerably. We also show that the convergence of the solutions found, with respect to number of reciprocal lattice vectors kept in the Fourier expansion of the crystal singlet density, is slow. These conclusions imply great sensitivity to the quality of the pair direct correlation function employed in the theory. (ii) We show the direct correlation function based and the pair correlation function based theories of freezing can be cast into a form which requires solution of isomorphous nonlinear integral equations. However, in the pair correlation function theory the usual neglect of the influence of inhomogeneity of the density distribution on the pair correlation function is shown to be inconsistent to the lowest order in the change of density on freezing, and to lead to erroneous predictions. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

Natural convection from a vertical cone in a porous medium due to the combined effects of heat and mass diffusion with non-uniform wall temperature/concentration or heat/mass flux and suction/injection

An analysis has been carried out to study the non-Darcy natural convention flow of Newtonian fluids on a vertical cone embedded in a saturated porous medium with power-law variation of the wall temperature/concentration or heat/mass flux and suction/injection with the streamwise distance x. Both non-similar and self-similar solutions have been obtained. The effects of non-Darcy parameter, ratio of the buoyancy forces due to mass and heat diffusion, variation of wall temperature/concentration or heat/mass flux and suction/injection on the Nusselt and Sherwood numbers have been studied.