Momentum advection on unstructured staggered quadrilateral meshes

The finite element method, as applied to problems in solid mechanics, typically uses a mesh with the velocities at the nodes and the remaining solution variables, including the density, located at the integration points. The arbitrary Lagrangian–Eulerian formulations used in solid mechanics are therefore faced with the challenge of transporting momentum, which is defined in terms of variables located at separate points in space, in a conservative manner. Two types of momentum transport methods have been developed over the years. The first constructs a dual mesh with the nodes as the integration points, a difficult task on an unstructured finite element mesh. The second uses the original mesh and constructs auxiliary variables for transport from which the final velocity may be recovered. An analysis demonstrates how the two methods are related. Simplified implementations of each type—dual mesh and element centered—are developed in detail and their performance is compared to verify the analysis. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Mesh motion techniques for the ALE formulation in 3D large deformation problems

Efficient mesh motion techniques are a key issue to achieve satisfactory results in the arbitrary Lagrangian–Eulerian (ALE) finite element formulation when simulating large deformation problems such as metal‐forming. In the updated Lagrangian (UL) formulation, mesh and material movement are attached and an excessive mesh distortion usually appears. By uncoupling mesh movement from material movement, the ALE formulation can relocate the mesh to avoid distortion. To facilitate the calculation process, the ALE operator is split into two steps at each analysis time step: UL step (where deformation due to loading is calculated without convective terms) and Eulerian step (where mesh motion is applied). In this work, mesh motion is performed by new nodal relocation methods, developed for eight‐node hexahedral elements, which can move internal and boundary nodes, improving and concentrating the mesh in critical zones. After mesh motion, data is transferred from the UL mesh to the relocated mesh using an expansion of stresses in a Taylor's series. Two numerical applications are presented, comparing results of UL and ALE formulation with results found in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Stress integration and mesh refinement for large deformation in geomechanics

This paper first discusses alternative stress integration schemes in numerical solutions to large‐ deformation problems in hardening materials. Three common numerical methods, i.e. the total‐Lagrangian (TL), the updated‐Lagrangian (UL) and the arbitrary Lagrangian–Eulerian (ALE) methods, are discussed. The UL and the ALE methods are further complicated with three different stress integration schemes. The objectivity of these schemes is discussed. The ALE method presented in this paper is based on the operator‐split technique where the analysis is carried out in two steps; an UL step followed by an Eulerian step. This paper also introduces a new method for mesh refinement in the ALE method. Using the known displacements at domain boundaries and material interfaces as prescribed displacements, the problem is re‐analysed by assuming linear elasticity and the deformed mesh resulting from such an analysis is then used as the new mesh in the second step of the ALE method. It is shown that this repeated elastic analysis is actually more efficient than mesh generation and it can be used for general cases regardless of problem dimension and problem topology. The relative performance of the TL, UL and ALE methods is investigated through the analyses of some classic geotechnical problems. Copyright © 2005 John Wiley & Sons, Ltd.     

An extended finite element formulation for contact in multi‐material arbitrary Lagrangian–Eulerian calculations

Multi‐material Eulerian and arbitrary Lagrangian–Eulerian methods were originally developed for solving hypervelocity impact problems, but they are attractive for solving a broad range of problems having large deformations, the evolution of new free surfaces, and chemical reactions. The contact, separation, and slip between two surfaces have traditionally been addressed by the mixture theory, however the accuracy of this approach is severely limited. To improve the accuracy, an extended finite element formulation is developed and example calculations are presented. As a side benefit, the mixture theory is eliminated from the multi‐material formulation, eliminating the issues associated with the equilibration time between adjacent materials. By design, the new formulation is relatively simple to implement in existing multi‐material codes, parallelizes without difficulty, and has a low memory burden. Copyright © 2006 John Wiley & Sons, Ltd.     

Updated Lagrangian/Arbitrary Lagrangian–Eulerian framework for interaction between a compressible neo‐Hookean structure and an incompressible fluid

We propose a numerical method for a fluid–structure interaction problem. The material of the structure is homogeneous, isotropic, and it can be described by the compressible neo‐Hookean constitutive equation, while the fluid is governed by the Navier–Stokes equations. Our study does not use turbulence model. Updated Lagrangian method is used for the structure and fluid equations are written in Arbitrary Lagrangian–Eulerian coordinates. One global moving mesh is employed for the fluid–structure domain, where the fluid–structure interface is an ‘interior boundary’ of the global mesh. At each time step, we solve a monolithic system of unknown velocity and pressure defined on the global mesh. The continuity of velocity at the interface is automatically satisfied, while the continuity of stress does not appear explicitly in the monolithic fluid–structure system. This method is very fast because at each time step, we solve only one linear system. This linear system was obtained by the linearization of the structure around the previous position in the updated Lagrangian formulation and by the employment of a linear convection term for the fluid. Numerical results are presented. Copyright © 2016 John Wiley & Sons, Ltd.     

Universal meshes for smooth surfaces with no boundary in three dimensions

We introduce a method to mesh the boundary Γ of a smooth, open domain in immersed in a mesh of tetrahedra. The mesh follows by mapping a specific collection of triangular faces in the mesh to Γ. Two types of surface meshes follow: (a) a mesh that exactly meshes Γ, and (b) meshes that approximate Γ to any order, by interpolating the map over the selected faces; that is, an isoparametric approximation to Γ. The map we use to deform the faces is the closest point projection to Γ. We formulate conditions for the closest point projection to define a homeomorphism between each face and its image. These are conditions on some of the tetrahedra intersected by the boundary, and they essentially state that each such tetrahedron should (a) have a small enough diameter, and (b) have two of its dihedral angles be acute. We provide explicit upper bounds on the mesh size, and these can be computed on the fly. We showcase the quality of the resulting meshes with several numerical examples. More importantly, all surfaces in these examples were meshed with a single background mesh. This is an important feature for problems in which the geometry evolves or changes, because it could be possible for the background mesh to never change as the geometry does. In this case, the background mesh would be a universal mesh 1 for all these geometries. We expect the method introduced here to be the basis for the construction of universal meshes for domains in three dimensions. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation of fluid–structure interactions using a velocity‐linked P2/P1 finite element method and the generalized‐ α method

A velocity‐linked algorithm for solving unsteady fluid–structure interaction (FSI) problems in a fully coupled manner is developed using the arbitrary Lagrangian–Eulerian method. The P2/P1 finite element is used to spatially discretize the incompressible Navier–Stokes equations and structural equations, and the generalized‐ α method is adopted for temporal discretization. Common velocity variables are employed at the fluid–structure interface for the strong coupling of both equations. Because of the velocity‐linked formulation, kinematic compatibility is automatically satisfied and forcing terms do not need to be calculated explicitly. Both the numerical stability and the convergence characteristics of an iterative solver for the coupled algorithm are investigated by solving the FSI problem of flexible tube flows. It is noteworthy that the generalized‐ α method with small damping is free from unstable velocity fields. However, the convergence characteristics of the coupled system deteriorate greatly for certain Poisson's ratios so that direct solvers are essential for these cases. Furthermore, the proposed method is shown to clearly display the advantage of considering FSI in the simulation of flexible tube flows, while enabling much larger time‐steps than those adopted in some previous studies. This is possible through the strong coupling of the fluid and structural equations by employing common primitive variables. Copyright © 2012 John Wiley & Sons, Ltd.     

Contact with friction in multi‐material arbitrary Lagrangian‐Eulerian formulations using X‐FEM

A contact method with friction for the multi‐dimensional Lagrangian step in multi‐material arbitrary Lagrangian–Eulerian (ALE) formulations is presented. In our previous research, the extended finite element method (X‐FEM) was used to create independent fields (i.e. velocity, strain rate, force, mass, etc.) for each material in the problem to model contact without friction. The research presented here includes the extension to friction and improvements to the accuracy and robustness of our previous study. The accelerations of the multi‐material nodes are obtained by coupling the material force and mass fields as a function of the prescribed contact; similarly, the velocities of the multi‐material nodes are recalculated using the conservation of momentum when the prescribed contact requires it. The coupling procedures impose the same nodal velocity on the coupled materials in the direction normal to their interface during the time step update. As a result, the overlap of materials is prevented and unwanted separation does not occur. Three different types of contacts are treated: perfectly bonded, frictionless slip, and slip with friction. Example impact problems are solved and the numerical solutions are presented. Copyright © 2008 John Wiley & Sons, Ltd.     

An efficient solution method for finite element ring‐rolling simulation

Finite element ring‐rolling simulation gives rise to poor conditioned non‐linear equations that require repeated solution. The associated computational costs are extreme making analysis impracticable in industry. This paper is concerned with a solution strategy that addresses this problem and involves the combined use of an arbitrary Lagrangian–Eulerian (ALE) formulation and a successive preconditioned conjugate gradient method (SPCGM). This approach, coupled to a finite element flow formulation, is shown to offer considerable computational savings. Through the combined use of the ALE flow formulation and the SPCGM the stability and condition of the non‐linear systems is enhanced. This purely iterative approach takes advantage of the slowly evolving velocity field and the self‐preconditioning offered by the SPCGM. The performance of the solver is compared against well‐known alternatives for varying problem sizes. The approach is shown to be generic but in particular makes ring‐rolling simulation a more practicable proposition. Copyright © 2000 John Wiley & Sons, Ltd.     

A three‐dimensional hybrid meshfree‐Cartesian scheme for fluid–body interaction

A numerical method based on a hybrid meshfree‐Cartesian grid is developed for solving three‐dimensional fluid–solid interaction (FSI) problems involving solid bodies undergoing large motion. The body is discretized and enveloped by a cloud of meshfree nodes. The motion of the body is tracked by convecting the meshfree nodes against a background of Cartesian grid points. Spatial discretization of second‐order accuracy is accomplished by the combination of a generalized finite difference (GFD) method and conventional finite difference (FD) method, which are applied to the meshfree and Cartesian nodes, respectively. Error minimization in GFD is carried out by singular value decomposition. The discretized equations are integrated in time via a second‐order fractional step projection method. A time‐implicit iterative procedure is employed to compute the new/evolving position of the immersed bodies together with the dynamically coupled solution of the flow field. The present method is applied on problems of free falling spheres and tori in quiescent flow and freely rotating spheres in simple shear flow. Good agreement with published results shows the ability of the present hybrid meshfree‐Cartesian grid scheme to achieve good accuracy. An application of the method to the self‐induced propulsion of a deforming fish‐like swimmer further demonstrates the capability and potential of the present approach for solving complex FSI problems in 3D. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of flapping wings using a panel method and a high‐order Navier–Stokes solver

The design of efficient flapping wings for human engineered micro aerial vehicles (MAVs) has long been an elusive goal, in part because of the large size of the design space. One strategy for overcoming this difficulty is to use a multifidelity simulation strategy that appropriately balances computation time and accuracy. We compare two models with different geometric and physical fidelity. The low‐fidelity model is an inviscid doublet lattice method with infinitely thin lifting surfaces. The high‐fidelity model is a high‐order accurate discontinuous Galerkin Navier–Stokes solver, which uses an accurate representation of the flapping wing geometry. To compare the performance of the two methods, we consider a model flapping wing with an elliptical planform and an analytically prescribed spanwise wing twist, at size scales relevant to MAVs. Our results show that in many cases, including those with mild separation, low‐fidelity simulations can accurately predict integrated forces, provide insight into the flow structure, indicate regions of likely separation, and shed light on design–relevant quantities. But for problems with significant levels of separation, higher‐fidelity methods are required to capture the details of the flow field. Inevitably high‐fidelity simulations are needed to establish the limits of validity of the lower fidelity simulations.Copyright © 2012 John Wiley & Sons, Ltd.     

A coupled BEM and arbitrary Lagrangian–Eulerian FEM model for the solution of two‐dimensional laminar flows in external flow fields

This paper describes a new computational model developed to solve two‐dimensional incompressible viscous flow problems in external flow fields. The model based on the Navier–Stokes equations in primitive variables is able to solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the pressure projection method. The external flow field is simulated using the boundary element method by solving a pressure Poisson equation that assumes the pressure as zero at the infinite boundary. The momentum equation of the flow motion is solved using the three‐step finite element method. The arbitrary Lagrangian–Eulerian method is incorporated into the model, to solve the moving boundary problems. The present model is applied to simulate various external flow problems like flow across circular cylinder, acceleration and deceleration of the circular cylinder moving in a still fluid and vibration of the circular cylinder induced by the vortex shedding. The simulation results are found to be very reasonable and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

On the numerical implementation of variational arbitrary Lagrangian–Eulerian (VALE) formulations

This paper is concerned with the implementation of variational arbitrary Lagrangian–Eulerian formulations, also known as variational r‐adaption methods. These methods seek to minimize the energy function with respect to the finite‐element mesh over the reference configuration of the body. We propose a solution strategy based on a viscous regularization of the configurational forces. This procedure eliminates the ill‐posedness of the problem without changing its solutions, i.e. the minimizers of the regularized problems are also minimizers of the original functional. We also develop strategies for optimizing the triangulation, or mesh connectivity, and for allowing nodes to migrate in and out of the boundary of the domain. Selected numerical examples demonstrate the robustness of the solution procedures and their ability to produce highly anisotropic mesh refinement in regions of high energy density. Copyright © 2006 John Wiley & Sons, Ltd.     

Two stress update algorithms for large strains: accuracy analysis and numerical implementation

Two algorithms for the stress update (i.e., time integration of the constitutive equation) in large-strain solid mechanics are compared from an analytical point of view. The order of the truncation error associated to the numerical integration is deduced for each algorithm a priori, using standard numerical analysis. This accuracy analysis has been performed by means of a convected frame formalism, which also allows a unified derivation of both algorithms in spite of their inherent differences. Then the two algorithms are adapted from convected frames to a fixed Cartesian frame and implemented in a small-strain finite element code. The implementation is validated by means of a set of simple deformation paths (simple shear, extension, extension and compression, extension and rotation) and two benchmark tests in non-linear mechanics (the necking of a circular bar and a shell under ring loads). In these numerical tests, the observed order of convergence is in very good agreement with the theoretical order of convergence, thus corroborating the accuracy analysis. © 1997 John Wiley & Sons, Ltd.     

A monolithic approach for interaction of incompressible viscous fluid and an elastic body based on fluid pressure Poisson equation

This paper describes a new monolithic approach based on the fluid pressure Poisson equation (PPE) to solve an interaction problem of incompressible viscous fluid and an elastic body. The PPE is derived so as to be consistent with the coupled equation system for the fluid‐structure interaction (FSI). Based on this approach, we develop two kinds of efficient monolithic methods. In both methods, the fluid pressure is derived implicitly so as to satisfy the incompressibility constraint, and all other unknown variables are derived fully explicitly or partially explicitly. The coefficient matrix of the PPE for the FSI becomes symmetric and positive definite and its condition is insensitive to inhomogeneity of material properties. The arbitrary Lagrangian–Eulerian (ALE) method is employed for the fluid part in order to take into account the deformable fluid‐structure interface. To demonstrate fundamental performances of the proposed approach, the developed two monolithic methods are applied to evaluate the added mass and the added damping of a circular cylinder as well as to simulate the vibration of a rectangular cylinder induced by vortex shedding. Copyright © 2005 John Wiley & Sons, Ltd.     

An implicit finite wear contact formulation based on dual mortar methods

Finite deformation contact problems with frictional effects and finite shape changes due to wear are investigated. To capture the finite shape changes, a third configuration besides the well‐known reference and spatial configurations is introduced, which represents the time‐dependent worn state. Consistent interconnections between these states are realized by an arbitrary Lagrangean–Eulerian formulation. The newly developed partitioned and fully implicit algorithm is based on a Lagrangean step and a shape evolution step. Within the Lagrangean step, contact constraints as well as the wear equations are weakly enforced following the well‐established mortar framework. Additional unknowns due to the employed Lagrange multiplier method for contact constraint enforcement and due to wear itself are eliminated by condensation procedures based on the concept of biorthogonality and the so‐called dual shape functions. Several numerical examples in both 2D and 3D are provided to demonstrate the performance and accuracy of the proposed numerical algorithm. Copyright © 2016 John Wiley & Sons, Ltd.