A PARTICLE TRACKING TECHNIQUE FOR THE LAGRANGIAN–EULERIAN FINITE ELEMENT METHOD IN MULTI-DIMENSIONS

This paper presents a multi-dimensional particle tracking technique for applying the Lagrangian–Eulerian finite element method to solve transport equations in transient-state simulations. In the Lagrangian– Eulerian approach, the advection term is handled in the Lagrangian step so that the associated numerical errors can be considerably reduced. It is important to have an adequate particle tracking technique for computing advection accurately in the Lagrangian step. The particle tracking technique presented here is designed to trace fictitious particles in the real-world flow field where the flow velocity is either measured or computed at a limited number of locations. The technique, named ‘in-element’ particle tracking, traces fictitious particles on an element-by-element basis. Given a velocity field, a fictitious particle is traced one element by one element until either a boundary is encountered or the available time is completely consumed. For the tracking within an element, the element is divided into a desired number of subelements with the interpolated velocity computed at all nodes of the subelements. A fictitious particle, thus, is traced one subelement by one subelement within the element. The desired number of subelements can be determined based on the complexity of the flow field being considered. The more complicated the flow field is, the more subelements are needed to achieve accurate particle tracking results. A single-velocity approach can be used to efficiently perform particle tracking in a smooth flow field, while an average-velocity approach can be employed to increase the tracking accuracy for more complex flow fields.     

Arbitrary Lagrangian–Eulerian finite element analysis of strain localization in transient problems

Non-local models guaranty that finite element computations on strain softening materials remain sound up to failure from a theoretical and computational viewpoint. The non-locality prevents strain localization with zero global dissipation of energy, and consequently finite element calculations converge upon mesh refinements to non-zero width localization zones. One of the major drawbacks of these models is that the element size needed in order to capture the localization zone must be smaller than the internal length. Hence, the total number of degrees of freedom becomes rapidly prohibitive for most engineering applications and there is an obvious need for mesh adaptivity. This paper deals with the application of the arbitrary Lagrangian–Eulerian (ALE) formulation, well known in hydrodynamics and fluid–structure interaction problems, to transient strain localization in a non-local damageable material. It is shown that the ALE formulation which is employed in large boundary motion problems can also be well suited for non-linear transient analysis of softening materials where localization bands appear. The remeshing strategy is based on the equidistribution of an indicator that quantifies the interelement jump of a selected state variable. Two well known one-dimensional examples illustrate the capabilities of this technique: the first one deals with localization due to a propagating wave in a bar, and the second one studies the wave propagation in a hollow sphere.