Error analysis of explicit TSERKN methods for highly oscillatory systems

In this paper, we are concerned with the error analysis for the two-step extended Runge-Kutta-Nyström-type (TSERKN) methods [Comput. Phys. Comm. 182 (2011) 2486–2507] for multi-frequency and multidimensional oscillatory systems y″(t) + My(t) = f(t, y(t)), where high-frequency oscillations in the solutions are generated by the linear part My(t). TSERKN methods extend the two-step hybrid methods [IMA J. Numer. Anal. 23 (2003) 197–220] by reforming both the internal stages and the updates so that they are adapted to the oscillatory properties of the exact solutions. However, the global error analysis for the TSERKN methods has not been investigated. In this paper we construct a new three-stage explicit TSERKN method of order four and present the global error bound for the new method, which is proved to be independent of ∥M∥ under suitable assumptions. This property of our new method is very important for solving highly oscillatory systems (1), where ∥M∥ may be arbitrarily large. We also analyze the stability and phase properties for the new method. Numerical experiments are included and the numerical results show that the new method is very competitive and promising compared with the well-known high quality methods proposed in the scientific literature.     

Envelope viscosity solutions of first- and second-order PDEs with u-dependence

General envelope methods are introduced which may be used to embed equations with u-dependence into equations without solution dependence. Furthermore, these methods present a rigorous way to consider so-called nodal solutions. That is, if w(t,x,z) is the viscosity solution of some pde, the nodal solution of an associated pde is a function u(t,x) so that w(t,x,u(t,x)) = 0. Examples are given to first- and second-order pdes arising in optimal control, differential games, minimal time problems, scalar conservation laws, geometric-type equations, and forward backward stochastic control.     

Analytic-numerical solutions with a priori error bounds for a class of strongly coupled mixed partial differential systems

This paper deals with the construction of analytic-numerical solutions with a priori error bounds for systems of the type u_t = Au_xx, u(0,t) + u_x(0,t) = 0, Bu(1,t) + Cu_x(1,t) = 0, 0 < x < 1, t > 0, u(x,0) = f(x). Here A, B, C are matrices for which no diagonalizable hypothesis is assumed. First an exact series solution is obtained after solving appropriate vector Sturm-Liouville-type problems. Given an admissible error ε and a bounded subdomain D, after appropriate truncation an approximate solution constructed in terms of data and approximate eigenvalues is given so that the error is less than the prefixed accuracy ε, uniformly in D.     

Numerical solutions with a priori error bounds for coupled self-adjoint time dependent partial differential systems

This paper is concerned with the construction of accurate continuous numerical solutions for partial self-adjoint differential systems of the type (P(t) u_t)_t = Q(t)u_xx, u(0, t) = u(d, t) = 0, u(x, 0) = f(x), u_t(x, 0) = g(x), 0 ≤ x ≤ d, t >- 0, where P(t), Q(t) are positive definite oR^r×r-valued functions such that P′(t) and Q′(t) are simultaneously semidefinite (positive or negative) for all t ≥ 0. First, an exact theoretical series solution of the problem is obtained using a separation of variables technique. After appropriate truncation strategy and the numerical solution of certain matrix differential initial value problems the following question is addressed. Given T > 0 and an admissible error ϵ > 0 how to construct a continuous numerical solution whose error with respect to the exact series solution is smaller than ϵ, uniformly in D(T) = {(x, t); 0 ≤ x ≤ d, 0 ≤ t ≤ T}. Uniqueness of solutions is also studied.     

Multiple solutions of some boundary value problems with parameters

In this paper, we study the existence and multiplicity of nontrivial solutions for the following second-order Dirichlet nonlinear boundary value problem with odd order derivative: −u^″(t)+au^′(t)+bu(t)=f(t,u(t)) for all t∈[0,1] with u(0)=u(1)=0, where a,b∈R¹, f∈C¹([0,1]×R¹,R¹). By using the Morse theory, we impose certain conditions on f which are able to guarantee that the problem has at least one nontrivial solution, two nontrivial solutions and infinitely many solutions, separately.     

Positive solutions for systems of nonlinear eigenvalue problems for functional differential equations

Values of?λ?are determined for which there exist positive solutions of the system of functional differential equations, u″?+?λa(t)f(v _t )?=?0,v″?+?λb(t)g(u _t )?=?0, for 0?t?u(s)?=?v(s)?=?φ(s), ?r?≤?s?≤?0, and the boundary conditions u(0)?=?v(0)?=?φ(0)?=?u(1)?=?v(1)?=?0. A Guo–Krasnosel'skii fixed point theorem is applied.     

A note on the optimal temporal decay estimates of solutions to the Cahn-Hilliard equation

This paper is concerned with the optimal temporal decay estimates on the solutions of the Cauchy problem of the Cahn-Hilliard equation. It is shown in Liu, Wang and Zhao (2007) [11] that such a Cauchy problem admits a unique global smooth solution u(t,x) provided that the smooth nonlinear function φ(u) satisfies a local growth condition. Furthermore if φ(u) satisfies a somewhat stronger local growth condition, the optimal temporal decay estimates on u(t,x) are also obtained in Liu, Wang and Zhao (2007) [11]. Thus a natural question is how to deduce the optimal temporal decay estimates on u(t,x) only under the local growth condition which is sufficient to guarantee the global solvability of the corresponding Cauchy problem and the main purpose of this paper is devoted to this problem. Our analysis is motivated by the technique developed recently in Ukai, Yang and Zhao (2006) [15] with a slight modification.     

Eigenvalue problems for second-order nonlinear dynamic equations on time scales

This paper is concerned with the existence and nonexistence of positive solutions of the second-order nonlinear dynamic equation uΔΔ(t)+λa(t)f(u(σ(t)))=0, t∈[0,1], satisfying either the conjugate boundary conditions u(0)=u(σ(1))=0 or the right focal boundary conditions u(0)=uΔ(σ(1))=0, where a and f are positive. We show that there exists a λ^∗>0 such that the above boundary value problem has at least two, one and no positive solutions for 0<λ<λ^∗, λ=λ^∗ and λ>λ^∗, respectively. Furthermore, by using the semiorder method on cones of the Banach space, we establish an existence and uniqueness criterion for positive solution of the problem. In particular, such a positive solution uλ(t) of the problem depends continuously on the parameter λ, i.e., uλ(t) is nondecreasing in λ, limλ→⁰⁺‖uλ‖=0 and limλ→+∞‖uλ‖=+∞.     

Global existence and nonexistence of solution for Cauchy problem of multidimensional double dispersion equations

In this paper we consider the Cauchy problem of multidimensional generalized double dispersion equations utt−Δu−Δutt+Δ²u=Δf(u), where f(u)=ap|u|. By potential well method we prove the existence and nonexistence of global weak solution without establishing the local existence theory. And we derive some sharp conditions for global existence and lack of global existence solution.     

Estimates of solutions of impulsive parabolic equations under Neumann boundary condition

In this paper, we prove the relation v(t)?u(t,x)?w(t), where u(t,x) is the solution of an impulsive parabolic equations under Neumann boundary condition ∂u(t,x)/∂ν=0, and v(t) and w(t) are solutions of two impulsive ordinary equations. We also apply these estimates to investigate the asymptotic behavior of a model in the population dynamics, and it is shown that there exists a unique solution of the model which converges to the periodic solution of an impulsive ordinary equation asymptotically.     

Mixed problems for the time-dependent telegraph equation: Continuous numerical solutions with a priori error bounds

This paper deals with the construction of continuous numerical solutions of mixed problems described by the time-dependent telegraph equation u_tt + c(t)u_t + b(t)u = a(t)u_xx, 0 < x < d, t > 0. Here a(t), b(t), and c(t) are positive functions with appropiate additional alternative hypotheses. First, using the separation of variables technique a theoretical series solution is obtained. Then, after truncation using one-step matrix methods and interpolating functions, a continuous numerical solution with a prefixed accuracy in a bounded subdomain is constructed.     

Globally smooth solutions of quasilinear hyperbolic systems in diagonal form

We prove the existence of a large class of globally smooth solutions of the Cauchy problem for the system of n equations u_t + Λ(x, t, u)u_x = 0, where Λ is a diagonal matrix. We show that, under certain monotonicity conditions on both Λ and the initial data u₀, the solution u will be locally Lipschitz continuous at positive times. Since u₀ is a function of locally bounded variation, our result thus provides both for the smoothing of discontinuities in u₀ as well as for the global preservation of smoothness. The global existence results from an a priori estimate of $\text{?u}\text{?x}$, which we obtain by a device which enables us to effectively uncouple the system of equations for $\text{?u}\text{?x}$. Finally, we prove a partial converse which demonstrates that our hypotheses are not overly restrictive.     

Blow-up and global solutions for a class of nonlinear parabolic equations with different kinds of boundary conditions

By constructing different auxiliary functions and using Hopf’s maximum principle, the sufficient conditions for the blow-up and global solutions are presented for nonlinear parabolic equation u_t = ∇(a(u)b(x)c(t)∇u) + f(x, u, q, t) with different kinds of boundary conditions. The upper bounds of the “blow-up time” and the “upper estimates” of global solutions are provided. Finally, some examples are presented as the application of the obtained results.     

On solutions to the Degasperis-Procesi equation

In this paper we consider a new integrable equation (the Degasperis-Procesi equation) derived recently by Degasperis and Procesi (1999) [3]. Analogous to the Camassa-Holm equation, this new equation admits blow-up phenomenon and infinite propagation speed. First, we give a proof for the blow-up criterion established by Zhou (2004) in [12]. Then, infinite propagation speed for the Degasperis-Procesi equation is proved in the following sense: the corresponding solution u(x,t) with compactly supported initial datum u₀(x) does not have compact x-support any longer in its lifespan. Moreover, we show that for any fixed time t>0 in its lifespan, the corresponding solution u(x,t) behaves as: u(x,t)=L(t)e−x for x?1, and u(x,t)=l(t)ex for x?−1, with a strictly increasing function L(t)>0 and a strictly decreasing function l(t)<0 respectively.     

Locally Lipschitz solutions of a single conservation law in several space variables

We consider the Cauchy problem for a single conservation law in several space variables. Letting u(x, t) denote the solution with initial data u₀, we state necessary and sufficient conditions on u₀ so that u(x, t) is locally Lipschitz continuous in the half space {t > 0}. These conditions allow for the preservation of smoothness of u₀ as well as for the smooth resolution of discontinuities in u₀. One consequence of our result is that u(x, t) cannot be locally Lipschitz unless u₀ has locally bounded variation. Another is that solutions which are bounded and locally Lipschitz continuous in {t > 0} automatically have boundary values u₀ at t = 0 in the sense that u(·, t) → u₀ in L_loc¹. Finally, we give an elementary proof that locally Lipschitz solutions satisfy Kruzkov's uniqueness condition.     

On the generalized Burgers equation

The present paper is concerned with the initial boundary value problem for the generalized Burgers equation u _t + g(t, u)u _x + f(t, u) = εu _xx which arises in many applications. We formulate a condition guaranteeing the a priori estimate of max |u _x | independent of ε and t and give an example demonstrating the optimality of this condition. Based on this estimate we prove the global existence of a unique classical solution of the problem and investigate the behavior of this solution for ε → 0 and t → + ∞. The Cauchy problem for this equation is considered as well.     

Stability of the Absolutely Continuous Spectrum of Random Schrödinger Operators on Tree Graphs

According to the Smolukowski-Kramers approximation, we show that the solution of the semi-linear stochastic damped wave equations μ u _tt (t,x)=Δu(t,x)?u _t (t,x)+b(x,u(t,x))+Q (t),u(0)=u ₀, u _t (0)=v ₀, endowed with Dirichlet boundary conditions, converges as μ goes to zero to the solution of the semi-linear stochastic heat equation u _t (t,x)=Δ u(t,x)+b(x,u(t,x))+Q (t),u(0)=u ₀, endowed with Dirichlet boundary conditions. Moreover we consider relations between asymptotics for the heat and for the wave equation. More precisely we show that in the gradient case the invariant measure of the heat equation coincides with the stationary distributions of the wave equation, for any μ>0.     

Nonlinear abstract differential equations with deviated argument

In this paper we consider the nonlinear differential equation with deviated argument u^′(t)=Au(t)+f(t,u(t),u[φ(u(t),t)]), t∈R₊, in a Banach space (X,‖⋅‖), where A is the infinitesimal generator of an analytic semigroup. Under suitable conditions on the functions f and φ, we prove a global existence and uniqueness result for the above equation.     

Two-sided estimates of a solution of the Neumann problem as t ← ∠ for a second-order Quasilinear Parabolic Equation

We establish exact upper and lower bounds as t ← ∠ for the norm ‖u(·, t)‖ _{L ∞(Ω)} of a solution of the Neumann problem for a second-order quasilinear parabolic equation in the region D=Ω×{>0}, where Ω is a region with noncompact boundary.     

Existence, uniqueness and stability of travelling waves in a discrete reaction-diffusion monostable equation with delay

In this paper, we study the existence, uniqueness and asymptotic stability of travelling wavefronts of the following equation:

u_t(x,t)=D[u(x+1,t)+u(x-1,t)-2u(x,t)]-du(x,t)+b(u(x,t-r)),