TY - JOUR
T1 - Stress recovery for the particle-in-cell finite element method
AU - Yang, Haibin
AU - Moresi, Louis N.
AU - Mansour, John
N1 - Publisher Copyright:
© 2020
PY - 2021/2
Y1 - 2021/2
N2 - The particle-in-cell finite element (PIC-FE) method has been widely used in geodynamic numerical modelling due to its efficiency in dealing with large deformations without the requirement of remeshing. However, material deformation within a Eulerian mesh frame will mix particles of contrasting strength properties (e.g., viscosity in Stokes problems) in a single element requiring some form of averaging to project particle properties to integration points. The numerical solutions are thus dependent on the way how the particle properties are projected to the integration points. An intra-element property discontinuity may introduce severe stress oscillations along the interfaces. In this study, we assess three preprocessing methods to smooth the viscosity contrast within one element. For simplified models with analytical solutions, the accuracy and convergence rate in L2 norm are systematically studied with ensembles. It is found that using higher-order quadrature elements does not improve the convergence rate for either the velocity or stress solution, both close to one. Additionally, the convergence rate of the maximum stress error, which exists adjacent to the mixed-material elements, is much less than one for all cases studied here. Comparing each component of the stress tensor, we find that the stress tensor component with the highest strain rate gradient across the material interface produces the maximum stress error. Such errors can be reduced by averaging the particle properties to the Gaussian quadrature point with an inverse-distance-weighted harmonic mean.
AB - The particle-in-cell finite element (PIC-FE) method has been widely used in geodynamic numerical modelling due to its efficiency in dealing with large deformations without the requirement of remeshing. However, material deformation within a Eulerian mesh frame will mix particles of contrasting strength properties (e.g., viscosity in Stokes problems) in a single element requiring some form of averaging to project particle properties to integration points. The numerical solutions are thus dependent on the way how the particle properties are projected to the integration points. An intra-element property discontinuity may introduce severe stress oscillations along the interfaces. In this study, we assess three preprocessing methods to smooth the viscosity contrast within one element. For simplified models with analytical solutions, the accuracy and convergence rate in L2 norm are systematically studied with ensembles. It is found that using higher-order quadrature elements does not improve the convergence rate for either the velocity or stress solution, both close to one. Additionally, the convergence rate of the maximum stress error, which exists adjacent to the mixed-material elements, is much less than one for all cases studied here. Comparing each component of the stress tensor, we find that the stress tensor component with the highest strain rate gradient across the material interface produces the maximum stress error. Such errors can be reduced by averaging the particle properties to the Gaussian quadrature point with an inverse-distance-weighted harmonic mean.
KW - Finite element method
KW - Numerical geodynamic modelling
KW - Particle-in-cell
KW - Stress fluctuation
KW - Stress smoothing
UR - http://www.scopus.com/inward/record.url?scp=85098935529&partnerID=8YFLogxK
U2 - 10.1016/j.pepi.2020.106637
DO - 10.1016/j.pepi.2020.106637
M3 - Article
SN - 0031-9201
VL - 311
JO - Physics of the Earth and Planetary Interiors
JF - Physics of the Earth and Planetary Interiors
M1 - 106637
ER -