Table of Contents
Journal of Computational Engineering
Volume 2014 (2014), Article ID 201958, 5 pages
Research Article

Navier-Stokes Equation and Computational Scheme for Non-Newtonian Debris Flow

1Institute for Scientific Methodology (ISEM), Via Ugo La Malfa 153, 90146 Palermo, Italy
2Department of Mathematics and Computer Science, University of Salerno, Via Ponte Don Melillo, 84084 Fisciano, Italy
3Department of Engineering, University of Sannio, Piazza Roma 21, 82100 Benevento, Italy

Received 26 November 2013; Revised 26 February 2014; Accepted 27 February 2014; Published 27 March 2014

Academic Editor: Pierre-yves Manach

Copyright © 2014 Ignazio Licata and Elmo Benedetto. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


This paper proposes a computational approach to debris flow model. In recent years, the theoretical activity on the classical Herschel-Bulkley model (1926) has been very intense, but it was in the early 80s that the opportunity to explore the computational model has enabled considerable progress in identifying the subclasses of applicability of different sets of boundary conditions and their approximations. Here we investigate analytically the problem of the simulation of uniform motion for heterogeneous debris flow laterally confined taking into account mainly the geological data and methodological suggestions derived from simulation with cellular automata and grid systems, in order to propose a computational scheme able to operate a compromise between “global” predictive capacities and computing effort.

1. Introduction

The mobility of granular clusters in the upper part of the mountain basins may cause a sliding similar to fluids currents. Obviously this happens in particular conditions of slope and if we have a particular solid-fluid concentration ratio. Inside the mixture that moves, there are several resistances. Despite the strong nonstationarity of currents, it is important to study the conditions of uniform or nearly uniform motion [1]. In fact, the phenomenon is so complex as to require the study of the problem in the simplest possible terms. Moreover, the equation of motion obtained in conditions of uniformity is also used in the simulations of variable motion and the uniform motion is an asymptotic condition where the natural currents tend to it in the absence of geometry variations of the contour and of supply conditions. In the case of currents in equilibrium with the bottom it was observed that the velocity profile along the wall has a zero gradient at the bottom.

The observations in correspondence with the free surface as well as the first velocity measurements made within the mixture showed that the transverse velocity profile exhibits a maximum at the centerline of the section and minimum values in correspondence with the side walls. Concentration measurements have been obtained by some authors for the case of granular mixtures devoid of fine material, in correspondence with the side walls of the channel. These measurements have shown that the concentration increases with the depth of the current, with a roughly linear trend, reaching a maximum value equal to the concentration at the bottom. In [2] the authors study the motion of a solid-liquid mixture under conditions of uniformity and, above all, they analyze the effect of the walls and the concentration distribution within the cross section. They have implemented a numerical code useful for calculating the local speed in the case of a nonhomogeneous Herschel-Bulkley fluid flowing in uniform motion within a channel. The numerical model describes the tendency of the material to sedimentation and erosion, through an iterative mechanism, for which the calculation in cells with speed lower than a critical threshold concentration is equal to that of the material at rest, while in cells with speed higher than the threshold concentration follows an increasing trend with depth.

The application of the model therefore allows simulating the different conditions that may occur for different combinations of the significant variables (bottom slope and concentration of the mixture) both in the presence of equilibrium with the erodible bed and in the presence of an interface with a nonerodible surface. The estimates of the average speed of the current obtained with the model described in [2] and those obtained assuming approximate homogeneous mixture were compared with the experimental measurements relating to a mud flow to the free surface, in uniform motion, with Herschel-Bulkley behavior. For details you can read [2]. The Herschel-Bulkley model, as is known, characterizes an entire class of strongly non-Newtonian fluid behaviors where the relationship between shear stress and shear rate is highly variable from point to point and it is impossible to define a viscosity constant. Traditionally, analytical efforts have focused on approximations able to return to the behavior described in that type of generalized Newtonian fluids. Unfortunately, the results, although mathematically flawless, have proved useful in very few cases of mainly biological interest (blood plasma). The rise of experimental mathematics has shifted the problem. The inability to find an analytical solution has been transformed in the search for useful “approximations classes” [3]. A number of papers [412] have fixed a solution for the Herschel-Bulkley model that has been very valuable for geological simulations. In it appears a constant, about 10−8 sec−1, on which depend the tangent tensions along the direction of flow. The physical reason is fairly straightforward: this value empirically fits a typical situation of the change in viscosity and the convective motions in the natural landslides. We wondered if it is possible to generalize this solution, eliminating the dependence on epsilon, so you can get the same useful evaluations also for situations more and more distant from those in which it developed the original model. We used a suggestion that comes from the simulation with cellular automata and “grid” methods [13, 14], where the global behavior is obtained from simple local rules that describe the interactions between elementary and homogeneous volumes of fluid. This method requires careful consideration of the portion to be considered “primary” and the simulation critically depends on these choices. Change grid is equivalent to a mathematical zoom on the phenomenon, a sort of “renormalization” procedure. Therefore we worked on the idea that an analytical generalization of the model should contain a threshold value associated with the velocity of concentration, used as a key parameter on which depends much of the physical content of phenomenon and it is easily measurable. The results seem to encourage this approach, in other respects “classic” and analytical.

2. Herschel-Bulkley Solution (with Epsilon)

In the study of debris flows the use of a numerical code is crucial, since the Navier-Stokes equation is nonlinear partial differential equation and it is often impossible to solve analytically. Therefore the goal is to find solutions computationally “complacent.” This departs substantially from the traditional mathematical physics approach. One of the key physical assumptions for the calculation code allows a numerical integration of the following conservation of momentum: where is the tensor of tensions, is the velocity, and is the gravitational acceleration. By projecting the equation along the motion axle, in the hypothesis of uniform flow, we get where and are the tangent tensions along the flow direction (x axle).

Since the tangent tension is not defined in the region where the gradient tensor is zero, we have used a relation proposed by Bercovier and Engleman [15] and we have approximated the model of Herschel-Bulkley through the following continuous relations [7]: where is about 10−8 sec−1, empirical value calibrated on geophysical data of landslides, is generalized viscosity coefficient, and is the flow behavior index.

3. The Attempt to Eliminate

As first step, the idea is to consider the following relations: as functions of and ; that is, and , defined in . For simplicity we set and , getting therefore the two functions:

We set in the neighborhood of the origin , so as to get a neighborhood of the origin varying the angular coefficient . We note that this hypothesis is very binding for the stress tensor. In this case we obtain and it is an indeterminate form. Now we suppose that and tend to zero with the same speed and that they are in the first positive quadrant on the plane ; this supposition is permissible since both the partial derivatives and are positive numbers; in fact the velocity is increasing along the directions and . Therefore in the neighborhood of the origin we have and the previous limit becomes

In practice we have modified the numerical code eliminating the parameter and we have set the tensions equal to , when and are smaller than a tolerance value equal to . Therefore we have identified the conditions that allow making assumptions of homogeneity on the behavior of the flow.

4. A New Concentration Law

In [15] the authors have used a linear distribution of concentration

The fact that the concentration grows with the depth has been observed experimentally but it is very difficult to find a relation or . We can consider a Taylor expansion about the point or written in a dimensional form where and are the height and the width of the river bed where the flow is confined. In the previous relation we have to calculate the constants using some conditions that the distribution of concentration has to satisfy. By truncating the series to the fourth order we have where

Global reasonable assumption is that the concentration is constant along the edges of the river bed and equal to maximum packing concentration and, instead, in the points of coordinates , the distribution is linear in accordance with (9). Therefore we have the following conditions:

These conditions allow finding the links among the coefficients . In fact we have

By uniting the two systems, we obtain where we have eliminated all the linearly independent equations!

By resolving the previous system we get

This relation has been modified introducing two new variables and that represent the height and the width of the river bed that are modified in the case that the flow velocity in the points of coordinates and is smaller than a value called velocity of threshold. Therefore the distribution of concentration is defined by

5. Final Results

Initially we have used the new distribution of concentration to integrate the problem through (3) without the parameter following the procedure described in Section 2. Unfortunately we have noted the difficulty of the code to integrate the problem in some cases and the next step has been to resolve the problem through the Bercovier and Engleman equations but using a concentration of the form suggested by relation (22). We have observed that, for very great widths of the river bed, the numerical solution of velocity, to the center, has the tendency to overlap with the velocity numerically calculated for the heterogeneous and not laterally confined flow (we call it analytic solution). By comparing with numerical solution obtained through relation (9), we can see that the new distribution is more neighbour to the analytic solution and therefore this relation gives us more interesting solutions. Besides we underline that, using the new concentration function , it is possible to integrate the problem with values of of the order of 10−8 sec−1 and instead, for such value, the linear distribution does not allow finding solutions. Let us notice that, in the simulation of the experiments, with widths of the river bed of 0.203 m we have observed that the variation of engraves few on the calculation of central velocity. Therefore it is possible to use values of that are not too small and we can get a solution in the smaller possible time. Finally we want to remember that in the precedent numerical codes there were many problems caused by the difference among the mass flows calculated in laboratory and those numerically calculated. Thanks to the new numerical code the differences are reduced, even if we are still distant from the real mass flows.

6. Conclusions

We have proposed a numerical code in the context of uniform motion for heterogeneous debris flow laterally confined. This computational scheme reflects a global description whose value critically depends on (8), which are hypothesis of “local regularity” of viscous convective motions, and it also depends on the distribution of concentration (22), which is here the variable macroscopic key. In other words, the idea is that it is precisely the variation of the resistance to deformation, not linearly proportional to the speed of angular deformation, that characterizes the typical inverse concentration distribution. “Local” hypothesis quite “elastic” allows a simple global description, easy to test computationally.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


  1. T. R. H. Davies, “Large debris flows: a macro-viscous phenomenon,” Acta Mechanica, vol. 63, no. 1–4, pp. 161–178, 1986. View at Publisher · View at Google Scholar · View at Scopus
  2. R. Martino and M. N. Papa, “Simulazione del moto uniforme per una corrente detritica eterogenea in presenza di pareti laterali,” in Proceedings of the 30th Convegno di Idraulica e Costruzioni Idrauliche (IDRA '06), 2006.
  3. M. Jakob and O. Hungr, Debris-Flow Hazards and Related Phenomena, Springer, 2005.
  4. A. Armanini, H. Capart, L. Fraccarollo, and M. Larcher, “Rheological stratification in experimental free-surface flows of granular-liquid mixtures,” Journal of Fluid Mechanics, vol. 532, pp. 269–319, 2005. View at Publisher · View at Google Scholar · View at Scopus
  5. A. Armanini, L. Fraccarollo, L. Guarino, R. Martino, and Y. Bin, “Experimental analysis of the general features of uniform flows over a loose bed,” in Proceedings of the 2nd International Conference on Debris Flow Hazards Mitigation: Mechanics, Prediction, and Assessment, Taipei, Taiwan, August 2000.
  6. R. Martino, “Experimental analysis on the rheological properties of a debris-flow deposit,” in Proceedings of the 3rd International Conference on Debris-Flow Hazards Mitigation: Mechanics, Prediction, and Assessment, Davose, Switzerland, 2003.
  7. R. Martino and M. N. Papa, “Effetto delle pareti nelle correnti detritiche: primi risultati,” in Proceedings of the 29th Convegno Nazionale di Idraulica e Costruzioni Idrauliche, Trento, Italy, 2004.
  8. R. Martino, C. Sabatino, and L. Taglialatela, “An experimental analysis on the rheology of a granular debris flow, part I: collisional stresses,” in Proceedings of the 29th of International Association of Hydraulic Engineering and Research Congress (IAHR '01), Beijing, China, September 2001.
  9. T. C. Papanastasiou, “Flows of materials with yield,” Journal of Rheology, vol. 31, no. 5, pp. 385–404, 1987. View at Publisher · View at Google Scholar · View at Scopus
  10. T. Takahashi, H. Nakagawa, and Y. Satofuka, “Newtonian fluid model for viscous debris-flow,” in Debris Flow Hazard Mitigation: Mechanics, Prediction, and Assessment, Balkema, 2000. View at Google Scholar
  11. Z. Wan and Z. Wang, Hyperconcentrated Flow, Balkema, Rotterdam, The Netherlands, 1994.
  12. K. X. Whipple, “Open-channel flow of Bingham fluids: applications in debris-flow research,” Journal of Geology, vol. 105, no. 2, pp. 243–262, 1997. View at Google Scholar · View at Scopus
  13. G. Iovine, S. Di Gregorio, and V. Lupiano, “Debris-flow susceptibility assessment through cellular automata modeling: an example from 15-16 December 1999 disaster at Cervinara and San Martino Valle Caudina (Campania, southern Italy),” Natural Hazards and Earth System Science, vol. 3, no. 5, pp. 457–468, 2003. View at Publisher · View at Google Scholar · View at Scopus
  14. Y. V. Vassilevski, K. D. Nikitin, M. A. Olshanskii, and K. M. Terekhov, “CFD technology for 3D simulation of large-scale hydrodynamic events and disasters,” Russian Journal of Numerical Analysis and Mathematical Modelling, vol. 27, no. 4, pp. 399–412, 2012. View at Publisher · View at Google Scholar
  15. M. Bercovier and M. Engleman, “A finite element method for incompressible non-Newtonian flows,” Journal Computational Physic, vol. 313, pp. 313–326, 1980. View at Google Scholar