Special exact soliton solutions for the K(2, 2) equation with non-zero constant pedestal

Special exact solutions of the K(2, 2) equation, u_t + (u²)_x + (u²)_xxx = 0, are investigated by employing the qualitative theory of differential equations. Our procedure shows that the K(2, 2) equation either has loop soliton, cusped soliton and smooth soliton solutions when sitting on the non-zero constant pedestal lim_x→±∞u = A ≠ 0, or possesses compacton solutions only when lim_x→±∞u = 0. Mathematical analysis and numerical simulations are provided for these soliton solutions of the K(2, 2) equation.     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

Control parameter estimation in a semi-linear parabolic inverse problem using a high accurate method

Parabolic inverse problems have an important role in many branches of science and technology. The aim of this research work is to solve these classes of equations using a high order compact finite difference scheme. We consider the following inverse problem for finding u(x, t) and p(t) governed by u_t = u_xx + p(t)u + φ(x, t) with an over specified condition inside the domain. Spatial derivatives are approximated using central difference scheme. The time advancement of the simulation is performed using a “third order compact Runge-Kutta method”. The convergence orders for the approximation of both u and p are of o(k³ + h²) which improves the results obtained in the literature. An exact test case is used to evaluate the validity of our numerical analysis. We found that the accuracy of the results is better than that of previous works in the literature.     

Generalized three-point difference schemes of high-order accuracy for systems of second-order nonlinear ordinary differential equations

For systems of second-order nonlinear ordinary differential equations with the Dirichlet boundary conditions, we develop generalized three-point difference schemes of high-order accuracy on a nonuniform grid. The construction of the suggested schemes requires solving four auxiliary Cauchy problems (two problems for systems of nonlinear ordinary differential equations and two problems for matrix linear ordinary differential equations) on the intervals [x _j−1, x _j ] (forward) and [x _j , x _j+1] (backward) at each grid point; this is done at each step by any single-step method of accuracy order $ \bar m $ \bar m = 2[(m+1)/2]. (Here m is a given positive integer, and [·] is the integer part of a number.) We prove that such three-point difference schemes have the accuracy order $ \bar m $ \bar m for the approximation to both the solution u of the boundary value problem and the flux K(x)d u/dx at the grid points.     

Critical Exponents of Quasilinear Parabolic Equations

In this paper we study the critical exponents of the Cauchy problem in Rⁿ of the quasilinear singular parabolic equations: u_t = div(|∇u|^m − 1∇u) + t^s|x|^σu^p, with non-negative initial data. Here s ≥ 0, (n − 1)/(n + 1) < m < 1, p > 1 and σ > n(1 − m) − (1 + m + 2s). We prove that p_c ≡ m + (1 + m + 2s + σ)/n > 1 is the critical exponent. That is, if 1 < p ≤ p_c then every non-trivial solution blows up in finite time, but for p > p_c, a small positive global solution exists.     

The Legendre-Burgers equation: When artificial dissipation fails

Artificial viscosity is a common device for stabilizing flows with shocks and fronts. The computational diffusion smears the frontal zone over a small distance μ where μ is chosen so that the discretization has a couple of grid points in the front, and thus is able to resolve the shock. Spectral element methods use a Legendre spectral viscosity whose effect is to damp the coefficient of P_n(x) by some amount that depends only on the degree n of the Legendre polynomial. Legendre viscosity is better than ordinary diffusion because it does not require spurious boundary conditions, does not increase the temporal stiffness of the differential equations, and can be applied locally on an element-by-element basis. Unfortunately, Legendre diffusion is equivalent to a diffusion with a spatially-varying coefficient that goes to zero at the boundaries. Using the simplest example, one in which the second derivative of Burgers equation is replaced by the Legendre operator to give the “Legendre-Burgers” equation, u_t + uu_x = ν[(1 − x²)u_x]_x, we show that the width of the computational front can similarly tend to zero at the endpoints, causing a numerical catastrophe.     

The improved linear multistep methods for differential equations with piecewise continuous arguments

This paper deals with the convergence of the linear multistep methods for the equation x′(t) = ax(t) + a₀x([t]). Numerical experiments demonstrate that the 2-step Adams-Bashforth method is only of order p = 0 when applied to the given equation. An improved linear multistep methods is constructed. It is proved that these methods preserve their original convergence order for ordinary differential equations (ODEs) and some numerical experiments are given.     

On some special codes over Fp + uFp + uFp

We define a new map between codes over F_p + uF_p + u²F_p and F_p which is different to that defined in [2]. It is proved that the image of the linear cyclic code over the commutative ring F_p + uF_p + u²F_p with length n under this map is a distance-invariant quasi-cyclic code of index p² with length p²n over F_p. Moreover, it is proved that, if (n, p) = 1, then every code with length p²n over F_p which is the image of a linear (1 − u²)-cyclic code with length n over F_p + uF_p + u²F_p under this map is permutation equivalent to a quasi-cyclic code of index p².     

A finite difference solution to a two-dimensional parabolic inverse problem

In this paper, we consider the problem of finding u = u(x, y, t) and p = p(t) which satisfy u_t = u_xx + u_yy + p(t)u + ? in R × [0, T], u(x, y, 0) = f(x, y), (x, y) ∈ R = [0, 1] × [0, 1], u is known on the boundary of R and u(x^∗, y^∗, t) = E(t), 0 < t ? T, where E(t) is known and (x^∗, y^∗) is a given point of R. Through a function transformation, the nonlinear two-dimensional diffusion problem is transformed into a linear problem, and a backward Euler scheme is constructed. It is proved by the maximum principle that the scheme is uniquely solvable, unconditionally stable and convergent in L_∞ norm. The convergence orders of u and p are of O(τ + h²). The impact of initial data errors on the numerical solution is also considered. Numerical experiments are presented to illustrate the validity of the theoretical results.     

Peakon, loop and solitary travelling wave solutions for the general modified CH-DP equation

Travelling wave solutions for the general modified CH-DP equation u_t − u_xxt + αu²u_x − βu_xu_xx = uu_xxx are developed. By using the dynamical system method, a peakon and a dark soliton are found to coexist for the same wave speed. Exact explicit blow-up solutions are given. By using numerical simulation, a loop solution for a special case is discussed.     

Differential invariants of the one-dimensional quasi-linear second-order evolution equation

We consider evolution equations of the form u_t = f(x, u, u_x)u_xx + g(x, u, u_x) and u_t = u_xx + g(x, u, u_x). In the spirit of the recent work of Ibragimov [Ibragimov NH. Laplace type invariants for parabolic equations. Nonlinear Dynam 2002;28:125–33] who adopted the infinitesimal method for calculating invariants of families of differential equations using the equivalence groups, we apply the method to these equations. We show that the first class admits one differential invariant of order two, while the second class admits three functional independent differential invariants of order three. We use these invariants to determine equations that can be transformed into the linear diffusion equation.     

Eigenvalues and the One-Dimensional p-Laplacian

We consider the boundary value problem (?_p(u′))′ + λF(t, u) = 0, with p > 1, t ∈ (0, 1), u(0) = u(1) = 0, and with λ > 0. The value of λ is chosen so that the boundary value problem has a positive solution. In addition, we derive an explicit interval for λ such that, for any λ in this interval, the existence of a positive solution to the boundary value problem is guaranteed. In addition, the existence of two positive solutions for λ in an appropriate interval is also discussed.     

An efficient algorithm for the multivariable Adomian polynomials

In this article the sum of the series of multivariable Adomian polynomials is demonstrated to be identical to a rearrangement of the multivariable Taylor expansion of an analytic function of the decomposition series of solutions u₁, u₂, … , u_m about the initial solution components u_1,0, u_2,0, … , u_m,0; of course the multivariable Adomian polynomials were developed and are eminently practical for the solution of coupled nonlinear differential equations. The index matrices and their simplified forms of the multivariable Adomian polynomials are introduced. We obtain the recurrence relations for the simplified index matrices, which provide a convenient algorithm for rapid generation of the multivariable Adomian polynomials. Another alternative algorithm for term recurrence is established. In these algorithms recurrence processes do not require complicated operations such as parametrization, expanding and regrouping, derivatives, etc. as practiced in prior art. The MATHEMATICA program generating the Adomian polynomials based on the algorithm in this article is designed.     

Properties of non-simultaneous blow-up in heat equations coupled via different localized sources

This paper deals with u_t = Δu + u^m(x, t)e^pv(0,t), v_t = Δv + u^q(0, t)e^nv(x,t), subject to homogeneous Dirichlet boundary conditions. The complete classification on non-simultaneous and simultaneous blow-up is obtained by four sufficient and necessary conditions. It is interesting that, in some exponent region, large initial data u₀(v₀) leads to the blow-up of u(v), and in some betweenness, simultaneous blow-up occurs. For all of the nonnegative exponents, we find that u(v) blows up only at a single point if m > 1(n > 0), while u(v) blows up everywhere for 0 ? m ? 1 (n = 0). Moreover, blow-up rates are considered for both non-simultaneous and simultaneous blow-up solutions.     

One-point pseudospectral collocation for the one-dimensional Bratu equation

The one-dimensional planar Bratu problem is u_xx + λ exp(u) = 0 subject to u(±1) = 0. Because there is an analytical solution, this problem has been widely used to test numerical and perturbative schemes. We show that over the entire lower branch, and most of the upper branch, the solution is well approximated by a parabola, u(x) ≈ u₀ (1 − x²) where u₀ is determined by collocation at a single point x = ξ. The collocation equation can be solved explicitly in terms of the Lambert W-function as u(0) ≈ −W(−λ(1 − ξ²)/2)/(1 − ξ²) where both real-valued branches of the W-function yield good approximations to the two branches of the Bratu function. We carefully analyze the consequences of the choice of ξ. We also analyze the rate of convergence of a series of even Chebyshev polynomials which extends the one-point approximation to arbitrary accuracy. The Bratu function is so smooth that it is actually poor for comparing methods because even a bad, inefficient algorithm is successful. It is, however, a solution so smooth that a numerical scheme (the collocation or pseudospectral method) yields an explicit, analytical approximation. We also fill some gaps in theory of the Bratu equation. We prove that the general solution can be written in terms of a single, parameter-free β(x) without knowledge of the explicit solution. The analytical solution can only be evaluated by solving a transcendental eigenrelation whose solution is not known explicitly. We give three overlapping perturbative approximations to the eigenrelation, allowing the analytical solution to be easily evaluated throughout the entire parameter space.     

A Variational Approach to Discontinuous Problems with Critical Sobolev Exponents

We employ variational techniques to study the existence and multiplicity of positive solutions of semilinear equations of the form − Δu = λh(x)H(u − a)u^q + u^2* − 1 in R^N, where λ, a > 0 are parameters, h(x) is both nonnegative and integrable on R^N, H is the Heaviside function, 2* is the critical Sobolev exponent, and 0 ≤ q < 2* − 1. We obtain existence, multiplicity and regularity of solutions by distinguishing the cases 0 ≤ q ≤ 1 and 1 < q < 2* − 1.     

Nonoscillation and oscillation of second order half-linear differential equations

We study the oscillation problems for the second order half-linear differential equation ^′[p(t)Φ(x^′)]+q(t)Φ(x)=0, where Φ(u)=|u|^r−1u with r>0, 1/p and q are locally integrable on R₊; p>0, q?0 a.e. on R₊, and . We establish new criteria for this equation to be nonoscillatory and oscillatory, respectively. When p≡1, our results are complete extensions of work by Huang [C. Huang, Oscillation and nonoscillation for second order linear differential equations, J. Math. Anal. Appl. 210 (1997) 712-723] and by Wong [J.S.W. Wong, Remarks on a paper of C. Huang, J. Math. Anal. Appl. 291 (2004) 180-188] on linear equations to the half-linear case for all r>0. These results provide corrections to the wrongly established results in [J. Jiang, Oscillation and nonoscillation for second order quasilinear differential equations, Math. Sci. Res. Hot-Line 4 (6) (2000) 39-47] on nonoscillation when 0<r<1 and on oscillation when r>1. The approach in this paper can also be used to fully extend Elbert's criteria on linear equations to half-linear equations which will cover and improve a partial extension by Yang [X. Yang, Oscillation/nonoscillation criteria for quasilinear differential equations, J. Math. Anal. Appl. 298 (2004) 363-373].