scholarly journals Channelized free-surface flow of cohesionless granular avalanches in a chute with shallow lateral curvature

1999 ◽  
Vol 392 ◽  
pp. 73-100 ◽  
Author(s):  
M. WIELAND ◽  
J. M. N. T. GRAY ◽  
K. HUTTER

A series of laboratory experiments and numerical simulations have been performed to investigate the rapid fluid-like flow of a finite mass of granular material down a chute with partial lateral confinement. The chute consists of a section inclined at 40° to the horizontal, which is connected to a plane run-out zone by a smooth transition. The flow is confined on the inclined section by a shallow parabolic cross-slope profile. Photogrammetric techniques have been used to determine the position of the evolving boundary during the flow, and the free-surface height of the stationary granular deposit in the run-out zone. The results of three experiments with different granular materials are presented and shown to be in very good agreement with numerical simulations based on the Savage–Hutter theory for granular avalanches. The basal topography over which the avalanche flows has a strong channelizing effect on the inclined section of the chute. As the avalanche reaches the run-out zone, where the lateral confinement ceases, the head spreads out to give the avalanche a characteristic ‘tadpole’ shape. Sharp gradients in the avalanche thickness and velocity began to develop at the interface between the nose and tail of the avalanche as it came to rest, indicating that a shock wave develops close to the end of the experiments.

2005 ◽  
Vol 127 (3) ◽  
pp. 572-582 ◽  
Author(s):  
Shin Hyung Rhee

The present study is concerned with liquid tank sloshing at low filling level conditions. The volume of fluid method implemented in a Navier–Stokes computational fluid dynamics code is employed to handle the free-surface flow of liquid sloshing. The geometric reconstruction scheme for the interface representation is employed to ensure sharpness at the free surface. The governing equations are discretized by second order accurate schemes on unstructured grids. Several different computational approaches are verified and numerical uncertainties are assessed. The computational results are validated against existing experimental data, showing good agreement. The capability is demonstrated for a generic membrane-type liquefied natural gas carrier tank with a simplified pump tower inside. The validation results suggest that the present computational approach is both easy to apply and accurate enough for more realistic problems.


2015 ◽  
Vol 767 ◽  
pp. 811-841 ◽  
Author(s):  
C. Frederik Brasz ◽  
Craig B. Arnold ◽  
Howard A. Stone ◽  
John R. Lister

AbstractWhen a solid boundary deforms rapidly into a quiescent liquid layer, a flow is induced that can lead to jet formation. An asymptotic analytical solution is presented for this flow, driven by a solid boundary deforming with dimensionless vertical velocity $V_{b}(x,t)={\it\epsilon}(1+\cos x)\,f(t)$, where the amplitude ${\it\epsilon}$ is small relative to the wavelength and the time dependence $f(t)$ approaches 0 for large $t$. Initially, the flow is directed outwards from the crest of the deformation and slows with the slowing of the boundary motion. A domain-perturbation method is used to reveal that, when the boundary stops moving, nonlinear interactions with the free surface leave a remnant momentum directed back towards the crest, and this momentum can be a precursor to jet formation. This scenario arises in a laser-induced printing technique in which an expanding blister imparts momentum into a liquid film to form a jet. The analysis provides insight into the physics underlying the interaction between the deforming boundary and free surface, in particular, the dependence of the remnant flow on the thickness of the liquid layer and the deformation amplitude and wavelength. Numerical simulations are used to show the range of validity of the analytical results, and the domain-perturbation solution is extended to an axisymmetric domain with a Gaussian boundary deformation to compare with previous numerical simulations of blister-actuated laser-induced forward transfer.


Author(s):  
Shin Hyung Rhee

The present study is concerned with the liquid tank sloshing at low filling level conditions. The volume of fluid method implemented in a Navier-Stokes computational fluid dynamics code is employed to handle the free-surface flow of liquid sloshing. A geometric reconstruction scheme for the interface representation is employed to ensure sharpness at the free-surface. The governing equations are discretized by second order accurate schemes on unstructured grids. Several different computational approaches are verified and numerical uncertainties are assessed. The computational results are validated against existing experimental data, showing good agreement. The capability is demonstrated for a generic membrane type LNG carrier tank with a simplified pump tower inside. The validation results suggest that the present computational approach is both easy to apply and accurate enough for more realistic problems.


2016 ◽  
Vol 808 ◽  
pp. 441-468 ◽  
Author(s):  
S. L. Gavrilyuk ◽  
V. Yu. Liapidevskii ◽  
A. A. Chesnokov

A two-layer long-wave approximation of the homogeneous Euler equations for a free-surface flow evolving over mild slopes is derived. The upper layer is turbulent and is described by depth-averaged equations for the layer thickness, average fluid velocity and fluid turbulent energy. The lower layer is almost potential and can be described by Serre–Su–Gardner–Green–Naghdi equations (a second-order shallow water approximation with respect to the parameter $H/L$, where $H$ is a characteristic water depth and $L$ is a characteristic wavelength). A simple model for vertical turbulent mixing is proposed governing the interaction between these layers. Stationary supercritical solutions to this model are first constructed, containing, in particular, a local turbulent subcritical zone at the forward slope of the wave. The non-stationary model was then numerically solved and compared with experimental data for the following two problems. The first one is the study of surface waves resulting from the interaction of a uniform free-surface flow with an immobile wall (the water hammer problem with a free surface). These waves are sometimes called ‘Favre waves’ in homage to Henry Favre and his contribution to the study of this phenomenon. When the Froude number is between 1 and approximately 1.3, an undular bore appears. The characteristics of the leading wave in an undular bore are in good agreement with experimental data by Favre (Ondes de Translation dans les Canaux Découverts, 1935, Dunod) and Treske (J. Hydraul Res., vol. 32 (3), 1994, pp. 355–370). When the Froude number is between 1.3 and 1.4, the transition from an undular bore to a breaking (monotone) bore occurs. The shoaling and breaking of a solitary wave propagating in a long channel (300 m) of mild slope (1/60) was then studied. Good agreement with experimental data by Hsiao et al. (Coast. Engng, vol. 55, 2008, pp. 975–988) for the wave profile evolution was found.


Author(s):  
Aggelos S. Dimakopoulos ◽  
Athanassios A. Dimas

A numerical model is presented for the simulation of the two-dimensional, inviscid, free-surface flow developing by the propagation and breaking of water waves over a flat bottom of steep slope. The simulation is based on the numerical solution of the unsteady, two-dimensional, Euler equations subject to the fully-nonlinear free-surface boundary conditions, the non-penetration condition at the bottom and appropriate inflow and outflow conditions. A boundary-fitted transformation, which includes both the time-dependent free surface and the arbitrary bottom shape, is applied. For the numerical solution of the Euler equations, a two-stage fractional time-step method is employed for the temporal discretization, while a hybrid scheme is used for the spatial discretization. Finite differences are used in the streamwise direction and a pseudo-spectral method in the vertical direction. An absorption zone is placed at the outflow region in order to minimize wave reflection by the outflow boundary. Wave breaking is modeled by a surface roller breaking model, which modifies the dynamic free-surface condition. The simulation results are in very good agreement with available experimental results for the wave propagation and breaking over bottom with slope 1:35. Results, from the simulations over bottom with steeper slopes of 1:15 and 1:10, which generate strong spilling and mild plunging breakers, respectively, are also in very good agreement with available predictions for the breaking depth and wave height. In all cases, a vortex is formed under the breaking wave front and convected in the surf zone.


1995 ◽  
Vol 297 ◽  
pp. 37-47 ◽  
Author(s):  
G. C. Hocking

Accurate numerical solutions to the problem of finding the location of the interface between two unconfined regions of fluid of different density during the withdrawal process are presented. Supercritical flows are considered, in which the interface is drawn directly into the sink. As the flow rate is reduced, the interface enters the sink more steeply, until the solution method breaks down just before the interface enters the sink vertically from above, and becomes flow from the lower layer only. This lower bound on supercritical flow is compared with the upper bound on single-layer (free surface) flow with good agreement.


Author(s):  
K. J. Bai ◽  
J. H. Kyoung ◽  
J. W. Kim

This paper describes a finite element method applied to a nonlinear free surface flow problem for a ship moving in three dimensions. The physical model is taken to simulate the towing tank experimental conditions. The exact nonlinear free-surface flow problem formulated by an initial/boundary value problem is replaced by an equivalent weak formulation. The same problem was considered earlier by Bai, et. al. [1] where some irregularities were observed in the downstream waves and a transom stern ship geometry could not be treated. In the present paper, specifically, three improvements are made from the earlier work. The first improvement is the introduction of the 5-point Chebyshev filtering scheme which eliminates the irregular and saw-toothed waves in the downstream. The second improvement is that now we can treat a transom stern ship geometry. The third improvement is the introduction of a new boundary condition to simulate a dry bottom behind a transom stern ship which is stretched from the free surface to the bottom at a high Froude number. Computations are made for two models. The first model is tested for the generation of the solitons in the upstream and smooth waves in the downstream. The second model is used to compute the generation of a dry bottom behind a transom stern which is one of highly nonlinear phenomena. The results of the first model show a good agreement with previous results for the generation of the solitons. The results of the second model also show a good agreement with the preliminary experimental observation for a dry-bottom, which will be reported in near future. The numerical simulation of the second model can be applied to the local flow behind a sail of a submarine in cruise, a sloshing problem in LNG tankers, and a dam breaking problem. Both computed models are assumed to be in shallow water for simplicity. However, the present numerical method can treat arbitrary water-depth and practical ship geometries.


2013 ◽  
Vol 720 ◽  
pp. 314-337 ◽  
Author(s):  
X. Cui ◽  
J. M. N. T. Gray

AbstractSnow avalanches and other hazardous geophysical granular flows, such as debris flows, lahars and pyroclastic flows, often impact on obstacles as they flow down a slope, generating rapid changes in the flow height and velocity in their vicinity. It is important to understand how a granular material flows around such obstacles to improve the design of deflecting and catching dams, and to correctly interpret field observations. In this paper small-scale experiments and numerical simulations are used to investigate the supercritical gravity-driven free-surface flow of a granular avalanche around a circular cylinder. Our experiments show that a very sharp bow shock wave and a stagnation point are generated in front of the cylinder. The shock standoff distance is accurately reproduced by shock-capturing numerical simulations and is approximately equal to the reciprocal of the Froude number, consistent with previous approximate results for shallow-water flows. As the grains move around the cylinder the flow expands and the pressure gradients rapidly accelerate the particles up to supercritical speeds again. The internal pressure is not strong enough to immediately push the grains into the space behind the cylinder and instead a grain-free region, or granular vacuum, forms on the lee side. For moderate upstream Froude numbers and slope inclinations, the granular vacuum closes up rapidly to form a triangular region, but on steeper slopes both experiments and numerical simulations show that the pinch-off distance moves far downstream.


Sign in / Sign up

Export Citation Format

Share Document