The development of travelling waves in a simple isothermal chemical system III. Cubic and mixed autocatalysis

Autocatalytic chemical reactions can support isothermal travelling waves of constant speed and form. This paper extends previous studies to cubic autocatalysis and to mixed systems where quadratic and cubic autocatalyses occur concurrently. A + B → 2B, rate = k q ab , (1) A + 2B → 3B, rate = k c ab 2 . (2) For pure cubic autocatalysis the wave has, at large times, a constant asymptotic speed v 0 (where v 0 = 1/√2 in the appropriate dimensionless units). This result is confirmed by numerical investigation of the initial-value problem. Perturbations to this stable wave-speed decay at long times as t -3/2 e -1/8 t . The mixed system is governed by a non-dimensional parameter μ = k q / k c a 0 which measures the relative rates of transformation by quadratic and cubic modes. In the mixed case ( μ ≠ 0) the reaction-diffusion wave has a form appropriate to a purely cubic autocatalysis so long as μ lies between ½ and 0. When μ exceeds ½, the reaction wave loses its symmetrical form, and all its properties steadily approach those of quadratic autocatalysis. The value μ = ½ is the value at which rates of conversion by the two paths are equal.

The initiation of travelling reaction-diffusion waves in the chemical system governed by the quadratic autocatalytic or branching reaction A + B → 2B (rate k 1 ab) coupled with the decay or termination step B + B → C (rate k 4 b 2 ) is discussed. The system is described by the non-dimensional parameter K - k 4 / k 1 and parameters representing the local initial input of B. It is shown that a travelling wave of permanent form will develop for all K (and no matter how small the initial input of B). Bounds on the solution of the initial-value problem are obtained as well as numerical integrations of the governing equations. The structure of the permanent form travelling waves that arise is discussed in some detail, as well as the asymptotic limits K → 0 and K → ∞. The behaviour of the solution for this problem is compared with solutions found previously for other related simple autocatalytic systems with autocatalyst decay.


The possibility of travelling reaction-diffusion waves developing in the isothermal chemical system governed by the cubic autocatalytic reaction A + 2B → 3B (rate k 3 ab 2 ) coupled with either the linear decay step B → C (rate k 2 b ) or the quadratic decay step B + B → C (rate k 4 b 2 ) is examined. Two simple solutions are obtained,namely the well-stirred analogue of the spatially inhomogeneous problem and the solution for small input of the autocatalyst B. Both of these suggest that, for the quadratic decay case, a wave will develop only if the non-dimensional parameter k ═ k 4 / k 3 a 0 < 1 (where a 0 is the initial concentration of the reactant A), with there being no restriction on the initial input of the autocatalyst B. However, for the linear decay case the initiation of a travelling wave depends on the parameter v ═ k 2 / k 3 a 2 0 and that, in addition, there is an input threshold on B before the formation of a wave will occur. The equations governing the fully developed travelling waves are then considered and it is shown that for the quadratic decay case the situation is similar to previous work in quadratic autocatalysis with linear decay, with a necessary condition for the existence of a travelling-wave solution being that K < 1. However, the case of linear decay is quite different, with a necessary condition for the existence of a travelling wave solution now found to be v < 1/4 Numerical solutions of the equations governing this case reveal further that a solution exists only for v < v c , with v c ≈ 0.0465, and that there are two branches of solution for 0 < v < v c . The behaviour of these lower branch solutions as v → 0 is discussed. The initial-value problem is then considered. For the quadratic decay case it is shown that the uniform state a ═ a 0 , b ═ 0 is globally asymptotically stable (i. e. a → a 0 , b → 0 uniformly for large times) for all k > 1. For the linear decay case it is shown that the development of a travelling wave requires β 0 > v (where β 0 is a measure of the initial input of B) for v < v c . These theoretical results are then complemented by numerical solutions of the initial-value problem for both cases, which confirm the various predictions of the theory. The behaviour of the solution of the equations governing the travelling waves is then discussed in the limits K → 0, v → 0 and K → 1. In the first case the solution approaches the solution for K ═ 0 (or v =0) on the length scale of the reaction-diffusion front, with there being a long tail region of length scale O ( K -1 ) (or O ( v -1 )) in which the autocatalyst B decays to zero. In the latter case we find that the concentration of reactant A is 1 + O [(1 - k )] and autocatalyst B is O[(1 - k 2 ] with the thickness of the reaction-diffusion front becoming large, of thickness O [(1- k ) -3/2 ].


We examine the possibility of generating a propagating chemical wave front when a local input of an autocatalyst B is introduced into a uniform concentration of a reactant A. The autocatalysis is assumed to be of order m so that A -> B at a rate k 1 [A][B] m , while the autocatalyst decays to an inert product C at a rate of order n , B -> C, rate k 2 [B] n . The situation is examined for m, n > 1 and emphasis placed on obtaining threshold criteria for the development of reaction-diffusion fronts.


The possibility of travelling reaction–diffusion waves developing in the chemical system governed by the quadratic autocatalytic or branching reaction A + B → 2B (rate k 1 ab ) coupled with the decay or termination step B → C (rate k 2 b ) is examined. Two simple solutions are obtained first, namely the well-stirred analogue of the spatially inhomogeneous problem and the solution for small input of the reactant B. Both of these indicate that the criterion for the existence of a travelling wave is that k 2 < k 1 a 0 , where a 0 is the initial concentration of reactant A. The equations governing the fully developed travelling waves are then discussed and it is shown that these possess a solution only if this criterion is satisfied, i. e. only if k = k 2 / k 1 a 0 < 1. Further properties of these waves are also established and, in particular, it is shown that the concentration of A increases monotonically from its fully reacted state at the rear of the wave to its unreacted state at the front, while the concentration of B has a single hump form. Numerical solutions of the full initial value problem are then obtained and these do confirm that travelling waves are possible only if k < 1 and suggest that, when this condition holds, these waves travel with the uniform speed v 0 = 2√ (1 – k ). This last result is established by a large time analysis of the full initial value problem that reveals that ahead of the reaction–diffusion front is a very weak diffusion-controlled region into which an exponentially small amount of B must diffuse before the reaction can be initiated. Finally, the behaviour of the travelling waves in the two asymptotic limits k → 0 and k → 1 are treated. In the first case the solution approaches that for the previously discussed k = 0 case on the length scale associated with the reaction–diffusion front, with the difference being seen on a much longer, O ( k –1 ), length scale. In the latter case we find that the concentration of A is 1 + O (1 – k ) and that of B is O ((1 – k ) 2 ), with the thickness of the reaction–diffusion front being of O ((1 – k ) ½ ).


We study the isothermal autocatalytic system , A + n B → ( n + 1)B , where n = 1 or 2 for quadratic or cubic autocatalysis respectively. In addition, we allow the chemical species, A and B, to have unequal diffusion rates. The propagating reaction-diffusion waves that may develop from a local initial input of the autocatalyst, B, are considered in one spatial dimension. We find that travelling wave solutions exist for all propagation speeds v ≥ v * n ,where v * n is a function of the ratio of the diffusion rates of the species A and B and represents the minimum propagation speed. It is also shown that the concentration of the autocatalyst, B, decays exponentially ahead of the wavefront for quadratic autocatalysis. However, for cubic autocatalysis, although the concentration of the autocatalyst decays exponentially ahead of the wavefront for travelling waves which propagate at speed v = v * 2 , this rate of decay is only algebraic for faster travelling waves with v > v * 2 . This difference in decay rate has implications for the selection of the long time wave speed when such travelling waves are generated from an initial-value problem.


Author(s):  
Kate Fang Zhang ◽  
Xiao-Qiang Zhao

This paper is devoted to the investigation of the asymptotic behaviour for a reaction–diffusion model with a quiescent stage. We first establish the existence of the asymptotic speed of spread and show that it coincides with the minimal wave speed for monotone travelling waves. Then we obtain a threshold result on the global attractivity of either zero or positive steady state in the case where the spatial domain is bounded.


1985 ◽  
Vol 17 (11-12) ◽  
pp. 325-326 ◽  
Author(s):  
H J. G. W. Donker ◽  
P. Opic ◽  
H. P. de Vries

Ca. 60 % of the Dutch activated sludge plants consist of completely mixed systems, experiments have been carried out in completely mixed pilot plants to study the biological P-removal. The research was carried out in two pilot plants. The pilot plants consisted of: anaerobic reactor, anoxic reactor, aerobic reactor and a clarifier. All the reactors were completely mixed. Both plants were fed with settled domestic waste water at a sludge loading of 400 and 250 g COD/kg sludge.day respectively. The results are given below:sludge loading (g COD/kg sludge.day)400400250ratio Anaerobic : Anoxic : Aerobic1: 1:2,71:1:4,11:1:2,7P-removal (%)802875N-removal (%)505065COD-removal (%)858585 It has been shown that there is no significant difference between the results at the two different sludge loadings. Remarkable is the difference between the ratio 1:1:2,7 in combination with the internal recirculation flow anoxic-anaerobic of 160 % and the ratio 1:1:4,1 with a recirculation flow of 30 %. During the start-up at a sludge loading of 250 g COD/kg sludge.day and an internal recirculation flow of 30 %, bulking sludge developed almost immediately. The Premoval was completely disturbed. Increasing the internal recirculation flow to 160% had a positive effect on settling properties and P-removal. This investigation has pointed out that a completely mixed system is suitable for biological P-removal, without negatively affecting the nitrification. Important factors in the process are the ratio anaerobic:anoxic:aerobic and the recirculation flows.


1989 ◽  
Vol 4 (2) ◽  
pp. 141-167 ◽  
Author(s):  
J. H. Merkin ◽  
D. J. Needham

Sign in / Sign up

Export Citation Format

Share Document