New generalized Jacobi elliptic function rational expansion method

In this work, a new generalized Jacobi elliptic function rational expansion method is based upon twenty-four Jacobi elliptic functions and eight double periodic Weierstrass elliptic functions, which solve the elliptic equation ?^′2=r+p?²+q?⁴, is described. As a consequence abundant new Jacobi-Weierstrass double periodic elliptic functions solutions for (3+1)-dimensional Kadmtsev-Petviashvili (KP) equation are obtained by using this method. We show that the new method can be also used to solve other nonlinear partial differential equations (NPDEs) in mathematical physics.     

Cubic superconvergent finite volume element method for one-dimensional elliptic and parabolic equations

In this paper, a cubic superconvergent finite volume element method based on optimal stress points is presented for one-dimensional elliptic and parabolic equations. For elliptic problem, it is proved that the method has optimal third order accuracy with respect to H¹ norm and fourth order accuracy with respect to L² norm. We also obtain that the scheme has fourth order superconvergence for derivatives at optimal stress points. For parabolic problem, the scheme is given and error estimate is obtained with respect to L² norm. Finally, numerical examples are provided to show the effectiveness of the method.     

Average implicit linear difference scheme for generalized Rosenau-Burgers equation

In this paper, a linear three-level average implicit finite difference scheme for the numerical solution of the initial-boundary value problem of Generalized Rosenau-Burgers equation is presented. Existence and uniqueness of numerical solutions are discussed. It is proved that the finite difference scheme is convergent in the order of O(τ² + h²) and stable. Numerical simulations show that the method is efficient.     

The method of boundary integral equations for a mixed scattering problem

In this paper we consider the scattering of an electromagnetic time-harmonic plane wave by an infinite cylinder having an open arc and a bounded domain in R² as cross section. To this end, we solve a scattering problem for the Helmholtz equation in R² where the scattering object is a combination of a crack Γ and a bounded obstacle D, and we have Dirichlet-impedance type boundary condition on Γ and Dirichlet boundary condition on ∂D (∂D∈C²). Applying potential theory, the problem can be reformulated as a boundary integral system. We establish the existence and uniqueness of a solution to the system by using the Fredholm theory.     

On reflection principles supported on a finite set

We investigate the existence of reflection formulas supported on a finite set. It is found that for solutions of the Laplace and Helmholtz equation there are no finitely supported reflection principles unless the support is a single point. This confirms that in order to construct a reflection formula that is not ‘point to point’, it is necessary to consider a continuous support. For solutions of the wave equation ^∂2u/∂x∂y=0, there exist finitely supported reflection principles that can be constructed explicitly. For solutions of the telegraph equation ^∂2u/∂x∂y+λ²u=0, we show that if a reflection principle is supported on less than five points then it is a point to point reflection principle.     

Estimates for wave and Klein-Gordon equations on modulation spaces

We prove that the fundamental semi-group eit(m 2I+|Δ|)1/2(m = 0) of the Klein-Gordon equation is bounded on the modulation space M ps,q(Rn) for all 0 < p,q ∞ and s ∈ R.Similarly,we prove that the wave semi-group eit|Δ|1/2 is bounded on the Hardy type modulation spaces μsp,q(Rn) for all 0 < p,q ∞,and s ∈ R.All the bounds have an asymptotic factor tn|1/p 1/2| as t goes to the infinity.These results extend some known results for the case of p 1.Also,some applications for the Cauchy problems related to the semi-group eit(m2I+|Δ|)1/2 are obtained.Finally we discuss the optimum of the factor tn|1/p 1/2| and raise some unsolved problems.     

The cubic fourth-order Schrödinger equation

Fourth-order Schrödinger equations have been introduced by Karpman and Shagalov to take into account the role of small fourth-order dispersion terms in the propagation of intense laser beams in a bulk medium with Kerr nonlinearity. In this paper we investigate the cubic defocusing fourth-order Schrödinger equation

it_∂u+Δ²u+²|u|u=0     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.     

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.     

On the use of high order difference methods for the system of one space second order nonlinear hyperbolic equations with variable coefficients

Implicit difference schemes of O(k⁴ + k²h² + h⁴), where k0, h 0 are grid sizes in time and space coordinates respectively, are developed for the efficient numerical integration of the system of one space second order nonlinear hyperbolic equations with variable coefficients subject to appropriate initial and Dirichlet boundary conditions. The proposed difference method for a scalar equation is applied for the wave equation in cylindrical and spherical symmetry. The numerical examples are given to illustrate the fourth order convergence of the methods.     

A new high accuracy non-polynomial tension spline method for the solution of one dimensional wave equation in polar co-ordinates

In this paper, we propose a new three-level implicit nine point compact finite difference formulation of O(k² + h⁴) based on non-polynomial tension spline approximation in r-direction and finite difference approximation in t-direction for the numerical solution of one dimensional wave equation in polar co-ordinates. We describe the mathematical formulation procedure in details and also discuss the stability of the method. Numerical results are provided to justify the usefulness of the proposed method.     

A uniformly convergent continuous-discontinuous Galerkin method for singularly perturbed problems of convection-diffusion type

In this paper, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving singularly perturbed convection-diffusion problems. On Shishkin mesh with linear elements, a rate O(N^-1lnN) in an associated norm is established, where N is the number of elements. Numerical experiments complement the theoretical results. Moreover, a rate O(N^-2ln²N) in a discrete L^∞ norm, and O(N^-2) in L² norm, are observed numerically on the Shishkin mesh.     

Trigonometric wavelet method for some elliptic boundary value problems

In this paper, we apply trigonometric wavelet method to investigate the numerical solution of elliptic boundary value problem in an exterior unit disk. The simple computation formulae of the entries in the stiffness matrix are obtained. It shows that we only need to compute 2(j²+1) elements of one ²j+2ײj+2 stiffness matrix. Moreover, the error estimates of the approximation solutions are given and some test examples are presented.     

Lie symmetries for the orbital linearization of smooth planar vector fields around singular points

This work is a generalization of the method proposed in [I.A. García, S. Maza, Linearization of analytic isochronous centers from a given commutator, J. Math. Anal. Appl. 339 (1) (2008) 740-745] of linearization of analytic isochronous centers from a given commutator. In this paper we propose a constructive procedure to get the change of variables that orbitally linearizes a smooth planar vector field on C² around an elementary singular point (i.e., a singular point with associated eigenvalues λ,μ∈C satisfying λ≠0) or a nilpotent singular point from a given infinitesimal generator of a Lie symmetry.     

On non-autonomous strongly damped wave equations with a uniform attractor and some averaging

In this paper we study the existence of a uniform attractor for strongly damped wave equations with a time-dependent driving force. If the time-dependent function is translation compact, then in a certain parameter region, the uniform attractor of the system has a simple structure: it is the closure of all the values of the unique, bounded complete trajectory of the wave equation. And it attracts any bounded set exponentially. At the same time, we consider the strongly damped wave equations with rapidly oscillating external force gε(x,t)=g(x,t,t/ε) having the average g⁰(x,t) as ε→⁰⁺. We prove that the Hausdorff distance between the uniform attractor Aε of the original equation and the uniform attractor A₀ of the averaged equation is less than O(ε1/2). We mention, in particular, that the obtained results can be used to study the usual damped wave equations.     

The behavior of solutions of second order delay differential equations

In this paper, we study the behavior of solutions of second order delay differential equation

y^″(t)=p₁y^′(t)+p₂y^′(t−τ)+q₁y(t)+q₂y(t−τ),     

Solutions of polynomial Pell's equation

Let D=F²+2G be a monic quartic polynomial in Z[x], where . Then for F/G∈Q[x], a necessary and sufficient condition for the solution of the polynomial Pell's equation X²−DY²=1 in Z[x] has been shown. Also, the polynomial Pell's equation X²−DY²=1 has nontrivial solutions X,Y∈Q[x] if and only if the values of period of the continued fraction of are 2, 4, 6, 8, 10, 14, 18, and 22 has been shown. In this paper, for the period of the continued fraction of is 4, we show that the polynomial Pell's equation has no nontrivial solutions X,Y∈Z[x].     

A second order finite difference-spectral method for space fractional diffusion equations

A high order finite difference-spectral method is derived for solving space fractional diffusion equations,by combining the second order finite difference method in time and the spectral Galerkin method in space.The stability and error estimates of the temporal semidiscrete scheme are rigorously discussed,and the convergence order of the proposed method is proved to be O(τ2+Nα-m)in L2-norm,whereτ,N,αand m are the time step size,polynomial degree,fractional derivative index and regularity of the exact solution,respectively.Numerical experiments are carried out to demonstrate the theoretical analysis.     

Evolution inclusions governed by the difference of two subdifferentials in reflexive Banach spaces

The existence of strong solutions of Cauchy problem for the following evolution equation du(t)/dt+∂?¹(u(t))-∂?²(u(t))∋f(t) is considered in a real reflexive Banach space V, where ∂?¹ and ∂?² are subdifferential operators from V into its dual V^*. The study for this type of problems has been done by several authors in the Hilbert space setting.The scope of our study is extended to the V-V^* setting. The main tool employed here is a certain approximation argument in a Hilbert space and for this purpose we need to assume that there exists a Hilbert space H such that V⊂H≡H^*⊂V^* with densely defined continuous injections.The applicability of our abstract framework will be exemplified in discussing the existence of solutions for the nonlinear heat equation: where Ω is a bounded domain in R^N. In particular, the existence of local (in time) weak solution is shown under the subcritical growth condition q<p^* (Sobolev's critical exponent) for all initial data This fact has been conjectured but left as an open problem through many years.     

The Smoothness of Scattering Operators for Sinh-Gordon and Nonlinear Schrödinger Equations

We show that the scattering operator carries a band in H ^s (? ⁿ ) ×H ^s?1 (? ⁿ ) into H ^s (? ⁿ ) ×H ^s?1(? ⁿ ) for the sinh-Gordon equation and an analogous result also holds true for the nonlinear Schrödinger equation with an exponential nonlinearity, where s≥n/2 is arbitrary and n≥ 2. Therefore, the scattering operators are infinitely smooth for the above two equations.