Destabilization of Immersed Dense Granular Material Submitted to Localized Fluidization: An Experimental and Numerical Study

An alternative experimental approach and a numerical analysis for the study of destabilization by localized fluidization of an immersed dense granular material are presented. To visualize the evolutions of the internal structure of the granular medium, the hydrogel beads, composed of about 99% of water and having substantially the same refraction indexes, are used as solid phase. A LED lighting system is used in place of a laser lighting system. As a result, the optical access restriction of porous structure is removed. A real economic alternative for the experimental study of fluid-grain coupling during destabilization by localized fluidization of a granular material is created. /e experimental phenomenology presented in the literature is verified: the system passes successively through three different stationary regimes: static regime, fluidized cavity regime, and fluidized chimney regime. Some restrictions of using hydrogel beads as particles in the study of liquid-solid interaction are also discussed.


Introduction
In general, fluidization consists of the setting in movement of a package of particles under the effect of a fluid flow. Several techniques have been used in the past for the experimental and numerical study of the hydrodynamic coupling between the phases present in the fluidization. Among them, we can single out the experimental works, by Zoueshthiagh and Merlen [1] and Philippe and Badiane [2], and some numerical works, including for Cui et al.  [3][4][5][6][7]. e study of the destabilization by localized fluidization of a granular material as well as the study of the behavior of the flows in porous media necessitated a wide use of several optical techniques: the iso-index technique and the laser-induced fluorescence (PLIF), the laser Doppler velocimetry (LDA), the particle tracking velocimetry (PTV), and the particle image velocimetry (PIV) have been extensively used [2,[8][9][10].
ese investigation methods strongly depend on the quality of the refractive index adequacy of the solid and liquid phases [11,12]. Many researchers have reported that, in order to visualize the evolutions of the internal structure of the granular medium without having recourse to a two-dimensional geometry, the interstitial fluid must be adapted to examples of mixtures of silicon oil, glycerol, or any other mineral oil [9,10]. However, it is possible that the optical properties of these fluid mixtures change over time during the experiment (sensitivity of their rheological properties to temperature) [13]. Although the use of these materials has given remarkable results in the different works, the maintenance of a constant visualization through the granular materials remains the most difficult task in the experimental works. In the ultimate concern to visualize the movement of particles of the granular structure inside the porous medium, Philippe and Badiane used the iso-index technique and laser-induced fluorescence (PLIF) [2]. A vertical laser sheet centered on the injection zone (s) illuminates the iso-index medium thus created. A possible replacement, of the laser system, which is expensive and cumbersome and for security problems, by LED lighting, has also been the subject of several studies in the past [14][15][16]. Moreover, some researchers have reported that hydrogel beads are also suitable as a solid phase for visualization through granular materials. Byron and Variano [17] in their studies of the interaction of small solid bodies with a fluid medium used a solid hydrogel material which had the same index of refraction as water. Also, to measure flow velocities around solids, Weitzman et al. [18] used solid hydrogel materials. Finally, Jay et al. [19] have successfully used aqua beads to model multiphase flows. e success of these different studies took us to try the use of hydrogel beads and LED lighting for the experimental study of destabilization by localized fluidization of an immersed dense granular medium. ese hydrogel beads are cross-linked materials that absorb large amounts of solvent without dissolving.
ese abilities make them unique materials [20,21]. According to the design of their matrices, the hydrogels can, by absorption, change their volume in response to stimuli, such as the temperature, the quality of the solvent, the pH, and the electric field [22,23].
An alternative experimental approach is proposed in this work. To visualize the evolutions of the internal structure of the granular medium, the hydrogels, composed of approximately 99% of water and having substantially the same refractive index as water, are used as solid phase. A threedimensional model has been developed using a numerical simulation based on the commercial software of discrete elements method (EDEM) to describe the mechanical behavior of the grain stack and a continuous model, based on the CFD ANSYS-FLUENT code, to solve the equation of fluid dynamics in confined geometry, more precisely in the pore space between grains.

Experimental Setup
e Porous medium's arrangement was done by a set of criteria (adjustment of the optical indices) with the ultimate goal of visualizing changes in the internal structure of the granular medium. We used hydrogel beads, which are denser than water when they are swollen, available commercially. e advantage of these hydrogel beads on the glass beads is that the optical access restrictions are eliminated.
is also eliminates the use of complex optical techniques (need for mineral oil blends to match the refractive indices) to visualize grain motions within the granular medium. is possibility of replacing a fixed medium of transparent glass beads immersed in liquid containing mixtures of mineral oils by simple hydrogel beads in water constitutes one of the innovations of this work. e hydrogel beads, with an initial diameter of 1 mm, contain absorbent polymers that expand when hydrated with a solvent and can reach a diameter of 8 to 9 mm if the solvent is purified water. e diameter can also be of the order of 4 to 5 mm if the solvent is salted water since their final size is very dependent on the amount of ions in the salted water. ey can contain about 99% of water if the solvent is purified water and therefore have the same index of refraction as water.
In addition, Harshani et al. [13] have shown that the diameter of these beads is kept constant for low flow rates, and that they can be used for at least three months without size and shape degradation when stored in the solvent.
A parallelepipedic cell (1000 × 30 × 600 mm) made in acrylic with transparent and good optical clarity walls was used as the main device. A series of five circular openings (5 mm, 10 mm, 15 mm, 10 mm, and 5 mm of diameter), spaced 45 mm each other, are drilled at the base of the cell. In order to generate a uniform, constant, and vertical flow into the cell, a single gear pump (VGS 24V OEM), a flowmeter, and a water tank were used as indicated in Figure 1. To visualize the grain movements inside the granular structure, we opted for an environment enlightened by a thin light plate from a high power LED module, to create a vertical sheet centered on the zone or zones of injection, which helps to improve the quality of captured images. A fast camera, placed in front of the cell, perpendicular to the plane of the LED sheet was used to capture the illuminated zones.
At first, the cell is filled with liquid. e hydrogel beads, introduced into the cell, are arranged in bulk and packed at the bottom of the cell according to the desired stacking height. It is therefore a starting sample that is well reproducible for other initial arrangements of the experimental phase.

Simulation Procedures: Methodology and Simulation Details
A three-dimensional model, based on the coupling of two particulate methods EDEM and the CFD-FLUENT code, has been developed to simulate particle-fluid interactions in the destabilization of immersed dense granular media. EDEM, used for describing the solid particle dynamics, is a commercial software of Discrete Element Method (DEM) expanded for taking into account the bulk material simulation. e continuous model, based on the Eulerian multiphase models code in CFD-FLUENT, is used for solving the fluid dynamics equation. Figure 2 provides an overview of EDEM-FLUENT coupling process flow. e present study, which concerns the case of dense granular materials immersed, takes into account the approach of deformable bodies in contact. is approach allows for multiple particle interactions (including possible local deformation at the point of contact that is frequent in dense fluidized beds).
Particles movement is tracked individually according to Newton's motion laws: where m p is the particle mass, u p is the particle velocity, F P is the particle gravity force, F D is the fluid-particle drag force, F C is the contact force, ω p is the particle angular velocity, I p 2 Journal of Spectroscopy is the particle inertia momentum, and T p is the applied resulting torque on the particle. e drag force represents the energy dissipation resulting from the effects of the fluid on the particles during the movement. It is expressed by where V p is the volume of particle and β is the interphase transfer term. e particle gravity force is expressed by e soft sphere model, proposed by Cundal and Strack [24] and used for the first time by Tsuji et al. [25], was adopted in this study to solve all the problems of particle deformation and collision of each particle. Forces are studied as the three models (spring, dashpots, and a slider). ese particle-particle and particle-wall contact forces are calculated as a function of spring constants and damping coefficients.
e following enumerated equations are used to account for particle overlaps in the model of soft spheres during particle collisions: F Ci is the net contact force acting on a particle i because of its neighbors, T pi is the couple acting on a particle i because of its neighbors, and f n,ij and f t,ij are the the normal and tangential forces on particle i by particle j.
where δ n,ij and δ t,ij express the normal and tangential interpenetration of the two particles i and j, κ n and κ t are, respectively, the stiffness at normal and tangential contact, σ n and σ t are, respectively, the normal and tangential viscous damping coefficients, and the normal and tangential relative velocities v n,ij and v t,ij are given by the following relations:  e friction slider limits the tangential force governed by Coulomb's friction law [26].
Water phase is modeled in CFD-FLUENT as a continuous fluid phase. e equations of continuity and momentum conservation are solved for each time step related to the fluid phase: e continuum equation for the q phase is e momentum balance for the q phase is Expression ρ q is the q phase density, g is the acceleration of gravity, u → p and u → q are the phase velocities, and τ q is the q phase stress-strain tensor given by where μ q and λ q the q phase's shear and bulk viscosity. k sl is the interphase momentum exchange coefficient. For granular flows, momentum exchange between phases is based on the fluid-solid exchange coefficients value. In CFD-ANSYS Fluent, the interphase exchange coefficient models are empirically based. e Gidaspow model, which is a combination of the Wen and Yu model [27] and the Ergun equation [28], is used in this study. is model is recommended for dense fluidized beds.
where μ l is the water dynamic viscosity, d s is the particle diameter, and C D is the drag force coefficient of a single particle. C D is expressed by Reynolds number is as follows: More detailed information on Gidaspow model can be found in [29]. e fluid phase is, at the start, analyzed by the commercial CFD ANSYS-FLUENT code using the solutions of the Navier-Stokes equation and the standard k-epsilon turbulence model. For a given time step, iteration of the simulation is performed in FLUENT until convergence is reached. en, EDEM calculates body forces acting on each particle (gravity, fluid drag, electrostatic) using data extracted from fluid mesh cells and eventually determines updated particle position, acceleration, and velocity. ese different coordinates of the particles thus updated are transferred to the CFD solver ANSYS-FLUENT and the forces exerted on the fluid by the particles are introduced into the fluid through a series of momentum source terms.
e simulations were carried out on a granular of silicon-si assembly which constituted the discrete phase provided by EDEM whereas the liquid water of kinematic viscosity ] f � 2.0 10-6 m 2 s −1 provided by FLUENT constituted the continuous phase. e laminar flow option and the double precision mode were selected in CFD-FLUENT. e particles of spherical form had a density of 2500 kg m −3 and an average size of 1.610 −3 m while the density of the water is 1000 kg m −3 . In order to certify the experimental observations, a parallelepipedic cell of 600 × 30 × 1000 mm constituted the computation domain for the simulations work. Figure 3 shows the computing domain meshed using the Integrated Computer Engineering and Manufacturing (ICEM) code and the simulation domain with EDEM particles. Sides of the mesh model were specified as wall boundary conditions, holes in the bottom were (Figure 3(c)) specified as velocity inlet according to corresponding simulation phase, and the upper surface was specified as an outflow boundary. e selection of the turbulence model affects the analysis because each model uses a different set of boundary conditions. e Reynolds-Averaged Naiver-Stokes (RANS) simulates turbulence in flow as a result of averaging the steady-state and dynamic flow variables. We adopted this model because of its simplicity and ease of understanding. Accurate fine grid resolution is essential to detect the smallest vortex formation in the turbulent flow; therefore high computing capacity is desired. And it is shown that, with generally more fine grids, all iterative solution methods converge with a good accuracy. e computation may require significantly more time [30]. Moreover, to capture the particles behavior accurately, a smaller integration time step is also essential. e computational domain consists of a total of 2979507 finer grids with a total of 594217 nodes. In addition, a time step of 10 −4 s is chosen in this work, which is supposed to be certainly small enough to produce reasonably accurate results. e validation of this numerical model is also based on the visual observations of dynamic process during the experimental work. Moreover, CFD-DEM two-way coupling have been more used in recent years to deal with the gas-particle coupling in two-phase flow systems with accurate results [31].
Furthermore, the direction of acceleration of gravity (g � 9.81 m.s −2 ) was evaluated as being opposite to the direction of flux entering the domain entrance. en, we imposed a flow velocity U greater than the fluidization threshold, at the base of the granular layer, and kept constant during the duration of the simulation.

Results and Discussion
In this work, the collapse or partial rupture of a civil engineering structure in the presence of water infiltration by fluidization is studied in the laboratory through a model based on a granular stack subjected to an ascending liquid flow. Also, a three-dimensional model, based on the coupling of two particulate methods EDEM and the CFD-FLUENT code, has been developed to realistically report the predominant physical mechanisms involved in our experimental studies. is study of the fluid-grain interaction is motivated by the fact that the resulting phenomena are sometimes at the origin of the appearance of hydromechanical instabilities. e fluid exerts aerodynamic or hydrodynamic forces on the grain interface, which can then move or deform. is displacement or deformation of the grains also affects the flow field locally. A better understanding of fluid-grain coupling and its role in the development of hydro mechanical instabilities can help to respond to a concrete problem of the safety of hydraulic structures, in which the fluidization located at the foot appears as an initiating mechanism of regressive erosion. We looked at the case of localized fluidization of an immersed grain layer with the single nozzle located at the center of the cell. For different flow rates through an injection diameter of d i � 15 mm, we observe the growth of the hydromechanical instability which is created at the level of the injection nozzle and propagates within the granular structure up to the surface of the bed.
When the injection speed U is regularly increased, the phenomenology demonstrated in the experiments [1,2] is perfectly highlighted. e bed is described successively by three different regimes: (i) e static regime, which corresponds to the low injection velocities; there is no deformation of the granular stack or significant movement of grains all along the bed. e stack behaves like a fixed porous medium.
(ii) e regime of the fluidized cavity corresponding to the sufficient injection velocities: we observe not only a significant grains movement just above the injection nozzle but also a local deformation of the stack at the surface of the injection zone.
(iii) e fluidized chimney regime corresponding to higher injection velocities: we observe a local fluidization throughout the thickness of the bed and the presence of a vertical chimney above the injection zone.
Indeed, an imbalance of the contact forces induced by the flow on the grains causes the destabilization phenomenon. When the induced forces are more significant than the apparent weight of the granular stack, there can be fluidization with large relative movements of grains in interior and a more or less regular expansion of the whole material. e lower surface was gradually destabilized thus releasing a shower of grains that gradually fed the fluidized zone. e size of the ascending grain layer in solid translation motion decreased regularly while the fluidized layer increased. ere is also fracturing of the bed's injection zone surface. In fact, the internal structure of the granular stack is reorganized by retroacting on the interstitial flow thus leading to the opening of a preferential path of low hydraulic resistance for the liquid flow and the modification of the dynamic behavior of the grains.
More information about the physical phenomenon involved in the three cited regimes can be found in [1].
To report more on the situation in which the stresses exerted by a liquid flow in a submerged granular medium are sufficient to put in motion certain grains, simulation works based on the coupling of two particular methods (EDEM and CFD ANSYS FLUENT) have been carried out. e same dynamic of growth of the fluidized zone in an initially static immersed granular stack was found with the simulations work. Figures 4-6 illustrate the three different stationary regimes observed during the experimental and the numerical works for a stack of initial height H 0 � 400 mm and with an injection diameter d i � 15: the static regime, the fluidized cavity regime, and the fluidized chimney regime. Moreover, by observing the behavior of the bed during the experimental phase, it is clear that the numerical observations show a good agreement with the experimental observations. However, other more in-depth studies devoted to these numerical and experimental observations are in progress and will be the subject of further submissions in the near future.

Conclusion
In natural environments, it is common to observe particular cases of vertical localized fluidization, a phenomenon which manifests itself as a hydraulic flow of infiltration into the foundations of the structure in Earth and can cause an overall lifting and possibly a hydraulic fracturing. e present work consists of experimental and numerical study of the destabilization by localized fluidization of an immersed dense granular material through the development of a new PIV approach able to correctly reproduce the interactions between fluid phase and solid phase. e aim was to force an ascending fluid through a bed of particles, which remain fixed until the flow reaches a certain flow rate, called the minimum flow of fluidization. Beyond this flow, when the drag force of the fluid completely offsets the effective weight of the particles, they are then fluidized. It seeks to understand the concrete problem of the safety of hydraulic structures, in which the localized fluidization at the foot appears as a mechanism initiating regressive erosion, one of the internal erosion modes. Restriction of the optical access to the porous structure was removed using hydrogel beads as solid particles and water as the fluid having substantially the same refractive index. A LED lighting system was used instead of a laser lighting system, thus avoiding a lot of risk of manipulation. A real economic alternative for the experimental study of fluid-grain coupling during destabilization by localized fluidization of a granular material was created.
e results of this study of destabilization by fluidization of an immersed grain layer by a variable size basal injection showed the observation of experimental phenomenology presented in literature (the fluidized bed passed successively by the static regime, the fluidized cavity regime, and the fluidized chimney regime). e committed effort must continue; a low-cost PIV system and hydrogel beads as a solid phase, reducing health and safety requirements at work, may well replace an expensive lazer lighting system and a mixture of liquids. e use of such type of inexpensive PIV systems in the study of the fluidization destabilization of a granular material may be tainted by certain difficulties such as the presence at the market of hydrogel beads which have densities similar to that of water when they are swollen (possibility of floating of balls in the water). is is avoidable by communicating to the seller all the characteristics of the desired hydrogel beads. e hydrogel beads should also be washed or changed after a series of tests for good visibility Journal of Spectroscopy through the porous medium. We must also take advantage of the information provided to better understand qualitatively, but especially quantitatively, the mechanisms leading to the destabilization of grains. In perspective, more detailed analysis should be carried out on the mechanisms leading to the destabilization of grains and the mode of expansion of the instabilities according to the hydrodynamic regime.

Data Availability
e data used to support the findings of the article are available on reasonable request to the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest.