An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. ...An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangu- lar 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 Neu- mann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet bound- ary 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 be- tween the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.展开更多
Performance comparisons are composed of two parts: the first part contains the systematically investigation of six difference schemes including CDS, FUDS, HDS, PLDS, SUDS and QUICK for convection terms in numerical fl...Performance comparisons are composed of two parts: the first part contains the systematically investigation of six difference schemes including CDS, FUDS, HDS, PLDS, SUDS and QUICK for convection terms in numerical fluid flow and heat transfer based on the finite volume method using staggered and Rhie-Chow’s momentum interpolation collocated grids, the second part contains the comparative computations being conducted on Rhie-Chow’s momentum interpolation collocated grid and Thiart’s finite difference scheme based nonstaggered grid. Three 3-D cases that have analytical or benchmark solutions are adopted. For the first part, the results of computations indicate that, all the six schemes have the same numerical accuracy when the diffusion term is predominant. With the increase of convection, the FUDS, HDS and PLDS almost have the same accuracy in two of those grid systems, while the SUDS and QUICK have higher accuracy than the former. The accuracy of CDS is something in between. For the same under-relaxation factors and convergence criterion, the convergence rate of each scheme on those two grid systems are nearly equal with that on the staggered grid being a little bit faster. For QUICK and CDS, smooth, non-oscillating solutions can be obtained even when local Peclet number may be as large as 31.2-31.3. For the second part, it is concluded that simplified collocated grid system is preferable from numerical accuracy, grid Peclet number limit, sensitivity to the underrelaxation factor and the freedom in choosing finite difference scheme for convection term.展开更多
基金Project supported by the National Natural Science Foundation of China(No.51176026)the Fundamental Research Funds for the Central Universities(No.DUT14RC(3)029)
文摘An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangu- lar 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 Neu- mann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet bound- ary 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 be- tween the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.
文摘Performance comparisons are composed of two parts: the first part contains the systematically investigation of six difference schemes including CDS, FUDS, HDS, PLDS, SUDS and QUICK for convection terms in numerical fluid flow and heat transfer based on the finite volume method using staggered and Rhie-Chow’s momentum interpolation collocated grids, the second part contains the comparative computations being conducted on Rhie-Chow’s momentum interpolation collocated grid and Thiart’s finite difference scheme based nonstaggered grid. Three 3-D cases that have analytical or benchmark solutions are adopted. For the first part, the results of computations indicate that, all the six schemes have the same numerical accuracy when the diffusion term is predominant. With the increase of convection, the FUDS, HDS and PLDS almost have the same accuracy in two of those grid systems, while the SUDS and QUICK have higher accuracy than the former. The accuracy of CDS is something in between. For the same under-relaxation factors and convergence criterion, the convergence rate of each scheme on those two grid systems are nearly equal with that on the staggered grid being a little bit faster. For QUICK and CDS, smooth, non-oscillating solutions can be obtained even when local Peclet number may be as large as 31.2-31.3. For the second part, it is concluded that simplified collocated grid system is preferable from numerical accuracy, grid Peclet number limit, sensitivity to the underrelaxation factor and the freedom in choosing finite difference scheme for convection term.