Medusa: A C++ Library for solving PDEs using Strong Form Mesh-Free methods

Abstract

Medusa, a novel library for implementation of strong form mesh-free methods, is described. We identify and present common parts and patterns among many such methods reported in the literature, such as node positioning, stencil selection and stencil weight computation. Many different algorithms exist for each part and the possible combinations offer a plethora of possibilities for improvements of solution procedures that are far from fully understood. As a consequence there are still many unanswered questions in mesh-free community resulting in vivid ongoing research in the field. Medusa implements the core mesh-free elements as independent blocks, which offers users great flexibility in experimenting with the method they are developing, as well as easily comparing it with other existing methods. The paper describes the chosen abstractions and their usage, illustrates aspects of the philosophy and design, offers some executions time benchmarks and demonstrates the application of the library on cases from linear elasticity and fluid flow in irregular 2D and 3D domains.

Cited by in corpus (14)

• Some observations regarding the RBF-FD approximation accuracy dependence on stencil size
• RBF-FD analysis of 2D time-domain acoustic wave propagation in heterogeneous media
• Dynamic mode decomposition as an analysis tool for time-dependent partial differential equations
• Discretization of non-uniform rational B-spline (NURBS) models for meshless isogeometric analysis
• Meshless method stencil evaluation with machine learning
• Refined radial basis function-generated finite difference analysis of non-Newtonian natural convection
• Stability analysis of RBF-FD and WLS based local strong form meshless methods on scattered nodes
• GPU accelerated RBF-FD solution of Poisson's equation
• A sharp-interface mesoscopic model for dendritic growth
• An RBF-based method for computational electromagnetics with reduced numerical dispersion
• Performance Trade-offs of High Order Meshless Approximation on Distributed Memory Systems
• Spatially-varying meshless approximation method for enhanced computational efficiency
• A superconvergence result in the RBF-FD method
• Adaptive hyperviscosity stabilisation for the RBF-FD method in solving advection-dominated transport equations