2024. Vol. 19. Issue 3, Pp. 119–124

URL: http://mfs.uimech.org/mfs2024.3.017,en

DOI: 10.21662/mfs2024.3.017

URL: http://mfs.uimech.org/mfs2024.3.017,en

DOI: 10.21662/mfs2024.3.017

Numerical study of the interaction of a shock wave with a layer of particles using the Baer–Nunziato model

P.A. Chuprov

ICAD RAS

The work is devoted to the application of the Baer–Nunziato model to study shock wave processes in particle beds. The numerical algorithm is based on an HLLC-like solver and takes into account the processes of establishing equilibrium at the interface, taking into account the effect of particle compaction. The problem statement corresponds to the full-scale experiment of B.C. Fan et al. As part of the work, qualitative and quantitative agreement was obtained with the results of natural experiments and calculations of other authors, and an explanation was proposed for the experimentally observed phenomena.

numerical modeling,

Baer–Nunziato equations,

compaction

Wave processes in granular media occur in various situations, e.g. during the gunpowder combustion or catastrophic events in coal mines and factories, known as dust explosions. Dust explosions are dangerous the fatal events, occurring in the presence of the chemically active dust. Shock waves are generated during such events, enhancing the explosion, even transitioning deflagration to detonation in some cases. Numerical modeling is an extremely valuable tool for understanding of the nature of such processes. The goal of this paper is to present a numerical study of a shock wave interaction with a layer of inert dust.

In this paper the numerical algorithm for modelling high energy wave processes in two-phase media is presented. This algorithm is based on the 2D Baer-Nunziato model and HLLC-like solver. Numerical algorithm includes pressure relaxation substep, which takes into account intergranular stress. This phenomenon is extremely important for qualitatively and quantitatively correct description of the two-phase wave processes in granular media. Interphase momentum transfer is also taken into account. HLLC-like solver is developed for the 2D BN equations for the particular case of phases being near mechanical equilibrium. HLLC solver allows for the better volume fraction discontinuity resolution that HLL method, but simple enough to be extrapolated for the three-phase case in the future work. Overall numerical algorithm is based on the Strang splitting technique with three stages: hyperbolic, pressure relaxation and momentum exchange. For the momentum exchange substep the explicit Euler method is used.

Described numerical algorithm was verified using data for the Rogue shock tube test. The algorithm achieved great qualitative and quantitative agreement with experimental data.
After verification, the algorithm was used to simulate the experiment of B.C. Fan on interaction of the passing shock wave with the dust layer. Within the framework of this paper
only wave effects inside the dust layer are considered. The effect of dust dispersion are neglected. It is shown in experiment that the layer is deformed by the shock wave,
deflecting from the initial level by some angle _{1}. Moreover, under this deformed surface is present the compacted layer of the dust.
This layer propagates at an angle _{2}. It is shown that _{1} increases with increase of Mach number of the passing wave, but _{2} decreases. Numerical experiments have shown that the compacted layer of dust acts as a plug, preventing penetration of the pressure wave from the region behind the shock wave inside the dust layer. Obviously with increase of the Mach number, layer is more and more compressed, increasing the effectiveness of the compaction plug. The parametric study was conducted, which confirmed the leading role of the plug in the formation of the flow inside the dust layer, and shed some light onto the dust compaction mechanism itself.

In the current paper, the robust and simple numerical algorithm for multiphase flow modeling is presented. This algorithm is verified and used to analyze inner structure of the compressed dust layer. The obtained results show agreement with the experimental data and are used to suggest an explanation to the experimentally known phenomenon.

- Shimura K., Matsuo A. Using an extended CFD-DEM for the two-dimensional simulation of shock-induced layered coal-dust combustion in a narrow
channel // Proc. Comb. Inst. 2019. V. 37, Issue 3. Pp. 3677–3684.

DOI: 10.1016/j.proci.2018.07.066 - Song Y., Zhang Q. Criterion and propagation process of spark-induced dust layered explosion // Fuel. 2020. V. 267. Pp. 117205.

DOI: 10.1016/j.fuel.2020.117205 - Guhathakurta S., Houim R.W. Impact of particle diameter and thermal radiation on the explosion of dust layers // Proc. Comb. Inst. 2023. V. 39, Issue 3.
Pp. 2905–2914.

DOI: 10.1016/j.proci.2022.10.011 - Baer M.R., Nunziato J.W. A two-phase mixture theory for the deflagration-to-detonation transition in reactive granular materials // International
Journal of Multiphase Flow. 1986. V. 12, Issue 6. Pp. 861–889.

DOI: 10.1016/0301-9322(86)90033-9 - Chuprov P.A. Numerical simulation of combustion wave propagation in a black powder charge using a two-fluid model // Adv. Theory and Practice
of Comp. Mech. Smart Innovation, Sys. Tech. Springer. 2022. V. 274. Pp. 167–178.

DOI: 10.1007/978-981-16-8926-0_12 - Saurel R., Favrie N., Petitpas F., Lallemand M.-H., Gavrilyuk S.K. Modeling dynamic and irreversible powder compaction // J. Fluid Mech. 2010. V. 664.
Pp. 348–396.

DOI: 10.1017/S0022112010003794 - Houim R.W., Oran E.S. A multiphase model for compressible granular-gaseous flows: formulation and initial tests // J. Fluid Mech. 2016. V. 789.
Pp. 166–220.

DOI: 10.1017/jfm.2015.728 - Li Q., Feng J., Cai T., Hu C. Difference scheme for two-phase flow // App. Math. Mech. 2004. V. 25. Pp. 536–545.

DOI: 10.1007/BF02437602 - Liang S., Liu W., Yuan L. Solving seven-equation model for compressible two-phase glow using multiple GPUs // Comp. & Fluids. 2014. V. 99. Pp. 156–
171.

DOI: 10.1016/j.compfluid.2014.04.021 - Tokareva S.A., Toro E.F. HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow // J. Comp. Phys. 2010. V. 229,
Issue 10. Pp. 3573–3604.

DOI: 10.1016/j.jcp.2010.01.016 - Lochon H., Daude F., Galon , HerardJ.-M. HLLC-type Riemann solver with approximated two-phase contact for the computation of the Baer-Nunziato
two-fluid model // J. Comp. Phys. 2016. V. 326. Pp. 733–762.

DOI: 10.1016/j.jcp.2016.09.015 - Fan B.C., Chen Z.H., JiangX.H., Li H.Z. Interaction of a shock wave with a loose dusty bulk layer // Shock Waves. 2007. V. 16. Pp. 179–187.

DOI: 10.1007/s00193-006-0059-5