2019. Vol. 14. Issue 2, Pp. 108–114

URL: http://mfs.uimech.org/mfs2019.2.015

DOI: 10.21662/mfs2019.2.015

URL: http://mfs.uimech.org/mfs2019.2.015

DOI: 10.21662/mfs2019.2.015

Modeling of the spherical explosion attenuation process using aqueous foam

Bolotnova R.Kh.^{∗}, Gainullina E.F.^{∗}, Nurislamova E.A.^{∗∗}

The two-phase model of dry aqueous foam dynamic behavior under the strong shock wave influence is presented under assumption that the foam structure under shock loading is destroyed into a suspension of monodispersed microdrops with the formation of a gas-droplet mixture. The system of equations for the model of aqueous foam includes the laws of conservation of mass, momentum and energy for each phase in accordance with the single-pressure, two-speed, two-temperature approximations in a three-dimensional formulation, taking into account the Schiller–Naumann interfacial drag force and the Ranz–Marshall interfacial contact heat transfer. The thermodynamic properties of air and water forming a gas-droplet mixture are described by the Peng–Robinson and Mie–Grueneisen equations of state. The presence of non-uniform process in height of aqueous foam syneresis, which is due to gravitational forces, is taken into account by setting the distribution of the liquid volume fraction in the foam. An additional consideration of the syneresis process during calculating the intensity of interphase drag forces according to the Schiller–Naumann model was controlled by introducing the parameter depending on the spatial distribution of the initial liquid volume fraction of the foam. The spherical explosion is modeled in the form of the shock wave pulse whose energy coincided with the charge energy of the HE used in the experiments. The problem numerical solution is implemented using the OpenFOAM free software package based on the two-step PIMPLE computational algorithm. The numerical solution of the problem, obtained on the basis of the proposed gas-droplet mixture model, is in satisfactory agreement with the experimental data on a spherical explosion in aqueous foam. The analysis of the spherical shock wave dynamics while its propagation through aqueous foam is given. The causes of the significant decrease in the amplitude and velocity shock waves propagation in the medium under study are investigated.

spherical shock wave,

aqueous foam,

OpenFOAM package,

numerical modeling

*Purpose.* Investigation of dynamics a spherical shock wave propagation in dry aqueous foam, taking into account the influence of interfacial drag forces, contact heat transfer and foam syneresis process.

*Methodology.* The two-phase model of gas-liquid mixture in the two-dimensional axisymmetric formulation under the conditions of single-pressure, two-speed, two-temperature approximations is proposed. While describing the behavior of aqueous foam under shock loading, it is assumed that in the process of its interaction with the strong shock wave, the foam structure is being destroyed into monodisperse microdroplets, and model of the gas-droplet mixture is applicable. The system of model equations includes the laws of conservation of mass, momentum and energy for each phase, taking into account the Schiller-Naumann interfacial drag forces and the Ranz-Marshall interphase heat transfer processes. The non-uniformity of foam deposition process is taken into account in the Schiller-Naumann drag model by introducing the parameter that depends on the spatial distribution of the initial liquid volume fraction of the foam. The Peng-Robinson and Mie-Grueneisen equations of state are used to describe the properties of gas and water. The proposed gas-droplet mixture model is implemented in the solver created in OpenFOAM numerical simulation open package. Discretization of the model equations system is carried out by the finite volume method. The two-step computational algorithm PIMPLE is chosen as the solution method.

*Findings.* The numerical simulation of a spherical shock wave evolution is performed under conditions of experiments on the HE explosion in the dry aqueous foam with the initial liquid volume fraction of 0.0083. The spherical region with the radius of 2 m is considered. The explosion is simulated in the center of region at initial time in the form of pressure pulse approximated by the smooth exponential function. The generated shock wave pulse with the maximum amplitude of 3000 MPa has the energy equal to the explosion energy in the experiment under consideration. The numerical modeling results of the task are compared with experimental data obtained from sensors attached at different heights and distances from the explosion center. To account for the foam syneresis process for each of the gauges, a separate problem was solved with the individual initial distribution of liquid fraction in the foam. Both in calculations and in experiments, the shock wave pulse at the short distance from the center of the explosion has two-wave structure, which, as shock wave propagates deep into the foam, is “spread” under the influence of interfacial drag forces and contact heat transfer. The comparative analysis of the calculated and experimental pressure profiles at the locations of the gauges demonstrates their good agreement.

*Value.* The numerical simulation and study of process of spherical shock wave propagation formed as a result of the HE explosion in aqueous foam was performed for the conditions of the selected experiments. While studying the shock wave dynamics in aqueous foam, the proposed model of the gas-droplet mixture that takes into account the forces of interfacial interaction, contact heat transfer and the effect of the aqueous foam syneresis due to gravitational forces was used. The numerical simulation of the investigated problem is based on the free OpenFOAM package. The comparison of the numerical solutions obtained for the pressure time dependences and experimental pressure oscillograms at the points of the sensors installation showed their good agreement on the velocity of propagation and the amplitude of shock wave pulse. The analysis the causes of shock wave energy dissipation in aqueous foam, leading to its significant attenuation was given.

- Bolotnova R.Kh., Galimzianov M.N., Agisheva U.O. [Simulation of a strong shock wave interaction process in gas liquid mixtures]
[University proceedings. Volga region. Physical and mathematical sciences]. 2011. No. 2. Pp. 3–14 (in Russian).

https://izvuz_fmn.pnzgu.ru/files/izvuz_fmn.pnzgu.ru/1211.pdf - Bolotnova R.Kh., Galimzianov M.N., Topolnikov A.S., Buzina V.A., Agisheva U.O. Nonlinear Effects In Bubbly Liquid With Shock Waves // World Academy of Science, Engineering and Technology. 2012. Vol. 6, No. 8. Pp. 1095–1102.

https://waset.org/publications/14192/nonlinear-effects-in-bubbly-liquid-with-shock-waves - Bolotnova R.Kh., Gainullina E.F. [The numerical modeling of spherical explosion in the foam] Trudy Instituta Mexaniki im. R.R. Mavlyutova Ufimskogo Nauchnogo Centra Rossijskoj Akademii Nauk
[Proceedings of the Mavlyutov Institute of Mechanics]. 2016. V. 11, No. 1. Pp. 60–65 (in Russian).

DOI: 10.21662/uim2016.1.009 - Bolotnova R.Kh., Gainullina E.F. [A research of damping properties of aqueous foam under the impact of spherical shock waves] Izvestiya VUZov. Povolzhskij Region.
Fiziko-Matematicheskie Nauki [University proceedings. Volga region. Physical and mathematical sciences]. 2017. No. 2. Pp. 108–121 (in Russian).

DOI: 10.21685/2072-3040-2017-2-9 - Bolotnova R.Kh., Gainullina E.F. Influence of Heat-exchange Processes on Decreasing an Intensity of a Spherical Explosion in Aqueous Foam // Fluid Dynamics. 2019. V. 54, Suppl. 7. Pp. S71–S78.

DOI: 10.1134/S0032823519030020 - Del Prete E., Chinnayya A., Domergue L., et al. Blast Wave Mitigation by Dry Aqueous Foams // Shock Waves. 2013. Vol. 23, No. 1 Pp. 39–53.

DOI: 10.1007/s00193-012-0400-0 - Bolotnova R.Kh., Agisheva U.O. [Spatial modeling of water foam dynamics with moving Lagrangian grids under shock wave impact] Vychislitel’nye metody i programmirovanie: novye
vychislitel’nye texnologii [Numerical methods and Programming]. 2014. No. 3. Pp. 427–440 (in Russian).

http://num-meth.srcc.msu.ru/zhurnal/tom_2014/pdf/v15r138.pdf - Agisheva U.O., Bolotnova R.Kh., Gainullina E.F., Korobchinskaya V.A. Features of Vortex Formation under the Impact of a Pressure Pulse on a Gas Region Bounded by the Foam Layer //
Fluid Dynamics. 2016. V. 51, No. 6. P. 757–766.

DOI: 10.1134/S0015462816060053 - Bolotnova R.Kh., Gainullina E.F. [Investigation of axisymmetric wave flows under interaction of a spherical impact pulse with a barrier of aqueous foam]
Trudy Instituta Mexaniki im. R.R. Mavlyutova Ufimskogo Nauchnogo Centra Rossijskoj Akademii Nauk [Proceedings of the Mavlyutov Institute of Mechanics]. 2017. V. 12, No. 2. Pp. 238—243 (in Russian).Болотнова Р.Х., Гайнуллина Э.Ф. Исследование осесимметричных волновых течений при взаимодействии сферического ударного импульса с барьером из водной пены //

DOI: 10.21662/uim2017.2.036 - Bolotnova R.Kh., Gainullina E.F. Features of the impact of spherical shock impulse on the boundary of gas with aqueous foam Vestnik Udmurtskogo Universiteta. Matematika. Mexanika.
Komp’yuternye Nauki [The Bulletin of Udmurt University. Mathematics. Mechanics. Computer Science]. 2018. V. 28, No. 3. Pp. 364–372 (in Russian).

DOI: 10.20537/vm180307 - Bolotnova R.Kh., Gainullina E.F. Wave dynamics and vortex formation under the impact of a spherical impulse on the boundary between gas and aqueous foam // J. Phys.: Conf. Ser. 2019. Vol. 1268. 012015.

DOI: 10.1088/1742-6596/1268/1/012015 - OpenFOAM. The Open Source Computational Fluid Dynamics (CFD) Toolbox.

http://www.openfoam.com - Zeno Tacconi. Feasibility analysis of a two-fluid solver for cavitation and interface capturing as implemented in OpenFOAM // Tesi di Laurea Magistrale in Ingegneria Energetica, Politecnico di Milano. 2018. 134 p.

https://www.politesi.polimi.it/handle/10589/139684 - Ranz W.E., Marshall W.R. Evaporation from Drops // Chem. Eng. Prog. 1952. Vol. 48, No. 22. P. 141–146.

http://dns2.asia.edu.tw/~ysho/YSHO-English/1000%20CE/PDF/Che%20Eng%20Pro48,%20141.pdf - Schiller L., Naumann Z. A Drag Coefficient Correlation // Z. Ver. Deutsch. Ing. 1935. Vol. 77. P. 40—65.
- Peng D.Y., Robinson D.B. A new two-constant equation of state // Industrial and Engineering Chemistry: Fundamentals. 1976. Vol. 15. P. 59–64.

DOI: 10.1021/i160057a011 - Jolgam S., Ballil A., Nowakowski A., Nicolleau F. On Equations of State for Simulations of Multiphase Flows // Proceedings of the World congress on Engineering. 2012. Vol. III.

http://www.iaeng.org/publication/WCE2012/WCE2012_pp1963-1968.pdf - Zhdan S.A. Numerical modeling of the explosion of a high explosive (HE) charge in foam // Combust., Explos., Shock Waves. 1990. V. 26. No. 2. Pp. 221–227.

https://elibrary.ru/item.asp?id=30555610