Smooth, cusped and sharp shock waves in a one-dimensional model of a microfluidic drop ensemble
Authors: Ramos J.I. 
Published: May 2021

Purpose
The purpose of this paper is to determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of microfluidic droplet ensembles, water flow in unsaturated flows, infiltration, etc., as functions of the powers of the convection and diffusion fluxes and upstream boundary condition; to study numerically the evolution of the wave for two different initial conditions; and to assess the accuracy of several finite difference methods for the solution of the degenerate, nonlinear, advection--diffusion equation that governs the model.

Design/methodology/approach
The theory of ordinary differential equations and several explicit, finite difference methods that use first- and second-order, accurate upwind, central and compact discretizations for the convection terms are used to determine the analytical solution for steadily propagating waves and the evolution of the wave fronts from hyperbolic tangent and piecewise linear initial conditions to steadily propagating waves, respectively. The amplitude and phase errors of the semi-discrete schemes are determined analytically and the accuracy of the discrete methods is assessed.

Findings
For non-zero upstream boundary conditions, it has been found both analytically and numerically that the shock wave is smooth and its steepness increases as the power of the diffusion term is increased and as the upstream boundary value is decreased. For zero upstream boundary conditions, smooth, cusped and sharp shock waves may be encountered depending on the powers of the convection and diffusion terms. For a linear diffusion flux, the shock wave is smooth, whereas, for a quadratic diffusion flux, the wave exhibits a cusped front whose left spatial derivative decreases as the power of the convection term is increased. For higher nonlinear diffusion fluxes, a sharp shock wave is observed. The wave speed decreases as the powers of both the convection and the diffusion terms are increased. The evolution of the solution from hyperbolic tangent and piecewise linear initial conditions shows that the wave back adapts rapidly to its final steady value, whereas the wave front takes much longer, especially for piecewise linear initial conditions, but the steady wave profile and speed are independent of the initial conditions. It is also shown that discretization of the nonlinear diffusion flux plays a more important role in the accuracy of first- and second-order upwind discretizations of the convection term than either a conservative or a non-conservative discretization of the latter. Second-order upwind and compact discretizations of the convection terms are shown to exhibit oscillations at the foot of the wave’s front where the solution is nil but its left spatial derivative is largest. The results obtained with a conservative, centered second--order accurate finite difference method are found to be in good agreement with those of the second-order accurate, central-upwind Kurganov--Tadmor method which is a non-oscillatory high-resolution shock-capturing procedure, but differ greatly from those obtained with a non-conservative, centered, second-order accurate scheme, where the gradients are largest.

Originality/value
A new, one-dimensional model for microfluidic droplet transport, water flow in unsaturated flows, infiltration, etc., that includes high-order convection fluxes and degenerate diffusion, is proposed and studied both analytically and numerically. Its smooth, cusped and sharp shock wave solutions have been determined analytically as functions of the powers of the nonlinear convection and diffusion fluxes and the boundary conditions. These solutions are used to assess the accuracy of several finite difference methods that use different orders of accuracy in space, and different discretizations of the convection and diffusion fluxes, and can be used to assess the accuracy of other numerical procedures for one-dimensional, degenerate, convection--diffusion equations.