A finite element method for a singularly perturbed boundary value problem

Summary We examine the problem:u+a(x)u–b(x)u=f(x) for 0<x<1,a(x)>0,b(x)>, ² = 4>0,a, b andf inC ² [0, 1], in (0, 1],u(0) andu(1) given. Using finite elements and a discretized Green's function, we show that the El-Mistikawy and Werle difference scheme on an equidistant mesh of widthh is uniformly second order accurate for this problem (i.e., the nodal errors are bounded byCh ², whereC is independent ofh and ). With a natural choice of trial functions, uniform first order accuracy is obtained in theL (0, 1) norm. On choosing piecewise linear trial functions (hat functions), uniform first order accuracy is obtained in theL ¹ (0, 1) norm.     

A method for solving the linear boundary value problem for an integro-differential equation

A method for solving the linear boundary value problem for an integro-differential equation is proposed that is based on interval partition and the introduction of additional parameters. Necessary and sufficient conditions for the solvability of the problem are obtained.     

Discontinuous Galerkin finite element method for a nonlinear boundary value problem

In this paper, we investigate the a priori and a posteriori error estimates for the discontinuous Galerkin finite element approximation to a regularization version of the variational inequality of the second kind. We show the optimal error estimates in the DG-norm （stronger than the H1 norm） and the L2 norm, respectively. Furthermore, some residual-based a posteriori error estimators are established which provide global upper bounds and local lower bounds on the discretization error. These a posteriori analysis results can be applied to develop the adaptive DG methods.     

A method for solving an exterior three-dimensional boundary value problem for the Laplace equation

We develop and experimentally study the algorithms for solving three-dimensionalmixed boundary value problems for the Laplace equation in unbounded domains. These algorithms are based on the combined use of the finite elementmethod and an integral representation of the solution in a homogeneous space. The proposed approach consists in the use of the Schwarz alternating method with consecutive solution of the interior and exterior boundary value problems in the intersecting subdomains on whose adjoining boundaries the iterated interface conditions are imposed. The convergence of the iterative method is proved. The convergence rate of the iterative process is studied analytically in the case when the subdomains are spherical layers with the known exact representations of all consecutive approximations. In this model case, the influence of the algorithm parameters on the method efficiency is analyzed. The approach under study is implemented for solving a problem with a sophisticated configuration of boundaries while using a high precision finite elementmethod to solve the interior boundary value problems. The convergence rate of the iterations and the achieved accuracy of the computations are illustrated with some numerical experiments.     

A numerical method for solving one nonlocal boundary value problem for a third-order hyperbolic equation

A nonlocal boundary value problem for a third-order hyperbolic equation with variable coefficients is considered in the one- and multidimensional cases. A priori estimates for the nonlocal problem are obtained in the differential and difference formulations. The estimates imply the stability of the solution with respect to the initial data and the right-hand side on a layer and the convergence of the difference solution to the solution of the differential problem.     

The finite element method for a boundary value problem with strong singularity

The existence and uniqueness of the R_ν-generalized solution for the third-boundary-value problem and the non-self-adjoint second-order elliptic equation with strong singularity are established. We construct a finite element method with a basis containing singular functions. The rate of convergence of the approximate solution to the R_ν-generalized solution in the norm of the Sobolev weighted space is established and, finally, results of numerical experiments are presented.     

Iterative method for solving the Neumann boundary value problem for biharmonic type equation

The solution of boundary value problems (BVP) for fourth order differential equations by their reduction to BVP for second order equations, with the aim to use the achievements for the latter ones attracts attention from many researchers. In this paper, using the technique developed by ourselves in recent works, we construct iterative method for the Neumann BVP for biharmonic type equation. The convergence rate of the method is proved and some numerical experiments are performed for testing it in dependence on the choice of an iterative parameter.     

The initial boundary value problem for a nonlocal singularly perturbed reaction-diffusion equation

The initial boundary value problem for a nonlinear singularly perturbed integro-parabolic equation is examined. An asymptotic expansion of the solution to this problem containing the temporal, spatial and corner boundary layers is constructed. The existence and local uniqueness of the solution is justified by using the asymptotic method of differential inequalities.     

A finite element method with a large mesh-width for a stiff two-point boundary value problem

We describe a finite element method for computation of numerical approximations of the solution of the second order singularly perturbed two-point boundary value problem on [?1, 1]

$\text{? u″ + pu′ = f, u(?1) = u(1) = 0, 0 < ? ∠ 1, (′ =}\text{d}\text{dx}\text{)}$

On a quasi-uniform mesh we construct exponentially fitted trial spaces which consist of piece-wise polynomials and of exponentials which fit locally to the singular solution of the equation or its adjoint. We discretise the Galerkin form for the boundary problem using such exponentially fitted trial spaces. We derive rigorous bounds for the error of discretisation with respect to the energy norm and we obtain superconvergence at the mesh-points, the error depending on ?, the mesh-width and the degree of the piece-wise polynomials.     

A method of solving a nonlinear boundary value problem with a parameter for a loaded differential equation

A nonlinear loaded differential equation with a parameter on a finite interval is studied. The interval is partitioned by the load points, at which the values of the solution to the equation are set as additional parameters. A nonlinear boundary value problem for the considered equation is reduced to a nonlinear multipoint boundary value problem for the system of nonlinear ordinary differential equations with parameters. For fixed parameters, we obtain the Cauchy problems for ordinary differential equations on the subintervals. Substituting the values of the solutions to these problems into the boundary condition and continuity conditions at the partition points, we compose a system of nonlinear algebraic equations in parameters. A method of solving the boundary value problem with a parameter is proposed. The method is based on finding the solution to the system of nonlinear algebraic equations composed.     

The p-version of the finite element method for a singularly perturbed boundary value problem

The p-version of the finite element method is applied to solve the singularly perturbed two-point boundary value problem with or without turning point. With the special choice of mesh points, global error estimates are derived. In some cases, the exponential rate of convergence is obtained. Some numerical results are given to show the performance of the proposed method.     

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in an L-shaped domain

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in an L-shaped domain is considered for when the solution has singularities at the corners of the domain. The densification of the Shishkin mesh near the inner corner where different boundary conditions meet is such that the solution obtained by the classical five-point difference scheme converges to the solution of the initial problem in the mesh norm L _?? ^h uniformly with respect to the small parameter with almost second order, i.e., as a smooth solution. Numerical analysis confirms the theoretical result.     

A mixed finite element method for solving the nonstationary stokes equation

Summary This paper deals with a mixed finite element method for approximating a fourth order initial value problem arising from the nonstationary Stokes problem. For piecewise linear shape functions error estimates are given with convergence rates similar to the elliptic case. Some numerical computations will illustrate the theoretical results.     

Iterative method for solving a nonlinear boundary value problem

In this paper, a boundary value problem for a nonlinear second-order ordinary differential equation is studied. By means of the maximum principle we established the existence and the uniqueness of a solution of the problem. Then for finding the solution an iterative method is proposed. It is proved that this method converges much faster than the Picar successive approximations and in a particular case it gives two-sided monotone approximations to the exact solution of the problem. Finally, some illustrative examples are considered to confirm the efficiency of the method.     

Bilateral difference method for solving the boundary value problem for an ordinary differential equation

A method is proposed for calculating the bilateral approximations of the solution of the boundary value problem on [0, 1] for the equation y+p(x)y-q(x)y=f(x) and the derivative of the solution having the maximum deviation O(h² (h)+h³) on {kh} _k ^N =0, where(t) is the sum of the continuity moduli of the functions p, q,f, on the set of points {kh} _k ^N =0, h=1/N by means of O(N) operations. The data obtained for fairly smooth p, q,f allow interpolation to be used for calculating the bilateral approximations of the solution and its higher derivatives having the maximum deviation O(h³) on [0, 1].Translated from Matematicheskie Zametkii, Vol. 11, No. 4, pp. 421–430, April, 1972.     

A least-squares finite element method for solving the polygonal-line arc-scattering problem

This paper is concerned with the scattering problem of a polygonal-line arc. We solve this polygonal-line arc-scattering problem by a least-squares finite element method. In the method, Fourier–Bessel functions is used to capture the singularities around tips and corners. A combination of fundamental solutions is used to represent the scattered field towards infinity. We also analyse the convergence and give an error estimate of the method. Numerical experiments are also presented to show the effectiveness of our method.     

A weighted adaptive least-squares finite element method for the Poisson-Boltzmann equation

The finite element methodology has become a standard framework for approximating the solution to the Poisson-Boltzmann equation in many biological applications. In this article, we examine the numerical efficacy of least-squares finite element methods for the linearized form of the equations. In particular, we highlight the utility of a first-order form, noting optimality, control of the flux variables, and flexibility in the formulation, including the choice of elements. We explore the impact of weighting and the choice of elements on conditioning and adaptive refinement. In a series of numerical experiments, we compare the finite element methods when applied to the problem of computing the solvation free energy for realistic molecules of varying size.