Direct numerical simulation of laminar-turbulent transition on grids with local refinement

D. I. Okhotnikov

Učënye zapiski Kazanskogo universiteta. Seriâ Fiziko-matematičeskie nauki, Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-Matematicheskie Nauki, Tome 159 (2017) no. 2, pp. 216-230

Résumé

This paper considers a method of building composite grids with local refinement for direct numerical simulation of flows in a channel with transverse cylindrical rib at moderate Reynolds numbers. The specific characteristic of the flows is LTT (laminar-turbulent transition) in a trace behind the rib. Grids characterized by resolution of the whole spectrum of vortices are needed for numerical simulation of LTT, thereby leading to significant computational costs and a need for supercomputers. The purpose of the study is to generate a cost-efficient computational grid based on the above-described method, as well as to test this grid for flows in a channel with the following two types of ribs: a semicircular ridge on the lower wall of the channel and cylinder near it. The grid has been built using the HybMesh generator, which is based on the composite approach, when the grid is built from a set of structured subdomains, in which grid generation and refinement is defined by several control parameters. As a result of superposition of the grids built in subdomains, a composite grid from structured parts superposed with unstructured inserts has been obtained. The results of the direct numerical simulation on a grid for the problem with a semicircular ridge have turned out to be in a good agreement with the results of the laboratory experiment, both along the profiles of average velocity and the profiles of pulsations. The results obtained on the grid for the problem with a circular cylinder have demonstrated good agreement with the results of the numerical simulation on high-definition grids, predicting in the right way the transition point to three-dimensionality in a trace behind the cylinder. The method of building composite grids with the use of the HybMesh generator allows to generate cost-efficient composite grids taking into account all flow peculiarities and allowing for calculations using medium-powered computers.

Détail

Keywords: flow in channel, direct numerical simulation, composite grids, local refinement.
Mots-clés : transverse rib, laminar-turbulent transition

Bibliographie

[1] Molochnikov V. M., Mazo A. B., Malyukov A. V., Kalinin E. I., Mikheev N. I., Dushina O. A., Paereliy A. A., “Distinctive features of vortical structures generation in separated channel flow behind a rib under transition to turbulence”, Thermophys. Aeromech., 21:3 (2014), 309–317 | DOI

[2] Dushina O. A., Molochnikov V. M., Mikheev N. I., Paerely A. A., “Evolution of kinematic structure of the flow behind a transverse rib for transitional flow regimes”, Thermophys. Aeromech., 19:2 (2012), 259–266 | DOI

[3] Mazo A. B., Okhotnikov D. I., “Local transition to turbulence behind an obstacle for a nominally laminar flow”, Lobachevskii J. Math., 37:3 (2015), 360–367 | DOI | MR

[4] Mazo A. B., Simulation of Turbulent Flows of Incompressible Fluid, Kazan. Gos. Univ., Kazan, 2007, 106 pp. (In Russian)

[5] Garbaruk A. V., Strelez M. H., Travin A. K., Shur M. L., Modern Approaches to Turbulence Modeling, Polytekh. Univ., St. Petersburg, 2016, 234 pp. (In Russian)

[6] Kalinin E. I., Mazo A. B., Isaev S. A., “Composite mesh generator for CFD problems”, IOP Conf. Ser.: Mater. Sci. Eng., 158:1 (2016), Art. 012047, 6 pp. | DOI

[7] Schlichting H., Gersten K., Boundary-Layer Theory, Springer, Berlin–Heidelberg, 2017, 805 pp. | MR | Zbl

[8] Fletcher C., Computational techniques for fluid dynamics, v. 2, Springer-Verlag, N.Y., 1988, 409 pp. | MR | Zbl

[9] ANSYS Fluent Tutorial Guide, Ansys Inc., Southpointe, 2011, 1146 pp.

[10] Harri H., Tri Q., Antti R., “Conjugate function method for numerical conformal mappings”, J. Comput. Appl. Math., 237:1 (2013), 340–353 | DOI | MR | Zbl

[11] Brown P. R., Porter R. M., “Conformal mapping of circular quadrilaterals and weierstrass elliptic functions”, Comp. Methods Funct. Theory, 11:2 (2011), 463–486 | DOI | MR | Zbl

[12] Liseikin V. D., Grid Generation Methods, Ed. 2, Springer, 2010, xviii+390 pp. | DOI | MR | Zbl

[13] Liseikin V. D., A computational differential geometry approach to grid generation, Springer-Verlag, Berlin–Heidelberg, 2007, xiv+294 pp. | DOI | MR | Zbl

[14] Geuzaine C., Remacle J.-F., “Gmsh: a three-dimensional finite element mesh generator with built-in pre- and post-processing facilities”, Int. J. Numer. Methods Eng., 79:11 (2009), 1309–1331 | DOI | MR | Zbl

[15] Remacle J.-F., Geuzaine C., Compere G., Marchandise E., “High-quality surface remeshing using harmonic maps”, Int. J. Numer. Methods Eng., 83:4 (2010), 403–425 | DOI | MR | Zbl

[16] Mikheev N. I., Dushin N. S., “A method for measuring the dynamics of velocity vector fields in a turbulent flow using smoke image-visualization videos”, Instrum. Exp. Tech., 59:6 (2016), 882–889 | DOI

[17] Zaripov D. I., Aslaev A. K., Miheev N. I., Dushin N. S., “Estimation of accuracy of new method of instantaneous flow velocity field measurement”, Tr. Academenergo, 2016, no. 1, 42–52 (In Russian)

[18] Nicolas K., Dimokratis G., Starvos K., “Three dimensional flow around a circular cylinder confined in a plane channel”, Phys. Fluids, 23:6 (2016), Art. 064106, 14 pp. | DOI

[19] HybMesh Generator, URL: https://github.com/kalininei/HybMesh

[20] Singha S., Sinhamahapatra K. P., “Flow past a circular cylinder between parallel walls at low Reynolds numbers”, Ocean Eng., 37:8–9 (2010), 757–769 | DOI