A novel implicit scheme based on free element method (FECM) is proposed to solve Navier–Stokes equations. The proposed method has the mesh-free features but elements are essential instead of point clouds. The isoparametric transformations, which are usually used in finite element method, are employed to obtain the coefficients in the spatial derivative approximation. The method satisfies the geometric conservation law and first-order consistency. An approximate Riemann solver is employed to compute convective fluxes at the midpoints between collocated points and field points. Meanwhile, a reconstruction approach is used to reach higher order discretization. Viscous terms are approximated with a corrected average progress of gradients. Furthermore, a matrix-free Lower-Upper Symmetric Gauss–Seidel (LU-SGS) scheme is introduced to speed up the convergence. In order to validate the proposed method, three test cases varying from subsonic to supersonic flows are simulated and discussed in detail. The test cases validate that the proposed method is as accurate as second order finite volume solver and is robust with different node sets. What's more, LU-SGS algorithm can significantly improve convergence speed compared with the explicit scheme.
Read full abstract