一维对流扩散方程非均匀网格上的指数型高阶紧致差分格式

基于非均匀网格上函数的泰勒级数展开,结合残参量修正法,推导了非均匀网格上对流扩散方程的高阶指数型紧致差分格式,选取的算例表明,格式兼有高精度和高分辨率的优点,能够很好的适用于大梯度变化,计算区域中含边界层和对流占优区域中的流动问题的求解.     

Composite Finite Elements for Elliptic Interface Problems

We present a Composite Finite Element Method for the approximation of linear elliptic boundary value problems of Dirichlet type with discontinuous coefficients. The challenge is the discontinuity of the coefficient (interface) which is not necessarily resolved by the underlying finite element mesh. The method is non-conforming in the sense that shape functions preserve continuity across the interface only in an approximative way. However, the construction allows to balance the non-conformity and the overall discretization error. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

用有限元方法求解界面特征值问题

用有限元方法计算椭圆型界面特征值问题,实验数据显示近似特征值的变化规律:界面特征值问题中系数的间断性对协调和非协调Crouzeix-Raviart有限元特征值的收敛性并无影响,而且对协调有限元特征值外推以后得到高精度的解,相应的外推值还提供特征值下界;Crouzeix-Raviart元特征值提供特征值下界,这对一般有界区域如"镂空"型区域也成立.另外,还展示近似特征函数的图形.     

The Automorphism Group of a Compact Generalized Polygon has Finite Dimension

 For every compact generalized polygon whose point space has finite dimension, we show that the automorphism group has finite dimension as well. This opens the field for the application of Lie theoretic methods. Received 12 November 1997; in revised form 10 March 1998     

Explicit Solutions of Nonconvex Variational Problems in Dimension One

We propose a method of finding the generalized solutions of nonconvex variational problems by solving an appropriate differential inclusion that is motivated by necessary conditions of optimality for such generalized minimizers. Accepted 28 September 1998     

Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

<正>In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.     

and Norms Error Estimates in Finite Element Method for Linear Parabolic Interface Problems

We derive new a priori error estimates for linear parabolic equations with discontinuous coefficients. Due to low global regularity of the solutions the error analysis of the standard finite element method for parabolic problems is difficult to adopt for parabolic interface problems. A finite element procedure is, therefore, proposed and analyzed in this paper. We are able to show that the standard energy technique of finite element method for non-interface parabolic problems can be extended to parabolic interface problems if we allow interface triangles to be curved triangles. Optimal pointwise-in-time error estimates in the L ²(Ω) and H ¹(Ω) norms are shown to hold for the semidiscrete scheme. A fully discrete scheme based on backward Euler method is analyzed and pointwise-in-time error estimates are derived. The interfaces are assumed to be arbitrary shape but smooth for our purpose.     

一种新型的高次多结点薄壁杆件单元

应用Hermite插值多项式以及最小势能原理导出了一种新型的高次多结点薄壁杆件单元,这种单元的特点是精度高,并可应用于桥梁、高层建筑巨型结构的有限元分析中.     

High-Order Multioperator Compact Schemes for Numerical Simulation of Unsteady Subsonic Airfoil Flow

On the basis of high-order schemes, the viscous gas flow over the NACA2212 airfoil is numerically simulated at a free-stream Mach number of 0.3 and Reynolds numbers ranging from 10³ to 10⁷. Flow regimes sequentially varying due to variations in the free-stream viscosity are considered. Vortex structures developing on the airfoil surface are investigated, and a physical interpretation of this phenomenon is given.     

On the Convergence of Finite Element Method for Second Order Elliptic Interface Problems

The purpose of this paper is to study the convergence of finite element approximation to the exact solution of general self-adjoint elliptic equations with discontinuous coefficients. Due to low global regularity of the solution, it is difficult to achieve optimal order of convergence with classical finite element methods [Numer. Math. 1998; 79:175–202]. In this paper, an isoparametric type of discretization is used to prove optimal order error estimates in L ² and H ¹ norms when the global regularity of the solution is low. The interface is assumed to be of arbitrary shape and is smooth for our purpose. Further, for the purpose of numerical computations, we discuss the effect of numerical quadrature on finite element solution, and the related optimal order estimates are also established.     

A Well-Conditioned,Nonconforming Nitsche's Extended Finite Element Method for Elliptic Interface Problems

In this paper, we introduce a nonconforming Nitsche's extended finite element method (NXFEM) for elliptic interface problems on unfitted triangulation elements. The solution on each side of the interface is separately expanded in the standard nonconforming piecewise linear polynomials with the edge averages as degrees of freedom. The jump conditions on the interface and the discontinuities on the cut edges (the segment of edges cut by the interface) are weakly enforced by the Nitsche's approach. In the method, the harmonic weighted fluxes are used and the extra stabilization terms on the interface edges and cut edges are added to guarantee the stability and the well conditioning. We prove that the convergence order of the errors in energy and $L^2$ norms are optimal. Moreover, the errors are independent of the position of the interface relative to the mesh and the ratio of the discontinuous coefficients. Furthermore, we prove that the condition number of the system matrix is independent of the interface position. Numerical examples are given to confirm the theoretical results.     

一维Markov过程的Nash不等式

本文得到了生灭过程和一维扩散过程满足Nash不等式的判别准则,并证明 了对此二类过程,非常返性蕴含相应半群如下收敛速度||P(t)||1→∞≤Ct-1.同时也给 出一般马氏链满足Nash不等式的充分条件.     

Nash Inequalities for Markov Processes in Dimension One

In this paper, we give characterizations of Nash inequalities for birth-death process and diffusion process on the line. As a by-product, we prove that for these processes, transience implies that the semigroups P(t) decay as ∥P(t)∥_1→∞≥Ct ₋₁. Sufficient conditions for general Markov chains are also obtained. Received November 29, 2000, Revised February 21, 2001, Accepted March 16, 2001     

High-Order Time Integration in the p-Version of Finite Elements

In isothermal solid mechanics, the space and time domain need to be discretized. This paper presents an approach using high-order finite elements based on Legendre polynomials for the spatial discretization and time-adaptive diagonal-implicit Runge-Kutta methods for the temporal discretization. The superiority of this approach over a traditional one using linear shape functions for the space domain and a Backward-Euler scheme for the time domain is demonstrated. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Inhomogeneous Voter Models in One Dimension

We establish conditions for survival and extinction of types of one-dimensional voter models, and show that increasing the flip rates at a finite number of sites typically does not affect survival, unless the flipping mechanism is altered. We provide an example of a modified voter model that does not survive but can be made to survive simply by altering the flip mechanism at one site. We also show that a rather general class of such models have clustering behavior.     

Almost Split Sequences in Dimension One

Let R be a complete discrete valuation ring with quotient field K, and let Λ be an R-order in a semisimple K-algebra. For an indecomposable Λ-lattice E, a sublattice Bi E satisfying Rad E ? Bi E is defined, and it is shown that the middle term H of an almost split sequence τE ? H ? E can be obtained by an amalgamation of E/Bi E with E′/τE for a suitable overlattice E′ of τE. The method is bound to dim R = 1.     

Some Stability Inequalities for Compact Finite Difference Schemes

For finite difference schemes of compact form on nonuniform grids approximating m-th order two-point boundary value problems stability inequalities are proved which use a norm analogous to the Spijker-norm in the case of multistep methods. The results are applied to a number of finite difference schemes for which they establish a higher order of convergence than naively expected.     

关于高阶精度WCNS格式的无粘通量分裂方法

高阶精度加权紧致非线性格式(WCNS)越来越广泛地应用于复杂流动数值模拟.WCNS可以与多种无粘通量分裂方法结合起来使用.但是,常见的通量分裂方法都是基于低阶格式发展起来的,目前还不清楚哪些通量分裂方法最适合WCNS,也不知道这些方法与高阶格式结合时将会产生什么效果.表面热流计算是高超声速流动数值模拟的难点之一,为了在热流计算时选择合适的通量,研究了多种通量分裂方法的耗散大小.每种通量都可以表示成中心部分与耗散部分之和.这些通量的中心部分相同且非常简单,但是耗散部分较为复杂,且不同的通量分裂方法可导致不同的耗散表达式.通过对通量耗散进行分析可以发现耗散大小与网格界面两侧的物理量跳跃近似线性正相关.数值计算表明高阶格式得到的网格界面左右两侧的物理量跳跃通常远比低阶格式小,因而带来的通量耗散小.通过3个典型算例考察了通量耗散对热流计算的影响,其中包括高超激波/边界层干扰算例.基于对van Leer通量、Steger-Warming通量、KFVS通量、Roe通量、AUSM类通量和HLL类通量的考察,给出了通量选择建议.