A stable SPH with adaptive B-spline kernel

01/04/2020
by   Saptarshi Kumar Lahiri, et al.
0

Tensile instability, often observed in smoothed particle hydrodynamics (SPH), is a numerical artifact that manifests itself by unphysical clustering or separation of particles. The instability originates in estimating the derivatives of the smoothing functions which, when interact with material constitution may result in negative stiffness in the discretized system. In the present study, a stable formulation of SPH is developed where the kernel function is continuously adapted at every material point depending on its state of stress. Bspline basis function with a variable intermediate knot is used as the kernel function. The shape of the kernel function is then modified by changing the intermediate knot position such that the condition associated with instability does not arise. While implementing the algorithm the simplicity and computational efficiency of SPH are not compromised. One-dimensional dispersion analysis is performed to understand the effect adaptive kernel on the stability. Finally, the efficacy of the algorithm is demonstrated through some benchmark elastic dynamics problems.

READ FULL TEXT

page 14

page 27

page 30

research
07/12/2023

An adaptive approach to remove tensile instability in SPH for weakly compressible fluids

Smoothed Particle Hydrodynamics (SPH) is plagued by the phenomenon of te...
research
08/01/2021

Fast and numerically stable particle-based online additive smoothing: the AdaSmooth algorithm

We present a novel sequential Monte Carlo approach to online smoothing o...
research
06/22/2014

On the Maximum Entropy Property of the First-Order Stable Spline Kernel and its Implications

A new nonparametric approach for system identification has been recently...
research
10/28/2022

Equally spaced points are optimal for Brownian Bridge kernel interpolation

In this paper we show how ideas from spline theory can be used to constr...
research
08/10/2021

Numerical simulation of hydraulic fracturing: a hybrid FEM based algorithm

In this paper a problem of numerical simulation of hydraulic fractures i...
research
05/06/2020

Mathematical foundations of stable RKHSs

Reproducing kernel Hilbert spaces (RKHSs) are key spaces for machine lea...
research
04/30/2015

On the estimation of initial conditions in kernel-based system identification

Recent developments in system identification have brought attention to r...

Please sign up or login with your details

Forgot password? Click here to reset