一维高精度离散GDQ方法

GDQ method is a kind of high order accurate numerical methods developed several years ago, which have been successfully used to simulate the solution of smooth engineering problems such as structure mechanics and incompressible fluid dynamics. In this paper, extending the traditional GDQ method, we develop a new kind of discontinuous GDQ methods to solve compressible flow problems of which solutions may be discontinuous. In order to capture the local features of fluid flows, firstly, the computational domain is divided into many small pieces of subdomains. Then, in each small subdomain, the GDQ method is implementedand some kinds of numerical flux limitation conditions will be required to keep the correct flow direction. At the boundary interface between subdomains, we also use some kind of flux conditions according to the flow direction. The numerical method obtained by the above steps has the advantages of high order accuracy and easy to treat boundary conditions. It can simulate perfectly nonlinear waves such as shock, rarefaction wave and contact discontinuity. Finally, the numerical experiments on one dimensional Burgers equation and Euler equations are given.The numerical results verify the validation of the method.     

CONVERGENCE ANALYSIS ON ITERATIVE METHODS FOR SEMIDEFINITE SYSTEMS

The convergence analysis on the general iterative methods for the symmetric and positive semidefinite problems is presented in this paper. First, formulated are refined necessary and sumcient conditions for the energy norm convergence for iterative methods. Some illustrative examples for the conditions are also provided. The sharp convergence rate identity for the Gauss-Seidel method for the semidefinite system is obtained relying only on the pure matrix manipulations which guides us to obtain the convergence rate identity for the general successive subspace correction methods. The convergence rate identity for the successive subspace correction methods is obtained under the new conditions that the local correction schemes possess the local energy norm convergence. A convergence rate estimate is then derived in terms of the exact subspace solvers and the parameters that appear in the conditions. The uniform convergence of multigrid method for a model problem is proved by the convergence rate identity. The work can be regradled as unified and simplified analysis on the convergence of iteration methods for semidefinite problems [8, 9].     

Exact Boundary Controllability and Exact Boundary Observability for a Coupled System of Quasilinear Wave Equations

Based on the theory of semi-global classical solutions to quasilinear hyperbolic systems, the authors apply a unified constructive method to establish the local exact boundary(null) controllability and the local boundary(weak) observability for a coupled system of 1-D quasilinear wave equations with various types of boundary conditions.     

On Local Nonreflecting Boundary Conditions for Time Dependent Wave Propagation

The simulation of wave phenomena in unbounded domains generally requires an artificial boundary to truncate the unbounded exterior and limit the computation to a finite region. At the artificial boundary a boundary condition is then needed, which allows the propagating waves to exit the computational domain without spurious reflection. In 1977, Engquist and Majda proposed the first hierarchy of absorbing boundary conditions, which allows a systematic reduction of spurious reflection without moving the artificial boundary farther away from the scatterer. Their pioneering work, which initiated an entire research area, is reviewed here from a modern perspective. Recent developments such as high-order local conditions and their extension to multiple scattering are also presented. Finally, the accuracy of high-order local conditions is demonstrated through numerical experiments.     

Natural Boundary Integral Method and its New Development

In this paper, the natural boundary integral method, and some related methods, including coupling method of the natural boundary elements and finite elements, which is also called DtN method or the method with exact artificial boundary conditions, domain decomposition methods based on the natural boundary reduction, and the adaptive boundary element method with hyper-singular a posteriori error estimates, are discussed.     

THE GENERALIZED LOCAL HERMITIAN AND SKEW-HERMITIAN SPLITTING ITERATION METHODS FOR THE NON-HERMITIAN GENERALIZED SADDLE POINT PROBLEMS

For large and sparse saddle point problems, Zhu studied a class of generalized local Hermitian and skew-Hermitian splitting iteration methods for non-Hermitian saddle point problem [M.-Z. Zhu, Appl. Math. Comput. 218 （2012） 8816-8824 ]. In this paper, we further investigate the generalized local Hermitian and skew-Hermitian splitting （GLHSS） iteration methods for solving non-Hermitian generalized saddle point problems. With different choices of the parameter matrices, we derive conditions for guaranteeing the con- vergence of these iterative methods. Numerical experiments are presented to illustrate the effectiveness of our GLHSS iteration methods as well as the preconditioners.     

A NONOVERLAPPING DOMAIN DECOMPOSITION METHOD FOR EXTERIOR 3-D PROBLEM

1. IntroductionIn recent y6ars, the elliptic boUndaly value problems ill unbounded domains have dlawnmore and more attention. TO solve an equation in an unbounded domain numerically, a basicidea is to licit the computation to a bounded domain by introducing an artWial boundary.Based on this idea, many numerical methods, such as the coupling of BEM and FEM, the FEMwith boundary conditions at atilicial boundary) the coupled finite-~ie elemellt ndhodthe DDM(domain decomposition method)(cf.,…     

Nonlinear semigroup approach to transport equations with delayed neutrons

This paper deal with a nonlinear transport equation with delayed neutron and general boundary conditions. We establish, via the nonlinear semigroups approach, the existence and uniqueness of the mild solution, weak solution, strong solution and local solution on L~p-spaces(1 ≤ p +∞). Local and non local evolution problems are discussed.     

A POSTERIORI ESTIMATOR OF NONCONFORMING FINITE ELEMENT METHOD FOR FOURTH ORDER ELLIPTIC PERTURBATION PROBLEMS

In this paper, we consider the nonconforming finite element approximations of fourth order elliptic perturbation problems in two dimensions. We present an a posteriori error estimator under certain conditions, and give an h-version adaptive algorithm based on the error estimation. The local behavior of the estimator is analyzed as well. This estimator works for several nonconforming methods, such as the modified Morley method and the modified Zienkiewicz method, and under some assumptions, it is an optimal one. Numerical examples are reported, with a linear stationary Cahn-HiUiard-type equation as a model problem.     

Fold completeness of a system of root vectors of a system of unbounded polynomial operator pencils in Banach spaces. II. Application

This paper is a continuation of the author’s paper in 2009,where the abstract theory of fold completeness in Banach spaces has been presented.Using obtained there abstract results,we consider now very general boundary value problems for ODEs and PDEs which polynomially depend on the spectral parameter in both the equation and the boundary conditions.Moreover,equations and boundary conditions may contain abstract operators as well.So,we deal,generally,with integro-differential equations,functional-differential equations,nonlocal boundary conditions,multipoint boundary conditions,integro-differential boundary conditions.We prove n-fold completeness of a system of root functions of considered problems in the corresponding direct sum of Sobolev spaces in the Banach Lq-framework,in contrast to previously known results in the Hilbert L 2-framework.Some concrete mechanical problems are also presented.     

THE NAVIER-STOKES EQUATIONS WITH THE KINEMATIC AND VORTICITY BOUNDARY CONDITIONS ON NON-FLAT BOUNDARIES

We study the initial-boundary value problem of the Navier-Stokes equations for incompressible fluids in a general domain in R^n with compact and smooth boundary, subject to the kinematic and vorticity boundary conditions on the non-flat boundary. We observe that, under the nonhomogeneous boundary conditions, the pressure p can be still recovered by solving the Neumann problem for the Poisson equation. Then we establish the well-posedness of the unsteady Stokes equations and employ the solution to reduce our initial-boundary value problem into an initial-boundary value problem with absolute boundary conditions. Based on this, we first establish the well-posedness for an appropriate local linearized problem with the absolute boundary conditions and the initial condition （without the incompressibility condition）, which establishes a velocity mapping. Then we develop apriori estimates for the velocity mapping, especially involving the Sobolev norm for the time-derivative of the mapping to deal with the complicated boundary conditions, which leads to the existence of the fixed point of the mapping and the existence of solutions to our initial-boundary value problem. Finally, we establish that, when the viscosity coefficient tends zero, the strong solutions of the initial-boundary value problem in R^n（n ≥ 3） with nonhomogeneous vorticity boundary condition converge in L^2 to the corresponding Euler equations satisfying the kinematic condition.     

A PRECONDITIONER FOR COUPLING SYSTEM OF NATURAL BOUNDARY ELEMENT AND COMPOSITE GRID FINITE ELEMENT

1. IntroductionThe coupling of boundary elemellts and finite elements is of great imPortance for the nu-mercal treatment of boundary value problems posed on unbounded domains. It permits us tocombine the advanages of boundary elements for treating domains extended to infinity withthose of finite elemenis in treating the comP1icated bounded domains.The standard procedure of coupling the boundary elemeni and finite elemeni methods isdescribed as follows. First, the (unbounded) domain is divided…     

Global Bounded Solutions for Reaction-diffusion Systems with a Full Matrix of Diffusion Coefficients and a Balance Law

This paper deals with an initial boundary value problem for the strongly coupled reaction-diffusion systems with a full matrix of diffusion coefficients. The global existence of solutions is proved by using the techniques based on invariant regions, Lyapunov functional methods, and local Lp prior estimates independent of time.     

NUMERICAL BOUNDARY CONDITIONS FOR THE FAST SWEEPING HIGH ORDER WENO METHODS FOR SOLVING THE EIKONAL EQUATION

High order fast sweeping methods have been developed recently in the literature to solve static Hamilton-Jacobi equations efficiently. Comparing with the first order fast sweeping methods, the high order fast sweeping methods are more accurate, but they often require additional numerical boundary treatment for several grid points near the boundary because of the wider numerical stencil. It is particularly important to treat the points near the inflow boundary accurately, as the information would flow into the computational domain and would affect global accuracy. In the literature, the numerical solution at these boundary points are either fixed with the exact solution, which is not always feasible, or computed with a first order discretization, which could reduce the global accuracy. In this paper, we discuss two strategies to handle the inflow boundary conditions. One is based on the numerical solutions of a first order fast sweeping method with several different mesh sizes near the boundary and a Richardson extrapolation, the other is based on a Lax-Wendroff type procedure to repeatedly utilizing the PDE to write the normal spatial derivatives to the inflow boundary in terms of the tangential derivatives, thereby obtaining high order solution values at the grid points near the inflow boundary. We explore these two approaches using the fast sweeping high order WENO scheme in [18] for solving the static Eikonal equation as a representative example. Numerical examples are given to demonstrate the performance of these two approaches.     

A PRACTICAL PARALLEL DIFFERENCE SCHEME FOR PARABOLIC EQUATIONS

A practical parallel difference scheme for parabolic equations is constructed as follows: to decompose the domain Ω into some overlapping subdomains, take flux of the last time layer as Neumann boundary conditions for the time layer on inner boundary points of subdomains, solve it with the fully implicit scheme on each subdomain, then take correspondent values of its neighbor subdomains as its values for inner boundary points of each subdomain and mean of its neighbor subdomain and itself at overlapping points. The scheme is unconditionally convergent. Though its truncation error is O(τ h), the convergent order for the solution can be improved to O(τ h2).     

关于非定常不可压Navier-Stokes方程的时间高精度隐式差分方法

The incompressible Navier-Stokes equations,upon spatial discretization,become a system of differential algebraic equations,formally of index2.But due to the special forms of the discrete gradient and disrete divergence,its index can be regarded as 1.Thus,in this paper,a systematic approach following the ODE theory and methods is presented for the construction of high-order time-accurate implicit schemes for the incompressible Navier-Stokes equations,with projection methods for efficiency of numerical solution.The 3rd order 3-step BDF with componentconsistent pressure-correction projection method is a first attempt in this direction;the related iterative solution of the auxiliary velocyty,the boundary conditions and the stability of the algorithm are discussed.Results of numerical tests on the incompressible Navier-Stokes equations with an exact solution are presented,confirming the accureacy,stability and component-consistency of the proposed method.     

EXISTENCE AND UNIQUENESS OF THE SOLUTION FOR NONLINEAR THIRD-ORDER ROBIN BOUNDARY VALUE PROBLEM WITH A TURNING POINT

The paper aims to obtain existence and uniqueness of the solution as well as asymptotic estimate of the solution for singularly perturbed nonlinear third- order Robin boundary value problem with a turning point.In order to achieve this aim,existence and uniqueness of the solution for third-order nonlinear Robin boundary value problem is derived first based on the upper and lower solutions method under relatively weaker conditions.In this manner,the goal of this paper is gained by applying the existence and uniqueness results mentioned above.     

SPECTRAL AND SPECTRAL ELEMENT METHODS FOR HIGH ORDER PROBLEMS WITH MIXED BOUNDARY CONDITIONS

In this paper, we investigate numerical methods for high order differential equations. We propose new spectral and spectral element methods for high order problems with mixed inhomogeneous boundary conditions, and prove their spectral accuracy by using the recent results on the Jacobi quasi-orthogonal approximation. Numerical results demonstrate the high accuracy of suggested algorithm, which also works well even for oscillating solutions.