Diffusion coefficients systems

Figure C2.3.8. Self-diffusion coefficients at 45°C for AOT ( ), water ( ) and decane ( ) in ternary AOT, brine (0.6% aqueous NaCl) and decane microemulsion system as a function of composition, a. This compositional parameter, a, is tire weight fraction of decane relative to decane and brine. Reproduced by pennission from figure 3 of [46].

The first thing to notice about these results is that the influence of the micropores reduces the effective diffusion coefficient below the value of the bulk diffusion coefficient for the macropore system. This is also clear in general from the forms of equations (10.44) and (10.48). As increases from zero, corresponding to the introduction of micropores, the variance of the response pulse Increases, and this corresponds to a reduction in the effective diffusion coefficient. The second important point is that the influence of the micropores on the results is quite small-Indeed it seems unlikely that measurements of this type will be able to realize their promise to provide information about diffusion in dead-end pores. 

The shear viscosity is a tensor quantity, with components T] y, t],cz, T)yx> Vyz> Vzx> Vzy If property of the whole sample rather than of individual atoms and so cannot be calculat< with the same accuracy as the self-diffusion coefficient. For a homogeneous fluid the cor ponents of the shear viscosity should all be equal and so the statistical error can be reducf by averaging over the six components. An estimate of the precision of the calculation c then be determined by evaluating the standard deviation of these components from tl average. Unfortunately, Equation (7.89) cannot be directly used in periodic systems, evi if the positions have been unfolded, because the unfolded distance between two particl may not correspond to the distance of the minimum image that is used to calculate the fore For this reason alternative approaches are required. 

Monte Carlo simulations require less computer time to execute each iteration than a molecular dynamics simulation on the same system. However, Monte Carlo simulations are more limited in that they cannot yield time-dependent information, such as diffusion coefficients or viscosity. As with molecular dynamics, constant NVT simulations are most common, but constant NPT simulations are possible using a coordinate scaling step. Calculations that are not constant N can be constructed by including probabilities for particle creation and annihilation. These calculations present technical difficulties due to having very low probabilities for creation and annihilation, thus requiring very large collections of molecules and long simulation times. 

In this expression, called Pick s first law, the proportionality constant D is the diffusion coefficient of the solute. Since J = (l/A)(dQ/dt) and c = Q/V, where Q signifies the quantity of solute in unspecified units, it follows that D has the units length time", or m sec in the SI system. The minus sign in Eq. (9.69)... 

Multicomponent Diffusion. In multicomponent systems, the binary diffusion coefficient has to be replaced by an effective or mean diffusivity Although its rigorous computation from the binary coefficients is difficult, it may be estimated by one of several methods (27—29). Any degree of counterdiffusion, including the two special cases "equimolar counterdiffusion" and "no counterdiffusion" treated above, may arise in multicomponent gas absorption. The influence of bulk flow of material through the films is corrected for by the film factor concept (28). It is based on a slightly different form of equation 13 ... 

The main conclusion to be drawn from these studies is that for most practical purposes the linear rate model provides an adequate approximation and the use of the more cumbersome and computationally time consuming diffusing models is generally not necessary. The Glueckauf approximation provides the required estimate of the effective mass transfer coefficient for a diffusion controlled system. More detailed analysis shows that when more than one mass transfer resistance is significant the overall rate coefficient may be estimated simply from the sum of the resistances (7) ... 

Most distillation systems ia commercial columns have Murphree plate efficiencies of 70% or higher. Lower efficiencies are found under system conditions of a high slope of the equiHbrium curve (Fig. lb), of high Hquid viscosity, and of large molecules having characteristically low diffusion coefficients. FiaaHy, most experimental efficiencies have been for biaary systems where by definition the efficiency of one component is equal to that of the other component. For multicomponent systems it is possible for each component to have a different efficiency. Practice has been to use a pseudo-biaary approach involving the two key components. However, a theory for multicomponent efficiency prediction has been developed (66,67) and is amenable to computational analysis. 

These three terms represent contributions to the flux from migration, diffusion, and convection, respectively. The bulk fluid velocity is determined from the equations of motion. Equation 25, with the convection term neglected, is frequently referred to as the Nemst-Planck equation. In systems containing charged species, ions experience a force from the electric field. This effect is called migration. The charge number of the ion is Eis Faraday s constant, is the ionic mobiUty, and O is the electric potential. The ionic mobiUty and the diffusion coefficient are related ... 

For prediction of gas phase diffusion coefficients in multicomponent hydi ocarbon/nonKydi ocai bon gas systems, the method of Wilke shown in Eq. (2-154) is used. 

Fuller-Schettler-Giddings The parameters and constants for this correlation were determined by regression analysis of 340 experimental diffusion coefficient values of 153 binary systems. Values of X Vj used in this equation are in Table 5-16. 

Cussler studied diffusion in concentrated associating systems and has shown that, in associating systems, it is the size of diffusing clusters rather than diffusing solutes that controls diffusion. is a reference diffusion coefficient discussed hereafter is the activity of component A and iC is a constant. By assuming that could be predicted by Eq. (5-223) with P = 1, iC was found to be equal to 0.5 based on five binaiy systems and vahdated with a sixth binaiy mixture. The limitations of Eq. (5-225) using and K defined previously have not been explored, so caution is warranted. Gurkan showed that K shoiild actually be closer to 0.3 (rather than 0.5) and discussed the overall results. 

At any plane in a Raoultian alloy system parallel to die original interface, the so-called chemical diffusion coefficient Dchem. which determines the flux of atoms at any given point, and is usually a function of the local composition so that according to Darken (1948), Dchem is given by... 

An experimental technique for die determination of Dchem in a binary alloy system in which die diffusion coefficient is a function of composition was originally developed by Matano (1932), based on a mathematical development... 

In part the parabolic law may also apply to multilayer oxide systems where the cation diffusion coefficient is much higher in the lower oxide tlran in the higher oxide, which, growing as a thin layer, undergoes plastic deformation at high temperatures, thus retaining the overall oxide layer as impervious to enuy of tire gas. 

An equation has been developed to relate the diffusion coefficients of two systems at the same temperature, as follows ... 

Bailey, R. G. and Chen. H. T., Predicting Diffusion Coefficients In Binary Gas Systems," Chemical Engineering, March 17, 1975, p. 86. 

The diffusivity coefficient for a liquid-liquid system can be estimated from the Wilke correlation ... 

Here the index 0 at Z) refers to the diffusion coefficient in a system at rest (no bias), and Vq denotes the zero approximation i.e., this relation is applicable for slow drift in dilute medium only. 

Equations (4-5) and (4-7) are alternative expressions for the estimation of the diffusion-limited rate constant, but these equations are not equivalent, because Eq. (4-7) includes the assumption that the Stokes-Einstein equation is applicable. Olea and Thomas" measured the kinetics of quenching of pyrene fluorescence in several solvents and also measured diffusion coefficients. The diffusion coefficients did not vary as t) [as predicted by Eq. (4-6)], but roughly as Tf. Thus Eq. (4-7) is not valid, in this system, whereas Eq. (4-5), used with the experimentally measured diffusion coefficients, gave reasonable agreement with measured rate constants. 

---