Wavelet Least Squares Methods for Boundary Value Problems

We consider the number of the weak solutions for some fourth order elliptic boundary value problem with bounded nonlinear term decaying at the origin. We get a theorem, which shows the existence of the bounded solution for this problem. We obtain this result by approaching the variational method and using the generalized mountain pass theorem for the fourth order elliptic problem with bounded nonlinear term.MSC:35J30, 35J40.

Scaling function and wavelet preconditioners for second order elliptic problems

(1981). Second order elliptic boundary value problems in a domain with edges. Communications in Partial Differential Equations: Vol. 6, No. 9, pp. 1023-1042.

We consider second order elliptic boundary value problems when essential boundary conditions are enforced with the aid of Lagrange multipliers. This is combined with a fictitious domain approach into which the physical domain is embedded. The resulting saddle point problem will be discretized in terms of wavelets, resulting in an operator equation in l2. Stability of the discretization and consequently the uniform boundedness of the condition number of the finite-dimensional operator independent of the discretization is guaranteed by an appropriate LBB condition. For the iterative solution of the saddle point system, an incomplete Uzawa algorithm is employed. It can be shown that the iterative scheme combined with a nested iteration strategy is asymptotically optimal in the sense that it provides the solution up to discretization error on discretization level J in an overall amount of iterations of order O(N J ), where N J is the number of unknowns on level J. Finally, numerical results are provided.

The Variational Problems and Equivalent Forms for a Second Order Elliptic Boundary Value Problem

This paper provides an unconditional optimal convergence of a fractional-step method for solving the Boussinesq equations. In this method, the convection is treated by the Lagrange-Galerkin technique, whereas the diffusion and the incompressibility are treated by the projection method. There are lots of authors who worked on this method, and some authors gave the error estimate of this method. But, to our best knowledge, the error estimate for this method is under certain time-step restrictions. In this paper, we prove that the methods are stable almost unconditionally, i.e., when τ and h are smaller than a given constant. The basic idea of our analysis is splitting the error function into three terms, one term between the finite element solution and the projection, the other term between the projection and the time-discrete solution, the third term between the time-discrete solution and the exact solution, and giving the error estimates for each term respectively. Then, we obtain the optimal error estimates in L2 and H1-norm for the velocity and L2-norm for the pressure. In order to show the efficiency of our method, some numerical results are presented.

Second order elliptic boundary value problems which are allowed to degenerate into zero order equations are considered. The behavior of the ordinary Galerkin finite element method without special arrangements to treat singularities is studied as the problem ranges from true second order to singularly perturbed.

This part contains new pointwise error estimates for the finite element method for second order elliptic boundary value problems on smooth bounded domains in R N \mathbb {R}^{N} . In a sense to be discussed below these sharpen known quasi–optimal L ∞ L_{\infty } and W ∞ 1 W^{1}_{\infty } estimates for the error on irregular quasi–uniform meshes in that they indicate a more local dependence of the error at a point on the derivatives of the solution u u . We note that in general the higher order finite element spaces exhibit more local behavior than lower order spaces. As a consequence of these estimates new types of error expansions will be derived which are in the form of inequalities. These expansion inequalities are valid for large classes of finite elements defined on irregular grids in R N \mathbb {R}^{N} and have applications to superconvergence and extrapolation and a posteriori estimates. Part II of this series will contain local estimates applicable to non–smooth problems.

Boundedness and continuity of solutions to second order elliptic boundary value problems

We derive an optimal multilevel preconditioner for nonconforming P1 element discretizations of second order elliptic boundary value problems.

We study the question of approximability for the inverse of the FEM stiffness matrix for (scalar) second order elliptic boundary value problems by blockwise low rank matrices such as those given by the H-matrix format. We show that exponential convergence in the local block rank can be achieved. We also show that exponentially accurate LU-decompositions in the H-matrix format are possible for the stiffness matrices arising in the FEM. Unlike prior works, our analysis avoids any coupling of the block rank r and the mesh width h and also covers mixed Dirichlet-Neumann-Robin boundary conditions.

SummaryThe main result in this paper is that the solution operator for the bi-Laplace problem with zero Dirichlet boundary conditions on a bounded smooth 2d-domain can be split in a positive part and a possibly negative part which both satisfy the zero boundary condition. Moreover, the positive part contains the singularity and the negative part inherits the full regularity of the boundary. Such a splitting allows one to find a priori estimates for fourth order problems similar to the one proved via the maximum principle in second order elliptic boundary value problems. The proof depends on a careful approximative fill-up of the domain by a finite collection of limaçons. The limaçons involved are such that the Green function for the Dirichlet bi-Laplacian on each of these domains is strictly positive.

We present an a posteriori error analysis of adaptive finite element approximations of distributed control problems for second order elliptic boundary value problems under bound constraints on the control. The error analysis is based on a residual-type a posteriori error estimator that consists of edge and element residuals. Since we do not assume any regularity of the data of the problem, the error analysis further invokes data oscillations. We prove reliability and efficiency of the error estimator and provide a bulk criterion for mesh refinement that also takes into account data oscillations and is realized by a greedy algorithm. A detailed documentation of numerical results for selected test problems illustrates the convergence of the adaptive finite element method.

The original additive average Schwarz method and its variants were introduced for solving second order elliptic boundary value problems. In this paper, we extend the original idea to fourth order problems by designing a variant of the additive average Schwarz method for the plate bending problem using the nonconforming Morley finite element for the discretization. Like the original average method this method uses only nonoverlapping subdomains, and its coarse space does not need any explicit coarse mesh. An analysis of the method, and some numerical experiments to validate the theory are presented.

We consider the interior and Dirichiet problems and problems with first order boundary conditions, for a second order homogeneous elliptic partial differential operator with constant coefficients. Under natural conditions on the operators, these problems give rise to isomorphisms between the appropriate spaces with homogeneous norms. From there we obtain a priori estimates and regularity results for boundary value problems in Sobolev spaces.