This paper presents a new version of the upwind compact finite difference scheme for solving the incompressible Navier-Stokes equations in generalized curvilinear coordinates.The artificial compressibility approach is...This paper presents a new version of the upwind compact finite difference scheme for solving the incompressible Navier-Stokes equations in generalized curvilinear coordinates.The artificial compressibility approach is used,which transforms the elliptic-parabolic equations into the hyperbolic-parabolic ones so that flux difference splitting can be applied.The convective terms are approximated by a third-order upwind compact scheme implemented with flux difference splitting,and the viscous terms are approximated by a fourth-order central compact scheme.The solution algorithm used is the Beam-Warming approximate factorization scheme.Numerical solutions to benchmark problems of the steady plane Couette-Poiseuille flow,the liddriven cavity flow,and the constricting channel flow with varying geometry are presented.The computed results are found in good agreement with established analytical and numerical results.The third-order accuracy of the scheme is verified on uniform rectangular meshes.展开更多
基金This work was supported by Natural Science Foundation of China(G10476032,G10531080)state key program for developing basic sciences(2005CB321703).
文摘This paper presents a new version of the upwind compact finite difference scheme for solving the incompressible Navier-Stokes equations in generalized curvilinear coordinates.The artificial compressibility approach is used,which transforms the elliptic-parabolic equations into the hyperbolic-parabolic ones so that flux difference splitting can be applied.The convective terms are approximated by a third-order upwind compact scheme implemented with flux difference splitting,and the viscous terms are approximated by a fourth-order central compact scheme.The solution algorithm used is the Beam-Warming approximate factorization scheme.Numerical solutions to benchmark problems of the steady plane Couette-Poiseuille flow,the liddriven cavity flow,and the constricting channel flow with varying geometry are presented.The computed results are found in good agreement with established analytical and numerical results.The third-order accuracy of the scheme is verified on uniform rectangular meshes.