Collocation Trefftz methods for the stokes equations with singularity

This study explores the boundary methods for the two‐dimensional homogeneous Stokes equations and investigates the particular solutions (PS) satisfying the Stokes equations. Smooth solutions for the Stokes equations are provided by explicit fundamental solutions (FS) and PS in this study, and singular corner solutions are also provided from linear elastostatics given in Li et al. ([Eng. Anal. Bound. Elem. 34 (2009), 533‐648, 2009). A new singularity model with an interior crack is proposed and solved by the collocation Trefftz method (CTM). The proposed method achieves highly accurate solutions with the first leading coefficient having 10 significant digits. These solutions may be used as a benchmark for testing results obtained by other numerical methods. Error bounds are derived for the CTM solutions using the PS. For a general corner, the exponent ν_k in r can only be obtained by numerical solutions of a system of nonlinear algebraic equations. Therefore, the combined method using many FS plus a few singular solutions is inevitable in most applications. For singularity problems, combining a few singular solutions with the FS is an advanced topic and is successfully implemented in Lee et al. (Eng. Anal. Bound. Elem. 24 (2010), 632–654); however, combining a few singular solutions with the smooth PS fails to converge in the first leading coefficient. As a result, the aforementioned method is not applicable to the singularity problems addressed in this article. With the help of particular and singular solutions, the hybrid Trefftz method with Lagrange multipliers can be developed for the Stokes equations. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A posteriori error estimates of discontinuous Galerkin methods for non-standard Volterra integro-differential equations

^** Email: jingtang{at}lsec.cc.ac.cn^*** Email: hermann{at}math.mun.ca In this paper we establish a posteriori error estimates forthe discontinuous Galerkin (DG) method applied to linear, semilinearand non-standard (non-linear) Volterra integro-differentialequations. We also present an analysis of the DG method withquadrature for the memory term. Numerical experiments basedon three integro-differential equations are used to illustratevarious aspects of the error analysis.     

Analytical Behaviour of Solutions of Boundary Integral Equations for a Nonsmooth Region

The Dirichlet problem for Helmholtz's equation in a domain exterior to some bounded smooth boundary in two dimensionsmay be solved by means of a combined potential of the singleand double layers. In this paper, the problem arising from allowingcorner points on the boundary is investigated. The resultingnoncompact operator is effectively split into singular and compactparts. By using the Mellin transforms, the equation can be convertedinto some Cauchy-type singular integral equations. Consequently,the singular form of the solution is found in terms of r^ßat a corner with 0>ß>1. As a first step towarddeveloping new numerical methods for the problem, one typicalexample is presented to demonstrate the slow convergence ofexisting methods without any modifications. Then the mesh-gradingtechnique designed for singular equations is successfully implementedto restore the order of convergence.     

Enlargement procedure for resolution of singularities at simple singular solutions of nonlinear equations

Summary A solution of a nonlinear equation in Hilbert spaces is said to be a simple singular solution if the Fréchet derivative at the solution has one-dimensional kernel and cokernel. In this paper we present the enlargement procedure for resolution of singularities at simple singular solutions of nonlinear equations. Once singularities are resolved, we can compute accurately the singular solution by Newton's method. Conditions for which the procedure terminates in finite steps are given. In particular, if the equation defined in ⁿ is analytic and the simple singular solution is geometrically isolated, the procedure stops in finite steps, and we obtain the enlarged problem with an isolated solution. Numerical examples are given.This research is partially supported by Grant-in-Aid for Encouragment of Young Scientist No. 60740119, the Ministry of EducationDedicated to Professor Seiiti Huzino on his 60th birthday     

Midpoint collocation for Cauchy singular integral equations

Summary A Cauchy singular integral equation on a smooth closed curve may be solved numerically using continuous piecewise linear functions and collocation at the midpoints of the underlying grid. Even if the grid is non-uniform, suboptimal rates of convergence are proved using a discrete maximum principle for a modified form of the collocation equations. The same techniques prove negative norm estimates when midpoint collocation is used to determine piecewise constant approximations to the solution of first kind equations with the logarithmic potential.This work was supported by the Australian Research Council through the program grant Numerical analysis for integrals, integral equations and boundary value problems     

On the convergence of entropy consistent schemes for lubrication type equations in multiple space dimensions

We present nonnegativity-preserving finite element schemes for a general class of thin film equations in multiple space dimensions. The equations are fourth order degenerate parabolic, and may contain singular terms of second order which are to model van der Waals interactions. A subtle discretization of the arising nonlinearities allows us to prove discrete counterparts of the essential estimates found in the continuous setting. By use of the entropy estimate, strong convergence results for discrete solutions are obtained. In particular, the limit of discrete fluxes will be identified with the flux in the continuous setting. As a by-product, first results on existence and positivity almost everywhere of solutions to equations with singular lower order terms can be established in the continuous setting.

     

Convergence in probability for perturbed stochastic integral equations

Summary In this work, one considers two stochastic integral equations indexed by some parameter and one studies the contiguity of their solutions when the parameter converges to some ₀. Two types of behaviour are described; they lead to the notion of regular and singular perturbations. The method which is used also enables a study of the rate of convergence. Applications to time discretization of equations are given.     

On multi-dimensional integral equations of convolution type with shift

The criterion of invertibility or Fredholmness of some multi-dimensional integral equations with Carleman type shifts are given. The investigation is based on some Banach space approach to equations with an involutive operator. A modified version of this approach is also presented in the paper.This approach is applied to multi-dimensional convolution type equations when the kernels may be integrable or of singular Calderon-Zygmund-Mikhlin type and shift generated by a linear transformation in the Euclidean space satisfying the generalized Carleman condition. The convolution type equations are also specially considered in the two-dimensional case in a sector on the plane symmetric with respect to one of the axes and the corresponding reflection shift. Another application deals with multi-dimensional equations with homogeneous kernels and the shift .     

Numerical analysis of a Stokes interface problem based on formulation using the characteristic function

Numerical analysis of a model Stokes interface problem with the homogeneous Dirichlet boundary condition is considered. The interface condition is interpreted as an additional singular force field to the Stokes equations using the characteristic function. The finite element method is applied after introducing a regularization of the singular source term. Consequently, the error is divided into the regularization and discretization parts which are studied separately. As a result, error estimates of order h^1/2 in H¹ × L² norm for the velocity and pressure, and of order h in L² norm for the velocity are derived. Those theoretical results are also verified by numerical examples.     

A convergent boundary integral method for three-dimensional water waves

We design a boundary integral method for time-dependent, three-dimensional, doubly periodic water waves and prove that it converges with accuracy, without restriction on amplitude. The moving surface is represented by grid points which are transported according to a computed velocity. An integral equation arising from potential theory is solved for the normal velocity. A new method is developed for the integration of singular integrals, in which the Green's function is regularized and an efficient local correction to the trapezoidal rule is computed. The sums replacing the singular integrals are treated as discrete versions of pseudodifferential operators and are shown to have mapping properties like the exact operators. The scheme is designed so that the error is governed by evolution equations which mimic the structure of the original problem, and in this way stability can be assured. The wavelike character of the exact equations of motion depends on the positivity of the operator which assigns to a function on the surface the normal derivative of its harmonic extension; similarly, the stability of the scheme depends on maintaining this property for the discrete operator. With grid points, the scheme can be implemented with essentially operations per time step.

     

Convergence behavior of Gauss–Newton’s method and extensions of the Smale point estimate theory

The notions of Lipschitz conditions with L$L$ average are introduced to the study of convergence analysis of Gauss–Newton’s method for singular systems of equations. Unified convergence criteria ensuring the convergence of Gauss–Newton’s method for one kind of singular systems of equations with constant rank derivatives are established and unified estimates of radii of convergence balls are also obtained. Applications to some special cases such as the Kantorovich type conditions, γ$\gamma$-conditions and the Smale point estimate theory are provided and some important known results are extended and/or improved.     

On the approximate solution of weakly singular integro-differential equations of Volterra type

Recently, the convergence rate of the collocation method for integral and integro-differential equations with weakly singular kernels has been studied in a series of papers [1–7]. The present paper belongs to the same series. We analyze the possibility of constructing approximate solutions of high-order accuracy on a uniform or almost uniform grid for weakly singular integro-differential equations of Volterra type.Translated from Differentsialnye Uravneniya, Vol. 40, No. 9, 2004, pp. 1271–1279.Original Russian Text Copyright © 2004 by Pedas.     

Monomial summability and doubly singular differential equations

In this work, we consider systems of differential equations that are doubly singular, i.e. that are both singularly perturbed and exhibit an irregular singular point. If the irregular singular point is at the origin, they have the form

     

A New Integral Equation Form of Integrable Reductions of the Einstein Equations

We further develop the monodromy transformation method for analyzing hyperbolic and elliptic integrable reductions of the Einstein equations. The compatibility conditions for alternative representations of solutions of the associated linear systems with a spectral parameter in terms of a pair of dressing (scattering) matrices yield a new set of linear (quasi-Fredholm) integral equations that are equivalent to the symmetry-reduced Einstein equations. In contrast to the previously derived singular integral equations constructed using conserved (nonevolving) monodromy data for fundamental solutions of the associated linear systems, the scalar kernels of the new equations involve functional parameters of a different type, the evolving (dynamic) monodromy data for scattering matrices. In the context of the Goursat problem, these data are completely determined for hyperbolic reductions by the characteristic initial data for the fields. The field components are expressed in quadratures in terms of solutions of the new integral equations.     

Collocation Methods and Their Modifications for Cauchy Singular Integral Equations on the Interval

In this paper we present polynomial collocation methods and their modi.cations for the numerical solution of Cauchy singular integral equations over the interval [-1, 1]. More precisely, the operators of the integral equations have the form with piecewise continuous coefficients a and b, and with a Jacobi weight . Using the splitting property of the singular values of the collocation methods, we obtain enough stable approximate methods to .nd the least square solution of our integral equation. Moreover, the modifications of the collocation methods enable us to compute kernel and cokernel dimensions of operators from a C*-algebra, which is generated by operators of the Cauchy singular integral equations.     

On the uniform convergence of a collocation method for a class of singular integral equations

We establish the uniform convergence of a collocation method for solving a class of singular integral equations. This method uses the Jacobi polynomials {P _n ^{(, )} } as basis elements and the zeros of a Chebyshev polynomial of the first kind as collocation points. Uniform convergence is shown to hold under the weak assumption that the kernel and the right-hand side are Hölder-continous functions. Convergence rates are also given.     

Quadrature method for singular integral equations on closed curves

Summary A method which combines quadrature with trigonometric interpolation is proposed for singular integral equations on closed curves. For the case of the circle, the present method is shown to be equivalent to the trigonometric -collocation method together with numerical quadrature for the compact term, and is shown to be stable inL ² provided the operatorA is invertible inL ². The results are extended to arbitraryC curves, to give a complete error analysis in the scale of Sobolev spacesH ^s . In the final section the case of a non-invertible operatorA is considered.     

Convergence of Step-by-step Methods for Non-linear Integro-Differential Equations

The theory of consistent step-by-step methods for solving Volterraintegral equations is extended to non-singular Volterra integro-differentialequations. It is shown that standard step-by-step algorithmsfor these more general equations are convergent. Several numericalexamples are included.     

A discontinuous Galerkin method on refined meshes for the two-dimensional time-harmonic Maxwell equations in composite materials

In this paper, a discontinuous Galerkin method for the two-dimensional time-harmonic Maxwell equations in composite materials is presented. The divergence constraint is taken into account by a regularized variational formulation and the tangential and normal jumps of the discrete solution at the element interfaces are penalized. Due to an appropriate mesh refinement near exterior and interior corners, the singular behaviour of the electromagnetic field is taken into account. Optimal error estimates in a discrete energy norm and in the L²$L^2$-norm are proved in the case where the exact solution is singular.     

On the monotonicity of some singular integral operators

The paper deals with the monotonicity of singular integral operators of the form where is the Cauchy singular integral operator on the interval (0,1) of the real axis and q is a power or logarithmic function. Under suitable assumptions, such singular integral operators are proved to be monotone and maximal monotone in spaces with power weights. Moreover, two related integral equations with weakly singular kernels of logarithmic type are studied. Copyright © 2012 John Wiley & Sons, Ltd.