zbMATH — the first resource for mathematics

Higher-order compact scheme for high-performance computing of stratified rotating flows. (English) Zbl 1410.76278
Summary: To take advantage of modern generation computing hardware, a scalable numerical method, based on higher-order compact scheme, is described to solve rotating stratified flows in cylindrical annular domains. An original approach combining 2d-pencil decomposition and reduced parallel diagonal dominant is proposed to improve the parallelization performance during the computation of Poisson/Helmholtz solvers and time explicit terms. The developed technique is validated with respect to analytical solutions, using the method of manufactured solutions, and available data for two specific configurations. The purpose is to demonstrate its ability to correctly capture the flow characteristics in strato-rotational instability and in baroclinic instability with associated small-scale features. Moreover, this code is found to drastically reduce the huge execution times often preventing detailed numerical investigations of these complex phenomena. Strong scaling test is carried out to assess the performance for up to 1024 cores using grid up to \(128 \times 568 \times 568\) in radial, axial and azimuthal directions.

76M20 Finite difference methods applied to problems in fluid mechanics
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
76D05 Navier-Stokes equations for incompressible viscous fluids
76D50 Stratification effects in viscous fluids
76U05 General theory of rotating fluids
Full Text: DOI
[1] Yakovlev, S.; Moxey, D.; Kirby, R. M.; Sherwin, S. J., To CG or to HDG: a comparative study in 3d, J Sci Comput, 67, 1, 192-220, (2016) · Zbl 1339.65225
[2] Witherden, F. D.; Farrington, A. M.; Vincent, P. E., PYFR: an open source framework for solving advection-diffusion type problems on streaming architectures using the flux reconstruction approach, Comput Phys Commun, 185, 11, 3028-3040, (2014) · Zbl 1348.65005
[3] King, J.; Yakovlev, S.; Fu, Z.; Kirby, R. M.; Sherwin, S. J., Exploiting batch processing on streaming architectures to solve 2d elliptic finite element problems: a hybridized discontinuous Galerkin (HDG) case study, J Sci Comput, 60, 2, 457-482, (2014) · Zbl 1304.65250
[4] Moxey, D.; Cantwell, C.; Kirby, R.; Sherwin, S., Optimising the performance of the spectral/hp element method with collective linear algebra operations, Comput Methods Appl Mech Eng, 310, 628-645, (2016)
[5] Karniadakis, G.; Sherwin, S., Spectral/hp element methods for computational fluid dynamics, (2013), Oxford University Press · Zbl 1256.76003
[6] Cantwell, C. D.; Moxey, D.; Comerford, A.; Bolis, A.; Rocco, G.; Mengaldo, G., Nektar++: an open-source spectral/hp element framework, Comput Phys Commun, 192, 205-219, (2015) · Zbl 1380.65465
[7] Fischer, P.; Kruse, J.; Mullen, J.; Tufo, H.; Lottes, J.; Kerkemeier, S., Nek5000: open source spectral element CFD solver, (2008), Argonne National Laboratory, Mathematics and Computer Science Division, Argonne, IL
[8] Lombard, J.-E. W.; Moxey, D.; Sherwin, S. J.; Hoessler, J. F.; Dhandapani, S.; Taylor, M. J., Implicit large-eddy simulation of a wingtip vortex, AIAA J, 54, 2, 506-518, (2015)
[9] Kerkemeier, S.; Parker, S.; Fischer, P., Scalability of the nek5000 spectral element code, Technical Report, (2010), Jülich blue gene/p extreme scaling workshop
[10] Abide, S.; Binous, M. S.; Zeghmati, B., An efficient parallel high-order compact scheme for the 3d incompressible Navier-Stokes equations, Int J Comut Fluid Dyn, 31, 4-5, 214-229, (2017)
[11] Gellert, M.; Rüdiger, G., Stratorotational instability in Taylor-Couette flow heated from above, J Fluid Mech, 623, 375-385, (2009) · Zbl 1157.76325
[12] Hide, R., An experimental study of thermal convection in a rotating liquid, Phil Trans R Soc Lond A, 250, 983, 441-478, (1958)
[13] Randriamampianina, A.; Crespo del Arco, E., Inertia-gravity waves in a liquid-filled, differentially heated, rotating annulus, J Fluid Mech, 782, 144-177, (2015) · Zbl 1381.76389
[14] Schornbaum, F.; Rude, U., Massively parallel algorithms for the lattice Boltzmann method on nonuniform grids, SIAM J Sci Comput, 38, 2, C96-C126, (2016) · Zbl 1381.76301
[15] Vàzquez, M.; Houzeaux, G.; Koric, S.; Artigues, A.; Aguado-Sierra, J.; Ars, R., Alya: multiphysics engineering simulation toward exascale, J Comput Sci, 14, 15-27, (2016)
[16] Baker, A. H.; Falgout, R. D.; Gamblin, T.; Kolev, T. V.; Schulz, M.; Yang, U. M., Scaling algebraic multigrid solvers: on the road to exascale, Competence in high performance computing, 215-226, (2011), Springer
[17] Notay, Y., An aggregation-based algebraic multigrid method, Electron Trans Numer Anal, 37, 6, 123-146, (2010) · Zbl 1206.65133
[18] Laizet, S.; Lamballais, E.; Vassilicos, J., A numerical strategy to combine high-order schemes, complex geometry and parallel computing for high resolution DNS of fractal generated turbulence, Comput Fluids, 39, 3, 471-484, (2010) · Zbl 1242.76078
[19] He, P., A high order finite difference solver for massively parallel simulations of stably stratified turbulent channel flows, Comput Fluids, 127, 161-173, (2016) · Zbl 1390.76572
[20] Thais, L.; Tejada-Martnez, A. E.; Gatski, T. B.; Mompean, G., A massively parallel hybrid scheme for direct numerical simulation of turbulent viscoelastic channel flow, Computers & Fluids, 43, 1, 134-142, (2011) · Zbl 1452.76075
[21] Chatterjee, A. G.; Verma, M. K.; Kumar, A.; Samtaney, R.; Hadri, B.; Khurram, R., Scaling of a fast Fourier transform and a pseudo-spectral fluid solver up to 196608 cores, J Parallel Distrib Comput, 113, 77-91, (2018)
[22] Shi, L.; Rampp, M.; Hof, B.; Avila, M., A hybrid mpi-openmp parallel implementation for pseudospectral simulations with application to taylorcouette flow, Comput Fluids, 106, 1-11, (2015) · Zbl 1390.76623
[23] Mininni, P. D.; Rosenberg, D.; Reddy, R.; Pouquet, A., A hybrid MPI-openmp scheme for scalable parallel pseudospectral computations for fluid turbulence, Parallel Comput, 37, 6, 316-326, (2011)
[24] Lele, S. K., Compact finite difference schemes with spectral-like resolution, J Comput Phys, 103, 1, 16-42, (1992) · Zbl 0759.65006
[25] Sun, X.-H., Application and accuracy of the parallel diagonal dominant algorithm, Parallel Comput, 21, 8, 1241-1267, (1995) · Zbl 0873.65016
[26] Haldenwang, P.; Labrosse, G.; Abboudi, S.; Deville, M., Chebyshev 3d spectral and 2d pseudospectral solvers for the Helmholtz equation, J Comput Phys, 55, 115-128, (1984) · Zbl 0544.65071
[27] Hugues, S.; Randriamampianina, A., An improved projection scheme applied to pseudospectral methods for the incompressible Navier-Stokes equations, Int J Numer Methods Fluids, 28, 501-521, (1998) · Zbl 0932.76065
[28] Le Bars, M.; Le Gal, P., Experimental analysis of the stratorotational instability in a cylindrical Couette flow, Phys Rev Lett, 99, 6, 064502, (2007)
[29] Rüdiger, G.; Hollerbach, R.; Kitchatinov, L. L., Magnetic processes in astrophysics: theory, simulations, experiments, (2013), John Wiley & Sons
[30] Ibanez, R.; Swinney, H. L.; Rodenborn, B., Observations of the stratorotational instability in rotating concentric cylinders, Phys Rev Fluids, 1, 053601, (2016)
[31] Rüdiger, G.; Seelig, T.; Schultz, M.; Gellert, M.; Egbers, C.; Harlander, U., The stratorotational instability of Taylor-Couette flows with moderate Reynolds numbers, Geophys Astrophys Fluid Dyn, 111, 6, 429-447, (2017)
[32] Pierrehumbert, R.; Swanson, K., Baroclinic instability, Annu Rev Fluid Mech, 27, 1, 419-467, (1995)
[33] Orlanski, I.; Cox, M. D., Baroclinic instability in Ocean currents, Geophys Astrophys Fluid Dyn, 4, 1, 297-332, (1972)
[34] Read, P. L.; Collins, M.; Früh, W.-G.; Lewis, S. R.; Lovegrove, A. F., Wave interactions and baroclinic chaos: a paradigm for long timescale variability in planetary atmospheres, Chaos Solitons Fractals, 9, 1-2, 231-249, (1998)
[35] Read, P. L., Transition to geostrophic turbulence in the laboratory, and as a paradigm in atmospheres and oceans, Surv Geophys, 22, 3, 265-317, (2001)
[36] Fritts, D. C.; Alexander, M. J., Gravity wave dynamics and effects in the middle atmosphere, Rev Geophys, 41, 1, (2003)
[37] Borchert, S.; Achatz, U.; Fruman, M. D., Gravity wave emission in an atmosphere-like configuration of the differentially heated rotating annulus experiment, J Fluid Mech, 758, 287-311, (2014)
[38] Hien, S.; Rolland, J.; Borchert, S.; Schoon, L.; Zülicke, C.; Achatz, U., Spontaneous inertia-gravity wave emission in the differentially heated rotating annulus experiment, J Fluid Mech, 838, 5-41, (2018)
[39] Von Larcher, T.; Viazzo, S.; Harlander, U.; Vincze, M.; Randriamampianina, A., Instabilities and small-scale waves within the Stewartson layers of a thermally driven rotating annulus, J Fluid Mech, 841, 380-407, (2018)
[40] Raspo, I.; Hugues, S.; Serre, E.; Randriamampianina, A.; Bontoux, P., A spectral projection method for the simulation of complex three - dimensional rotating flows, Comput Fluids, 31, 4-7, 745-767, (2002) · Zbl 1028.76036
[41] Abide, S.; Viazzo, S., A 2d compact fourth-order projection decomposition method, J Comput Phys, 206, 1, 252-276, (2005) · Zbl 1087.76083
[42] Knikker, R., Study of a staggered fourth-order compact scheme for unsteady incompressible viscous flows, Int J Numer Methods Fluids, 59, 10, 1063-1092, (2009) · Zbl 1158.76380
[43] Schiestel, R.; Viazzo, S., A Hermitian-Fourier numerical method for solving the incompressible Navier-Stokes equations, Comput fluids, 24, 6, 739-752, (1995) · Zbl 0845.76060
[44] Abide, S.; Chesneau, X.; Zeghmati, B., A fourth-order iterative solver for the singular Poisson equation, (Dimov, I.; Faragó, I.; Vulkov, L., Numerical analysis and its applications, (2013), Springer Berlin Heidelberg Berlin, Heidelberg), 143-150 · Zbl 1352.65416
[45] Pekurovsky, D., P3DFFT: a framework for parallel computations of Fourier transforms in three dimensions, SIAM J Sci Comput, 34, 4, C192-C209, (2012) · Zbl 1253.65205
[46] Li, N.; Laizet, S., 2DECOMP&FFT-a highly scalable 2d decomposition library and FFT interface, Proceedings of the Cray User Group conference, 1-13, (2010)
[47] Roache, P. J., Code verification by the method of manufactured solutions, J Fluids Eng, 124, 1, 4-10, (2002)
[48] Harlander, U.; von Larcher, T.; Wang, Y.; Egbers, C., PIV- and LDV-measurements of baroclinic wave interactions in a thermally driven rotating annulus, Exp Fluids, 51, 1, 37-49, (2011)
[49] von Larcher, T.; Dörnbrack, A., Numerical simulations of baroclinic driven flows in a thermally driven rotating annulus using the immersed boundary method, Meteorol Z, 1-12, (2014)
[50] Vincze, M.; Borchert, S.; Achatz, U.; von Larcher, T.; Baumann, M.; Liersch, C., Benchmarking in a rotating annulus: a comparative experimental and numerical study of baroclinic wave dynamics, Meteorol Z, 23, 6, 611-635, (2015)
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.