Parallel, second-order and consistent remeshing transfer operators for evolving meshes with superconvergence property on surface and volume

This paper investigates several field transfer techniques that can be used to remap data between three-dimensional unstructured meshes, either after full remeshing of the computational domain or after mesh regularization resulting from an ALE (Arbitrary Lagrangian or Eulerian) formulation. The transfer is focused on state (or secondary) variables that are piecewise discontinuous and consequently only defined at integration points. The proposed methods are derived from recovery techniques that have initially been developed by Zienkiwicz et al. in the frame of error estimation. Obtaining a higher order interpolation with the recovered fields allows reducing the inescapable diffusion error resulting from the projection on the new mesh. Several variants of the method are investigated: (a) either based on nodal patches or on element patches, (b) by enforcing the balance equation in a weak sense or in a strong sense or not, (c) by using first or second interpolation orders. A special attention is paid to the accuracy of the transfer operators for surface values, which can play a first order role in several mechanical problems. In order to take into account the constraint due to parallel calculations, a new iterative approach is proposed. All methods are evaluated and compared on analytical tests functions, both for the ALE formulation and for full remeshings, before being applied to an actual metal forming problem. In all studied examples, in addition to improved accuracy, higher order convergence rates are observed both for volume and surface values, so providing quite accurate transfer operators for various applications.

[1]  David Dureisseix,et al.  Information transfer between incompatible finite element meshes: Application to coupled thermo-viscoelasticity , 2006 .

[2]  Thierry Coupez,et al.  Adaptive remeshing based on a posteriori error estimation for forging simulation , 2006 .

[3]  Zhi Zong,et al.  A modified superconvergent patch recovery method and its application to large deformation problems , 2004 .

[4]  Pierre Villon,et al.  Diffuse approximation for field transfer in non linear mechanics , 2006 .

[5]  T. Liszka An interpolation method for an irregular net of nodes , 1984 .

[6]  O. C. Zienkiewicz,et al.  Recovery procedures in error estimation and adaptivity. Part II: Adaptivity in nonlinear problems of elasto-plasticity behaviour , 1999 .

[7]  G. Touzot,et al.  Continuous Stress Fields in Finite Element Analysis , 1977 .

[8]  O. C. Zienkiewicz,et al.  Recovery procedures in error estimation and adaptivity: Adaptivity in non-linear problems of elasto-plasticity behaviour , 1998 .

[10]  Nils-Erik Wiberg,et al.  Error estimation and adaptive procedures based on superconvergent patch recovery (SPR) techniques , 1997 .

[11]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[12]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[13]  O. C. Zienkiewicz,et al.  The superconvergent patch recovery (SPR) and adaptive finite element refinement , 1992 .

[15]  Michael T. Heath,et al.  Common‐refinement‐based data transfer between non‐matching meshes in multiphysics simulations , 2004 .

[16]  M. M. Rashid,et al.  Material state remapping in computational solid mechanics , 2002 .

[17]  Zhiwei Lin,et al.  A local rezoning and remapping method for unstructured mesh , 2011, Comput. Phys. Commun..

[18]  Nils-Erik Wiberg,et al.  Improved element stresses for node and element patches using superconvergent patch recovery , 1995 .

[19]  John S. Campbell,et al.  Local and global smoothing of discontinuous finite element functions using a least squares method , 1974 .

[20]  O. C. Zienkiewicz,et al.  A simple error estimator and adaptive procedure for practical engineerng analysis , 1987 .

[21]  Len G. Margolin,et al.  Second-order sign-preserving conservative interpolation (remapping) on general grids , 2003 .

[22]  M. Ortiz,et al.  Adaptive Lagrangian modelling of ballistic penetration of metallic targets , 1997 .

[23]  Lionel Fourment,et al.  Error estimators for viscoplastic materials: application to forming processes , 1995 .

[24]  Nicholas Zabaras,et al.  Shape optimization and preform design in metal forming processes , 2000 .

[25]  O. C. Zienkiewicz,et al.  Recovery procedures in error estimation and adaptivity Part I: Adaptivity in linear problems , 1999 .

[26]  O. C. Zienkiewicz,et al.  Superconvergence and the superconvergent patch recovery , 1995 .

[27]  Lionel Fourment,et al.  Friction model for friction stir welding process simulation: Calibrations from welding experiments , 2010 .

[28]  David R. Owen,et al.  Transfer operators for evolving meshes in small strain elasto-placticity , 1996 .

[29]  C. Pavanachand,et al.  Remeshing issues in the finite element analysis of metal forming problems , 1998 .

[30]  Amir R. Khoei,et al.  The superconvergence patch recovery technique and data transfer operators in 3D plasticity problems , 2007 .

[31]  Pierre Villon,et al.  On a consistent field transfer in non linear inelastic analysis and ultimate load computation , 2008 .

[32]  N.-E. Wiberg Superconvergent patch recovery—a key to quality assessed FE solutions , 1997 .

[33]  T. Liszka,et al.  The finite difference method at arbitrary irregular grids and its application in applied mechanics , 1980 .

[34]  Bijan Boroomand,et al.  RECOVERY BY EQUILIBRIUM IN PATCHES (REP) , 1997 .

[35]  Marc G. D. Geers,et al.  A robust and consistent remeshing-transfer operator for ductile fracture simulations , 2006 .

[36]  Klaus-Jürgen Bathe,et al.  Error indicators and adaptive remeshing in large deformation finite element analysis , 1994 .

[37]  J. Z. Zhu,et al.  The superconvergent patch recovery and a posteriori error estimates. Part 2: Error estimates and adaptivity , 1992 .

[38]  S. Guerdoux,et al.  A 3D numerical simulation of different phases of friction stir welding , 2009 .