• Abstract

• Abstract

In this paper we review recent developments in the analysis of finite elementmethods for incompressible flow problems with local projection stabilization(LPS). These methods preserve the favourable stability and approximationproperties of classical residual-based stabilization (RBS) techniques but avoidthe strong coupling of velocity and pressure in the stabilizationterms. LPS-methods belong to the class of symmetric stabilization techniquesand may be characterized as variational multiscale methods. In this workwe summarize the most important a priori estimates of thisclass of stabilization schemes developed in the past 6 years.We consider the Stokes equations, the Oseen linearizationand the Navier-Stokes equations. Furthermore, we apply it tooptimal control problems with linear(ized) flow problems, since thesymmetry of the stabilization leads to the nice featurethat the operations "discretize" and"optimize" commute.

• Abstract

We consider an adaptive finite element method (AFEM) for obstacle problems associated with linear second order elliptic boundary value problems and prove a reduction in the energy norm of the discretization error which leads to R-linear convergence. This result is shown to hold up to a consistency error due to the extension of the discrete multipliers (point functionals) to H −1 and a possible mismatch between the continuous and discrete coincidence and noncoincidence sets. The AFEM is based on a residual-type error estimator consisting of element and edge residuals. The a posteriori error analysis reveals that the significant difference to the unconstrained case lies in the fact that these residuals only have to be taken into account within the discrete noncoincidence set. The proof of the error reduction property uses the reliability and the discrete local efficiency of the estimator as well as a perturbed Galerkin orthogonality. Numerical results are given illustrating the performance of the AFEM.

• Abstract

It is well understood that a good way to discretize a pointwise length constraint in partial differential equations or variational problems is to impose it at the nodes of a triangulation that defines a lowest order finite element space. This article pursues this approach and discusses the iterative solution of the resulting discrete nonlinear system of equations for a simple model problem which defines harmonic maps into spheres. An iterative scheme that is globally convergent and energy decreasing is combined with a locally rapidly convergent approximation scheme. An explicit example proves that the local approach alone may lead to ill-posed problems; numerical experiments show that it may diverge or lead to highly irregular solutions with large energy if the starting value is not chosen carefully. The combination of the global and local method defines a reliable algorithm that performs very efficiently in practice and provides numerical approximations with low energy.

• Abstract

The initial value problem for the spatially homogeneous Boltzmann equation is considered. A deterministic numerical scheme for this problem is developed by the use of the three way decomposition of the unknown function as well as of the collision integral. On this way, almost linear complexity of the algorithm is achieved. Some numerical examples are presented.

• Abstract

The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.

• Abstract

A unified a posteriori error analysis has been developed in [18, 21–23] to analyze the finite element error a posteriori under a universal roof. This paper contributes to the finite element meshes with hanging nodes which are required for local mesh-refining. The twodimensional 1−irregular triangulations into triangles and parallelograms and their combinations are considered with conforming and nonconforming finite element methods named after or by Courant, Q₁, Crouzeix-Raviart, Han, Rannacher-Turek, and others for the Poisson, Stokes and Navier-Lamé equations. The paper provides a unified a priori and a posteriori error analysis for triangulations with hanging nodes of degree ≤ 1 which are fundamental for local mesh refinement in self-adaptive finite element discretisations.

• Abstract

In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization Galerkin method. A priori error estimates are derived for the state, the adjoint state and the control. Moreover, residual type a posteriori error estimates in the L²-norm are obtained. Finally, two numerical experiments are presented to illustrate the theoretical results.

• Abstract

We propose a novel numerical approach for delay differential equations with vanishing proportional delays based on spectral methods. A Legendre-collocation method is employed to obtain highly accurate numerical approximations to the exact solution. It is proved theoretically and demonstrated numerically that the proposed method converges exponentially provided that the data in the given pantograph delay differential equation are smooth.

• Abstract

The application of a standard Galerkin finite element method for convection-diffusion problems leads to oscillations in the discrete solution, therefore stabilization seems to be necessary. We discuss several recent stabilization methods, especially its combination with a Galerkin method on layer-adapted meshes. Supercloseness results obtained allow an improvement of the discrete solution using recovery techniques.

• Abstract

In this paper, we consider the local discontinuous Galerkin method (LDG) for solving singularly perturbed convection-diffusion problems in one- and two-dimensional settings. The existence and uniqueness of the LDG solutions are verified. Numerical experiments demonstrate that it seems impossible to obtain uniform superconvergence for numerical fluxes under uniform meshes. Thanks to the implementation of two-type different anisotropic meshes, i.e., the Shishkin and an improved grade meshes, the uniform 2p + 1-order superconvergence is observed numerically for both one-dimensional and two-dimensional cases.

• Abstract

An anisotropic nonconforming finite element method is presented for a class of nonlinear Sobolev equations. The optimal error estimates and supercloseness are obtained for both semi-discrete and fully-discrete approximate schemes, which are the same as the traditional finite element methods. In addition, the global superconvergence is derived through the postprocessing technique. Numerical experiments are included to illustrate the feasibility of the proposed method.

• Abstract

In this paper, a two-scale higher-order finite elementdiscretization scheme is proposed and analyzed for a Schrodinger equation on tensor product domains. With thescheme, the solution of the eigenvalue problem on a fine grid can bereduced to an eigenvalue problem on a much coarser grid together with someeigenvalue problems on partially fine grids. It is showntheoretically and numerically that the proposed two-scale higher-order schemenot only significantly reduces the number of degrees of freedom butalso produces very accurate approximations.

• Abstract

Projection methods are efficient operator-splitting schemes toapproximatesolutions of the incompressible Navier-Stokes equations.As a major drawback, they introduce spurious layers, both in space and time.In this work, we survey convergence results forhigher order projection methods, in the presence ofonly strong solutions of the limiting problem; in particular, wehighlight concomitantdifficulties in the construction process of accurate higher orderschemes, suchas limitedregularities of the limiting solution, and a lack of accurate initialdata for thepressure.Computational experiments are includedto compare the presented schemes, and illustrate the difficulties mentioned.

• Abstract

In this survey paper we report on recent developments of the hp-version of the boundary element method (BEM). As model problems we consider weakly singular and hypersingular integral equations of the first kind on a planar, open surface. We show that the Galerkin solutions (computed with the hp-version on geometric meshes) converge exponentially fast towards the exact solutions of the integral equations. An hp-adaptive algorithm is given and the implementation of the hp-version BEM is discussed together with the choice of efficient preconditioners for the ill-conditioned boundary element stiffness matrices. We also comment on the use of the hp-version BEM for solving Signorini contact problems in linear elasticity where the contact conditions are enforced only on the discrete set of Gauss-Lobatto points. Numerical results are presented which underline the theoretical results.

• Abstract

In this article, we analyse three related preconditioned steepestdescent algorithms, which are partially popular in Hartree-Fock andKohn-Sham theory as well as invariant subspace computations, fromthe viewpoint of minimization of the corresponding functionals,constrained by orthogonality conditions. We exploit the geometry ofthe admissible manifold, i.e., the invariance with respect tounitary transformations, to reformulate the problem on the Grassmannmanifold as the admissible set. We then prove asymptotical linearconvergence of the algorithms under the condition that the Hessianof the corresponding Lagrangian is elliptic on the tangent space ofthe Grassmann manifold at the minimizer.

• Abstract

This paper presents an efficient moving mesh method to solve anonlinear singular problem with an optimal control constrainedcondition. The physical problem is governed by a new model ofturbulent flow in circular tubes proposed by Luo et al. usingPrandtl's mixing-length theory. Our algorithm is formed by an outeriterative algorithm for handling the optimal control condition andan inner adaptive mesh redistribution algorithm for solving thesingular governing equations. We discretize the nonlinear problem byusing a upwinding approach, and the resulting nonlinear equationsare solved by using the Newton-Raphson method. The mesh is generatedand the grid points are moved by using the arc-lengthequidistribution principle. The numerical results demonstrate thatproposed algorithm is effective in capturing the boundary layersassociated with the turbulent model.

• Abstract

In this paper, we consider lower order rectangular finite element methods for the singularly perturbed Stokes problem. The model problem reduces to a linear Stokes problem when the perturbation parameter is large and degenerates to a mixed formulation of Poisson’s equation as the perturbation parameter tends to zero. We propose two 2D and two 3D nonconforming rectangular finite elements, and derive robust discretization error estimates. Numerical experiments are carried out to verify the theoretical results.