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.

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] 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] 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] 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] 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] Sticky particles and scalar conservation laws, SIAM J. Numer. Anal., Volume 35 (1998), pp. 2317-2328 | Article | MR 1655848 | Zbl 0924.35080

[6] 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] 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] 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] Well-balanced central-upwind schemes for the Euler equations with gravitation (Submitted)

[10] 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] 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] On a practical implementation of particle methods, Appl. Numer. Math., Volume 56 (2006), pp. 1418-1431 | Article | MR 2245465 | Zbl 1103.65103

[13] 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 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] Vortex methods, Cambridge University Press, Cambridge, 2000

[16] 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] 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] Numerical simulation of collapsing volcanic columns, J. Geophys. Res., Volume 98 (1993), pp. 4231-4259 | Article

[19] 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] 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] On Godunov-type methods for gas dynamics, SIAM J. Numer. Anal., Volume 25 (1988), pp. 294-318 | Article | MR 933726

[22] 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] 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] Strong stability-preserving high-order time discretization methods, SIAM Rev., Volume 43 (2001), pp. 89-112 | Article | MR 1854647 | Zbl 0967.65098

[25] 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] A hyperbolic Eulerian model for dilute two-phase suspensions, Journal of Modern Physics, Volume 2 (2011), pp. 997-1011 | Article

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

[28] 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] 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] 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 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] Multiparticle simulation of collapsing volcanic columns and pyroclastic flow, J. Geophys. Res, Volume 108 (2003), pp. 1-22 | Article

[42] 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] 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] 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] 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] The immersed boundary method, Acta Numer., Volume 11 (2002), pp. 479-517 | Article | MR 2009378 | Zbl 1123.74309

[47] 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] 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] 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] Numerical analysis of dusty-gas flows, J. Comput. Phys., Volume 176 (2002), pp. 129-144 | Article | Zbl 1120.76350

[51] An Eulerian model for particles nonisothermally carried by a compressible fluid, Int. J. Heat Mass Transfer, Volume 65 (2013), pp. 845-854 | Article

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

[53] 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] 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] 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] 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] Numerical models of Plinian eruption columns and pyroclastic, J. Geophys. Res, Volume 94 (1989), pp. 1867-1887 | Article

[58] 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] Hydrodynamic aspects of caldera-forming eruptions: numerical models, J. Geophys. Res, Volume 89 (1984), pp. 8269-8285 | Article

[60] 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] 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] Some exact blowup solutions to the pressureless Euler equations in ${\mathbb{R}}^{N}$, Commun. Nonlinear Sci. Numer. Simul., Volume 16 (2011) no. 8, pp. 2993-2998 | Article | MR 2782989 | Zbl 1419.76542