High order difference schemes using the local anisotropic basis function method

DOI10.1016/J.JCP.2020.109549MaRDI QIDQ2194355zbMATHOpenAlexFDO

Authors J. R. C. King, S. J. Lind, A. M. A. Nasar

Publication date 25 August 2020

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1912.06441

zbMATH Keywords

consistencyhigh orderradial basis functionsmoothed particle hydrodynamicsmesh-free

Mathematics Subject Classification ID

Probabilistic methods, particle methods, etc. for boundary value problems involving PDEs (65N75) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Probabilistic methods, particle methods, etc. for initial value and initial-boundary value problems involving PDEs (65M75) Particle methods and lattice-gas methods (76M28)

Abstract: Mesh-free methods have significant potential for simulations in complex geometries, as the time consuming process of mesh-generation is avoided. Smoothed Particle Hydrodynamics (SPH) is the most widely used mesh-free method, but suffers from a lack of consistency. High order, consistent, and local (using compact computational stencils) mesh-free methods are particularly desirable. Here we present a novel framework for generating local high order difference operators for arbitrary node distributions, referred to as the Local Anisotropic Basis Function Method (LABFM). Weights are constructed from linear sums of anisotropic basis functions (ABFs), chosen to ensure exact reproduction of polynomial fields up to a given order. The ABFs are based on a fundamental Radial Basis Function (RBF), and the choice of fundamental RBF has small effect on accuracy, but influences stability. LABFM is able to generate high order difference operators with compact computational stencils (4th order with 25 nodes, 8th order with 60 nodes in two dimensions). At domain boundaries (with incomplete support) LABFM automatically provides one-sided differences of the same order as the internal scheme, up to 4th order. We use the method to solve elliptic, parabolic and mixed hyperbolic-parabolic PDEs, showing up to 8th order convergence. The inclusion of hyperviscosity is straightforward, and can effectively provide stability when solving hyperbolic problems. LABFM is a promising new mesh-free method for the numerical solution of PDEs in complex geometries. The method is highly scalable, and for Eulerian schemes, the computational efficiency is competitive with RBF-FD for a given accuracy. A particularly attractive feature is that in the low order limit, LABFM collapses to SPH, and there is potential for Arbitrary Lagrangian-Eulerian schemes with natural adaptivity of resolution and accuracy.

Recommendations

• Solving PDEs with radial basis functions
• A high order meshless method with compact support
• Construction and analysis of meshless finite difference methods
• Minisymposium “Meshfree Methods for the Solution of PDEs”
• Enhancing SPH using moving least-squares and radial basis functions

Cites work

• scientific article; zbMATH DE number 3640828 (Why is no real title available?)
• scientific article; zbMATH DE number 1396373 (Why is no real title available?)
• A compatible high-order meshless method for the Stokes equations with applications to suspension flows
• A conservative mesh-free scheme and generalized framework for conservation laws
• A conservative, consistent, and scalable meshfree mimetic method
• A generalized smoothed particle hydrodynamics method for nonlinear dynamic problems
• A high-order staggered meshless method for elliptic problems
• A stabilized radial basis-finite difference (RBF-FD) method with hybrid kernels
• An Eulerian-Lagrangian incompressible SPH formulation (ELI-SPH) connected with a sharp interface
• An algorithm to improve consistency in smoothed particle hydrodynamics
• Compact finite difference schemes with spectral-like resolution
• Enhancing finite differences with radial basis functions: experiments on the Navier-Stokes equations
• High-order Eulerian incompressible smoothed particle hydrodynamics with transition to Lagrangian free-surface motion
• Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves
• Local radial basis function-based differential quadrature method and its application to solve two-dimensional incompressible Navier--Stokes equations
• Meshfree Methods: A Comprehensive Review of Applications
• Modeling incompressible flows using a finite particle method
• Modified smoothed particle hydrodynamics method and its application to transient problems
• Multiquadrics - a scattered data approximation scheme with applications to computational fluid-dynamics. I: Surface approximations and partial derivative estimates
• Multiquadrics -- a scattered data approximation scheme with applications to computational fluid-dynamics. II: Solutions to parabolic, hyperbolic and elliptic partial differential equations
• Novel finite particle formulations based on projection methodologies
• Numerical methods for high dimensional Hamilton-Jacobi equations using radial basis functions
• On a quasi-linear parabolic equation occurring in aerodynamics
• On stability and convergence of projection methods based on pressure Poisson equation
• On the role of polynomials in RBF-FD approximations. I: Interpolation and accuracy
• On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs
• Particle Methods for a 1D Elastic Model Problem: Error Analysis and Development of a Second-Order Accurate Formulation
• Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree
• Reproducing kernel particle methods
• Reproducing kernel particle methods for structural dynamics
• Restoring particle consistency in smoothed particle hydrodynamics
• Smoothed particle hydrodynamics and its diverse applications
• Smoothed particle hydrodynamics: theory and application to non-spherical stars
• Solving PDEs with radial basis functions
• Solving parabolic and hyperbolic equations by the generalized finite difference method
• Solving second order non-linear elliptic partial differential equations using generalized finite difference method
• Stabilization of RBF-generated finite difference methods for convective PDEs
• The exponentially convergent trapezoidal rule
• Truncation error in mesh-free particle methods
• Variational and momentum preservation aspects of smooth particle hydrodynamic formulations

Cited in

(6)

• High-order simulations of isothermal flows using the local anisotropic basis function method (LABFM)
• High-order consistent SPH with the pressure projection method in 2-D and 3-D
• Coupling of an SPH-based solver with a multiphysics library
• Towards high-order consistency and convergence of conservative SPH approximations
• An alternative SPH formulation: ADER-WENO-SPH
• A mesh-free framework for high-order direct numerical simulations of combustion in complex geometries

This page was built for publication: High order difference schemes using the local anisotropic basis function method

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2194355)