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

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.

[1]  Joerg Kuhnert,et al.  A meshfree generalized finite difference method for surface PDEs , 2018, Comput. Math. Appl..

[2]  Bengt Fornberg,et al.  Fast variable density 3-D node generation , 2019, SIAM J. Sci. Comput..

[3]  Varun Shankar,et al.  Hyperviscosity-Based Stabilization for Radial Basis Function-Finite Difference (RBF-FD) Discretizations of Advection-Diffusion Equations , 2018, J. Comput. Phys..

[4]  Steven J. Ruuth,et al.  An RBF-FD closest point method for solving PDEs on surfaces , 2018, J. Comput. Phys..

[5]  Bengt Fornberg,et al.  Fast generation of 2-D node distributions for mesh-free PDE discretizations , 2015, Comput. Math. Appl..

[6]  Mario Mulansky,et al.  Odeint - Solving ordinary differential equations in C++ , 2011, ArXiv.

[7]  Gregor Kosec,et al.  A local numerical solution of a fluid-flow problem on an irregular domain , 2016, Adv. Eng. Softw..

[8]  Frédéric Hecht,et al.  New development in freefem++ , 2012, J. Num. Math..

[9]  Holger Wendland,et al.  Scattered Data Approximation: Conditionally positive definite functions , 2004 .

[10]  Bakhtier Farouk,et al.  A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure , 1991 .

[11]  W. Benz Smooth Particle Hydrodynamics: A Review , 1990 .

[12]  P. Draper,et al.  SWIFT: SPH With Inter-dependent Fine-grained Tasking , 2018 .

[13]  Gregor Kosec,et al.  Solution of thermo‐fluid problems by collocation with local pressure correction , 2008 .

[14]  Bengt Fornberg,et al.  On the role of polynomials in RBF-FD approximations: I. Interpolation and accuracy , 2016, J. Comput. Phys..

[15]  Andrew T. T. McRae,et al.  Firedrake: automating the finite element method by composing abstractions , 2015, ACM Trans. Math. Softw..

[16]  Bengt Fornberg,et al.  On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs , 2017, J. Comput. Phys..

[17]  Hua Li,et al.  Meshless Methods and Their Numerical Properties , 2013 .

[18]  Jure Slak,et al.  On generation of node distributions for meshless PDE discretizations , 2018, SIAM J. Sci. Comput..

[19]  Konstantinos Poulios,et al.  GetFEM , 2020, ACM Trans. Math. Softw..

[20]  H. Sadat,et al.  PERFORMANCE AND ACCURACY OF A MESHLESS METHOD FOR LAMINAR NATURAL CONVECTION , 2000 .

[21]  H. Munthe-Kaas,et al.  Coordinate free numerics : Closing the gap between 'pure' and 'applied' mathematics ? , 1996 .

[22]  Stephen M. Longshaw,et al.  DualSPHysics: Open-source parallel CFD solver based on Smoothed Particle Hydrodynamics (SPH) , 2015, Comput. Phys. Commun..

[23]  Jure Slak,et al.  Parallel coordinate free implementation of local meshless method , 2018, 2018 41st International Convention on Information and Communication Technology, Electronics and Microelectronics (MIPRO).

[24]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[25]  Anders Logg,et al.  DOLFIN: Automated finite element computing , 2010, TOMS.

[26]  Zhaoli Guo,et al.  Numerical study of three-dimensional natural convection in a cubical cavity at high Rayleigh numbers , 2017 .

[27]  G. Kosec,et al.  Cooling of overhead power lines due to the natural convection , 2019 .

[28]  W. Bangerth,et al.  deal.II—A general-purpose object-oriented finite element library , 2007, TOMS.

[29]  Hamou Sadat,et al.  A new meshless approach for three dimensional fluid flow and related heat transfer problems , 2012, 1809.02502.

[30]  A. I. Tolstykh,et al.  On using radial basis functions in a “finite difference mode” with applications to elasticity problems , 2003 .

[31]  Alessandro Veneziani,et al.  The LifeV library: engineering mathematics beyond the proof of concept , 2017, 1710.06596.

[32]  G. Kosec,et al.  Refined Meshless Local Strong Form solution of Cauchy–Navier equation on an irregular domain , 2019, Engineering Analysis with Boundary Elements.

[34]  J. Kuhnert,et al.  Finite Pointset Method Based on the Projection Method for Simulations of the Incompressible Navier-Stokes Equations , 2003 .

[35]  Eugenio Oñate,et al.  A finite point method for elasticity problems , 2001 .

[36]  Marc Duflot,et al.  Meshless methods: A review and computer implementation aspects , 2008, Math. Comput. Simul..

[37]  Sudarshan Tiwari,et al.  Finite pointset method for simulation of the liquid-liquid flow field in an extractor , 2008, Comput. Chem. Eng..

[38]  Martin Robinson,et al.  Particle-based and Meshless Methods with Aboria , 2017, SoftwareX.

[39]  Luis Gavete,et al.  Improvements of generalized finite difference method and comparison with other meshless method , 2003 .

[40]  Manuel Kindelan,et al.  RBF-FD formulas and convergence properties , 2010, J. Comput. Phys..

[41]  David Levin,et al.  The approximation power of moving least-squares , 1998, Math. Comput..

[42]  Oleg Davydov,et al.  Adaptive RBF-FD method for elliptic problems with point singularities in 2D , 2016, Appl. Math. Comput..

[43]  Rainald Löhner,et al.  A general advancing front technique for filling space with arbitrary objects , 2004 .

[44]  R. Trobec,et al.  Weak and strong from meshless methods for linear elastic problem under fretting contact conditions , 2019, Tribology International.

[45]  Yo-Ming Hsieh,et al.  ESFM: An Essential Software Framework for Meshfree Methods , 2014, Adv. Eng. Softw..

[46]  J. Robert Buchler,et al.  The Numerical Modelling of Nonlinear Stellar Pulsations , 1990 .

[47]  Natasha Flyer,et al.  A 3-D RBF-FD solver for modeling the atmospheric global electric circuit with topography (GEC-RBFFD v1.0) , 2015 .

[48]  Gerd Heber,et al.  An overview of the HDF5 technology suite and its applications , 2011, AD '11.

[49]  Guirong Liu Mesh Free Methods: Moving Beyond the Finite Element Method , 2002 .

[50]  Arpit Agarwal,et al.  PySPH: A Python-based Framework for Smoothed Particle Hydrodynamics , 2019, ACM Trans. Math. Softw..

[51]  Riccardo Zamolo,et al.  Two algorithms for fast 2D node generation: Application to RBF meshless discretization of diffusion problems and image halftoning , 2018, Comput. Math. Appl..

[52]  Slawomir Milewski Selected computational aspects of the meshless finite difference method , 2012, Numerical Algorithms.

[53]  Ying Liu,et al.  Node Placement Method by Bubble Simulation and Its Application , 2010 .

[54]  Robert Michael Kirby,et al.  Robust Node Generation for Meshfree Discretizations on Irregular Domains and Surfaces , 2018, SIAM J. Sci. Comput..