Consistency and Convergence of a High Order Accurate Meshless Method for Solution of Incompressible Fluid Flows

by   Shantanu Shahane, et al.

Computations of incompressible flows with velocity boundary conditions require solution of a Poisson equation for pressure with all Neumann boundary conditions. Discretization of such a Poisson equation results in a rank-deficient matrix of coefficients. When a non-conservative discretization method such as finite difference, finite element, or spectral scheme is used, such a matrix also generates an inconsistency which makes the residuals in the iterative solution to saturate at a threshold level that depends on the spatial resolution and order of the discretization scheme. In this paper, we examine inconsistency for a high-order meshless discretization scheme suitable for solving the equations on a complex domain. The high order meshless method uses polyharmonic spline radial basis functions (PHS-RBF) with appended polynomials to interpolate scattered data and constructs the discrete equations by collocation. The PHS-RBF provides the flexibility to vary the order of discretization by increasing the degree of the appended polynomial. In this study, we examine the convergence of the inconsistency for different spatial resolutions and for different degrees of the appended polynomials by solving the Poisson equation for a manufactured solution as well as the Navier-Stokes equations for several fluid flows. We observe that the inconsistency decreases faster than the error in the final solution, and eventually becomes vanishing small at sufficient spatial resolution. The rate of convergence of the inconsistency is observed to be similar or better than the rate of convergence of the discretization errors. This beneficial observation makes it unnecessary to regularize the Poisson equation by fixing either the mean pressure or pressure at an arbitrary point. A simple point solver such as the SOR is seen to be well-convergent, although it can be further accelerated using multilevel methods.


page 36

page 40


A Non-Nested Multilevel Method for Meshless Solution of the Poisson Equation in Heat Transfer and Fluid Flow

We present a non-nested multilevel algorithm for solving the Poisson equ...

Parallel implementation of a compatible high-order meshless method for the Stokes' equations

A parallel implementation of a compatible discretization scheme for stea...

High-order Methods for a Pressure Poisson Equation Reformulation of the Navier-Stokes Equations with Electric Boundary Conditions

Pressure Poisson equation (PPE) reformulations of the incompressible Nav...

A Semi-Implicit Meshless Method for Incompressible Flows in Complex Geometries

We present an exponentially convergent semi-implicit meshless algorithm ...

On the efficient preconditioning of the Stokes equations in tight geometries

If the Stokes equations are properly discretized, it is well-known that ...

A High-Order Accurate Meshless Method for Solution of Incompressible Fluid Flow Problems

Meshless solution to differential equations using radial basis functions...

A time-spectral Stokes solver for simulation of time-periodic flows in complex geometries

Simulation of unsteady creeping flows in complex geometries has traditio...

Please sign up or login with your details

Forgot password? Click here to reset