A comparison of time adaptive integration methods for small and large strain viscoelasticity

In quasistatic solid mechanics the initial boundary value problem has to be solved in the space and time domain. The spatial discretization is done using finite elements. For the temporal discretization three different classes of Runge-Kutta methods are compared. These methods are diagonally implicit Runge-Kutta schemes (DIRK), linear implicit Runge-Kutta methods (Rosenbrock type methods) and half-explicit Runge-Kutta schemes (HERK). It will be shown that the application of half-explicit or linear-implicit Runge-Kutta methods can enormously reduce the computational time in particular situations. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Second-order time-integration in inelastic dynamical systems

In structural dynamics the initial boundary value problem has to be solved in the space and time domain. The spatial discretization is done using finite elements. For the temporal discretization two classes of Runge-Kutta methods and the generalized-α method are compared. The representatives for the Runge-Kutta methods are diagonally implicit Runge-Kutta schemes (DIRK) and diagonally implicit Runge-Kutta-Nyström schemes (DIRKN). (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Characterization of allA-stable methods of order 2m-4

The characterization ofA-stable methods is often considered as a very difficult task (see e.g. [1]). In recent years, simple proofs have been found for methods of orderp2m-2 (see [2], [3], [7]). In this paper, we characterize theA-acceptable approximations of orderp 2m-4 and apply the result to 12-parameter families of implicit Runge-Kutta methods.     

Optimal m-stage Runge-Kutta schemes for steady-state solutions of hyperbolic equations

A numerical scheme is developed to find optimal parameters and time step of m-stage Runge-Kutta (RK) schemes for accelerating the convergence to -steady-state solutions of hyperbolic equations. These optimal RK schemes can be applied to a spatial discretization over nonuniform grids such as Chebyshev spectral discretization. For each m given either a set of all eigenvalues or a geometric closure of all eigenvalues of the discretization matrix, a specially structured nonlinear minimax problem is formulated to find the optimal parameters and time step. It will be shown that each local solution of the minimax problem is also a global solution and therefore the obtained m-stage RK scheme is optimal. A numerical scheme based on a modified version of the projected Lagrangian method is designed to solve the nonlinear minimax problem. The scheme is generally applicable to any stage number m. Applications in solving nonsymmetric systems of linear equations are also discussed. © 1993 John Wiley & Sons, Inc.     

Cardiovascular fluid-structure interaction: A partitioned approach utilizing the p-FEM

In this paper, we propose a technique for simulating the fluid-structure interaction in blood vessels. A partitioned approach is used, which allows for an independent discretization of the fluid domain and the structural domain. We choose the finite volume method to solve the Navier-Stokes equations and the p-version of the finite element method (p-FEM) to solve the equations of geometrically nonlinear structural dynamics. The solution strategy can be seen as a first approach towards a comprehensive study of the hemodynamics in vascular substitutes with the goal of improving their long term functionality. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

p-version finite elements and finite cells for finite strain elastoplastic problems

In this paper, the p-version finite element method and its fictitious domain extension, the finite cell method, are extended to the finite strain J₂ plasticity. High-order shape functions are used for the finite element approximation of volume-preserving plastic dominated deformations. The accuracy and efficiency of p-version elements and cells in the finite plastic strain range are evaluated by the computation of two benchmark problems. It is shown that they provide locking free behavior and simplified meshing. These results are verified in comparison with the results of h-version elements in F-bar formulation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of a mixed least-squares formulation using different approximation spaces

The main goal of the present work is the comparison of the performance of a least-squares mixed finite element formulation where the solution variables (displacements and stresses) are interpolated using different approximation spaces. Basis for the formulation is a weak form resulting from the minimization of a least-squares functional, compare e.g. [1]. As suitable functions for standard interpolation polynomials of Lagrangian type are chosen. For the conforming discretization of the Sobolev space vector-valued Raviart-Thomas interpolation functions, see also [2], are used. The resulting elements are named as P_mP_k and RT_mP_k. Here m (stresses) and k (displacements) denote the approximation order of the particular interpolation function. For the comparison we consider a two-dimensional cantilever beam under plain strain conditions and small strain assumptions. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Small idempotent clones I

G. Grätzer and A. Kisielewicz devoted one section of their survey paper concerning p _n-sequences and free spectra of algebras to the topic Small idempotent clones (see Section 6 of [18]). Many authors, e.g., [8], [14, 15], [22], [25] and [29, 30] were interested in p _n-sequences of idempotent algebras with small rates of growth. In this paper we continue this topic and characterize all idempotent groupoids (G, ·) with p ₂(G, ·) 2 (see Section 7). Such groupoids appear in many papers see, e.g. [1], [4], [21], [26, 27], [25], [28, 30, 31, 32] and [34].     

Numerical homogenization of foam-like structures based on the FE2-approach

The computation of foam–like structures is still a topic of research. There are two basic approaches: the microscopic model where the foam–like structure is entirely resolved by a discretization (e.g. with Timoshenko beams) on a micro level, and the macroscopic approach which is based on a higher order continuum theory. A combination of both of them is the FE²-approach where the mechanical parameters of the macroscopic scale are obtained by solving a Dirichlet boundary value problem for a representative microstructure at each integration point. In this contribution, we present a two–dimensional geometrically nonlinear FE²-framework of first order (classical continuum theories on both scales) where the microstructures are discretized by continuum finite elements based on the p-version. The p-version elements have turned out to be highly efficient for many problems in structural mechanics. Further, a continuum–based approach affords two additional advantages: the formulation of geometrical and material nonlinearities is easier, and there is no problem when dealing with thicker beam–like structures. In our numerical example we will investigate a simple macroscopic shear test. Both the macroscopic load displacement behavior and the evolving anisotropy of the microstructures will be discussed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The approximation theory for the p-version finite element method and application to non-linear elliptic PDEs

Approximation theoretic results are obtained for approximation using continuous piecewise polynomials of degree p on meshes of triangular and quadrilateral elements. Estimates for the rate of convergence in Sobolev spaces , are given. The results are applied to estimate the rate of convergence when the p-version finite element method is used to approximate the -Laplacian. It is shown that the rate of convergence of the p-version is always at least that of the h-version (measured in terms of number of degrees of freedom used). If the solution is very smooth then the p-version attains an exponential rate of convergence. If the solution has certain types of singularity, the rate of convergence of the p-version is twice that of the h-version. The analysis generalises the work of Babuska and others to the case . In addition, the approximation theoretic results find immediate application for some types of spectral and spectral element methods. Received August 2, 1995 / Revised version received January 26, 1998     

Optimal error estimate of two linear and momentum-preserving Fourier pseudo-spectral schemes for the RLW equation

In this paper, two novel linear-implicit and momentum-preserving Fourier pseudo-spectral schemes are proposed and analyzed for the regularized long-wave equation. The numerical methods are based on the blend of the Fourier pseudo-spectral method in space and the linear-implicit Crank–Nicolson method or the leap-frog scheme in time. The two fully discrete linear schemes are shown to possess the discrete momentum conservation law, and the linear systems resulting from the schemes are proved uniquely solvable. Due to the momentum conservative property of the proposed schemes, the Fourier pseudo-spectral solution is proved to be bounded in the discrete L^∞ norm. Then by using the standard energy method, both the linear-implicit Crank–Nicolson momentum-preserving scheme and the linear-implicit leap-frog momentum-preserving scheme are shown to have the accuracy of in the discrete L^∞ norm without any restrictions on the grid ratio, where N is the number of nodes and τ is the time step size. Numerical examples are carried out to verify the correction of the theory analysis and the efficiency of the proposed schemes.     

Comparison principles applied to obstacle problems with penalties

Penalty methods form a well known technique to embed elliptic variational inequality problems into a family of variational equations (cf. [6], [13], [17]). Using the specific inverse monotonicity properties of these problems L _∞-bounds for the convergence can be derived by means of comparison solutions. Lagrange duality is applied to estimate parameters involved.

For piecewise linear finite elements applied on weakly acute triangulations in combination with mass lumping the inverse monotonicity of the obstacle problems can be transferred to its discretization. This forms the base of similar error estimations in the maximum norm for the penalty method applied to the discrete problem.

The technique of comparison solutions combined with the uniform boundedness of the Lagrange multipliers leads to decoupled convergence estimations with respect to the discretization and penalization parameters.     

求解带刚性源项标量双曲型守恒律方程的保有界WCNS格式

带刚性源项的双曲守恒律方程是很多物理问题,特别是化学反应流的数学模型.本文考虑带刚性源项的标量双曲型守恒律方程,通过时空分离的方式,发展了一类保有界的WCNS格式.对于空间离散,我们将参数化的通量限制器推广到WCNS框架,使得方程对流项离散后满足极值原理.对于时间离散,我们将半离散的WCNS改写成指数形式,采用三阶修正...     

Error estimates for the combinedh andp versions of the finite element method

Summary In theh-version of the finite element method, convergence is achieved by refining the mesh while keeping the degree of the elements fixed. On the other hand, thep-version keeps the mesh fixed and increases the degree of the elements. In this paper, we prove estimates showing the simultaneous dependence of the order of approximation on both the element degrees and the mesh. In addition, it is shown that a proper design of the mesh and distribution of element degrees lead to a better than polynomial rate of convergence with respect to the number of degrees of freedom, even in the presence of corner singularities. Numerical results comparing theh-version,p-version, and combinedh-p-version for a one dimensional problem are presented.     

Stagnation in the p-version of the finite element method

A special feature of the p-version of the finite element method for solving a differential boundary value problem stated in the form of minimizing a quadratic functional on a certain set is studied. This special feature results in approximate solutions remaining unchanged on finite numbers of increasing finite-dimensional subsets of increasing dimension, in which solutions are sought. Necessary and sufficient conditions for the existence of this feature are found, and the stagnation effect is interpreted for a specially constructed example. For the adaptive p-version of the finite element approach, a modified strategy is proposed that takes this feature into account and thus improves the reliability of the method.     

Local error estimation for singly-implicit formulas by two-step Runge-Kutta methods

In this paper it is shown that the local discretization error ofs-stage singly-implicit methods of orderp can be estimated by embedding these methods intos-stage two-step Runge-Kutta methods of orderp+1, wherep=s orp=s+1. These error estimates do not require any extra evaluations of the right hand side of the differential equations. This is in contrast with the error estimation schemes based on embedded pairs of two singly-implicit methods proposed by Burrage.The work of A. Bellen and M. Zennaro was supported by the CNR and MPI. The work of Z. Jackiewicz was supported by the CNR and by the NSF under grant DMS-8520900.     

Orthogonal Systems of Singular Functions and Numerical Treatment of Problems with Degeneration of Data

The orthogonal systems of singular functions are considered. They are applied to the error analysis of the p-version of the finite element method for elliptic problems with degeneration of data and strong singularity of solution.     

Additive Schwarz methods for thep-version finite element method

Summary In this paper, we study some additive Schwarz methods (ASM) for thep-version finite element method. We consider linear, scalar, self adjoint, second order elliptic problems and quadrilateral elements in the finite element discretization. We prove a constant bound independent of the degreep and the number of subdomainsN, for the condition number of the ASM iteration operator. This optimal result is obtained first in dimension two. It is then generalized to dimensionn and to a variant of the method on the interface. Numerical experiments confirming these results are reported. As is the case for other additive Schwarz methods, our algorithms are highly parallel and scalable.This work was supported in part by the Applied Math. Sci. Program of the U.S. Department of Energy under contract DE-FG02-88ER25053 and, in part, by the National Science Foundation under Grant NSF-CCR-9204255     

Two-Grid hp-Version DGFEMs for Strongly Monotone Second-Order Quasilinear Elliptic PDEs

In this article we develop the a priori error analysis of so-called two-grid hp-version discontinuous Galerkin finite element methods for the numerical approximation of strongly monotone second-order quasilinear partial differential equations. In this setting, the fully nonlinear problem is first approximated on a coarse finite element space V(𝒯_H, P ). The resulting ‘coarse’ numerical solution is then exploited to provide the necessary data needed to linearize the underlying discretization on the finer space V(𝒯_h, p ); thereby, only a linear system of equations is solved on the richer space V(𝒯_h, p ). Numerical experiments confirming the theoretical results are presented. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical solution of laminar incompressible generalized Newtonian fluids flow

This paper deals with the numerical solution of laminar viscous incompressible flows for generalized Newtonian fluids in the branching channel. The generalized Newtonian fluids contain Newtonian fluids, shear thickening and shear thinning non-Newtonian fluids. The mathematical model is the generalized system of Navier-Stokes equations. The finite volume method combined with an artificial compressibility method is used for spatial discretization. For time discretization the explicit multistage Runge-Kutta numerical scheme is considered. Steady state solution is achieved for t → ∞ using steady boundary conditions and followed by steady residual behavior. For unsteady solution a dual-time stepping method is considered. Numerical results for flows in two dimensional and three dimensional branching channel are presented.