Transport equations formulation

The TDE solute module is formulated with one equation describing pollutant mass balance of the species in a representative soil volume dV = dxdydz. The solute module is frequently known as the dispersive, convective differential mass transport equation, in porous media, because of the wide employment of this equation, that may also contain an adsorptive, a decay and a source or sink term. The one dimensional formulation of the module is ... 

Barkey, Tobias and Muller formulated the stability analysis for deposition from well-supported solution in the Tafel regime at constant current [48], They used dilute-solution theory to solve the transport equations in a Nernst diffusion layer of thickness S. The concentration and electrostatic potential are given in this approximation... 

Besides the resuspension of particles, the perfect sink model also neglects the effect of deposited particles on incoming particles. To overcome these limitations, recent models [72, 97-99] assume that particles accumulate within a thin adsorption layer adjacent to the collector surface, and replace the perfect sink conditions with the boundary condition that particles cannot penetrate the collector. General continuity equations are formulated both for the mobile phase and for the immobilized particles in which the immobilization reaction term is decomposed in an accumulation and a removal term, respectively. Through such equations, one can keep track of the particles which arrive at the primary minimum distance and account for their normal and tangential motion. These equations were solved both approximately, and by numerical integration of the governing non-stationary transport equations. 

Note that although the density is constant, we have included it in the transport equations to be consistent with the formulation used in commercial CFD codes. 

The simplest aggregation and breakage models can be formulated in terms of the NDF (o), which uses volume as the independent variable.6 The microscopic transport equation for the NDF has the form (Wang et al., 2005a,b)... 

Section II deals with the general formalism of Prigogine and his co-workers. Starting from the Liouville equation, we derive an exact transport equation for the one-particle distribution function of an arbitrary fluid subject to a weak external field. This equation is valid in the so-called "thermodynamic limit , i.e. when the number of particles N —> oo, the volume of the system 2-> oo, with Nj 2 = C finite. As a by-product, we obtain very easily a formulation for the equilibrium pair distribution function of the fluid as well as a general expression for the conductivity tensor. 

The first term on the right-hand side of (2.61) is the spectral transfer function, and involves two-point correlations between three components of the velocity vector (see McComb (1990) for the exact form). The spectral transfer function is thus unclosed, and models must be formulated in order to proceed in finding solutions to (2.61). However, some useful properties of T (k, t) can be deduced from the spectral transport equation. For example, integrating (2.61) over all wavenumbers yields the transport equation for the turbulent kinetic energy ... 

In (5.297), the interpolation parameter is defined separately for each component. Note, however, that unlike the earlier examples, there is no guarantee that the interpolation parameters will be bounded between zero and one. For example, the equilibrium concentration of intermediate species may be negligible despite the fact that these species can be abundant in flows dominated by finite-rate chemistry. Thus, although (5.297) provides a convenient closure for the chemical source term, it is by no means guaranteed to produce accurate predictions A more reliable method for determining the conditional moments is the formulation of a transport equation that depends explicitly on turbulent transport and chemical reactions. We will look at this method for both homogeneous and inhomogeneous flows below. 

It is now necessary to formulate transport equations for p and (s). However, by making the following substitutions, the same transport equations as are used in the multienvironment PDF model can be used for the multi-environment LES model ... 

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. 

The remaining challenge is then to formulate and solve transport equations for the mapping functions g(z x, 0 (Gao and O Brien 1991 Pope 1991b). Note that if g(z x, 0 is known, then the FP model can be used to describe Z(0, and 0(0 will follow from (6.121). Since the PDF of Z is stationary and homogeneous, the FP model needed to describe it will be particularly simple. With the mapping closure, the difficulties associated with the chemical source terms are thus shifted to the model for g(z x, 0. 

Note that the RANS formulation used in (B.44) and (B.45) can easily be extended to the LES, as outlined in Section 5.10. Moreover, by following the same steps as outlined above, DQMOM can be used with the joint velocity, composition PDF transport equation. Finally, the reader can observe that the same methodology is applicable to more general distribution functions than probability density functions. Indeed, DQMOM can be applied to general population balance equations such as those used to describe multi-phase flows. 

Divisek et al. presented a similar two-phase, two-dimensional model of DMFC. Two-phase flow and capillary effects in backing layers were considered using a quantitatively different but qualitatively similar function of capillary pressure vs liquid saturation. In practice, this capillary pressure function must be experimentally obtained for realistic DMFC backing materials in a methanol solution. Note that methanol in the anode solution significantly alters the interfacial tension characteristics. In addition, Divisek et al. developed detailed, multistep reaction models for both ORR and methanol oxidation as well as used the Stefan—Maxwell formulation for gas diffusion. Murgia et al. described a one-dimensional, two-phase, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymer membrane for a liquid-feed DMFC. 

A variety of specific mathematical formulations of the CTRW approach have been considered to date, and network models have also been applied (Bijeljic and Blunt 2006). A key result in development of the CTRW approach is a transport equation that represents a strong generalization of the advection-dispersion equation. As shown by Berkowitz et al. (2006), an extremely broad range of transport patterns can be described with the (ensemble-averaged) equation... 

The set of partial differential equations developed for the simultaneous transfer of moisture, hear, and reactive chemicals under saturated/unsaturated soil conditions has been solved by the Galerkin finite element method. The chemical transport equations are formulated in terms of the total analytical concentration of each component species, and can be solved sequentially (Wu and Chieng, 1995). 

The transport equations for laminar motion can be formulated, in general, easily and difficulties may lie only in their solution. On the other hand, for turbulent motion the formulation of the basic equations for the time-averaged local quantities constitutes a major physical difficulty. In recent developments, one considers that turbulence (chaos) is predictable from the time-dependent transport equations. However, this point of view is beyond the scope of the present treatment. For the present, some simple procedures based on physical models and scaling will be employed to obtain useful results concerning turbulent heat or mass transfer. 

Finally, let us briefly point out some essential features of the stability analysis for a more general transport problem. It can be exemplified by the moving a//9 phase boundary in the ternary system of Figure 11-12. Referring to Figure 11-7 and Eqn. (11.10), it was a single independent (vacancy) flux that caused the motion of the boundary. In the case of two or more independent components, we have to formulate the transport equation (Fick s second law) for each component, both in the a- and /9-phase. Each of the fluxes jf couples at the boundary b with jf, i = A,B,... (see, for example, Eqn. (11.2)). Furthermore, in the bulk, the fluxes are also coupled (e.g., by electroneutrality or site conservation). 

In this section, we consider the influence of inhomogeneous stress on matter transport. Tb illustrate the problem, let us formulate a simple transport equation diffusion of an interstitial component i in an otherwise immobile solid (e.g., H in Pd). Furthermore, we neglect cross effects. For an electrically neutral species i (i.e., H) we then have... 

Qualitatively, to model the chemically reacting plasma, neutral-species chemistry and transport and the electrical or plasma physics aspects of the discharge must be considered. The chemical and physical properties included in the model of the discharge must at least reflect the current understanding of the most important processes. The conventional view of discharge behavior has been described previously, but will now be reviewed briefly to set the stage for equation formulation. 

This was the idea behind concept (14b) in Ref.4). The corresponding formulation of Eq. (15i) of Ref.4), however, was unhappily chosen Eq. (15f)ofthis paper should have been used. If the kinetics of separation were explicitly introduced into the transport Equation of PDC (instead of the implicit concept (I7a-c) of the flow-equilibrium), an integrodifferential equation more complicated than (41 a-b) would be obtained, which could hardly be solved analytically... 

The elution and migration effects are also found in PDC, but not the compression effect, since no precipitant and temperature gradients are applied. Unlike BWF, no simple transport Equation is applied in PDC because the replacement of the linear Eq. (3b) by the non-linear (17c) would lead to an integrodifferential equation similar to (41 a-b), but more complicated, if some explicit formulation were used instead of the implicit one, based on a flow-equilibrium and on a perturbation calculus, applied to an integrated transport Equation... 

Formulation of the transport equations in the terms of the aforementioned mechanism leads to a set of integro-differential and differential equations. The numerical solution of these equations is associated with serious computational difficulties (57). 

Extension of the equilibrium model to column or field conditions requires coupling the ion-exchange equations with the transport equations for the 5 aqueous species (Eq. 1). To accomplish this coupling, we have adopted the split-operator approach (e.g., Miller and Rabideau, 1993), which provides considerable flexibility in adjusting the sorption submodel. In addition to the above conceptual model, we are pursuing more complex formulations that couple cation exchange with pore diffusion, surface diffusion, or combined pore/surface diffusion (e.g., Robinson et al., 1994 DePaoli and Perona, 1996 Ma et al., 1996). However, the currently available data are inadequate to parameterize such models, and the need for a kinetic formulation for the low-flow conditions expected for sorbing barriers has not been established. These issues will be addressed in a future publication. 

The proposed model consists of a biphasic mechanical description of the tissue engineered construct. The resulting fluid velocity and displacement fields are used for evaluating solute transport. Solute concentrations determine biosynthetic behavior. A finite deformation biphasic displacement-velocity-pressure (u-v-p) formulation is implemented [12, 7], Compared to the more standard u-p element the mixed treatment of the Darcy problem enables an increased accuracy for the fluid velocity field which is of primary interest here. The system to be solved increases however considerably and for multidimensional flow the use of either stabilized methods or Raviart-Thomas type elements is required [15, 10]. To model solute transport the input features of a standard convection-diffusion element for compressible flows are employed [20], For flexibility (non-linear) solute uptake is included using Strang operator splitting, decoupling the transport equations [9],... 

It is worth noting at this point that the various scientific theories that quantitatively and mathematically formulate natural phenomena are in fact mathematical models of nature. Such, for example, are the kinetic theory of gases and rubber elasticity, Bohr s atomic model, molecular theories of polymer solutions, and even the equations of transport phenomena cited earlier in this chapter. Not unlike the engineering mathematical models, they contain simplifying assumptions. For example, the transport equations involve the assumption that matter can be viewed as a continuum and that even in fast, irreversible processes, local equilibrium can be achieved. The paramount difference between a mathematical model of a natural process and that of an engineering system is the required level of accuracy and, of course, the generality of the phenomena involved. 

The derivation of the transport equations starts with the formulation of the entropy production rate. A differential change of the entropy of the isolated system dSsys is... 

The two-dimensional Navier-Stokes equation is solved in stream function-vorticity formulation, as reported variously in Sengupta et al. (2001, 2003), Sengupta Dipankar (2005). Brinckman Walker (2001) also simulated the burst sequence of turbulent boundary layer excited by streamwise vortices (in X- direction) using the same formulation for which a stream function was defined in the y — z) -plane only. To resolve various small scale events inside the shear layer, the vorticity transport equation (VTE) and the stream function equation (SFE) are solved in the transformed — rj) —... 

---