Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. © 2013 AIP Publishing LLC.

Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.

Models for acoustic propagation through turbofan exhaust flows-lined ducts

Turbomachinery noise radiating into the rearward arc is an important problem. This noise is scattered by the trailing edges of the nacelle and the jet exhaust, and interacts with the shear layers between the external flow, bypass stream and jet, en route to the far field. In the past a range of relevant model problems involving semi-infinite cylinders have been solved. However, one limitation of these previous solutions is that they do not allow for the jet nozzle protruding a finite distance beyond the end of the nacelle (or in certain configurations being buried a finite distance upstream). With this in mind, we have used the matrix Wiener-Hopf technique to allow precisely this finite nacelle-jet nozzle separation to be included. We have previously reported results for the case of hard-walled ducts, which requires factorisation of a 2 × 2 matrix. In this paper we extend this work by allowing one of the duct walls, in this case the outer wall of the jet pipe, to be acoustically lined. This results in the need to factorise a 3 × 3 matrix, which is completed by use of a combination of pole-removal and Pad́e approximant techniques. Sample results are presented, investigating in particular the effects of exit plane stagger and liner impedance. Here we take the mean flow to be zero, but extension to nonzero Mach numbers in the core and bypass flow has also been completed. Copyright © 2009 by Nigel Peake & Ben Veitch.