Molecular diffusion models

Parameter for molecular diffusion model In a moving zone, equivalent to the reciprocal of Peclet number, dispersion number Reynolds number, Re Prandtle number, Pr Schmidt number Sc... 

Section IA summarizes the molecular model of diffusion of Pace and Datyner (1 2) which proposes that the diffusion of gases in a polymeric matrix is determined by the cooperative main-chain motions of the polymer. In Section IB we report carbon-13 nmr relaxation measurement which show that the diffusion of gases in poly(vinyl chloride) (PVC) - tricresyl phosphate (TCP) systems is controlled by the cooperative motions of the polymer chains. The correlation of the phenomenological diffusion coefficients with the cooperative main-chain motions of the polymer provides an experimental verification for the molecular diffusion model. 

The models above may be useful for predicting mass fluxes in MD however, each of these models has its limitations. The Knudsen and Poiseuille model require knowledge of r, 8, and e, which in general can be estimated by applying the models to experimental gas fluxes through the given membrane. The molecular diffusion model is inadequate at low-partial pressures of air, as it predicts infinite flux since, in totally deaerated membrane 7 tends to zero. 

One apparent advantage of this model is that it predicts a dependence of flux on the one-half power of the diffusion coefficient, which is closer to reality than the first power for the molecular diffusion model. And it is more realistic, in that eddy-type packets are considered. However, it does not have an application advantage, since there is little physical basis for estimating the exposure... 

Film Theory. Many theories have been put forth to explain and correlate experimentally measured mass transfer coefficients. The classical model has been the film theory (13,26) that proposes to approximate the real situation at the interface by hypothetical "effective" gas and Hquid films. The fluid is assumed to be essentially stagnant within these effective films making a sharp change to totally turbulent flow where the film is in contact with the bulk of the fluid. As a result, mass is transferred through the effective films only by steady-state molecular diffusion and it is possible to compute the concentration profile through the films by integrating Fick s law ... 

Fig. 9. Uptake curves for N2 in two samples of carbon molecular sieve showing conformity with diffusion model (eq. 24) for sample 1 (A), and with surface resistance model (eq. 26) for example 2 (0)j LDF = linear driving force. Data from ref. 18.

Dispersion Model An impulse input to a stream flowing through a vessel may spread axially because of a combination of molecular diffusion and eddy currents that together are called dispersion. Mathematically, the process can be represented by Fick s equation with a dispersion coefficient replacing the diffusion coefficient. The dispersion coefficient is associated with a linear dimension L and a linear velocity in the Peclet number, Pe = uL/D. In plug flow, = 0 and Pe oq and in a CSTR, oa and Pe = 0. 

Models for description of liquids should provide us with an understanding of the dynamic behavior of the molecules, and thus of the routes of chemical reactions in the liquids. While it is often relatively easy to describe the molecular structure and dynamics of the gaseous or the solid state, this is not true for the liquid state. Molecules in liquids can perform vibrations, rotations, and translations. A successful model often used for the description of molecular rotational processes in liquids is the rotational diffusion model, in which it is assumed that the molecules rotate by small angular steps about the molecular rotation axes. One quantity to describe the rotational speed of molecules is the reorientational correlation time T, which is a measure for the average time elapsed when a molecule has rotated through an angle of the order of 1 radian, or approximately 60°. It is indirectly proportional to the velocity of rotational motion. 

Glaser and Litt (G4) have proposed, in an extension of the above study, a model for gas-liquid flow through a b d of porous particles. The bed is assumed to consist of two basic structures which influence the fluid flow patterns (1) Void channels external to the packing, with which are associated dead-ended pockets that can hold stagnant pools of liquid and (2) pore channels and pockets, i.e., continuous and dead-ended pockets in the interior of the particles. On this basis, a theoretical model of liquid-phase dispersion in mixed-phase flow is developed. The model uses three bed parameters for the description of axial dispersion (1) Dispersion due to the mixing of streams from various channels of different residence times (2) dispersion from axial diffusion in the void channels and (3) dispersion from diffusion into the pores. The model is not applicable to turbulent flow nor to such low flow rates that molecular diffusion is comparable to Taylor diffusion. The latter region is unlikely to be of practical interest. The model predicts that the reciprocal Peclet number should be directly proportional to nominal liquid velocity, a prediction that has been confirmed by a few determinations of residence-time distribution for a wax desulfurization pilot reactor of 1-in. diameter packed with 10-14 mesh particles. 

In the film-penetration model (H19), it is assumed that the reactant A penetrates through the surface element by one-dimensional unsteady-state molecular diffusion. Convective transport is assumed to be insignificant. The diffusing stream of the reactant A is depleted along the path of diffusion by its reversible reaction with the reactant B, which is an existing component of the liquid surface element. If such a reaction can be represented as... 

Marchello and Toor (M2) proposed a mixing model for transfer near a boundary which assumes that localized mixing occurs rather than gross displacement of the fluid elements. This model can be said to be a modified penetration-type model. Kishinevsky (K6-K8) assumed a surface-renewal mechanism with eddy diffusion rather than molecular diffusion controlling the transfer at the interface. 

When a polymer film is exposed to a gas or vapour at one side and to vacuum or low pressure at the other, the mechanism generally accepted for the penetrant transport is an activated solution-diffusion model. The gas dissolved in the film surface diffuses through the film by a series of activated steps and evaporates at the lower pressure side. It is clear that both solubility and diffusivity are involved and that the polymer molecular and morphological features will affect the penetrant transport behaviour. Some of the chemical and morphological modification that have been observed for some epoxy-water systems to induce changes of the solubility and diffusivity will be briefly reviewed. 

Kishinev ski/23 has developed a model for mass transfer across an interface in which molecular diffusion is assumed to play no part. In this, fresh material is continuously brought to the interface as a result of turbulence within the fluid and, after exposure to the second phase, the fluid element attains equilibrium with it and then becomes mixed again with the bulk of the phase. The model thus presupposes surface renewal without penetration by diffusion and therefore the effect of diffusivity should not be important. No reliable experimental results are available to test the theory adequately. 

At lower Reynolds numbers, the axial velocity profile will not be flat and it might seem that another correction must be added to Equation (9.14). It turns out, however, that Equation (9.14) remains a good model for real turbulent reactors (and even some laminar ones) given suitable values for D. The model lumps the combined effects of fluctuating velocity components, nonflat velocity profiles, and molecular diffusion into the single parameter D. 

At a close level of scrutiny, real systems behave differently than predicted by the axial dispersion model but the model is useful for many purposes. Values for Pe can be determined experimentally using transient experiments with nonreac-tive tracers. See Chapter 15. A correlation for D that combines experimental and theoretical results is shown in Figure 9.6. The dimensionless number, udt/D, depends on the Reynolds number and on molecular diffusivity as measured by the Schmidt number, Sc = but the dependence on Sc is weak for... 

In the absence of diffusion, all hydrodynamic models show infinite variances. This is a consequence of the zero-slip condition of hydrodynamics that forces Vz = 0 at the walls of a vessel. In real systems, molecular diffusion will ultimately remove molecules from the stagnant regions near walls. For real systems, W t) will asymptotically approach an exponential distribution and will have finite moments of all orders. However, molecular diffusivities are low for liquids, and may be large indeed. This fact suggests the general inappropriateness of using to characterize the residence time distribution in a laminar flow system. Turbulent flow is less of a problem due to eddy diffusion that typically results in an exponentially decreasing tail at fairly low multiples of the mean residence time. 

The devolatilization of a component in an internal mixer can be described by a model based on the penetration theory [27,28]. The main characteristic of this model is the separation of the bulk of material into two parts A layer periodically wiped onto the wall of the mixing chamber, and a pool of material rotating in front of the rotor flights, as shown in Figure 29.15. This flow pattern results in a constant exposure time of the interface between the material and the vapor phase in the void space of the internal mixer. Devolatilization occurs according to two different mechanisms Molecular diffusion between the fluid elements in the surface layer of the wall film and the pool, and mass transport between the rubber phase and the vapor phase due to evaporation of the volatile component. As the diffusion rate of a liquid or a gas in a polymeric matrix is rather low, the main contribution to devolatilization is based on the mass transport between the surface layer of the polymeric material and the vapor phase. 

This simple relaxation theory becomes invalid, however, if motional anisotropy, or internal motions, or both, are involved. Then, the rotational correlation-time in Eq. 30 is an effective correlation-time, containing contributions from reorientation about the principal axes of the rotational-diffusion tensor. In order to separate these contributions, a physical model to describe the manner by which a molecule tumbles is required. Complete expressions for intramolecular, dipolar relaxation-rates for the three classes of spherical, axially symmetric, and asymmetric top molecules have been evaluated by Werbelow and Grant, in order to incorporate into the relaxation theory the appropriate rotational-diffusion model developed by Woess-ner. Methyl internal motion has been treated in a few instances, by using the equations of Woessner and coworkers to describe internal rotation superimposed on the overall, molecular tumbling. Nevertheless, if motional anisotropy is present, it is wiser not to attempt a quantitative determination of interproton distances from measured, proton relaxation-rates, although semiquantitative conclusions are probably justified by neglecting motional anisotropy, as will be seen in the following Section. 

Axial and radial dispersion or non-ideal flow in tubular reactors is usually characterised by analogy to molecular diffusion, in which the molecular diffusivity is replaced by eddy dispersion coefficients, characterising both radial and longitudinal dispersion effects. In this text, however, the discussion will be limited to that of tubular reactors with axial dispersion only. Otherwise the model equations become too complicated and beyond the capability of a simple digital simulation language. 

The modeling of mass transport from the bulk fluid to the interface in capillary flow typically applies an empirical mass transfer coefficient approach. The mass transfer coefficient is defined in terms of the flux and driving force J = kc(cbuik-c). For non-reactive steady state laminar flow in a square conduit with constant molecular diffusion D, the mass balance in the fluid takes the form... 

Steps (i) and (ii) are controlled by molecular diffusion. Higher operating temperatures can improve the kinetics of mass transfer in all three steps. Vandenburg et al. [37] described the kinetics of PFE extraction using the hot-ball model [286] derived for SFE extractions. 

The basic biofilm model149,150 idealizes a biofilm as a homogeneous matrix of bacteria and the extracellular polymers that bind the bacteria together and to the surface. A Monod equation describes substrate use molecular diffusion within the biofilm is described by Fick s second law and mass transfer from the solution to the biofilm surface is modeled with a solute-diffusion layer. Six kinetic parameters (several of which can be estimated from theoretical considerations and others of which must be derived empirically) and the biofilm thickness must be known to calculate the movement of substrate into the biofilm. 

The flux vector accounts for mass transport by both convection (i.e., blood flow, interstitial fluid flow) and conduction (i.e., molecular diffusion), whereas S describes membrane transport between adjacent compartments and irreversible elimination processes. For the three-subcompartment organ model presented in Figure 2, with concentration both space- and time-dependent, the conservation equations are... 

---