Modified Schwarz–Christoffel mappings using approximate curve factors

The Schwarz–Christoffel mapping from the upper half-plane to a polygonal region in the complex plane is an integral of a product with several factors, where each factor corresponds to a certain vertex in the polygon. Different modifications of the Schwarz–Christoffel mapping in which factors are replaced with the so-called curve factors to achieve polygons with rounded corners are known since long times. Among other requisites, the arguments of a curve factor and its correspondent scl factor must be equal outside some closed interval on the real axis.In this paper, the term approximate curve factor is defined such that many of the already known curve factors are included as special cases. Additionally, by alleviating the requisite on the argument from exact to asymptotic equality, new types of curve factors are introduced. While traditional curve factors have a C¹ regularity, C^∞ regular approximate curve factors can be constructed, resulting in smooth boundary curves when used in conformal mappings.Applications include modelling of wave scattering in waveguides. When using approximate curve factors in modified Schwarz–Christoffel mappings, numerical conformal mappings can be constructed that preserve two important properties in the waveguides. First, the direction of the boundary curve can be well controlled, especially towards infinity, where the application requires two straight parallel walls. Second, a smooth (C^∞) boundary curve can be achieved.     

Interpolation approximations based on Gauss–Lobatto–Legendre–Birkhoff quadrature

We derive in this paper the asymptotic estimates of the nodes and weights of the Gauss–Lobatto–Legendre–Birkhoff (GLLB) quadrature formula, and obtain optimal error estimates for the associated GLLB interpolation in Jacobi weighted Sobolev spaces. We also present a user-oriented implementation of the pseudospectral methods based on the GLLB quadrature nodes for Neumann problems. This approach allows an exact imposition of Neumann boundary conditions, and is as efficient as the pseudospectral methods based on Gauss–Lobatto quadrature for PDEs with Dirichlet boundary conditions.     

On numerical improvement of Gauss–Legendre quadrature rules

Among all integration rules with n points, it is well-known that n-point Gauss–Legendre quadrature rule∫₋₁¹f(x) dx∑_i=1ⁿw_if(x_i)has the highest possible precision degree and is analytically exact for polynomials of degree at most 2n−1, where nodes x_i are zeros of Legendre polynomial P_n(x), and w_i's are corresponding weights.In this paper we are going to estimate numerical values of nodes x_i and weights w_i so that the absolute error of introduced quadrature rule is less than a preassigned tolerance ε₀, say ε₀=10⁻⁸, for monomial functionsf(x)=x^j, j=0,1,…,2n+1.(Two monomials more than precision degree of Gauss–Legendre quadrature rules.) We also consider some conditions under which the new rules act, numerically, more accurate than the corresponding Gauss–Legendre rules. Some examples are given to show the numerical superiority of presented rules.     

Riemann–Hilbert analysis for Jacobi polynomials orthogonal on a single contour

Classical Jacobi polynomials , with α,β>-1, have a number of well-known properties, in particular the location of their zeros in the open interval (-1,1). This property is no longer valid for other values of the parameters; in general, zeros are complex. In this paper we study the strong asymptotics of Jacobi polynomials where the real parameters α_n,β_n depend on n in such a way that

with . We restrict our attention to the case where the limits A,B are not both positive and take values outside of the triangle bounded by the straight lines A=0, B=0 and A+B+2=0. As a corollary, we show that in the limit the zeros distribute along certain curves that constitute trajectories of a quadratic differential.The non-hermitian orthogonality relations for Jacobi polynomials with varying parameters lie in the core of our approach; in the cases we consider, these relations hold on a single contour of the complex plane. The asymptotic analysis is performed using the Deift–Zhou steepest descent method based on the Riemann–Hilbert reformulation of Jacobi polynomials.     

Uniform Convergence of Fourier–Jacobi Series

Necessary and sufficient conditions which imply the uniform convergence of the Fourier–Jacobi series of a continuous function are obtained under an assumption that the Fourier–Jacobi series is convergent at the end points of the segment of orthogonality [−1,1]. The conditions are in terms of the modulus of continuity, Λ-variation, and the modulus of variation of a function.     

Composite generalized Laguerre–Legendre pseudospectral method for Fokker–Planck equation in an infinite channel

In this paper, we propose a composite generalized Laguerre–Legendre pseudospectral method for the Fokker–Planck equation in an infinite channel, which behaves like a parabolic equation in one direction, and behaves like a hyperbolic equation in other direction. We establish some approximation results on the composite generalized Laguerre–Legendre–Gauss–Radau interpolation, with which the convergence of proposed composite scheme follows. An efficient implementation is provided. Numerical results show the spectral accuracy in space of this approach and coincide well with theoretical analysis. The approximation results and techniques developed in this paper are also very appropriate for many other problems on multiple-dimensional unbounded domains, which are not of standard types.     

Higher order representations of the Robbins–Monro process

For quasi-linear regression functions, the Robbins–Monro process X_n is decomposed in a sum of a linear form and a quadratic form both defined in the observation errors. Under regularity conditions, the remainder term is of order O(n^−3/2) with respect to the L_p-norm. If a cubic form is added, the remainder term can be improved up to an order of O(n⁻²). As a corollary the expectation of X_n is expanded up to an error of order O(n⁻²). This is used to correct the bias of X_n up to an error of order O(n^−3/2 log n).     

On the spectral norms of Cauchy–Toeplitz and Cauchy–Hankel matrices

In this study, we have found upper and lower bounds for the spectral norm of Cauchy–Toeplitz and Cauchy–Hankel matrices in the forms T_n=[1/(a+(i−j)b)]ⁿ_i,j=1, H_n=[1/(a+(i+j)b)]ⁿ_i,j=1.     

An a Posteriori Error Estimate for a Semi-Lagrangian Scheme for Hamilton–Jacobi Equations

We present an a posteriori estimate for a first order semi-Lagrangian method for Hamilton–Jacobi equations. The result requires piecewise C ^1,1 regularity of the viscosity solution and is stated for the Bellman equation related to the infinite horizon problem, although it can be applied to more general Hamilton–Jacobi equations with convex Hamiltonians. This estimate suggests different numerical indicators that can be used to construct an adaptive algorithm for the approximation of the viscosity solution.     

Quadratic Hermite–Padé polynomials associated with the exponential function

The asymptotic behavior of quadratic Hermite–Padé polynomials associated with the exponential function is studied for n→∞. These polynomials are defined by the relation

(*)

p_n(z)+q_n(z)e^z+r_n(z)e^2z=O(z³ⁿ⁺²) as z→0,

where O(·) denotes Landau's symbol. In the investigation analytic expressions are proved for the asymptotics of the polynomials, for the asymptotics of the remainder term in (*), and also for the arcs on which the zeros of the polynomials and of the remainder term cluster if the independent variable z is rescaled in an appropriate way. The asymptotic expressions are defined with the help of an algebraic function of third degree and its associated Riemann surface. Among other possible applications, the results form the basis for the investigation of the convergence of quadratic Hermite–Padé approximants, which will be done in a follow-up paper.     

Sphericity test in a GMANOVA–MANOVA model with normal error

For the GMANOVA–MANOVA model with normal error: , , we study in this paper the sphericity hypothesis test problem with respect to covariance matrix: Σ=λI_q (λ is unknown). It is shown that, as a function of the likelihood ratio statistic Λ, the null distribution of Λ^2/n can be expressed by Meijer’s function, and the asymptotic null distribution of −2logΛ is (as n→∞). In addition, the Bartlett type correction −2ρlogΛ for logΛ is indicated to be asymptotically distributed as with order n⁻² for an appropriate Bartlett adjustment factor −2ρ under null hypothesis.     

The Bishop–Phelps–Bollobás theorem fails for bilinear forms on

In this paper we show that the Bishop–Phelps–Bollobás theorem fails for bilinear forms on l₁×l₁, while it holds for linear operators from l₁ to l_∞.     

Global smoothness and uniform convergence of smooth Gauss–Weierstrass singular operators

In this article we continue with the study of smooth Gauss–Weierstrass singular integral operators over the real line regarding their simultaneous global smoothness preservation property with respect to the L_p norm, 1≤p≤∞, by involving higher order moduli of smoothness. Also we study their simultaneous approximation to the unit operator with rates involving the modulus of continuity with respect to the uniform norm. The produced Jackson type inequalities are almost sharp containing elegant constants, and they reflect the high order of differentiability of the engaged function.     

Group classification of the generalized Emden–Fowler-type equation

We perform the Lie group classification of the Emden–Fowler-type equation xu^″+nu^′+x^νF(u)=0, which arises in several applications. These include the theory of stellar structure, the thermal behaviour of a spherical cloud of gas, isothermal gas spheres and the theory of thermionic currents. Seven cases arise for the possible extension of the principal Lie algebra, which in this case is trivial. Three new cases occur for which we have non-trivial Lie point symmetry algebra. We compare these cases with the Noether symmetry cases. Moreover, we also make comparisons with the partial Noether operators. Finally for three cases we reduce the Emden–Fowler-type equation to quadratures.     

Proof of the Loebl–Komlós–Sós conjecture for large, dense graphs

The Loebl–Komlós–Sós conjecture states that for any integers k and n, if a graph G on n vertices has at least n/2 vertices of degree at least k, then G contains as subgraphs all trees on k+1 vertices. We prove this conjecture in the case when k is linear in n, and n is sufficiently large.     

Asymptotic behavior for the Navier–Stokes equations in 2D exterior domains

We show that the L^p spatial–temporal decay rates of solutions of incompressible flow in an 2D exterior domain. When a domain has a boundary, pressure term makes an obstacle since we do not have enough information on the pressure term near the boundary. To overcome the difficulty, we adopt the ideas in He, Xin [C. He, Z. Xin, Weighted estimates for nonstationary Navier–Stokes equations in exterior domain, Methods Appl. Anal. 7 (3) (2000) 443–458], and our previous results [H.-O. Bae, B.J. Jin, Asymptotic behavior of Stokes solutions in 2D exterior domains, J. Math. Fluid Mech., in press; H.-O. Bae, B.J. Jin, Temporal and spatial decay rates of Navier–Stokes solutions in exterior domains, submitted for publication]. For the spatial decay rate estimate, we first extend temporal decay rate result of the Navier–Stokes solutions for general L^p space when the initial velocity is in , 1<rq<∞ (1<r<q=∞).     

On a problem of Duke–Erdős–Rödl on cycle-connected subgraphs

In this short note, we prove that for β<1/5 every graph G with n vertices and n^2−β edges contains a subgraph G^′ with at least cn^2−2β edges such that every pair of edges in G^′ lie together on a cycle of length at most 8. Moreover edges in G^′ which share a vertex lie together on a cycle of length at most 6. This result is best possible up to the constant factor and settles a conjecture of Duke, Erdős, and Rödl.     

A rotation method which gives linear estimates for powers of the Ahlfors–Beurling operator

In [O. Dragičević, A. Volberg, Sharp estimate of the Ahlfors–Beurling operator via averaging martingale transforms, Michigan Math. J. 51 (2) (2003) 415–435] the Ahlfors–Beurling operator T was represented as an average of two-dimensional martingale transforms. The same result can be proven for powers Tⁿ. Motivated by [T. Iwaniec, G. Martin, Riesz transforms and related singular integrals, J. Reine Angew. Math. 473 (1996) 25–57], we deduce from here that Tⁿ_p are bounded from above by Cnp^*, . We further improve this estimate to obtain the optimal behaviour of the L^p norms in question.     

A robust semi-explicit difference scheme for the Kuramoto–Tsuzuki equation

In this paper, we propose a robust semi-explicit difference scheme for solving the Kuramoto–Tsuzuki equation with homogeneous boundary conditions. Because the prior estimate in L_∞-norm of the numerical solutions is very hard to obtain directly, the proofs of convergence and stability are difficult for the difference scheme. In this paper, we first prove the second-order convergence in L₂-norm of the difference scheme by an induction argument, then obtain the estimate in L_∞-norm of the numerical solutions. Furthermore, based on the estimate in L_∞-norm, we prove that the scheme is also convergent with second order in L_∞-norm. Numerical examples verify the correction of the theoretical analysis.     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.