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.     

Multiple Positive Solutions of Semilinear Differential Equations with Singularities

The existence of positive solutions of a second order differentialequation of the form z'+g(t)f(z)=0 (1.1) with the separated boundary conditions: z(0) – ßz'(0)= 0 and z(1)+z'(1) = 0 has proved to be important in physicsand applied mathematics. For example, the Thomas–Fermiequation, where f = z^3/2 and g = t^–1/2 (see [12, 13, 24]),so g has a singularity at 0, was developed in studies of atomicstructures (see for example, [24]) and atomic calculations [6].The separated boundary conditions are obtained from the usualThomas–Fermi boundary conditions by a change of variableand a normalization (see [22, 24]). The generalized Emden–Fowlerequation, where f = z^p, p > 0 and g is continuous (see [24,28]) arises in the fields of gas dynamics, nuclear physics,chemically reacting systems [28] and in the study of multipoletoroidal plasmas [4]. In most of these applications, the physicalinterest lies in the existence and uniqueness of positive solutions.     

Approximate Monotonicity: Theory and Applications

The ideas of value distribution for measurable functions fromR to R are applied to functions which are approximately monotonicon sets of positive measure. (For definitions see 1.) A functionp(x) is introduced, describing the local relative value distributionin the neighbourhood of a point x, and it is shown that almosteverywhere p(x) = 0 or wherever p(x) exists, implying approximate differentiability, with thefunction approximately oscillatory elsewhere. These resultsare applied to the analysis of angular boundary behaviour forHerglotz functions, where they have implications for the spectralanalysis of differential and other operators.     

Singularities and Limit Functions in Iteration of Meromorphic Functions

Let f(z) be a transcendental meromorphic function. The paperinvestigates, using the hyperbolic metric, the relation betweenthe forward orbit P(f) of singularities of f^–1 and limitfunctions of iterations of f in its Fatou components. It ismainly proved, among other things, that for a wandering domainU, all the limit functions of {fⁿ|_U} lie in the derived setof P(f) and that if f^np|_V q(n +) for a Fatou component V, theneither q is in the derived set of S_p (f) or f^p(q) = q. As applicationsof main theorems, some sufficient conditions of the non-existenceof wandering domains and Baker domains are given.     

Singularities of Jacobi series on C2 and the Poisson process equation

A classical theorem of Gabor Szego relates the singularities of real zonal harmonic expansions with those of associated analytic functions of a single complex variable. Zeev Nehari developed the counterpart for Legendre series on the C-plane by generalizing Szego's theorem. This paper function theretically identifies the singularities of analytic symmetric Jacobi series on C² with those of analytic functions on the C-plane. One feature is that information about the singularities of solutions of Solomon Bochner's Poisson process equation flow from the expansion coefficients. Others are that the Szego and Nehari theorems appear on characteristic subspaces. And, that this PDE, unlike those normally encountered in function theory, is hyperbolic in the real domain.     

A transform method for linear evolution PDEs on a finite interval

We study initial boundary value problems for linear scalar evolutionpartial differential equations, with spatial derivatives ofarbitrary order, posed on the domain {t > 0, 0 < x <L}. We show that the solution can be expressed as an integralin the complex k-plane. This integral is defined in terms ofan x-transform of the initial condition and a t-transform ofthe boundary conditions. The derivation of this integral representationrelies on the analysis of the global relation, which is an algebraicrelation defined in the complex k-plane coupling all boundaryvalues of the solution. For particular cases, such as the case of periodic boundaryconditions, or the case of boundary value problems for even-orderPDEs, it is possible to obtain directly from the global relationan alternative representation for the solution, in the formof an infinite series. We stress, however, that there existinitial boundary value problems for which the only representationis an integral which cannot be written as an infinite series.An example of such a problem is provided by the linearized versionof the KdV equation. Similarly, in general the solution of odd-orderlinear initial boundary value problems on a finite intervalcannot be expressed in terms of an infinite series.     

On the mixed Cauchy problem with data on singular conics

We consider a problem of mixed Cauchy type for certain holomorphicpartial differential operators with the principal part Q_2p(D)essentially being the (complex) Laplace operator to a power,^p. We provide inital data on a singular conic divisor givenby P = 0, where P is a homogeneous polynomial of degree 2p.We show that this problem is uniquely solvable if the polynomialP is elliptic, in a certain sense, with respect to the principalpart Q_2p(D).     

Properties of Removable Singularities for Hardy Spaces of Analytic Functions

Removable singularities for Hardy spaces H^p() = {f Hol(): |f|^p u in for some harmonic u}, 0 < p < are studied. A setE = is a weakly removable singularity for H^p(\E) if H^p(\E) Hol(), and a strongly removable singularity for H^p(\E) if H^p(\E)= H^p(). The two types of singularities coincide for compactE, and weak removability is independent of the domain . The paper looks at differences between weak and strong removability,the domain dependence of strong removability, and when removabilityis preserved under unions. In particular, a domain and a setE that is weakly removable for all H^p, but not strongly removablefor any H^p(\E), 0 < p < , are found. It is easy to show that if E is weakly removable for H^p(\E)and q > p, then E is also weakly removable for H^q(\E). Itis shown that the corresponding implication for strong removabilityholds if and only if q/p is an integer. Finally, the theory of Hardy space capacities is extended, anda comparison is made with the similar situation for weightedBergman spaces.     

Asymptotic simplification for a reaction-diffusion problem with a nonlinear boundary condition

We study non-negative solutions of the porous medium equationwith a source and a nonlinear flux boundary condition, u_t =(u^m)_xx + u^p in (0, ), x (0, T); – (u^m)_x (0, t) = u^q (0,t) for t (0, T); u (x, 0) = u₀ (x) in (0, ), where m > 1,p, q > 0 are parameters. For every fixed m we prove thatthere are two critical curves in the (p, q-plane: (i) the criticalexistence curve, separating the region where every solutionis global from the region where there exist blowing-up solutions,and (ii) the Fujita curve, separating a region of parametersin which all solutions blow up from a region where both globalin time solutions and blowing-up solutions exist. In the caseof blow up we find the blow-up rates, the blow-up sets and theblow-up profiles, showing that there is a phenomenon of asymptoticsimplification. If 2q < p + m the asymptotics are governedby the source term. On the other hand, if 2q > p + m theevolution close to blow up is ruled by the boundary flux. If2q = p + m both terms are of the same order.     

Singularities of Centre Symmetry Sets

The center symmetry set (CSS) of a smooth hypersurface S inan affine space Rⁿ is the envelope of lines joining pairs ofpoints where S has parallel tangent hyperplanes. The idea stemsfrom a definition of Janeczko, in an alternative version dueto Giblin and Holtom. For n = 2 the envelope is always real,while for n > 3 the existence of a real envelope dependson the geometry of the hypersurface. In this paper we make alocal study of the CSS, some results applying to n 5 and othersto the cases n = 2,3. The method is to construct a generatingfunction whose bifurcation set contains the CSS and possiblysome other redundant components. Focal sets of smooth hypersurfacesare a special case of the construction, but the CSS is an affineand not a euclidean invariant. Besides the familiar local formsof focal sets there are other local forms corresponding to boundarysingularities, and yet others which do not appear to have arisenelsewhere in a geometrical context. There are connections withFinsler geometry. This paper concentrates on the theory andthe proof of the local normal forms for the CSS. 2000 MathematicsSubject Classification 57R45, 58K40, 32S25, 58B20.     

Scattering Matrix for Manifolds with Conical Ends

Let M be a manifold with conical ends. (For precise definitionssee the next section; we only mention here that the cross-sectionK can have a nonempty boundary.) We study the scattering forthe Laplace operator on M. The first question that we are interestedin is the structure of the absolute scattering matrix S(s).If M is a compact perturbation of Rⁿ, then it is well-knownthat S(s) is a smooth perturbation of the antipodal map on asphere, that is, S(s)f(·)=f(–·) (mod C) On the other hand, if M is a manifold with a scattering metric(see [8] for the exact definition), it has been proved in [9]that S(s) is a Fourier integral operator on K, of order 0, associatedto the canonical diffeomorphism given by the geodesic flow atdistance . In our case it is possible to prove that S(s) isin fact equal to the wave operator at a time t = plus C terms.See Theorem 3.1 for the precise formulation. This result isnot too difficult and is obtained using only the separationof variables and the asymptotics of the Bessel functions. Our second result is deeper and concerns the scattering phasep(s) (the logarithm of the determinant of the (relative) scatteringmatrix).     

Triple-deck Flows over Wavy Surfaces

Sinusoidally perturbed laminar flow over a flat plate is considered,in which the amplitude and wave number bear the triple-deckrelation to the Reynolds number R. To find the resulting flowfield and in particular the surface shear stress and pressureit is necessary to solve the non-linear lower-deck equation.Analytic solutions of the linearized equations are obtainedand the behaviour of the stresses near the leading edge of thewaviness is examined. As regions further downstream are consideredit is found that the linearized results predict the same surfacestress phase shifts as those derived by Benjamin (1959) forboundary-layer flows over wavy surfaces of amplitude and wavenumberindependent of R. For larger amplitudes, the full non-linearequations are solved numerically via a modification of the spectralmethod of Burggraf & Duck (1981). This modification is shownto significantly increase both the speed and accuracy of theoriginal method for lower-deck problems in general. Specialattention is paid to the problems associated with the semi-infiniteextent of the surface perturbation. The results of these non-linearcalculations of the stresses are examined to determine the effectsof increasing amplitude h and of increasing distance downstream.It is found that the surface stress extrema are phase-shifteddownstream with respect to the surface perturbation both asthe amplitude is increased and as regions further downstreamof the leading edge are considered. Moreover, the stress profilesbecome considerably distorted from sinusoidal as h increases.It is found that separation occurs at higher values of h thanpredicted by linear theory. Finally, differences between resultsfor positive and negative h are examined.     

On the Number of Singularities in Generic Deformations of Map Germs

Let f:Cⁿ, 0C^p, 0 be a K-finite map germ, and let i=(i₁, ...,i_k) be a Boardman symbol such that ⁱ has codimension n in thecorresponding jet space J^k(n, p). When its iterated successorshave codimension larger than n, the paper gives a list of situationsin which the number of ⁱ points that appear in a generic deformationof f can be computed algebraically by means of Jacobian idealsof f. This list can be summarised in the following way: f musthave rank n–i₁ and, in addition, in the case p=6, f mustbe a singularity of type ^i₁,i₂.     

Riemann problem for kinematical conservation laws and geometrical features of nonlinear wavefronts

A pair of kinematical conservation laws (KCL) in a ray coordinatesystem (,t) are the basic equations governing the evolutionof a moving curve in two space dimensions. We first study elementarywave solutions and then the Riemann problem for KCL when themetric g, associated with the coordinate designating differentrays, is an arbitrary function of the velocity of propagationm of the moving curve. We assume that m>1 (m is appropriatelynormalized), for which the system of KCL becomes hyperbolic.We interpret the images of the elementary wave solutions inthe (,t)-plane to the (x,y)-plane as elementary shapes of themoving curve (or a nonlinear wavefront when interpreted in aphysical system) and then describe their geometrical properties.Solutions of the Riemann problem with different initial datagive the shapes of the nonlinear wavefront with different combinationsof elementary shapes. Finally, we study all possible interactionsof elementary shapes.     

Symmetries of Surface Singularities

The study of reductive group actions on a normal surface singularityX is facilitated by the fact that the group Aut X of automorphismsof X has a maximal reductive algebraic subgroup G which containsevery reductive algebraic subgroup of Aut X up to conjugation.If X is not weighted homogeneous then this maximal group G isfinite (Scheja, Wiebe). It has been determined for cusp singularitiesby Wall. On the other hand, if X is weighted homogeneous butnot a cyclic quotient singularity then the connected componentG₁ of the unit coincides with the C* defining the weighted homogeneousstructure (Scheja, Wiebe, Wahl). Thus the main interest liesin the finite group G/G₁. Not much is known about G/G₁. Ganterhas given a bound on its order valid for Gorenstein singularitieswhich are not log-canonical. Aumann-Körber has determinedG/G₁ for all quotient singularities. We propose to study G/G₁ through the action of G on the minimalgood resolution of X. If X is weightedhomogeneous but not a cyclic quotient singularity, let E₀ bethe central curve of the exceptional divisor of . We show that the natural homomorphism GAut E₀ haskernel C* and finite image. In particular, this re-proves therest of Scheja, Wiebe and Wahl mentioned above. Moreover, itallows us to view G/G₁ as a subgroup of Aut E₀. For simple ellipticsingularities it equals (Z_bxZ_b)Aut₀ E₀ where –b is theself-intersection number of E₀, Z_bxZ_b is the group of b-torsionpoints of the elliptic curve E₀ acting by translations, andAut₀ E₀ is the group of automorphisms fixing the zero elementof E₀. If E₀ is rational then G/G₁ is the group of automorphismsof E₀ which permute the intersection points with the branchesof the exceptional divisor while preserving the Seifert invariantsof these branches. When there are exactly three branches weconclude that G/G₁ is isomorphic to the group of automorphismsof the weighted resolution graph. This applies to all non-cyclicquotient singularities as well as to triangle singularities.We also investigate whether the maximal reductive automorphismgroup is a direct product GG₁xG/G₁. This is the case, for instance,if the central curve E₀ is rational of even self-intersectionnumber or if X is Gorenstein such that its nowhere-zero 2-form has degree ±1. In the latter case there is a ‘natural’section G/G₁G of GG/G₁ given by the group of automorphisms inG which fix . For a simple elliptic singularity one has GG₁xG/G₁if and only if –E₀ · E₀ = 1.     

The Characterization of the Regularity of the Jacobian Determinant in the Framework of Bessel Potential Spaces on Domains

Let 2 m n. The paper gives necessary and sufficient conditionson the parameters s1, s2, ..., sm, p1, p2, ..., pm such thatthe Jacobian determinant extends to a bounded operator fromHs1p1 x Hs2p2 x ... x Hsmpm into S'. Here all spaces are definedon Rn or on domains Rn. In almost all cases the regularity ofthe Jacobian determinant is calculated exactly.     

On the Accuracy of Quadrature Laplace Transform Inversion and the Solution of State Space Equations

A quadrature formula is shown to be an approximation of thepower-series method of inverting Laplace transforms. This togetherwith the properties of the constants derived from a power-seriesexpansion of the Pad? approximation to exp (s) yield an importantupper limit on t which is quite sharp in determining the breakdownpoint up to and after which the approximation is accurate andinaccurate respectively. The solution of state space equationsusing the quadrature inversion formula is also discussed.     

On a boundary value problem for a p‐Dirac equation

The p‐Laplace equation is a nonlinear generalization of the Laplace equation. This generalization is often used as a model problem for special types of nonlinearities. The p‐Laplace equation can be seen as a bridge between very general nonlinear equations and the linear Laplace equation. The aim of this paper is to solve the p‐Laplace equation for 1 < p < 2 and to find strong solutions. The idea is to apply a hypercomplex integral operator and spatial function theoretic methods to transform the p‐Laplace equation into the p‐Dirac equation. This equation will be solved iteratively by using a fixed‐point theorem. Applying operator‐theoretical methods for the p‐Dirac equation and p‐Laplace equation, the existence and uniqueness of solutions in certain Sobolev spaces will be proved. Copyright © 2016 John Wiley & Sons, Ltd.     

The norm of the Riemann-Liouville operator on Lp[0,1]: a probabilistic approach

We obtain explicit lower and upper bounds for the norm of theRiemann–Liouville operator V^s on L^p[0, 1] which are asymptoticallysharp, thus completing previous results by Eveson. Similar statementsare shown with respect to the norms ||V^s f||_p, whenever f satisfiescertain smoothness properties. It turns out that the correctrate of convergence of ||V^s f||_p as s depends both on theinfimum of the support of f and on the degree of smoothnessof f. We use a probabilistic approach which allows us to giveunified proofs.     

(k, l)-Algebraic Stability of Gauss Methods

It is well known that the s-stage Gauss Runge-Kutta methodsof order 2s are algebraically stable, or equivalently (1, 0)-algebraicallystable. In this paper, we show that there exists some l_s >0 such that the Gauss methods are (k, l) algebraically stablefor l [0, l_s) with k(l)=e^2l+O(l^p+1, where p=2s if s=1 or s=2,and p=2 if s>3.