Hybrid Finite-Volume-Particle Method for Dusty Gas Flows
The SMAI journal of computational mathematics, Volume 3 (2017) , pp. 139-180.

We first study the one-dimensional dusty gas flow modeled by the two-phase system composed of a gaseous carrier (gas phase) and a particulate suspended phase (dust phase). The gas phase is modeled by the compressible Euler equations of gas dynamics and the dust phase is modeled by the pressureless gas dynamics equations. These two sets of conservation laws are coupled through source terms that model momentum and heat transfers between the phases. When an Eulerian method is adopted for this model, one can notice the obtained numerical results are typically significantly affected by numerical diffusion. This phenomenon occurs since the pressureless gas equations are nonstrictly hyperbolic and have degenerate structure in which singular delta shocks are formed, and these strong singularities are vulnerable to the numerical diffusion.

We introduce a low dissipative hybrid finite-volume-particle method in which the compressible Euler equations for the gas phase are solved by a central-upwind scheme, while the pressureless gas dynamics equations for the dust phase are solved by a sticky particle method. The obtained numerical results demonstrate that our hybrid method provides a sharp resolution even when a relatively small number of particle is used.

We then extend the hybrid finite-volume-particle method to the three-dimensional dusty gas flows with axial symmetry. In the studied model, gravitational effects are taken into account. This brings an additional level of complexity to the development of the finite-volume-particle method since a delicate balance between the flux and gravitational source terms should be respected at the discrete level. We test the proposed method on a number of numerical examples including the one that models volcanic eruptions.

Published online:
DOI: https://doi.org/10.5802/smai-jcm.23
Classification: 65M08,  76M12,  76M28,  86-08,  76M25,  35L65
Keywords: Two-phase dusty gas flow model, three-dimensional axial symmetry, compressible Euler equations, pressureless gas dynamics, finite-volume-particle methods, central-upwind schemes, sticky particle methods, operator splitting methods
@article{SMAI-JCM_2017__3__139_0,
author = {Alina Chertock and Shumo Cui and Alexander Kurganov},
title = {Hybrid {Finite-Volume-Particle} {Method} for {Dusty} {Gas} {Flows}},
journal = {The SMAI journal of computational mathematics},
pages = {139--180},
publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
volume = {3},
year = {2017},
doi = {10.5802/smai-jcm.23},
mrnumber = {3707729},
zbl = {1416.76143},
language = {en},
url = {https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.23/}
}
Alina Chertock; Shumo Cui; Alexander Kurganov. Hybrid Finite-Volume-Particle Method for Dusty Gas Flows. The SMAI journal of computational mathematics, Volume 3 (2017) , pp. 139-180. doi : 10.5802/smai-jcm.23. https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.23/

[1] N. Botta; R. Klein; S. Langenberg; S. Lützenkirchen Well balanced finite volume methods for nearly hydrostatic flows, J. Comput. Phys., Volume 196 (2004) no. 2, pp. 539-565 | Article | MR 2054350 | Zbl 1109.86304

[2] F. Bouchut On zero pressure gas dynamics, Advances in kinetic theory and computing (Ser. Adv. Math. Appl. Sci.) Volume 22, World Sci. Publ., River Edge, NJ, 1994, pp. 171-190 | MR 1323183 | Zbl 0863.76068

[3] F. Bouchut; F. James Duality solutions for pressureless gases, monotone scalar conservation laws, and uniqueness, Comm. Partial Differential Equations, Volume 24 (1999) no. 11-12, pp. 2173-2189 | Article | MR 1720754 | Zbl 0937.35098

[4] F. Bouchut; S. Jin; X. Li Numerical approximations of pressureless and isothermal gas dynamics, SIAM J. Numer. Anal., Volume 41 (2003), pp. 135-158 | Article | MR 1974495 | Zbl 1060.76080

[5] Y. Brenier; E. Grenier Sticky particles and scalar conservation laws, SIAM J. Numer. Anal., Volume 35 (1998), pp. 2317-2328 | Article | MR 1655848 | Zbl 0924.35080

[6] S. Carcano; L. Bonaventura; T. Esposti Ongaro; A. Neri A semi-implicit, second order accurate numerical model for multiphase underexpanded volcanic jets, Geosci. Model Dev. Discuss., Volume 6 (2013) no. 1, pp. 399-452 | Article

[7] P. Chandrashekar; C. Klingenberg A second order well-balanced finite volume scheme for Euler equations with gravity, SIAM J. Sci. Comput., Volume 37 (2015) no. 3, p. B382-B402 | Article | MR 3348134 | Zbl 1320.76078

[8] G.-Q. Chen; H. Liu Formation of $\delta$-shocks and vacuum states in the vanishing pressure limit of solutions to the Euler equations for isentropic fluids, SIAM J. Math. Anal., Volume 34 (2003), pp. 925-938 | Article | MR 1969608 | Zbl 1038.35035

[9] A. Chertock; S. Cui; A. Kurganov; Ş. N. Özcan; E. Tadmor Well-balanced central-upwind schemes for the Euler equations with gravitation (Submitted)

[10] A. Chertock; S. Cui; A. Kurganov; T. Wu Well-balanced positivity preserving central-upwind scheme for the shallow water system with friction terms, Internat. J. Numer. Meth. Fluids, Volume 78 (2015), pp. 355-383 | Article | MR 3350178

[11] A. Chertock; A. Kurganov On a hybrid finite-volume particle method, M2AN Math. Model. Numer. Anal, Volume 38 (2004), pp. 1071-1091 | Article | Numdam | MR 2108945 | Zbl 1077.65091

[12] A. Chertock; A. Kurganov On a practical implementation of particle methods, Appl. Numer. Math., Volume 56 (2006), pp. 1418-1431 | Article | MR 2245465 | Zbl 1103.65103

[13] A. Chertock; A. Kurganov; G. Petrova Finite-volume-particle methods for models of transport of pollutant in shallow water, J. Sci. Comput., Volume 27 (2006), pp. 189-199 | Article | MR 2285775 | Zbl 1101.76036

[14] A. Chertock; A. Kurganov; Yu. Rykov A new sticky particle method for pressureless gas dynamics, SIAM J. Numer. Anal., Volume 45 (2007), pp. 2408-2441 | Article | MR 2361896 | Zbl 1234.35145

[15] G.-H. Cottet; P. D. Koumoutsakos Vortex methods, Cambridge University Press, Cambridge, 2000

[16] S. Dartevelle; W. Rose; J. Stix; K. Kelfoun; J.W. Vallance Numerical modeling of geophysical granular flows: 2. Computer simulations of plinian clouds and pyroclastic flows and surges, Geochem. Geophys. Geosyst., Volume 5 (2004) no. 8 | Article

[17] V. Desveaux; M. Zenk; C. Berthon; C. Klingenberg A well-balanced scheme for the Euler equation with a gravitational potential, Finite volumes for complex applications. VII. Methods and theoretical aspects (Springer Proc. Math. Stat.) Volume 77, Springer, Cham, 2014, pp. 217-226 | MR 3213352 | Zbl 1304.76027

[18] F. Dobran; A. Neri; G. Macedonio Numerical simulation of collapsing volcanic columns, J. Geophys. Res., Volume 98 (1993), pp. 4231-4259 | Article

[19] J. Dufek; G. W. Bergantz Dynamics and deposits generated by the Kos Plateau Tuff eruption: Controls of basal particle loss on pyroclastic flow transport, Geochem. Geophys. Geosyst., Volume 8 (2007) no. 12 | Article

[20] W. E; Yu. G. Rykov; Ya. G. Sinai Generalized variational principles, global weak solutions and behavior with random initial data for systems of conservation laws arising in adhesion particle dynamics, Comm. Math. Phys., Volume 177 (1996), pp. 349-380 | Article | MR 1384139 | Zbl 0852.35097

[21] B. Einfeld On Godunov-type methods for gas dynamics, SIAM J. Numer. Anal., Volume 25 (1988), pp. 294-318 | Article | MR 933726

[22] P. Glaister Flux difference splitting for the Euler equations with axial symmetry, J. Engrg. Math., Volume 22 (1988) no. 2, pp. 107-121 | Article | MR 953171 | Zbl 0663.76070

[23] S. Gottlieb; D. Ketcheson; C.-W. Shu Strong stability preserving Runge-Kutta and multistep time discretizations, World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2011, xii+176 pages | Zbl 1241.65064

[24] S. Gottlieb; C.-W. Shu; E. Tadmor Strong stability-preserving high-order time discretization methods, SIAM Rev., Volume 43 (2001), pp. 89-112 | Article | MR 1854647 | Zbl 0967.65098

[25] M. Gurris; D. Kuzmin; S. Turek Finite element simulation of compressible particle-laden gas flows, J. Comput. Appl. Math., Volume 233 (2010) no. 12, pp. 3121-3129 | Article | MR 2595283 | Zbl 1252.76041

[26] S. Hank; R. Saurel; O. Le Metayer A hyperbolic Eulerian model for dilute two-phase suspensions, Journal of Modern Physics, Volume 2 (2011), pp. 997-1011 | Article

[27] A. Harlow; A. A. Amsden Numerical calculation of multiphase fluid flow, J. Comput. Phys., Volume 17 (1975), pp. 19-52 | Article | Zbl 0297.76079

[28] A. Harten; P. Lax; B. van Leer On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM Rev., Volume 25 (1983), pp. 35-61 | Article | MR 693713 | Zbl 0565.65051

[29] A. Kurganov; C.-T. Lin On the reduction of numerical dissipation in central-upwind schemes, Commun. Comput. Phys., Volume 2 (2007), pp. 141-163 | MR 2305919 | Zbl 1164.65455

[30] A. Kurganov; S. Noelle; G. Petrova Semi-discrete central-upwind scheme for hyperbolic conservation laws and Hamilton-Jacobi equations, SIAM J. Sci. Comput., Volume 23 (2001), pp. 707-740 | Article

[31] A. Kurganov; G. Petrova A second-order well-balanced positivity preserving central-upwind scheme for the Saint-Venant system, Commun. Math. Sci., Volume 5 (2007), pp. 133-160 | Article | MR 2310637 | Zbl 1226.76008

[32] A. Kurganov; E. Tadmor New high resolution central schemes for nonlinear conservation laws and convection-diffusion equations, J. Comput. Phys., Volume 160 (2000), pp. 241-282 | Article | MR 1756766 | Zbl 0987.65085

[33] A. Kurganov; E. Tadmor Solution of two-dimensional Riemann problems for gas dynamics without Riemann problem solvers, Numer. Methods Partial Differential Equations, Volume 18 (2002), pp. 584-608 | Article | MR 1919599 | Zbl 1058.76046

[34] M.-C. Lai; C. S. Peskin An immersed boundary method with formal second-order accuracy and reduced numerical viscosity, J. Comput. Phys., Volume 160 (2000), pp. 705-719 | Article | MR 1763828 | Zbl 0954.76066

[35] R. J. LeVeque The dynamics of pressureless dust clouds and delta waves, J. Hyperbolic Differ. Equ., Volume 1 (2004), pp. 315-327 | Article | MR 2070130 | Zbl 1079.76074

[36] R. J. LeVeque; D. S. Bale Wave propagation methods for conservation laws with source terms, Hyperbolic problems: theory, numerics, applications, Vol. II (Zürich, 1998) (Internat. Ser. Numer. Math.) Volume 130, Birkhäuser, Basel, 1999, pp. 609-618 | MR 1717233 | Zbl 0927.35062

[37] G. Li; Y. Xing High order finite volume WENO schemes for the Euler equations under gravitational fields, J. Comput. Phys., Volume 316 (2016), pp. 145-163 | Article | MR 3494348 | Zbl 1349.76356

[38] K.-A. Lie; S. Noelle On the artificial compression method for second-order nonoscillatory central difference schemes for systems of conservation laws, SIAM J. Sci. Comput., Volume 24 (2003) no. 4, pp. 1157-1174 | Article | MR 1976211 | Zbl 1038.65078

[39] J. Luo; K. Xu; N. Liu A well-balanced symplecticity-preserving gas-kinetic scheme for hydrodynamic equations under gravitational field, SIAM J. Sci. Comput., Volume 33 (2011) no. 5, pp. 2356-2381 | Article | MR 2837536 | Zbl 1232.76044

[40] H. Miura; I. I. Glass On a dusty-gas shock tube, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci., Volume 382 (1982) no. 1783, pp. 373-388 | Article

[41] A. Neri; T. E. Ongaro; G. Macedonio; D Gidaspow Multiparticle simulation of collapsing volcanic columns and pyroclastic flow, J. Geophys. Res, Volume 108 (2003), pp. 1-22 | Article

[42] H. Nessyahu; E. Tadmor Nonoscillatory central differencing for hyperbolic conservation laws, J. Comput. Phys., Volume 87 (1990) no. 2, pp. 408-463 | Article | MR 1047564 | Zbl 0697.65068

[43] B. Nilsson; O. S. Rozanova; V. M. Shelkovich Mass, momentum and energy conservation laws in zero-pressure gas dynamics and $\delta$-shocks: II, Appl. Anal., Volume 90 (2011) no. 5, pp. 831-842 | Article | MR 2786660 | Zbl 1222.35122

[44] B. Nilsson; V. M. Shelkovich Mass, momentum and energy conservation laws in zero-pressure gas dynamics and delta-shocks, Appl. Anal., Volume 90 (2011) no. 11, pp. 1677-1689 | Article | MR 2842606 | Zbl 1237.35111

[45] M. Pelanti; R. J. Leveque High-resolution finite volume methods for dusty gas jets and plumes, SIAM J. Sci. Comput., Volume 28 (2006), pp. 1335-1360 | Article | MR 2255460 | Zbl 1123.76042

[46] C. S. Peskin The immersed boundary method, Acta Numer., Volume 11 (2002), pp. 479-517 | Article | MR 2009378 | Zbl 1123.74309

[47] P.-A. Raviart An analysis of particle methods, Numerical methods in fluid dynamics (Como, 1983) (Lecture Notes in Math.) Volume 1127, Springer, Berlin, 1985, pp. 243-324 | MR 802214

[48] Yu. G. Rykov Propagation of singularities of shock wave type in a system of equations of two-dimensional pressureless gas dynamics, Mat. Zametki, Volume 66 (1999), p. 760-769 (Russian); translation in Math. Notes 66 (1999), pp. 628–635 (2000) | MR 1754475 | Zbl 1003.76044

[49] Yu. G. Rykov On the nonhamiltonian character of shocks in 2-D pressureless gas, Boll. Unione Mat. Ital. Sez. B Artic. Ric. Mat., Volume (8) 5 (2002), pp. 55-78 | MR 1881444 | Zbl 1096.35117

[50] T. Saito Numerical analysis of dusty-gas flows, J. Comput. Phys., Volume 176 (2002), pp. 129-144 | Article | Zbl 1120.76350

[51] B. Shotorban; G. B. Jacobs; O. Ortiz; Q. Truong An Eulerian model for particles nonisothermally carried by a compressible fluid, Int. J. Heat Mass Transfer, Volume 65 (2013), pp. 845-854 | Article

[52] C.-W. Shu; S. Osher Efficient implementation of essentially non-oscillatory shock-capturing schemes, J. Comput. Phys., Volume 77 (1988), pp. 439-471 | Article

[53] G. Strang On the construction and comparison of difference schemes, SIAM J. Numer. Anal., Volume 5 (1968), pp. 506-517 | Article | MR 235754 | Zbl 0184.38503

[54] P. K. Sweby High resolution schemes using flux limiters for hyperbolic conservation laws, SIAM J. Numer. Anal., Volume 21 (1984) no. 5, pp. 995-1011 | Article | MR 760628 | Zbl 0565.65048

[55] C. T. Tian; K. Xu; K. L. Chan; L. C. Deng A three-dimensional multidimensional gas-kinetic scheme for the Navier-Stokes equations under gravitational fields, J. Comput. Phys., Volume 226 (2007) no. 2, pp. 2003-2027 | Article | MR 2356404 | Zbl 1388.76313

[56] R. Touma; U. Koley; C. Klingenberg Well-balanced unstaggered central schemes for the Euler equations with gravitation, SIAM J. Sci. Comput., Volume 38 (2016) no. 5, p. B773-B807 | Article | MR 3546981 | Zbl 1346.76100

[57] G. Valentine; K. Wohletz Numerical models of Plinian eruption columns and pyroclastic, J. Geophys. Res, Volume 94 (1989), pp. 1867-1887 | Article

[58] B. van Leer Towards the ultimate conservative difference scheme. V. A second-order sequel to Godunov’s method, J. Comput. Phys., Volume 32 (1979) no. 1, pp. 101-136 | Article | Zbl 1364.65223

[59] K. H. Wohletz; T. R. McGetchin; M. T. Sandford II; E. M. Jones Hydrodynamic aspects of caldera-forming eruptions: numerical models, J. Geophys. Res, Volume 89 (1984), pp. 8269-8285 | Article

[60] Y. Xing; C.-W. Shu High order well-balanced WENO scheme for the gas dynamics equations under gravitational fields, J. Sci. Comput., Volume 54 (2013) no. 2-3, pp. 645-662 | Article | MR 3011375 | Zbl 1260.76022

[61] K. Xu; J. Luo; S. Chen A well-balanced kinetic scheme for gas dynamic equations under gravitational field, Adv. Appl. Math. Mech., Volume 2 (2010), pp. 200-210 | MR 2606447

[62] M. Yuen Some exact blowup solutions to the pressureless Euler equations in ${ℝ}^{N}$, Commun. Nonlinear Sci. Numer. Simul., Volume 16 (2011) no. 8, pp. 2993-2998 | Article | MR 2782989 | Zbl 1419.76542