Recent advances in the simulation of particle-laden flows

Marco Ramaioli
Springer Heidelberg, 2014
Article

Résumé

A substantial number of algorithms exists for the simulation of moving particles suspended in fluids. However, finding the best method to address a particular physical problem is often highly non-trivial and depends on the properties of the particles and the involved fluid(s) together. In this report, we provide a short overview on a number of existing simulation methods and provide two state of the art examples in more detail. In both cases, the particles are described using a Discrete Element Method (DEM). The DEM solver is usually coupled to a fluid-solver, which can be classified as grid-based or mesh-free (one example for each is given). Fluid solvers feature different resolutions relative to the particle size and separation. First, a multicomponent lattice Boltzmann algorithm (mesh-based and with rather fine resolution) is presented to study the behavior of particle stabilized fluid interfaces and second, a Smoothed Particle Hydrodynamics implementation (mesh-free, meso-scale resolution, similar to the particle size) is introduced to highlight a new player in the field, which is expected to be particularly suited for flows including free surfaces.

Official source

https://infoscience.epfl.ch/record/204174?ln=fr

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

, ,

A numerical model is presented for the simulation of complex fluid flows with free surfaces. The unknowns are the velocity and pressure fields in the liquid region, together with a function defining the volume fraction of liquid. Although the mathematical formulation of the model is similar to the volume of fluid (VOF) method, the numerical schemes used to solve the problem are different. A splitting method is used for the time discretization. At each time step, two advection problems and a generalized Stokes problem are to be solved. Two different grids are used for the space discretization. The two advection problems are solved on a fixed, structured grid made out of small rectangular cells, using a forward characteristic method. The generalized Stokes problem is solved using a finite element method on a fixed, unstructured mesh. Numerical results are presented for several test cases: the filling of an S-shaped channel, the filling of a disk with core, the broken dam in a confined domain. (C) 1999 Academic Press.

1999

Chargement

The research work reported in the present dissertation is aimed at the analysis of complex physical phenomena involving instabilities and nonlinearities occurring in fluids through state-of-the-art numerical modeling. Solutions of intricate fluid physics problems are devised in two particularly arduous situations: fluid domains with moving boundaries and the high-Reynolds-number regime dominated by nonlinear convective effects. Shear-driven flows of incompressible Newtonian fluids enclosed in cavities of varying geometries are thoroughly investigated in the two following frameworks: transition with a free surface and confined turbulence. The physical system we consider is made of an incompressible Newtonian fluid filling a bounded, or partially bounded cavity. A series of shear-driven flows are easily generated by setting in motion some part of the container boundary. These driven-cavity flows are not only technologically important, they are of great scientific interest because they display almost all physical fluid phenomena that can possibly occur in incompressible flows, and this in the simplest geometrical settings. Thus corner eddies, secondary flows, longitudinal vortices, complex three-dimensional patterns, chaotic particle motions, nonuniqueness, transition, and turbulence all occur naturally and can be studied in the same geometry. This facilitates the comparison of results from experiments, analysis, and computation over the whole range of Reynolds numbers. The flows under consideration are part of a larger class of confined flows driven by linear or angular momentum gradients. This dissertation reports a detailed study of a novel numerical method developed for the simulation of an unsteady free-surface flow in three-space-dimensions. This method relies on a moving-grid technique to solve the Navier-Stokes equations expressed in the arbitrary Lagrangian-Eulerian (ALE) kinematics and discretized by the spectral element method. A comprehensive analysis of the continuous and discretized formulations of the general problem in the ALE frame, with nonlinear, non-homogeneous and unsteady boundary conditions is presented. In this dissertation, we first consider in the internal turbulent flow of a fluid enclosed in a bounded cubical cavity driven by the constant translation of its lid. The solution of this flow relied on large-eddy simulations, which served to improve our physical understanding of this complex flow dynamics. Subsequently, a novel subgrid model based on approximate deconvolution methods coupled with a dynamic mixed scale model was devised. The large-eddy simulation of the lid-driven cubical cavity flow based on this novel subgrid model has shown improvements over traditional subgrid-viscosity type of models. Finally a new interpretation of approximate deconvolution models when used with implicit filtering as a way to approximate the projective grid filter was given. This led to the introduction of the grid filter models. Through the use of a newly-developed method of numerical simulation, in this dissertation we solve unsteady flows with a flat and moving free-surface in the transitional regime. These flows are the incompressible flow of a viscous fluid enclosed in a cylindrical container with an open top surface and driven by the steady rotation of the bottom wall. New flow states are investigated based on the fully three-dimensional solution of the Navier-Stokes equations for these free-surface cylindrical swirling flows, without resorting to any symmetry properties unlike all other results available in the literature. To our knowledge, this study delivers the most general available results for this free-surface problem due to its original mathematical treatment. This second part of the dissertation is a basic research task directed at increasing our understanding of the influence of the presence of a free surface on the intricate transitional flow dynamics of shear-driven flows.

EPFL2007

Chargement

Large Eddy Simulation of Particle Dynamics Inside a Differentially Heated Cavity

Christoph Bosshard

In nuclear safety, most severe accident scenarios lead to the presence of fission products in aerosol form in the closed containment atmosphere. It is important to understand the particle depletion process to estimate the risk of a release of radioactivity to the environment should a containment break occur. As a model for the containment, we use the three-dimensional differentially heated cavity (DHC) problem. DHC is a cubical box with a hot wall and a cold wall on vertical opposite sides. On the other walls of the cube we have adiabatic boundary conditions. For the velocity field the no-slip boundary condition is valid. The flow of the air in the cavity is described by the Boussinesq equations. Complex flow patterns develop and the flow characteristics depend on the non-dimensional Rayleigh and Prandtl numbers. The predominant flow type in the DHC is a turbulent natural convection flow. This work aims at reaching Rayleigh numbers and turbulent levels as high as possible given the available computational resources. The method used to simulate the turbulent flow is the large eddy simulation (LES) where the dynamics of the large eddies is resolved by the computational grid and the small eddies are modelled by the introduction of subgrid scale quantities using a filter function. Numerically, the LES equations are discretized by the spectral element method. Particle trajectories are computed using the Lagrangian particle tracking method, including the relevant forces (drag, gravity, thermophoresis). Four different particle sets with each set containing one million particles and diameters of 10 μm, 15 μm, 25 μm and 35 μm are simulated. The complexity and the size of the large three-dimensional problem requires the use of massively parallel supercomputers. Spectral element methods are naturally suitable for parallelisation by distributing the elements among the processors. For the Lagrangian particle tracking we use a method where equal numbers of particles are assigned to every processor. The flow field is broadcast and every particle processor tracks the assigned particles, a procedure which leads to a perfect load balancing. Simulation results for the flow field and particle sizes from 15 μm to 35 μm at a Rayleigh number of 109 are compared to previous results from a direct numerical simulation. First order statistics of the LES flow fields are in very good agreement with the direct numerical simulation while the agreement of second order moments is fair. Also the turbulent structures associated to the maximum of turbulent kinetic energy production are correctly reproduced. Particle statistics in the LES and the direct numerical simulation were similar and the settling rates practically identical. Contrary to previous particle simulations in LES, it was found that no model was necessary for the influence of the unresolved flow scales on the particle motions. This can be explained, because the important settling mechanism is through gravity and particle deposition at the walls by turbophoresis is negligible.

EPFL2012

Chargement

EPFL2007

Large Eddy Simulation of Particle Dynamics Inside a Differentially Heated Cavity

Christoph Bosshard

EPFL2012