## Abstract

The surface of a turbulent liquid is visualized as consisting of a large number of chaotic eddies or liquid elements. Assuming that surface elements of a particular age have renewal frequencies that are integral multiples of a fundamental frequency quantum, and further assuming that the renewal frequency distribution is of the Boltzmann type, performing a population balance for these elements leads to the Danckwerts surface age distribution. The basic quantum is what has been traditionally called the rate of surface renewal. The Higbie surface age distribution follows if the renewal frequency distribution of such elements is assumed to be continuous. Four age distributions, which reflect different start-up conditions of the absorption process, are then used to analyse transient physical gas absorption into a large volume of liquid, assuming negligible gas-side mass-transfer resistance. The first two are different versions of the Danckwerts model, the third one is based on the uniform and Higbie distributions, while the fourth one is a mixed distribution. For the four cases, theoretical expressions are derived for the rates of gas absorption and dissolved-gas transfer to the bulk liquid. Under transient conditions, these two rates are not equal and have an inverse relationship. However, with the progress of absorption towards steady state, they approach one another. Assuming steady-state conditions, the conventional one-parameter Danckwerts age distribution is generalized to a two-parameter age distribution. Like the two-parameter logarithmic normal distribution, this distribution can also capture the bell-shaped nature of the distribution of the ages of surface elements observed experimentally in air–sea gas and heat exchange. Estimates of the liquid-side mass-transfer coefficient made using these two distributions for the absorption of hydrogen and oxygen in water are very close to one another and are comparable to experimental values reported in the literature.

## 1. Introduction

At the turn of the twentieth century, physics was faced with the ‘ultraviolet catastrophe’. The classical Rayleigh–Jeans formula was unable to explain the bell-shaped nature of the observed spectrum of blackbody radiation. The problem was resolved by Max Planck, who thereby laid the foundation of quantum mechanics. Planck invoked the now celebrated assumption that the energy of a particle (e.g. a gas molecule in a box or radiation in a cavity) could not take arbitrary values but only discrete multiples of some fundamental energy. With this energy quantum concept, Planck was able to derive his famous formula for the energy density of blackbody radiation, which reduces to the Rayleigh–Jeans formula in the limit of low frequencies [1]. In 1951, Danckwerts [2] published his classic paper on gas absorption in a turbulent liquid in which he presented the surface renewal model of mass transfer. This model visualizes the gas–liquid interface, where the absorption occurs, to be continuously replenished by fresh liquid elements that arrive from the bulk liquid. In contrast with the earlier model of Higbie [3] which had assumed that all surface elements had the same residence time at the gas–liquid interface (i.e. uniform age distribution), Danckwerts derived an exponential age distribution by using the postulate that all surface elements, irrespective of their individual ages, had equal probability of being replaced by fresh elements arriving from the bulk liquid. The only parameter that appears in this distribution is the frequency or rate of surface renewal *S*, which is the mean rate of production of fresh surface according to Danckwerts [2], and which depends on the prevailing hydrodynamic conditions. This well-known distribution, which carries his name and an early experimental verification of which was provided by Lamb *et al.* [4] in the case of gas absorption into a stirred liquid, has since found wide use in chemical engineering, and over the years, many applications and extensions of the surface renewal model have appeared in the literature [5–17] (table 1).

In mass-transfer studies, the surface renewal model is considered to be more realistic than the film model because it predicts that the liquid-side mass-transfer coefficient *k*_{L} is proportional to the square root of the diffusion coefficient *D* of the dissolved gas in the liquid—a result that has often been confirmed experimentally. For example, Kuthan & Broz [18] obtained experimental values of *k*_{L} for the absorption of helium, nitrogen and propane by a liquid film of aqueous ethylene glycol flowing over a smooth wetted wall and an expanded metal sheet. For the case of the wetted wall, they found *k*_{L} ∝ *D*^{0.5}, which is in agreement with the surface renewal model. For the expanded metal sheet, *k*_{L} ∝ *D*^{0.64} and the film-penetration model was found to be more applicable. According to Astarita [19], in the case of a liquid in contact with a solid or a more viscous liquid phase, *k*_{L} ∝ *D*^{2/3}. All of this is in sharp contrast to the film model which predicts that *k*_{L} ∝ *D*.

The surface of a turbulent liquid is characterized by bursting and chaotic movements of eddies that emerge from beneath the surface, and by the presence of turbulent sweeps, upwellings, downwellings and vortices that profoundly influence the interfacial mass-transfer process [20–22]. According to Komori *et al.* [20], mass transfer across the gas–liquid interface is dominated by large-scale, surface renewal eddies with the liquid-side mass-transfer coefficient being proportional to the square root of the surface renewal frequency. And according to Banerjee [21], liquid-side controlled gas exchange in the case of clean, non-breaking interfaces is well predicted by surface renewal models with the renewal frequency being that of the turbulent sweeps that impinge on the interface. Metzger & Dobbins [23] presented the following theoretical expression for *S*:
*C*_{1}–*C*_{3} are constants, *ρ* is the density of the fluid, *ν* is the kinematic viscosity of the fluid, *M*_{s} is the modulus of surface elasticity of the fluid and *E* is the specific rate of energy dissipation in the fluid as a whole (due to turbulent mixing). They deduced this relation by equating the resisting pressure (which opposes eddy motion) at the surface of the fluid to the product of the fluid density and the square of the eddy velocity, and by assuming that this velocity is proportional to the Kolomogoroff velocity factor (which is equal to the product of the kinematic viscosity of the fluid and the energy dissipation rate per unit mass of fluid raised to a power of 0.25). Equation (1.1) shows that there is a direct relation between *S* and *E*, and, keeping the random and erratic movements of the surface eddies in mind, the possibility arises that this motion is a result of the specific rate of energy dissipation in the fluid being quantized, which would be manifested in a quantization of the eddy renewal frequency. This raises the question of whether the Danckwerts age distribution has a deeper theoretical basis and can be derived by arguments similar to those used by Planck, with the surface renewal frequency *S* playing some kind of a fundamental role.

The answer to this question is provided in §2 of this manuscript. Specifically, the surface of a turbulent liquid is visualized as consisting of a large number of chaotic eddies or liquid elements. Assuming that surface elements of a particular age have renewal frequencies that are integral multiples of a fundamental frequency quantum, and further assuming that the renewal frequency distribution is of the Boltzmann type, performing a population balance for these elements leads to the Danckwerts surface age distribution. The basic quantum is what has been traditionally called the rate of surface renewal, *S*. Higbie surface age distribution follows if the renewal frequency distribution of such elements is assumed to be continuous. The development to be presented in what follows will lead to the derivation of four different unsteady-state age distributions. These four distributions correspond to four different hypotheses about the behaviour of liquid elements at the gas–liquid interface and reflect different initial states of the interface when the absorption commences. The first two age distributions (Cases 1 and 2) are two different versions of the Danckwerts model and the third one (Case 3) is based on the uniform and Higbie distributions, whereas the fourth (Case 4) is a mixed distribution. These distributions should be useful in modelling unsteady-state gas absorption in a gas–liquid reactor or a packed tower with the surface renewal model. Although these age distributions have appeared in previous and more recent works, as will be indicated later, none of them have traced their basis to the frequency quantum hypothesis that is proposed in this paper.

The frequency quantum hypothesis used herein is to be treated as an ansatz. The parameter *S*, which depends upon the level of turbulence or flow instability and fluid properties, is postulated to be a fundamental quantity which governs both the rate of gas absorption at the gas–liquid interface and the rate of dissolved-gas transfer to the bulk liquid from the interface. It is this basic parameter which determines the mean eddy renewal time—the relation between them is explored later.

A portion of the energy supplied for fluid mixing is transmitted to eddies of finer and finer scales (i.e. the Richardson cascade) and is ultimately dissipated by fluid friction. It may be speculated that another portion of this energy causes groups of fluid molecules to reach different excited states, and it is the transition of molecules between these states that leads to the release of discrete bursts of energy, which is reflected in the chaotic or discontinuous motion observed in turbulent flow. For example, Llaguno & Muriel [24] explained the multilevel turbulent data of the flow of nitrogen at room temperature obtained in a specially constructed experimental apparatus by postulating transitions between different rotational energy levels of the nitrogen molecule, which are caused by intermolecular collisions. By introducing a modification (postulated to arise from inelastic interactions of quantum origin among fluid molecules, which may be in different excited states) to the incompressible Navier–Stokes equation, Jirkovsky & Muriel [25] were able to derive the flattened paraboloid velocity profile in the turbulent flow of a fluid in a pipe and between two parallel plates, and also to fit the experimentally observed velocity profile with a single exponential-type equation in the near wall (viscous and buffer sublayers), intermediate and far away region in the turbulent flow of water in a channel. We note that the universal velocity profile for turbulent flow in a pipe is conventionally represented by three separate expressions depending upon the distance from the wall [26].

There is a vital aspect of the surface renewal model that has received scant attention in the literature where the focus has been on the rate of interphase transfer, i.e. the rate of gas absorption at the gas–liquid interface. Owing to surface renewal, there is a constant movement of liquid elements between the bulk liquid and the gas–liquid interface with a consequent net transfer of dissolved gas from the interface to the bulk liquid. In the case of physical absorption, as will be shown later, except at steady state, these two rates are not equal because of the accumulation of absorbed gas in the surface elements, while for absorption with chemical reaction, these two rates will not be equal even under steady-state conditions as demonstrated by Chatterjee [27] for the case of a first-order reaction. This is because of the consumption of dissolved gas by the chemical reaction in the surface elements. It should be noted that it is the rate of transfer of dissolved gas to the bulk liquid and not the rate of gas absorption at the gas–liquid interface that forms a part of the mass balance for dissolved gas in the bulk liquid, without which a rational design of gas–liquid contacting equipment is not possible. As pointed out by Chatterjee [27], this conceptual error has been made by some investigators and even by Danckwerts himself in Chapter 8 of his well-known book on gas–liquid reactions [28] in his general mass-balance equations for the dissolved-gas and liquid phase reactant in a packed column, although he presented an extensive discussion of this problem in the context of the film model in Chapter 6 of this book. To our knowledge, Loffler & Merchuk [29] and Merchuk [30], using the film penetration model of mass transfer, were the first to present the general form of the equation for the rate of dissolved-gas transfer for the case of absorption with a first-order reaction in a continuous flow, stirred tank reactor, assuming no gas phase resistance to mass transfer. The case of a finite gas-phase mass-transfer resistance was considered by Chatterjee [27] who, using the surface renewal and film penetration models, derived analytical expressions for the dissolved-gas transfer rate for the case of absorption with a first-order reaction under steady-state conditions. This development made possible an analysis with the surface renewal model of a gas–liquid reactor in which the dissolved gas is consumed by a first-order reaction [27], analysis of which has conventionally been made with the film model [28,31].

The use of the four age distributions derived in this work is illustrated in §3 of this manuscript where they are used in the theoretical modelling of transient physical absorption of a gas into a large volume of liquid (i.e. constant bulk liquid concentration of dissolved gas), assuming negligible gas-side mass-transfer resistance. In particular, explicit mathematical expressions for the transient rates of absorption of the gas at the gas–liquid interface and transfer of dissolved gas to the bulk liquid are derived. Section 4 of the manuscript presents numerical comparisons among the four age distribution cases.

Although experimental confirmation of the Danckwerts age distribution has been provided by Lamb *et al.* [4] for the case of a stirred liquid and by Lesage *et al.* [16] in the case of pipe flow, some studies by the physical oceanography community have reported that the age distribution of surface elements in air–sea heat and gas exchange, as measured both directly and indirectly, does not follow the exponential age distribution of Danckwerts but rather the logarithmic normal (LN) or Chi distributions [32–34]. This is addressed in §5, where it is shown that the generalized Danckwerts (GD) age distribution (derived in the electronic supplementary material, S1), which follows from an extension of the conventional Danckwerts age distribution, gives a shape which is very close to that of the LN distribution. Using both of these distributions, estimates of the liquid-side mass-transfer coefficient in the case of the absorption of hydrogen and oxygen in water at three different wind speeds, which were used by Garbe *et al.* [32] in their experiments, are presented. Section 6 of the manuscript offers some concluding remarks.

We note that the case of transient physical gas absorption in which the gas-phase resistance is finite and the bulk liquid concentration of the dissolved gas is a function of process time is a more difficult theoretical problem. For such a situation, equations for the rates of gas absorption and dissolved-gas transfer can be expressed in the form of convolution integrals; these have been provided by Chatterjee [35] using the age distribution for Case 1. However, only general conclusions were inferred and no numerical results were presented by him.

## 2. Theoretical development

### 2.1. Population balance

The surface (assumed to be of unit area) of a turbulent liquid is visualized as being populated by a mosaic of chaotic eddies or liquid elements. At any time *t*_{p} since the start of the process (e.g. absorption of a gas in a well-mixed liquid, crossflow membrane filtration, etc.), such elements are assumed to have different surface ages, which will range from zero to *t*_{p}. Let *F*(*t*, *t*_{p}) be the cumulative fraction of surface elements that have ages between 0 and *t* at process (or ‘clock’) time *t*_{p} and let *t*_{p} for elements with ages lying within the time interval from *t* to *t* + d*t* then gives
*average* renewal frequency of surface elements having ages lying in the time interval from *t* to *t* + d*t*.

### 2.2. Expression for R ¯ ( t )

*individual* renewal frequencies *R* of surface elements. This quantity, which is a reflection of the energy content of individual or groups of surface elements, is assumed to come in discrete packets (due to the random, discontinuous motions of the elements) which are integral multiples of a basic or *fundamental* frequency quantum. This basic quantum, as mentioned before, is taken to be equal to what has been conventionally called the rate of renewal of liquid elements at the surface, which is denoted by *S*. Thus, *R* = 0, *S*, 2*S*, 3*S*, …. If we now assume a Boltzmann-type distribution for the individual renewal frequency *R*, the average renewal frequency *t* at the interface is given by

In the theoretical development up to this point, as one conceptually moves from *S* (fundamental frequency quantum) through *R* (individual renewal frequency) to

### 2.3. General expression for *F*(*t*, *t*_{p})

Substituting equation (2.3) into equation (2.2) yields
*A* is a parameter whose value depends upon the initial state of the surface when the process was started at *t*_{p} = 0. Equation (2.5) contains the essence of the age distributions derived in this work.

Differentiating equation (2.5) with respect to *t* yields the density age distribution *f*(*t*,*t*_{p}), i.e.
*f*, e.g. that given by Danckwerts [2], equation (2.7) is the starting point and the cumulative age distribution *F* can then be obtained by integrating the density age distribution *f* with respect to *t*. However, in this method, the essential nature of the parameter *S*, that it is a fundamental frequency quantum, is obscured. This is in contrast with the procedure adopted in our work in which equation (2.1) for the cumulative age distribution *F* is the starting point whose time derivative is the density age distribution *f*.

In the development that follows, we will consider the four age distributions mentioned previously. The first two are unsteady-state versions of the Danckwerts (steady-state) model. Both have appeared previously in the literature and, as will be shown subsequently, reduce to the famous steady-state age distribution proposed by Danckwerts as the process time *t*_{p} → ∞.

#### 2.3.1. Case 1 (Danckwerts distribution)

The surface (assumed to be of unit area, as mentioned earlier) is instantaneously and completely formed at *t*_{p} = 0 with liquid elements flowing into it from the bulk liquid and departing from it to the bulk liquid at a constant rate for *t*_{p} ≥ 0. As
*A* in equation (2.5) is given by
*t* yields the following age distribution:

The above equation has been used in the modelling of crossflow ultrafiltration and microfiltration [36–40]. As *t*_{p} → ∞, it reduces to the well-known, steady-state age distribution *S*e^{−St}, which was proposed by Danckwerts [2].

#### 2.3.2. Case 2 (Danckwerts distribution)

In this case, it is assumed that there are liquid elements, visualized as being ‘blue’, that are already present on the surface at *t*_{p} = 0 when the process starts and ‘red’ elements begin displacing the blue elements by the mechanism of surface renewal. At any time *t*_{p}, the surface will therefore consist of a mixture of both red and blue elements, the population of the latter decreasing as *t*_{p} increases. The red elements will have ages lying within 0 ≤ *t* < *t*_{p} while the blue elements will have ages exactly equal to *t*_{p.} From this mental picture and equations (2.5) and (2.8), it follows that
*t* = *t*_{p} that becomes vanishingly small as *t*_{p} → ∞. This discontinuity is due to the fraction of blue elements present on the surface at process time *t*_{p}, which is equal to *u*(*t*) is the unit step function. Differentiating this equation with respect to *t* yields
*δ*(*t*) is the delta function.

Chung *et al.* [8] and Sada *et al.* [9] have previously presented equation (2.15) while a derivation of this equation, based on a stochastic population balance of interfacial fluid elements, has been provided by Fan *et al.* [14]. We note that, in the general equation for the unsteady-state age distribution provided by Chung *et al.* [8] of which equation (2.15) is a special case, the form of the steady-state age distribution has to be known or assumed *a priori* in order to derive the specific form of the unsteady-state age distribution, unlike in this work. In a recent paper, Zhang & Chatterjee [41] used this age distribution in the derivation of analytical expressions for the permeate flux decline and cake build-up in constant pressure, crossflow microfiltration. It is once again observed that, as *t*_{p} → ∞, equation (2.15) reduces to the steady-state Danckwerts age distribution.

#### 2.3.3. Case 3 (uniform and Higbie distributions)

We again assume a Boltzmann-type distribution for *R* but now assume that this distribution is continuous. The average renewal frequency *t* at the surface is obtained from
*t* yields the uniform age distribution:

The age distribution represented by equation (2.19) corresponds to the situation when the surface (initially empty) is ‘filling up’ with liquid elements that arrive from the bulk liquid, and is being formed; i.e. there is no outflow of elements from the surface back to the bulk liquid until a certain value of *t*_{p} = *τ* (Higbie time). Once the surface is completely formed at *t*_{p} = *τ*, liquid elements which flow into it from the bulk liquid also start departing from it to the bulk liquid at a constant rate. If it is assumed that thereafter (i.e. for *t*_{p} > *τ*) all liquid elements have the same residence time at the surface with elements flowing into and out of the surface in a *plug-flow* fashion, equation (2.19) reduces to the Higbie distribution:
*τ* is the residence time of liquid elements at the interface, which is a hydrodynamic parameter. Thus, the Higbie distribution (equation (2.20)) is a special case of the uniform distribution (equation (2.19)).

It can be observed that, as *St*_{p} → 0, i.e. as *S* → 0 (low surface renewal) or as *t*_{p} → 0 (near the start of the process), which implies that *St* → 0, equation (2.11) reduces to equation (2.19), i.e. the Danckwerts distribution (Case 1) reduces to the uniform distribution. This is analogous to the reduction of the Planck formula to the classical Rayleigh–Jeans formula in the limit of low frequencies in the energy spectrum radiated by a blackbody [1].

#### 2.3.4. Case 4 (mixed distribution)

This distribution can be thought of as being a combination of ideas contained in the distributions discussed earlier, and has also been used by Zhang & Chatterjee [41] to model permeate flux decline and cake build-up in constant pressure, crossflow microfiltration. For the benefit of the reader, we recapitulate their derivation. Initially, i.e. at *t*_{p} = 0, the surface is assumed to be empty of liquid elements. At *t*_{p} = 0+ when the process starts, the surface starts filling up with such elements, with no outflow of elements to the bulk liquid until the surface renewal mechanism is assumed to trigger at *t*_{p} = 1/*S* after which fresh liquid elements from the bulk begin to displace those already occupying the surface, which start flowing to the bulk liquid. The age distribution can be derived as follows.

Within the time interval 0 < *t*_{p} ≤ 1/*S*, the age distribution of the surface elements will be uniform and given by equation (2.19). Thus
*t*_{p} = 1/*S*, let all liquid elements that already occupy the surface be thought of as being blue, while those elements from the bulk that start displacing the blue elements from *t*_{p} = 1/*S* onwards be thought of as having the colour red. Similar to Case 2 of the Danckwerts distribution, at any time *t*_{p} > 1/*S* the surface will consist of both red and blue elements with ages lying from 0 to *t*_{p} − 1/*S* and *t*_{p} − 1/*S* to *t*_{p}, respectively. The age distribution for the red elements, which have ages in 0 ≤ *t* ≤ *t*_{p} − 1/*S*, will be given by *S*e^{−St} (see equation (2.15)), while that for the blue elements, which have ages in *t*_{p} − 1/*S* < *t* ≤ *t*_{p}, will be equal to

## 3. Unsteady-state gas absorption in a large volume of liquid

As an example of the application of the different age distributions, consider the physical, unsteady-state absorption of a gas into a large volume of liquid (i.e. constant bulk liquid concentration of dissolved gas). If the gas-side mass-transfer resistance is negligible, the concentration profile and instantaneous rate of absorption of the gas per unit area in a surface element (of infinite depth) having an age of *t* are given by [19,28]
*c*(*x*,*t*) is the dissolved-gas concentration in the element at location *x* (measured from the gas–liquid interface) and time *t*, *c*_{s} and *c*_{b} are the dissolved-gas surface and bulk liquid concentrations (assumed to be constants), respectively, and *D*, as mentioned earlier, is the diffusion coefficient of the dissolved gas in the liquid. The average rate of absorption of the gas at process time *t*_{p} [*R*_{abs}(*t*_{p})] can then be obtained from

We will now present an expression for the transient net rate of transfer *R*_{trans}(*t*_{p}) of the dissolved gas to the bulk liquid. The expression for *R*_{trans} for surface elements of infinite depth is given by

For later use, we define the following dimensionless quantities:
*z*, *y*) is the extended Euler gamma function.

### 3.1. Case 1 (Danckwerts distribution)

Using equation (2.11) and equations (3.1)–(3.11) yields the following:
*St*_{p} → 0 implies that *S* → 0 (low surface renewal) or *t*_{p} → 0 (near the start of absorption). Under either of these situations, equation (3.18) shows that the ratio of the rates of transfer to absorption varies linearly with the dimensionless process time

### 3.2. Case 2 (Danckwerts distribution)

From equation (2.15) and equations (3.1)–(3.11), it can be shown that

Equation (3.20) has been presented before by Chung *et al.* [8] who used it to derive an expression for the time-averaged mass-transfer coefficient. They found good agreement of this expression with experimental data on mass transfer in liquid films in wetted wall columns and obtained values of *S*, which increases with the liquid mass flow rate, varying from 0.294 to 2.42 s^{−1} as the gas–liquid contact time decreased from 3.4 to 0.41 s. As

### 3.3. Case 3 (uniform and Higbie distributions)

For this case, there are two time periods which have to be considered separately; these are 0 < *t*_{p} ≤ *τ* and *t*_{p} > *τ*. During the first period, the surface is being formed by the arrival of fresh liquid elements from the bulk liquid with no outflow of such elements back to the bulk liquid (i.e. *R*_{trans}(*t*_{p}) = 0). The surface is completely formed at *t*_{p} = *τ* after which, as indicated earlier, liquid elements flow to it from the bulk liquid and depart from it to the bulk liquid at a constant rate, with all elements residing on the surface for the same amount of time *τ*. Thus, *S* = 1/*τ*. For 0 < *t*_{p} ≤ *τ*, we use equations (2.19), (3.2), (3.3) and (3.5)–(3.9) to obtain

For *t*_{p} > *τ*, the age distribution, which is given by equation (2.20), is independent of *t*_{p} and the dimensionless absorption rate can be obtained by setting *t*_{p} > *τ*, all elements reside on the surface for the same amount of time *τ*. This rate is given by
*τ* = 1/*S*, equation (3.30), after using equations (3.1) and (3.8), is transformed into
*r* undergoes a step change from 0 to 1 at *t*_{p} = *τ*) and thereafter remains constant at this value.

### 3.4. Case 4 (mixed distribution)

From equation (2.23) and equations (3.1)–(3.11), it can be deduced that

We note from equations (3.35) and (3.38) that *r* undergoes a sudden transition from 0 to 2/3 at *r* asymptotically approaches a value of 1.

## 4. Numerical comparison among the four different age distribution cases

Table 2 presents a summary of the results for the four cases that correspond to the different age distributions.

Figures 1–4 show the cumulative age distributions (equations (2.10), (2.13), (2.18) and (2.24)) as functions of the dimensionless age *t**( = *St*) and dimensionless process time

It is observed from table 3 (*t** = 0.6, the values of *F* = 0.76, 0.45, 0.67 and 0.67 for Cases 1, 2, 3 and 4, respectively. Thus, the population of younger elements at the gas–liquid interface is greatest for Case 1, intermediate for Cases 3 and 4, and smallest for Case 2. As younger elements have a higher absorption rate (equation (3.2)) than older elements, the age-averaged absorption rate will be greatest for Case 1, intermediate for Cases 3 and 4, and lowest for Case 2. By contrast, younger elements will have a lower dissolved-gas concentration (equation (3.1)); therefore, the rate of dissolved-gas transfer to the bulk liquid will have an inverse relation to the rate of gas absorption at the gas–liquid interface. Thus, as illustrated in table 3, the average rate of transfer is higher for Case 2 compared to that for Case 1. The rate of transfer is zero for Cases 3 and 4 because the interface has not yet been completely filled up with liquid elements.

Table 4 shows the values of *F*, *t** = 0.6, the population of younger elements at the gas–liquid interface is greatest for Case 1, intermediate for Cases 2 and 4, and smallest for Case 3. At *t** = 0.9, the population of younger elements decreases in the order of cases 1, 4, 2 and 3, respectively. Not considering Case 3, which is an outlier because of the (imposed) sudden transition from a uniform distribution to the Higbie distribution at *t*_{p }= *τ* (i.e.

As *F* for Cases 1, 2 and 4 approaches the cumulative, steady-state age distribution:

Figures 5–8 show the age-averaged rates of gas absorption and dissolved-gas transfer to the bulk liquid as a function of process time in dimensionless coordinates. For the cases analysed in this work, the rate of absorption shows an exponential decay-like behaviour, decreasing from an initially high value to an asymptotic steady-state value as

Figure 9 presents the behaviour of the ratio *r* of the rate of transfer to that of absorption as a function of *r* = 0 at *r* near the origin is much greater for Case 2 when compared with Case 1 (see equations (3.24) and (3.18)). Also, *r* is in general considerably smaller for Case 1 compared with Case 2, with the difference becoming vanishingly small as *r* = 0 for *r* = 0 for

## 5. Comparison of the logarithmic normal and generalized Danckwerts age distributions

For simplicity, we will assume steady-state conditions (i.e. *m* and *σ* are the mean and standard deviation of the distribution. The mean eddy renewal or burst time *t*_{ren} can be obtained from

The GD age distribution, derived in the electronic supplementary material, S1, is an extension of the conventional Danckwerts age distribution. The extension follows from the introduction of a time lag into the basic differential equation that describes the conventional Danckwerts model (see equation (A.2) in electronic supplementary material). We note that the time delay concept has been used in the modelling of various phenomena such as biogas production by anaerobic digestion of organic wastes [43,44], population dynamics, forest fires and wavefronts in reaction–diffusion systems [45,46], and circadian rhythms [47]. The GD age distribution is given by *a* is a parameter of the distribution with the other parameter being *S*. We note that age distributions similar to equation (A.18) that involve the gamma function have been previously proposed by Harriott [6] and Bullin & Dukler [48]. However, these authors did not provide any theoretical justification, as done in the present work. We see that, for *a* = 0, equation (A.18) reduces to the standard or conventional Danckwerts age distribution, i.e. *f*(*t*) = *S*e^{−St}. Substitution of equation (A.18) into equation (5.2) yields
*a* = 0 (i.e. the standard Danckwerts age distribution), equation (5.4) gives *t*_{ren} = 1/*S*, i.e. the average eddy renewal time *t*_{ren} is equal to the reciprocal of the fundamental frequency quantum *S*, which is a well-known result. However, this is not true for other values of *a*. As the parameter *a* → ∞, *t*_{ren} → 0.5/*S*; thus for finite values of *a*, 0.5/S < *t*_{ren} < 1/*S* for the GD distribution.

Garbe *et al*. [32] experimentally measured the age distribution of surface renewal events on the Heidelberg Aeolotron in air–water heat exchange. They found a good fit of the LN distribution (equation (5.1)) to these data and reported values of *m*, *σ* and *t*_{ren} (calculated by equation (5.3)) at wind speeds of 2, 4.2 and 8 m s^{−1}; the values of these quantities are given in table 5. Figures 10–12 show plots of the LN age distribution using these values of *m* and *σ* for the three different wind speeds. The figures also show the GD age distribution (equation (A.18)), with the values of *a*, *S* and *t*_{ren} (calculated by equation (5.4)) given in table 6. It is evident that the LN and GD age distributions match one another quite closely. For the LN distribution, as the wind speed increases, the values of *m, σ* and *t*_{ren} decrease (table 5), whereas for the GD distribution, *a* and *S* increase, while *t*_{ren} decreases with wind speed (table 6). As the wind speed increases by a factor of 4, the mean eddy renewal or burst time decreases by a factor of 21 (LN distribution) and 16 (GD distribution), revealing the dramatic effect of turbulence. There is also a significant narrowing in the distribution of ages with a rise in wind speed as can be observed by comparing figures 10–12 among each other (see also in electronic supplementary material, figure A1 ).

The liquid side mass transfer coefficient *k*_{L} for the GD distribution is given by (see equations (A.16) and (A.17) in electronic supplementary material)
*a* = 0, the above equation reduces to

Hutchinson & Sherwood [49] investigated the absorption of eight different pure gases at 25°C in a stirred flask containing water whose surface was exposed to the gas. Table 7 shows values of the liquid-side mass-transfer coefficient *k*_{L}, which they obtained from dissolved-gas concentration measurements, for hydrogen and oxygen at two different stirrer speeds. This table also gives values of *D* (diffusion coefficient of dissolved gas in the liquid) that we estimated from values of the diffusion coefficient of the dissolved gas in water at 20°C reported in table III of [49] and by using the Wilke–Chang correlation [50] to correct them to 25°C. As the stirrer speed increases from 171 r.p.m. (revolutions per minute) to 1025 r.p.m., *k*_{L} varies from 1.06 × 10^{−5} to 3.25 × 10^{−5} m s^{−1} for hydrogen and 8.33 × 10^{−6} to 2.12 × 10^{−5} m s^{−1} for oxygen. Thus, a sixfold increase in the stirrer speed increases the value of *k*_{L} by a factor of 3.1 for hydrogen and 2.5 for oxygen.

Using the values of the diffusion coefficients of hydrogen and oxygen in water at 25°C (table 7), we calculated values of *k*_{L} for both the GD and LN distributions using equations (5.5) and (5.6), respectively, for the three different wind speeds used by Garbe *et al*. [32] in their experiments. The values of the parameters of the LN distribution (*m* and *σ*) and GD distribution (*a* and *S*) required for the calculations were obtained from tables 5 and 6, respectively. The results, shown in table 8 and figures 13 and 14, clearly demonstrate that both distributions predict very similar values of *k*_{L} at each wind speed, which, moreover, are comparable to those reported by Hutchinson & Sherwood [49] as can be observed from table 7. A fourfold increase in wind speed increases the value of *k*_{L} (for both gases) by factors of 3.4 and 3.1 for the LN and GD distributions, respectively. We also note that the increase in *k*_{L} with wind speed is not linear.

Although we have assumed steady-state conditions in this section, an unsteady-state version of the GD age distribution can also be derived whose details are not presented here. It is not obvious to the authors how to extend the LN age distribution to transient conditions.

## 6. Concluding remarks

This work is based on the primary assumption that the individual renewal frequency of eddies or elements at the surface of a turbulent liquid obeys a Boltzmann-type distribution. With a further assumption that the individual renewal frequency is either discrete or continuous, a population balance leads to the Danckwerts age distribution or a uniform age distribution, a special case of which is the Higbie age distribution. Thus, the Higbie and Danckwerts distributions are analogous to the Rayleigh–Jeans and Planck formulae, respectively, which describe the energy spectrum radiated by a blackbody.

The analysis of transient physical absorption of a gas into a large volume of liquid using the four age distributions revealed the crucial difference between the rate of absorption of the gas at the gas–liquid interface and the rate of its transfer to the bulk liquid. Under transient conditions, these two rates are not equal and have an inverse relationship. However, with the progress of absorption towards steady state, they approach one another.

The conventional, one-parameter Danckwerts age distribution can be generalized to a two-parameter age distribution which is able to capture the bell-shaped nature of the distribution of ages of surface elements that has been observed experimentally in air–sea heat and gas exchange. This distribution, called the GD age distribution, is equivalent to the two-parameter LN age distribution that has often been used by the physical oceanography community. Estimates of the liquid-side mass-transfer coefficient made using these two age distributions for the absorption of hydrogen and oxygen in water are very close to one another and are comparable to experimental values reported in the literature.

## Data accessibility

All theoretical calculations presented herein were performed with the equations developed in this paper. The only experimental data referred to or used in this work are those in the prior manuscripts of Chung *et al.* [8], Garbe *et al*. [32] and Hutchinson & Sherwood [49].

## Authors' contributions

C.M. was responsible for suggesting critical enhancements to the manuscript and its reorganization, helping with the literature search and checking the results. S.G.C. was responsible for developing the theoretical framework, performing the calculations and drafting the manuscript. Both authors gave their final approval for publication.

## Competing interests

We declare we have no competing interests.

## Funding

We received no funding for this study.

## Footnotes

Electronic supplementary material is available online at https://dx.doi.org/10.6084/m9.figshare.c.3810556.

- Received February 6, 2017.
- Accepted June 7, 2017.

- © 2017 The Authors.

Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.