Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

doi:10.1007/s10483-015-1964-7

Applied Mathematics and Mechanics (English Edition) ›› 2015, Vol. 36 ›› Issue (8): 1073-1090.doi: https://doi.org/10.1007/s10483-015-1964-7

Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

Benwen LI¹, Shangshang CHEN²

1. Institute of Thermal Engineering, School of Energy and Power Engineering, Dalian University of Technology, Dalian 116024, Liaoning Province, China;
2. Key Laboratory of Electromagnetic Processing of Materials(Ministry of Education), Northeastern University, Shenyang 110819, China

收稿日期:2014-12-24 修回日期:2015-05-22 出版日期:2015-08-01 发布日期:2015-08-01
通讯作者: Benwen LI E-mail:heatli@dlut.edu.cn
基金资助:
Project supported by the National Natural Science Foundation of China (No. 51176026) and the Fundamental Research Funds for the Central Universities (No.DUT14RC(3)029)

Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

Benwen LI¹, Shangshang CHEN²

1. Institute of Thermal Engineering, School of Energy and Power Engineering, Dalian University of Technology, Dalian 116024, Liaoning Province, China;
2. Key Laboratory of Electromagnetic Processing of Materials(Ministry of Education), Northeastern University, Shenyang 110819, China

Received:2014-12-24 Revised:2015-05-22 Online:2015-08-01 Published:2015-08-01
Contact: Benwen LI E-mail:heatli@dlut.edu.cn
Supported by:
Project supported by the National Natural Science Foundation of China (No. 51176026) and the Fundamental Research Funds for the Central Universities (No.DUT14RC(3)029)

摘要/Abstract

摘要： An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and incompressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangular sub-domains. In each sub-domain, an improved projection scheme with second-order accuracy is used to deal with the coupling of velocity and pressure, and the Chebyshev collocation spectral method (CSM) is adopted to execute the spatial discretization. The influence matrix technique is employed to enforce the continuities of both variables and their normal derivatives between the adjacent sub-domains. The imposing of the Neumann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet boundary conditions on interface and using the first-order normal derivatives as transmission conditions to keep the continuities of variables is proposed to overcome this trouble. Three test cases are used to verify the accuracy and efficiency, and the detailed comparison between the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.

关键词: domain decomposition, influence matrix technique, Chebyshev collocation spectral method, incompressible Navier-Stokes equation

Abstract: An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and incompressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangular sub-domains. In each sub-domain, an improved projection scheme with second-order accuracy is used to deal with the coupling of velocity and pressure, and the Chebyshev collocation spectral method (CSM) is adopted to execute the spatial discretization. The influence matrix technique is employed to enforce the continuities of both variables and their normal derivatives between the adjacent sub-domains. The imposing of the Neumann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet boundary conditions on interface and using the first-order normal derivatives as transmission conditions to keep the continuities of variables is proposed to overcome this trouble. Three test cases are used to verify the accuracy and efficiency, and the detailed comparison between the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.

Key words: Chebyshev collocation spectral method, domain decomposition, influence matrix technique, incompressible Navier-Stokes equation

中图分类号:

Benwen LI, Shangshang CHEN. Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations[J]. Applied Mathematics and Mechanics (English Edition), 2015, 36(8): 1073-1090.

参考文献

[1] Canuto, C., Hussaini, M. Y., Quarteroni, A., and Zang, T. A. Spectral Methods in Fluid Dynamics, Springer, New York (1988)
[2] Peyret, R. Spectral Methods for Incompressible Viscous Flow, Springer, New York (2002)
[3] Orszag, S. A. Spectral methods for problems in complex geometries. Journal of Computational Physics, 37(1), 70-92(1980)
[4] Canuto, C. and Funaro, D. The Schwarz algorithm for spectral methods. SIAM Journal on Numerical Analysis, 25(1), 24-40(1988)
[5] Pateraa, A. T. A spectral element method for fluid dynamics:laminar flow in a channel expansion. Journal of Computational Physics, 54(3), 468-488(1984)
[6] Phillips, T. N. and Karageorghis, A. Efficient direct method for solving the spectral collocation equations for stokes flow in rectangularly decomposable domains. SIAM Journal on Scientific and Statistical Computing, 10(1), 89-103(1989)
[7] Macaraeg, M. G. and Streett, C. L. Improvement in spectral collocation discretization through a multiple domain technique. Applied Numerical Mathematics, 2(2), 95-108(1986)
[8] Zanolli, P. Domain decomposition algorithms for spectral methods. Calcolo, 24, 201-240(1987)
[9] Funaro, D., Quarteroni, A., and Zanolli, P. An iterative procedure with interface relaxation for domain decomposition methods. SIAM Journal on Numerical Analysis, 25(6), 1213-1236(1988)
[10] Louchart, O., Randriamampianina, A., and Leonardi, E. Spectral domain decomposition technique for the incompressible Navier-Stokes equations. Numerical Heat Transfer, Part A, 34(5), 495-518(1998)
[11] Louchart, O. and Randriamampianina, A. A spectral iterative domain decomposition technique for the incompressible Navier-Stokes equations. Applied Numerical Mathematics, 33(1), 233-240(2000)
[12] Pulicani, J. P. A spectral multi-domain method for the solution of 1-D-Helmholtz and Stokes-type equations. Computers & Fluids, 16(2), 207-215(1988)
[13] Lacroix, J. M., Peyret, R., and Pulicani, J. P. A pseudospectral multi-domain method for the Navier-Stokes equations with application to double-diffusive convection. GAMM Conference on Numerical Methods in Fluid Mechanics, 167-174(1988)
[14] Raspo, I. and Bontoux, P. Direct multidomain spectral method for the computation of various fluid dynamic problems. Journal of Applied Mathematics and Mechanics, 79(1), 33-36(1999)
[15] Raspo, I., Ouazzani, J., and Peyret, R. A spectral multidomain technique for the computation of the Czochralskimelt configuration. International Journal of Numerical Methods for Heat & Fluid Flow, 6(1), 31-58(1996)
[16] Bwemba, R. and Pasquetti, R. On the influence matrix used in the spectral solution of the 2D Stokes problem (vorticity-stream function formulation). Applied Numerical Mathematics, 16(3), 299-315(1995)
[17] Boguslawski, A. and Kubacki, S. An influence matrix technique for multi-domain solution of the Navier-Stokes equations in a vorticity-streamfunction formulation. Journal of Theoretical and Applied Mechanics, 47(1), 17-40(2009)
[18] Raspo, I. A direct spectral domain decomposition method for the computation of rotating flows in a T-shape geometry. Computers & Fluids, 32, 431-456(2003)
[19] Danabasoglu, G., Biringen, S., and Streett, C. L. Application of the spectral multidomain method to the Navier-Stokes equations. Journal of Computational Physics, 113(2), 155-164(1994)
[20] Sabbah, C. and Pasquetti, R. A divergence-free multidomain spectral solver of the Navier-Stokes equations in geometries of high aspect ratio. Journal of Computational Physic, 139(2), 359-379(1998)
[21] Droll, P. and Schäfer, M. A pseudospectral multi-domain method for the incompressible Navier-Stokes equations. Journal of Scientific Computing, 17, 365-374(2002)
[22] Chen, S. S. and Li, B. W. Application of collocation spectral domain decomposition method to solve radiative heat transfer in 2D partitioned domains. Journal of Quantitative Spectroscopy & Radiative Transfer, 149, 275-284(2014)
[23] Vay, J. L., Haber, I., and Godfrey, B. B. A domain decomposition method for pseudo-spectral electromagnetic simulations of plasmas. Journal of Computational Physic, 243, 260-268(2013)
[24] Grinberg, L. and Karniadakis, G. E. A new domain decomposition method with overlapping patches for ultrascale simulations:application to biological flows. Journal of Computational Physic, 229(15), 5541-5563(2010)
[25] Sidler, R., Carcione, J. M., and Holliger, K. A pseudo-spectral method for the simulation of poroelastic seismic wave propagation in 2D polar coordinates using domain decomposition. Journal of Computational Physic, 235, 846-864(2013)
[26] Golbabai, A. and Javidi, M. A spectral domain decomposition approach for the generalized Burger's-Fisher equation. Chaos, Solitons and Fractals, 39(1), 385-392(2009)
[27] Taleei, A. and Dehghan, M. Time-splitting pseudo-spectral domain decomposition method for the soliton solutions of the one- and multi-dimensional nonlinear Schrödinger equations. Computer Physics Communications, 185(6), 1515-1528(2014)
[28] Hameed, M. and Ellahi, R. Numerical and analytical solutions of an Oldroyd 8-constant MHD fluid with nonlinear slip conditions. International Journal for Numerical Methods in Fluids, 67(10), 1234-1246(2011)
[29] Hameed, M. and Ellahi, R. Thin film flow of non-Newtonian MHD fluid on a vertically moving belt. International Journal for Numerical Methods in Fluids, 66(11), 1409-1419(2011)
[30] Ellahi, R. Numerical analysis of steady non-Newtonian flows with heat transfer analysis, MHD and nonlinear slip effects. International Journal of Numerical Methods for Heat & Fluid Flow, 22(1), 24-28(2012)
[31] Ellahi, R. The effects of MHD and temperature dependent viscosity on the flow of non-Newtonian nanofluid in a pipe:analytical solutions. Applied Mathematical Modelling, 37(3), 1451-1467(2013)
[32] Luo, X. H., Li, B. W., Zhang, J. K., and Hu, Z. M. Simulation of thermal radiation effects on MHD free convection in a square cavity using Chebyshev collocation spectral method. Numerical Heat Transfer, Part A, 66(7), 792-815(2014)
[33] Hugues, S. and Randriamampianina, A. An improved projection scheme applied to pseudospectral methods for the incompressible Navier-Stokes equations. International Journal for Numerical Methods in Fluids, 28(3), 501-521(1998)
[34] Tuckerman, L. S. Divergence-free velocity fields in nonperiodic geometries. Journal of Computational Physics, 80(2), 403-441(1989)
[35] Chen, H. L., Su, Y. H., and Shizgal, B. D. A direct spectral collocation Poisson solver in polar and cylindrical coordinates. Journal of Computational Physics, 160(2), 453-469(2000)
[36] Johnston, H. and Liu, J. G. Accurate, stable and efficient Navier-Stokes solvers based on explicit treatment of the pressure term. Journal of Computational Physics, 199(1), 221-259(2004)
[37] Burggraf, O. R. Analytical and numerical studies of the structure of steady separated flows. Journal of Fluid Mechanics, 24(1), 113-151(1966)
[38] Ghia, U., Ghia, K. N., and Shin, C. T. High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method. Journal of Computational Physics, 48(3), 387-411(1982)
[39] Wahba, E. M. Steady flow simulations inside a driven cavity up to Reynolds number 35000. Computers & Fluids, 66, 85-97(2012)
[40] Erturk, E., Corke, T. C., and Gokcol, C. Numerical solutions of 2D steady incompressible driven cavity flow at high Reynolds numbers. International Journal for Numerical Methods in Fluids, 48(7), 747-774(2005)
[41] Moffat, H. K. Viscous and resistive eddies near a sharp corner. Journal of Fluid Mechanics, 18(1), 1-18(1964)
[42] Ehrenstein, U. and Peyret, R. A Chebyshev collocation method for the Navier-Stokes equations with application to double-diffusive convection. International Journal for Numerical Methods in Fluids, 9(4), 427-452(1989)

Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

Direct spectral domain decomposition method for 2D incompressible Navier-Stokes equations

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 5

编辑推荐

Metrics

本文评价

[1]	Chao XU, Dongyang SHI, Xin LIAO. Low order nonconforming mixed finite element method for nonstationary incompressible Navier-Stokes equations[J]. Applied Mathematics and Mechanics (English Edition), 2016, 37(8): 1095-1112.
[2]	陈光华陈光明戴智华. Modified domain decomposition method for Hamilton-Jacobi-Bellman equations[J]. Applied Mathematics and Mechanics (English Edition), 2010, 31(12): 1585-1592.
[3]	尚月强罗振东. A parallel two-level finite element method for the Navier-Stokes equations[J]. Applied Mathematics and Mechanics (English Edition), 2010, 31(11): 1429-1438.
[4]	尚月强何银年. Fourier analysis of Schwarz domain decomposition methods for the biharmonic equation[J]. Applied Mathematics and Mechanics (English Edition), 2009, 30(9): 1177-1182.
[5]	温功碧;陈作斌. UNSTEADY/STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY[J]. Applied Mathematics and Mechanics (English Edition), 2004, 25(1): 59-72.