1 Introduction
In a recent paper Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) conducted experiments that showed that in a thin liquid film, on scales typical of those found in a kitchen sink, the circular jump produced by the normal impact of a round laminar jet onto an infinite plane is independent of the orientation of the surface. These experiments conclusively showed that gravity does not play a significant role in the origin/formation of these jumps – in sharp contrast with previous theoretical analyses. They also used conservation of energy, including both surface tension and gravity, to determine the radius of the jump. As we discuss below, on the scale of a kitchen sink the predicted radius is found to be almost independent of gravity, and is in excellent agreement with experiments, which cover fluids with a range of surface tension values.
Since the conclusions of Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) are markedly different from the accepted view of the role of gravity in these thin-film jumps the paper has attracted criticism, most particularly in a recent paper by Duchesne, Andersen & Bohr (Reference Duchesne, Andersen and Bohr2019) who, despite the experimental evidence, continue to question the role of surface tension in these flows. The underlying issue appears to be the current understanding of the interfacial flow which is that, in both in hydrostatics and in hydrodynamics, the influence of surface tension (when surface tension is uniform) is fully expressed as a pressure, commonly known as the Laplace pressure. The argument is that the force due to surface tension, which acts a tensile force at the line boundary of the interfacial surface, is mathematically equivalent to a pressure normal to the interface and, since the liquid velocity is tangential to the interface, this force can do no work. For hydrostatics, we have no argument with this viewpoint, but in hydrodynamics, we will show that the interfacial surface energy is not conserved, except when the liquid film is completely flat and surface tension force is trivially balanced.
Before proceeding further with this discussion we note that the importance of surface tension had previously been noted by Mohajer & Li (Reference Mohajer and Li2015) who measured the jump radius and the height of the liquid downstream of the jump for water and a water–surfactant solution (see their figures 4 and 13). They observed a significant increase in the jump radius when the surface tension was reduced by the addition of the surfactant. They also reported that for a range of flow rates the jump height remained constant and depended only on the surface tension of the liquid.
Despite this experimental evidence several other recent papers (Fernandez-Feria, Sanmiguel-Rojas & Benilov Reference Fernandez-Feria, Sanmiguel-Rojas and Benilov2019; Scheichl Reference Scheichl2019; Sen et al. Reference Sen, Chatterjee, Crockett, Ganguly, Yu and Megaridis2019; Wang & Khayat Reference Wang and Khayat2019) continue to support the previous gravity-based theory. For example, Fernandez-Feria et al. (Reference Fernandez-Feria, Sanmiguel-Rojas and Benilov2019) performed numerical simulations of the downward vertical impingement of a liquid jet onto a horizontal plate and compared their results with the experimental data from Button et al. (Reference Button, Davidson, Jameson and Sader2010), which were obtained for an upward directed jet onto a ceiling. Despite the opposite vertical orientations the authors concluded gravity is the dominant force.
The role of surface tension is subtle, and one objective of this paper is to clarify the contribution of surface tension to the dynamics. The other objective is to show that the jump radius can be predicted using momentum conservation and show that this is consistent with the energy-based approach used by Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018). In order to focus on the role of surface tension we predict the jump location using conservation of momentum, ignoring gravity completely.
Conventionally, the circular jump has been studied in an experimental arrangement where a vertical jet impinges onto a horizontal plate and either flows over a weir or off the edge of the plate. In this paper we are concerned with the situation in which the plate is effectively an infinite plane and we consider the jump before the spreading liquid film encounters either the weir or the edge. The observations of Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) show (for sufficiently high Reynolds numbers) that the jump, once formed, has a constant radius until the spreading film reaches the edge of the plate. In practice, once the film reaches the edge of the plate a different boundary condition is imposed that changes the depth of the downstream subcritical film, which, in turn, changes the position of the previously formed jump. Usually, in this later adjustment, since the subcritical region is significantly deeper than the supercritical region upstream of the jump, gravity can play a significant role. However, we will show that until the further downstream condition influences the jump, the initial jump location is determined by surface tension alone.
The paper is organised as follows. We begin in § 2 by a consideration of dimensional analysis which will set the parameter ranges that are appropriate for our theoretical analysis. The role of surface tension at the liquid–gas interface in a flowing, as distinct from stationary, liquid is derived and applied to a radially spreading thin film in § 3. A prediction for the jump radius is then obtained from application of conservation of momentum in both the radial and film-normal directions in § 4. The energy-based approach presented in Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) is revisited and compared with the momentum based approach in § 5. In § 6, we discuss the relative importance of surface tension and re-visit the analysis of Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) and discuss the flaw in their theory. Finally, our conclusions are given in § 7.
2 Dimensional analysis
We begin by considering the implications of dimensional analysis. The relative importance of gravity 
$g$ and surface tension depends on the fluid properties, the value of the surface tension 
$\unicode[STIX]{x1D6FE}$, the density 
$\unicode[STIX]{x1D70C}$ and kinematic viscosity 
$\unicode[STIX]{x1D708}$ and the flow characterised by the jet flow rate 
$Q$. Dimensional analysis shows that when either gravity or surface tension is ignored the jump radius scales, respectively, as
$$\begin{eqnarray}\displaystyle R_{ST}\sim \frac{Q^{3/4}\unicode[STIX]{x1D70C}^{1/4}}{\unicode[STIX]{x1D708}^{1/4}\unicode[STIX]{x1D6FE}^{1/4}}\quad \text{or}\quad R_{G}\sim \frac{Q^{5/8}}{\unicode[STIX]{x1D708}^{3/8}g^{1/8}}. & & \displaystyle\end{eqnarray}$$ In a particular flow the jump will either be caused by surface tension if 
$R_{ST}<R_{G}$ and gravity if 
$R_{ST}>R_{G}$. This criterion leads to a critical flow rate 
$Q_{C}$ given by
$$\begin{eqnarray}\displaystyle Q_{C}=\frac{\unicode[STIX]{x1D6FE}^{2}}{\unicode[STIX]{x1D708}\unicode[STIX]{x1D70C}^{2}g}, & & \displaystyle\end{eqnarray}$$ below which surface tension is the dominant force and above which gravity is important. For water 
$Q_{C}\,\approx \,500~\text{cm}^{3}~\text{s}^{-1}$ (i.e. ≈30 l min-1), which is significantly larger than the flow in the standard kitchen sink.
In appendix A (see (A 2)) we further show that including gravity provides a correction to the radius predicted by surface tension. This correction to the pure surface tension radius 
$R_{ST}$ takes the values 0.95 and 0.77, for 
$Q=2Q_{C}$ and 
$10Q_{C}$, respectively, which shows that even for 
$Q\sim 10Q_{C}$, the jump is dominated by the surface tension of the liquid. Table 1 gives values of 
$Q_{C}$ for liquids used in experiments (Bohr, Dimon & Putkaradze Reference Bohr, Dimon and Putkaradze1993; Bhagat et al. Reference Bhagat, Jha, Linden and Wilson2018; Duchesne et al. Reference Duchesne, Andersen and Bohr2019) and we see that 
$Q\ll Q_{C}$ in the experiments in water, and that the flow rates are only close to 
$Q_{c}$ in the experiments with ethylene glycol and silicone oil. However, as we note from (A 2), even in these latter cases this scaling analysis implies that surface tension is still dominant. Consequently, we conclude that all the experiments listed in table 1 are expected to fall within the range of parameters to which the theory developed below that ignores gravity is expected to apply.
It is also worth noting that apart from surface tension the only parameters that have been varied in experiments are the viscosity and jet flow rate. The scaling relations (2.1) have very similar dependences on these two parameters, which makes it very hard to distinguish between them on experimental grounds alone. This is possibly why some controversy persists about the interpretation of the experimental observations.
3 Theory
3.1 Surface forces and energy in hydrodynamics
Surface forces at a liquid interface arise from short range intermolecular interactions and the total short range force acting on a fluid element is determined by the surface area of the element, while the volume of the element is not relevant for such forces (Batchelor Reference Batchelor2000). The molecular origin of surface tension and surface energy is associated with intermolecular cohesive forces. The average free energy associated with a fluid molecule in the liquid bulk, where it is surrounded by similar molecules is independent of its position. However, within a distance less than the range of the cohesive forces 
$(10^{-9}~\text{m})$ from the interface, the liquid molecules have an additional free energy which is proportional to the interfacial surface area. A net exchange of liquid molecules between the interface and the bulk that necessarily accompanies a change in surface area will be associated with a net energy exchange.
Table 1. Parameters used in published experiments. The jet flow rate 
$Q$, the critical flow rate 
$Q_{C}$ at which gravity begins to play a role, the surface tension 
$\unicode[STIX]{x1D6FE}$, the kinematic viscosity 
$\unicode[STIX]{x1D708}$ and the density 
$\unicode[STIX]{x1D70C}$ of the fluid.
Now consider the case of an impacting liquid jet and follow the small material volume of fluid, 
$\unicode[STIX]{x03C0}a^{2}l$ in the jet (indicated as 0 in figure 1), where 
$a$ is the radius of the jet and 
$l$ is the length of the cylindrical fluid element. The surface area of a liquid interface changes as it spreads on the plane as a thin film of thickness 
$h$, and it is straightforward to show that the ratio of the final to initial surface area is 
$a/2h$. In the spreading film when 
$h$ decreases with radius and 
$a/2h\gg 1$, implying that there is an energy penalty as liquid molecules leave the bulk to increase the surface area of the spreading film. Whereas, at the hydraulic jump, an abrupt increase of liquid film thickness implies a release of the surface energy. This clearly demonstrates that even in a steady flow a spatial variation of liquid film thickness implies an active exchange of mechanical and surface energy. We now discuss the way in which this exchange of energy is accounted for.
Figure 1. An illustrative sketch of a circular hydraulic jump on an infinite plate that was reported experimentally by Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018), along with a history of the interfacial surface area of a material fluid volume, 
$V=\unicode[STIX]{x03C0}a^{2}l$ (indicated as 0), where 
$a$ is the radius of the jet. The annular surfaces are the top views of the surface of the material volume at different instances when it passes through locations 1 and 2 in the supercritical region and 3 in subcritical region. The surface area of the material volume is inversely proportional to the liquid film thickness 
$h$, which in the supercritical region where 
$h$ decreases implies an increase in surface area. At the location of hydraulic jump an abrupt change of liquid film thickness implies a release of the surface energy.
Figure 2. A schematic of the differential control volume showing the slope of the thin liquid film, and a contour, 
$C$ enclosing an interfacial surface, 
$S$.
3.2 Surface tension and the interfacial boundary condition
Fluid motion is governed by the Navier–Stokes equations which express conservation of momentum in a fluid continuum. For flows with an interface between two fluids, the Navier–Stokes equations do not express the surface tension force acting on the interface. This force is introduced as a normal stress boundary condition at the interface.
Consider a fixed surface 
$S$, with unit normal 
$\boldsymbol{n}$, bounded by a closed contour 
$C$ with the line element 
$\text{d}\boldsymbol{l}$ around the contour in the interface between two immiscible fluids, taken here for simplicity to be the common case of a liquid and a gas denoted by the subscripts 
$L$ and 
$G$, respectively, with constant surface tension 
$\unicode[STIX]{x1D6FE}$ (see figure 2). Since the surface tension force acts in a direction perpendicular to 
$\boldsymbol{n}$ and the contour 
$C$, continuity of the normal stress is expressed as
$$\begin{eqnarray}\displaystyle \int _{S}(\boldsymbol{T}_{G}-\boldsymbol{T}_{L})\,\boldsymbol{\cdot }\,\boldsymbol{n}\,\text{d}S+\unicode[STIX]{x1D6FE}\int _{C}\,\text{d}\boldsymbol{l}\times \boldsymbol{n}=0, & & \displaystyle\end{eqnarray}$$ where 
$\boldsymbol{T}=-p\boldsymbol{I}+\unicode[STIX]{x1D707}[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}]$ is the total stress, with pressure 
$p$ and velocity 
$\boldsymbol{u}$, and 
$\unicode[STIX]{x1D707}$ is the dynamic viscosity of the fluid. Using the vector identity and 
$\int _{C}\,\text{d}\boldsymbol{l}\times \boldsymbol{n}\equiv -\int _{S}(\unicode[STIX]{x1D735}_{s}\,\boldsymbol{\cdot }\,\boldsymbol{n})\boldsymbol{n}\,\text{d}S$,
$$\begin{eqnarray}\displaystyle \int _{S}(\boldsymbol{T}_{G}-\boldsymbol{T}_{L})\,\boldsymbol{\cdot }\,\boldsymbol{n}\,\text{d}S=\unicode[STIX]{x1D6FE}\int _{S}(\unicode[STIX]{x1D735}_{s}\,\boldsymbol{\cdot }\,\boldsymbol{n})\boldsymbol{n}\,\text{d}S. & & \displaystyle\end{eqnarray}$$ Since the surface 
$S$ is arbitrary, we then obtain the conventional form of the dynamic boundary condition
$$\begin{eqnarray}\displaystyle (\boldsymbol{T}_{G}-\boldsymbol{T}_{L})\,\boldsymbol{\cdot }\,\boldsymbol{n}=\unicode[STIX]{x1D6FE}(\unicode[STIX]{x1D735}_{s}\,\boldsymbol{\cdot }\,\boldsymbol{n})\boldsymbol{n}, & & \displaystyle\end{eqnarray}$$ where 
$\unicode[STIX]{x1D735}_{s}=[\boldsymbol{I}-\boldsymbol{n}\boldsymbol{n}]\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735}$ is the surface gradient, relating the jump in the normal stress to the curvature of the surface in agreement with Bush & Aristoff (Reference Bush and Aristoff2003) and valid pointwise in space. However, it is important to note that the surface tension forces act along the local tangents at different points of the closed contour 
$C$, and it is the net resultant of these tensile forces that acts along the normal to the interface. In the case of a stationary liquid and gas, the viscous term is zero, and this equation gives the usual Laplace pressure in the liquid associated with the curvature of the surface. In the case of a flowing liquid, (3.3) will give the correct net resultant force on an interfacial surface 
$S$, but does not give a detailed account of the tensile forces along the contour 
$C$. The velocity of the liquid may vary around the contour 
$C$ implying that liquid molecules crossing the boundary may not be conserved and need to be accounted for in the surface energy exchange. We consider this effect in the next section.
3.3 Interfacial surface energy conservation
Assuming the dynamic viscosity of the gas (for air and water the ratio 
$\unicode[STIX]{x1D707}_{a}/\unicode[STIX]{x1D707}_{w}\sim 10^{-2}$) is negligible compared to that of the liquid and that the velocities in the liquid and gas are comparable, and denoting pressure in the gas as 
$p_{G}$ and in the liquid as 
$p_{L}$, (3.1) can be written as
$$\begin{eqnarray}\int _{S}(-p_{G}+p_{L})\boldsymbol{n}\,\text{d}S-\unicode[STIX]{x1D707}\int _{S}\boldsymbol{n}\,\boldsymbol{\cdot }\,[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}]\,\text{d}S+\unicode[STIX]{x1D6FE}\int _{C}\,\text{d}\boldsymbol{l}\times \boldsymbol{n}=0,\end{eqnarray}$$ where 
$\boldsymbol{u}$ is the velocity in the liquid at the interface. The surface energy flux can be calculated by taking a dot product between the interfacial velocity at the contour 
$C$ and the forces to obtain (since 
$\boldsymbol{u}\,\boldsymbol{\cdot }\,\boldsymbol{n}=0$ on the interface)
$$\begin{eqnarray}\unicode[STIX]{x1D707}\int _{S}\boldsymbol{u}\,\boldsymbol{\cdot }\,(\boldsymbol{n}\,\boldsymbol{\cdot }\,[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}])\,\text{d}S=\unicode[STIX]{x1D6FE}\int _{C}\boldsymbol{u}\,\boldsymbol{\cdot }\,(\text{d}\boldsymbol{l}\times \boldsymbol{n}),\end{eqnarray}$$which can be re-expressed as
$$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D707}\int _{S}\{\boldsymbol{u}\,\boldsymbol{\cdot }\,[(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}]+\boldsymbol{n}\,\boldsymbol{\cdot }\,[(\boldsymbol{u}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}]\}\,\text{d}S & = & \displaystyle \unicode[STIX]{x1D6FE}\int _{S}\{(\unicode[STIX]{x1D735}_{I}\boldsymbol{\cdot }\boldsymbol{u})-[(\boldsymbol{n}\boldsymbol{\cdot }\unicode[STIX]{x1D735}_{I})\boldsymbol{u}]\boldsymbol{\cdot }\boldsymbol{n}\}\,\text{d}S,\nonumber\\ \displaystyle & = & \displaystyle \unicode[STIX]{x1D6FE}\frac{\text{D}S}{\text{D}t},\end{eqnarray}$$ the rate of change of the material surface area 
$S$ (see (3.1.5) of Batchelor (Reference Batchelor2000)), where the subscript 
$I$ of the operator 
$\unicode[STIX]{x1D735}_{I}$ is to emphasise that the operator is applied at the interface only.
Thus the increase in surface energy is supplied by the kinetic energy of the flow, which, of course, is zero in the case of stationary liquid. We also emphasise that the velocity field here is the surface velocity field, for which 
$(\unicode[STIX]{x1D735}_{I}\boldsymbol{\cdot }\boldsymbol{u})$ may or may not be equal to zero due to the addition (or loss) of fluid elements to the liquid interface as discussed above. The surface energy will only be conserved when the velocity field at the interface satisfies the condition 
$\{(\unicode[STIX]{x1D735}_{I}\boldsymbol{\cdot }\boldsymbol{u})-[(\boldsymbol{n}\boldsymbol{\cdot }\unicode[STIX]{x1D735}_{I})\boldsymbol{u}]\boldsymbol{\cdot }\boldsymbol{n}\}=0$ (which is satisfied for a flat surface for which net force due to surface tension is zero), otherwise there will be an active energy exchange with the bulk.
It is important to re-emphasise that the pointwise dynamic boundary condition in its conventional form (3.3), which if used to calculate the surface energy over a surface 
$S$ does not account for the variation in velocity along the contour 
$C$, will give a flawed result implying that the net surface energy exchange is zero. We emphasise that (3.6) is the correct interfacial boundary condition to be applied at the liquid interface and will be used in the analysis below.
4 Momentum conservation
We now apply conservation of momentum to the axisymmetric flow spreading radially from the point of impact of the jet on the plane, ignoring gravity. In cylindrical coordinates 
$(r,\unicode[STIX]{x1D703})$ with the origin at the point of jet impact, let 
$u,w$ be the radial and vertical velocity components, respectively.
First we determine the surface boundary condition. Following Bush & Aristoff (Reference Bush and Aristoff2003), we write the equation of the axisymmetric surface in implicit form
$$\begin{eqnarray}J(r,z)=z-h(r)=0,\end{eqnarray}$$which yields the vector normal to the film surface
$$\begin{eqnarray}\boldsymbol{n}=\frac{\unicode[STIX]{x1D735}J}{|\unicode[STIX]{x1D735}J|}=\frac{\hat{\boldsymbol{z}}-h^{\prime }\hat{\boldsymbol{r}}}{(1+h^{\prime 2})^{1/2}},\end{eqnarray}$$ where 
$\hat{\boldsymbol{r}}$ and 
$\hat{\boldsymbol{z}}$ are unit vectors in the radial and wall-normal directions, respectively, and 
$h^{\prime }=\text{d}h/\text{d}r$. We define the angle 
$\unicode[STIX]{x1D6FC}$ as the tangent to the surface defined by 
$h^{\prime }=\tan \unicode[STIX]{x1D6FC}$. Then 
$\cos \unicode[STIX]{x1D6FC}=1/(1+h^{\prime 2})^{1/2}$ and 
$\sin \unicode[STIX]{x1D6FC}=h^{\prime }/(1+h^{\prime 2})^{1/2}$, and (4.2) can also be written as
$$\begin{eqnarray}\boldsymbol{n}=\hat{\boldsymbol{z}}\cos \unicode[STIX]{x1D6FC}-\hat{\boldsymbol{r}}\sin \unicode[STIX]{x1D6FC}.\end{eqnarray}$$ We denote the radial velocity component at the surface 
$z=h$ as 
$u_{s}(r)$. From the kinematic boundary condition the vector surface velocity can be written as
$$\begin{eqnarray}\boldsymbol{u}=u_{s}(r)\hat{\boldsymbol{r}}+u_{s}(r)h^{\prime }(r)\hat{\boldsymbol{z}},\end{eqnarray}$$which implies that, at the interface,
$$\begin{eqnarray}\{(\unicode[STIX]{x1D735}_{I}\boldsymbol{\cdot }\boldsymbol{u})-[(\boldsymbol{n}\boldsymbol{\cdot }\unicode[STIX]{x1D735}_{I})\boldsymbol{u}]\boldsymbol{\cdot }\boldsymbol{n}\}=\frac{1}{r}\frac{\text{d}(u_{s}r)}{\text{d}r}+\frac{u_{s}h^{\prime }h^{\prime \prime }}{(1+h^{\prime ^{2}})}=\frac{1}{r(1+{h^{\prime }}^{2})^{(1/2)}}\frac{\text{d}\{u_{s}r(1+{h^{\prime }}^{2})^{(1/2)}\}}{\text{d}r}.\end{eqnarray}$$ Substituting (4.5) in (3.6) and recognising 
$\text{d}S=r(1+{h^{\prime }}^{2})^{(1/2)}\,\text{d}\unicode[STIX]{x1D703}\,\text{d}r$ yields
$$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D707}\int _{S}\{\boldsymbol{u}\,\boldsymbol{\cdot }\,[(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}]+\boldsymbol{n}\,\boldsymbol{\cdot }\,[(\boldsymbol{u}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}]\}\,\text{d}S=\unicode[STIX]{x1D6FE}\int _{S}\frac{\text{d}\{u_{s}r(1+{h^{\prime }}^{2})^{(1/2)}\}}{\text{d}r}\,\text{d}r\,d\unicode[STIX]{x1D703}. & & \displaystyle\end{eqnarray}$$ Note that these results depend only on the surface velocity and hold independent of the wall-normal velocity profile. However, in order to obtain an explicit momentum equation, we need to calculate the shape factors arising due to the velocity profile. We recognise that a similarity velocity profile in wall-normal direction will not satisfy the zero interfacial stress condition 
$(\boldsymbol{n}\boldsymbol{\cdot }\unicode[STIX]{x1D735})\boldsymbol{u}\neq 0$ (see appendix B). Nevertheless, the liquid film upstream of the hydraulic jump remains almost flat, and the shape factor is not expected to change much in the radial direction. Therefore, we use the Watson (Reference Watson1964) velocity in similarity form to calculate the shape factor
$$\begin{eqnarray}u(r,z)=u_{s}(r)f(\unicode[STIX]{x1D702}),\quad \unicode[STIX]{x1D702}\equiv \frac{z}{h(r)},\quad 0\leqslant \unicode[STIX]{x1D702}\leqslant 1,\end{eqnarray}$$ where 
$u_{s}$ as above is the surface velocity and 
$f(0)=0,f(1)=1$. Conservation of mass implies
$$\begin{eqnarray}\int _{0}^{h}ur\,\text{d}z=ru_{s}h\int _{0}^{1}f(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}\equiv C_{1}u_{s}rh=Q/2\unicode[STIX]{x03C0},\end{eqnarray}$$ and 
$C_{1}\equiv \int _{0}^{1}f(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}$ is an integration constant arising from the velocity profile. Incompressibility gives
$$\begin{eqnarray}w=-\int _{0}^{\unicode[STIX]{x1D702}}\frac{1}{r}\frac{\unicode[STIX]{x2202}ru}{\unicode[STIX]{x2202}r}\,\text{d}z=-\frac{1}{r}\frac{\text{d}}{\text{d}r}\left[(ru_{s}h)\int _{0}^{\unicode[STIX]{x1D702}}f(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}\right].\end{eqnarray}$$Then, using (4.8),
$$\begin{eqnarray}\displaystyle w=-u_{s}h\frac{\text{d}}{\text{d}r}\int _{0}^{\unicode[STIX]{x1D702}}f(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}=uh^{\prime }\unicode[STIX]{x1D702}=u_{s}h^{\prime }\unicode[STIX]{x1D702}f(\unicode[STIX]{x1D702}), & & \displaystyle\end{eqnarray}$$which automatically satisfies the kinematic boundary conditions at the wall and interface. We now use these expressions for the velocity in equations expressing momentum conservation in the radial and wall-normal directions.
4.1 Force on an axisymmetric thin film
For the control volume shown in figure 2, with the interface bounded by the closed contour 
$C$, the differential arc length 
$\text{d}\boldsymbol{l}=(r\,\text{d}\unicode[STIX]{x1D703})\hat{\unicode[STIX]{x1D73D}}|_{r}^{r+\text{d}r}-(\text{d}r)\hat{\boldsymbol{r}}|_{\unicode[STIX]{x1D703}}^{\unicode[STIX]{x1D703}+\text{d}\unicode[STIX]{x1D703}}-(\text{d}h)\hat{\boldsymbol{z}}|_{\unicode[STIX]{x1D703}}^{\unicode[STIX]{x1D703}+\text{d}\unicode[STIX]{x1D703}}$. Consequently, the force due to surface tension on this interface is
$$\begin{eqnarray}\text{d}F_{\unicode[STIX]{x1D6FE}}=\unicode[STIX]{x1D6FE}(\text{d}\boldsymbol{l}\times \boldsymbol{n})=\unicode[STIX]{x1D6FE}(\text{d}\unicode[STIX]{x1D703}r\cos \unicode[STIX]{x1D6FC}\hat{\boldsymbol{r}})|_{r}^{r+\text{d}r}+\unicode[STIX]{x1D6FE}(\text{d}\unicode[STIX]{x1D703}r\sin \unicode[STIX]{x1D6FC}\hat{\boldsymbol{z}})|_{r}^{r+\text{d}r}+\left.\unicode[STIX]{x1D6FE}\left(\frac{\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\hat{\unicode[STIX]{x1D73D}}\right)\right|_{\unicode[STIX]{x1D703}}^{\unicode[STIX]{x1D703}+\text{d}\unicode[STIX]{x1D703}}.\end{eqnarray}$$ Considering the circular symmetry, in the limit of 
$\text{d}r$ and 
$\text{d}\unicode[STIX]{x1D703}\rightarrow 0$
$$\begin{eqnarray}\text{d}F_{\unicode[STIX]{x1D6FE}}=\left(\unicode[STIX]{x1D6FE}\text{d}(r\cos \unicode[STIX]{x1D6FC})-\unicode[STIX]{x1D6FE}\frac{\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\right)\text{d}\unicode[STIX]{x1D703}\hat{\boldsymbol{r}}+\unicode[STIX]{x1D6FE}\,\text{d}(r\sin \unicode[STIX]{x1D6FC})\,\text{d}\unicode[STIX]{x1D703}\hat{\boldsymbol{z}}.\end{eqnarray}$$Therefore, from (3.5) and (4.11), the radial and vertical components of the normal stress at the free surface can be written as
$$\begin{eqnarray}\displaystyle F_{\unicode[STIX]{x1D6FE},r} & \equiv & \displaystyle \int _{r}^{r+dr}\unicode[STIX]{x1D707}(\boldsymbol{n}\,\boldsymbol{\cdot }\,[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}])\,\boldsymbol{\cdot }\,\hat{\boldsymbol{r}}\left.\left(r\,\text{d}\unicode[STIX]{x1D703}\frac{\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\right)\right|_{h}=\int _{r}^{r+\text{d}r}P\sin \unicode[STIX]{x1D6FC}r\,\text{d}\unicode[STIX]{x1D703}\frac{\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\nonumber\\ \displaystyle & & \displaystyle +\,\unicode[STIX]{x1D6FE}r\,\text{d}\unicode[STIX]{x1D703}\cos \unicode[STIX]{x1D6FC}|_{r}^{r+\text{d}r}-\unicode[STIX]{x1D6FE}\frac{\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\,\text{d}\unicode[STIX]{x1D703},\end{eqnarray}$$and
$$\begin{eqnarray}\displaystyle F_{\unicode[STIX]{x1D6FE},z} & \equiv & \displaystyle \int _{r}^{r+\text{d}r}\unicode[STIX]{x1D707}(\boldsymbol{n}\,\boldsymbol{\cdot }\,[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}])\,\boldsymbol{\cdot }\,\hat{\boldsymbol{z}}\left.\left(r\,\text{d}\unicode[STIX]{x1D703}\frac{\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\right)\right|_{h}\nonumber\\ \displaystyle & = & \displaystyle \left.\int _{r}^{r+dr}Pr\,\text{d}\unicode[STIX]{x1D703}\,\text{d}r+r\,\text{d}\unicode[STIX]{x1D703}\unicode[STIX]{x1D6FE}\sin \unicode[STIX]{x1D6FC}\right|_{r}^{r+\text{d}r},\end{eqnarray}$$ where 
$P=p_{L}-p_{G}$. As we will see below, application of momentum conservation requires expressions for the radial gradients of these forces, and as axisymmetry implies we can drop 
$\text{d}\unicode[STIX]{x1D703}$, these are given by
$$\begin{eqnarray}\displaystyle \frac{\text{d}F_{\unicode[STIX]{x1D6FE},r}}{\text{d}r} & = & \displaystyle \frac{\text{d}}{\text{d}r}\left(\int _{r}^{r+\text{d}r}\unicode[STIX]{x1D707}(\boldsymbol{n}\,\boldsymbol{\cdot }\,[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}])\,\boldsymbol{\cdot }\,\hat{\boldsymbol{r}}\left.\left(\frac{r\,\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\right)\right|_{h}\right)\nonumber\\ \displaystyle & = & \displaystyle rP\tan \unicode[STIX]{x1D6FC}-\frac{\unicode[STIX]{x1D6FE}}{\cos \unicode[STIX]{x1D6FC}}+\unicode[STIX]{x1D6FE}\frac{\text{d}(r\cos \unicode[STIX]{x1D6FC})}{\text{d}r}\end{eqnarray}$$and
$$\begin{eqnarray}\frac{\text{d}F_{\unicode[STIX]{x1D6FE},z}}{\text{d}r}=\frac{\text{d}}{\text{d}r}\left(\int _{r}^{r+\text{d}r}\unicode[STIX]{x1D707}(\boldsymbol{n}\,\boldsymbol{\cdot }\,[\unicode[STIX]{x1D735}\boldsymbol{u}+(\unicode[STIX]{x1D735}\boldsymbol{u})^{\text{T}}])\,\boldsymbol{\cdot }\,\hat{\boldsymbol{z}}\left.\left(\frac{r\,\text{d}r}{\cos \unicode[STIX]{x1D6FC}}\right)\right|_{h}\right)=rP+\unicode[STIX]{x1D6FE}\frac{\text{d}(r\sin \unicode[STIX]{x1D6FC})}{\text{d}r}.\end{eqnarray}$$4.2 Momentum balance in the radial direction
In the absence of gravity, the momentum equation in the radial direction is
$$\begin{eqnarray}\unicode[STIX]{x1D70C}\left(u\frac{\unicode[STIX]{x2202}u}{\unicode[STIX]{x2202}r}+w\frac{\unicode[STIX]{x2202}u}{\unicode[STIX]{x2202}z}\right)=\frac{1}{r}\frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rr})}{\unicode[STIX]{x2202}r}+\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D70F}_{rz}}{\unicode[STIX]{x2202}z},\end{eqnarray}$$ where 
$\unicode[STIX]{x1D70F}=\unicode[STIX]{x1D61B}+p\boldsymbol{I}$ is the total stress consisting of the deviator stress 
$\unicode[STIX]{x1D61B}$ and the pressure 
$p$ (in cylindrical coordinates). Since the pressure in the gas at the surface is constant (atmospheric) and the film is thin, the radial pressure gradient in the liquid is only a function of radius caused by the curvature of the surface.
We use incompressibility (
$\unicode[STIX]{x1D735}\,\boldsymbol{\cdot }\,\boldsymbol{u}=0$ in the bulk fluid), integrate (4.16) across the film from the wall to the interface, use axisymmetry to eliminate 
$\text{d}\unicode[STIX]{x1D703}$, substitute for the velocity from (4.7) and (4.9) and apply the surface tension boundary condition (4.5), to obtain
$$\begin{eqnarray}\frac{\text{d}}{\text{d}r}\int _{0}^{h}\unicode[STIX]{x1D70C}ru^{2}\,\text{d}z=-\int _{0}^{h}\frac{\text{d}p}{\text{d}r}r\,\text{d}z+rp\tan \unicode[STIX]{x1D6FC}-\frac{\unicode[STIX]{x1D6FE}}{\cos \unicode[STIX]{x1D6FC}}+\unicode[STIX]{x1D6FE}\frac{\text{d}(r\cos \unicode[STIX]{x1D6FC})}{\text{d}r}-\unicode[STIX]{x1D707}r\left.\left(\frac{\unicode[STIX]{x2202}u}{\unicode[STIX]{x2202}z}+\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}r}\right)\right|_{0}.\end{eqnarray}$$ Noting that conservation of mass (4.8) implies 
$u_{s}rh=\text{const.}$, we obtain
$$\begin{eqnarray}C_{2}\left[\unicode[STIX]{x1D70C}u_{s}rh\frac{\text{d}u_{s}}{\text{d}r}\right]=-\frac{\text{d}p}{\text{d}r}rh+rp\tan \unicode[STIX]{x1D6FC}-\frac{\unicode[STIX]{x1D6FE}}{\cos \unicode[STIX]{x1D6FC}}+\unicode[STIX]{x1D6FE}\frac{\text{d}(r\cos \unicode[STIX]{x1D6FC})}{\text{d}r}-\unicode[STIX]{x1D707}r\left.\left(\frac{\unicode[STIX]{x2202}u}{\unicode[STIX]{x2202}z}+\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}r}\right)\right|_{0},\end{eqnarray}$$ where 
$C_{2}=\int _{0}^{1}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}$ is a second integration constant. Equation (4.18) can be written in terms of the interface slope
$$\begin{eqnarray}\displaystyle C_{2}\left[\unicode[STIX]{x1D70C}u_{s}rh\frac{\text{d}u_{s}}{\text{d}r}\right] & = & \displaystyle -\frac{\text{d}p}{\text{d}r}rh+rph^{\prime }-\unicode[STIX]{x1D6FE}(1+h^{\prime 2})^{1/2}+\unicode[STIX]{x1D6FE}\frac{1}{(1+h^{\prime 2})^{1/2}}\nonumber\\ \displaystyle & & \displaystyle -\,\unicode[STIX]{x1D6FE}\frac{rh^{\prime }h^{\prime \prime }}{(1+h^{\prime 2})^{3/2}}-\unicode[STIX]{x1D70F}_{w}r,\end{eqnarray}$$ where 
$\unicode[STIX]{x1D70F}_{w}=-\unicode[STIX]{x1D707}(\unicode[STIX]{x2202}u/\unicode[STIX]{x2202}z+\unicode[STIX]{x2202}w/\unicode[STIX]{x2202}r)|_{0}$ is the wall shear stress.
4.3 Momentum balance in the wall-normal direction
We now apply conservation of momentum in the wall-normal 
$z$ direction in the differential control volume shown in figure 2. Since the film is thin, the pressure is independent of 
$z$ and for an axisymmetric flow
$$\begin{eqnarray}\unicode[STIX]{x1D70C}ur\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}r}+\unicode[STIX]{x1D70C}wr\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}=\frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}+r\frac{\unicode[STIX]{x2202}\unicode[STIX]{x1D70F}_{zz}}{\unicode[STIX]{x2202}z}.\end{eqnarray}$$As before, we integrate across the film to obtain
$$\begin{eqnarray}\int _{0}^{h}\unicode[STIX]{x1D70C}ur\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}r}\,\text{d}z+\int _{0}^{h}\unicode[STIX]{x1D70C}wr\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\,\text{d}z=\int _{0}^{h}\frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h}.\end{eqnarray}$$Substituting for the velocity from (4.7) and (4.9) gives
$$\begin{eqnarray}\displaystyle & & \displaystyle \unicode[STIX]{x1D70C}u_{s}rh\frac{\text{d}(u_{s}h^{\prime })}{\text{d}r}\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}-\unicode[STIX]{x1D70C}h^{\prime 2}u_{s}^{2}r\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}-\unicode[STIX]{x1D70C}\left.\frac{u_{s}^{2}rh^{\prime 2}\unicode[STIX]{x1D702}^{2}f^{2}(\unicode[STIX]{x1D702})}{2}\right|_{0}^{1}\nonumber\\ \displaystyle & & \displaystyle \quad +\,\unicode[STIX]{x1D70C}h^{\prime 2}u_{s}^{2}r\left.\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}+r\frac{\unicode[STIX]{x1D70C}w^{2}}{2}\right|_{0}^{h}=\frac{\text{d}F_{\unicode[STIX]{x1D6FE},z}}{\text{d}r}-\left.\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}|_{h}.\qquad\end{eqnarray}$$Substituting (4.15) into (4.22) yields
$$\begin{eqnarray}\displaystyle & & \displaystyle \unicode[STIX]{x1D70C}u_{s}rh\frac{\text{d}(u_{s}h^{\prime })}{\text{d}r}\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}\nonumber\\ \displaystyle & & \displaystyle \quad =\left.rP+\unicode[STIX]{x1D6FE}\frac{\text{d}(r\sin \unicode[STIX]{x1D6FC})}{\text{d}r}-\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h},\nonumber\\ \displaystyle & & \displaystyle \quad =rP+\unicode[STIX]{x1D6FE}\sin \unicode[STIX]{x1D6FC}+\unicode[STIX]{x1D6FE}r\cos \unicode[STIX]{x1D6FC}\frac{\text{d}\unicode[STIX]{x1D6FC}}{\text{d}r}-\left.\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h},\nonumber\\ \displaystyle & & \displaystyle \quad =rP+\frac{\unicode[STIX]{x1D6FE}h^{\prime }}{(1+h^{\prime 2})^{1/2}}+\frac{\unicode[STIX]{x1D6FE}rh^{\prime \prime }}{(1+{h^{\prime }}^{3/2})}-\left.\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h}.\end{eqnarray}$$ In the thin liquid film upstream of the hydraulic jump, the interface slope remains small and we will ignore the higher-order terms in 
$\text{d}h/\text{d}r$. Applying this approximation and re-arranging (4.23) gives an expression for the curvature of the film
$$\begin{eqnarray}\displaystyle & & \displaystyle h^{\prime \prime }(\unicode[STIX]{x1D70C}u_{s}^{2}rh\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}-\unicode[STIX]{x1D6FE}r)+h^{\prime }\left(\unicode[STIX]{x1D70C}u_{s}rh\frac{\text{d}u_{s}}{\text{d}r}\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}-\unicode[STIX]{x1D6FE}\right)\nonumber\\ \displaystyle & & \displaystyle \quad =\left.rP-\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h}.\end{eqnarray}$$ Finally, substituting 
$\unicode[STIX]{x1D70C}u_{s}rh(\text{d}u_{s}/\text{d}r)$ from the radial momentum balance (4.19) gives
$$\begin{eqnarray}\displaystyle & & \displaystyle [C_{3}\unicode[STIX]{x1D70C}u_{s}^{2}rh-\unicode[STIX]{x1D6FE}r]h^{\prime \prime }+\frac{C_{3}}{C_{2}}h^{\prime }\left(-\frac{\text{d}P}{\text{d}r}rh+rPh^{\prime }-\unicode[STIX]{x1D6FE}(1+h^{\prime 2})^{1/2}+\unicode[STIX]{x1D6FE}\frac{1}{(1+h^{\prime 2})^{1/2}}\right.\nonumber\\ \displaystyle & & \displaystyle \quad -\left.\unicode[STIX]{x1D6FE}\frac{h^{\prime }h^{\prime \prime }}{(1+h^{\prime 2})^{3/2}}-\unicode[STIX]{x1D70F}_{w}r-\frac{C_{2}}{C_{3}}\unicode[STIX]{x1D6FE}\right)=rP-\left.\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h},\end{eqnarray}$$ where 
$C_{3}=\int _{0}^{1}\unicode[STIX]{x1D702}f^{2}(\unicode[STIX]{x1D702})\,\text{d}\unicode[STIX]{x1D702}$ is a third integration constant.
Ignoring higher-order terms in 
$h^{\prime }$ and re-arranging gives an expression for the curvature of the film
$$\begin{eqnarray}h^{\prime \prime }=\frac{\displaystyle -h^{\prime }\frac{C_{3}}{C_{2}}\left[-\frac{\text{d}P}{\text{d}r}rh+r\unicode[STIX]{x1D70F}_{w}-\frac{C_{2}}{C_{3}}\unicode[STIX]{x1D6FE}\right]+rP+\left.\int \frac{\unicode[STIX]{x2202}(r\unicode[STIX]{x1D70F}_{rz})}{\unicode[STIX]{x2202}r}\,\text{d}z\right|_{0}-2\unicode[STIX]{x1D707}r\left.\frac{\unicode[STIX]{x2202}w}{\unicode[STIX]{x2202}z}\right|_{h}}{C_{3}\unicode[STIX]{x1D70C}u_{s}^{2}rh(1-We^{-1})},\end{eqnarray}$$ where the Weber number 
$We$ is defined by
$$\begin{eqnarray}We\equiv \frac{C_{3}\unicode[STIX]{x1D70C}u_{s}^{2}h}{\unicode[STIX]{x1D6FE}}.\end{eqnarray}$$ Consequently, we predict a singularity in the curvature of the film at a critical radius where the film thickness is such that 
$We=1$. This criterion gives the location of the hydraulic jump.
4.4 Revisiting the radial momentum balance
We now revisit the radial momentum balance (4.19), which can also be written as
$$\begin{eqnarray}\displaystyle C_{2}[\unicode[STIX]{x1D70C}u_{s}rh\,\text{d}u_{s}] & = & \displaystyle -rh\,\text{d}p+rp\,\text{d}h-\unicode[STIX]{x1D6FE}(1+h^{\prime 2})^{1/2}\,\text{d}r+\unicode[STIX]{x1D6FE}\frac{\text{d}r}{(1+h^{\prime 2})^{1/2}}\nonumber\\ \displaystyle & & \displaystyle -\,\unicode[STIX]{x1D6FE}\frac{rh^{\prime }h^{\prime \prime }\,\text{d}r}{(1+h^{\prime 2})^{3/2}}-\unicode[STIX]{x1D70F}_{w}r\,\text{d}r,\end{eqnarray}$$ and apply it at the jump radius. We note that the jump is a singularity (see figure 3b), where 
$\text{d}h/\text{d}r=\tan \unicode[STIX]{x1D6FC}\rightarrow \infty$, 
$\text{d}r\rightarrow 0$, 
$\text{d}h=H$, a finite quantity, and 
$\unicode[STIX]{x1D6FC}$ changes from 
$0$ to 
$\unicode[STIX]{x03C0}/2$, and the radial velocity changes from 
$u_{s}$ to 
$0$. Substituting the trigonometric forms of the functions and integrating (4.28) at the jump location 
$r=R$ gives
$$\begin{eqnarray}\displaystyle \int C_{2}\unicode[STIX]{x1D70C}u_{s}Rh\,\text{d}u_{s} & = & \displaystyle -\int Rh\,\text{d}p+\int Rp\,\text{d}h-\left.\int \unicode[STIX]{x1D6FE}\sin \unicode[STIX]{x1D6FC}\,\text{d}h+\unicode[STIX]{x1D6FE}R\cos \unicode[STIX]{x1D6FC}\right|_{0}^{\unicode[STIX]{x03C0}/2}\nonumber\\ \displaystyle & & \displaystyle +\,\int \unicode[STIX]{x1D6FE}\cos \unicode[STIX]{x1D6FC}\,\text{d}r-\int \unicode[STIX]{x1D70F}_{w}R\,\text{d}r.\end{eqnarray}$$Then in scaled terms (4.29) can be written
$$\begin{eqnarray}-C_{2}\unicode[STIX]{x1D70C}{u_{s}}^{2}Rh\approx RH\left(p-\frac{\unicode[STIX]{x1D6FE}}{R}\right)-\unicode[STIX]{x1D6FE}R.\end{eqnarray}$$ Since, at the jump, the pressure 
$p$ scales as 
$\unicode[STIX]{x1D6FE}/R$ (see Bush & Aristoff Reference Bush and Aristoff2003), the first term on the right-hand side of (4.29) is zero, which gives 
$We=1$ as the condition for the hydraulic jump. Thus conservation of radial momentum gives a similar result for the jump condition.
Figure 3. (a) Schematic velocity profiles in a flow with conventional flat surface assumption with zero radial viscous stress (black) and liquid interface with zero tangential stress but non-zero radial stress (red). The surface tension force retards the flow in the radial direction and accelerates it in the wall-normal direction, giving non-zero radial and wall-normal viscous stresses at the surface (3.1). (b) An idealised hydraulic jump.
5 Relation to energy conservation
For a better physical understanding we can also get the interfacial energy from a consideration of the interfacial force. From (4.11), the energy flux associated with surface tension force on the control volume (figure 2) is
$$\begin{eqnarray}\text{d}F_{\unicode[STIX]{x1D6FE}}\boldsymbol{\cdot }(u\hat{\boldsymbol{r}}+w\hat{\boldsymbol{z}})=\unicode[STIX]{x1D6FE}(\text{d}\unicode[STIX]{x1D703}ur\cos \unicode[STIX]{x1D6FC})|_{r}^{r+\text{d}r}+\unicode[STIX]{x1D6FE}(\text{d}\unicode[STIX]{x1D703}wr\sin \unicode[STIX]{x1D6FC})|_{r}^{r+\text{d}r}.\end{eqnarray}$$Circular symmetry implies that there is no net flux of fluid in the azimuthal direction and so substituting (4.4) into (5.1) yields
$$\begin{eqnarray}\text{d}F_{\unicode[STIX]{x1D6FE}}\boldsymbol{\cdot }(u\hat{\boldsymbol{r}}+w\hat{\boldsymbol{z}})=\unicode[STIX]{x1D6FE}(\text{d}\unicode[STIX]{x1D703}u_{s}r(1+h^{\prime 2})^{1/2})|_{r}^{r+\text{d}r},\end{eqnarray}$$which is the equivalent result to (4.6) derived by considering the rate of change of interfacial surface area and work done by the surface tension force. Consequently, in a control volume approach, the force due to surface tension, which appears through the normal stress boundary condition, can be incorporated as a surface force on the circumference of the control volume consistent with the analysis in Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018).
6 Relative importance of surface tension
We now return to the question of the relative importance of gravity and surface tension in these flows. Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) argued that the analysis of Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) is wrong and that surface tension does not play a significant role in the formation of these thin-film jumps. Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) wrote the energy equation over a volume 
$V$ bounded by a closed surface 
$A$ (their equation (13)) as
$$\begin{eqnarray}-\int _{A}\left[v_{j}\left({\textstyle \frac{1}{2}}\unicode[STIX]{x1D70C}v_{i}v_{i}+p\right)-\unicode[STIX]{x1D707}v_{i}v_{i,j}\right]n_{j}\,\text{d}A-{\textstyle \frac{1}{2}}\unicode[STIX]{x1D707}\int _{V}v_{i,j}^{2}\,\text{d}V=0,\end{eqnarray}$$ where the pressure 
$p$ is taken to be sum of a Laplace term and a hydrostatic pressure term. Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) initially ignore all the viscous terms in (6.1) and integrated the first term over a cylindrical surface spanning the film to obtain (in our notation)
$$\begin{eqnarray}\unicode[STIX]{x1D712}(r)=\int _{0}^{h}u\left({\textstyle \frac{1}{2}}\unicode[STIX]{x1D70C}(u^{2}+w^{2})+p\right)r\,\text{d}z\approx {\textstyle \frac{1}{2}}C_{3}\unicode[STIX]{x1D70C}ru_{s}^{3}h+C_{1}u_{s}rhp.\end{eqnarray}$$ They then differentiate this expression with respect to 
$r$, and set the derivative 
$\unicode[STIX]{x1D712}^{\prime }$ equal to the viscous terms which they denote as 
$\unicode[STIX]{x1D709}$. This leads to their equation (20), where the viscous terms remain unspecified as 
$\unicode[STIX]{x1D709}$. However, the surface energy term 
$\int \unicode[STIX]{x1D707}v_{i}v_{i,j}n_{j}\,\text{d}S$ evaluated at the free surface 
$S$ is, from (3.6),
$$\begin{eqnarray}\int _{S}\!\unicode[STIX]{x1D707}v_{i}v_{i,j}n_{j}\,\text{d}S\equiv \unicode[STIX]{x1D707}\!\!\int _{S}\{\boldsymbol{u}\,\boldsymbol{\cdot }\,[(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}]+\boldsymbol{n}\,\boldsymbol{\cdot }\,[(\boldsymbol{u}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}]\}\,\text{d}S=\unicode[STIX]{x1D6FE}\!\!\int _{S}\{(\unicode[STIX]{x1D735}_{I}\,\boldsymbol{\cdot }\,\boldsymbol{u})-[(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735}_{I})\boldsymbol{u}]\,\boldsymbol{\cdot }\,\boldsymbol{n}\}\,\text{d}S,\end{eqnarray}$$and includes the effect of surface tension and depends on the surface velocity. In their analysis, Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) did not include the surface energy term.
Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) integrated (6.1) over the inner cylinder of the control volume. Instead, in our analysis we integrate (6.1) for the complete closed annular control volume from 0 to 
$2\unicode[STIX]{x03C0}$, which yields
$$\begin{eqnarray}-{\textstyle \frac{1}{2}}C_{3}\unicode[STIX]{x1D70C}ru_{s}^{3}h|_{r}^{r+\text{d}r}-C_{1}u_{s}rhp|_{r}^{r+\text{d}r}+\unicode[STIX]{x1D6FE}\text{d}\{u_{s}r(1+{h^{\prime }}^{2})^{(1/2)}\}-\unicode[STIX]{x1D70F}_{w}r\,\text{d}r=0.\end{eqnarray}$$ Recognising that 
$u_{s}rh=$ constant, 
$\unicode[STIX]{x1D70F}_{w}$ is the wall friction and dividing (6.4) by 
$\text{d}r$, the limit of 
$h^{\prime }\rightarrow 0$ gives
$$\begin{eqnarray}-C_{3}\unicode[STIX]{x1D70C}ru_{s}^{2}h\frac{\text{d}u_{s}}{\text{d}r}+\unicode[STIX]{x1D6FE}r\frac{\text{d}u_{s}}{\text{d}r}+\unicode[STIX]{x1D6FE}u_{s}-C_{1}u_{s}rh\frac{\text{d}p}{\text{d}r}+\unicode[STIX]{x1D70F}_{w}r\,\text{d}r=0.\end{eqnarray}$$ This equation is equivalent to (in the absence of gravity) (5.4) in Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) which was obtained from a control volume approach and leads to a critical Weber number of order one at the jump (and not 
$O(\unicode[STIX]{x1D6FC})$ as Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) claim).
Including this surface term (6.3) explicitly in 
$\unicode[STIX]{x1D709}$, since it involves 
$u_{s}$ and 
$\unicode[STIX]{x1D6FE}$, will change the result in Duchesne et al. (Reference Duchesne, Andersen and Bohr2019) and leads to the conclusion that surface tension is the dominant force, consistent with the increase in surface area documented in § 3.3.
7 Conclusions
Applying conservation of radial and wall-normal momentum to the flow in an expanding axisymmetric thin film shows that the curvature of the film is singular at a finite radius determined by a critical value of the Weber number. This singularity arises from the wall-normal momentum conservation which implies that 
$\text{d}^{2}h/\text{d}r^{2}\rightarrow \infty$ whereas 
$(\text{d}h/\text{d}r)^{2}\rightarrow 0$ at a finite radius. In § 3.1, we presented a physical mechanism to demonstrate the mechanical and surface energy exchange in steady flow. Physically, in the thin-film limit, 
$h^{\prime \prime }=(1+h^{\prime 2})(\text{d}\unicode[STIX]{x1D6FC}/\text{d}r)\approx \text{d}\unicode[STIX]{x1D6FC}/\text{d}r$, and the singularity implies that this change in interface slope 
$\text{d}\unicode[STIX]{x1D6FC}/\text{d}r$ coincides with the hydraulic jump. Therefore, the mathematical expression in tandem with the physical interpretation indicates a jump. We also show that a similar jump condition can be obtained by applying conservation of radial momentum.
Using conservation of energy Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) showed that the radial velocity gradient is also singular, in this case 
$\text{d}h/\text{d}r\rightarrow \infty$, at a critical Weber number which is numerically slightly different. This radius was identified in experiments as a jump in the flow depth to a thicker and slower flow downstream, and excellent quantitative agreement was found in the predicted and observed values of the jump radius.
There is a small numerical difference between the two predictions of the jump radius 
$R$ given by conservation of energy and momentum, respectively, namely
$$\begin{eqnarray}\frac{R}{R_{ST}}=\left(\frac{1}{f^{\prime }(0)(2\unicode[STIX]{x03C0})^{3}}\frac{C_{2}}{C_{1}^{3}}\right)^{1/4}=0.2705\quad \text{and}\quad \left(\frac{1}{f^{\prime }(0)(2\unicode[STIX]{x03C0})^{3}}\frac{C_{3}}{C_{1}^{3}}\right)^{1/4}=0.2481.\end{eqnarray}$$ These numerical values are obtained from Watson’s similarity profile which, as acknowledged, is only an approximation to the flow in the film, and the predictions are both smaller than the experimentally measured values of 0.289 
$\pm$ 0.015. Since 
$C_{1}$ is the area under the curve 
$f(\unicode[STIX]{x1D702})$ then this will be smaller for the real profile (figure 3), leading to a larger prediction. Also, since 
$0\leqslant \unicode[STIX]{x1D702}\leqslant 1$, 
$C_{2}>C_{3}$, the jump radius estimate from momentum conservation is always smaller than that obtained from energy conservation, suggesting that energy is dissipated in the jump.
Acknowledgements
R.K.B. wishes to thank his PhD supervisor Professor D. I. Wilson for his advice and support.
Declaration of interests
The authors report no conflict of interest.
Appendix A. Scaling relation including gravity and surface tension
Bhagat et al. (Reference Bhagat, Jha, Linden and Wilson2018) provided a scaling relationship for the jump radius 
$R_{ST}$ from the following three conditions: (i) that the radial flow velocity 
$u$ and depth 
$h$ is balanced by viscous drag, 
$u/R\sim \unicode[STIX]{x1D708}/h^{2}$, (ii) continuity 
$uRh\sim Q$ and (iii) the jump is surface tension dominated, which implies at the jump, the Weber number 
$We=\unicode[STIX]{x1D70C}u^{2}h/\unicode[STIX]{x1D6FE}\approx 1$. Further, their theoretical analysis including gravity gave the condition of hydraulic jump to be
$$\begin{eqnarray}\frac{1}{We}+\frac{1}{Fr^{2}}=1,\end{eqnarray}$$ where the Froude number 
$Fr=u/\sqrt{gh}$. Incorporating both gravity and surface tension through the use of (A 1) modifies the scaling relation (2.1) in the form
$$\begin{eqnarray}R\sim R_{ST}\left[\sqrt{\left(\frac{Q_{C}}{Q}\right)^{2}+2\left(\frac{Q_{C}}{Q}\right)}-\left(\frac{Q_{C}}{Q}\right)\right]^{1/4}.\end{eqnarray}$$ As discussed in § 2, this gravitational correction to the radius 
$R_{ST}$ determined by surface tension alone is 5 % for 
$Q=2Q_{C}$ and 23 % for 
$Q=10Q_{C}$.
Appendix B. Self-similar velocity profile
Following Watson’s (Reference Watson1964) analysis, the conventional literature on the hydraulic jump (e.g. Bohr et al. Reference Bohr, Dimon and Putkaradze1993; Bush & Aristoff Reference Bush and Aristoff2003; Kasimov Reference Kasimov2008) assumes a self-similar velocity profile of the form
$$\begin{eqnarray}u(r,z)=u_{s}(r)f(\unicode[STIX]{x1D702}),\quad \unicode[STIX]{x1D702}\equiv \frac{z}{h(r)},\quad 0\leqslant \unicode[STIX]{x1D702}\leqslant 1,\end{eqnarray}$$ where 
$u_{s}$ is the surface velocity and 
$f(0)=0,f(1)=1,f^{\prime }(1)=0$. which further implies (see Watson Reference Watson1964) that
$$\begin{eqnarray}\boldsymbol{u}=u_{s}(r)f(\unicode[STIX]{x1D702})\hat{r}+u_{s}(r)h^{\prime }(r)\unicode[STIX]{x1D702}f(\unicode[STIX]{x1D702})\hat{z}.\end{eqnarray}$$ Here, we will show that the conventional self-similar velocity profiles do not satisfy the zero interfacial shear-stress condition or 
$(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}\neq 0$. We write
$$\begin{eqnarray}(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}=\left(-\text{sin}\,\unicode[STIX]{x1D6FC}\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}r}+\cos \unicode[STIX]{x1D6FC}\frac{\unicode[STIX]{x2202}}{\unicode[STIX]{x2202}z}\right)(u_{s}(r)f(\unicode[STIX]{x1D702})\hat{r}+u_{s}(r)h^{\prime }(r)\unicode[STIX]{x1D702}f(\unicode[STIX]{x1D702})\hat{z})\end{eqnarray}$$or
$$\begin{eqnarray}\displaystyle & & \displaystyle \left[-\text{sin}\,\unicode[STIX]{x1D6FC}\frac{\text{d}u_{s}}{\text{d}r}f(\unicode[STIX]{x1D702})-\sin \unicode[STIX]{x1D6FC}u_{s}\frac{h^{\prime }}{h}\unicode[STIX]{x1D702}f^{\prime }(\unicode[STIX]{x1D702})+\frac{1}{h}\cos \unicode[STIX]{x1D6FC}u_{s}\,f^{\prime }(\unicode[STIX]{x1D702})\right]\hat{r}\nonumber\\ \displaystyle & & \displaystyle \quad +\left[-\text{sin}\,\unicode[STIX]{x1D6FC}h^{\prime }\frac{\text{d}u_{s}}{\text{d}r}\unicode[STIX]{x1D702}f(\unicode[STIX]{x1D702})-\sin \unicode[STIX]{x1D6FC}h^{\prime \prime }u_{s}\unicode[STIX]{x1D702}f(\unicode[STIX]{x1D702})+\sin \unicode[STIX]{x1D6FC}h^{\prime }u_{s}\unicode[STIX]{x1D702}f(\unicode[STIX]{x1D702})\frac{h^{\prime }}{h}\right.\nonumber\\ \displaystyle & & \displaystyle \quad +\left.\sin \unicode[STIX]{x1D6FC}h^{\prime }u_{s}\frac{h^{\prime }}{h}\unicode[STIX]{x1D702}^{2}f^{\prime }(\unicode[STIX]{x1D702})+\cos \unicode[STIX]{x1D6FC}u_{s}\frac{h^{\prime }}{h}f(\unicode[STIX]{x1D702})+\cos \unicode[STIX]{x1D6FC}u_{s}\frac{h^{\prime }}{h}\unicode[STIX]{x1D702}f^{\prime }(\unicode[STIX]{x1D702})\right]\hat{z}.\end{eqnarray}$$ Evaluating 
$(\boldsymbol{n}\boldsymbol{\cdot }\unicode[STIX]{x1D735})\boldsymbol{u}$, equation (B 3) at the free surface with boundary conditions, 
$f(1)=1$ and 
$f^{\prime }(1)=0$ yields
$$\begin{eqnarray}\left[-\text{sin}\,\unicode[STIX]{x1D6FC}\frac{\text{d}u_{s}}{\text{d}r}\right]\hat{r}+\left(-\text{sin}\,\unicode[STIX]{x1D6FC}h^{\prime }\frac{\text{d}u_{s}}{\text{d}r}-\sin \unicode[STIX]{x1D6FC}h^{\prime \prime }u_{s}+\sin \unicode[STIX]{x1D6FC}h^{\prime }u_{s}\frac{h^{\prime }}{h}+\cos \unicode[STIX]{x1D6FC}u_{s}\frac{h^{\prime }}{h}\right)\hat{z}\end{eqnarray}$$or
$$\begin{eqnarray}\displaystyle & & \displaystyle -\frac{h^{\prime }}{(1+h^{\prime 2})^{1/2}}\frac{\text{d}u_{s}}{\text{d}r}\hat{r}+\left(-\frac{h^{\prime 2}}{(1+h^{\prime 2})^{1/2}}\frac{\text{d}u_{s}}{\text{d}r}-\frac{h^{\prime }h^{\prime \prime }u_{s}}{(1+h^{\prime 2})^{1/2}}+\frac{1}{h}\frac{h^{\prime 3}u_{s}}{(1+h^{\prime 2})^{1/2}}\right.\nonumber\\ \displaystyle & & \displaystyle \quad +\left.\frac{1}{h}\frac{h^{\prime }}{(1+h^{\prime 2})^{1/2}}u_{s}\right)\hat{z}.\end{eqnarray}$$ The analysis implies that only for a completely flat film for which 
$h^{\prime }=0$ and
$\boldsymbol{n}=\hat{z}$, 
$(\boldsymbol{n}\,\boldsymbol{\cdot }\,\unicode[STIX]{x1D735})\boldsymbol{u}=0$ can be trivially satisfied. In all other cases the tangential stress is non-zero.
Open access
