Iterative operator-splitting methods with higher-order time integration methods and applications for parabolic partial differential equations

In this paper we design higher-order time integrators for systems of stiff ordinary differential equations. We combine implicit Runge–Kutta and BDF methods with iterative operator-splitting methods to obtain higher-order methods. The idea of decoupling each complicated operator in simpler operators with an adapted time scale allows to solve the problems more efficiently. We compare our new methods with the higher-order fractional-stepping Runge–Kutta methods, developed for stiff ordinary differential equations. The benefit is the individual handling of each operator with adapted standard higher-order time integrators. The methods are applied to equations for convection–diffusion reactions and we obtain higher-order results. Finally we discuss the applications of the iterative operator-splitting methods to multi-dimensional and multi-physical problems.     

EXISTENCE AND UNIQUENESS RESULTS FOR VISCOUS, HEAT-CONDUCTING 3-D FLUID WITH VACUUM

We prove the local existence and uniqueness of the strong solution to the 3-D full Navier-Stokes equations whose the viscosity coefficients and the thermal conductivity coefficient depend on the density and the temperature. The initial density may vanish in an open set and the domain could be bounded or unbounded. Finally, we show the blow-up of the smooth solution to the compressible Navier-Stokes equations in R^n （n 〉 1） when the initial density has compactly support and the initial total momentum is nonzero.     

解非线性约束方程的拉格朗日全局投影方法

基于最优化方法求解约束非线性方程组的一个突出困难是计算 得到的仅是该优化问题的稳定点或局部极小点,而非方程组的解点.由此引出的问题是如何从一个稳定点出发得到一个相对于方程组解更好的点. 该文采用投影型算法,推广了Nazareth-Qi$^{[8,9]}$ 求解无约束非线性方程组的拉格朗日全局算法(Lagrangian Global-LG)于约束方程上; 理论上证明了从优化问题的稳定点出发,投影LG方法可寻找到一个更好的点. 数值试验证明了LG方法的有效性.     

TREANOR算法有病态吗?

文在计算“HL-1”两台电机并车过程时,、用了Treanor方法,且指出算法有病态现象。我们的分析是Treanor算法不存在病态现象,实际计算也表明Treanor算法可正确地描述并车过程。 描述两台双Y30°电机并车同步过程的微分方程为     

A fast algorithm for finding small solutions of <Emphasis Type=”Italic”>F(X,Y)=G(X,Y)</Emphasis> over number fields

Summary In the definition of semi-open sets and other kinds of generalized open sets, the role of the starting topology is given to a generalized topology. Some properties of the classes of sets so obtained is examined.     

THE EXPONENTIAL STABILIZATION FOR A SEMILINEAR WAVE EQUATION WITH LOCALLY DISTRIBUTED FEEDBACK

This paper considers the exponential decay of the solution to a damped semilinear wave equation with variable coefficients in the principal part by Riemannian multiplier method. A differential geometric condition that ensures the exponential decay is obtained.     

Brownian Chen series and Atiyah-Singer theorem

The purpose of this work is to give a new and short proof of the Atiyah-Singer local index theorem for the Dirac operator on the spin bundle. This proof is obtained by using heat semigroups approximations based on the truncation of Brownian Chen series.     

Euler schemes and large deviations for stochastic Volterra equations with singular kernels

In this paper, an Euler type approximation is constructed for stochastic Volterra equation with singular kernels, which provides an algorithm for numerical calculation. Then, the large deviation estimates of small perturbation to equations of this type are obtained. We finally apply them to SDEs with the kernel of fractional Brownian motion with Hurst parameter H∈(0,1).     

Parameter estimation of ordinary differential equations

This paper addresses the development of a new algorithm forparameter estimation of ordinary differential equations. Here,we show that (1) the simultaneous approach combined with orthogonalcyclic reduction can be used to reduce the estimation problemto an optimization problem subject to a fixed number of equalityconstraints without the need for structural information to devisea stable embedding in the case of non-trivial dichotomy and(2) the Newton approximation of the Hessian information of theLagrangian function of the estimation problem should be usedin cases where hypothesized models are incorrect or only a limitedamount of sample data is available. A new algorithm is proposedwhich includes the use of the sequential quadratic programming(SQP) Gauss–Newton approximation but also encompassesthe SQP Newton approximation along with tests of when to usethis approximation. This composite approach relaxes the restrictionson the SQP Gauss–Newton approximation that the hypothesizedmodel should be correct and the sample data set large enough.This new algorithm has been tested on two standard problems.     

An efficient transient Navier–Stokes solver on compact nonuniform space grids

In this paper, we propose an implicit higher-order compact (HOC) finite difference scheme for solving the two-dimensional (2D) unsteady Navier–Stokes (N–S) equations on nonuniform space grids. This temporally second-order accurate scheme which requires no transformation from the physical to the computational plane is at least third-order accurate in space, which has been demonstrated with numerical experiments. It efficiently captures both transient and steady-state solutions of the N–S equations with Dirichlet as well as Neumann boundary conditions. The proposed scheme is likely to be very useful for the computation of transient viscous flows involving free and wall bounded shear layers which invariably contain spatial scale variation. Numerical results are presented and compared with analytical as well as established numerical data. Excellent comparison is obtained in all the cases.