Remapping method for transferring data between two meshes using a modified iterative SPR approach for parallel resolution
Résumé
In metal forming simulations, the transfer of data from one mesh to another can be very often required such as, in ALE formulation at each time step or in Lagrangian formulation at each remeshing step. Its accuracy is one of the main concerns for researchers, since the accumulation of generated diffusion can lead to larger numerical error and create convergence problems. Since 1992, the Super convergent Patch Recovery method (SPR) introduced by Zienkiewicz & Zhu [1], was a major breakthrough for the methods to estimate errors in FE solution. Later, it has also been used to recover nodal fields from integration points, in order to transfer data between two meshes [10]. The original method raises some difficulties to treat the domain boundaries. They have not quite properly been handled, in particular in the frame of parallel implementation, despite the fact that, surface phenomena, such as contact and friction, play such an important role in metal forming applications. In the present paper, it is presented a modified iterative SPR method which deals with boundary points with the same order of accuracy as the interior points. It does not require increasing the patch size and it is easier to implement in parallel environment. Also, when interpolating the field on new mesh, a new and consistent technique (P0+ transport) of enriched field has been used. It involves building a P1+ field by mixing the recovered SPR nodal field with the known field at integration points. Results are presented in form of convergence rate and L2 error norm, for several analytical functions for a SPR based transfer operator against a simple volumetric based transfer operator, before being applied to an actual metal forming problem. All computations presented here have been done by using four processors in parallel environment.