Shock formation in two-layer equal-density viscous gravity currents. (29th January 2019)

• Record Type: Journal Article
• Authors: Dauck, Tim-Frederik; Box, Finn; Gell, Laura; Neufeld, Jerome A.; Lister, John R.
• Is Part Of: Journal of fluid mechanics. Volume 863(2019)
• Page Start: 730
• Page End: 756
• Publication Date: 2019-01-29
• Subjects: gravity currents, -- lubrication theory; Fluid mechanics -- Periodicals; 532.005
• DOI: 10.1017/jfm.2018.1015
• Languages: English
• ISSNs: 0022-1120

Abstract

Abstract : The flow of a viscous gravity current over a lubricating layer of fluid is modelled using lubrication theory. We study the case of an axisymmetric current with constant influx which allows for a similarity solution, which depends on three parameters: a non-dimensional influx rate ${\mathcal{Q}}$ ; a viscosity ratio $m$ between the lower and upper layer fluid; and a relative density difference $\unicode[STIX]{x1D700}$ . The limit of equal densities $\unicode[STIX]{x1D700}=0$ is singular, as the interfacial evolution equation changes nature from parabolic to hyperbolic. Theoretical analysis of this limit reveals that a discontinuity, or shock, in the interfacial height forms above a critical viscosity ratio $m_{crit}=3/2$, i.e. for a sufficiently less viscous upper-layer fluid. The physical mechanism for shock formation is described, which is based on advective steepening of the interface between the two fluids and relies on the lack of a contribution to the pressure gradient from the interfacial slope for equal-density fluids. In the limit of small but non-zero density differences, local travelling-wave solutions are found which regularise the singular structure of a potential shock and lead to a constraint on the possible shock heights in the form of an Oleinik entropy condition. Calculation of a simplified time-dependent system reveals the appropriate boundary conditions for the late-time similarity solution, which includes a shock at the nose of the current for $m>3/2$ . The numerically calculated similarity solutions compare well to experimental measurements with respect to the predictions of self-similarity, the radial extent and the self-similar top-surface shapes of the current.