The finite element method for parabolic equations

Summary In this first of two papers, computable a posteriori estimates of the space discretization error in the finite element method of lines solution of parabolic equations are analyzed for time-independent space meshes. The effectiveness of the error estimator is related to conditions on the solution regularity, mesh family type, and asymptotic range for the mesh size. For clarity the results are limited to a model problem in which piecewise linear elements in one space dimension are used. The results extend straight-forwardly to systems of equations and higher order elements in one space dimension, while the higher dimensional case requires additional considerations. The theory presented here provides the basis for the analysis and adaptive construction of time-dependent space meshes, which is the subject of the second paper. Computational results show that the approach is practically very effective and suggest that it can be used for solving more general problems.The work was partially supported by ONR Contract N00014-77-C-0623     

A priori error estimates for a coupled finite element method and mixed finite element method for a fluid-solid interaction problem

This paper presents a heterogeneous finite element method fora fluid–solid interaction problem. The method, which combinesa standard finite element discretization in the fluid regionand a mixed finite element discretization in the solid region,allows the use of different meshes in fluid and solid regions.Both semi-discrete and fully discrete approximations are formulatedand analysed. Optimal order a priori error estimates in theenergy norm are shown. The main difficulty in the analysis iscaused by the two interface conditions which describe the interactionbetween the fluid and the solid. This is overcome by explicitlybuilding one of the interface conditions into the finite elementspaces. Iterative substructuring algorithms are also proposedfor effectively solving the discrete finite element equations.     

半线性抛物方程的时空有限元方法

讨论半线性抛物方程的连续Galerkin时空有限元方法,利用有限元方法和有限差分方法相结合的技巧,证明了弱解的存在唯一性,给出了时间最大模,空间L~2模,即L~∞(L~2)模误差估计.并给出数值算例证明了连续时空有限元方法对于半线性抛物方程的有效性.     

Discontinuous finite volume element method for parabolic problems

In this article, we consider the semidiscrete and the backward Euler fully discrete discontinuous finite volume element methods for the second‐order parabolic problems and obtain the optimal order error estimates in a mesh dependent norm and in the L²‐norm. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

A compact finite element method for elastic bodies

A nonconforming finite element method is described for treating linear equilibrium problems, and a convergence proof showing second order accuracy is given. The close relationship to a related compact finite difference scheme due to Phillips and Rose [1] is examined. A Condensation technique is shown to preserve the compactness property and suggests an approach to a certain type of homogenization.     

A finite element method for the Tricomi problem

Summary We consider the numerical solution of the Tricomi problem. Using a weak formulation based on different spaces of test and trial functions, we construct a new Galerkin procedure for the Tricomi problem. Existence, uniqueness, and uniform stability of the approximate solution is proven, and a priori error bounds are given.Research supported in part by the Department of Energy under contract DOE E(40-1)3443     

A finite element method for cohesive crack modelling

The present contribution is concerned with the computational modelling of cohesive cracks, whereby the discontinuity is not limited to interelement boundaries, but is allowed to propagate freely through the elements. Inelastic material behaviour is described by a discrete constitutive law, formulated in terms of tractions and displacements at the surface. Details on the implementation and numerical examples are given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Augmented hybrid finite element method for domain decomposition

We present an Augmented Hybrid Finite Element Method for Domain Decompositon. In this method, finite element approximations are defined independently on each subdomain and do not match at interface. This dows the user to mda local change of design, of meshes on one aubdomain without modifying other subdomains. Optimal reaults are obtained for a second-order model problem.     

Pre- and post-processing for the finite element method

The finite element method provides a powerful procedure to mathematically model physical phenomena. The technique is numerically formulated and is effectively used on a broad range of computers. The method has increased in both popularity and functionality with the development of user friendly pre- and post-processing software. Pre-processing software is used to create the model, generate an appropriate finite element grid, apply the appropriate boundary conditions, and view the total model. Post-processing provides visualization of the computed results. This paper addresses the pertinent issues of pre- and post-processing for finite element analysis. It reviews the capabilities that are provided by pre- and post-processors and suggests enhancements and new features that will likely be developed in the near future.     

Parallel finite element splitting-up method for parabolic problems

An efficient method for solving parabolic systems is presented. The proposed method is based on the splitting-up principle in which the problem is reduced to a series of independent 1D problems. This enables it to be used with parallel processors. We can solve multidimensional problems by applying only the 1D method and consequently avoid the difficulties in constructing a finite element space for multidimensional problems. The method is suitable for general domains as well as rectangular domains. Every 1D subproblem is solved by applying cubic B-splines. Several numerical examples are presented.     

A finite element method for the Sivashinsky equation

The Sivashinsky equation is a nonlinear evolutionary equation of fourth order in space. In this paper we have analyzed a semidiscrete finite element method and completely discrete scheme based on the backward Euler method and Crank–Nicolson–Galerkin scheme. A linearized backward Euler method have been developed and error bounds are derived for an L² projection.     

Nonconforming finite element method for stochastic Stokes equations

In this paper, we study nonconforming finite element method for stochastic Stokes equation driven by white noise. We apply “green function framework” and standard duality technique to study the error estimate for velocity in L²$L^2$-norm and for pressure in H^-1$H^{-1}$-norm. Finally, numerical experiment proves our theoretical results.     

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     

NONCONFORMING FINITE ELEMENT PENALTY METHOD FOR STOKES EQUATION

a special penalty method is presented to improve the accuracy of the standard penaltymethod (or solving Stokes equation with nonconforming finite element, It is shown that thismethod with a larger penalty parameter can achieve the same accuracy as the staodaxd methodwith a smaller penalty parameter. The convergence rate of the standard method is just hall order of this penalty method when using the same penalty parameter, while the extrapolationmethod proposed by Faik et al can not yield so high accuracy of convergence. At last, we alsoget the super-convergence estimates for total flux.     

A finite element method for contact/impact

Ideas from the analysis of differential-algebraic equations are applied to the numerical solution of frictionless contact/impact problems in solid mechanics. Index-one and two formulations for dynamic contact–impact within the context of the finite element method are derived. The resulting equations are shown to stabilize the kinematic fields at the contact interface, at the expense of a small energy loss, which is shown to decrease consistently with mesh refinement. This energy dissipation is shown to be necessary for the establishment of persistent contact. A Newmark-type time integration scheme is derived from the proposed formulation, and shown to yield excellent results in modeling the transition to contact/impact.     

A time finite element method for structural dynamics

A time (Galerkin) finite element method (time FEM) for structural dynamics is proposed in this paper. The key lies in a variational formulation that is well-posed and equivalent to the conventional strong form of governing equations of structural dynamics. Based on the variational formulation, a time finite element formulation is naturally established and its convergence property is easily derived through an a priori error analysis. Technical details on practical implementation of the time FEM are presented. Numerical examples are studied to verify the proposed time FEM.     

Higher-order conformal decomposition finite element method

Numerical simulations of structures using higher-order finite elements is still a challenging task, in particular for domains with curved boundaries. A new higher-order accurate approach is proposed, combining the advantages of the classical p-FEM with embedded domain methods. Boundaries and/or interfaces are described implicitly using the level set method. In the elements cut by the zero level set, an automatic decomposition into interface aligned, i. e. conforming sub-elements is realized. Transfinite mappings are utilized to construct higher-order sub-elements by mappings of reference elements to the two sides of the boundary or interface. It is shown that although the resulting sub-elements are not always well-shaped, optimal convergence rates are possible. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A finite element time relaxation method

We discuss a finite element time-relaxation method for high Reynolds number flows. The method uses local projections on polynomials defined on macroelements of each pair of two elements sharing a face. We prove that this method shares the optimal stability and convergence properties of the continuous interior penalty (CIP) method. We give the formulation both for the scalar convection–diffusion equation and the time-dependent incompressible Euler equations and the associated convergence results. This note finishes with some numerical illustrations.