CONSTRUCTION OF GEOMETRIC PARTIAL DIFFERENTIAL EQUATIONS FOR LEVEL SETS

Geometric partial differential equations of level-set form are usually constructed by a variational method using either Dirac delta function or co-area formula in the energy functional to be minimized. However, the equations derived by these two approaches are not consistent. In this paper, we present a third approach for constructing the level-set form equations. By representing various differential geometry quantities and differential geometry operators in terms of the implicit surface, we are able to reformulate three classes of parametric geometric partial differential equations （second-order, fourth-order and sixth- order） into the level-set forms. The reformulation of the equations is generic and simple, and the resulting equations are consistent with their parametric form counterparts. We further prove that the equations derived using co-area formula are also consistent with the parametric forms. However, these equations are of much complicated forms than these given by the equations we derived.     

EFFICIENT NUMERICAL ALGORITHMS FOR THREE-DIMENSIONAL FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS

This paper detailedly discusses the locally one-dimensional numerical methods for ef- ficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional diffusion equation. The second order finite difference scheme is used to discretize the space fractional derivative and the Crank-Nicolson procedure to the time derivative. We theoretically prove and numerically verify that the presented numerical methods are unconditionally stable and second order convergent in both space and time directions. In particular, for the Riesz fractional dif- fusion equation, the idea of reducing the splitting error is used to further improve the algorithm, and the unconditional stability and convergency are also strictly proved and numerically verified for the improved scheme.     

HOMOTOPY CONTINUATION METHODS FOR STOCHASTIC TWO-POINT BOUNDARY VALUE PROBLEMS DRIVEN BY ADDITIVE NOISES

In this paper the homotopy continuation method for stochastic two-point boundary value problems driven by additive noises is studied. The existence of the solution of the homotopy equation is proved. Numerical schemes are constructed and error estimates are obtained. Numerical experiments demonstrate the effectiveness of the homotopy continu- ation method over other commonly used methods such as the shooting method.     

REDUCED BASIS METHOD FOR PARAMETRIZED ELLIPTIC ADVECTION-REACTION PROBLEMS

In this work we consider the Reduced Basis method for the solution of parametrized advection-reaction partial differential equations. For the generation of the basis we adopt a stabilized finite element method and we define the Reduced Basis method in the ＂primal- dual＂ formulation for this stabilized problem. We provide a priori Reduced Basis error estimates and we discuss the effects of the finite element approximation on the Reduced Basis error. We propose an adaptive algorithm, based on the a posteriori Reduced Basis error estimate, for the selection of the sample sets upon which the basis are built; the idea leading this algorithm is the minimization of the computational costs associated with the solution of the Reduced Basis problem. Numerical tests demonstrate the efficiency, in terms of computational costs, of the ＂primal-dual＂ Reduced Basis approach with respect to an ＂only primal＂ one. Parametrized advection-reaction partial differential equations, Reduced Basis method, ＂primal-dual＂ reduced basis approach, Stabilized finite element method, a posteriori error estimation.     

A SPECTRAL METHOD FOR PANTOGRAPH-TYPE DELAY DIFFERENTIAL EQUATIONS AND ITS CONVERGENCE ANALYSIS

We propose a novel numerical approach for delay differential equations with vanishing proportional delays based on spectral methods. A Legendre-collocation method is employed to obtain highly accurate numerical approximations to the exact solution. It is proved theoretically and demonstrated numerically that the proposed method converges exponentially provided that the data in the are smooth. given pantograph delay differential equation     

NUMERICAL RECONSTRUCTION OF SMALL PERTURBATIONS IN THE ELECTROMAGNETIC COEFFICIENTS OF A DIELECTRIC MATERIAL

The aim of this paper is to solve numerically the inverse problem of reconstructing small amplitude perturbations in the magnetic permeability of a dielectric material from partial or total dynamic boundary measurements. Our numerical algorithm is based on the resolution of the time-dependent Maxwell equations, an exact controllability method and Fourier inversion for localizing the perturbations. Two-dimensional numerical experiments illustrate the performance of the reconstruction method for different configurations even in the case of limited-view data.     

OPTIMAL POINT-WISE ERROR ESTIMATE OF A COMPACT FINITE DIFFERENCE SCHEME FOR THE COUPLED NONLINEAR SCHRODINGER EQUATIONS

In this paper, we analyze a compact finite difference scheme for computing a coupled nonlinear SchrSdinger equation. The proposed scheme not only conserves the totM mass and energy in the discrete level but also is decoupled and linearized in practical computa- tion. Due to the difficulty caused by compact difference on the nonlinear term, it is very hard to obtain the optimal error estimate without any restriction on the grid ratio. In order to overcome the difficulty, we transform the compact difference scheme into a special and equivalent vector form, then use the energy method and some important lemmas to obtain the optimal convergent rate, without any restriction on the grid ratio, at the order of O（h4 ＋r2） in the discrete L∞ -norm with time step - and mesh size h. Finally, numerical results are reported to test our theoretical results of the proposed scheme.     

A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS

We propose an effective stopping criterion for higher-order fast sweeping schemes for static Hamilton-Jacobi equations based on ratios of three consecutive iterations. To design the new stopping criterion we analyze the convergence of the first-order Lax-Friedrichs sweeping scheme by using the theory of nonlinear iteration. In addition, we propose a fifth-order Weighted PowerENO sweeping scheme for static Hamilton-Jacobi equations with convex Hamiltonians and present numerical examples that validate the effectiveness of the new stopping criterion.     

A NUMERICAL STUDY OF THE GAUSSIAN BEAM METHODS FOR SCHRODINGER-POISSON EQUATIONS

As an important model in quantum semiconductor devices, the SchrSdinger-Poisson equations have generated widespread interests in both analysis and numerical simulations in recent years. In this paper, we present Gaussian beam methods for the numerical simulation of the one-dimensional Schrodinger-Poisson equations. The Gaussian beam methods for high frequency waves outperform the geometrical optics method in that the former are accurate even around caustics. The purposes of the paper are first to develop the Gaussian beam methods, based on our previous methods for the linear SchrSdinger equation, for the Schrodinger-Poisson equations, and then check their validity for this weakly-nonlinear system.     

FINITE ELEMENT METHODS FOR THE NAVIER-STOKES EQUATIONS BY H（div） ELEMENTS

We derived and analyzed a new numerical scheme for the Navier-Stokes equations by using H（div） conforming finite elements. A great deal of effort was given to an establishment of some Sobolev-type inequalities for piecewise smooth functions. In particular, the newly derived Sobolev inequalities were employed to provide a mathematical theory for the H（div） finite element scheme. For example, it was proved that the new finite element scheme has solutions which admit a certain boundedness in terms of the input data. A solution uniqueness was also possible when the input data satisfies a certain smallness condition. Optimal-order error estimates for the corresponding finite element solutions were established in various Sobolev norms. The finite element solutions from the new scheme feature a full satisfaction of the continuity equation which is highly demanded in scientific computing.     

A NUMERICAL STUDY FOR THE PERFORMANCE OF THE WENO SCHEMES BASED ON DIFFERENT NUMERICAL FLUXES FOR THE SHALLOW WATER EQUATIONS

In this paper we investigate the performance of the weighted essential non-oscillatory （WENO） methods based on different numerical fluxes, with the objective of obtaining better performance for the shallow water equations by choosing suitable numerical fluxes. We consider six numerical fluxes, i.e., Lax-Friedrichs, local Lax-Friedrichs, Engquist-Osher, Harten-Lax-van Leer, HLLC and the first-order centered fluxes, with the WENO finite volume method and TVD Runge-Kutta time discretization for the shallow water equations. The detailed numerical study is performed for both one-dimensional and two-dimensional shallow water equations by addressing the property, and resolution of discontinuities. issues of CPU cost, accuracy, non-oscillatory     

Regularity condition of solutions to the quasi-geostrophic equations in Besov spaces with negative indices

With a Hǒlder type inequality in Besov spaces, we show that every strong solution θ（t, x） on （0, T ） of the dissipative quasi-geostrophic equations can be continued beyond T provided that ⊥θ（t, x） ∈L 2γ/γ-2δ （（0, T ）; B^δ-γ/2 ∞∞（R^2）） for 0 〈 δ 〈 γ/2 .     

Standing Waves of the Coupled Nonlinear Schr ¨odinger Equations

In this paper, we study the existence of standIng waves of the coupled nonlinear Schrodinger equations. The proofs of which rely on the Lyapunov-Schmidt methods and contraction mapping principle are due to F. Weinstein in .     

ON SOLVABILITY AND WAVEFORM RELAXATION METHODS FOR LINEAR VARIABLE-COEFFICIENT DIFFERENTIAL-ALGEBRAIC EQUATIONS

This paper is concerned with the solvability and waveform relaxation methods of linear variable-coefficient differential-algebraic equations （DAEs）. Most of the previous works have been focused on linear variable-coefficient DAEs with smooth coefficients and data, yet no results related to the convergence rate of the corresponding waveform relaxation methods has been obtained. In this paper, we develope the solvability theory for the linear variable-coefficient DAEs on Legesgue square-integrable function space in both traditional and least squares senses, and determine the convergence rate of the waveform relaxation methods for solving linear variable-coefficient DAEs.     

A MIXED FINITE ELEMENT METHOD ON A STAGGERED MESH FOR NAVIER-STOKES EQUATIONS

In this paper, we introduce a mixed finite element method on a staggered mesh for the numerical solution of the steady state Navier-Stokes equations in which the two components of the velocity and the pressure are defined on three different meshes. This method is a conforming quadrilateral Q1 × Q1 - P0 element approximation for the Navier-Stokes equations. First-order error estimates are obtained for both the velocity and the pressure. Numerical examples are presented to illustrate the effectiveness of the proposed method.     

偶数阶中立型偏微分方程系统的振动性定理

讨论了一类偶数阶中立型偏泛函微分方程系统在Robin边界条件下解的振动性,获得了所有解振动的若干充分条件,同时也给出了实际应用例子.     

ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS

The main purpose of this work is to provide a novel numerical approach for the Volterra integral equations based on a spectral approach. A Legendre-collocation method is proposed to solve the Volterra integral equations of the second kind. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical results confirm the theoretical prediction of the exponential rate of convergence. The result in this work seems to be the first successful spectral approach （with theoretical justification） for the Volterra type equations.     

Hilbert-Space-Valued Super-Brownian Motion and Related Evolution Equations

A stochastic partial differential equation in which the square root of the solution appears as the diffusion coefficient is studied as a particular case of stochastic evolution equations. Weak existence of a solution is proved by the Euler approximation scheme. The super-Brownian motion on [0, 1] is also studied as a Hilbert-space-valued equation. In this set up, weak existence, pathwise uniqueness, and positivity of solutions are obtained in any dimension d . Accepted 23 October 1998     

THE FINITE DIFFERENCE METHOD FOR DISSIPATIVE KLEIN-GORDON-SCHRODINGER EQUATIONS IN THREE SPACE DIMENSIONS

A fully discrete finite difference scheme for dissipative Klein-Gordon-SchrSdinger equations in three space dimensions is analyzed. On the basis of a series of the time-uniform priori estimates of the difference solutions and discrete version of Sobolev embedding the- orems, the stability of the difference scheme and the error bounds of optimal order for the difference solutions are obtained in H2 × H2 ×H1 over a finite time interval. Moreover, the existence of a maximal attractor is proved for a discrete dynamical system associated with the fully discrete finite difference scheme.     

具有连续时滞的中立型双曲偏微分方程系统解的振动性

研究了一类具有连续时滞变量的非线性中立型双曲型偏微分方程系统,解的振动性.获得了该方程组在Robin边值条件和Dirichlet边值条件下解振动的充分条件.