Physica A: Statistical Mechanics and its Applications
Reaction and diffusion in a lamellar structure: the effect of the lamellar arrangement upon yield
Introduction
Determining the yield of a chemical reaction in a fluid flow is a problem that has been of interest to many researchers [1], [2], [3], [4]. For many reactions, fluid mixing can have a significant effect on both the rate of reaction and the quality of the product [5] – the ratio between desired products and waste, formed in secondary, simultaneous (and unwanted) reactions. Clearly these issues are important for industrial chemical reaction processes, such as the Reaction Injection Moulding polymerisation process, where the reaction may fail to go to completion even if there is a stoichiometric balance of the reactants [6], [7].
In this paper we consider a model for chemical reactions between two initially segregated chemical species in a two-dimensional laminar chaotic fluid flow. This flow generates a complicated structure of thin striations of the two species [8]. The existence of these striations provides the basis for a one-dimensional ‘lamellar’ model of the fully two-dimensional problem, and so much work has been directed towards elucidating the structure of the striations [9], [10]. The interface between the chemical species remains sharp if they react with one another, and the reaction rate is taken to be infinite [11], [12]. The reaction is much harder to simulate numerically if instead the rate of reaction is finite, in which case the sharp interface between the species is replaced by a reaction zone in which the reactants are converted into products. Although there are then no sharply defined striations as such, it is useful, at least before the system becomes nearly homogeneous, to continue to use the terminology of striations and interfaces. At the early stage, it is sufficient to model the system as a series of reaction fronts at the interfaces [13], [14], [15], but as the species diffuse further across the striations, the reaction fronts meet, and interact: the full problem must then be considered. The time for this interaction to occur is governed by the rates of diffusion of the chemical species, and by the width of the striations.
Since the flow is incompressible, the striation widths are related to the length of the interface between the reactants: the thinner the striations, the longer the interface. It has been shown that the mean striation width is proportional to at large time, , where is the positive Lyapunov exponent of the chaotic flow [1]. The thinner the striation width, the faster the rate of reaction, and for the system we consider subsequently, the higher the yield. A detailed description of the effect of the striation width on the yield is given in [16].
Striation widths are an important ingredient in the one-dimensional lamellar model of a reaction system, in which the curved striations of the two-dimensional system are replaced by an array of exactly parallel lamellae. However, the striations produced by typical chaotic flows in both experimental and numerical studies have a wide range of widths [17], [18]. This range is a consequence of the tremendous variety of length stretches experienced at different points in the flow. For example, striation widths calculated in a recent numerical study varied over three orders of magnitude [19]. As noted in [20], the lamellar model does not provide an adequate approximation based on the mean striation width alone, and more complete information about the distribution of striation widths is needed.
Another factor that influences yield is the order in which striations of different widths are arranged [5], [20]. This aspect of the problem has received little attention hitherto. Previously, numerical simulations have assumed that the widths of adjacent striations are uncorrelated [18]. This assumption may not hold in realistic fluid flows; its influence is investigated here.
Section snippets
The reaction system
We consider the competitive–consecutive (or two-step) reactionwhere and are the rate constants for the two reactions, andAs an initial state for the system we assume that the reactants and are segregated in a one-dimensional array of lamellae (see Fig. 1) and that the products (desired product) and (unwanted by-product) are initially absent.
We let , , and denote the concentrations of the corresponding chemical species in which case the system is
Numerical simulations of the one-dimensional reaction–diffusion system
In this section we describe numerical simulations of the one-dimensional reaction–diffusion system – in the intervalwhich initially contains lamellae of species and species of two different widths, . The detailed distribution of and is described for each simulation by listing the lamellae in order, with a capital letter indicating a thick lamella and a small letter indicating a thin lamella. The arrangement in Fig. 1 is thus . The aim is to determine the effect upon
Discussion
We have considered a competitive–consecutive chemical reaction taking place in a one-dimensional lamellar array, and in particular we have shown that when the lamellae are of different widths, the order in which they lie can significantly affect the yield of the product. As the number of lamellae increases, the range of possible yields grows with the number of permutations of the lamellae. Approximating an arrangement of lamellae by a periodic array of lamellae of mean width tends to
Acknowledgements
This work was supported by the United Kingdom Engineering and Physical Sciences Research Council (EPSRC) under its Applied Nonlinear Mathematics Programme.
References (21)
- et al.
Chemical reactions in chaotic flows
Chem. Eng. J.
(1996) - et al.
The influence of segregation on the yield for a series–parallel reaction
Chem. Eng. Sci.
(1998) - et al.
Effect of stretching path on the course of polymerizations: applications to idealized unpremixed reactors
Chem. Eng. Sci.
(1987) - et al.
Lamellar modelling of reaction, diffusion and mixing in a two-dimensional flow
Chem. Eng. J.
(1998) - et al.
A two-stage reaction with initially separated reactants
Physica A
(1998) - et al.
Effect of striation thickness distribution on the course of an unpremixed polymerization
Chem. Eng. Sci.
(1987) - J.M. Ottino, The Kinematics of Mixing: Stretching, Chaos, and Transport, Cambridge University Press, Cambridge,...
- et al.
A fluid mechanical approach to turbulent mixing and chemical reaction. Part III: Computational and experimental results for the new micromixing model
Chem. Eng. Commun.
(1984) - et al.
Modelling of rapidly mixed fast-crosslinking exothermic polymerizations
A.I.Ch.E.J.
(1983) - et al.
Analysis of impingement mixing–reaction data: use of a lamellar model to generate fluid mixing information
Ind. Eng. Chem. Res.
(1989)
Cited by (24)
Quantifying magma mixing with the Shannon entropy: Application to simulations and experiments
2015, LithosCitation Excerpt :This advection/diffusion numerical system was proposed as a proxy to simulate magma mixing processes between degassed and crystal-free melts (Morgavi et al., 2013a; Perugini et al., 2003). The stretching and folding dynamics were triggered by applying the “sine-flow” chaotic mixing system, which is a well known prototypical chaotic flow (Clifford et al., 1998, 1999; Liu et al., 1994). As shown by Perugini et al. (2003), in order to make the system more similar to a natural scenario, the flow is continuously re-modulated by changing randomly the control parameter of the flow (i.e. parameter k in the works of Morgavi et al., 2013a; Perugini et al., 2004).
Shear-flow suppression of hotspot growth
2011, Physica D: Nonlinear PhenomenaDetermination of the local micromixing structure in laminar flows
2010, Chemical Engineering JournalChaotic mixing of a competitive-consecutive reaction
2004, Physica D: Nonlinear PhenomenaCitation Excerpt :Here, we adopt a simple form of the lamellar model, which follows the evolution of a single pair of lamellae. A more sophisticated variant of the model represents a more realistic set of lamellae, which have a range of different initial widths matching the striation width distribution measured from full two-dimensional simulations [14,19,20,21]. This latter approach clearly builds in more information from the full simulations, and might therefore be expected to provide a more accurate lamellar model, but is not considered here.
Kinematic significance of morphological structures generated by mixing of magmas: A case study from Salina Island (southern Italy)
2004, Earth and Planetary Science Letters