Material point method enhanced by modified gradient of shape function

Duan Zhang

doi:10.1016/j.jcp.2011.04.032

2011-07-01

158

PRE

Duan Zhang *

Xia Ma

P. T. Giguere

DOI：10.1016/j.jcp.2011.04.032

原文链接

原文求助

摘要

En 中文

A numerical scheme of computing quantities involving gradients of shape functions is introduced for the material point method (MPM), so that the quantities are continuous as material points move across cell boundaries. The noise and instability caused by cell crossing of the material points are then eliminated. In this scheme, the formulas used to compute these quantities can be expressed in the same forms as in the original material point method, but with the gradient of the shape function modified. For one-dimensional cases, the gradient of the shape function used in the generalized interpolation material point (GIMP) method is a special case of the modified gradient if the characteristic function of a material point is introduced. The characteristic function of a material point is not otherwise needed in this scheme, therefore difficulties in tracking its evolution are avoided. Although the support of the modified gradient of a shape function is enlarged from the cell containing the material point to also include the immediate neighbor cells, all the non-local effects of a material point can be accounted for by two consecutive local operations. Therefore this scheme can be used in calculations with unstructured grids. This scheme is proved to satisfy mass and momentum conservations exactly. The error in energy conservation is shown to be second order on both spatial and temporal discretizations. Although the error in energy conservation is the same order as that in the original material point method, numerical examples show that this scheme has significantly better energy conservation properties than those of the original material point method. (C) 2011 Elsevier Inc. All rights reserved.

Keyword:

Particle-in-cell

Material point method

Multiphase flow

Fluid-structure interaction