×

Efficient algorithms and implementations of entropy-based moment closures for rarefied gases. (English) Zbl 1376.76059

Summary: We present efficient algorithms and implementations of the 35-moment system equipped with the maximum-entropy closure in the context of rarefied gases. While closures based on the principle of entropy maximization have been shown to yield very promising results for moderately rarefied gas flows, the computational cost of these closures is in general much higher than for closure theories with explicit closed-form expressions of the closing fluxes, such as Grad’s classical closure. Following a similar approach as of C. K. Garrett et al. [ibid. 302, 573–590 (2015; Zbl 1349.65580)], we investigate efficient implementations of the computationally expensive numerical quadrature method used for the moment evaluations of the maximum-entropy distribution by exploiting its inherent fine-grained parallelism with the parallelism offered by multi-core processors and graphics cards. We show that using a single graphics card as an accelerator allows speed-ups of two orders of magnitude when compared to a serial CPU implementation. To accelerate the time-to-solution for steady-state problems, we propose a new semi-implicit time discretization scheme. The resulting nonlinear system of equations is solved with a Newton type method in the Lagrange multipliers of the dual optimization problem in order to reduce the computational cost. Additionally, fully explicit time-stepping schemes of first and second order accuracy are presented. We investigate the accuracy and efficiency of the numerical schemes for several numerical test cases, including a steady-state shock-structure problem.

MSC:

76P05 Rarefied gas flows, Boltzmann equation in fluid mechanics
76M12 Finite volume methods applied to problems in fluid mechanics

Citations:

Zbl 1349.65580

Software:

MKL; CUDA
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Grad, H., On the kinetic theory of rarefied gases, Commun. Pure Appl. Math., 2, 4, 331-407 (1949) · Zbl 0037.13104
[2] Grad, H., The profile of a steady plane shock wave, Commun. Pure Appl. Math., 5, 3, 257-300 (1952) · Zbl 0047.18801
[3] Brini, F., Hyperbolicity region in extended thermodynamics with 14 moments, Contin. Mech. Thermodyn., 13, 1, 1-8 (2001) · Zbl 0980.76076
[4] Torrilhon, M., Hyperbolic moment equations in kinetic gas theory based on multi-variate Pearson-IV-distributions, Commun. Comput. Phys., 7, 4, 639-673 (2010) · Zbl 1364.82050
[5] Cai, Z.; Fan, Y.; Li, R., On hyperbolicity of 13-moment system, Kinet. Relat. Models, 7, 3, 415-432 (2014) · Zbl 1304.82063
[6] Dreyer, W., Maximisation of the entropy in non-equilibrium, J. Phys. A, Math. Gen., 20, 18, 6505-6517 (1987) · Zbl 0633.76081
[7] Müller, I.; Ruggeri, T., Rational Extended Thermodynamics, Springer Tracts Nat. Philos., vol. 37 (1993), Springer-Verlag: Springer-Verlag New York · Zbl 0801.35141
[8] Levermore, C. D., Moment closure hierarchies for kinetic theories, J. Stat. Phys., 83, 5, 1021-1065 (1996) · Zbl 1081.82619
[9] \(M^c\) Donald, J. G.; Torrilhon, M., Affordable robust moment closures for CFD based on the maximum-entropy hierarchy, J. Comput. Phys., 251, 500-523 (2013) · Zbl 1349.82057
[10] Schaerer, R. P.; Torrilhon, M., A new closure mitigating the sub-shock phenomenon in the continuous shock-structure problem, PAMM, 13, 499-500 (2013)
[11] Schaerer, R. P.; Torrilhon, M., On singular closures for the 5-moment system in kinetic gas theory, Commun. Comput. Phys., 17, 2, 371-400 (2015) · Zbl 1373.76287
[12] Schaerer, R. P.; Torrilhon, M., The 35-moment system with the maximum-entropy closure for rarefied gas flows, Eur. J. Mech. B, Fluids (2017), in press · Zbl 1408.76459
[13] Garrett, C. K.; Hauck, C. D.; Hill, J., Optimization and large scale computation of an entropy-based moment closure, J. Comput. Phys., 302, 573-590 (2015) · Zbl 1349.65580
[14] Bhatnagar, P. L.; Gross, E. P.; Krook, M., A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems, Phys. Rev., 94, 3, 511-525 (1954) · Zbl 0055.23609
[15] Holway, L. H., New statistical models for kinetic theory: methods of construction, Phys. Fluids, 9, 1658-1673 (1966)
[16] Shakov, E. M., Generalization of the Krook kinetic relaxation equation, Fluid Dyn., 3, 95-96 (1968)
[17] Junk, M., Maximum entropy for reduced moment problems, Math. Models Methods Appl. Sci., 10, 7, 1001-1025 (2000) · Zbl 1012.44005
[18] Jaynes, E. T., Information theory and statistical mechanics, Phys. Rev., 106, 4, 620-630 (1957) · Zbl 0084.43701
[19] Jaynes, E. T., Information theory and statistical mechanics. II, Phys. Rev., 108, 2, 171-190 (1957) · Zbl 0084.43701
[20] Alldredge, G. W.; Hauck, C. D.; O’Leary, D. P.; Tits, A. L., Adaptive change of basis in entropy-based moment closures for linear kinetic equations, J. Comput. Phys., 258, 489-508 (2014) · Zbl 1349.82066
[21] \(Intel^®\) Corporation, \(Intel^®\) Math Kernel Library
[22] Lindholm, E.; Nickolls, J.; Oberman, S.; Montrym, J., NVIDIA Tesla: a unified graphics and computing architecture, IEEE MICRO, 28, 2, 39-55 (2008)
[23] Nickolls, J.; Buck, I.; Garland, M.; Skadron, K., Scalable parallel programming with CUDA, ACM Queue, 6, 2, 40-53 (2008)
[24] NVIDIA Corporation, CUDA C programming guide (2015)
[25] Saini, S.; Chang, J.; Jin, H., Performance evaluation of the Intel Sandy Bridge based NASA Pleiades using scientific and engineering applications, (High Performance Computing SystemsPerformance Modeling, Benchmarking and Simulation. High Performance Computing SystemsPerformance Modeling, Benchmarking and Simulation, 4th International Workshop (2014), Springer International Publishing), 25-51
[26] Williams, S.; Waterman, A.; Patterson, D., Roofline: an insightful visual performance model for multicore architectures, Commun. ACM, 52, 4, 65-76 (2009)
[27] Torrilhon, M.; Struchtrup, H., Regularized 13-moment-equations: shock structure calculations and comparison to Burnett models, J. Fluid Mech., 513, 171-198 (2004) · Zbl 1107.76069
[28] Groth, C. P.T.; McDonald, J. G., Towards physically realizable and hyperbolic moment closures for kinetic theory, Contin. Mech. Thermodyn., 21, 6, 467-493 (2009) · Zbl 1184.76839
[29] Alldredge, G. W.; Hauck, C. D.; Tits, A. L., High-order entropy-based closures for linear transport in slab geometry II: a computational study of the optimization problem, SIAM J. Sci. Comput., 34, 4, B361-B391 (2012) · Zbl 1297.82032
[30] Boillat, G.; Ruggeri, T., On the shock structure problem for hyperbolic system of balance laws and convex entropy, Contin. Mech. Thermodyn., 10, 5, 285-292 (1998) · Zbl 0922.76237
[31] Vuyst, F. D., Stable and accurate hybrid finite volume methods based on pure convexity arguments for hyperbolic systems of conservation law, J. Comput. Phys., 193, 2, 426-468 (2004) · Zbl 1039.35010
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.