Bakhvalov-Shishkin网格上求解边界层问题的差分进化算法

该文在Bakhvalov-Shishkin网格上求解具有左边界层或右边界层的对流扩散方程,并采用差分进化算法对Bakhvalov-Shishkin网格中的参数进行优化,获得了该网格上具有最优精度的数值解.对三个算例进行了数值模拟,数值结果表明:采用差分进化算法求解具有较高的计算精度和收敛性,特别是边界层的数值解精度明显...     

奇异摄动问题内罚间断有限方法的最优阶一致收敛性分析

本文在Bakhvalov-Shishkin网格上分析了采用高次元的内罚间断有限元方法求解一维对流扩散型奇异摄动问题的最优阶一致收敛性. 取k(k≥1)次分片多项式和网格剖分单元数为N时,在能量范数度量下, Bakhvalov-Shishkin网格上可获得O(N^-k)的一致误差估计. 在数值算例部分对理论分析结果进行了验证.     

椭圆型奇异摄动问题差分格式的一致收敛性分析

本文研究了一类椭圆型奇异摄动问题.利用Bakhvalov-Shishkin网格上的差分方法,获得了数值解一致一阶收敛于真解的结果.     

奇异摄动问题最优阶一致收敛的间断有限元分析

采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在Bakhvalov-Shishkin网格上具有最优阶的一致收敛性,即在能量范数度量下其误差估计为O(N~(-1)),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.     

Convergence analysis of sectional methods for solving breakage population balance equations-I: the fixed pivot technique

In this work we study the convergence of the fixed pivot techniques (Kumar and Ramkrishna Chem. Eng. Sci. 51, 1311–1332, 1996) for breakage problems. In particular, the convergence is investigated on four different types of uniform and non-uniform meshes. It is shown that the fixed pivot technique is second order convergent on a uniform and non-uniform smooth meshes. Furthermore, it gives first order convergence on a locally uniform mesh. Finally the analysis shows that the method does not converge on a non-uniform random mesh. The mathematical results of convergence analysis are also validated numerically.     

热传导-对流问题的回溯二重水平法

该文给出定常的热传导-对流问题的有限元逼近的一种二重水平方法. 这种二重水平方法包括解一个小的非线性的粗网格系统、一个细网格上的线性Oseen问题和一个粗网格上的线性校正问题. 同时,给出了这种近似解的存在性和收敛性分析.     

Postprocessing and higher order convergence for the mixed finite element approximations of the eigenvalue problem

In this paper, we propose a method to improve the convergence rate of the lowest order Raviart-Thomas mixed finite element approximations for the second order elliptic eigenvalue problem. Here, we prove a supercloseness result for the eigenfunction approximations and use a type of finite element postprocessing operator to construct an auxiliary source problem. Then solving the auxiliary additional source problem on an augmented mixed finite element space constructed by refining the mesh or by using the same mesh but increasing the order of corresponding mixed finite element space, we can increase the convergence order of the eigenpair approximation. This postprocessing method costs less computation than solving the eigenvalue problem on the finer mesh directly. Some numerical results are used to confirm the theoretical analysis.     

四阶奇异摄动边值问题在自适应网格上的一致收敛分析

we study a difference scheme for the fourth-order singular pertur-bation differential equation on the Bakhvalov-Shishkin grid by Green‘‘s function.The method is shown to be uniformly convergent with respect to the perturbation parameter,of order N^-2 in the maxmum norm on Bakhvalov-Shishkin meshes.Numerical results support our theoretical results.     

A non-monotone line search multidimensional filter-SQP method for general nonlinear programming

In this paper, we propose a non-monotone line search multidimensional filter-SQP method for general nonlinear programming based on the Wächter–Biegler methods for nonlinear equality constrained programming. Under mild conditions, the global convergence of the new method is proved. Furthermore, with the non-monotone technique and second order correction step, it is shown that the proposed method does not suffer from the Maratos effect, so that fast local convergence to second order sufficient local solutions is achieved. Numerical results show that the new approach is efficient.     

An additive Schwarz method for a stationary convection-diffusion problem

In this paper, we consider an additive Schwarz method applied to a linear, second order, nonsymmetric, indefinite problem. We discuss the solution of linear system of algebraic equations that arise from the streamline method for the above problem. An alternative linear system, which has the same solution as the system obtained by the streamline method, is derived and the GMRES method is used to solve this system. We show that the rate of convergence does not depend on the mesh size, nor on the number of local problems if the coarse mesh is fine enough.     

Numerical analysis of a second order algorithm for simplified magnetohydrodynamic flows

In this paper, we construct a second order algorithm based on the spectral deferred correction method for the time-dependent magnetohydrodynamics flows at a low magnetic Reynolds number. We present a complete theoretical analysis to prove that this algorithm is unconditionally stable, consistent and second order accuracy. Finally, two numerical examples are given to illustrate the convergence and effectiveness of our algorithm.     

A higher order difference method for singularly perturbed parabolic partial differential equations

This paper studies a higher order numerical method for the singularly perturbed parabolic convection-diffusion problems where the diffusion term is multiplied by a small perturbation parameter. In general, the solutions of these type of problems have a boundary layer. Here, we generate a spatial adaptive mesh based on the equidistribution of a positive monitor function. Implicit Euler method is used to discretize the time variable and an upwind scheme is considered in space direction. A higher order convergent solution with respect to space and time is obtained using the postprocessing based extrapolation approach. It is observed that the convergence is independent of perturbation parameter. This technique enhances the order of accuracy from first order uniform convergence to second order uniform convergence in space as well as in time. Comparative study with the existed meshes show the highly effective behavior of the present method.     

奇异摄动问题在Bakhvalov-Shishkin 网格上的流线扩散有限元逼近

本文采用线性插值的流线扩散有限元在Bakhvalov-Shishkin网格上求解一维对流扩散型的奇异摄动问题. 在ε ≤ N^-1的前提下,可以得到,关于扰动参数ε 是一致收敛的. 在离散的SD范数下,其u-u_I的误差阶提高到N^-2,u-u_h的误差阶达到N^-2（lnN）^0.5. 最后,通过数值算例,验证了理论分析.     

Explicit jump immersed interface method for virtual material design of the effective elastic moduli of composite materials

Virtual material design is the microscopic variation of materials in the computer, followed by the numerical evaluation of the effect of this variation on the material’s macroscopic properties. The goal of this procedure is an in some sense improved material. Here, we give examples regarding the dependence of the effective elastic moduli of a composite material on the geometry of the shape of an inclusion. A new approach on how to solve such interface problems avoids mesh generation and gives second order accurate results even in the vicinity of the interface. The Explicit Jump Immersed Interface Method is a finite difference method for elliptic partial differential equations that works on an equidistant Cartesian grid in spite of non-grid aligned discontinuities in equation parameters and solution. Near discontinuities, the standard finite difference approximations are modified by adding correction terms that involve jumps in the function and its derivatives. This work derives the correction terms for two dimensional linear elasticity with piecewise constant coefficients, i.e. for composite materials. It demonstrates numerically convergence and approximation properties of the method.      

A reduced Hessian algorithm with line search filter method for nonlinear programming

This paper proposes a line search filter reduced Hessian method for nonlinear equality constrained optimization. The feature of the presented algorithm is that the reduced Hessian method is used to produce a search direction, a backtracking line search procedure to generate step size, some filtered rules to determine step acceptance, second order correction technique to reduce infeasibility and overcome the Maratos effects. It is shown that this algorithm does not suffer from the Maratos effects by using second order correction step, and under mild assumptions fast convergence to second order sufficient local solutions is achieved. The numerical experiment is reported to show the effectiveness of the proposed algorithm.     

Mesh‐independent convergence of the modified inexact Newton method for a second order non‐linear problem

In this paper, we consider an inexact Newton method applied to a second order non‐linear problem with higher order non‐linearities. We provide conditions under which the method has a mesh‐independent rate of convergence. To do this, we are required, first, to set up the problem on a scale of Hilbert spaces and second, to devise a special iterative technique which converges in a higher than first order Sobolev norm. We show that the linear (Jacobian) system solved in Newton's method can be replaced with one iterative step provided that the initial non‐linear iterate is accurate enough. The closeness criteria can be taken independent of the mesh size. Finally, the results of numerical experiments are given to support the theory. Published in 2005 by John Wiley & Sons, Ltd.     

Improving the mass conservation of the level set method in a finite element context

In this Note, a new algorithm is proposed for improving the mass conservation of the level set method in the finite element context. Two kinds of Lagrange multipliers are introduced, associated respectively to the redistancing and advection equations. The first one, is located at the vicinity of the interface, while the second one is associated to a correction that is global to the domain. The performances of the proposed method are tested on the Zalesak test case, and the convergence rate versus the element mesh size are founded to be improved.     

A stabilized multi-level method for non-singular finite volume solutions of the stationary 3D Navier–Stokes equations

This paper proposes and analyzes a stabilized multi-level finite volume method (FVM) for solving the stationary 3D Navier?CStokes equations by using the lowest equal-order finite element pair without relying on any solution uniqueness condition. This multi-level stabilized FVM consists of solving the nonlinear problem on the coarsest mesh and then performing one Newton correction step on each subsequent mesh, thus only solving a large linear system. An optimal convergence rate for the finite volume approximations of nonsingular solutions is first obtained with the same order as that for the usual finite element solution by using a relationship between the stabilized FVM and a stabilized finite element method. Then the multi-level finite volume approximate solution is shown to have a convergence rate of the same order as that of the stabilized finite volume solution of the stationary Navier?CStokes equations on a fine mesh with an appropriate choice of the mesh size: ${ h_{j} ~ h_{j-1}^{2}, j = 1,\ldots, J}$ . Finally, numerical results presented validate our theoretical findings.     

An iteration formula for the simultaneous determination of the zeros of a polynomial

The need for efficient algorithms for determining zeros of given polynomials has been stressed in many applications. In this paper we give a new cubic iteration method for determining simultaneously all the zeros of a polynomial (assumed distinct) starting with ‘reasonably close’ initial approximations (also assumed distinct).The polynomial is expressed as an expansion in terms of the starting and their correction terms.A formula which gives cubic convergence without involving second derivatives is derived by retaining terms up to second order of the expansion in the correction terms.Numerical evidence is given to illustrate the cubic convergence of the process.     

Convergence of a continuous Galerkin method with mesh modification for nonlinear wave equations

We consider space-time continuous Galerkin methods with mesh modification in time for semilinear second order hyperbolic equations. We show a priori estimates in the energy norm without mesh conditions. Under reasonable assumptions on the choice of the spatial mesh in each time step we show optimal order convergence rates. Estimates of the jump in the Riesz projection in two successive time steps are also derived.