Singularities of solutions to linear,second order analytic elliptic equations in two independent variables I. The completely regular boundary

Two methods are described for the a priori location of singularities of solutions to exterior boundary value problems. One uses an expansion for the solution in a circle centered on a regular exterior point P. A singularity lies on the circle of convergence. The envelope of these circles, generated as P makes a circuit about the closed boundary, circumscribes the singularities. The radius of convergence depends on singularities of the solution u(s) and its normal derivative v(s) on the boundary. The second method employs complex characteristics to relate singularities of the boundary data to real singularities of the solution. Integral equations connecting (y), v(s) and the analytic boundary condition are used to continue the data into the complex s-plane and to locate their singularities. Explicit solution of the integral equations is unnecessary; some nonlinear boundary conditions can be handled.     

The hp-Version of the Boundary Element Method in ℝ3 The Basic Approximation Results

This paper deals with the basic approximation properties of the h–p version of the boundary element method (BEM) in ℝ³. We extend the results on the exponential convergence of the h–p version of the boundary element method on geometric meshes from problems in polygonal domains to problems in polyhedral domains. In 2D elliptic boundary value problems the solutions have only corner singularities whereas in 3D problems they contain additional edge and corner-edge singularities. The solutions of the corresponding boundary integral equations inherit those singularities. The detailed investigations in our analysis take care of the various types of those singularities. While edge singularities can be analysed using standard one-dimensional approximation results the corner-edge singularities demand a new analysis. © 1997 by B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

Boundary singularities with a simple decomposition

We define the decomposition of a boundary singularity as a pair (a singularity in the ambient space together with a singularity of the restriction to the boundary). We prove that the Lagrange transform is an involution on the set of boundary singularities that interchanges the singularities that occur in the decomposition of a boundary singularity. We classify the boundary singularities for which both of these singularities are simple. Bibliography: 8 titles.Translated fromTrudy Seminara imeni I. G. Petrovskogo, No. 15, pp. 55–69, 1991.     

Generalized boundary values of the solutions of semilinear elliptic equations from weighted functional spaces

In weighted C-spaces, we establish the solvability of a boundary-value problem for a semilinear elliptic equation of order 2m in a bounded domain with generalized functions given on its boundary, strong power singularities at some points of the boundary, and finite orders of singularities on the entire boundary. The behavior of the solution near the boundary of the domain is analyzed.     

Realizations of Differential Operators on Conic Manifolds with Boundary

We study the closed extensions (realizations) of differential operators subject to homogeneous boundary conditions on weighted L _p -Sobolev spaces over a manifold with boundary and conical singularities. Under natural ellipticity conditions we determine the domains of the minimal and the maximal extension. We show that both are Fredholm operators and give a formula for the relative index. Mathematics Subject Classifications (2000): Primary 58J32; Secondary 35G70, 35S15.     

Finite element solutions for three‐dimensional elliptic boundary value problems on unbounded domains

The finite element (FE) solutions of a general elliptic equation ?div([a_ij] ??u) + u = f in an exterior domain Ω, which is the complement of a bounded subset of R ³, is considered. The most common approach to deal with exterior domain problems is truncating an unbounded subdomain Ω_∞, so that the remaining part Ω_B = Ω\Ω _∞ is bounded, and imposing an artificial boundary condition on the resulted artificial boundary Γ_a = Ω _∞ ∩ Ω _B. In this article, instead of discarding an unbounded subdomain Ω_∞ and introducing an artificial boundary condition, the unbounded domain is mapped to a unit ball by an auxiliary mapping. Then, a similar technique to the method of auxiliary mapping, introduced by Babu?ka and Oh for handling the domain singularities, is applied to obtain an accurate FE solution of this problem at low cost. This method thus does have neither artificial boundary nor any restrictions on f. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

The hp‐version of the BEM with quasi‐uniform meshes for a three‐dimensional crack problem: The case of a smooth crack having smooth boundary curve

The article considers a three‐dimensional crack problem in linear elasticity with Dirichlet boundary conditions. The crack in this model problem is assumed to be a smooth open surface with smooth boundary curve. The hp‐version of the boundary element method with weakly singular operator is applied to approximate the unknown jump of the traction which is not L₂‐regular due to strong edge singularities. Assuming quasi‐uniform meshes and uniform distributions of polynomial degrees, we prove an a priori error estimate in the energy norm. The estimate gives an upper bound for the error in terms of the mesh size h and the polynomial degree p. It is optimal in h for any given data and quasi‐optimal in p for sufficiently smooth data. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

On the Solution of the Stokes Equation on Regions with Corners

The detailed behavior of solutions to Stokes equations on regions with corners has been historically difficult to characterize. The solutions to Stokes equations on regions with corners are known to develop singularities in the vicinity of corners; in particular, the solutions are known to have infinite oscillations along almost every ray that meet at the corner. While the nature of singularities for the differential equation have been analyzed in great detail, very little is known about the nature of singularities for the corresponding integral equations. In this paper, we observe that, when the Stokes equation is formulated as a boundary integral equation, the solutions are representable by rapidly convergent series of the form , where t is the distance from the corner and the parameters μ_j, β_j are real, and are determined via an explicit formula depending on the angle at the corner. In addition to being analytically perspicuous, these representations lend themselves to the construction of highly accurate and efficient numerical discretizations, significantly reducing the number of degrees of freedom required for the solution of the corresponding integral equations. The results are illustrated by several numerical examples. © 2020 Wiley Periodicals LLC     

Some applications of a galerkin-collocation method for boundary integral equations of the first kind

Here we apply the boundary integral method to several plane interior and exterior boundary value problems from conformal mapping, elasticity and fluid dynamics. These are reduced to equivalent boundary integral equations on the boundary curve which are Fredholm integral equations of the first kind having kernels with logarithmic singularities and defining strongly elliptic pseudodifferential operators of order - 1 which provide certain coercivity properties. The boundary integral equations are approximated by Galerkin's method using B-splines on the boundary curve in connection with an appropriate numerical quadrature, which yields a modified collocation scheme. We present a complete asymptotic error analysis for the fully discretized numerical equations which is based on superapproximation results for Galerkin's method, on consistency estimates and stability properties in connection with the illposedness of the first kind equations in L₂. We also present computational results of several numerical experiments revealing accuracy, efficiency and an amazing asymptotical agreement of the numerical with the theoretical errors. The method is used for computations of conformal mappings, exterior Stokes flows and slow viscous flows past elliptic obstacles.     

Correlation between pole location and asymptotic behavior for Painlevé I solutions

We extend the technique of asymptotic series matching to exponential asymptotics expansions (transseries) and show by using asymptotic information that the extension provides a method of finding singularities of solutions of nonlinear differential equations. This transasymptotic matching method is applied to Painlevé's first equation, P1. The solutions of P1 that are bounded in some direction towards infinity can be expressed as series of functions obtained by generalized Borel summation of formal transseries solutions; the series converge in a neighborhood of infinity. We prove (under certain restrictions) that the boundary of the region of convergence contains actual poles of the associated solution. As a consequence, the position of these exterior poles is derived from asymptotic data. In particular, we prove that the location of the outermost pole x_p(C) on ℝ⁺ of a solution is monotonic in a parameter C describing its asymptotics on anti‐Stokes lines and obtain rigorous bounds for x_p(C). We also derive the behavior of x_p(C) for large C ∈ ℂ. The appendix gives a detailed classical proof that the only singularities of solutions of P1 are poles. © 1999 John Wiley & Sons, Inc.     

Variation homological diagrams,¶and corner singularities

We describe the general homological framework (the variation arrays and variation homological diagrams) in which can be studied hypersurface isolated singularities as well as boundary singularities and corner singularities from the point of view of duality. We then show that any corner singularity is extension, in a sense which is defined, of the corner singularities of less dimension on which it is built. This framework is also used to rewrite Thom–Sebastiani type properties for isolated singularities and to establish them for boundary singularities. Received: 27 June 2000 / Revised version: 18 October 2000     

On the string-theoretic Euler number of a class of absolutely isolated singularities

An explicit computation of the so-called string-theoretic E-function of a normal complex variety X with at most log-terminal singularities can be achieved by constructing one snc-desingularization of X, accompanied with the intersection graph of the exceptional prime divisors, and with the precise knowledge of their structure. In the present paper, it is shown that this is feasible for the case in which X is the underlying space of a class of absolutely isolated singularities (including both usual ? _n -singularities and Fermat singularities of arbitrary dimension). As byproduct of the exact evaluation of lim, for this class of singularities, one gets counterexamples to a conjecture of Batyrev concerning the boundedness of the string-theoretic index. Finally, the string-theoretic Euler number is also computed for global complete intersections in ℙ_ℂ ^N with prescribed singularities of the above type. Received: 2 January 2001 / Revised version: 22 May 2001     

On the boundary behavior of solutions of the Beltrami equations

We show that any homeomorphic solution of the Beltrami equation v from the Sobolev class W _loc^1,1 is a so-called lower Q-homeomorphism with Q(z) = K _μ(z), where K _μ(z) is the dilatation ratio of this equation. On this basis, we develop the theory of boundary behavior and removing of singularities of these solutions.     

Pseudo‐reflection phenomena for singularities in thin elastic shells

We consider problems of statics of thin elastic shells with hyperbolic middle surface subjected to boundary conditions ensuring the geometric rigidity of the surface. The asymptotic behaviour of the solutions when the relative thickness tends to zero is then given by the membrane approximation. It is a hyperbolic problem propagating singularities along the characteristics. We address here the reflection phenomena when the propagated singularities arrive to a boundary. As the boundary conditions are not the classical ones for a hyperbolic system, there are various cases of reflection. Roughly speaking, singularities provoked elsewhere are not reflected at all at a free boundary, whereas at a fixed (or clamped) boundary the reflected singularity is less singular than the incident one. Reflection of singularities provoked along a non‐characteristic curve C are also considered. Copyright © 2003 John Wiley & Sons, Ltd.     

Saddle singularities of complexity 1 of integrable Hamiltonian systems

Properties of saddle singularities of rank 0 and complexity 1 for integrable Hamiltonian systems are studied. An invariant (f _n -graph) solving the problem of semi-local classification of saddle singularities of rank 0 for an arbitrary complexity was constructed earlier by the author. In this paper, a more simple form of the invariant for singularities of complexity 1 is obtained and some properties of such singularities are described in algebraic terms. In addition, the paper contains a list of saddle singularities of complexity 1 for systems with three degrees of freedom.     

Penalty-combined approaches to the Ritz-Galerkin and finite element methods for singularity problems of elliptic equations

Penlty coupling techniques on an interface boundary, artificial or material, are first presented for combining the Ritz–Galerkin and finite element methods. An optimal convergence rate first is proved in the Sobolev norms. Moreover, a significant coupling strategy, L + 1 = O(|ln h|), between these two methods are derived for the Laplace equation with singularities, where L + 1 is the total number of particular solutions used in the Ritz–Galerkin method, and h is the maximal boundary length of quasiuniform elements used in the linear finite element method. Numreical experiments have been carried out for solving the benchmark model: Motz's problem. Both theoretical analysis and numreical experiments clearly display the importance of penalty-combined methods is solving elliptic equations with singularities.     

On the integral equation method for the plane mixed boundary value problem of the Laplacian

Although the plane boundary value problem for the Laplacian with given Dirichlet data on one part Γ₂ and given Neumann data on the remaining part Γ₂ of the boundary is the simplest case of mixed boundary value problems, we present several applications in classical mathematical physics. Using Green's formula the problem is converted into a system of Fredholm integral equations for the yet unknown values of the solution u on Γ₂ and the also desired values of the normal derivatie on Γ₁. One of these equations has principal part of the second kind, whereas that one of the other is of the first kind. Since any improvement of constructive methods requires higher regularity of u but, on the other hand, grad u possesses singularities at the collision points Γ₁ ∩ Γ₂ even for C^∞ data, u is decomposed into special singular terms and a regular rest. This is incorporated into the integral equations and the modified system is solved in appropriate Sobolev spaces. The solution of the system requires to solve a Fredholm equation of the first kind on the arc Γ₂ providing an improvement of regularity for the smooth part of u. Since the integral equations form a strongly elliptic system of pseudodifferential operators, the Galerkin procedure converges. Using regular finite element functions on Γ₁ and Γ₂ augmented by the special singular functions we obtain optimal order of asymptotic convergence in the norm corresponding to the energy norm of u and also superconvergence as well as high orders in smoother norms if the given data are smooth (and not the solution).     

On the boundary value problem of the biharmonic operator on domains with angular corners

The paper is concerned with boundary singularities of weak solutions of boundary value problems governed by the biharmonic operator. The presence of angular corner points or points at which the type of boundary condition changes in general causes local singularities in the solution. For that case the general theory of V. A. Kondrat'ev provides a priori estimates in weighted Sobolev norms and asymptotic singular representations for the solution which essentially depend on the zeros of certain transcendental functions. The distribution of these zeros will be analysed in detail for the biharmonic operator under several boundary conditions. This leads to sharp a priori estimates in weighted Sobolev norms where the weight function is characterized by the inner angle of the boundary corner. Such estimates for “negative” Sobolev norms are used to analyse also weakly nonlinear perturbations of the biharmonic operator as, for instance, the von Kármán model in plate bending theory and the stream function formulation of the steady state Navier-Stokes problem. It turns out that here the structure of the corner singularities is essentially the same as in the corresponding linear problem.     

Exterior Stokes flows with stick-slip boundary conditions

Steady two-dimensional creeping flows induced by line singularities in the presence of an infinitely long circular cylinder with stick-slip boundary conditions are examined. The singularities considered here include a rotlet, a potential source and a stokeslet located outside a cylinder and lying in a plane containing the cylinder axis. The general exterior boundary value problem is formulated and solved in terms of a stream function by making use of the Fourier expansion method. The solutions for various singularity driven flows in the presence of a cylinder are derived from the general results. The stream function representation of the solutions involves a definite integral whose evaluation depends on a non-dimensional slip parameter l₁\lambda_1. For extremal values, l₁ = 0\lambda_1 = 0 and l₁ = 1\lambda_1 = 1, of the slip parameter our results reduce to solutions of boundary value problems with stick (no-slip) and perfect slip conditions, respectively.¶The slip parameter influences the flow patterns significantly. The plots of streamlines in each case show interesting flow patterns. In particular, in the case of a single rotlet/stokeslet (with axis along y-direction) flows, eddies are observed for various values of l₁\lambda_1. The flow fields for a pair of singularities located on either side of the cylinder are also presented. In these flows, eddies of different sizes and shapes exist for various values of l₁\lambda_1 and the singularity locations. Plots of the fluid velocity on the surface show locations of the stagnation points on the surface of the cylinder and their dependencies on l₁\lambda_1 and singularity locations.     

Singular solutions of the Poisson equation at edges of three‐dimensional domains and their treatment with a predictor–corrector finite element method

Solutions of boundary value problems in three‐dimensional domains with edges may exhibit singularities which are known to influence both the accuracy of the finite element solutions and the rate of convergence in the error estimates. This paper considers boundary value problems for the Poisson equation on typical domains Ω ? ?³ with edge singularities and presents, on the one hand, explicit computational formulas for the flux intensity functions. On the other hand, it proposes and analyzes a nonconforming finite element method on regular meshes for the efficient treatment of the singularities. The novelty of the present method is the use of the explicit formulas for the flux intensity functions in defining a postprocessing procedure in the finite element approximation of the solution. A priori error estimates in H¹(Ω) show that the present algorithm exhibits the same rate of convergence as it is known for problems with regular solutions.