Error analysis for hybrid Trefftz methods coupling traction conditions in linear elastostatics

For linear elastostatics, the Lagrange multiplier to couple the displacement (i.e., Dirichlet) condition is well known in mathematics community, but the Lagrange multiplier to couple the traction (i.e., Neumann) condition is popular for elasticity problems by the Trefftz method in engineering community, which is called the Hybrid Trefftz method (HTM). However, there has not been any analysis for these Lagrange multipliers to couple the traction condition so far. New error analysis of the HTM for elasticity problems is explored in this paper, to derive error bounds with the optimal convergence rates. Numerical experiments are reported to support this analysis. The error analysis of the HTM for linear elastostatics is the main aim of this paper. In this paper, the collocation Trefftz method (CTM) without a multiplier is also introduced, accompanied with error analysis. Numerical comparisons are made for HTM and CTM using fundamental solutions (FS) and particular solutions (PS). The error analysis and numerical computations show that the accuracy of the HTM is equivalent to that of the CTM, but the stability of the CTM is good. For elasticity and other complicated problems, the simplicity of algorithms and programming grants the CTM a remarkable advantage. More numerical comparisons show that using PS is more efficient than using FS in both HTM and CTM. However, since the optimal convergence rates are the most important criterion in evaluation of numerical methods, the global performance of the HTM is as good as that of the CTM. The comparisons of HTM and CTM using FS and PS are the next aim of this article. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Classification and refined singularity of positive solutions for nonlinear Maxwell equations arising in mesoscopic electromagnetism

This paper is a continuation of [13], where we studied the existence and other analytic properties of positive radial solutions for a system of nonlinear Maxwell equations in the plane R²${\mathbb{R}}^2$, which arises in the modelling of mesoscopic scale electromagnetic phenomena. In this paper we derive local estimates of singular positive solutions, based on which a classification theorem of general positive solutions is established. The refined singularity of general positive solutions is also investigated by employing the theory of infinite dimensional dynamical systems.     

Superconvergence of the iterated collocation methods for Hammerstein equations

In this paper, we analyse the iterated collocation method for Hammerstein equations with smooth and weakly singular kernels. The paper expands the study which began in [16] concerning the superconvergence of the iterated Galerkin method for Hammerstein equations. We obtain in this paper a similar superconvergence result for the iterated collocation method for Hammerstein equations. We also discuss the discrete collocation method for weakly singular Hammerstein equations. Some discrete collocation methods for Hammerstein equations with smooth kernels were given previously in [3, 18].     

分片常系数电导率问题的边界元配点法

用边界元方法讨论了具有分片常系数电导率方程Δ↓(γΔ↓u)=0的Dirichlet边值问题,由于方程的基本解无法显式写出,在应用通常边界元时存在很大的困难,基于这个电导率方程的解的积分表达式,导出一个在边界和交界面上的积分方程组,并讨论了这个方程组的性质,对于这个积分方程组,用配点法进行求解,且给出其误差分析．相应的数值例子证实了算法的有效性．应该指出的是本文所用的方法也适用于具有分片常系数椭圆方程的不同边界问题。     

求解二维Poisson方程的重心插值配点法

采用重心Lagrange插值配点法计算了二维Poisson方程.采用重心Lagrange插值法构造近似函数,由配点法离散Poisson方程及其边界条件.数值算例表明方法具有理论简单、计算精度高的特点.     

A Landesman-Lazer type condition for second-order differential equations with a singularity at resonance

In this paper, we are concerned with the existence of positive periodic solutions for second-order differential equations with a singularity at resonance. A Landesman-Lazer type condition is given to obtain the existence of positive periodic solutions using phase-plane analysis method and topological degree method.     

On the Divergence of Collocation Solutions in Smooth Piecewise Polynomial Spaces for Volterra Integral Equations

In this paper we present a characterization of those smooth piecewise polynomial collocation spaces that lead to divergent collocation solutions for Volterra integral equations of the second kind. The key to these results is an equivalence result between such collocation solutions and collocation solutions in slightly smoother spaces for initial-value problems for ordinary differential equations. For the latter problems Mülthei (1979/1980) established a complete divergence (and convergence) theory. Our analysis can be extended to furnish divergence results for smooth collocation solutions to Volterra integral equations of the first kind. AMS subject classification (2000) 65R20, 65L20, 65L60.Received May 2004. Accepted September 2004. Communicated by Tom Lyche.Hermann Brunner: This research was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC).     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

On certain modified boundary value problems for ordinary differential equations with irregular singularity

We consider a linear second-order differential equation with irregularly singular point at the beginning of the interval. For the corresponding homogeneous differential equation, we obtain the asymptotics of the solutions and their derivatives near the singular point. Using some modified Green functions and taking into account the asymptotics, we consider three boundary value problems with various boundary conditions (including a weighted one) at the singular point, proving theorems on the existence and uniqueness of the solutions and giving their structure. Lithuanian Mathematical Journal, Vol. 49, No. 1, 2009, pp. 109–121     

Low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions

In this paper, we consider low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions. The time discretization is based on the Euler implicit scheme, and the spatial discretization is based on the low‐order element (P₁−P₁ or P₁−P₀) for the approximation of the velocity and pressure. Moreover, some error estimates for the numerical solution of fully discrete stabilized finite element scheme are obtained. Finally, numerical experiments are performed to confirm our theoretical results.     

Collocation methods for solving multivariable integral equations of the second kind

In a recent paper (Allouch, in press) [5] on one dimensional integral equations of the second kind, we have introduced new collocation methods. These methods are based on an interpolatory projection at Gauss points onto a space of discontinuous piecewise polynomials of degree r$r$ which are inspired by Kulkarni’s methods (Kulkarni, 2003) [10], and have been shown to give a 4r+4$4r+4$ convergence for suitable smooth kernels. In this paper, these methods are extended to multi-dimensional second kind equations and are shown to have a convergence of order 2r+4$2r+4$. The size of the systems of equations that must be solved in implementing these methods remains the same as for Kulkarni’s methods. A two-grid iteration convergent method for solving the system of equations based on these new methods is also defined.     

Study on connections of the MFS, Trefftz method, indirect BIEM and invariant MFS in the three-dimensional Laplace problems containing spherical boundaries

Following the success of a study on the method of fundamental solutions using an image concept [13], we extend to solve the three-dimensional Laplace problems containing spherical boundaries by using the three approaches. The case of eccentric sphere for the Laplace problem is considered. The optimal locations for the source distribution to include the foci in the MFS are also examined by using the image concept in the 3D problems. Whether a free constant is required or not in the MFS is also studied. The error distribution is discussed after comparing with the analytical solution derived by using the bispherical coordinates. Besides, the relationship between the Trefftz bases and the singularity in the MFS for the three-dimensional Laplace problems is also addressed. It is found that one source of the MFS contains several interior and exterior Trefftz sets through a degenerate kernel. On the contrary, one single Trefftz base can be superimposed by some lumped sources in the MFS through an indirect BIEM. Based on this finding, the relationship between the fictitious boundary densities of the indirect BIEM and the singularity strength in the MFS can be constructed due to the fact that the MFS is a lumped version of an indirect BIEM.     

The Trefftz method using fundamental solutions for biharmonic equations

In this paper, the Trefftz method of fundamental solution (FS), called the method of fundamental solution (MFS), is used for biharmonic equations. The bounds of errors are derived for the MFS with Almansi’s fundamental solutions (denoted as the MAFS) in bounded simply connected domains. The exponential and polynomial convergence rates are obtained from highly and finitely smooth solutions, respectively. The stability analysis of the MAFS is also made for circular domains. Numerical experiments are carried out for both smooth and singularity problems. The numerical results coincide with the theoretical analysis made. When the particular solutions satisfying the biharmonic equation can be found, the method of particular solutions (MPS) is always superior to the MFS and the MAFS, based on numerical examples. However, if such singular particular solutions near the singular points do not exist, the local refinement of collocation nodes and the greedy adaptive techniques can be used for seeking better source points. Based on the computed results, the MFS using the greedy adaptive techniques may provide more accurate solutions for singularity problems. Moreover, the numerical solutions by the MAFS with Almansi’s FS are slightly better in accuracy and stability than those by the traditional MFS. Hence, the MAFS with the AFS is recommended for biharmonic equations due to its simplicity.     

Energy methods for stochastic differential equations

This article is devoted to the existence of strong solutions to stochastic differential equations (SDEs). Compared with Ito's theory, we relax the assumptions on the volatility term and replace the global Lipschitz continuity condition with a local Lipschitz continuity condition and a Hoelder continuity condition. In particular, our general SDE covers the Cox–Ingersoll–Ross SDE as a special case. We note that the general weak existence theory presumably extends to our general SDE (although the explicit time dependence of the drift term and the volatility term might require some extra considerations). However, avoiding weak existence theory we prove the existence of a strong solution directly using a priori estimates (the so-called energy estimates) derived from the SDE. The benefit of this approach is that the argument only requires some basic knowledge about stochastic and functional analysis. Moreover, the underlying principle has developed to become one of the cornerstones of the modern theory of partial differential equations (PDEs). In this sense, the general goal of this article is not just to establish the existence of a strong solution to the SDE under consideration but rather to introduce a new principle in the context of SDEs that has already proven to be successful in the context of PDEs.     

带非局部边界条件的热方程的LEGENDRE配置法

对带非局部边界条件的热方程的初边值问题提出了LEGENDRE配置法,并给出其半离散逼近和全离散逼近的稳定性和收敛性分析.数值试验验证了方法的有效性.     

One-bump positive solutions for the Kirchhoff equations with perturbation in

In this paper, we study the Kirchhoff equations with perturbation in . Applying the finite reduction method, we prove that the equation has one-bump positive solutions under some suitable conditions that are given in Section 1.     

Projection methods for the solution of one class of singular integrodifferential equations

In this paper we propose a computational scheme of the general projection method for the solution of singular integrodifferential equations in the theory of stream lines and thermal conductivity. We theoretically substantiate this scheme from the standpoint of the theory of approximate methods of functional analysis. We consider particular cases of the general projection method, namely, the method of moments and the collocation method.     

Projection stabilized nonconforming finite element methods for the stokes problem

In this article we study a projection‐stabilized nonconforming finite element discretization of the Stokes problem. We present a priori error analysis and give a recovery‐based a posteriori error estimator for the considered problem. Numerical results illustrate the theoretical performance of the error estimator. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 218–240, 2017