scholarly journals A multilayer shallow model for dry granular flows with the -rheology: application to granular collapse on erodible beds

2016 ◽  
Vol 798 ◽  
pp. 643-681 ◽  
Author(s):  
E. D. Fernández-Nieto ◽  
J. Garres-Díaz ◽  
A. Mangeney ◽  
G. Narbona-Reina

In this work we present a multilayer shallow model to approximate the Navier–Stokes equations with the ${\it\mu}(I)$-rheology through an asymptotic analysis. The main advantages of this approximation are (i) the low cost associated with the numerical treatment of the free surface of the modelled flows, (ii) the exact conservation of mass and (iii) the ability to compute two-dimensional profiles of the velocities in the directions along and normal to the slope. The derivation of the model follows Fernández-Nieto et al. (J. Comput. Phys., vol. 60, 2014, pp. 408–437) and introduces a dimensional analysis based on the shallow flow hypothesis. The proposed first-order multilayer model fully satisfies a dissipative energy equation. A comparison with steady uniform Bagnold flow – with and without the sidewall friction effect – and laboratory experiments with a non-constant normal profile of the downslope velocity demonstrates the accuracy of the numerical model. Finally, by comparing the numerical results with experimental data on granular collapses, we show that the proposed multilayer model with the ${\it\mu}(I)$-rheology qualitatively reproduces the effect of the erodible bed on granular flow dynamics and deposits, such as the increase of runout distance with increasing thickness of the erodible bed. We show that the use of a constant friction coefficient in the multilayer model leads to the opposite behaviour. This multilayer model captures the strong change in shape of the velocity profile (from S-shaped to Bagnold-like) observed during the different phases of the highly transient flow, including the presence of static and flowing zones within the granular column.

Author(s):  
Hisham Elsafti ◽  
Hocine Oumeraci

In this study, the fully-coupled and fully-dynamic Biot governing equations in the open-source geotechFoam solver are extended to account for pore fluid viscous stresses. Additionally, turbulent pore fluid flow in deformable porous media is modeled by means of the conventional eddy viscosity concept without the need to resolve all turbulence scales. A new approach is presented to account for porous media resistance to flow (solid-to-fluid coupling) by means of an effective viscosity, which accounts for tortuosity, grain shape and local turbulences induced by flow through porous media. The new model is compared to an implemented extended Darcy-Forchheimer model in the Navier-Stokes equations, which accounts for laminar, transitional, turbulent and transient flow regimes. Further, to account for skeleton deformation, the porosity and other model parameters are updated with regard to strain of geomaterials. The presented model is calibrated by means of available results of physical experiments of unidirectional and oscillatory flows.


1995 ◽  
Vol 303 ◽  
pp. 215-232 ◽  
Author(s):  
H. M. Badr ◽  
S. C. R. Dennis ◽  
S. Kocabiyik ◽  
P. Nguyen

The transient flow field caused by an infinitely long circular cylinder placed in an unbounded viscous fluid oscillating in a direction normal to the cylinder axis, which is at rest, is considered. The flow is assumed to be started suddenly from rest and to remain symmetrical about the direction of motion. The method of solution is based on an accurate procedure for integrating the unsteady Navier–Stokes equations numerically. The numerical method has been carried out for large values of time for both moderate and high Reynolds numbers. The effects of the Reynolds number and of the Strouhal number on the laminar symmetric wake evolution are studied and compared with previous numerical and experimental results. The time variation of the drag coefficients is also presented and compared with an inviscid flow solution for the same problem. The comparison between viscous and inviscid flow results shows a better agreement for higher values of Reynolds and a Strouhal numbers. The mean flow for large times is calculated and is found to be in good agreement with previous predictions based on boundary-layer theory.


1989 ◽  
Vol 111 (3) ◽  
pp. 545-552 ◽  
Author(s):  
R. Nordmann ◽  
F. J. Dietzen ◽  
H. P. Weiser

The compressible flow in a seal can be described by the Navier-Stokes equations in connection with a turbulence model (k–ε model) and an energy equation. By introducing a perturbation analysis in these differential equations we obtain zeroth order equations for the centered position and first order equations for small motions of the shaft about the centered position. These equations are solved by a finite difference technique. The zeroth order equations describe the leakage flow. Integrating the pressure solution of the first order equations yields the fluid forces and the rotordynamic coefficients, respectively.


2010 ◽  
Vol 644 ◽  
pp. 97-106 ◽  
Author(s):  
B. U. FELDERHOF

The flow of a viscous compressible fluid in a circular tube generated by a sudden impulse at a point on the axis is studied on the basis of the linearized Navier–Stokes equations. A no-slip boundary condition is assumed to hold on the wall of the tube. Owing to the finite velocity of sound the flow behaviour differs qualitatively from that of an incompressible fluid. The flow velocity and the pressure disturbance at any fixed point different from the source point vanish at short time and decay at long times with a t−3/2 power law.


2011 ◽  
Vol 9 (5) ◽  
pp. 1257-1283 ◽  
Author(s):  
S. C. Fu ◽  
R. M. C. So ◽  
W. W. F. Leung

AbstractThe objective of this paper is to seek an alternative to the numerical simulation of the Navier-Stokes equations by a method similar to solving the BGK-type modeled lattice Boltzmann equation. The proposed method is valid for both gas and liquid flows. A discrete flux scheme (DFS) is used to derive the governing equations for two distribution functions; one for mass and another for thermal energy. These equations are derived by considering an infinitesimally small control volume with a velocity lattice representation for the distribution functions. The zero-order moment equation of the mass distribution function is used to recover the continuity equation, while the first-order moment equation recovers the linear momentum equation. The recovered equations are correct to the first order of the Knudsen number(Kn);thus, satisfying the continuum assumption. Similarly, the zero-order moment equation of the thermal energy distribution function is used to recover the thermal energy equation. For aerodynamic flows, it is shown that the finite difference solution of the DFS is equivalent to solving the lattice Boltzmann equation (LBE) with a BGK-type model and a specified equation of state. Thus formulated, the DFS can be used to simulate a variety of aerodynamic and hydrodynamic flows. Examples of classical aeroacoustics, compressible flow with shocks, incompressible isothermal and non-isothermal Couette flows, stratified flow in a cavity, and double diffusive flow inside a rectangle are used to demonstrate the validity and extent of the DFS. Very good to excellent agreement with known analytical and/or numerical solutions is obtained; thus lending evidence to the DFS approach as an alternative to solving the Navier-Stokes equations for fluid flow simulations.


1997 ◽  
Vol 119 (3) ◽  
pp. 556-567 ◽  
Author(s):  
Mihai Arghir ◽  
Jean Freˆne

The paper presents a method to calculate the rotordynamic coefficients of circumferentially-grooved liquid seals operating in centered position and turbulent flow regimes. The method is based on the integration of the averaged Navier-Stokes equations and uses a coordinate transformation proposed by Dietzen and Nordmann (1987). The effect of the coordinate transformation on the components of the stress tensor is included in the first order transport equations. To ensure grid independent solutions, numerical boundary conditions for the first-order velocities were formulated using the logarithmic law. The perturbation of the turbulent viscosity was also considered. A pressure recovery effect at the exit section was included in the first order mathematical model. The method is validated by calculations for straight and circumferentially-grooved seals. Comparisons with experimental and theoretical results show a good agreement for straight seals and for seals with few grooves, and a reasonable agreement for severe industrial cases (high Reynolds numbers and large number of grooves).


Sign in / Sign up

Export Citation Format

Share Document