3.1. Lagrangian change of variables

We introduce the following change of independent variables which was first introduced by Plesset & Zwick (Reference Plesset and Zwick1952):

(3.1a,b)\begin{equation} x = \frac{Pe^{1/2}}{3} [r^{3} - \hat{R}(\tau)^{3}], \quad t = \tau, \end{equation}

in which the new dependent variables are $R(t) = \hat {R}(\tau )$, $c(x,t) = \hat {c}(r,\tau )$ and $m(t) = \hat {m}(\tau )$. The dependent variables without bars or hats correspond to dimensionless variables in the Lagrangian frame of reference. Equation (2.10) may be rewritten as

(3.2)\begin{equation} \frac{\partial {c}}{\partial {t}} = \frac{\partial}{\partial {x}} \left [ \left ( R^{3} + \frac{3 x}{Pe^{1/2}} \right)^{4/3} \frac{\partial {c}}{\partial {x}} \right ] \quad \mbox{for} \ x > 0. \end{equation}

The Rayleigh–Plesset equation (2.9) becomes

(3.3)\begin{equation} R \frac{\mbox{d}^{2} R}{\mbox{d} t^{2}} + \frac{3}{2} \left( \frac{\mbox{d} {R}}{\mbox{d} {t}} \right)^{2} = \frac{p_{g0} m^{\kappa}}{R^{3 \kappa}} - \frac{2}{We R} - 1 - \frac{4}{Re R} \frac{\mbox{d} {R}}{\mbox{d} {t}} - p_a \sin (2 {\rm \pi}\nu_a t) \end{equation}

and the conservation of the mass of gas in the bubble (2.13) becomes

(3.4)\begin{equation} \frac{\mbox{d} {m}}{\mbox{d} {t}} = \frac{4 {\rm \pi}M}{Pe^{1/2}} R^{4} \frac{\partial {c}}{\partial {x}} (0,t). \end{equation}

Equations (3.2)–(3.4) are singularly perturbed in $1/Pe^{1/2}$ having both an inner expansion when $x \ll Pe^{1/2}$ ($r=\hat {R}(\tau )+O(Pe^{-1/2})$) and an outer expansion when $x \sim Pe^{1/2}$ ($r=\hat {R}(\tau )+O(1)$). There are also multiple time scales. However, we will consider matched asymptotic expansions when $t \ll Pe^{1/2}$.

3.2. Inner region

We introduce asymptotic expansions taking the first two orders for the mass of gas in the bubble

(3.5a–c)\begin{equation} c \sim c_0,\quad R \sim R_0, \quad m \sim m_0 + \frac{m_1}{Pe^{1/2}}, \end{equation}

as $Pe \rightarrow \infty$. Therefore, the leading-order problem for the concentration of gas in the liquid on this length scale and time scale becomes

(3.6)\begin{equation} \frac{\partial {c_0}}{\partial {t}} = R_0^{4} \frac{\partial^{2} {c_0}}{\partial {x}^{2}} \quad \mbox{for} \ x > 0 \end{equation}

subject to the boundary condition at the bubble surface

(3.7)\begin{equation} c_0(0,t) = f(t) \equiv{-} 1 + {\rm d} p_{g0} \frac{m_0^{\kappa}}{R_0^{3 \kappa}}, \end{equation}

the initial condition

(3.8)\begin{equation} c_0(x,0)=0 \quad \mbox{for} \ x>0 \end{equation}

and a far-field condition which will follow from matching with the outer region.

At leading order, (3.4) becomes

(3.9)\begin{equation} \frac{\mbox{d} {m_0}}{\mbox{d} {t}} = 0, \end{equation}

with $m_0(0)=1$. We deduce that the mass of gas in the bubble does not change on this time scale at leading order, so $m_0 = 1$. At next order, we have

(3.10)\begin{equation} \frac{\mbox{d} {m_1}}{\mbox{d} {t}} = 4 {\rm \pi}M R_0^{4} \frac{\partial {c_0}}{\partial {x}} (0,t), \end{equation}

with the initial condition $m_1(0)=0$.

3.3. Outer region

We introduce the following change of variables ($y=O(1)$):

(3.11a,b)\begin{equation} y = \frac{x}{Pe^{1/2}} = \frac{1}{3} [r^{3} - \hat{R}(\tau)^{3}],\quad \tau = \tau. \end{equation}

The concentration of gas in the outer liquid region is denoted $C$. The equation for mass transport (2.10) for $C(y, t) = \hat {c}(r, \tau )$ becomes

(3.12)\begin{equation} \frac{\partial {C}}{\partial {t}} = \frac{1}{Pe}\frac{\partial}{\partial {y}} \left [ \left ( R^{3} + 3 y \right)^{4/3} \frac{\partial {C}}{\partial {y}} \right ] \quad \mbox{for} \ y > 0, \end{equation}

with far-field boundary condition

(3.13)\begin{equation} C \rightarrow 0 \quad \mbox{as} \ y \rightarrow \infty \end{equation}

and the initial condition

(3.14)\begin{equation} C = 0 \quad \mbox{on} \ t=0. \end{equation}

We introduce an expansion of the following form:

(3.15a,b)\begin{equation} C \sim C_0 + \frac{C_1}{Pe},\quad R \sim R_0, \end{equation}

as $Pe \rightarrow \infty$. At leading order in (3.12), we obtain

(3.16)\begin{equation} \frac{\partial {C_0}}{\partial {t}}=0, \end{equation}

with the initial condition $C_0(0)=0$ which is readily integrated to yield $C_0 = 0$.

3.4. Matching

We use limit matching technique in the form

(3.17)\begin{equation} \lim_{x \rightarrow \infty} c_0 = \lim_{y \rightarrow 0} C_0 \end{equation}

to obtain the matching condition on the inner region

(3.18)\begin{equation} c_0 \rightarrow 0 \quad\mbox{as} \ x \rightarrow \infty. \end{equation}

3.5. Leading-order inner solution

After the application of a Fourier sine transform to (3.6)–(3.8) and (3.18), the following initial value problem is obtained:

(3.19)\begin{equation} \frac{\partial {\tilde{c}}}{\partial {t}} + s^{2} R_0^{4} \tilde{c} = s R_0^{4} f \quad \mbox{with} \ \tilde{c}(s,0)=0, \end{equation}

where $f(t)$ was defined in (3.7) and

(3.20)\begin{equation} \tilde{c} (s,t) = \int_{x=0}^{\infty} c_0(x,t) \sin (s x) \,\mbox{d} x. \end{equation}

We solve this initial value problem to obtain

(3.21)\begin{equation} \tilde{c}(s,t) = \int_{\xi=0}^{t} f(\xi) R_0(\xi)^{4} s \exp \left({-}s^{2} \int_{\xi}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right) \mbox{d} \xi. \end{equation}

We then apply the inverse Fourier sine transform to find a solution for the concentration of gas in the liquid in the form

(3.22)\begin{equation} c_0(x,t) = \frac{2}{\rm \pi} \int_{s=0}^{\infty} \int_{\xi=0}^{t} f(\xi) R_0(\xi)^{4} \exp \left({-}s^{2} \int_{\xi}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right) \mbox{d} \xi s \sin( s x ) \,\mbox{d} s. \end{equation}

It is desirable to change the order of the integration in (3.22), but it is unclear whether the integral of the absolute value exists which is a requirement to apply Fubini's theorem. Therefore, integration by parts is applied to the inner integral which yields

(3.23)\begin{align} c_0(x,t) &= \frac{2}{\rm \pi} \int_{s=0}^{\infty} \left \{ f(t) - f(0) \exp \left({-}s^{2} \int_{0}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right) \right. \nonumber\\ &\quad \left. -\int_{\xi=0}^{t} f^{\prime}(\xi) \exp \left({-}s^{2} \int_{\xi}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right) \mbox{d} \xi \right \} \frac{\sin( s x )}{s}\, \mbox{d} s. \end{align}

Fubini's theorem may now be applied to change the order of the integration. Using tables of integral transforms (Erdélyi et al. Reference Erdélyi, Magnus, Oberhettinger and Tricomi1954), or otherwise, we determine the integrals in $s$ to yield

(3.24)\begin{align} c_0(x,t) &= f(t) - f(0) \,\mbox{erf} \left ( \frac{x}{2} \left \{ \int_{0}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right \}^{{-}1/2} \right ) \nonumber\\ &\quad -\int_{\xi=0}^{t} f^{\prime}(\xi)\, \mbox{erf} \left( \frac{x}{2} \left \{ \int_{\xi}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right \}^{{-}1/2} \right) \mbox{d} \xi, \end{align}

where

(3.25)\begin{equation} \mbox{erf} (\xi) = \frac{2}{\sqrt{\rm \pi}} \int_{y=0}^{\xi} \exp({-}y^{2}) \,\mbox{d} y. \end{equation}

The solution (3.24) may be seen to be an exact solution of the leading-order inner approximation by substituting this formula back into (3.6)–(3.8) and (3.18).

Fyrillas & Szeri (Reference Fyrillas and Szeri1994) adopted the idea of linear superposition when they split their model into an oscillatory problem and a smooth problem. Their solution was a linear sum of the solutions to the oscillatory and smooth problems. Linear sums of solutions are inappropriate for nonlinear problems, they are an engineering approximation at best. Fyrillas & Szeri (Reference Fyrillas and Szeri1994) then solved (3.6) with Fourier series, but this problem has a continuous spectrum due to the semi-infinite spatial domain, and Fourier series are more appropriate for problems with a discrete spectrum on a finite spatial domain. They concluded that the mass transport associated with the transient of the oscillatory problem is negligible, but this is clearly incorrect as we describe in the mass flux below. Furthermore, we note that their solution to the smooth problem makes sense only near the threshold when the bubble gas in the liquid averaged over a cycle of oscillation must neither increase nor decrease.

The concentration gradient at the bubble interface located at $x=0$ is determined by taking the partial derivative of this analytical solution with respect to $x$, as follows:

(3.26)\begin{align} \frac{\partial {c_0}}{\partial {x}}(0,t) ={-}\frac{f(0)}{\sqrt{\rm \pi}} \left \{ \int_{0}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right \}^{{-}1/2} - \frac{1}{\sqrt{\rm \pi}} \int_{\xi=0}^{t} f^{\prime}(\xi) \left \{ \int_{\xi}^{t} R_0(\eta)^{4} \,\mbox{d} \eta \right \}^{{-}1/2} \,\mbox{d} \xi. \end{align}

The first term on the right-hand side of (3.26) corresponds to a slowly fading memory of the initial condition. The second term is the history term as described by Peñas-López et al. (Reference Peñas-López, Parrales, Rodríguez-Rodríguez and van der Meer2016). It represents a memory of the entire preceding time history of the first time derivative of the concentration at the bubble surface. The distinction between this expression for the mass flux and the corresponding one in Peñas-López et al. (Reference Peñas-López, Parrales, Rodríguez-Rodríguez and van der Meer2016) is that the influence of convection associated with the bubble volume oscillations is accounted for here. Each contribution of the first derivative of the concentration of gas at the bubble surface is weighted by the integral of the fourth power of the bubble radius. Therefore, the flux of dissolved gas has been determined as an exact solution of the leading-order inner problem. The key theoretical difficulty in this topic has been overcome (Fyrillas & Szeri Reference Fyrillas and Szeri1994). The expression for the concentration gradient at the bubble interface (3.26) allows us to reduce the system of (2.9)–(2.15a–d) to a system of two ordinary differential equations for the growth and dissolution.

3.6. Simplified mathematical model

The simplified mathematical model is the dimensionless Rayleigh–Plesset equation

(3.27)\begin{equation} R \frac{\mbox{d}^{2} R}{\mbox{d} t^{2}} + \frac{3}{2} \left( \frac{\mbox{d} {R}}{\mbox{d} {t}} \right)^{2} = \frac{p_{g0} m^{\kappa}}{R^{3 \kappa}} - \frac{2}{We R} - 1 - \frac{4}{Re R} \frac{\mbox{d} {R}}{\mbox{d} {t}} - p_a \sin (2 {\rm \pi}\nu_a t) \end{equation}

and conservation of the mass of gas in the bubble

(3.28)\begin{align} \frac{\mbox{d} {m}}{\mbox{d} {t}} ={-}\frac{4 {\rm \pi}^{1/2} M}{Pe^{1/2}} R^{4} \left [ f(0) \left \{ \int_{0}^{t} R(\eta)^{4} \,\mbox{d} \eta \right \}^{{-}1/2} + \int_{\xi=0}^{t} f^{\prime}(\xi) \left \{ \int_{\xi}^{t} R(\eta)^{4} \,\mbox{d} \eta \right \}^{{-}1/2} \,\mbox{d} \xi \right], \end{align}

with the initial conditions

(3.29a–c)\begin{equation} R(0) = 1, \quad \frac{\mbox{d} {R}}{\mbox{d} {t}}(0) = 0, \quad m(0) = 1, \end{equation}

in which

(3.30a,b)\begin{equation} f(0) ={-} 1 + {\rm d} p_{g0},\quad f^{\prime}(\xi) ={-} 3 \kappa \,{\rm d} p_{g0} \frac{m^{\kappa}}{R^{3 \kappa + 1}} \frac{\mbox{d} {R}}{\mbox{d} {\xi}}. \end{equation}

We remain consistent by neglecting terms of the order of $1/ Pe^{1/2}$ in the expression for $f^{\prime }(\xi )$. Equation (3.28) has reverted to the form (3.4) as it is equivalent to (3.9) and (3.10). The fourth power of $R$ on the right-hand side of (3.28) encapsulates both the area and shell effects responsible for the growth of the bubble. In the remainder of this article, we study this simplified initial value problem (3.27)–(3.29a–c) using numerical methods.

4.1. Algorithm

The system of ordinary differential equations (3.27)–(3.29a–c) is integrated numerically with the Numerical Algorithms Group (NAG) routine D02EJF. The time integrals on the right-hand side of (3.28) are evaluated with the trapezoidal composite rule and Bode's composite rule. Bode's rule is exploited to integrate over longer time intervals, which reduces the computer memory requirements and speeds up the evaluation of history effects whilst retaining the accuracy of the trapezoidal rule on shorter time intervals.

4.2. Validation

Experimental results do not exist in the published literature on the time scale of a few hundred bubble oscillations on which the simplified model (3.27)–(3.29a–c) is valid. In order to compare with the experiments of Eller (Reference Eller1969), we note the following. The main change in the physical balances occurs when time is of the order of the Péclet number for mass transfer, the outer expansion in the above matched asymptotic expansions does not become active until this time. The introduction of multiple scales when time is of the order of the square root of the Péclet number for mass transfer is to maintain the same physics. The mass of gas in the bubble is being forced to retain a short memory of the concentration gradient on the order one time scale and not being allowed to develop a long memory. There is no new physical balance on this time scale of the order of the square root of the Péclet number for mass transfer. Therefore, the underlying physics on the time scale less than the Péclet number for mass transfer will be unchanged. If we are able to determine the underlying physics during the times when our simplified model is valid, then this same physics will be valid for all time less than the order of the Péclet number for mass transfer.

We consider the case of air-saturated water in a vertical pipe (Eller Reference Eller1969). A piston at the bottom of the pipe was used to excite a standing acoustic wave at $26.6$ kHz and a hypodermic needle was used to inject air bubbles into the pipe with a range of initial equilibrium radii. The following values for gas bubbles in water are adopted following the experimental conditions: $\varDelta = 10^{5}\ {\rm kg~m}^{-1}\ {\rm s}^{-2}$; $\sigma = 0.073\ {\rm Nm}^{-1}$; $\kappa = 1.4$, $\mu = 0.001\ {\rm Pas}$; $\rho = 10^{3}\ {\rm kg~m}^{-3}$; $\bar {\nu }_{a} = 26.6 \times 10^{3}$ Hz; $\bar {p}_a = 2 \times 10^{4}\ {\rm kg~m}^{-1}\ {\rm s}^{-2}$; $D = 2 \times 10^{-9}\ {\rm m}^{2}\ {\rm s}^{-1}$; and $k_{H} = 2.5 \times 10^{6}\ {\rm m}^{2}\ {\rm s}^{-2}$. The initial pressure of the bubble gas is given by $\bar {p}_{g0} = \varDelta + 2 \sigma / \bar {R}_{eqi}$ and the initial mass of gas in the bubble is given by

(4.1)\begin{equation} \bar{m}_i = \phi_{air} \rho_{air} \frac{4}{3} {\rm \pi}\bar{R}_{eqi}^{3}, \end{equation}

in which $\rho _{air}$ is the density and $\phi _{air}$ is the initial volume fraction of air in the bubble. As the solubility of oxygen in water is higher than the solubility of nitrogen, the air in water is assumed to be 36 % oxygen giving $\rho _{air}=1.3\ {\rm kg~m}^{-3}$. We adopt a reasonable estimate of the initial volume fraction of air in the bubble $\phi _{air}=0.9$. The saturation concentration of the gas in the liquid, $\bar {C}_{sat} = 4.15 \times 10^{-2}\ {\rm kg~m}^{-3}$, has been chosen to fit the threshold bubble radius of approximately $27\ \mathrm {\mu }\rm {m}$. We choose the initial uniform concentration of the gas in the liquid $\bar {C}_i = \bar {C}_{sat}$ (Eller Reference Eller1969). The Péclet number for mass transfer takes values in the range $8 \times 10^{4}$ to $3 \times 10^{5}$. In order to validate the simplified model (3.27)–(3.29a–c), numerical solutions were obtained for these parameter values.

Based on the bubble buoyancy, Eller (Reference Eller1969) measured the equilibrium bubble radius rather than the bubble radius itself. The dimensional equilibrium bubble radius, $\bar {R}_{eq}$, is evaluated as the solution to the following algebraic equation:

(4.2)\begin{equation} \bar{p}_{g0} \left ( \frac{\bar{m}}{\bar{m}_i} \right )^{\kappa} \left( \frac{\bar{R}_{eqi}}{\bar{R}_{eq}} \right)^{3 \kappa} - \frac{2 \sigma}{\bar{R}_{eq}} -(\bar{p}_{\infty} - \bar{p}_v)=0. \end{equation}

We firstly consider two typical situations where the equilibrium bubble radius is larger and smaller than the threshold bubble radius, respectively. An example of the bubble growth is shown in figure 1(a) for an initial equilibrium bubble radius of $\bar {R}_{eqi} = 30\ \mathrm {\mu }\rm {m}$. The equilibrium bubble radius oscillates around a rapidly increasing mean value. The oscillations initially exhibit beating but eventually the frequency of the oscillations settles down to the applied pressure frequency. An example of the bubble dissolution is shown in figure 1(b) for an initial equilibrium bubble radius of $\bar {R}_{eqi} = 20\ \mathrm {\mu }\rm {m}$. In this case, the equilibrium bubble radius oscillates around a rapidly decreasing mean value and the frequency of the oscillations almost immediately follows the applied pressure frequency. The mean values of the equilibrium bubble radius, $\bar {a}$, are well approximated by a square root dependence in time. The square root dependence indicates the relative significance of the diffusion of air in the water for the bubble growth and dissolution as it corresponds to the scaling in the Boltzmann variable ($x/t^{1/2}$). The approximation by a square root holds except for values of the initial equilibrium bubble radius of $\bar {R}_{eqi}$ between $55\ \mathrm {\mu }{\rm m}$ and $70\ \mathrm {\mu }\rm {m}$ which will be discussed below.

Figure 1. The equilibrium bubble radius, $\bar {R}_{eq}$, using a numerical solution of the simplified mathematical model (3.27)–(3.29a–c) and an initial equilibrium bubble radius of (a) $\bar {R}_{eqi}=30\ \mathrm {\mu }\rm {m}$ and (b) $\bar {R}_{eqi}=20\ \mathrm {\mu }\rm {m}$. Power law fits for the mean value of the equilibrium bubble radius of (a) $\bar {a}=\bar {R}_{eqi}+7 \times 10^{-7} (\bar {t} U \bar {R}_{eqi})^{1/2}$ and (b) $\bar {a}=\bar {R}_{eqi}-2.4 \times 10^{-6} (\bar {t} U \bar {R}_{eqi})^{1/2}$ are compared with the numerical solution.

Figure 2 shows the oscillations of the bubble radius corresponding to figure 1(a). The periods of oscillation in figure 2 agree with the period of the applied external pressure of approximately $3.8 \times 10^{-5}$ s. We note that the maxima are farther above the equilibrium bubble radius than the minima are below. This feature supports the area and shell mechanisms of rectified diffusion (Young Reference Young1989).

Figure 2. The bubble radius, $\bar {R}$, and the equilibrium bubble radius, $\bar {R}_{eq}$, using a numerical solution of the simplified mathematical model (3.27)–(3.29a–c) and an initial equilibrium bubble radius of $\bar {R}_{eqi}=30\ \mathrm {\mu }\rm {m}$ corresponding to figure 1(a).

In order to compare with the time-averaged growth rates in the experiments of Eller (Reference Eller1969), we require to average the growth rate of the mean value of the equilibrium bubble radius. As the mean values of the equilibrium bubble radius follow a square root dependence in time, the time-averaged growth rate varies with the reciprocal of the square root of the time period. The time period over which the growth rate is averaged has been chosen to be $0.05$ s to fit with the experiments. The change in the mean value of the equilibrium bubble radius over this time period is of the order of a tenth of a micron, which is small in comparison with the equilibrium bubble radius itself. Numerical simulations are undertaken for a range of initial equilibrium radii $\bar {R}_{eqi}$. Figure 3 compares the time-averaged rate of growth of the mean value of the equilibrium bubble radius with the experimental results, the agreement being excellent up to an initial equilibrium bubble radius of $\bar {R}_{eqi} = 55\ \mathrm {\mu }\rm {m}$.

Figure 3. Comparison of the time-averaged rate of growth of the mean value of the equilibrium bubble radius, $\mbox {d} \bar {a}/ \mbox {d} \bar {t}$, using a numerical solution of the simplified mathematical model (3.27)–(3.29a–c), the experimental results of Eller (Reference Eller1969) and the theoretical results of Fyrillas & Szeri (Reference Fyrillas and Szeri1994). In the experiment, a bubble in air-saturated water is subject to a sound field of 26.6 kHz.

The bubble radius has a resonance at approximately $\bar {R}_{eqi}=60\ \mathrm {\mu }{\rm m}$ in our results, which has also been observed in the theoretical results of Fyrillas & Szeri (Reference Fyrillas and Szeri1994). Bubble resonance corresponds to the appearance of a high frequency mode in the bubble oscillations. Due to the inaccuracy of the polytropic bubble models near resonance as discussed in Appendix C of Fyrillas & Szeri (Reference Fyrillas and Szeri1994) and the absence of the square root dependence in time, results are omitted for the initial equilibrium bubble radius $\bar {R}_{eqi}$ between $55\ \mathrm {\mu }\rm {m}$ and $70\ \mathrm {\mu }\rm {m}$. In comparison with the predictions of Fyrillas & Szeri (Reference Fyrillas and Szeri1994), which are also shown in figure 3, our resonance affects a larger range of the initial equilibrium bubble radius. Between $55\ \mathrm {\mu }\rm {m}$ and $70\ \mathrm {\mu }{\rm m}$, it is possible that microstreaming governs the bubble growth rate in the experimental results of Eller (Reference Eller1969). For values of the initial equilibrium bubble radius greater than $70\ \mathrm {\mu }\rm {m}$, the square root dependence in time returns and the bubble growth rate is again dominated by the diffusion of air in the water. The bubble growth rate for values above the resonance are also shown in figure 3. The present theory displays a much better agreement with the experimental data.

Despite the comparison in figure 3, we are not saying that Fyrillas & Szeri (Reference Fyrillas and Szeri1994) is without its applications. For example, stable single bubble sonoluminescence requires that the bubble is diffusively stable (Brenner, Hilgenfeldt & Lohse Reference Brenner, Hilgenfeldt and Lohse2002). The bubble gas averaged over a cycle of oscillation must neither increase nor decrease which corresponds to operating exactly on the threshold. Fyrillas & Szeri (Reference Fyrillas and Szeri1994) accurately predicts the threshold as shown in figure 3.

4.3. Parametric analysis

As discussed in the introduction, many authors have explained the underprediction of their theoretical growth rates in terms of the importance of surfactants. We now investigate how the time-averaged growth rates at $\bar {R}_{eqi}=40\ \mathrm {\mu }\rm {m}$ in figure 3 would vary in response to a reduction in surface tension. Figure 4(a) shows that the growth rate increases linearly in response to a reduction in surface tension. It is unclear how small quantities of surfactant could significantly increase the growth rate in terms of a small reduction in surface tension. Figure 5 of Lee et al. (Reference Lee, Kentish and Ashokkumar2005) has also reported an approximately linear decrease of the bubble growth rate with surface tension for the addition of the surfactant sodium dodecyl sulphate.

Figure 4. The time-averaged rate of growth of the mean value of the equilibrium bubble radius, $\mbox {d} \bar {a}/ \mbox {d} \bar {t}$, using a numerical solution of the simplified mathematical model (3.27)–(3.29a–c) with $\bar {R}_{eqi}=40\ \mathrm {\mu }\rm {m}$ for (a) a range of values for the surface tension and (b) a range of values of the initial uniform concentration of the gas in the liquid.

Figure 4(b) shows that the bubble growth rate increases linearly with the initial uniform concentration of the gas in the liquid. The bubble growth rate is significantly more sensitive to changes in the initial concentration of the gas than it is to the surface tension.