Nonlinear Steepest Descent and Numerical Solution of Riemann‐Hilbert Problems

The effective and efficient numerical solution of Riemann‐Hilbert problems has been demonstrated in recent work. With the aid of ideas from the method of nonlinear steepest descent for Riemann‐Hilbert problems, the resulting numerical methods have been shown, in practice, to retain accuracy as values of certain parameters become arbitrarily large. Remarkably, this numerical approach does not require knowledge of local parametrices; rather, the deformed contour is scaled near stationary points at a specific rate. The primary aim of this paper is to prove that this observed asymptotic accuracy is indeed achieved. To do so, we first construct a general theoretical framework for the numerical solution of Riemann‐Hilbert problems. Second, we demonstrate the precise link between nonlinear steepest descent and the success of numerics in asymptotic regimes. In particular, we prove sufficient conditions for numerical methods to retain accuracy. Finally, we compute solutions to the homogeneous Painlevé II equation and the modified Korteweg–de Vries equation to explicitly demonstrate the practical validity of the theory. © 2014 Wiley Periodicals, Inc.     

Riemann–Hilbert Problems,Matrix Orthogonal Polynomials and Discrete Matrix Equations with Singularity Confinement

In this paper, matrix orthogonal polynomials in the real line are described in terms of a Riemann–Hilbert problem. This approach provides an easy derivation of discrete equations for the corresponding matrix recursion coefficients. The discrete equation is explicitly derived in the matrix Freud case, associated with matrix quartic potentials. It is shown that, when the initial condition and the measure are simultaneously triangularizable, this matrix discrete equation possesses the singularity confinement property, independently if the solution under consideration is given by the recursion coefficients to quartic Freud matrix orthogonal polynomials or not.     

Riemann‐Hilbert problem and matrix discrete Painlevé II systems

Matrix Szeg? biorthogonal polynomials for quasi‐definite matrices of Hölder continuous weights are studied. A Riemann‐Hilbert problem is uniquely solved in terms of the matrix Szeg? polynomials and its Cauchy transforms. The Riemann‐Hilbert problem is given as an appropriate framework for the discussion of the Szeg? matrix and the associated Szeg? recursion relations for the matrix orthogonal polynomials and its Cauchy transforms. Pearson‐type differential systems characterizing the matrix of weights are studied. These are linear systems of ordinary differential equations that are required to have trivial monodromy. Linear ordinary differential equations for the matrix Szeg? polynomials and its Cauchy transforms are derived. It is shown how these Pearson systems lead to nonlinear difference equations for the Verblunsky matrices and two examples, of Fuchsian and non‐Fuchsian type, are considered. For both cases, a new matrix version of the discrete Painlevé II equation for the Verblunsky matrices is found. Reductions of these matrix discrete Painlevé II systems presenting locality are discussed.     

A remark on the Riemann‐Hilbert problem with a vanishing coefficient

This paper concerns the existence of nontrivial solutions of the Riemann‐Hilbert problem with a continuous coefficient whose values belong to two rays in the complex plane. Our results extend those recently obtained by E. Shargorodsky and J. F. Toland [6]. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

N solutions for a derivative nonlinear Schrödinger‐type equation via Riemann‐Hilbert approach

The inverse scattering transform for the derivative nonlinear Schrödinger‐type equation is studied via the Riemann‐Hilbert approach. In the direct scattering process, the spectral analysis of the Lax pair is performed, from which a Riemann‐Hilbert problem is established for the derivative nonlinear Schrödinger‐type equation. In the inverse scattering process, N‐soliton solutions of the derivative nonlinear Schrödinger‐type equation are obtained by solving Riemann‐Hilbert problems corresponding to the reflectionless cases. Moreover, the dynamics of the exact solutions are discussed.     

Double scaling limit in the random matrix model: The Riemann‐Hilbert approach

We prove the existence of the double scaling limit in the unitary matrix model with quartic interaction, and we show that the correlation functions in the double scaling limit are expressed in terms of the integrable kernel determined by the ψ function for the Hastings‐McLeod solution to the Painlevé II equation. The proof is based on the Riemann‐Hilbert approach, and the central point of the proof is an analysis of analytic semiclassical asymptotics for the ψ function at the critical point in the presence of four coalescing turning points. © 2003 Wiley Periodicals, Inc.     

Riemann‐Hilbert problems and soliton solutions of a multicomponent mKdV system and its reduction

An arbitrary order matrix spectral problem is introduced and its associated multicomponent AKNS integrable hierarchy is constructed. Based on this matrix spectral problem, a kind of Riemann‐Hilbert problems is formulated for a multicomponent mKdV system in the resulting AKNS integrable hierarchy. Through special corresponding Riemann‐Hilbert problems with an identity jump matrix, soliton solutions to the presented multicomponent mKdV system are explicitly worked out. A specific reduction of the multicomponent mKdV system is made, together with its reduced Lax pair and soliton solutions.     

On the Riemann‐Hilbert boundary value problem for generalized analytic functions in the framework of variable exponent spaces

Let Γ be a simple closed curve that bounds the finite domain D , z =z (ζ )=z (r e ^{i ?} ) be the conformal mapping of the circle {ζ :|ζ |<1} onto the domain D . Furthermore, let the functions A (z ), B (z ) be given on D and U ^{s ,2}(A ;B ;D ) be the set of regular solutions of the equation . We call the Smirnov class E ^{p (t )}(A ;B ;D ) the set of those generalized functions W in D for which where p (t ) is a positive measurable function on Γ. We consider the Riemann‐Hilbert problem: Define a function W (z ) from the class E ^{p (t )}(A ;B ;D ) for which the equality, is fulfilled almost everywhere on Γ. It is assumed that Γ is a piecewise‐smooth curve without external peaks; , p is Log Hölder continuous and the function belongs to the class A (p (t );Γ), which is the generalization of the well‐known Simonenko class A (p ;Γ), where p is constant. The solvability conditions are established, and solutions are constructed.     

Singular Value Decomposition of a Finite Hilbert Transform Defined on Several Intervals and the Interior Problem of Tomography: The Riemann‐Hilbert Problem Approach

We study the asymptotics of singular values and singular functions of a finite Hilbert transform (FHT), which is defined on several intervals. Transforms of this kind arise in the study of the interior problem of tomography. We suggest a novel approach based on the technique of the matrix Riemann‐Hilbert problem (RHP) and the steepest‐descent method of Deift‐Zhou. We obtain a family of matrix RHPs depending on the spectral parameter λ and show that the singular values of the FHT coincide with the values of λ for which the RHP is not solvable. Expressing the leading‐order solution as λ → 0 of the RHP in terms of the Riemann Theta functions, we prove that the asymptotics of the singular values can be obtained by studying the intersections of the locus of zeroes of a certain Theta function with a straight line. This line can be calculated explicitly, and it depends on the geometry of the intervals that define the FHT. The leading‐order asymptotics of the singular functions and singular values are explicitly expressed in terms of the Riemann Theta functions and of the period matrix of the corresponding normalized differentials, respectively. We also obtain the error estimates for our asymptotic results. © 2016 Wiley Periodicals, Inc.     

Cubic and Quartic Transformations of the Sixth Painlevé Equation in Terms of Riemann–Hilbert Correspondence

The starting point of this paper is a classification of quadratic polynomial transformations of the monodromy manifold for the 2 × 2 isomonodromic Fuchsian systems associated to the Painlevé VI equation. Up to birational automorphisms of the monodromy manifold, we find three transformations. Two of them are identified as the action of known quadratic or quartic transformations of the Painlevé VI equation. The third transformation of the monodromy manifold gives a new transformation of degree 3 of Picard’s solutions of Painlevé VI.     

Extremal solutions for p‐Laplacian boundary value problems with the right‐handed Riemann‐Liouville fractional derivative

The paper is concerned with the solvability for several nonlinear boundary value problems of fractional p‐Laplacian differential equation involving the right‐handed Riemann‐Liouville derivative. By applying monotone iterative technique, lower and upper solutions method and the Banach fixed point theorem, sufficient conditions for existence and uniqueness of extremal solutions are obtained and they extend existing results. At last, two examples are provided to illustrate the results.     

Scalar and matrix Riemann–Hilbert approach to the strong asymptotics of Padé approximants and complex orthogonal polynomials with varying weight

We describe methods for the derivation of strong asymptotics for the denominator polynomials and the remainder of Padé approximants for a Markov function with a complex and varying weight. Two approaches, both based on a Riemann–Hilbert problem, are presented. The first method uses a scalar Riemann–Hilbert boundary value problem on a two-sheeted Riemann surface, the second approach uses a matrix Riemann–Hilbert problem. The result for a varying weight is not with the most general conditions possible, but the loss of generality is compensated by an easier and transparent proof.     

Two‐Timing Hypothesis,Distinguished Limits,Drifts, and Pseudo‐Diffusion for Oscillating Flows

The aim of this paper is: using the two‐timing method to study and classify the multiplicity of distinguished limits and asymptotic solutions for the advection equation with a general oscillating velocity field. Our results are:

• (i) the dimensionless advection equation that contains two independent small parameters, which represent the ratio of two characteristic time‐scales and the spatial amplitude of oscillations; the related scaling of the variables and parameters uses the Strouhal number;
• (ii) an infinite sequence of distinguished limits has been identified; this sequence corresponds to the successive degenerations of a drift velocity;
• (iii) we have derived the averaged equations and the oscillatory equations for the first four distinguished limits; derivations are performed up to the fourth orders in small parameters;
• (v) we have shown, that each distinguished limit generates an infinite number of parametric solutions; these solutions differ from each other by the slow time‐scale and the amplitude of the prescribed velocity;
• (vi) we have discovered the inevitable presence of pseudo‐diffusion terms in the averaged equations, pseudo‐diffusion appears as a Lie derivative of the averaged tensor of quadratic displacements; we have analyzed the matrix of pseudo‐diffusion coefficients and have established its degenerated form and hyperbolic character; however, for one‐dimensional cases, the pseudo‐diffusion can appear as ordinary diffusion;
• (vii) the averaged equations for four different types of oscillating velocity fields have been considered as the examples of different drifts and pseudo‐diffusion;
• (viii) our main methodological result is the introduction of a logical order into the area and classification of an infinite number of asymptotic solutions; we hope that it can help in the study of the similar problems for more complex systems;
• (ix) our study can be used as a test for the validity of the two‐timing hypothesis, because in our calculations we do not employ any additional assumptions.

     

Some orthogonal polynomials on the finite interval and Gaussian quadrature rules for fractional Riemann‐Liouville integrals

Inspired with papers by Bokhari, Qadir, and Al‐Attas (2010) and by Rapai?, ?ekara, and Govedarica (2014), in this paper we investigate a few types of orthogonal polynomials on finite intervals and derive the corresponding quadrature formulas of Gaussian type for efficient numerical computation of the left and right fractional Riemann‐Liouville integrals. Several numerical examples are included to demonstrate the numerical efficiency of the proposed procedure.     

The Bernoulli Integral for a Certain Class of Non‐Stationary Viscous Vortical Flows of Incompressible Fluid

It has been shown in our paper [1] that there is a wide class of 3D motions of incompressible viscous fluid which can be described by one scalar function dabbed the quasi‐potential. This class of fluid flows is characterized by three‐component velocity field having two‐component vorticity field; both these fields can depend of all three spatial variables and time, in general. Governing equations for the quasi‐potential have been derived and simple illustrative example of 3D flow has been presented. Here, we derive the Bernoulli integral for that class of flows and compare it against the known Bernoulli integrals for the potential flows or 2D stationary vortical flows of inviscid fluid. We show that the Bernoulli integral for this class of fluid motion possesses unusual features: it is valid for the vortical nonstationary motions of a viscous incompressible fluid. We present a new very nontrivial analytical example of 3D flow with two‐component vorticity which hardly can be obtained by any of known methods. In the last section, we suggest a generalization of the developed concept which allows one to describe a certain class of 3D flows with the 3D vorticity.     

Weak solutions for fractional differential equations in nonreflexive Banach spaces via Riemann‐Pettis integrals

The aim of this paper is to develop fractional calculus for vector‐valued functions using the weak Riemann integral. Also, we establish the existence of weak solutions for a class of fractional differential equations with fractional weak derivatives.     

Hilbert空间中广义平衡问题和不动点问题的粘滞逼近法

在Hilbert空间,我们用粘滞逼近法建立了一迭代序列来逼近两个集合的公共点,这两个集合分别是广义平衡问题的解集和渐进非扩张映射的不动点集.我们表明这一迭代序列强收敛到这两个集合的公共点,而且这一公共点还是一变分不等式的解.用这一结果,还研究了三个强收敛问题和优化问题.     

一类二维流场模型的周期解与浑沌解的共存性*

本文研究Kelvin-Stuart猫眼流在周期扰动下的动力学行为,运用Melnikov方法确定出振动型周期轨产生偶数阶次谐分枝、旋转型周期轨产生任意阶次谐分枝的条件,并进一步发现周期解与浑沌解共存的复杂现象.