We find solutions that describe the levelling of a thin fluid film, comprising a non-Newtonian power-law fluid, that coats a substrate and evolves under the influence of surface tension. We consider the evolution from periodic and localized initial conditions as separate cases. The particular (similarity) solutions in each of these two cases exhibit the generic property that the profiles are weakly singular (that is, higher-order derivatives do not exist) at points where the pressure gradient vanishes. Numerical simulations of the thin film equation, with either periodic or localized initial condition, are shown to approach the appropriate particular solution.

This paper derives approximate ‘sunset’ similarity solutions for receding plane strain and radially symmetric hydraulic fractures in permeable elastic media close to the point of closure. Local analysis is used to show that a receding hydraulic fracture has a linear aperture asymptote $\hat {w}\sim \hat {s}$ in the fracture tip, where $\hat {s}$ is the distance from the fracture front. Due to the regularity of the linear asymptote, it is possible to determine similarity solutions in the form of power series expansions, which, for integers $N\ge 2$ and values of the radius decay exponent $\gamma =1/N$, can be shown to terminate to yield polynomial solutions for the fracture aperture of degree $N$. Of this countable infinity of polynomial solutions, the final aperture profile as the fracture approaches closure is associated with the second-degree polynomial with $\gamma =1/2$ called the sunset solution. For the reverse time $t^{\prime }$ measured from closure, the sunset solution is characterized by $w\sim t^{\prime }$ and $R\sim t^{\prime 1/2}$. Of all the admissible polynomial similarity solutions, the sunset solution is shown to form an attractor, as $t^{\prime }\rightarrow 0$, for receding hydraulic fractures associated with a wide variety of points in parametric space. Using the sunset solution, it is possible to estimate the duration of recession, assuming that the fracture aperture and radius at the start of recession are given, and determine how it scales with a dimensionless shut-in parameter. As the fracture approaches closure, the term responsible for coupling the elastic force balance and fluid conservation becomes subdominant to the other terms in the lubrication equation, which reduces to a local kinematic relation between the decaying fracture aperture and the leak-off velocity. This fundamental decoupling of dynamics from kinematics results in the sunset solution being dependent on only a single material parameter – namely the leak-off coefficient. This isolation of the leak-off coefficient by the sunset solution opens the possibility to determine this parameter from laboratory or field measurements.

In this paper, we extend the results in the literature for boundary layer flow over a horizontal plate, by considering the buoyancy force term in the momentum equation. Using a similarity transformation, we transform the partial differential equations of the problem into coupled nonlinear ordinary differential equations. We first analyse several special cases dealing with the properties of the exact and approximate solutions. Then, for the general problem, we construct series solutions for arbitrary values of the physical parameters. Furthermore, we obtain numerical solutions for several sets of values of the parameters. The numerical results thus obtained are presented through graphs and tables and the effects of the physical parameters on the flow and heat transfer characteristics are discussed. The results obtained reveal many interesting behaviours that warrant further study of the equations related to non-Newtonian fluid phenomena, especially the shear-thinning phenomena. Shear thinning reduces the wall shear stress.