Friction law for dense granular flows: application to the motion of a mass down a rough inclined plane

The problem of the spreading of a granular mass released at the top of a rough inclined plane was investigated. We experimentally measure the evolution of the avalanche from the initiation up to the deposit using a Moir\'e image processing technique. The results are quantitatively compared with the prediction of an hydrodynamic model based on depth averaged equations. In the model, the interaction between the flowing layer and the rough bottom is described by a non trivial friction force whose expression is derived from measurements on steady uniform flows. We show that the spreading of the mass is quantitatively predicted by the model when the mass is released on a plane free of particles. When an avalanche is triggered on an initially static layer, the model fails in quantitatively predicting the propagation but qualitatively captures the evolution.Comment: 19 pages, 10 figures, to be published in J. Fluid Mec

Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity

The aim of this paper is to compare a hyperelastic with a hypoelastic model describing the Eulerian dynamics of solids in the context of non-linear elastoplastic deformations. Specifically, we consider the well-known hypoelastic Wilkins model, which is compared against a hyperelastic model based on the work of Godunov and Romenski. First, we discuss some general conceptual differences between the two approaches. Second, a detailed study of both models is proposed, where differences are made evident at the aid of deriving a hypoelastic-type model corresponding to the hyperelastic model and a particular equation of state used in this paper. Third, using the same high order ADER Finite Volume and Discontinuous Galerkin methods on fixed and moving unstructured meshes for both models, a wide range of numerical benchmark test problems has been solved. The numerical solutions obtained for the two different models are directly compared with each other. For small elastic deformations, the two models produce very similar solutions that are close to each other. However, if large elastic or elastoplastic deformations occur, the solutions present larger differences.Comment: 14 figure

Modeling Shallow Water Flows on General Terrains

A formulation of the shallow water equations adapted to general complex terrains is proposed. Its derivation starts from the observation that the typical approach of depth integrating the Navier-Stokes equations along the direction of gravity forces is not exact in the general case of a tilted curved bottom. We claim that an integration path that better adapts to the shallow water hypotheses follows the "cross-flow" surface, i.e., a surface that is normal to the velocity field at any point of the domain. Because of the implicitness of this definition, we approximate this "cross-flow" path by performing depth integration along a local direction normal to the bottom surface, and propose a rigorous derivation of this approximation and its numerical solution as an essential step for the future development of the full "cross-flow" integration procedure. We start by defining a local coordinate system, anchored on the bottom surface to derive a covariant form of the Navier-Stokes equations. Depth integration along the local normals yields a covariant version of the shallow water equations, which is characterized by flux functions and source terms that vary in space because of the surface metric coefficients and related derivatives. The proposed model is discretized with a first order FORCE-type Godunov Finite Volume scheme that allows implementation of spatially variable fluxes. We investigate the validity of our SW model and the effects of the bottom geometry by means of three synthetic test cases that exhibit non negligible slopes and surface curvatures. The results show the importance of taking into consideration bottom geometry even for relatively mild and slowly varying curvatures

Three-dimensional numerical simulation of magnetohydrodynamic-gravity waves and vortices in the solar atmosphere

With the adaptation of the FLASH code we simulate magnetohydrodynamic-gravity waves and vortices as well as their response in the magnetized three-dimensional (3D) solar atmosphere at different heights to understand the localized energy transport processes. In the solar atmosphere strongly structured by gravitational and magnetic forces, we launch a localized velocity pulse (in horizontal and vertical components) within a bottom layer of 3D solar atmosphere modelled by initial VAL-IIIC conditions, which triggers waves and vortices. The rotation direction of vortices depends on the orientation of an initial perturbation. The vertical driver generates magnetoacoustic-gravity waves which result in oscillations of the transition region, and it leads to the eddies with their symmetry axis oriented vertically. The horizontal pulse excites all magnetohydrodynamic-gravity waves and horizontally oriented eddies. These waves propagate upwards, penetrate the transition region, and enter the solar corona. In the high-beta plasma regions the magnetic field lines move with the plasma and the temporal evolution show that they swirl with eddies. We estimate the energy fluxes carried out by the waves in the magnetized solar atmosphere and conclude that such wave dynamics and vortices may be significant in transporting the energy to sufficiently balance the energy losses in the localized corona. Moreover, the structure of the transition region highly affects such energy transports, and causes the channelling of the propagating waves into the inner corona.Comment: 11 Pages, 12 Figures, Accepted for the publication in MNRA

A high-order nonconservative approach for hyperbolic equations in fluid dynamics

It is well known, thanks to Lax-Wendroff theorem, that the local conservation of a numerical scheme for a conservative hyperbolic system is a simple and systematic way to guarantee that, if stable, a scheme will provide a sequence of solutions that will converge to a weak solution of the continuous problem. In [1], it is shown that a nonconservative scheme will not provide a good solution. The question of using, nevertheless, a nonconservative formulation of the system and getting the correct solution has been a long-standing debate. In this paper, we show how get a relevant weak solution from a pressure-based formulation of the Euler equations of fluid mechanics. This is useful when dealing with nonlinear equations of state because it is easier to compute the internal energy from the pressure than the opposite. This makes it possible to get oscillation free solutions, contrarily to classical conservative methods. An extension to multiphase flows is also discussed, as well as a multidimensional extension

From traffic and pedestrian follow-the-leader models with reaction time to first order convection-diffusion flow models

In this work, we derive first order continuum traffic flow models from a microscopic delayed follow-the-leader model. Those are applicable in the context of vehicular traffic flow as well as pedestrian traffic flow. The microscopic model is based on an optimal velocity function and a reaction time parameter. The corresponding macroscopic formulations in Eulerian or Lagrangian coordinates result in first order convection-diffusion equations. More precisely, the convection is described by the optimal velocity while the diffusion term depends on the reaction time. A linear stability analysis for homogeneous solutions of both continuous and discrete models are provided. The conditions match the ones of the car-following model for specific values of the space discretization. The behavior of the novel model is illustrated thanks to numerical simulations. Transitions to collision-free self-sustained stop-and-go dynamics are obtained if the reaction time is sufficiently large. The results show that the dynamics of the microscopic model can be well captured by the macroscopic equations. For non--zero reaction times we observe a scattered fundamental diagram. The scattering width is compared to real pedestrian and road traffic data

High-Order Unstructured Lagrangian One-Step WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems: Applications to Compressible Multi-Phase Flows

In this article we present the first better than second order accurate unstructured Lagrangian-type one-step WENO finite volume scheme for the solution of hyperbolic partial differential equations with non-conservative products. The method achieves high order of accuracy in space together with essentially non-oscillatory behavior using a nonlinear WENO reconstruction operator on unstructured triangular meshes. High order accuracy in time is obtained via a local Lagrangian space-time Galerkin predictor method that evolves the spatial reconstruction polynomials in time within each element. The final one-step finite volume scheme is derived by integration over a moving space-time control volume, where the non-conservative products are treated by a path-conservative approach that defines the jump terms on the element boundaries. The entire method is formulated as an Arbitrary-Lagrangian-Eulerian (ALE) method, where the mesh velocity can be chosen independently of the fluid velocity. The new scheme is applied to the full seven-equation Baer-Nunziato model of compressible multi-phase flows in two space dimensions. The use of a Lagrangian approach allows an excellent resolution of the solid contact and the resolution of jumps in the volume fraction. The high order of accuracy of the scheme in space and time is confirmed via a numerical convergence study. Finally, the proposed method is also applied to a reduced version of the compressible Baer-Nunziato model for the simulation of free surface water waves in moving domains. In particular, the phenomenon of sloshing is studied in a moving water tank and comparisons with experimental data are provided

The GeoClaw software for depth-averaged flows with adaptive refinement

Many geophysical flow or wave propagation problems can be modeled with two-dimensional depth-averaged equations, of which the shallow water equations are the simplest example. We describe the GeoClaw software that has been designed to solve problems of this nature, consisting of open source Fortran programs together with Python tools for the user interface and flow visualization. This software uses high-resolution shock-capturing finite volume methods on logically rectangular grids, including latitude--longitude grids on the sphere. Dry states are handled automatically to model inundation. The code incorporates adaptive mesh refinement to allow the efficient solution of large-scale geophysical problems. Examples are given illustrating its use for modeling tsunamis, dam break problems, and storm surge. Documentation and download information is available at www.clawpack.org/geoclawComment: 18 pages, 11 figures, Animations and source code for some examples at http://www.clawpack.org/links/awr10 Significantly modified from original posting to incorporate suggestions of referee