×

Fully decoupled monolithic projection method for natural convection problems. (English) Zbl 1375.76178

Summary: To solve time-dependent natural convection problems, we propose a fully decoupled monolithic projection method. The proposed method applies the Crank-Nicolson scheme in time and the second-order central finite difference in space. To obtain a non-iterative monolithic method from the fully discretized nonlinear system, we first adopt linearizations of the nonlinear convection terms and the general buoyancy term with incurring second-order errors in time. Approximate block lower-upper decompositions, along with an approximate factorization technique, are additionally employed to a global linearly coupled system, which leads to several decoupled subsystems, i.e., a fully decoupled monolithic procedure. We establish global error estimates to verify the second-order temporal accuracy of the proposed method for velocity, pressure, and temperature in terms of a discrete \(\ell^2\)-norm. Moreover, according to the energy evolution, the proposed method is proved to be stable if the time step is less than or equal to a constant. In addition, we provide numerical simulations of two-dimensional Rayleigh-Bénard convection and periodic forced flow. The results demonstrate that the proposed method significantly mitigates the time step limitation, reduces the computational cost because only one Poisson equation is required to be solved, and preserves the second-order temporal accuracy for velocity, pressure, and temperature. Finally, the proposed method reasonably predicts a three-dimensional Rayleigh-Bénard convection for different Rayleigh numbers.

MSC:

76R10 Free convection
76M20 Finite difference methods applied to problems in fluid mechanics
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Patterson, J.; Imberger, J., Unsteady natural convection in a rectangular cavity, J. Fluid Mech., 100, 65-86 (1980) · Zbl 0433.76071
[2] De Vahl Davis, G., Natural convection of air in a square cavity: a bench mark numerical solution, Int. J. Numer. Methods Fluids, 3, 249-264 (1983) · Zbl 0538.76075
[3] Le Quéré, P.; Alziary de Roquefort, T., Computation of natural convection in two-dimensional cavities with Chebyshev polynomials, J. Comput. Phys., 57, 210-228 (1985) · Zbl 0585.76128
[4] Le Quéré, P., Accurate solutions to the square thermally driven cavity at high Rayleigh number, Comput. Fluids, 20, 29-41 (1991) · Zbl 0731.76054
[5] Kerr, R. M., Rayleigh number scaling in numerical convection, J. Fluid Mech., 310, 139-179 (1996) · Zbl 0875.76584
[6] Wan, D. C.; Patnaik, B. S.V.; Wei, G. W., A new benchmark quality solution for the buoyancy-driven cavity by discrete singular convolution, Numer. Heat Transf., Part B, 40, 199-228 (2001)
[7] Armfield, S.; Street, R., An analysis and comparison of the time accuracy of fractional-step methods for the Navier-Stokes equations on staggered grids, Int. J. Numer. Methods Fluids, 38, 255-282 (2002) · Zbl 1027.76033
[8] Ouertatani, N.; Cheikh, N. B.; Beya, B. B.; Lili, T., Numerical simulation of two-dimensional Rayleigh-Bénard convection in an enclosure, C. R., Méc., 336, 464-470 (2008) · Zbl 1143.76483
[9] Löwe, J.; Lube, G., A projection-based variational multiscale method for large-eddy simulation with application to non-isothermal free convection problems, Math. Models Methods Appl. Sci., 22, Article 1150011 pp. (2012), (31 pages) · Zbl 1426.76287
[10] Zhang, K.; Yang, M.; Zhang, Y., A compact finite-difference scheme based on the projection method for natural-convection heat transfer, Numer. Heat Transf., Part B, 61, 259-278 (2012)
[11] Deteix, J.; Jendoubi, A.; Yakoubi, D., A coupled prediction scheme for solving the Navier-Stokes and convection-diffusion equations, SIAM J. Numer. Anal., 52, 2415-2439 (2014) · Zbl 1307.76072
[12] Zhang, Y.; Hou, Y.; Zhao, J., Error analysis of a fully discrete finite element variational multiscale method for the natural convection problem, Comput. Math. Appl., 68, 543-567 (2014) · Zbl 1362.76056
[13] Park, S.; Lee, C., Analysis of coherent structures in Rayleigh-Bénard convection, J. Turbul., 16, 1162-1178 (2015)
[14] Pan, X.; Kim, K.; Lee, C.; Choi, J.-I., A decoupled monolithic projection method for natural convection problems, J. Comput. Phys., 314, 160-166 (2016) · Zbl 1349.76518
[15] Qian, Y.; Zhang, T., The second order projection method in time for the time-dependent natural convection problem, Appl. Math., 61, 299-315 (2016) · Zbl 1424.76023
[16] Chorin, A. J., Numerical solution of the Navier-Stokes equations, Math. Comput., 22, 745-762 (1968) · Zbl 0198.50103
[17] Kim, J.; Moin, P., Application of a fractional-step method to incompressible Navier-Stokes equations, J. Comput. Phys., 59, 308-323 (1985) · Zbl 0582.76038
[18] Bell, J. B.; Colella, P.; Glaz, H. M., A second-order projection method for the incompressible Navier-Stokes equations, J. Comput. Phys., 85, 257-283 (1989) · Zbl 0681.76030
[19] Perot, J., An analysis of the fractional step method, J. Comput. Phys., 108, 51-58 (1993) · Zbl 0778.76064
[20] Choi, H.; Moin, P., Effects of the computational time step on numerical solutions of turbulent flow, J. Comput. Phys., 113, 1-4 (1994) · Zbl 0807.76051
[21] Shen, J., On error estimates of the projection methods for the Navier-Stokes equations: second-order schemes, Math. Comput., 65, 1039-1065 (1996) · Zbl 0855.76049
[22] Kim, K.; Baek, S. J.; Sung, H. J., An implicit velocity decoupling procedure for the incompressible Navier-Stokes equations, Int. J. Numer. Methods Fluids, 38, 125-138 (2002) · Zbl 1059.76046
[23] Pan, X.; Lee, C.; Kim, K.; Choi, J.-I., Analysis of velocity-components decoupled projection method for the incompressible Navier-Stokes equations, Comput. Math. Appl., 71, 1722-1743 (2016) · Zbl 1443.65137
[24] Chen, Y.; Shen, J., Efficient, adaptive energy stable schemes for the incompressible Cahn-Hilliard Navier-Stokes phase-field models, J. Comput. Phys., 308, 40-56 (2016) · Zbl 1352.65229
[25] Moin, P., Fundamentals of Engineering Numerical Analysis, 102-104 (2010), Cambridge University Press: Cambridge University Press New York · Zbl 1228.65003
[26] van Kan, J., A second-order accurate pressure-correction scheme for viscous incompressible flow, SIAM J. Sci. Stat. Comput., 7, 870-891 (1986) · Zbl 0594.76023
[27] Jangveladze, T.; Kiguradze, Z.; Neta, B., Large time asymptotic and numerical solution of a nonlinear diffusion model with memory, Comput. Math. Appl., 59, 254-273 (2010) · Zbl 1189.65195
[28] Wanschura, M.; Kuhlmann, H. C.; Rath, H. J., Three-dimensional instability of axisymmetric buoyant convection in cylinders heated from below, J. Fluid Mech., 326, 399-415 (1996) · Zbl 0884.76020
[29] Heywood, J. G.; Rannacher, R., Finite element approximation of the nonstationary Navier-Stokes problem. Part IV. Error analysis for second-order time discretization, SIAM J. Numer. Anal., 27, 353-384 (1990) · Zbl 0694.76014
[30] Marion, M.; Temam, R., Navier-Stokes equations: theory and approximation, (Handbook of Numerical Analysis, vol. 6 (1998), North-Holland: North-Holland Amsterdam), 503-688 · Zbl 0921.76040
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.