High order difference schemes using the local anisotropic basis function method
From MaRDI portal
Publication:2194355
DOI10.1016/J.JCP.2020.109549zbMATH Open1440.76113arXiv1912.06441OpenAlexW3025860596MaRDI QIDQ2194355FDOQ2194355
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)
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.
Full work available at URL: https://arxiv.org/abs/1912.06441
Recommendations
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)
Cites Work
- Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree
- On a quasi-linear parabolic equation occurring in aerodynamics
- Numerical methods for high dimensional Hamilton-Jacobi equations using radial basis functions
- Title not available (Why is that?)
- Compact finite difference schemes with spectral-like resolution
- Multiquadrics -- a scattered data approximation scheme with applications to computational fluid-dynamics. II: Solutions to parabolic, hyperbolic and elliptic partial differential equations
- Smoothed particle hydrodynamics: theory and application to non-spherical stars
- Reproducing kernel particle methods
- Multiquadrics - a scattered data approximation scheme with applications to computational fluid-dynamics. I: Surface approximations and partial derivative estimates
- Local radial basis function-based differential quadrature method and its application to solve two-dimensional incompressible Navier--Stokes equations
- The exponentially convergent trapezoidal rule
- Stabilization of RBF-generated finite difference methods for convective PDEs
- Restoring particle consistency in smoothed particle hydrodynamics
- Smoothed particle hydrodynamics and its diverse applications
- Reproducing kernel particle methods for structural dynamics
- A generalized smoothed particle hydrodynamics method for nonlinear dynamic problems
- Variational and momentum preservation aspects of smooth particle hydrodynamic formulations
- Solving parabolic and hyperbolic equations by the generalized finite difference method
- Modeling incompressible flows using a finite particle method
- Truncation error in mesh-free particle methods
- Title not available (Why is that?)
- Solving PDEs with radial basis functions
- Modified smoothed particle hydrodynamics method and its application to transient problems
- On stability and convergence of projection methods based on pressure Poisson equation
- Incompressible smoothed particle hydrodynamics for free-surface flows: A generalised diffusion-based algorithm for stability and validations for impulsive flows and propagating waves
- A conservative mesh-free scheme and generalized framework for conservation laws
- Solving second order non-linear elliptic partial differential equations using generalized finite difference method
- Novel finite particle formulations based on projection methodologies
- An algorithm to improve consistency in smoothed particle hydrodynamics
- Enhancing finite differences with radial basis functions: experiments on the Navier-Stokes equations
- On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs
- On the role of polynomials in RBF-FD approximations. I: Interpolation and accuracy
- A compatible high-order meshless method for the Stokes equations with applications to suspension flows
- A high-order staggered meshless method for elliptic problems
- A conservative, consistent, and scalable meshfree mimetic method
- A stabilized radial basis-finite difference (RBF-FD) method with hybrid kernels
- Meshfree Methods: A Comprehensive Review of Applications
- Particle Methods for a 1D Elastic Model Problem: Error Analysis and Development of a Second-Order Accurate Formulation
- High-order Eulerian incompressible smoothed particle hydrodynamics with transition to Lagrangian free-surface motion
- An Eulerian-Lagrangian incompressible SPH formulation (ELI-SPH) connected with a sharp interface
Cited In (6)
- High-order simulations of isothermal flows using the local anisotropic basis function method (LABFM)
- An alternative SPH formulation: ADER-WENO-SPH
- High-order consistent SPH with the pressure projection method in 2-D and 3-D
- Towards high-order consistency and convergence of conservative SPH approximations
- Coupling of an SPH-based solver with a multiphysics library
- 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)
