Current work and open problems in the numerical analysis of Volterra functional equations with vanishing delays

The aims of this paper are (i) to present a survey of recent advances in the analysis of superconvergence of collocation solutions for linear Volterra-type functional integral and integro-differential equations with delay functions θ(t) vanishing at the initial point of the interval of integration (with ia(t) = qt (0 < q < 1, t ⩾ 0) being an important special case), and (ii) to point, by means of a list of open problems, to areas in the numerical analysis of such Volterra functional equations where more research needs to be carried out.      

Oscillation of second order neutral delay differential equations of Emden-Fowler type

We present new oscillation criteria for the second order nonlinear neutral delay differential equation [y(t)-py(t-τ)]'+ q(t)y ^λ (g(t)) sgn y(g(t)) = 0, t ≧ t ₀. Our results solve an open problem posed by James S.W . Wong [24]. The relevance of our results becomes clear due to a carefully selected example. This revised version was published online in June 2006 with corrections to the Cover Date.     

Existence and uniqueness of periodic solutions for a class of nonlinear n ‐th order differential equations with delays

In this paper, we use the coincidence degree theory to establish new results on the existence and uniqueness of T ‐periodic solutions for a class of nonlinear n ‐th order differential equations with delays of the form x⁽ⁿ⁾(t) + f (x^{(n‐ 1)}(t)) + g (t, x (t ‐ τ (t))) = p (t).     

Geometric meshes in collocation methods for Volterra integral equations with proportional delays

In this paper we analyse the local superconvergence propertiesof iterated piecewise polynomial collocation solutions for linearsecond-kind Volterra integral equations with (vanishing) proportionaldelays qt (0 < q < 1). It is shown that on suitable geometricmeshes depending on q, collocation at the Gauss points leadsto almost optimal superconvergence at the mesh points. Thiscontrasts with collocation on uniform meshes where the problemregarding the attainable order of local superconvergence remainsopen.     

Limit-point type solutions of nonlinear differential equations

We are concerned with the nonexistence of L²-solutions of a nonlinear differential equation x″=a(t)x+f(t,x). By applying technique similar to that exploited by Hallam [SIAM J. Appl. Math. 19 (1970) 430-439] for the study of asymptotic behavior of solutions of this equation, we establish nonexistence of solutions from the class L²(t₀,∞) under milder conditions on the function a(t) which, as the examples show, can be even square integrable. Therefore, the equation under consideration can be classified as of limit-point type at infinity in the sense of the definition introduced by Graef and Spikes [Nonlinear Anal. 7 (1983) 851-871]. We compare our results to those reported in the literature and show how they can be extended to third order nonlinear differential equations.     

Some oscillation criteria for second order nonlinear functional ordinary differential equations

The main objective of this article is to study the oscillatory behavior of the solutions of the following nonlinear functional differential equations (a(t)x'(t))' δ1p(t)x'(t) δ2q(t)f(x(g(t))) = 0,for 0 ≤ t0 ≤ t, where δ1 = ±1 and δ2 = ±1. The functions p,q,g : [t0, ∞) → R, f :R → R are continuous, a(t) ＞ 0, p(t) ≥ 0,q(t) ≥ 0 for t ≥ t0, limt→∞ g(t) = ∞, and q is not identically zero on any subinterval of [t0, ∞). Moreover, the functions q(t),g(t), and a(t) are continuously differentiable.     

Singular solutions of Hessian fully nonlinear elliptic equations

We study Hessian fully nonlinear uniformly elliptic equations and show that the second derivatives of viscosity solutions of those equations (in 12 or more dimensions) can blow up in an interior point of the domain. We prove that the optimal interior regularity of such solutions is no more than C1+?, showing the optimality of the known interior regularity result. The same is proven for Isaacs equations. We prove the existence of non-smooth solutions to fully nonlinear Hessian uniformly elliptic equations in 11 dimensions. We study also the possible singularity of solutions of Hessian equations defined in a neighborhood of a point and prove that a homogeneous order 0<α<1 solution of a Hessian uniformly elliptic equation in a punctured ball should be radial.     

Generating orthogonal matrix polynomials satisfying second order differential equations from a trio of triangular matrices

The method developed in [A.J. Durán, F.A. Grünbaum, Orthogonal matrix polynomials satisfying second order differential equations, Int. Math. Res. Not. 10 (2004) 461–484] led us to consider matrix polynomials that are orthogonal with respect to weight matrices W(t) of the form , , and (1−t)^α(1+t)^βT(t)T^*(t), with T satisfying T^′=(2Bt+A)T, T(0)=I, T^′=(A+B/t)T, T(1)=I, and T^′(t)=(−A/(1−t)+B/(1+t))T, T(0)=I, respectively. Here A and B are in general two non-commuting matrices. We are interested in sequences of orthogonal polynomials (P_n)_n which also satisfy a second order differential equation with differential coefficients that are matrix polynomials F₂, F₁ and F₀ (independent of n) of degrees not bigger than 2, 1 and 0 respectively. To proceed further and find situations where these second order differential equations hold, we only dealt with the case when one of the matrices A or B vanishes.The purpose of this paper is to show a method which allows us to deal with the case when A, B and F₀ are simultaneously triangularizable (but without making any commutativity assumption).     

Remarks on a paper of C. Huang

We prove oscillation and nonoscillation theorems for the second order linear differential equation (E) y″+q(t)y=0, where q(t)?0 and locally integrable on These results are extensions of earlier results of Huang [J. Math. Anal. Appl. 210 (1997) 712-723]. Furthermore, we show that the oscillation criterion established for Eq. (E) can be extended to the delayed differential equation y″+q(t)y(σ(t))=0, where σ(t)?t and lim_t→∞σ(t)=∞.     

Superconvergence of the direct discontinuous Galerkin method for a time‐fractional initial‐boundary value problem

A time‐fractional reaction–diffusion initial‐boundary value problem with periodic boundary condition is considered on Q ? Ω × [0, T] , where Ω is the interval [0, l] . Typical solutions of such problem have a weak singularity at the initial time t = 0. The numerical method of the paper uses a direct discontinuous Galerkin (DDG) finite element method in space on a uniform mesh, with piecewise polynomials of degree k ≥ 2 . In the temporal direction we use the L1 approximation of the Caputo derivative on a suitably graded mesh. We prove that at each time level of the mesh, our L1‐DDG solution is superconvergent of order k + 2 in L²(Ω) to a particular projection of the exact solution. Moreover, the L1‐DDG solution achieves superconvergence of order (k + 2) in a discrete L²(Q) norm computed at the Lobatto points, and order (k + 1) superconvergence in a discrete H¹(Q) seminorm at the Gauss points; numerical results show that these estimates are sharp.     

On the convergence and superconvergence for a class of two-dimensional time fractional reaction-subdiffusion equations

This paper analyzes a class of two-dimensional (2-D) time fractional reaction-subdiffusion equations with variable coefficients. The high-order L2-1_σ time-stepping scheme on graded meshes is presented to deal with the weak singularity at the initial time t = 0, and the bilinear finite element method (FEM) on anisotropic meshes is used for spatial discretization. Using the modified discrete fractional Grönwall inequality, and combining the interpolation operator and the projection operator, the L²-norm error estimation and H¹-norm superclose results are rigorously proved. The superconvergence result in the H¹-norm is derived by applying the interpolation postprocessing technique. Finally, numerical examples are presented to verify the validation of our theoretical analysis.     

On the oscillation of Hill's equations under periodic forcing

We use the Floquet theory of the Hill's equation to prove the conjecture that all solutions of the second order forced linear differential equation y^″+c(sint)y=cost, are oscillatory on [0,∞) for all c≠0.     

The minimax approach for a class of variable order fractional differential equation

This paper introduces an approximate solution for Liouville‐Caputo variable order fractional differential equations with order 0 < α(t) ≤ 1 . The solution is adapted using a family of fractional‐order Chebyshev functions with unknown coefficients. These coefficients have been obtained by using an optimization approach based on minimax technique and the least pth optimization function. Several linear and nonlinear fractional‐order differential equations are discussed using the proposed technique for fixed and variable order fractional‐order derivatives. Moreover, the response of RC charging circuit with variable order fractional capacitor is studied for different cases. Several comparisons with related published techniques have been added to illustrate the accuracy of the proposed approach.     

On solvability and regularity of a parametrized version of optimality conditions

We investigate a linear homotopyF(·,t) connecting an appropriate smooth equationG=0 with Kojima's (nonsmooth) systemK=0 describing critical points (primal —dual) of a nonlinear optimization problem (NLP) in finite dimension.Fort=0, our system may be seen e.g. as a starting system for an embedding procedure to determine a critical point to NLP. Fort1, it may be regarded as a regularization ofK.Conditions for regularity (necessary and sufficient) and solvability (sufficient) are studied. Though, formally, they can be given in a unified way, we show that their meaning differs fort < 1 andt=1. Particularily, no MFCQ-like condition must be imposed in order to ensure regularity fort < 1.     

Spreading Properties and Complex Dynamics for Monostable Reaction–Diffusion Equations

This paper is concerned with the study of the large-time behavior of the solutions u of a class of one-dimensional reaction–diffusion equations with monostable reaction terms f, including in particular the classical Fisher-KPP nonlinearities. The nonnegative initial data u ₀(x) are chiefly assumed to be exponentially bounded as x tends to + ∞ and separated away from the unstable steady state 0 as x tends to ? ∞. On the one hand, we give some conditions on u ₀ which guarantee that, for some λ > 0, the quantity c _λ = λ +f′(0)/λ is the asymptotic spreading speed, in the sense that lim  _t→+∞ u(t, ct) = 1 (the stable steady state) if c < c _λ and lim  _t→+∞ u(t, ct) = 0 if c > c _λ. These conditions are fulfilled in particular when u ₀(x) e ^λx is asymptotically periodic as x → + ∞. On the other hand, we also construct examples where the spreading speed is not uniquely determined. Namely, we show the existence of classes of initial conditions u ₀ for which the ω-limit set of u(t, ct + x) as t tends to + ∞ is equal to the whole interval [0, 1] for all x ∈ ? and for all speeds c belonging to a given interval (γ₁, γ₂) with large enough γ₁ < γ₂ ≤ + ∞.     

ON THE DIFFUSION PHENOMENON OF QUASILINEAR HYPERBOLICWAVES

Introduction1.1.ConsiderthefollowingquasilinearhyperbolicCauchyproblemwithlineardamping{:;;!OTt=-:i<:,;>>L06,(11)wherexER",t20,anda(.)isasmoothfunctionsatisfyinga(y)~1 O(lyl")aslyl-0,orEN.(1.2)Thepurposeofthispaperistoshowthat,atleastwhenn53,theasymptoticprofileofthesolutionu(x,t)of(l.1)isgivenbythesolutionv(x,t)ofthecorrespondingparabolicproblem{:;.t>ivj:相似文献   

Boundary penalty finite element methods for blending surfaces, . Superconvergence and stability and examples

This paper is Part III of the study on blending surfaces by partial differential equations (PDEs). The blending surfaces in three dimensions (3D) are taken into account by three parametric functions, x(r,t),y(r,t) and z(r,t). The boundary penalty techniques are well suited to the complicated tangent (i.e., normal derivative) boundary conditions in engineering blending. By following the previous papers, Parts I and II in Li (J. Comput. Math. 16 (1998) 457–480; J. Comput. Appl. Math. 110 (1999) 155–176) the corresponding theoretical analysis is made to discover that when the penalty power σ=2, σ=3 (or 3.5) and 0<σ1.5 in the boundary penalty finite element methods (BP-FEMs), optimal convergence rates, superconvergence and optimal numerical stability can be achieved, respectively. Several interesting samples of 3D blending surfaces are provided, to display the remarkable advantages of the proposed approaches in this paper: unique solutions of blending surfaces, optimal blending surfaces in minimum energy, ease in handling the complicated boundary constraint conditions, and less CPU time and computer storage needed. This paper and Li (J. Comput. Math. 16 (1998) 457–480; J. Comput. Appl. Math.) provide a foundation of blending surfaces by PDE solutions, a new trend of computer geometric design.     

Superconvergence of finite volume element method for a nonlinear elliptic problem

In this article, we consider the finite volume element method for the second‐order nonlinear elliptic problem and obtain the H¹ and W^1,∞ superconvergence estimates between the solution of the finite volume element method and that of the finite element method, which reveal that the finite volume element method is in close relationship with the finite element method. With these superconvergence estimates, we establish the L^p and W^1,p (2 < p ≤ ∞) error estimates for the finite volume element method for the second‐order nonlinear elliptic problem. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

On Lagrange interpolation for |<Emphasis Type="Italic">X</Emphasis>|<Superscript><Emphasis Type="Italic">α</Emphasis></Superscript> (0 < <Emphasis Type="Italic">α</Emphasis> < 1)

We study the optimal order of approximation for |x| ^α (0 < α < 1) by Lagrange interpolation polynomials based on Chebyshev nodes of the first kind. It is proved that the Jackson order of approximation is attained. Supported by the NSFC, 10601065.     

A sufficient condition for the existence of a “dead zone” for solutions of degenerate semilinear parabolic equations and inequalities

We consider the solutions of degenerate parabolic equations and inequalities of the formLu-u _t = |u| ^q sgnu and sgnu(Lu−u _t )−|u| ^q ≥0, 0<q<1, with the elliptic operatorL in divergent or nondivergent form. We establish a dependence of the maximum modulus of the solution on the domain and on the equation (inequality) such that this dependence guarantees the existence of a “dead zone” of the solution. In this case, the character of degeneracy is unessential. Translated fromMatematicheskie Zametki, Vol. 60, No. 6, pp. 824–831, December, 1996.