^{1}

^{*}

^{2}

^{*}

Waves of finite amplitude on a thin layer of non-Newtonian fluid modelled as a power-law fluid are considered. In the long wave approximation, the system of equations taking into account the viscous and nonlinear effects has the hyper- bolic type. For the two-parameter family of periodic waves in the film flow on a vertical wall the modulation equations for nonlinear wave trains are derived and investigated. The stability criterium for roll waves based on the hyperbolicity of the modulation equations is suggested. It is shown that the evolution of stable roll waves can be described by self-similar solutions of the modulation equations.

Mud flows are frequently encountered in mountainous regions, especially after torrential rains, and often exhibit a series of breaking waves (roll waves). These type of waves can also be observed as an event following volcano eruptions. A report on roll waves can be found in extensive references quoted by Ng & Mei [

The discontinuous waves play an important role in engineering and geophysical processes. Roll waves consist of a periodic pattern of bores separated by continuous profiles of free boundary. The transition from uniform flow to intermittent flow regime is usually tackled by resorting to stability theory. When it is perturbed, a steady flow becomes unstable, if certain criteria are satisfied, and evolves towards wave breaking. Dressler has been the first, who gave the analytical solution for such waves in open channel flows. Based on long wave approximation, Dressler’s theory of roll waves was extended in [

The aim of the paper is to give a nonlinear study on stability of permanent roll waves on a shear thinning fluid in the frame of one-dimensional, unsteady, gradually varied, laminar mud flow with the shear stress being evaluated in a conventional manner. Starting from long waves equations averaged over the normal to the bed [

All results presented herein can be regarded as a generalization to a power law mud fluid in laminar flow regime of non linear stability method alredy applied to Newtonian turbulent flows in open channels (Boudlal & Liapidevskii [

Consider two-dimensional film flows of a non-Newtonian liquid on a vertical wall. The coordinate system is defined as follows: the axis is directed vertically down and axis is horizontal and directed outward of the liquid layer. The longitudinal velocity component is denoted by The boundary layer approximation is assumed to be valid, with a power-law shear stress relation for laminar flows taken in the form

Here is the viscosity coefficient of dimension

and is the flow index. The case n = 1 corresponds to the Newtonian fluid and is the ordinary dynamic viscosity [

the governing equations, which consist of the mass and momentum conservation equations averaged in the ordinate direction, are reduced to the system [

Here h is the layer thickness, is the depth-averaged velocity, is the bottom stress, g is the gravity acceleration, is the momentum flux factor. For a shear-thinning fluid, which will be considered bellow, we have and, consequently,.

In dimensionless variables introduced in [

Equation (2.3) take the form (asterisks are omitted)

Here the reference depth and the reference velocity are expressed through the given flow rate and the viscosity as follows:

Note that Equation (2.5) are hyperbolic with the characteristics

It is shown by linear analysis in [

is unstable. The instability of the uniform steady-state flow also can be easily checked by the Whitham method [

It is clear that (2.8) is satisfied for, since

and.

Analogously to unstable flow regimes in open channel flows, the roll waves or periodic discontinuous solutions, have been constructed for (2.5) [

Consider the travelling waves propagating with a constant velocity. Introducing the variable and assuming the flow being steady in the coordinate system moving with the velocity, Equation (2.5) are reduced to the ODE

where and

Let y and be the critical depth and the critical velocity with It follows from (3.1) that the necessary condition of roll wave existence is

In view of (3.2) we have

Note that

and, consequently, at the critical depth.

Let be fixed and. Equations (3.1)-(3.2) take the form

To find the function in (3.6), the following identity is used:

.

Furthermore, the function can be represented as follows

where, and is a root of the equation at the interval. The continuous function

is positive for. It means that the function vanishes only at two points and for [

.

Now we can construct the two-parameter family of roll waves as follows: for given and we put, ,. Note that and are the conjugate depths, since the Rankine-Hugoniot conditions for discontinuous solutions of (2.5), which are reduced to the relation

or

,(3.8)

are fulfilled. The stability condition for shocks satisfying (3.8) takes the form (Rozhdestvenskii and Janenko [

or

Thus the admissible values of the governing parameters, for which a roll wave exists, belong to the domain

For, a roll wave consists of the smooth part defined by the monotonous solution of (3.6) at the interval and of the jump with the conjugate depths and

It is shown in the previous section that analogously to the roll waves in open channel flows governed by the classic shallow water theory (Whitham [

All averaged quantities can be expressed as functions of and as follows :

Here we have used (3.1)-(3.2), (3.8) for periodic roll waves defined by parameters and.

In view of (4.2) and (3.4) the modulation equations take the form

where

The nonstationary evolution of the governing parameters for a periodic wave train is described by Equation (4.3). We say the roll waves are stable if the modulation Equation (4.3) for corresponding values are hyperbolic. Considering and as new dependent variables instead of and, we can find the characteristics of (4.3) from the quadratic equation

where

Here “·” denotes the differentiation on and.

Equation (4.3) are hyperbolic for. It means that the hyperbolicity domain depends only on the variable or Note that for. It can be shown by numerical calculations that the hyperbolicity interval is rather narrow and it lies in vicinity of. For n = 1 the values of corresponding to the hyperbolicity domain belong to the interval with Therefore, the boundary of the hyperbolicity domain can be found effectively by the following approximation of the function for long waves

Approximate values of characteristics are given by (4.5)-(4.6) with instead of and instead of. In this case the modulation equations simplify considerably since the functions and are linear on:

The dependence for Equations (4.3), (4.8) with n = 1 is shown in

The simple centered waves, i.e. the self-similar solutions of (4.3) depending on the variable, can be found analytically. Equations (4.3), (4.8) for simple

waves take the form

A solution of (4.9) can be represented in the form

Here is a solution of (4.3) for the corresponding family of characteristics. The solution (4.10) can be applied to explain some specific features of roll wave dynamics considered below.

In this section we compare the numerical calculations performed for the model derived by Meza and Balakotaiah (2008) for vertically falling films of Newtonian fluid with numerical solutions of (2.3) for n = 1. Both models are derived for moderate Reynolds numbers of flow, but in contrast to the former model, the surface tension in Equation (2.3) is ignored. Let us rescale (2.3) slightly according to (Meza and Balakotaiah, 2008). For that we introduce the Reynolds number and the Weber number as follows

Here is the surface tension. Consider two models of film flows on a vertical wall depending on and explicitly. The first one is the Shkadov model [

It is clear that for and. Equations (2.5) are equivalent to (5.1) and the systems can be transformed to each other by dilatation of independent variables. The momentum equation derived by Meza and Balakotaiah [

We will compare below numerical solutions of the model (5.2) presented in [

We restrict our attention to the case 1 in [

Numerical calculations using a variant of the Godunov standard scheme illustrate the development of roll waves on the free surface of a thin film flowing on a vertical wall in the frame of the hyperbolic model (5.1). In Figures 2-4 the perturbations of free surface of the liquid layer calculated by (5.1) are comparing with correspond-

ing calculations from [

This algorithm gives the reasonable values for the boundaries of the states “0” and “1” in numerical calculations of the nonstationary problem of roll wave evolution. Therefore, the transition “0” - “1” in Figures 4(a)- (c) can be described by the centered simple wave of (4.9) moving on the left and the roll waves shown in

We have investigated the roll wave generation on vertically falling films by the nonlinear hyperbolic model (2.3), in which the viscous effects are taken into account by assuming that the velocity profile of an exact steady-state solution of the two-dimensional problem can be used also in the one-dimensional wavy flows. The model is based on the first order long-wave approximation of non-Newtonian fluid flows, which shear rate is modeled by a power law [

the periodic discontinuous solutions of (2.3) (roll waves) can be described by two parameters analogously to the roll waves in open channel flows. Moreover, the nonlinear stability of finite amplitude roll waves can be expressed in the terms of the hyperbolicity of modulation Equations (4.3) for the governing parameters of roll waves. Comparison of numerical calculations of roll wave evolution for the models with and without surface tension effects reveal that in spite of the difference in the individual wave shapes the behavior of roll wave packets is very alike for the models. In particular, the transition from “saturated” to “tsunami” waves described in [

The work was supported by the Russian Foundation for Basic Research (Grant No. 10-01-00338) and by the Program for support of leading scientific schools of the Russian Federation (Grant No. NSh-6706.2012.1).