Energy conserving and dissipative time finite element schemes for N‐body stiff problems

Petrov–Galerkin finite element method is adopted to develop a family of temporal integrators, which preserves the feature of energy conservation or numerical dissipation for non‐linear N‐body dynamical systems. This leads to an enhancement of numerical stability and the integrators may therefore offer some advantage for the numerical solution of stiff systems in long‐term simulations. Dynamically tuneable numerical integration is exploited to improve the accuracy of the time‐stepping schemes. Representative simulations for simple non‐linear systems show the performance of the schemes in controlling over or damping out unresolved high frequencies. Copyright © 2004 John Wiley & Sons, Ltd.     

On the enforcing energy conservation of time finite elements for discrete elasto‐dynamics problems

In the context of the time‐finite element method, algorithmic stresses, which enable the conservation of energy, are designed for temporal integrators derived from the midpoint and trapezoidal schemes. This is achieved through an appropriate modification of the standard midpoint and trapezoidal quadrature rules used for the numerical integration of time integrals. Either scalar scaling or vectorial adjustments can be employed for the modification, and well‐designed simple tests allow to investigate the quality of these different strategies of energy‐conserving enforcements. Numerical examples with semi‐discrete elasto‐dynamics problems are presented to show the superior stability of energy‐conserving schemes. Copyright © 2006 John Wiley & Sons, Ltd.     

Conservation properties of a time FE method—part II: Time‐stepping schemes for non‐linear elastodynamics

In the present paper one‐step implicit integration algorithms for non‐linear elastodynamics are developed. The discretization process rests on Galerkin methods in space and time. In particular, the continuous Galerkin method applied to the Hamiltonian formulation of semidiscrete non‐linear elastodynamics lies at the heart of the time‐stepping schemes. Algorithmic conservation of energy and angular momentum are shown to be closely related to quadrature formulas that are required for the calculation of time integrals. We newly introduce the ‘assumed strain method in time’ which enables the design of energy–momentum conserving schemes and which can be interpreted as temporal counterpart of the well‐established assumed strain method for finite elements in space. The numerical examples deal with quasi‐rigid motion as well as large‐strain motion. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical integration of the stiff dynamics of geometrically exact shells: an energy‐dissipative momentum‐conserving scheme

This paper presents a new family of time‐stepping algorithms for the integration of the dynamics of non‐linear shells. We consider the geometrically exact shell theory involving an inextensible director field (the so‐called five‐parameter shell model). The main characteristic of this model is the presence of the group of finite rotations in the configuration manifold describing the deformation of the solid. In this context, we develop time‐stepping algorithms whose discrete solutions exhibit the same conservation laws of linear and angular momenta as the underlying physical system, and allow the introduction of a controllable non‐negative energy dissipation to handle the high numerical stiffness characteristic of these problems. A series of algorithmic parameters for the different components of the deformation of the shell (i.e. membrane, bending and transverse shear) fully control this numerical dissipation, recovering existing energy‐momentum schemes as a particular choice of these algorithmic parameters. We present rigorous proofs of the numerical properties of the resulting algorithms in the full non‐linear range. Furthermore, it is argued that the numerical dissipation is introduced in the high‐frequency range by considering the proposed algorithm in the context of a linear problem. The finite element implementation of the resulting methods is described in detail as well as considered in the final arguments proving the aforementioned conservation/dissipation properties. We present several representative numerical simulations illustrating the performance of the newly proposed methods. The robustness gained over existing methods in these stiff problems is confirmed in particular. Copyright © 2002 John Wiley & Sons, Ltd.     

Time finite element based Moreau‐type integrators

With the postulation of the principle of virtual action, we propose, in this paper, a variational framework for describing the dynamics of finite dimensional mechanical systems, which contain frictional contact interactions. Together with the contact and impact laws formulated as normal cone inclusions, the principle of virtual action directly leads to the measure differential inclusions commonly used in the dynamics of nonsmooth mechanical systems. The discretization of the principle of virtual action in its strong and weak variational form by local finite elements in time provides a structured way to derive various time‐stepping schemes. The constitutive laws for the impulsive and nonimpulsive contact forces, ie, the contact and impact laws, are treated on velocity‐level by using a discrete contact law for the percussion increments in the sense of Moreau. Using linear shape functions and different quadrature rules, we obtain three different stepping schemes. Besides the well‐established Moreau time‐stepping scheme, we can present two alternative integrators referred to as symmetric and variational Moreau‐type stepping schemes. A suitable benchmark example shows the superiority of the newly proposed integrators in terms of energy conservation properties, accuracy, and convergence.     

Conservation properties of a time FE method. Part IV: Higher order energy and momentum conserving schemes

In the present paper a systematic development of higher order accurate time stepping schemes which exactly conserve total energy as well as momentum maps of underlying finite‐dimensional Hamiltonian systems with symmetry is shown. The result of this development is the enhanced Galerkin (eG) finite element method in time. The conservation of the eG method is generally related to its collocation property. Total energy conservation, in particular, is obtained by a new projection technique. The eG method is, moreover, based on objective time discretization of the used strain measure. This paper is concerned with particle dynamics and semi‐discrete non‐linear elastodynamics. The related numerical examples show good performance in presence of stiffness as well as for calculating large‐strain motions. Copyright © 2005 John Wiley & Sons, Ltd.     

On the enforcing energy conservation of time‐finite elements for particle systems

In the framework of the time‐finite element method, energy‐conserving temporal integrators are developed through the design of non‐standard quadrature. The latter can be achieved by a global enforcement of the force field issued from the standard quadrature. This might, however, lead to an undesirable coupling of all degrees of freedom. A local enforcement is suggested to enhance the enforcement of energy conservation through a tailored modification to each component of acting forces. Representative examples for non‐linear particle systems are presented in order to examine the quality of various approaches in the enforcement of energy conservation. Copyright © 2006 John Wiley & Sons, Ltd.     

An energy momentum consistent integration scheme using a polyconvexity‐based framework for nonlinear thermo‐elastodynamics

The present contribution provides a new approach to the design of energy momentum consistent integration schemes in the field of nonlinear thermo‐elastodynamics. The method is inspired by the structure of polyconvex energy density functions and benefits from a tensor cross product that greatly simplifies the algebra. Furthermore, a temperature‐based weak form is used, which facilitates the design of a structure‐preserving time‐stepping scheme for coupled thermoelastic problems. This approach is motivated by the general equation for nonequilibrium reversible‐irreversible coupling (GENERIC) framework for open systems. In contrast to complex projection‐based discrete derivatives, a new form of an algorithmic stress formula is proposed. The spatial discretization relies on finite element interpolations for the displacements and the temperature. The superior performance of the proposed formulation is shown within representative quasi‐static and fully transient numerical examples.     

High‐accuracy differential quadrature finite element method and its application to free vibrations of thin plate with curvilinear domain

Based on the differential quadrature (DQ) rule, the Gauss Lobatto quadrature rule and the variational principle, a DQ finite element method (DQFEM) is proposed for the free vibration analysis of thin plates. The DQFEM is a highly accurate and rapidly converging approach, and is distinct from the differential quadrature element method (DQEM) and the quadrature element method (QEM) by employing the function values themselves in the trial function for the title problem. The DQFEM, without using shape functions, essentially combines the high accuracy of the differential quadrature method (DQM) with the generality of the standard finite element formulation, and has superior accuracy to the standard FEM and FDM, and superior efficiency to the p‐version FEM and QEM in calculating the stiffness and mass matrices. By incorporating the reformulated DQ rules for general curvilinear quadrilaterals domains into the DQFEM, a curvilinear quadrilateral DQ finite plate element is also proposed. The inter‐element compatibility conditions as well as multiple boundary conditions can be implemented, simply and conveniently as in FEM, through modifying the nodal parameters when required at boundary grid points using the DQ rules. Thus, the DQFEM is capable of constructing curvilinear quadrilateral elements with any degree of freedom and any order of inter‐element compatibilities. A series of frequency comparisons of thin isotropic plates with irregular and regular planforms validate the performance of the DQFEM. Copyright © 2009 John Wiley & Sons, Ltd.     

Energy‐conserving and decaying Algorithms in non‐linear structural dynamics

A generalized formulation of the Energy‐Momentum Methodwill be developed within the framework of the Generalized‐α Methodwhich allows at the same time guaranteed conservation or decay of total energy and controllable numerical dissipation of unwanted high frequency response. Furthermore, the latter algorithm will be extended by the consistently integrated constraints of energy and momentum conservation originally derived for the Constraint Energy‐Momentum Algorithm. The goal of this general approach of implicit energy‐conserving and decaying time integration schemes is, to compare these algorithms on the basis of an equivalent notation by the means of an overall algorithmic design and hence to investigate their numerical properties. Numerical stability and controllable numerical dissipation of high frequencies will be studied in application to non‐linear structural dynamics. Among the methods considered will be the Newmark Method, the classical α‐methods, the Energy‐Momentum Methodwith and without numerical dissipation, the Constraint Energy‐Momentum Algorithm and the Constraint Energy Method. Copyright © 1999 John Wiley & Sons, Ltd.     

An optimal high‐order non‐reflecting finite element scheme for wave scattering problems

A new finite element scheme is proposed for the numerical solution of time‐harmonic wave scattering problems in unbounded domains. The infinite domain in truncated via an artificial boundary ?? which encloses a finite computational domain Ω. On ?? a local high‐order non‐reflecting boundary condition (NRBC) is applied which is constructed to be optimal in a certain sense. This NRBC is implemented in a special way, by using auxiliary variables along the boundary ??, so that it involves no high‐order derivatives regardless of its order. The order of the scheme is simply an input parameter, and it may be arbitrarily high. This leads to a symmetric finite element formulation where standard C⁰ finite elements are used in Ω. The performance of the method is demonstrated via numerical examples, and it is compared to other NRBC‐based schemes. The method is shown to be highly accurate and stable, and to lead to a well‐conditioned matrix problem. Copyright © 2002 John Wiley & Sons, Ltd.     

Finite element formulation of exact non‐reflecting boundary conditions for the time‐dependent wave equation

A modified version of an exact Non‐reflecting Boundary Condition (NRBC) first derived by Grote and Keller is implemented in a finite element formulation for the scalar wave equation. The NRBC annihilate the first N wave harmonics on a spherical truncation boundary, and may be viewed as an extension of the second‐order local boundary condition derived by Bayliss and Turkel. Two alternative finite element formulations are given. In the first, the boundary operator is implemented directly as a ‘natural’ boundary condition in the weak form of the initial–boundary value problem. In the second, the operator is implemented indirectly by introducing auxiliary variables on the truncation boundary. Several versions of implicit and explicit time‐integration schemes are presented for solution of the finite element semidiscrete equations concurrently with the first‐order differential equations associated with the NRBC and an auxiliary variable. Numerical studies are performed to assess the accuracy and convergence properties of the NRBC when implemented in the finite element method. The results demonstrate that the finite element formulation of the (modified) NRBC is remarkably robust, and highly accurate. Copyright © 1999 John Wiley & Sons, Ltd.     

Assumed strain finite element methods for conserving temporal integrations in non‐linear solid dynamics

This paper presents a new class of assumed strain finite elements to use in combination with general energy‐momentum‐conserving time‐stepping algorithms so that these conservation properties in time are preserved by the fully discretized system in space and time. The case of interest corresponds to nearly incompressible material responses, in the fully non‐linear finite strain elastic and elastoplastic ranges. The new elements consider the classical scaling of the deformation gradient with an assumed Jacobian (its determinant) defined locally through a weighted averaging procedure at the element level. The key aspect of the newly proposed formulation is the definition of the associated linearized strain operator or B‐bar operator. The developments presented here start by identifying the conditions that this discrete operator must satisfy for the fully discrete system in time and space to inherit exactly the conservation laws of linear and angular momenta, and the conservation/dissipation law of energy for elastic and inelastic problems, respectively. Care is also taken of the preservation of the relative equilibria and the corresponding group motions associated with the momentum conservation laws, and characterized by purely rotational and translational motions superimposed to the equilibrium deformed configuration. With these developments at hand, a new general B‐bar operator is introduced that satisfies these conditions. The new operator not only accounts for the spatial interpolations (e.g. bilinear displacements with piece‐wise constant volume) but also depends on the discrete structure of the equations in time. The aforementioned conservation/dissipation properties of energy and momenta are then proven to hold rigorously for the final numerical schemes, unconditionally of the time step size and the material model (elastic or elastoplastic). Different finite elements are considered in this framework, including quadrilateral and triangular elements for plane problems and brick elements for three‐dimensional problems. Several representative numerical simulations are presented involving, in particular, the use of energy‐dissipating momentum‐conserving time‐stepping schemes recently developed by the author and co‐workers for general finite strain elastoplasticity in order to illustrate the properties of the new finite elements, including these conservation/dissipation properties in time and their locking‐free response in the quasi‐incompressible case. Copyright © 2007 John Wiley & Sons, Ltd.     

A stable energy‐conserving approach for frictional contact problems based on quadrature formulas

A common approach for the numerical simulation of non‐linear multi‐body contact problems is the use of Lagrange multipliers to model the contact conditions. The stability of standard algorithms is improved by introducing a modified mass matrix which assigns no mass to the potential contact nodes. By this, the spurious algorithmic oscillations in the multiplier do not occur any more, which facilitates the application of the primal–dual active set strategy to dynamical contact problems. The new mass matrix is calculated via a modified quadrature formula that needs no extra computational cost. In addition the conservation properties of the underlying algorithm are transferred to the modified mass version. Different numerical examples for frictional two‐body contact problems illustrate the improvement in the results for the contact stresses. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical analysis of dissolution processes in cementitious materials using discontinuous and continuous Galerkin time integration schemes

The present paper is concerned with the numerical integration of non‐linear reaction–diffusion problems by means of discontinuous and continuous Galerkin methods. The first‐order semidiscrete initial value problem of calcium leaching of cementitious materials, based on a phenomenological dissolution model, an electrolyte diffusion model and the spatial p‐finite element discretization, is used as a highly non‐linear model problem. A p‐finite element method is used for the spatial discretization. In the context of discontinuous Galerkin methods the semidiscrete mass balance and the continuity of the primary variables are weakly formulated within time steps and between time steps, respectively. Continuous Galerkin methods are obtained by the strong enforcement of the continuity condition as special cases. The introduction of a natural time co‐ordinate allows for the application of standard higher order temporal shape functions of the p‐Lagrange type and the well‐known Gauss–Legendre quadrature of associated time integrals. It is shown, that arbitrary order accurate integration schemes can be developed within the framework of the proposed temporal p‐Galerkin methods. Selected benchmark analyses of calcium dissolution demonstrate the robustness of these methods with respect to pronounced changes of the reaction term and non‐smooth changes of Dirichlet boundary conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Cyclic steady states of nonlinear electro‐mechanical devices excited at resonance

We present an efficient numerical method to solve for cyclic steady states of nonlinear electro‐mechanical devices excited at resonance. Many electro‐mechanical systems are designed to operate at resonance, where the ramp‐up simulation to steady state is computationally very expensive – especially when low damping is present. The proposed method relies on a Newton–Krylov shooting scheme for the direct calculation of the cyclic steady state, as opposed to a naïve transient time‐stepping from zero initial conditions. We use a recently developed high‐order Eulerian–Lagrangian finite element method in combination with an energy‐preserving dynamic contact algorithm in order to solve the coupled electro‐mechanical boundary value problem. The nonlinear coupled equations are evolved by means of an operator split of the mechanical and electrical problem with an explicit as well as implicit approach. The presented benchmark examples include the first three fundamental modes of a vibrating nanotube, as well as a micro‐electro‐mechanical disk resonator in dynamic steady contact. For the examples discussed, we observe power law computational speed‐ups of the form S  = 0.6·ξ  ^? 0.8, where ξ is the linear damping ratio of the corresponding resonance frequency. Copyright © 2016 John Wiley & Sons, Ltd.     

Solving the advection–diffusion equation on unstructured meshes with discontinuous/mixed finite elements and a local time stepping procedure

Explicit schemes are known to provide less numerical diffusion in solving the advection–diffusion equation, especially for advection‐dominated problems. Traditional explicit schemes use fixed time steps restricted by the global CFL condition in order to guarantee stability. This is known to slow down the computation especially for heterogeneous domains and/or unstructured meshes. To avoid this problem, local time stepping procedures where the time step is allowed to vary spatially in order to satisfy a local CFL condition have been developed. In this paper, a local time stepping approach is used with a numerical model based on discontinuous Galerkin/mixed finite element methods to solve the advection–diffusion equation. The developments are detailed for general unstructured triangular meshes. Numerical experiments are performed to show the efficiency of the numerical model for the simulation of (i) the transport of a solute on highly unstructured meshes and (ii) density‐driven flow, where the velocity field changes at each time step. The model gives stable results with significant reduction of the computational cost especially for the non‐linear problem. Moreover, numerical diffusion is also reduced for highly advective problems. Copyright © 2009 John Wiley & Sons, Ltd.     

On time integration in the XFEM

The extended finite element method (XFEM) is often used in applications that involve moving interfaces. Examples are the propagation of cracks or the movement of interfaces in two‐phase problems. This work focuses on time integration in the XFEM. The performance of the discontinuous Galerkin method in time (space–time finite elements (FEs)) and time‐stepping schemes are analyzed by convergence studies for different model problems. It is shown that space–time FE achieve optimal convergence rates. Special care is required for time stepping in the XFEM due to the time dependence of the enrichment functions. In each time step, the enrichment functions have to be evaluated at different time levels. This has important consequences in the quadrature used for the integration of the weak form. A time‐stepping scheme that leads to optimal or only slightly sub‐optimal convergence rates is systematically constructed in this work. Copyright © 2009 John Wiley & Sons, Ltd.