On the ellipticity of coupled finite element and one-equation boundary element methods for boundary value problems

In this paper we extend some recent results on the stability of the Johnson–Nédelec coupling of finite and boundary element methods in the case of boundary value problems. In Of and Steinbach (Z Angew Math Mech 93:476–484, 2013), Sayas (SIAM J Numer Anal 47:3451–3463, 2009) and Steinbach (SIAM J Numer Anal 49:1521–1531, 2011), the case of a free-space transmission problem was considered, and sufficient and necessary conditions are stated which ensure the ellipticity of the bilinear form for the coupled problem. The proof was based on considering the energies which are related to both the interior and exterior problem. In the case of boundary value problems for either interior or exterior problems, additional estimates are required to bound the energy for the solutions of related subproblems. Moreover, several techniques for the stabilization of the coupled formulations are analysed. Applications involve boundary value problems with either hard or soft inclusions, exterior boundary value problems, and macro-element techniques.     

Two-grid finite volume element method for linear and nonlinear elliptic problems

Two-grid finite volume element discretization techniques, based on two linear conforming finite element spaces on one coarse and one fine grid, are presented for the two-dimensional second-order non-selfadjoint and indefinite linear elliptic problems and the two-dimensional second-order nonlinear elliptic problems. With the proposed techniques, solving the non-selfadjoint and indefinite elliptic problem on the fine space is reduced into solving a symmetric and positive definite elliptic problem on the fine space and solving the non-selfadjoint and indefinite elliptic problem on a much smaller space; solving a nonlinear elliptic problem on the fine space is reduced into solving a linear problem on the fine space and solving the nonlinear elliptic problem on a much smaller space. Convergence estimates are derived to justify the efficiency of the proposed two-grid algorithms. A set of numerical examples are presented to confirm the estimates. The work is supported by the National Natural Science Foundation of China (Grant No: 10601045).     

An efficient boundary element method for two-dimensional transient wave propagation problems

The boundary element method (BEM) has been recognized by its unique feature of requiring neither internal cells nor their associated domain integrals in the computation. The method preserves its elegance for transient problems by means of a certain time-stepping scheme that initiates the time integration always from the initial time. Unfortunately, this time-marching scheme becomes rather difficult to apply, because the computation time and storage requirement grow dramatically with the increasing number of time steps. This paper shows that a reduction of one half of the computation time as well as the storage requirement can be achieved by an efficient truncation scheme for two-dimensional transient wave propagation problems. In particular, a guiding parameter for the determination of the truncation limit is proposed, and the overall measure of the error with respect to the truncation guide parameter is established.     

A direct discontinuous Galerkin finite element method for convection-dominated two-point boundary value problems

The direct discontinuous Galerkin (DDG) finite element method, using piecewise polynomials of degree k ≥ 1 on a Shishkin mesh, is applied to convection-dominated singularly perturbed two-point boundary value problems. Consistency, stability and convergence of order k (up to a logarithmic factor) are proved in an energy-type norm appropriate to the method and problem. The results are robust, i.e., they hold uniformly for all values of the singular perturbation parameter. Numerical experiments confirm the theoretical convergence rate.

     

A finite difference method for free boundary problems

Fornberg and Meyer-Spasche proposed some time ago a simple strategy to correct finite difference schemes in the presence of a free boundary that cuts across a Cartesian grid. We show here how this procedure can be combined with a minimax-based optimization procedure to rapidly solve a wide range of elliptic-type free boundary value problems.     

A cell boundary element method for elliptic problems

An elementary analysis on the cell boundary element (CBEM) was given by Jeon and Sheen. In this article we improve the previous results in various aspects. First of all, stability and convergence analysis on the rectangular grids are established. Moreover, error estimates are improved. Our improved analysis was possible by recasting of the CBEM in a Petrov‐Galerkin setting. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A boundary element formulation in linear piezoelectric problems

A new boundary element formulation is presented for linear piezoelectric problems under the assumption of electrodynamic quasi-static approximation. The domain, its boundary and the time region are fully discretized into a finite number of boundary-volume-time elements, and then a set of linear equations are derived from the governing integral equations. The solving scheme for the discretized linear equations in which the displacement vector, the electric potential and the corresponding fluxes are included as unknowns, is discussed.

---

Zusammenfassung Eine neue Randelement-Formulierung wird für piezoelektrische Probleme vorgeschlagen. Dabei wird eine elektrodynamische quasistatische Näherung angenommen. Das Gebiet, der Rand und der Zeitbereich werden in eine endliche Zahl von Rand-Volumen-Zeit-Elementen diskretisiert, und dann wird aus den beherrschenden Integralgleichungen eine Reihe linearer Gleichungen hergeleitet. Das Lösungsverfahren für die diskretisierten linearen Gleichungen, in denen der Verschiebungsvektor, das elektrische Potential und die zugehörigen Flüsse als Unbekannte eingehen, wird diskutiert.

     

An augmented mixed finite element method for 3D linear elasticity problems

In this paper we introduce and analyze a new augmented mixed finite element method for linear elasticity problems in 3D. Our approach is an extension of a technique developed recently for plane elasticity, which is based on the introduction of consistent terms of Galerkin least-squares type. We consider non-homogeneous and homogeneous Dirichlet boundary conditions and prove that the resulting augmented variational formulations lead to strongly coercive bilinear forms. In this way, the associated Galerkin schemes become well posed for arbitrary choices of the corresponding finite element subspaces. In particular, Raviart-Thomas spaces of order 0 for the stress tensor, continuous piecewise linear elements for the displacement, and piecewise constants for the rotation can be utilized. Moreover, we show that in this case the number of unknowns behaves approximately as 9.5 times the number of elements (tetrahedrons) of the triangulation, which is cheaper, by a factor of 3, than the classical PEERS in 3D. Several numerical results illustrating the good performance of the augmented schemes are provided.     

A finite element method scheme for boundary value problems with noncoordinated degeneration of input data

A finite element method scheme is constructed for boundary value problems with noncoordinated degeneration of input data and singularity of a solution. We look at a rate with which an approximate solution by the proposed finite element method converges toward an exact R _ν -generalized solution in the weight set W _2,ν*+β ^2+1/1 (Ω, δ), and establish estimates for the finite element approximation.     

Preconditioners for the p-version of the Galerkin method for a coupled finite element/boundary element system

We propose and analyze efficient preconditioners for solving systems of equations arising from the p-version for the finite element/boundary element coupling. The first preconditioner amounts to a block Jacobi method, whereas the second one is partly given by diagonal scaling. We use the generalized minimum residual method for the solution of the linear system. For our first preconditioner, the number of iterations of the GMRES necessary to obtain a given accuracy grows like log² p, where p is the polynomial degree of the ansatz functions. The second preconditioner, which is more easily implemented, leads to a number of iterations that behave like p log³ p. Computational results are presented to support this theory. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 47–61, 1998     

Solving parabolic Wick-stochastic boundary value problems using a finite element method

A class of linear parabolic stochastic boundary value problems of Wick-type is studied. The equations are understood in a weak sense on a suitable stochastic distribution space, and existence and uniqueness results are provided. The paper continues to discuss a numerical method for this type of problem, based on a Galerkin type of approximation. Estimates showing linear convergence in time and space are derived, and rate of convergence results for the stochastic dimension are reported.     

A multiscale cell boundary element method for elliptic problems

In this paper we introduce the multiscale cell boundary element method (MsCBE method). The method is obtained by applying the oversampling technique of the MsFEM by Hou and Wu [T.Y. Hou, X.H. Wu, A multiscale finite element method for elliptic problems in composite materials and porous media, J. Comput. Phys. 134 (1997) 169–189] to the newly developed numerical method, the cell boundary element(CBE) method by the author and his colleagues. The advantage of the MsCBE method is that it preserves flux exactly on arbitrary subdomain without needing the dual mesh. A complete H¹ convergence analysis and numerical examples confirming our analysis are presented.     

A space-time finite element method for fractional wave problems

Numerical Algorithms - In this paper, we analyze a space-time finite element method for fractional wave problems involving the time fractional derivative of order γ (1 < γ...     

An asymptotic finite element method for improvement of solutions of boundary layer problems

Summary A scheme that uses singular perturbation theory to improve the performance of existing finite element methods is presented. The proposed scheme improves the error bounds of the standard Galerkin finite element scheme by a factor of O(ⁿ⁺¹) (where is the small parameter andn is the order of the asymptotic approximation). Numerical results for linear second order O.D.E.'s are given and are compared with several other schemes.     

A finite element method for nearly incompressible elasticity problems

A finite element method is considered for dealing with nearly incompressible material. In the case of large deformations the nonlinear character of the volumetric contribution has to be taken into account. The proposed mixed method avoids volumetric locking also in this case and is robust for (with being the well-known Lamé constant). Error estimates for the -norm are crucial in the control of the nonlinear terms.

     

A finite element method for solving singular boundary-value problems

It is proved that under certain assumptions on the functions q(t) and f(t), there is one and only one function u₀(t) ∈ at which the functional

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.