Analysis of finite element method for one-dimensional time-dependent Schrödinger equation on unbounded domain

This paper addresses the theoretical analysis of a fully discrete scheme for the one-dimensional time-dependent Schrödinger equation on unbounded domain. We first reduce the original problem into an initial-boundary value problem in a bounded domain by introducing a transparent boundary condition, then fully discretize this reduced problem by applying Crank–Nicolson scheme in time and linear or quadratic finite element approximation in space. By a rigorous analysis, this scheme has been proved to be unconditionally stable and convergent, its convergence order has also be obtained. Finally, two numerical examples are performed to show the accuracy of the scheme.     

The coupling of bem and fem for the two-dimensional viscous flow problem

In this paper we propose a numerical scheme for treating the problem of sJow viscous flow past an obstacle in the plane. This scheme is a combination of boundary element and finite element methods. By introducing an auxiliary boundary curve, we divide the region under consideration into two subregions, an inner and an outer region. In the inner region, we employ a finite element method (FEM) for solving a system of simplified field equations with proper natural boundary conditions. In the outer region, the solution is expressed in the form of a simple-layer potential with density function satisfying a system of modified integral equations of the first kind. The latter are solved by a boundary element method (BEM). Both solutions are matched on the common auxiliary boundary curve. Error estimates in suitable function spaces are derived in terms of the mesh widths as well as the small parameters, the Reynolds numbers     

Thermo-mechanical modeling of crack propagation in dynamically loaded elastomer specimens using a scaled boundary finite element approach

Recently, a scaled boundary finite element (SBFE) formulation for geometrically and physically nonlinear materials has been developed using the scaled boundary finite element method (SBFEM). The SBFE formulation has been employed to describe plane stress problems of notched and unnotched hyperelastic elastomer specimens. In this contribution, the derived SBFE formulation is extended to nonlinear time- and temperature-dependent material behavior. Subsequently, the SBFE formulation is incorporated into a crack propagation scheme to model crack propagation in cyclically loaded elastomer specimens of the so-called tear fatigue analyzer (TFA). (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of a difference scheme for the heat equation in a long strip by artificial boundary conditions

The numerical solution of the heat equation on a strip in two dimensions is considered. An artificial boundary is introduced to make the computational domain finite. On the artificial boundary, an exact boundary condition is proposed to reduce the original problem to an initial‐boundary value problem in a finite computational domain. A difference scheme is constructed by the method of reduction of order to solve the problem in the finite computational domain. It is proved that the difference scheme is uniquely solvable, unconditionally stable and convergent with the convergence order 2 in space and order 3/2 in time in an energy norm. A numerical example demonstrates the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Eigenvibrations of a bar with elastically attached load

We study the problem on the eigenvibrations of a bar with an elastically attached load. The problem is reduced to finding the eigenvalues and eigenfunctions of an ordinary secondorder differential problem with a spectral parameter nonlinearly occurring in the boundary condition at the load attachment point. We prove the existence of countably many simple positive eigenvalues of the differential problem. The problem is approximated by a grid scheme of the finite element method. We study the convergence and accuracy of the approximate solutions.     

Numerical solutions of a two-phase membrane problem

In this paper different numerical methods for a two-phase free boundary problem are discussed. In the first method a novel iterative scheme for the two-phase membrane is considered. We study the regularization method and give an a posteriori error estimate which is needed for the implementation of the regularization method. Moreover, an efficient algorithm based on the finite element method is presented. It is shown that the sequence constructed by the algorithm is monotone and converges to the solution of the given free boundary problem. These methods can be applied for the one-phase obstacle problem as well.     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

Preconditioner for a discrete Laplace operator on a condensed grid

In this paper, it is proved that the discrete Laplace operator approximating a Dirichlet boundary value problem for a Poisson equation by a finite element method with piecewise-linear functions on an evenly condensed grid that is topologically equivalent to a rectangular grid (i.e., obtained by shifting the rectangular grid nodes) is equivalent, in the range, to the operator of a 5-point difference scheme on a uniform grid.     

Coupled FEM-BEM for elastoplastic contact problems using Lagrange multipliers

The coupling of the elastoplastic finite element and elastic boundary element methods for two-dimensional frictionless contact stress analysis is presented. Interface traction matching (boundary element approach), which involves the force terms in the finite element analysis being transformed to tractions, is chosen for the coupling method. The analysis at the contact region is performed by the finite element method, and the Lagrange multiplier approach is used to apply the contact constraints. Since the analyses of elastoplastic problems are non-linear and involve iterative solution, the reduced size of the final system of equations introduced by combining the two methods is very advantageous, especially for contact problems where the nature of the problem also involves an iterative scheme.     

On Two Iterative Least-Squares Finite Element Schemes for the Incompressible Navier–Stokes Problem

This paper is mainly devoted to a comparative study of two iterative least-squares finite element schemes for solving the stationary incompressible Navier–Stokes equations with velocity boundary condition. Introducing vorticity as an additional unknown variable, we recast the Navier–Stokes problem into a first-order quasilinear velocity–vorticity–pressure system. Two Picard-type iterative least-squares finite element schemes are proposed to approximate the solution to the nonlinear first-order problem. In each iteration, we adopt the usual L ² least-squares scheme or a weighted L ² least-squares scheme to solve the corresponding Oseen problem and provide error estimates. We concentrate on two-dimensional model problems using continuous piecewise polynomial finite elements on uniform meshes for both iterative least-squares schemes. Numerical evidences show that the iterative L ² least-squares scheme is somewhat suitable for low Reynolds number flow problems, whereas for flows with relatively higher Reynolds numbers the iterative weighted L ² least-squares scheme seems to be better than the iterative L ² least-squares scheme. Numerical simulations of the two-dimensional driven cavity flow are presented to demonstrate the effectiveness of the iterative least-squares finite element approach.     

Convergence of simple adaptive Galerkin schemes based on h − h/2 error estimators

We discuss several adaptive mesh-refinement strategies based on (h − h/2)-error estimation. This class of adaptive methods is particularly popular in practise since it is problem independent and requires virtually no implementational overhead. We prove that, under the saturation assumption, these adaptive algorithms are convergent. Our framework applies not only to finite element methods, but also yields a first convergence proof for adaptive boundary element schemes. For a finite element model problem, we extend the proposed adaptive scheme and prove convergence even if the saturation assumption fails to hold in general.     

An adaptive approach for modeling a fiber-matrix composite with the FE2 method

In materials with a complicated microstructre [1], the macroscopic material behaviour is unknown. In this work a Fiber-Matrix composite is considered with elasto-plastic fibers. A homogenization of the microscale leads to the macroscopic material properties. In the present work, this is realized in the frame of a FE² formulation. It combines two nested finite element simulations. On the macroscale, the boundary value problem is modelled by finite elements, at each integration point a second finite element simulation on the microscale is employed to calculate the stress response and the material tangent modulus. One huge disadvantage of the approach is the high computational effort. Certainly, an accompanying homogenization is not necessary if the material behaves linear elastic. This motivates the present approach to deal with an adaptive scheme. An indicator, which makes use of the different boundary conditions (BC) of the BVP on microscale, is suggested. The homogenization with the Dirichlet BC overestimates the material tangent modulus whereas the Neumann BC underestimates the modulus [2]. The idea for an adaptive modeling is to use both of the BCs during the loading process of the macrostructure. Starting initially with the Neumann BC leads to an overestimation of the displacement response and thus the strain state of the boundary value problem on the macroscale. An accompanying homogenization is performed after the strain reaches a limit strain. Dirichlet BCs are employed for the accompanying homogenization. Some numerical examples demonstrate the capability of the presented method. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A full multigrid method for nonlinear eigenvalue problems

We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations. Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We prove that the computational work of this new scheme is truly optimal, the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

A Multilevel Correction Method for Steklov Eigenvalue Problem by Nonconforming Finite Element Methods

In this paper, a multilevel correction scheme is proposed to solve the Steklov eigenvalue problem by nonconforming finite element methods. With this new scheme, the accuracy of eigenpair approximations can be improved after each correction step which only needs to solve a source problem on finer finite element space and an Steklov eigenvalue problem on the coarsest finite element space. This correction scheme can increase the overall efficiency of solving eigenvalue problems by the nonconforming finite element method. Furthermore, as same as the direct eigenvalue solving by nonconforming finite element methods, this multilevel correction method can also produce the lower-bound approximations of the eigenvalues.     

边界层问题的小波—有限元解

本文将小波分析与有限元法结合起来，建立了一种小波－有限元计算格式，并用该算法计算了一个典型的边界层问题，探讨了寻找边界层位置的过程以及计算边界层区的内部解及外部解的步骤。计算结果表明，用该法寻找的边界层位置以及所求得的内部解与真实结果完全符合。     

Approximation capability of a bilinear immersed finite element space

This article discusses a bilinear immersed finite element (IFE) space for solving second‐order elliptic boundary value problems with discontinuous coefficients (interface problem). This is a nonconforming finite element space and its partition can be independent of the interface. The error estimates for the interpolation of a Sobolev function indicate that this IFE space has the usual approximation capability expected from bilinear polynomials. Numerical examples of the related finite element method are provided. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

平面弹性问题及Stokes问题的一个矩形有限元方法

In this paper,a locking-free nonconforming rectangular finite element scheme is presented for the planar elasticity problem with pure displacement boundary condition.Meanwhile,we prove that this element is also convergent for stationary Stokes problem.     

无界区域上Stokes问题的自然边界元与有限元耦合法

§1.引言 对于用有限元方法求解平面有界区域上的Stokes问题,国内外已有大量工作,例如可见[2]、[9]及其所引文献.但对无界区域上的这一问题,由于区域的无界性给有限元方法带来了困难,边界元方法及边界元与有限元的耦合法便显示其优越性.本文提出用自然边界元与有限元的耦合法求解无界区域上的Stokes问题.这一耦合法早在作者以前的工作中被应用于求解调和问题、重调和问题和平面弹性问题,但将它用于求解     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．