A Meshfree Generalized Finite Difference Method for Surface PDEs

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

On Meshfree GFDM Solvers for the Incompressible Navier-Stokes Equations

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

A Meshfree Lagrangian Method for Flow on Manifolds

In this paper, we present a novel meshfree framework for fluid flow simulations on arbitrarily curved surfaces. First, we introduce a new meshfree Lagrangian framework to model flow on surfaces. Meshfree points or particles, which are used to discretize the domain, move in a Lagrangian sense along the given surface. This is done without discretizing the bulk around the surface, without parametrizing the surface, and without a background mesh. A key novelty that is introduced is the handling of flow with evolving free boundaries on a curved surface. The use of this framework to model flow on moving and deforming surfaces is also introduced. Then, we present the application of this framework to solve fluid flow problems defined on surfaces numerically. In combination with a meshfree Generalized Finite Difference Method (GFDM), we introduce a strong form meshfree collocation scheme to solve the Navier-Stokes equations posed on manifolds. Benchmark examples are proposed to validate the Lagrangian framework and the surface Navier-Stokes equations with the presence of free boundaries

Particle-based adaptive coupling of 3D and 2D fluid flow models

peer reviewedThis paper proposes the notion of model adaptivity for fluid flow modelling, where the under- lying model (the governing equations) is adaptively changed in space and time. Specifically, this work introduces a hybrid and adaptive coupling of a 3D bulk fluid flow model with a 2D thin film flow model. As a result, this work extends the applicability of existing thin film flow models to complex scenarios where, for example, bulk flow develops into thin films after striking a surface. At each location in space and time, the proposed framework automatically decides whether a 3D model or a 2D model must be applied. Using a meshless approach for both 3D and 2D models, at each particle, the decision to apply a 2D or 3D model is based on the user-prescribed resolution and a local principal component analysis. When a particle needs to be changed from a 3D model to 2D, or vice versa, the discretization is changed, and all relevant data mapping is done on-the-fly. Appropriate two-way coupling conditions and mass conservation considerations between the 3D and 2D models are also developed. Numerical results show that this model adaptive framework shows higher flexibility and compares well against finely resolved 3D simulations. In an actual application scenario, a 3 factor speed up is obtained, while maintaining the accuracy of the solution.U-AGR-6054 - IAS-AUDACITY ADONIS - BORDAS Stéphan

Higher-Order GFDM for Linear Elliptic Operators

We present a novel approach of discretizing diffusion operators of the form $\nabla\cdot(\lambda\nabla u)$ in the context of meshfree generalized finite difference methods. Our ansatz uses properties of derived operators and combines the discrete Laplace operator with reconstruction functions approximating the diffusion coefficient $\lambda$. Provided that the reconstructions are of a sufficiently high order, we prove that the order of accuracy of the discrete Laplace operator transfers to the derived diffusion operator. We show that the new discrete diffusion operator inherits the diagonal dominance property of the discrete Laplace operator and fulfills enrichment properties. Our numerical results for elliptic and parabolic partial differential equations show that even low-order reconstructions preserve the order of the underlying discrete Laplace operator for sufficiently smooth diffusion coefficients. In experiments, we demonstrate the applicability of the new discrete diffusion operator to interface problems with point clouds not aligning to the interface and numerically prove first-order convergence

A Meshfree Generalized Finite Difference Method for Solution Mining Processes

Experimental and field investigations for solution mining processes have improved intensely in recent years. Due to today's computing capacities, three-dimensional simulations of potential salt solution caverns can further enhance the understanding of these processes. They serve as a "virtual prototype" of a projected site and support planning in reasonable time. In this contribution, we present a meshfree Generalized Finite Difference Method (GFDM) based on a cloud of numerical points that is able to simulate solution mining processes on microscopic as well as macroscopic scales, which differ significantly in both the spatial and temporal scale. Focusing on anticipated industrial requirements, Lagrangian and Eulerian formulations including an Arbitrary Lagrangian-Eulerian (ALE) approach are considered

A Lagrangian meshfree model for solidification of liquid thin-films

peer reviewedIn this paper, a new method to model solidification of thin liquid films is proposed. This method is targeted at applications like aircraft icing and tablet coating where the formation of liquid films from impinging droplets on a surface form a critical part of the physics of the process. The proposed model takes into account the (i) unsteadiness in temperature distribution, (ii) heat transfer at the interface between the solid and the surface, (iii) volumetric expansion/contraction and (iv) the liquid thin-film behaviour, each of which are either partly or fully ignored in existing models. The liquid thin-film, modelled using the Discrete Droplet Method (DDM), is represented as a collection of discrete droplets that are tracked in a Lagrangian sense. The height of the liquid film is estimated as a summation of Gaussian kernel functions associated with each droplet. At each droplet location, a solid height is also computed. The evolution of the solid height is governed by the Stefan problem. The flow of the liquid thin-film is solved just as in the case of DDM, while also taking into consideration the shape of the solidified region lying beneath the droplet. The results presented in this work show the reliability of the proposed model in simulating solidification of thin-films and its applicability to complex problems such as ice-formation on aircraft wings. The model has been verified for canonical problems that have analytical solutions. For the more complex problems of icing, the results of the model are compared with data from literature, without considering a background air flow. The comparison can be improved by coupling this model with suitable air flow solvers, as shown in the final test case.U-AGR-6054 - IAS-AUDACITY ADONIS - BORDAS Stéphan