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.     

Clustered generalized finite element methods for mesh unrefinement,non‐matching and invalid meshes

In spite of significant advancements in automatic mesh generation during the past decade, the construction of quality finite element discretizations on complex three‐dimensional domains is still a difficult and time demanding task. In this paper, the partition of unity framework used in the generalized finite element method (GFEM) is exploited to create a very robust and flexible method capable of using meshes that are unacceptable for the finite element method, while retaining its accuracy and computational efficiency. This is accomplished not by changing the mesh but instead by clustering groups of nodes and elements. The clusters define a modified finite element partition of unity that is constant over part of the clusters. This so‐called clustered partition of unity is then enriched to the desired order using the framework of the GFEM. The proposed generalized finite element method can correctly and efficiently deal with: (i) elements with negative Jacobian; (ii) excessively fine meshes created by automatic mesh generators; (iii) meshes consisting of several sub‐domains with non‐matching interfaces. Under such relaxed requirements for an acceptable mesh, and for correctly defined geometries, today's automated tetrahedral mesh generators can practically guarantee successful volume meshing that can be entirely hidden from the user. A detailed technical discussion of the proposed generalized finite element method with clustering along with numerical experiments and some implementation details are presented. Copyright © 2006 John Wiley & Sons, Ltd.     

Energy-based r-adaptivity: a solution strategy and applications to fracture mechanics

This paper deals with energy based r-adaptivity in finite hyperelastostatics. The focus lies on the development of a numerical solution strategy. Although the concept of improving the accuracy of a finite element solution by minimizing the discrete potential energy with respect to the material node point positions is well-known, the numerical implementation of the underlying minimization problem is difficult. In this paper, energy based r-adaptivity is defined as a minimization problem with inequality constraints. The constraints are introduced to restrict the maximum distortion of the finite element mesh. As a solution strategy for the constrained problem, we use a classical barrier method. Beside the theoretical aspects and the implementation, a numerical experiment is presented. We illustrate the performance of the proposed r-adaptivity in the case of a cracked specimen.     

An error‐estimate‐free and remapping‐free variational mesh refinement and coarsening method for dissipative solids at finite strains

A variational h‐adaptive finite element formulation is proposed. The distinguishing feature of this method is that mesh refinement and coarsening are governed by the same minimization principle characterizing the underlying physical problem. Hence, no error estimates are invoked at any stage of the adaption procedure. As a consequence, linearity of the problem and a corresponding Hilbert‐space functional framework are not required and the proposed formulation can be applied to highly non‐linear phenomena. The basic strategy is to refine (respectively, unrefine) the spatial discretization locally if such refinement (respectively, unrefinement) results in a sufficiently large reduction (respectively, sufficiently small increase) in the energy. This strategy leads to an adaption algorithm having O(N) complexity. Local refinement is effected by edge‐bisection and local unrefinement by the deletion of terminal vertices. Dissipation is accounted for within a time‐discretized variational framework resulting in an incremental potential energy. In addition, the entire hierarchy of successive refinements is stored and the internal state of parent elements is updated so that no mesh‐transfer operator is required upon unrefinement. The versatility and robustness of the resulting variational adaptive finite element formulation is illustrated by means of selected numerical examples. Copyright © 2008 John Wiley & Sons, Ltd.     

hp‐Generalized FEM and crack surface representation for non‐planar 3‐D cracks

A high‐order generalized finite element method (GFEM) for non‐planar three‐dimensional crack surfaces is presented. Discontinuous p‐hierarchical enrichment functions are applied to strongly graded tetrahedral meshes automatically created around crack fronts. The GFEM is able to model a crack arbitrarily located within a finite element (FE) mesh and thus the proposed method allows fully automated fracture analysis using an existing FE discretization without cracks. We also propose a crack surface representation that is independent of the underlying GFEM discretization and controlled only by the physics of the problem. The representation preserves continuity of the crack surface while being able to represent non‐planar, non‐smooth, crack surfaces inside of elements of any size. The proposed representation also provides support for the implementation of accurate, robust, and computationally efficient numerical integration of the weak form over elements cut by the crack surface. Numerical simulations using the proposed GFEM show high convergence rates of extracted stress intensity factors along non‐planar curved crack fronts and the robustness of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

An element‐free method for in‐plane notch problems with anisotropic materials

This study developed an element‐free Galerkin method (EFGM) to simulate notched anisotropic plates containing stress singularities at the notch tip. Two‐dimensional theoretical complex displacement functions are first deduced into the moving least‐squares interpolation. The interpolation functions and their derivatives are then determined to calculate the nodal stiffness using the Galerkin method. In the numerical validation, an interface layer of the EFGM is used to combine the mesh between the traditional finite elements and the proposed singular notch EFGM. The H‐integral determined from finite element analyses with a very fine mesh is used to validate the numerical results of the proposed method. The comparisons indicate that the proposed method obtains more accurate results for the displacement, stress, and energy fields than those determined from the standard finite element method. Copyright © 2013 John Wiley & Sons, Ltd.     

A posteriori finite element bounds to linear functional outputs of the three‐dimensional Navier–Stokes equations

An implicit a posteriori finite element error estimation method is presented to inexpensively calculate lower and upper bounds for a linear functional output of the numerical solutions to the three‐dimensional Navier–Stokes (N–S) equations. The novelty of this research is to utilize an augmented Lagrangian based on a coarse mesh linearization of the N–S equations and the finite element tearing and interconnecting (FETI) procedure. The latter approach extends the a posteriori bound method to the three‐dimensional Crouzeix–Raviart space for N–S problems. The computational advantage of the bound procedure is that a single coupled non‐symmetric large problem can be decomposed into several uncoupled symmetric small problems. A simple model problem, which is selected here to illustrate the procedure, is to find the lower and upper bounds of the average velocity of a pressure driven, incompressible, steady Newtonian fluid flow moving at low Reynolds numbers through an endless square channel which has an array of rectangular obstacles. Numerical results show that the bounds for this output are rigorous, i.e. always in the asymptotic certainty regime, that they are sharp and that the required computational resources decrease significantly. Parallel implementation on a Beowulf cluster is also reported. Copyright © 2004 John Wiley & Sons, Ltd.     

Subdivision schemes for smooth contact surfaces of arbitrary mesh topology in 3D

This paper presents a strategy to parameterize contact surfaces of arbitrary mesh topology in 3D with at least C¹‐continuity for both quadrilateral and triangular meshes. In the regular mesh domain, four quadrilaterals or six triangles meet in one node, even C²‐continuity is attained. Therefore, we use subdivision surfaces, for which non‐physical pressure jumps are avoided for contact interactions. They are usually present when the contact kinematics is based on facet elements discretizing the interacting bodies. The properties of subdivision surfaces give rise to basically four different implementation strategies. Each strategy has specific features and requires more or less efforts for an implementation in a finite element program. One strategy is superior with respect to the others in the sense that it does not use nodal degrees of freedom of the finite element mesh at the contact surface. Instead, it directly uses the degrees of freedom of the smooth surface. Thereby, remarkably, it does not require an interpolation. We show how the proposed method can be used to parameterize adaptively refined meshes with hanging nodes. This is essential when dealing with finite element models whose geometry is generated by means of subdivision techniques. Three numerical 3D problems demonstrate the improved accuracy, robustness and performance of the proposed method over facet‐based contact surfaces. In particular, the third problem, adopted from biomechanics, shows the advantages when designing complex contact surfaces by means of subdivision techniques. Copyright © 2004 John Wiley & Sons, Ltd.     

Molecular mechanics in the context of the finite element method

In molecular mechanics, the formalism of the finite element method can be exploited in order to analyze the behavior of atomic structures in a computationally efficient way. Based on the atom‐related consideration of the atomic interactions, a direct correlation between the type of the underlying interatomic potential and the design of the related finite element is established. Each type of potential is represented by a specific finite element. A general formulation that unifies the various finite elements is proposed. Arbitrary diagonal‐ and cross‐terms dependent on bond length, valence angle, dihedral angle, improper dihedral angle and inversion angle can also be considered. The finite elements are formulated in a geometrically exact setting; the related formulas are stated in detail. The mesh generation can be performed using well‐known procedures typically used in molecular dynamics. Although adjacent elements overlap, a double counting of the element contributions (as a result of the assembly process) cannot occur a priori. As a consequence, the assembly process can be performed efficiently line by line. The presented formulation can easily be implemented in standard finite element codes; thus, already existing features (e.g. equation solver, visualization of the numerical results) can be employed. The formulation is applied to various interatomic potentials that are frequently used to describe the mechanical behavior of carbon nanotubes. The effectiveness and robustness of this method are demonstrated by means of several numerical examples. Copyright © 2008 John Wiley & Sons, Ltd.     

Effective volume‐fraction optimization for thermal stress reduction in FGMs using irregular h‐refinement

The volume fraction of constituent particles in functionally graded materials (FGM) varies continuously and functionally, and its optimal tailoring could be made by the numerical optimization incorporated with the finite element method. In such a case, the mesh density in finite element discretization of the volume fraction field influences the final design quality such that the further objective function reduction requires the refinement of volume fraction meshes. However, the uniform refinement of the volume fraction mesh is not effective, from the numerical point of view, particularly when the finite difference scheme is employed. This numerical inefficiency could be resolved by locally increasing the mesh density only where more volume‐fraction flexibility (i.e. more mesh density) is required. In this paper, we propose an effective volume‐fraction optimization procedure by applying the irregular h‐refinement to the volume fraction discretization. Copyright © 2003 John Wiley & Sons, Ltd.     

The maximum principle violations of the mixed‐hybrid finite‐element method applied to diffusion equations

The abundant literature of finite‐element methods applied to linear parabolic problems, generally, produces numerical procedures with satisfactory properties. However, some initial–boundary value problems may cause large gradients at some points and consequently jumps in the solution that usually needs a certain period of time to become more and more smooth. This intuitive fact of the diffusion process necessitates, when applying numerical methods, varying the mesh size (in time and space) according to the smoothness of the solution. In this work, the numerical behaviour of the time‐dependent solutions for such problems during small time duration obtained by using a non‐conforming mixed‐hybrid finite‐element method (MHFEM) is investigated. Numerical comparisons with the standard Galerkin finite element (FE) as well as the finite‐difference (FD) methods are checked. Owing to the fact that the mixed methods violate the discrete maximum principle, some numerical experiments showed that the MHFEM leads sometimes to non‐physical peaks in the solution. A diffusivity criterion relating the mesh steps for an artificial initial–boundary value problem will be presented. One of the propositions given to avoid any non‐physical oscillations is to use the mass‐lumping techniques. Copyright © 2002 John Wiley & Sons, Ltd.     

A p‐hierarchical adaptive procedure for the scaled boundary finite element method

This paper describes a p‐hierarchical adaptive procedure based on minimizing the classical energy norm for the scaled boundary finite element method. The reference solution, which is the solution of the fine mesh formed by uniformly refining the current mesh element‐wise one order higher, is used to represent the unknown exact solution. The optimum mesh is assumed to be obtained when each element contributes equally to the global error. The refinement criteria and the energy norm‐based error estimator are described and formulated for the scaled boundary finite element method. The effectivity index is derived and used to examine quality of the proposed error estimator. An algorithm for implementing the proposed p‐hierarchical adaptive procedure is developed. Numerical studies are performed on various bounded domain and unbounded domain problems. The results reflect a number of key points. Higher‐order elements are shown to be highly efficient. The effectivity index indicates that the proposed error estimator based on the classical energy norm works effectively and that the reference solution employed is a high‐quality approximation of the exact solution. The proposed p‐hierarchical adaptive strategy works efficiently. Copyright © 2007 John Wiley & Sons, Ltd.     

Investigations on boundary temperature control analysis considering a moving body based on the adjoint variable and the fictitious domain finite element methods

This paper presents an examination of moving‐boundary temperature control problems. With a moving‐boundary problem, a finite‐element mesh is generated at each time step to express the position of the boundary. On the other hand, if an overlapped domain, that is, comprising foreground and background meshes, is prepared, the moving boundary problem can be solved without mesh generation at each time step by using the fictitious domain method. In this study, boundary temperature control problems with a moving boundary are formulated using the finite element, the adjoint variable, and the fictitious domain methods, and several numerical experiments are carried out. Copyright © 2015 John Wiley & Sons, Ltd.     

Formulation and implementation of a high-order 3-D domain integral method for the extraction of energy release rates

This article presents a three dimensional (3-D) formulation and implementation of a high-order domain integral method for the computation of energy release rate. The method is derived using surface and domain formulations of the J-integral and the weighted residual method. The J-integral along 3-D crack fronts is approximated by high-order Legendre polynomials. The proposed implementation is tailored for the Generalized/eXtended Finite Element Method and can handle discontinuities arbitrarily located within a finite element mesh. The domain integral calculations are based on the same integration elements used for the computation of the stiffness matrix. Discontinuities of the integrands across crack surfaces and across computational element boundaries are fully accounted for. The proposed method is able to deliver smooth approximations and to capture the boundary layer behavior of the J-integral using tetrahedral meshes. Numerical simulations of mode-I and mixed mode benchmark fracture mechanics examples verify expected convergence rates for the computed energy release rates. The results are also in good agreement with other numerical solutions available in the literature.     

Numerical simulations of strain localization in inelastic solids using mesh‐free methods

In this paper, a comprehensive account on using mesh‐free methods to simulate strain localization in inelastic solids is presented. Using an explicit displacement‐based formulation in mesh‐free computations, high‐resolution shear‐band formations are obtained in both two‐dimensional (2‐D) and three‐dimensional (3‐D) simulations without recourse to any mixed formulation, discontinuous/incompatible element or special mesh design. The numerical solutions obtained here are insensitive to the orientation of the particle distributions if the local particle distribution is quasi‐uniform, which, to a large extent, relieves the mesh alignment sensitivity that finite element methods suffer. Moreover, a simple h‐adaptivity procedure is implemented in the explicit calculation, and by utilizing a mesh‐free hierarchical partition of unity a spectral (wavelet) adaptivity procedure is developed to seek high‐resolution shear‐band formations. Moreover, the phenomenon of multiple shear band and mode switching are observed in numerical computations with a relatively coarse particle distribution in contrast to the costly fine‐scale finite element simulations. Copyright © 2000 John Wiley & Sons, Ltd.     

Conservation properties of a time FE method. Part I: time‐stepping schemes for N‐body problems

In the present paper one‐step implicit integration algorithms for the N‐body problem are developed. The time‐stepping schemes are based on a Petrov–Galerkin finite element method applied to the Hamiltonian formulation of the N‐body problem. The approach furnishes algorithmic energy conservation in a natural way. The proposed time finite element method facilitates a systematic implementation of a family of time‐stepping schemes. A particular algorithm is specified by the associated quadrature rule employed for the evaluation of time integrals. The influence of various standard as well as non‐standard quadrature formulas on algorithmic energy conservation and conservation of angular momentum is examined in detail for linear and quadratic time elements. Copyright © 2000 John Wiley & Sons, Ltd.     

An adaptive non‐conforming finite‐element method for Reissner–Mindlin plates

Adaptive algorithms are important tools for efficient finite‐element mesh design. In this paper, an error controlled adaptive mesh‐refining algorithm is proposed for a non‐conforming low‐order finite‐element method for the Reissner–Mindlin plate model. The algorithm is controlled by a reliable and efficient residual‐based a posteriori error estimate, which is robust with respect to the plate's thickness. Numerical evidence for this and the efficiency of the new algorithm is provided in the sense that non‐optimal convergence rates are optimally improved in our numerical experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

Multi‐level hp‐adaptivity for cohesive fracture modeling

Discretization‐induced oscillations in the load–displacement curve are a well‐known problem for simulations of cohesive crack growth with finite elements. The problem results from an insufficient resolution of the complex stress state within the cohesive zone ahead of the crack tip. This work demonstrates that the hp‐version of the finite element method is ideally suited to resolve this complex and localized solution characteristic with high accuracy and low computational effort. To this end, we formulate a local and hierarchic mesh refinement scheme that follows dynamically the propagating crack tip. In this way, the usually applied static a priori mesh refinement along the complete potential crack path is avoided, which significantly reduces the size of the numerical problem. Studying systematically the influence of h‐refinement, p‐refinement, and hp‐refinement, we demonstrate why the suggested hp‐formulation allows to capture accurately the complex stress state at the crack front preventing artificial snap‐through and snap‐back effects. This allows to decrease significantly the number of degrees of freedom and the simulation runtime. Furthermore, we show that by combining this idea with the finite cell method, the crack propagation within complex domains can be simulated efficiently without resolving the geometry by the mesh. Copyright © 2016 John Wiley & Sons, Ltd.     

A parallel implicit/explicit hybrid time domain method for computational electromagnetics

The numerical solution of Maxwell's curl equations in the time domain is achieved by combining an unstructured mesh finite element algorithm with a cartesian finite difference method. The practical problem area selected to illustrate the application of the approach is the simulation of three‐dimensional electromagnetic wave scattering. The scattering obstacle and the free space region immediately adjacent to it are discretized using an unstructured mesh of linear tetrahedral elements. The remainder of the computational domain is filled with a regular cartesian mesh. These two meshes are overlapped to create a hybrid mesh for the numerical solution. On the cartesian mesh, an explicit finite difference method is adopted and an implicit/explicit finite element formulation is employed on the unstructured mesh. This approach ensures that computational efficiency is maintained if, for any reason, the generated unstructured mesh contains elements of a size much smaller than that required for accurate wave propagation. A perfectly matched layer is added at the artificial far field boundary, created by the truncation of the physical domain prior to the numerical solution. The complete solution approach is parallelized, to enable large‐scale simulations to be effectively performed. Examples are included to demonstrate the numerical performance that can be achieved. Copyright © 2009 John Wiley & Sons, Ltd.