NAJan 31, 2017
A Flux Conserving Meshfree Method for Conservation LawsPratik Suchde, Joerg Kuhnert, Simon Schroeder et al.
Lack of conservation has been the biggest drawback in meshfree generalized finite difference methods (GFDMs). In this paper, we present a novel modification of classical meshfree GFDMs to include local balances which produce an approximate conservation of numerical fluxes. This numerical flux conservation is done within the usual moving least squares framework. Unlike Finite Volume Methods, it is based on locally defined control cells, rather than a globally defined mesh. We present the application of this method to an advection diffusion equation and the incompressible Navier - Stokes equations. Our simulations show that the introduction of flux conservation significantly reduces the errors in conservation in meshfree GFDMs.
NAJan 12, 2017
On Meshfree GFDM Solvers for the Incompressible Navier-Stokes EquationsPratik Suchde, Joerg Kuhnert, Sudarshan Tiwari
Meshfree solution schemes for the incompressible Navier--Stokes equations are usually based on algorithms commonly used in finite volume methods, such as projection methods, SIMPLE and PISO algorithms. However, drawbacks of these algorithms that are specific to meshfree methods have often been overlooked. In this paper, we study the drawbacks of conventionally used meshfree Generalized Finite Difference Method~(GFDM) schemes for Lagrangian incompressible Navier-Stokes equations, both operator splitting schemes and monolithic schemes. The major drawback of most of these schemes is inaccurate local approximations to the mass conservation condition. Further, we propose a new modification of a commonly used monolithic scheme that overcomes these problems and shows a better approximation for the velocity divergence condition. We then perform a numerical comparison which shows the new monolithic scheme to be more accurate than existing schemes.
NAApr 3, 2017
Point Cloud Movement For Fully Lagrangian Meshfree MethodsPratik Suchde, Joerg Kuhnert
In Lagrangian meshfree methods, the underlying spatial discretization, referred to as a point cloud or a particle cloud, moves with the flow velocity. In this paper, we consider different numerical methods of performing this movement of points or particles. The movement is most commonly done by a first order method, which assumes the velocity to be constant within a time step. We show that this method is very inaccurate and that it introduces volume and mass conservation errors. We further propose new methods for the same which prescribe an additional ODE system that describes the characteristic velocity. Movement is then performed along this characteristic velocity. The first new way of moving points is an extension of mesh-based streamline tracing ideas to meshfree methods. In the second way, the movement is done based on the difference in approximated streamlines between two time levels, which approximates the pathlines in unsteady flow. Numerical comparisons show these method to be vastly superior to the conventionally used first order method.
NANov 15, 2018
A Meshfree Generalized Finite Difference Method for Surface PDEsPratik Suchde, Joerg Kuhnert
In this paper, we propose a novel meshfree Generalized Finite Difference Method (GFDM) approach to discretize PDEs defined on manifolds. Derivative approximations for the same are done directly on the tangent space, in a manner that mimics the procedure followed in volume-based meshfree GFDMs. As a result, the proposed method not only does not require a mesh, it also does not require an explicit reconstruction of the manifold. In contrast to existing methods, it avoids the complexities of dealing with a manifold metric, while also avoiding the need to solve a PDE in the embedding space. A major advantage of this method is that all developments in usual volume-based numerical methods can be directly ported over to surfaces using this framework. We propose discretizations of the surface gradient operator, the surface Laplacian and surface Diffusion operators. Possibilities to deal with anisotropic and discontinous surface properties (with large jumps) are also introduced, and a few practical applications are presented.