Continuum simulation of the discharge of the granular silo: a validation test for the mu(I)-visco-plastic flow law

Using both a continuum Navier-Stokes solver, with the mu(I)-flow-law implemented to model the viscous behavior, and the discrete Contact Dynamics algorithm, the discharge of granular silos is simulated in two dimensions from the early stages of the discharge until complete release of the material. In both cases, the Beverloo scaling is recovered. We first do not attempt quantitative comparison, but focus on the qualitative behavior of velocity and pressure at different locations in the flow. A good agreement is obtained in the regions of rapid flows, while areas of slow creep are not entirely captured by the continuum model. The pressure field shows a general good agreement. The evolution of the free surface implies differences, however, the bulk deformation is essentially identical in both approaches. The influence of the parameters of the mu(I)-flow-law is systematically investigated, showing the importance of the dependence on the inertial number I to achieve quantitative agreement between continuum and discrete discharge. The general ability of the continuum model to reproduce qualitatively the granular behavior is found to be very encouraging.Comment: 12 pages, 15 figure

Droplet migration in a Hele-Shaw cell: Effect of the lubrication film on the droplet dynamics

Droplet migration in a Hele--Shaw cell is a fundamental multiphase flow problem which is crucial for many microfluidics applications. We focus on the regime at low capillary number and three-dimensional direct numerical simulations are performed to investigate the problem. In order to reduce the computational cost, an adaptive mesh is employed and high mesh resolution is only used near the interface. Parametric studies are performed on the droplet horizontal radius and the capillary number. For droplets with an horizontal radius larger than half the channel height the droplet overfills the channel and exhibits a pancake shape. A lubrication film is formed between the droplet and the wall and particular attention is paid to the effect of the lubrication film on the droplet velocity. The computed velocity of the pancake droplet is shown to be lower than the average inflow velocity, which is in agreement with experimental measurements. The numerical results show that both the strong shear induced by the lubrication film and the three-dimensional flow structure contribute to the low mobility of the droplet. In this low-migration-velocity scenario the interfacial flow in the droplet reference frame moves toward the rear on the top and reverses direction moving to the front from the two side edges. The velocity of the pancake droplet and the thickness of the lubrication film are observed to decrease with capillary number. The droplet velocity and its dependence on capillary number cannot be captured by the classic Hele--Shaw equations, since the depth-averaged approximation neglects the effect of the lubrication film

Adaptive modelling of long-distance wave propagation and fine-scale flooding during the Tohoku tsunami

The 11 March 2011 Tohoku tsunami is simulated using the quadtree-adaptive Saint-Venant solver implemented within the Gerris Flow Solver. The spatial resolution is adapted dynamically from 250 m in flooded areas up to 250 km for the areas at rest. Wave fronts are tracked at a resolution of 1.8 km in deep water. The simulation domain extends over 73° of both latitude and longitude and covers a significant part of the north-west Pacific. The initial wave elevation is obtained from a source model derived using seismic data only. Accurate long-distance wave prediction is demonstrated through comparison with DART buoys timeseries and GLOSS tide gauges records. The model also accurately predicts fine-scale flooding compared to both satellite and survey data. Adaptive mesh refinement leads to orders-of-magnitude gains in computational efficiency compared to non-adaptive methods. The study confirms that consistent source models for tsunami initiation can be obtained from seismic data only. However, while the observed extreme wave elevations are reproduced by the model, they are located further south than in the surveyed data. Comparisons with inshore wave buoys data indicate that this may be due to an incomplete understanding of the local wave generation mechanisms

The granular silo as a continuum plastic flow: the hour-glass vs the clepsydra

The granular silo is one of the many interesting illustrations of the thixotropic property of granular matter: a rapid flow develops at the outlet, propagating upwards through a dense shear flow while material at the bottom corners of the container remains static. For large enough outlets, the discharge flow is continuous; however, by contrast with the clepsydra for which the flow velocity depends on the height of fluid left in the container, the discharge rate of granular silos is constant. Implementing a plastic rheology in a 2D Navier-Stokes solver (following the mu(I)-rheology or a constant friction), we simulate the continuum counterpart of the granular silo. Doing so, we obtain a constant flow rate during the discharge and recover the Beverloo scaling independently of the initial filling height of the silo. We show that lowering the value of the coefficient of friction leads to a transition toward a different behavior, similar to that of a viscous fluid, and where the filling height becomes active in the discharge process. The pressure field shows that large enough values of the coefficient of friction ($\simeq$ 0.3) allow for a low-pressure cavity to form above the outlet, and can thus explain the Beverloo scaling. In conclusion, the difference between the discharge of a hourglass and a clepsydra seems to reside in the existence or not of a plastic yield stress.Comment: 6 pages, 6 figure

Transition in a numerical model of contact line dynamics and forced dewetting

We investigate the transition to a Landau-Levich-Derjaguin film in forced dewetting using a quadtree adaptive solution to the Navier-Stokes equations with surface tension. We use a discretization of the capillary forces near the receding contact line that yields an equilibrium for a specified contact angle $\theta_\Delta$ called the numerical contact angle. Despite the well-known contact line singularity, dynamic simulations can proceed without any explicit additional numerical procedure. We investigate angles from $15^\circ$ to $110^\circ$ and capillary numbers from $0.00085$ to $0.2$ where the mesh size $\Delta$ is varied in the range of $0.0035$ to $0.06$ of the capillary length $l_c$. To interpret the results, we use Cox's theory which involves a microscopic distance $r_m$ and a microscopic angle $\theta_e$. In the numerical case, the equivalent of $\theta_e$ is the angle $\theta_\Delta$ and we find that Cox's theory also applies. We introduce the scaling factor or gauge function $\phi$ so that $r_m = \Delta/\phi$ and estimate this gauge function by comparing our numerics to Cox's theory. The comparison provides a direct assessment of the agreement of the numerics with Cox's theory and reveals a critical feature of the numerical treatment of contact line dynamics: agreement is poor at small angles while it is better at large angles. This scaling factor is shown to depend only on $\theta_\Delta$ and the viscosity ratio $q$. In the case of small $\theta_e$, we use the prediction by Eggers [Phys. Rev. Lett., vol. 93, pp 094502, 2004] of the critical capillary number for the Landau-Levich-Derjaguin forced dewetting transition. We generalize this prediction to large $\theta_e$ and arbitrary $q$ and express the critical capillary number as a function of $\theta_e$ and $r_m$. An analogy can be drawn between $r_m$ and the numerical slip length.Comment: This version of the paper includes the corrections indicated in Ref. [1

Numerical simulation of spreading drops

We consider a liquid drop that spreads on a wettable surface. Different time evolutions have been observed for the base radius r depending of the relative role played by inertia, viscosity, surface tension and the wetting condition. Numerical simulations were performed to discuss the relative effect of these parameters on the spreading described by the evolution of the base radius r(t) and the spreading time tS. Different power law evolutions r(t) ∝ tⁿ have been observed when varying the parameters. At the early stage of the spreading, the power law t½ (n = 1/2) is observed as long as capillarity is balanced by inertia at the contact line. When increasing the viscosity contribution, the exponent n is found to increase despite the increase of the spreading time. The effect of the surface wettability is observed for liquids more viscous than water. For a small contact angle, the power law t½ is then followed by the famous Tanner law t1/10 once the drop shape has reached a spherical cap

Instability regimes in the primary breakup region of planar coflowing sheets

International audienceThis article investigates the appearance of instabilities in two planar coflowing fluid sheets with different densities and viscosities via experiments, numerical simulation and linear stability analysis. At low dynamic pressure ratios a convective instability is shown to appear for which the frequency of the waves in the primary atomization region is influenced by both liquid and gas velocities. For large dynamic pressure ratios an asymptotic regime is obtained in which frequency is solely controlled by gas velocity and the instability becomes absolute. The transition from convective to absolute is shown to be influenced by the velocity defect induced by the presence of the separator plate. We show that in this regime the splitter plate thickness can also affect the nature of the instability if it is larger than the gas vorticity thickness. Computational and experimental results are in agreement with the predictions of a spatio-temporal stability analysis

Multifluid flows with weak and strong discontinuous interfaces using an elemental enriched space

In a previous paper, the authors presented an elemental enriched space to be used in a finite-element framework (EFEM) capable of reproducing kinks and jumps in an unknown function using a fixed mesh in which the jumps and kinks do not coincide with the interelement boundaries. In this previous publication, only scalar transport problems were solved (thermal problems). In the present work, these ideas are generalized to vectorial unknowns, in particular, the incompressible Navier-Stokes equations for multifluid flows presenting internal moving interfaces. The advantage of the EFEM compared with global enrichment is the significant reduction in computing time when the internal interface is moving. In the EFEM, the matrix to be solved at each time step has not only the same amount of degrees of freedom (DOFs) but also the same connectivity between the DOFs. This frozen matrix graph enormously improves the efficiency of the solver. Another characteristic of the elemental enriched space presented here is that it allows a linear variation of the jump, thus improving the convergence rate, compared with other enriched spaces that have a constant variation of the jump. Furthermore, the implementation in any existing finite-element code is extremely easy with the version presented here because the new shape functions are based on the usual finite-element method shape functions for triangles or tetrahedrals, and once the internal DOFs are statically condensed, the resulting elements have exactly the same number of unknowns as the nonenriched finite elements.Peer ReviewedPreprin