A particle-resolved direct numerical simulation method for the compressible gas flow and arbitrary shape solid moving with a uniform framework. (English) Zbl 07901317
Summary: Compressible particle-resolved direct numerical simulations (PR-DNS) are widely used in explosion-driven dispersion of particles simulations, multiphase turbulence modelling, and stage separation for two-stage-to-orbit vehicles. The direct forcing immersed boundary method (IBM) is a promising method and widely applied in low speed flow while there is few research regarding compressible flows. We developed a novel IBM to resolve supersonic and hypersonic gas flows interacting with irregularly shaped multi-body particle. The main innovation is that current method can solve the interaction of particles and high-speed fluids, particle translation and rotation, and collision among complex-shaped particles within a uniform framework. Specially, high conservation and computation consumption are strictly satisfied, which is critical for resolving the high speed compressible flow feature. To avoid the non-physical flow penetration around particle surface, an special iterative algorithm is specially derived to handle the coupling force between the gas and particles. The magnitude of the velocity difference error could be reduced by 6-8 orders compared to that of a previous method. Additionally, aerodynamic force integration was achieved using the momentum equation to ensure momentum conservation for two-phase coupling. A high-efficiency cell-type identification method for each step was proposed, and mapping among LPs and cells was used again to select the immersed cells. As for the collision force calculation, the complex shape of a particle was represented by a cloud of LPs and the mapping of LPs and cells was used to reduce the complexity of the algorithm for contact searching. The repetitive use of the mapping relationship could reduce the internal memory and improve the efficiency of the proposed algorithm. Moreover, various verification cases were conducted to evaluate the simulation performance of the proposed algorithm, including two- and three-dimensional moving and motionless particles with regular and complex shapes interacting with high-speed flow. Specifically, an experiment involving a shock passing through a sphere was designed and conducted to provide high-precision data. The corresponding results of the large-scale numerical simulation agree well with those obtained experimentally. The current method supports flow simulations at a particle-resolved scale in engineering.
Keywords:
compressible multiphase flow; particle-resolved direct numerical simulation; immersed boundary method; flux conservationSoftware:
CUDAReferences:
| [1] | Sundaresan, S.; Ozel, A.; Kolehmainen, J., Toward constitutive models for momentum, species, and energy transport in gas-particle flows, Annu. Rev. Chem. Biomol. Eng., 9, 1, 2018 |
| [2] | Zhong, W.; Yu, A.; Zhou, G., ChemInform abstract: CFD simulation of dense particulate reaction system: approaches, recent advances and applications, Cheminform, 140, 50, 16-43, 2016 |
| [3] | Falah, A.; Bernd, Improvement, validation and application of CFD]DEM model to dense gas-solid flow in a fluidized bed, 2013, Cambridge University Press |
| [4] | Baer, M.; Nunziato, J., A two-phase mixture theory for the deflagration-to-detonation transition (DDT) in reactive granular materials, Int. J. Multiph. Flow, 12, 6, 861-889, 1986 · Zbl 0609.76114 |
| [5] | Saurel, R.; Abgrall, R., A simple method for compressible multifluid flows, SIAM J. Sci. Comput., 21, 3, 1115-1145, 1999 · Zbl 0957.76057 |
| [6] | Poroshyna, Y. E.; Utkin, P., Numerical simulation of a normally incident shock wave-dense particles layer interaction using the Godunov solver for the Baer-Nunziato equations, Int. J. Multiph. Flow, 142, Article 103718 pp., 2021 |
| [7] | Snider, D. M.; Clark, S. M.; O’rourke, P. J., Eulerian-Lagrangian method for three-dimensional thermal reacting flow with application to coal gasifiers, Chem. Eng. Sci., 66, 6, 1285-1295, 2011 |
| [8] | Tenneti, S.; Subramaniam, S., Particle-resolved direct numerical simulation for gas-solid flow model development, Annu Rev. Fluid. Mech., 46, 1, 277-286, 2014 |
| [9] | Osnes, A. N.; Vartdal, M.; Omang, M. G., Computational analysis of shock-induced flow through stationary particle clouds, Int. J. Multiph. Flow, 114, 268-286, 2019 |
| [10] | Xia, Y.; Lin, Z.; Pan, D., Turbulence modulation by finite-size heavy particles in a downward turbulent channel flow, Phys. Fluids, 33, 6, 2021 |
| [11] | Fan, P.; Lin, Z.; Xu, J., A three-dimensional fictitious domain method for direct numerical simulations of particle-laden flows with heat transfer, Phys. Fluids, 35, 6, 2023 |
| [12] | Zhu, C.; Yu, Z.; Pan, D., Interface-resolved direct numerical simulations of the interactions between spheroidal particles and upward vertical turbulent channel flows, J. Fluid. Mech., 891, A6, 2020 · Zbl 1460.76868 |
| [13] | Souli, M.; Ouahsine, A.; Lewin, L., ALE formulation for fluid-structure interaction problems, Comput. Methods Appl. Mech. Eng., 190, 5-7, 659-675, 2000 · Zbl 1012.76051 |
| [14] | Shen, Y.; Ren, Y.; Ding, H., A 3D conservative sharp interface method for simulation of compressible two-phase flows, J. Comput. Phys., 403, Article 109107 pp., 2020 · Zbl 1453.76111 |
| [15] | Lin, J.-Y.; Shen, Y.; Ding, H., Simulation of compressible two-phase flows with topology change of fluid-fluid interface by a robust cut-cell method, J. Comput. Phys., 328, 140-159, 2017 |
| [16] | Geisenhofer, M.; Kummer, F.; Müller, B., A discontinuous Galerkin immersed boundary solver for compressible flows: adaptive local time stepping for artificial viscosity-based shock-capturing on cut cells, Int. J. Numer. Methods Fluids, 91, 9, 448-472, 2019 |
| [17] | Peskin, C. S., The immersed boundary method, Acta Numer., 11, 479-517, 2002 · Zbl 1123.74309 |
| [18] | Peskin, C. S., Flow patterns around heart valves: a numerical method, J. Comput. Phys., 10, 2, 252-271, 1972 · Zbl 0244.92002 |
| [19] | Hashemi, Z.; Abouali, O.; Ahmadi, G., Direct numerical simulation of particle-fluid interactions: a review, Iran. J. Sci. Technol., Trans. Mech. Eng., 41, 71-89, 2017 |
| [20] | Huang, W.-X.; Tian, F.-B., Recent trends and progress in the immersed boundary method, Proc. Instit. Mech. Eng. Part C: J. Mech. Eng. Sci., 233, 23-24, 7617-7636, 2019 |
| [21] | De Vanna, F.; Picano, F.; Benini, E., A sharp-interface immersed boundary method for moving objects in compressible viscous flows, Comput. Fluids., 201, Article 104415 pp., 2020 · Zbl 1519.76206 |
| [22] | Abalakin, I. V.; Vasilyev, O. V.; Zhdanova, N. Y.S., Characteristic based volume penalization method for numerical simulation of compressible flows on unstructured meshes, Comput. Math. Math. Phys, 61, 1315-1329, 2021 · Zbl 1479.76065 |
| [23] | Uddin, H.; Kramer, R. M.; Pantano, C., A Cartesian-based embedded geometry technique with adaptive high-order finite differences for compressible flow around complex geometries, J. Comput. Phys., 262, 379-407, 2014 · Zbl 1349.76542 |
| [24] | Mizuno, Y.; Takahashi, S.; Nonomura, T., A simple immersed boundary method for compressible flow simulation around a stationary and moving sphere, Math. Probl. Eng., 2015, 2015 · Zbl 1394.76113 |
| [25] | Zhou, K.; Balachandar, S., An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing, J. Comput. Phys., 424, Article 109862 pp., 2021 · Zbl 07508467 |
| [26] | Carmouze, Q.; Fraysse, F.; Saurel, R., Coupling rigid bodies motion with single phase and two-phase compressible flows on unstructured meshes, J. Comput. Phys., 375, 1314-1338, 2018 |
| [27] | Kapahi, A.; Mousel, J.; Sambasivan, S., Parallel, sharp interface Eulerian approach to high-speed multi-material flows, Comput. Fluids, 83, 144-156, 2013 · Zbl 1290.76134 |
| [28] | Van Noordt, W.; Ganju, S.; Brehm, C., Immersed-boundary wall-modeled large-eddy simulation of high mach number boundary layer flows, (Proceedings of the AIAA Aviation 2021 Forum, F, 2021), C |
| [29] | Chaudhuri, A.; Hadjadj, A.; Chinnayya, A., On the use of immersed boundary methods for shock/obstacle interactions, J. Comput. Phys., 230, 5, 1731-1748, 2011 · Zbl 1391.76531 |
| [30] | Pantano, C.; Karagiozis, K.; Kamakoti, R., A low-numerical dissipation immersed interface method for the compressible Navier-Stokes equations, (Proceedings of the APS Division of Fluid Dynamics Meeting Abstracts, F, 2008), [C] |
| [31] | Takahashi, S.; Nonomura, T.; Fukuda, K., A numerical scheme based on an immersed boundary method for compressible turbulent flows with shocks: application to two-dimensional flows around cylinders, J. Appl. Math., 2014, 2014 |
| [32] | Sambasivan, S. K.; Udaykumar, H., Ghost fluid method for strong shock interactions part 2: immersed solid boundaries, AIAA J., 47, 12, 2923-2937, 2009 |
| [33] | Wang, J.; Gu, X.; Wu, J., A sharp-interface immersed boundary method for simulating high-speed compressible inviscid flows, Adv. Aerodyn., 2, 1, 1-19, 2020 |
| [34] | Borazjani, I., Sharp-interface immersed-boundary method for compressible flows with shock-particle interaction, AIAA J., 59, 4, 1169-1183, 2021 |
| [35] | Liu, J.; Qiu, J.; Hu, O., Adaptive Runge-Kutta discontinuous Galerkin method for complex geometry problems on Cartesian grid, Int. J. Numer. Methods Fluids, 73, 10, 847-868, 2013 · Zbl 1455.65151 |
| [36] | Soni, V.; Hadjadj, A.; Roussel, O., On the use of adaptive multiresolution method with time-varying tolerance for compressible fluid flows, Shock Waves, 29, 37-50, 2019 |
| [37] | Reviznikov, D. L.; Sposobin, A. V.E.; Ivanov, I. E., Oscillatory flow regimes resulting from the shock layer-particle interaction, High Temp., 58, 278-283, 2020 |
| [38] | Boustani, J.; Barad, M. F.; Kiris, C. C., An immersed boundary fluid-structure interaction method for thin, highly compliant shell structures, J. Comput. Phys., 438, Article 110369 pp., 2021 · Zbl 07505964 |
| [39] | Wang, L.; Currao, G. M.; Han, F., An immersed boundary method for fluid-structure interaction with compressible multiphase flows, J. Comput. Phys., 346, 131-151, 2017 · Zbl 1378.74038 |
| [40] | Chaudhuri, A.; Hadjadj, A.; Sadot, O., Computational study of shock-wave interaction with solid obstacles using immersed boundary methods, Int. J. Numer. Methods Eng., 89, 8, 975-990, 2012 · Zbl 1242.76158 |
| [41] | Bailoor, S.; Annangi, A.; Seo, J. H., Fluid-structure interaction solver for compressible flows with applications to blast loading on thin elastic structures, Appl. Math. Model., 52, 470-492, 2017 · Zbl 1480.74062 |
| [42] | Kimata, H.; Taguchi, M.; Xue, X., Unsteady aerodynamic simulation of supersonic flexible parachute including fluid-structure interaction, (proceedings of the 23rd AIAA Aerodynamic Decelerator Systems Technology Conference, F, 2015), [C] |
| [43] | He, Y.; Li, D.; Liu, S., An immersed boundary method based on volume fraction, Proc. Eng., 99, 677-685, 2015 |
| [44] | Choung, H.; Saravanan, V.; Lee, S., Nonlinear weighting process in ghost-cell immersed boundary methods for compressible flow, J. Comput. Phys., 433, Article 110198 pp., 2021 · Zbl 1515.76112 |
| [45] | Zhang, Y.; Fang, X.; Zou, J., Numerical simulations of shock/obstacle interactions using an improved ghost-cell immersed boundary method, Comput. Fluids, 182, 128-143, 2019 · Zbl 1411.76109 |
| [46] | Mousavinia, S.; Emdad, H., Application of a new hybrid immersed boundary method for the simulation of viscous compressible flows, Int. J. Appl. Mech., 13, 05, Article 2150051 pp., 2021 |
| [47] | Keistler, P., An immersed boundary method for supersonic flow, (proceedings of the 46th AIAA Aerospace Sciences Meeting and Exhibit, F, 2008), [C] |
| [48] | Miyazaki, S., Cartesian grid-based CFD for shock/shock interference of two circular cylinders, (Proceedings of the AIAA AVIATION 2021 FORUM, F, 2021), C |
| [49] | Sha, S.; Chen, Z.; Jiang, X., Influences of obstacle geometries on shock wave attenuation, Shock Waves, 24, 573-582, 2014 |
| [50] | Seshadri, P. K.; De, A., Investigation of shock wave interactions involving stationary and moving wedges, Phys. Fluids, 32, 9, 2020 |
| [51] | Mao, C.; Jin, T.; Luo, K., Investigation of supersonic turbulent flows over a sphere by fully resolved direct numerical simulation, Phys. Fluids, 31, 5, 2019 |
| [52] | Luo, K.; Luo, Y.; Jin, T., Numerical analysis on shock-cylinder interaction using immersed boundary method, Sci. China Technol. Sci., 60, 9, 1423-1432, 2017 |
| [53] | Xue, X.; Nakamura, Y., Numerical simulation of a three-dimensional flexible parachute system under supersonic conditions, Transactions of the Japan Society for Aeronautical and Space Sciences, 11, 99-108, 2013, Aerospace Technology: Aerospace Technology Japan |
| [54] | Mizuno, Y.; Takahashi, S.; Nonomura, T., Direct numerical simulation of shock waves passed by multiple particles by using immersed boundary method, (Proceedings of the 54th AIAA Aerospace Sciences Meeting, F, 2016), [C] |
| [55] | Xiao, W.; Mao, C.; Jin, T., Fully resolved simulation of a shockwave interacting with randomly clustered particles via a ghost-cell immersed boundary method, Phys. Fluids, 32, 6, Article 066105 pp., 2020 |
| [56] | Kiris, C. C.; Housman, J. A.; Barad, M. F., Computational framework for launch, ascent, and vehicle aerodynamics (LAVA), Aerosp. Sci. Technol., 55, 189-219, 2016 |
| [57] | Kamatsuchi, T., Flow simulation around complex geometries with solution adaptive cartesian grid method, (Proceedings of the 18th AIAA Computational Fluid Dynamics Conference, F, 2007), [C] |
| [58] | Ghosh, S.; Choi, J.-I.; Edwards, J. R., Simulation of shock/boundary-layer interactions with bleed using immersed-boundary methods, J. Propuls. Power., 26, 2, 203-214, 2010 |
| [59] | Shallcross, G. S.; Fox, R. O.; Capecelatro, J., A volume-filtered description of compressible particle-laden flows, Int. J. Multiph. Flow, 122, Article 103138 pp., 2020 |
| [60] | Boukharfane, R.; Ribeiro, F. H.E.; Bouali, Z., A combined ghost-point-forcing/direct-forcing immersed boundary method (IBM) for compressible flow simulations, Comput. Fluids., 162, 91-112, 2018 · Zbl 1390.76554 |
| [61] | Yu, Z.; Shao, X., A direct-forcing fictitious domain method for particulate flows, J. Comput. Phys., 227, 1, 292-314, 2007 · Zbl 1280.76052 |
| [62] | Yu, Z.; Lin, Z.; Shao, X., A parallel fictitious domain method for the interface-resolved simulation of particle-laden flows and its application to the turbulent channel flow, Eng. Appl. Comput. Fluid Mech., 10, 1, 160-170, 2016 |
| [63] | Uhlmann, M., An immersed boundary method with direct forcing for the simulation of particulate flows, J. Comput. Phys., 209, 2, 448-476, 2005 · Zbl 1138.76398 |
| [64] | Silva, A. L.E.; Silveira-Neto, A.; Damasceno, J., Numerical simulation of two-dimensional flows over a circular cylinder using the immersed boundary method, J. Comput. Phys., 189, 2, 351-370, 2003 · Zbl 1061.76046 |
| [65] | Breugem, W.-P., A second-order accurate immersed boundary method for fully resolved simulations of particle-laden flows, J. Comput. Phys., 231, 13, 4469-4498, 2012 · Zbl 1245.76064 |
| [66] | Crowe, C.; Sommerfeld, M.; Tsuji, Y., Multiphase flows with droplets and particles, Ž, 1998 |
| [67] | Hoffmann, A. C.; Stein, L. E., Computational fluid dynamics, Gas Cyclones and Swirl Tubes, 123-135, 2002, Springer |
| [68] | Toro, E. F., The HLL and HLLC Riemann solvers, Riemann solvers and numerical methods for fluid dynamics, 315-344, 2009, Springer · Zbl 1227.76006 |
| [69] | Harten, A., High resolution schemes for hyperbolic conservation laws, J. Comput. Phys., 135, 2, 260-278, 1997 · Zbl 0890.65096 |
| [70] | Liu, X.-D.; Osher, S.; Chan, T., Weighted essentially non-oscillatory schemes, J. Comput. Phys., 115, 1, 200-212, 1994 · Zbl 0811.65076 |
| [71] | Tian, B.; Zeng, J.; Meng, B., Compressible multiphase particle-in-cell method (CMP-PIC) for full pattern flows of gas-particle system, J. Comput. Phys., 418, Article 109602 pp., 2020 · Zbl 07506167 |
| [72] | Kruggel-Emden, H.; Sturm, M.; Wirtz, S., Selection of an appropriate time integration scheme for the discrete element method (DEM), Comput. Chem. Eng., 32, 10, 2263-2279, 2008 |
| [73] | Kruggel-Emden, H.; Rickelt, S.; Wirtz, S., A study on the validity of the multi-sphere discrete element method, Powder Technol., 188, 2, 153-165, 2008 |
| [74] | Mo, H.; Lien, F. S.; Zhang, F., An immersed boundary method for solving compressible flow with arbitrarily irregular and moving geometry, Int. J. Numer. Methods Fluids, 88, 5, 239-263, 2018 |
| [75] | Luo, K.; Wang, Z.; Fan, J., Full-scale solutions to particle-laden flows: multidirect forcing and immersed boundary method, Phys. Rev. E, 76, 6, Article 066709 pp., 2007 |
| [76] | Matas, J.-P.; Morris, J. F., Guazzelli inertial migration of rigid spherical particles in Poiseuille flow, J. Fluid. Mech., 515, 171-195, 2004 · Zbl 1130.76301 |
| [77] | Li, J.; He, J.; Meng, B., Investigation of dust lifting by a moving shock wave based on compressible multiphase particle-in-cell method, Phys. Fluids, 34, 10, 2022 |
| [78] | Garg, R.; Galvin, J.; Li, T., Open-source MFIX-DEM software for gas-solids flows: Part I—verification studies, Powder Technol., 220, 122-137, 2012 |
| [79] | Chen, F.; Drumm, E. C.; Guiochon, G., Prediction/verification of particle motion in one dimension with the discrete-element method, Int. J. Geomech., 7, 5, 344-352, 2007 |
| [80] | Kaca, J., An interferometric investigation of the diffraction of a planar shock wave over a semicircular cylinder, Nasa Sti/Recon Tech. Rep. N, 89, 16126, 1988 |
| [81] | Tanno, H.; Itoh, K.; Saito, T., Interaction of a shock with a sphere suspended in a vertical shock tube, Shock Waves, 13, 3, 191-200, 2003 |
| [82] | Mo H., Lien F.-S., Zhang F., et al. A simple field function for solving complex and dynamic fluid-solid system on Cartesian grid. arXiv preprint arXiv:170202474, 2017. |
| [83] | Yoshihara, H.; Sacher, P., Test cases for inviscid flow field methods, AGARD advisory report AR-211, 1985 |
| [84] | Wei, F.; Jin, L.; Liu, J., GPU acceleration of a 2D compressible Euler solver on CUDA-based block-structured Cartesian meshes, J. Braz. Soc. Mech. Sci. Eng., 42, 1-12, 2020 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
