Reaction-diffusion-advection equation

The problem of reacting transport is different due to the presence of the production term/ (0, ..., 0N) [see Eq. (1)] that makes it to be nonlinear. Here we consider the simplest nontrivial case of Eq. (1) a unique scalar field 0(x, t) evolving according to the advection-reaction-diffusion equation... 

We consider the case of a unique scalar reactive field 0(jc, t). This model is appropriate in aqueous autocatalytic premixed reactions, as well as in gaseous combustion with a large flow intensity but low value of gas expansion across the flame [7]. The field 0 evolves according to the advection-reaction-diffusion (ARD) equation ... 

If the solute imdergoes any chemical changes, a reaction term must be added to Eq. 12.4. In the absence of specific rate law information, diagenetic reactions are generally assumed to be first-order with respect to the solute concentration. Thus, the one-dimensional advection-diffusion equation far a nonconservative solute is given by... 

The advection—diffusion equation with a source term can be solved by CFD algorithms in general. Patankar provided an excellent introduction to numerical fluid flow and heat transfer. Oran and Boris discussed numerical solutions of diffusion—convection problems with chemical reactions. Since fuel cells feature an aspect ratio of the order of 100, 0(100), the upwind scheme for the flow-field solution is applicable and proves to be very effective. Unstructured meshes are commonly employed in commercial CFD codes. 

Transport Processes and Gauss Theorem One-Dimensional Diffusion/Advection/Reaction Equation Box 22.1 One-Dimensional Diffusion/Advection/Reaction Equation at Steady-State... 

Table 22.2 Solution of Diffusion/Advection/Reaction Equation at Steady-State A, and A2 of Eq. 22-9 for Different Types of Boundary Conditions...

The one-dimensional diffusion/advection/reaction equation at steady-state is (22-7) ... 

Situations in which either Da or Pe are much larger or much smaller than 1 indicate that in the diffusion-advection-reaction equation some of these processes are dominant while others can be disregarded. Figure 22.3 gives an overview of such cases. A first distinction is made according to the size of Da ... 

The model (Fig. 23.6) consists of three compartments, (a) the surface mixed water layer (SMWL) or epilimnion, (b) the remaining open water column (OP), and (c) the surface mixed sediment layer (SMSL). SMWL and OP are assumed to be completely mixed their mass balance equations correspond to the expressions derived in Box 23.1, although the different terms are not necessarily linear. The open water column is modeled as a spatially continuous system described by a diffusion/advection/ reaction... 

It was evident in the earliest profiles obtained for dissolved °Th that this prediction does not hold, in that concentrations of dissolved °Th were observed to increase with depth in a nearly linear fashion (Bacon and Anderson, 1982 Nozaki et al., 1981). These results were interpreted to indicate that sorption processes are a reversible reaction in which particulate thorium may be released back into solution, as represented by k-i in Figure 1 (Bacon and Anderson, 1982 Nozaki et al., 1981). Neglecting advection and diffusion, mass balance equations for dissolved and particulate thorium can be written, respectively, as... 

Equation 2.20 is the advection-dispersion (AD) equation. In the petroleum literature, the term convection-diffusion (CD) equation is used, or simply diffusion equation (Brigham, 1974). When a reaction term is included, the term advection-reaction-dispersion (ARD) equation is used elsewhere. When the adsorption term is expressed as a reaction term, the ARD equation is as discussed later in Section 2.4. Several solutions of Eq. 2.20 have been presented in the literature, depending on the boundary conditions imposed. In general, they are various combinations of the error function. When the porous medium is long compared with the length of the mixed zone, they all give virtually identical results. 

A natural response to the limitations of both geochemical equilibrium models and the solute transport models (see 10.3 for a discussion) is to couple the two. Over the last two decades, a number of models that couple advective-dispersive-diffusive transport with fully speciated chemical reactions have been developed (see reviews by Engesgaard and Christensen, 1988 Grove and Stollenwerk, 1987 Mangold and Tsang, 1991). In the coupled models, the solute transport and chemical equilibrium equations are simultaneously evaluated. 

In a reaction-diffusion-advection formulation, the model considers a one-dimensional slice transverse to concentration filaments, and models its evolution by an equation of the type... 

Consider the binary chemical reaction A + B —> C. The reaction-diffusion-advection equations read, in the case of equal diffusion co-... 

In the case of the fast binary reaction we could eliminate the reaction term from the reaction-diffusion-advection equation. But in general this is not possible. In this chapter we consider another class of chemical and biological activity for which some explicit analysis is still feasible. We consider the case in which the local-reaction dynamics has a unique stable steady state at every point in space. If this steady state concentration was the same everywhere, then it would be a trivial spatially uniform solution of the full reaction-diffusion-advection problem. However, when the local chemical equilibrium is not uniform in space, due to an imposed external inhomogeneity, the competition between the chemical and transport dynamics may lead to a complex spatial structure of the concentration field. As we will see in this chapter, for this class of chemical or biological systems the dominant processes that determine the main characteristics of the solutions are the advection and the reaction dynamics, while diffusion does not play a major role in the large Peclet number limit considered here. Thus diffusion can be neglected in a first approximation. 

In the previous chapters we discussed various aspects of chemical and biological activity in fluid flows presenting certain classes of dynamical behavior that can be described by reaction-diffusion-advection equations and analyzed using dynamical systems approaches. However, there are many research areas of chemical or biological processes taking place in fluid environments that were not covered in the previous chapters. Here we briefly discuss some of these areas and point the reader to the relevant literature for further reading. Apart from classical well-studied topics, here we also focus on more recent developments and active areas of current research. 

The species-B balance equation includes advective transport, Fickian diffusion, and depletion by chemical reaction. The binary diffusion coefficient D represents downstream diffusion of reactant species-B relative to upstream diffusion of product species-C. The expression for Ys, the surface mass fraction of B (gas side), is obtained from a species balance at the surface on B which includes advective transport of pure B to the interface on the condensed phase side and both advective and diffusive transport of B away from the surface on the gas side. The downstream condition K(oo)=0 represents the assumption of complete conversion... 

We use Model C, given by (5.27), for the mean-field equation for p x,t) with z) = (p(t)w z). The standard diffusion approximation of (5.27), i.e., taking the limit of small jump lengths and small waiting times, yields the reaction-diffusion-advection equation... 

Here Npe > 1 means that transport in the chemical isolation layer is dominated by advection while Wpe < 1 implies that transport is dominated by diffusion. Advection and diffusion in either the cap isolation layer or bioturbation layer are not independent because advection tends to reduce diffusion gradients and diffusion tends to reduce the advective flux. In the cap isolation layer, a reasonable approximation is to assume that the flux is well-estimated by the dominant flux (either advection or diffusion). Solutions to the steady-state transport equations considering both diffusion and advection with and without reaction are feasible, but are algebraically more complicated and deviate significantly from solutions assuming only the dominant process in the relatively narrow range of approximately 0.3 < Npe < 3. Even within this range, the dominant process correctly estimates the flux within a factor of 2. 

The conservation of species is actually the law of conservation of mass applied to each species in a mixture of various species. The fluid element, as described in Sections 6.2.1.1 through 6.2.1.3, does not comprise of pure fluid with only one species, such as water, but of many species forming a multicomponent mixture. This law is mathematically described by the continuity equation for species, also known as the species equation, advection-diffusion equation, or convection-diffusion equation. If the species equation additionally includes a reaction term, it is known as the reaction-diffusion-advection equation. 

Above-mentioned reaction, diffusion and advection influence mass transfer in rock-water system. It is generally difficult to solve the differential equation including all these mechanisms. Thus, the two coupled models at constant temperature and pressure will be explained below. They are (1) reaction-fluid flow model, (2) reaction-diffusion model, (3) diffusion-fluid flow model. In addition to these coupled models, model taking into account the change in temperature will be considered. 

Yaimacopoulos, A.N., Tomlin, A.S., Brindley, J., Merkin, J.H., Pilling, MJ. Error propagation in approximations to reaction-diffusion-advection equations. Phys. Lett. A 223, 82-90 (1996b) Yarwood, G., Rao, S., Yocke, M., Whitten, G. Updates to the Carbon Bimd chemical mechanism CB05. Final Report to the US EPA, RT-0400675 (2005)... 

Consider the instantaneous release of a fixed mass of material, Q, into an infinite expanse of air (a ground surface will be added later). The coordinate system is fixed at the source. Assuming no reaction or molecular diffusion, the concentration C of material resulting from this release is given by the advection equation... 

The ratio vJD can then be used to calculate a chemical reaction rate for a nonconservative solute, S. To do this, the one-dimensional advection-diffusion model is modified to include a chemical reaction term, J. This new equation is called the one-dimensional advection-diffusion-reaction model and has the following form ... 

As we saw with the steady-state water-column application of the one-dimensional advection-diffusion-reaction equation (Eq. 4.14), the basic shapes of the vertical concentration profiles can be predicted from the relative rates of the chemical and physical processes. Figure 4.21 provided examples of profiles that exhibit curvatures whose shapes reflected differences in the direction and relative rates of these processes. Some generalized scenarios for sedimentary pore water profiles are presented in Figure 12.7 for the most commonly observed shapes. 

---