Hostname: page-component-cd9895bd7-7cvxr Total loading time: 0 Render date: 2024-12-30T21:03:40.841Z Has data issue: false hasContentIssue false

Instability of a thin viscous film flowing under an inclined substrate: the emergence and stability of rivulets – CORRIGENDUM

Published online by Cambridge University Press:  17 September 2021

Abstract

Type
Corrigendum
Copyright
© The Author(s), 2021. Published by Cambridge University Press

Recent developments revealed an error in Lerisson et al. (Reference Lerisson, Ledda, Balestra and Gallaire2020), which propagates into Ledda et al. (Reference Ledda, Lerisson, Balestra and Gallaire2020). The mistake, reconducted from Lerisson et al. (Reference Lerisson, Ledda, Balestra and Gallaire2020), stems from the fact that the normalization of the curvature involves two independent length scales, which we erroneously assumed to be identical in the rest of this series of papers. The presence of the parameter $\frac {h_N}{l_c}\sqrt {\sin \theta }=\frac {h_N}{\ell _c^{*}}=1/\tilde {\ell }_c^{*}$, with $\ell_c=\sqrt {{\gamma }/{(\rho g)}}$, was overlooked and assumed to be $\tilde {\ell }_c^{*}=1$. The correct expression of the full curvature reads

\[ \kappa=\frac{\dfrac{\partial^{2} h}{\partial x^{2}}\left(1+\left(\dfrac{h_N}{\ell_c^{*}}\dfrac{\partial h}{\partial y}\right)^{2}\right)+\dfrac{\partial^{2} h}{\partial y^{2}}\left(1+\left(\dfrac{h_N}{\ell_c^{*}}\dfrac{\partial h}{\partial x}\right)^{2}\right)-2 \left(\dfrac{h_N}{\ell_c^{*}}\right)^{2} \dfrac{\partial h}{\partial x} \dfrac{\partial h}{\partial y} \dfrac{\partial^{2} h}{\partial x \partial y}}{\left(1+\left(\dfrac{h_N}{\ell_c^{*}}\dfrac{\partial h}{\partial x}\right)^{2}+\left(\dfrac{h_N}{\ell_c^{*}}\dfrac{\partial h}{\partial y}\right)^{2}\right)^{3 / 2}}. \]

Since the incorrect expression was employed in Ledda et al. (Reference Ledda, Lerisson, Balestra and Gallaire2020) (equations (3.3) and (A1)), some clarifications should be made about the validity of the results.

The numerical simulations of § 3.2 (and in particular figure3) are valid only for the case $\tilde {\ell }_c^{*}=1$. However, the discussion about the effect of the linear advection velocity $u=\cot \theta \tilde {\ell }_c^{*}$ remains valid.

In § 4, the performed linear stability analyses are valid again only for the case ${\tilde {\ell }_c^{*}=1}$ (Ledda et al. (Reference Ledda, Lerisson, Balestra and Gallaire2020) figures 4, 5, 6, 7). For each value of $\tilde {\ell _c^{*}}$, a different rivulet profile is identified (see corrigendum of Lerisson et al. Reference Lerisson, Ledda, Balestra and Gallaire2020), and thus the statement about the uniqueness of the rivulet in § 4.1 needs to be rectified. We complete the analysis by reporting in figure 1 of this corrigendum the results of the linear temporal and spatial stability analyses for different values of $\tilde {\ell }_c^{*}$. We recall that the perturbation with respect to the rivulet profile is assumed of the form $\eta = \hat {\eta }(y) \exp [\mathrm {i}(k_x x -\omega t)]$. The temporal stability analysis (i.e. $k_x \in \mathbb {R}$ is fixed and one looks for $\omega \in \mathbb {C}$) shows a reduction of the temporal growth rate $\mathrm {Im}(\omega )$ as $\tilde {\ell }_c^{*}$ increases. This reduction yields to the complete quenching of the secondary instability for sufficiently large $u$ and $\tilde {\ell }_c^{*}$ (figure 1). Similarly, the spatial stability analysis (i.e. $\omega \in \mathbb {R}$ is fixed and one looks for $k_x \in \mathbb {C}$) shows an analogous reduction of the spatial growth rate $-\mathrm {Im}(k_x)$ with $\tilde {\ell }_c^{*}$. In contrast, $\mathrm {Re}(\omega )$ in the temporal approach and $\mathrm {Re}(k_x)$ in the spatial one do not vary significantly with increasing $\tilde {\ell }_c^{*}$.

Figure 1. Results of the linear stability analysis for varying $\tilde {\ell }_c^{*}$. (a) Temporal stability analysis: temporal growth rate as a function of $k_x$. (b) Spatial growth rate from the spatial stability analysis (lines) and from the Gaster transformation (circles) as a function of $\omega$. (c) Real part of the complex frequency from the temporal approach as a function of $k_x$. (d) Variation of $\mathrm {Re}(k_x)$ with $\omega$, for the spatial approach. (e) Most unstable mode for $u=0.7$ and varying $\tilde {\ell }_c^{*}$, real (solid lines) and imaginary (dashed lines) parts, for the temporal approach. The different colours correspond to $\tilde {\ell }_c^{*}=1$ (blue), $\tilde {\ell }_c^{*}=1.1$ (black), $\tilde {\ell }_c^{*}=1.25$ (red), $\tilde {\ell }_c^{*}=1.43$ (maroon), $\tilde {\ell }_c^{*}=1.67$ (cyan), $\tilde {\ell }_c^{*}=2$ (green), $\tilde {\ell }_c^{*}=2.5$ (purple), $\tilde {\ell }_c^{*}=3.33$ (yellow), $\tilde {\ell }_c^{*}=5$ (orange), and $\tilde {\ell }_c^{*}=10$ (light blue).

An increase of $\tilde {\ell }_c^{*}$ leads to a decrease of the film thickness and thus of the hydro-static pressure gradients due to the gravity component normal to the substrate. As a consequence, the growth rates are reduced and the flow is thus stabilized. However, the quantities $\mathrm {Re}(\omega )$ and $\mathrm {Re}(k_x)$ are directly related to the advection of perturbations, which is dominated by $u$.

The results of § 5 are updated with the corrected theoretical spatial amplifications from the inlet. Following the procedure outlined in the paper, the theoretical amplification follows the same trend as the experimental one, without the need to fit the initial amplitude at $x=0$, that is now simply assumed to be the saturation value of the optical sensor, varying slightly from one set of experiments to the other (see figure 2 of this corrigendum). However, we have now included errorbars in figure 2 of this corrigendum (which replaces figure9 of the original paper), which reflect the relative tolerance in the measurement of $\varDelta$ and $h_N/\ell _c$ of $15\,\%$. This relates to the experimental error due to the undersampling in the identification of the position of the maximum thickness of the rivulet, since the resolution in the spanwise location was $\approx 1\,$mm, while the rivulet thickness varies from $h=1.5$ to $h=1.7$ in a region of $\approx1.5\ {\rm mm}$ extension. Besides, the relation (5.1)b of the original manuscript is now an approximation, since $1.65<\max (h)<1.71$ as $\tilde {\ell }_c^{*}$ varies, which gives a supplementary uncertainty of $\approx4$ % on $\max (h)$ and therefore on $h_N$.

Figure 2. Values of $\varDelta$ (blue dots) as a function of $h_{N} / l_{c}$, for different values of $\theta$. The black horizontal line denotes the plateau value due to the resolution of the optical sensor. The red lines denote the amplification estimated using the spatial stability analysis of § 4.4 and the size of the plate, i.e. $\varDelta =\varDelta _{0} \exp (-\mathrm {Im}(k_{x}) L)$, with initial amplitude the saturation value of the optical sensor. See text for the definition of the error bars.

Despite the assumptions reported in the paper and these uncertainties, the trends are very similar and thus the correlation between theoretical and experimental amplification is confirmed. We have also updated figure10 (figure 3 of this corrigendum) to include the marginal stability threshold below which perturbations are damped. The experimental values of $\varDelta$ and $h_N/\ell _c$ are mean values in the $15\,\%$ uncertainty region mentioned earlier. Agood agreement is observed between the theoretical onset of the instability and the growth of lenses on the rivulets. Note that the correlation between experiments and theory and the trends are more convincing in this corrected version than in the original manuscript. The main result of this corrigendum is the existence of a theoretical threshold below which rivulets are stable, in complement to the previous conclusions.

Figure 3. Results of the analysis in the $(\theta , h_N/l_c)$ plane: experimental measurements of $\varDelta$ (coloured dots) and inlet disturbance amplification $\varDelta /\varDelta _0 = \exp (-\mathrm {Im}(k_x)L)$ evaluated by the spatial stability analysis of § 4.4 (blue iso-contours). The red solid line denotes the iso-contour $\mathrm {Im}(\omega )=0$, which identifies the marginal stability threshold.

The results of § 6 remain valid since in the linear and weakly nonlinear models the curvature is linearized. While this does not modify the discussion, it should be specified in figure13 that the fully nonlinear simulation was performed with $\tilde {\ell }_c^{*}=1$.

We acknowledge the Swiss National Science Foundation under grant 200021_178971.

Declaration of interests

The authors report no conflict of interest.

References

REFERENCES

Ledda, P.G., Lerisson, G., Balestra, G. & Gallaire, F. 2020 Instability of a thin viscous film flowing under an inclined substrate: the emergence and stability of rivulets. J. Fluid Mech. 904, A23.CrossRefGoogle Scholar
Lerisson, G., Ledda, P.G., Balestra, G. & Gallaire, F. 2020 Instability of a thin viscous film flowing under an inclined substrate: steady patterns. J. Fluid Mech. 898, A6.CrossRefGoogle Scholar
Figure 0

Figure 1. Results of the linear stability analysis for varying $\tilde {\ell }_c^{*}$. (a) Temporal stability analysis: temporal growth rate as a function of $k_x$. (b) Spatial growth rate from the spatial stability analysis (lines) and from the Gaster transformation (circles) as a function of $\omega$. (c) Real part of the complex frequency from the temporal approach as a function of $k_x$. (d) Variation of $\mathrm {Re}(k_x)$ with $\omega$, for the spatial approach. (e) Most unstable mode for $u=0.7$ and varying $\tilde {\ell }_c^{*}$, real (solid lines) and imaginary (dashed lines) parts, for the temporal approach. The different colours correspond to $\tilde {\ell }_c^{*}=1$ (blue), $\tilde {\ell }_c^{*}=1.1$ (black), $\tilde {\ell }_c^{*}=1.25$ (red), $\tilde {\ell }_c^{*}=1.43$ (maroon), $\tilde {\ell }_c^{*}=1.67$ (cyan), $\tilde {\ell }_c^{*}=2$ (green), $\tilde {\ell }_c^{*}=2.5$ (purple), $\tilde {\ell }_c^{*}=3.33$ (yellow), $\tilde {\ell }_c^{*}=5$ (orange), and $\tilde {\ell }_c^{*}=10$ (light blue).

Figure 1

Figure 2. Values of $\varDelta$ (blue dots) as a function of $h_{N} / l_{c}$, for different values of $\theta$. The black horizontal line denotes the plateau value due to the resolution of the optical sensor. The red lines denote the amplification estimated using the spatial stability analysis of § 4.4 and the size of the plate, i.e. $\varDelta =\varDelta _{0} \exp (-\mathrm {Im}(k_{x}) L)$, with initial amplitude the saturation value of the optical sensor. See text for the definition of the error bars.

Figure 2

Figure 3. Results of the analysis in the $(\theta , h_N/l_c)$ plane: experimental measurements of $\varDelta$ (coloured dots) and inlet disturbance amplification $\varDelta /\varDelta _0 = \exp (-\mathrm {Im}(k_x)L)$ evaluated by the spatial stability analysis of § 4.4 (blue iso-contours). The red solid line denotes the iso-contour $\mathrm {Im}(\omega )=0$, which identifies the marginal stability threshold.