SPACE-TIME DISCONTINUOUS GALERKIN METHOD FOR MAXWELL EQUATIONS IN DISPERSIVE MEDIA

In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the underlying problem. Unconditional L2-stability and error estimate of order O τr+1+ hk+1/2 are obtained when polynomials of degree at most r and k are used for the temporal discretization and spatial discretization respectively. 2-D and 3-D numerical examples are given to validate the theoretical results. Moreover, numerical results show an ultra-convergence of order 2r + 1 in temporal variable t.     

FULL DISCRETE TWO-LEVEL CORRECTION SCHEME FOR NAVIER-STOKES EQUATIONS

In this paper, a full discrete two-level scheme for the unsteady Navier-Stokes equations based on a time dependent projection approach is proposed. In the sense of the new projection and its related space splitting, non-linearity is treated only on the coarse level subspace at each time step by solving exactly the standard Galerkin equation while a linear equation has to be solved on the fine level subspace to get the final approximation at this time step. Thus, it is a two-level based correction scheme for the standard Galerkin approximation. Stability and error estimate for this scheme are investigated in the paper.     

GLOBAL ATTRACTIVITY OF A DELAYED RATIO-DEPENDENT COOPERATIVE SYSTEM WITH STAGE STRUCTURE

An autonomous stage-structured ratio-dependent cooperative system, which was proposed by Muhammadhaji, Teng and Abdurahman, is revisited in this paper. By introducing a new lemma and using the iterative method, a set of sufficient conditions which guarantee the global attractivity of the positive equilibrium is obtained. It is shown that the conditions which guarantee the permanence of the system are enough to ensure the global attractivity of the system. Our result not only complements but also supplements one of the main results of Muhammadhaji, Teng and Abdurahman （Permanence and extinction analysis for a delayed ratio-dependent cooperative system with stage structure, Afr. Mat.（2013）DOI 10.1007/s13370-013-0162-6）.     

ON THE OPTIMIZATION OF EXTRAPOLATION METHODS FOR SINGULAR LINEAR SYSTEMS

We discuss semiconvergence of the extrapolated iterative methods for solving singular linear systems. We obtain the upper bounds and the optimum convergence factor of the extrapolation method as well as its associated optimum extrapolation parameter. Numerical examples are given to illustrate the theoretical results.     

A DIAGONAL SPLIT-CELL MODEL FOR THE OVERLAPPING YEE FDTD METHOD

In this paper, we present a nonorthogonal overlapping Yee method for solv- ing Maxwell＇s equations using the diagonal split-cell model. When material interface is presented, the diagonal split-cell model does not require permittivity averaging so that better accuracy can be achieved. Our numerical results on optical force computation show that the standard FDTD method converges linearly, while the proposed method achieves quadratic convergence and better accuracy.     

ERROR ANALYSIS FOR A FAST NUMERICAL METHOD TO A BOUNDARY INTEGRAL EQUATION OF THE FIRST KIND

For two-dimensional boundary integral equations of the first kind with logarithmic kernels, the use of the conventional boundary element methods gives linear systems with dense matrix. In a recent work [J. Comput. Math., 22 （2004）, pp. 287-298], it is demonstrated that the dense matrix can be replaced by a sparse one if appropriate graded meshes are used in the quadrature rules. The numerical experiments also indicate that the proposed numerical methods require less computational time than the conventional ones while the formal rate of convergence can be preserved. The purpose of this work is to establish a stability and convergence theory for this fast numerical method. The stability analysis depends on a decomposition of the coefficient matrix for the collocation equation. The formal orders of convergence observed in the numerical experiments are proved rigorously.     

NUMERICAL LOCALIZATION OF ELECTROMAGNETIC IMPERFECTIONS FROM A PERTURBATION FORMULA IN THREE DIMENSIONS

This work deals with the numerical localization of small electromagnetic inhomogeneities. The underlying inverse problem considers, in a three-dimensional bounded domain, the time-harmonic Maxwell equations formulated in electric field. Typically, the domain contains a finite number of unknown inhomogeneities of small volume and the inverse problem attempts to localize these inhomogeneities from a finite number of boundary measurements. Our localization approach is based on a recent framework that uses an asymptotic expansion for the perturbations in the tangential boundary trace of the curl of the electric field. We present three numerical localization procedures resulting from the combination of this asymptotic expansion with each of the following inversion algorithms： the Current Projection method, the MUltiple Signal Classification （MUSIC） algorithm, and an Inverse Fourier method. We perform a numerical study of the asymptotic expansion and compare the numerical results obtained from the three localization procedures in different settings.     

A NEW NONMONOTONE TRUST REGION ALGORITHM FOR SOLVING UNCONSTRAINED OPTIMIZATION PROBLEMS

Based on the nonmonotone line search technique proposed by Gu and Mo （Appl. Math. Comput. 55, （2008） pp. 2158-2172）, a new nonmonotone trust region algorithm is proposed for solving unconstrained optimization problems in this paper. The new algorithm is developed by resetting the ratio ρk for evaluating the trial step dk whenever acceptable. The global and superlinear convergence of the algorithm are proved under suitable conditions. Numerical results show that the new algorithm is effective for solving unconstrained optimization problems.     

The exceptional set for sums of unlike powers of primes

It is established that all even positive integers up to N but at most O(N15/16+ε) exceptions can be expressed in the form p21+ p32+ p43+ p54,where p1,p2,p3 and p4 are prime numbers.     

ASYMPTOTIC DECAY TOWARD RAREFACTION WAVE FOR A HYPERBOLIC-ELLIPTIC COUPLED SYSTEM ON HALF SPACE

We consider the asymptotic behavior of solutions to a model of hyperbolicelliptic coupled system on the half-line R＋ = （0, ∞）,with the Dirichlet boundary condition u（0, t） = 0. S. Kawashima and Y. Tanaka [Kyushu J. Math., 58（2004）, 211-250] have shown that the solution to the corresponding Cauchy problem behaviors like rarefaction waves and obtained its convergence rate when u_ u＋. Our main concern in this paper is the boundary effect. In the case of null-Dirichlet boundary condition on u, asymptotic behavior of the solution （u, q） is proved to be rarefaction wave as t tends to infinity. Its convergence rate is also obtained by the standard L2-energy method and Ll-estimate. It decays much lower than that of the corresponding Cauchy problem.     

THE RESTRICTIVELY PRECONDITIONED CONJUGATE GRADIENT METHODS ON NORMAL RESIDUAL FOR BLOCK TWO-BY-TWO LINEAR SYSTEMS

The restrictively preconditioned conjugate gradient （RPCG） method is further developed to solve large sparse system of linear equations of a block two-by-two structure. The basic idea of this new approach is that we apply the RPCG method to the normal-residual equation of the block two-by-two linear system and construct each required approximate matrix by making use of the incomplete orthogonal factorization of the involved matrix blocks. Numerical experiments show that the new method, called the restrictively preconditioned conjugate gradient on normal residual （RPCGNR）, is more robust and effective than either the known RPCG method or the standard conjugate gradient on normal residual （CGNR） method when being used for solving the large sparse saddle point problems.     

Bifurcations and chaos in Duffing equation with damping and external excitations

Duffng equation with damping and external excitations is investigated.By using Melnikov method and bifurcation theory,the criterions of existence of chaos under periodic perturbations are obtained.By using second-order averaging method,the criterions of existence of chaos in averaged system under quasi-periodic perturbations forff=nω+εσ,n=2,4,6(whereσis not rational toω)are investigated.However,the criterions of existence of chaos for n=1,3,5,7-20 can not be given.The numerical simulations verify the theoretical analysis,show the occurrence of chaos in the averaged system and original system under quasiperiodic perturbation for n=1,2,3,5,and expose some new complex dynamical behaviors which can not be given by theoretical analysis.In particular,the dynamical behaviors under quasi-periodic perturbations are different from that under periodic perturbations,and the period-doubling bifurcations to chaos has not been found under quasi-periodic perturbations.     

HOLDER ESTIMATES FOR A CLASS OF DEGENERATE ELLIPTIC EQUATIONS

In this paper we study the regularity theory of the solutions of a class of degenerate elliptic equations in divergence form. By introducing a proper distance and applying the compactness method we establish the HSlder type estimates for the weak solutions.     

A GHOST FLUID BASED FRONT TRACKING METHOD FOR MULTIMEDIUM COMPRESSIBLE FLOWS

Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.     

CONVERGENCE ANALYSIS OF THE JACOBI SPECTRAL-COLLOCATION METHOD FOR FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS

We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method,which shows that the errors of the approximate solution decay exponentially in L∞norm and weighted L2-norm. The numerical examples are given to illustrate the theoretical results.     

A REACTIVE DYNAMIC CONTINUUM USER EQUILIBRIUM MODEL FOR BI-DIRECTIONAL PEDESTRIAN FLOWS

In this paper, a reactive dynamic user equilibrium model is extended to simulate two groups of pedestrians traveling on crossing paths in a continuous walking facility. Each group makes path choices to minimize the travel cost to its destination in a reactive manner based on instantaneous information. The model consists of a conservation law equation coupled with an Eikonal-type equation for each group. The velocity-density relationship of pedestrian movement is obtained via an experimental method. The model is solved using a finite volume method for the conservation law equation and a fast-marching method for the Eikonal-type equation on unstructured grids. The numerical results verify the rationality of the model and the validity of the numerical method. Based on this continuum model, a number of results, e.g., the formation of strips or moving clusters composed of pedestrians walking to the same destination, are also observed.     

A REDUCED-ORDER MFE FORMULATION BASED ON POD METHOD FOR PARABOLIC EQUATIONS

In this paper, we extend the applications of proper orthogonal decomposition (POD) method, i.e., apply POD method to a mixed finite element (MFE) formulation naturally satisfied Brezz-Babu?ska for para...     

INFINITELY MANY SOLUTIONS FOR AN ELLIPTIC PROBLEM INVOLVING CRITICAL NONLINEARITY

We study the following elliptic problem:{-div(a(x)Du)=Q(x)|u|2-2u+λu x ∈Ω ,u=0 on ΩUnder certain assumptions on a and Q, we obtain existence of infinitely many solutions by variational method.     

CHEBYSHEV METHODS WITH DISCRETE NOISE： THE τ-ROCK METHODS

Stabilized or Chebyshev explicit methods have been widely used in the past to solve stiff ordinary differential equations. Making use of special properties of Chebyshev-like polynomials, these methods have favorable stability properties compared to standard explicit methods while remaining explicit. A new class of such methods, called ROCK, introduced in [Numer. Math., 90, 1-18, 2001] has recently been extended to stiff stochastic differential equations under the name S-ROCK [C. R. Acad. Sci. Paris, 345（10）, 2007 and Commun. Math. Sci, 6（4）, 2008]. In this paper we discuss the extension of the S-ROCK methods to systems with discrete noise and propose a new class of methods for such problems, the τ-ROCK methods. One motivation for such methods is the simulation of multi-scale or stiff chemical kinetic systems and such systems are the focus of this paper, but our new methods could potentially be interesting for other stiff systems with discrete noise. Two versions of the τ-ROCK methods are discussed and their stability behavior is analyzed on a test problem. Compared to the τ-leaping method, a significant speed-up can be achieved for some stiff kinetic systems. The behavior of the proposed methods are tested on several numerical experiments.     

The Tanh Method for Kink Solution of Some Modified Nonlinear Equation

The tanh method is a very powerful technique for computation of exact traveling wave, in this paper this method has been employed for special modified states of Burger, Klein-Gordon and Fisher-Burger equations and the solitary solution of these equations are derived.