Non-linear transport

As an illustration, consider the isothermal, isobaric diffusional mixing of two elemental crystals, A and B, by a vacancy mechanism. Initially, A and B possess different vacancy concentrations Cy(A) and Cy(B). During interdiffusion, these concentrations have to change locally towards the new equilibrium values Cy(A,B), which depend on the local (A, B) composition. Vacancy relaxation will be slow if the external surfaces of the crystal, which act as the only sinks and sources, are far away. This is true for large samples. Although linear transport theory may apply for all structure elements, the (local) vacancy equilibrium is not fully established during the interdiffusion process. Consequently, the (local) transport coefficients (DA,DB), which are proportional to the vacancy concentration, are no longer functions of state (Le., dependent on composition only) but explicitly dependent on the diffusion time and the space coordinate. Non-linear transport equations are the result. 

In non-linear transport models, see e.g. [72], the chemical reaction rates depend on more parameters not only on concentrations and temperature but also on deformation rates, gradients of concentrations, etc. for a possible generalization of the presented procedure see [108, 164] and the end of this section. 

Electro-osmotic effects are of intrinsic interest also, since they give rise to a number of steady-state phenomena which can be conveniently studied experimentally. These afford examples where the power of non-equilibrium thermodynamics can be easily demonstrated. This is a unique phenomenon where non-linear transport processes have been experimentally observed and studied. 

The electro-osmotic data in the experimentally studied non-linear region are found to satisfy the following non-linear transport equation [12, 65]... 

As expected for a charge density wave (CDW) system, below the M-I transition, and for electrical fields high enough to depin the CDW, non-linear transport phenomena are observed in many a-phase compounds, namely those with M = Au and Pt [92—94]. In spite of the non-linear transport being a recent aspect of the research in this family of compounds, the experimental data already existing provides by far the strongest and more complete evidence for non-linear transport as due to a CDW motion in molecular conductors. The Au and Pt compounds with low transition temperatures at 8.2 and 12.2 K, have been more studied because their... 

Fig. 16. A. Plot of log iNa as a function of T 1 (°K) using the experimental values of the rate constants and the location of the binding sites in Eq. 4. The Gibbs free energy of activation is calculated from Eq. 3 the AS are taken to be zero, and the current is calculated by means of Eq. 4. The purpose is to demonstrate that multibarrier channel transport can be seen as single rate process with average values for the enthalpies of activation. Non-linearity of such a plot is then taken to arise form the dynamic nature of the channel.

The above mechanistic aspect of electron transport in electroactive polymer films has been an active and chemically rich research topic (13-18) in polymer coated electrodes. We have called (19) the process "redox conduction", since it is a non-ohmic form of electrical conductivity that is intrinsically different from that in metals or semiconductors. Some of the special characteristics of redox conductivity are non-linear current-voltage relations and a narrow band of conductivity centered around electrode potentials that yield the necessary mixture of oxidized and reduced states of the redox sites in the polymer (mixed valent form). Electron hopping in redox conductivity is obviously also peculiar to polymers whose sites comprise spatially localized electronic states. 

WASP/TOXIWASP/WASTOX. The Water Quality Analysis Simulation Program (WASP, 3)is a generalized finite-difference code designed to accept user-specified kinetic models as subroutines. It can be applied to one, two, and three-dimensional descriptions of water bodies, and process models can be structured to include linear and non-linear kinetics. Two versions of WASP designed specifically for synthetic organic chemicals exist at this time. TOXIWASP (54) was developed at the Athens Environmental Research Laboratory of U.S. E.P.A. WASTOX (55) was developed at HydroQual, with participation from the group responsible for WASP. Both codes include process models for hydrolysis, biolysis, oxidations, volatilization, and photolysis. Both treat sorption/desorption as local equilibria. These codes allow the user to specify either constant or time-variable transport and reaction processes. 

The last term on the right-hand side of (1.28) is the chemical source term. As will be seen in Chapter 5, the chemical source term is often a complex, non-linear function of the scalar fields , and thus solutions to (1.28) are very different than those for the z nm-scalar transport equation wherein S is null. 

All non-linear terms involving spatial gradients require transported PDF closures. Examples of such terms are viscous dissipation, pressure fluctuations, and scalar dissipation. 

We shall see that transported PDF closures forthe velocity field are usually linear in V. Thus (/ D) will depend only on the first two moments of U. In general, non-linear velocity models could be formulated, in which case arbitrary moments of U would appear in the Reynolds-stress transport equation. 

Because the conditional scalar Laplacian is approximated in the FP model by a non-linear diffusion process (6.91), (6.145) will not agree exactly with CMC. Nevertheless, since transported PDF methods can be easily extended to inhomogeneous flows,113 which are problematic for the CMC, the FP model offers distinct advantages. 

In a transported PDF simulation, the chemical source term, (6.249), is integrated over and over again with each new set of initial conditions. For fixed inlet flow conditions, it is often the case that, for most of the time, the initial conditions that occur in a particular simulation occupy only a small sub-volume of composition space. This is especially true with fast chemical kinetics, where many of the reactions attain a quasi-steady state within the small time step At. Since solving the stiff ODE system is computationally expensive, this observation suggests that it would be more efficient first to solve the chemical source term for a set of representative initial conditions in composition space,156 and then to store the results in a pre-computed chemical lookup table. This operation can be described mathematically by a non-linear reaction map ... 

An example of a smart tabulation method is the intrinsic, low-dimensional manifold (ILDM) approach (Maas and Pope 1992). This method attempts to reduce the number of dimensions that must be tabulated by projecting the composition vectors onto the nonlinear manifold defined by the slowest chemical time scales.162 In combusting systems far from extinction, the number of slow chemical time scales is typically very small (i.e, one to three). Thus the resulting non-linear slow manifold ILDM will be low-dimensional (see Fig. 6.7), and can be accurately tabulated. However, because the ILDM is non-linear, it is usually difficult to find and to parameterize for a detailed kinetic scheme (especially if the number of slow dimensions is greater than three ). In addition, the shape, location in composition space, and dimension of the ILDM will depend on the inlet flow conditions (i.e., temperature, pressure, species concentrations, etc.). Since the time and computational effort required to construct an ILDM is relatively large, the ILDM approach has yet to find widespread use in transported PDF simulations outside combustion. 

The initial concentration distribution is therefore simply translated at the velocity of the liquid steady flow and full equilibrium between the liquid and its matrix require that the amount of element transported by the concentration wave is constant. In more realistic cases, either the flow is non-steady due to abrupt changes in fluid advection rate or porosity, or solid-liquid equilibrium is not achieved. These cases may lead to non-linear terms in the chromatographic equation (9.4.35) and unstable behavior. The rather complicated theory of these processes is beyond the scope of the present book. 

At sufficiently high concentrations of the transported solute particles, the surface coverage becomes important and non-linear laws for the rate of adsorption should be used. 

Both a simplified continuous and discrete model, describing the behaviour of single component mass transport in chromatographic columns with non-linear distribution isotherm, were developed and simulated by Smit et al. Studies of more complex but still relatively simple (multicomponent) transport models have been published (see e.g. ... 

---