Optimal error estimate of the penalty finite element method for the time-dependent Navier-Stokes equations

A fully discrete penalty finite element method is presented for the two-dimensional time-dependent Navier-Stokes equations. The time discretization of the penalty Navier-Stokes equations is based on the backward Euler scheme; the spatial discretization of the time discretized penalty Navier-Stokes equations is based on a finite element space pair which satisfies some approximate assumption. An optimal error estimate of the numerical velocity and pressure is provided for the fully discrete penalty finite element method when the parameters and are sufficiently small.

     

Noninterpolatory Hermite subdivision schemes

Bivariate interpolatory Hermite subdivision schemes have recently been applied to build free-form subdivision surfaces. It is well known to geometric modelling practitioners that interpolatory schemes typically lead to ``unfair" surfaces--surfaces with unwanted wiggles or undulations--and noninterpolatory (a.k.a. approximating in the CAGD community) schemes are much preferred in geometric modelling applications. In this article, we introduce, analyze and construct noninterpolatory Hermite subdivision schemes, a class of vector subdivision schemes which can be applied to iteratively refine Hermite data in a not necessarily interpolatory fashion. We also study symmetry properties of such subdivision schemes which are crucial for application in free-form subdivision surfaces.

A key step in our mathematical analysis of Hermite type subdivision schemes is that we make use of the strong convergence theory of refinement equations to convert a prescribed geometric condition on the subdivision scheme--namely, the subdivision scheme is of Hermite type--to an algebraic condition on the subdivision mask. The latter algebraic condition can then be used in a computational framework to construct specific schemes.

     

Robust Contour Tracking Model Using a Variational Level-Set Algorithm

In this work, we provide a novel variational level-set based object contour tracking approach. Thus, a mathematically rigorous variant of the Chan-Vese algorithm for image segmentation via geometric active contour model is proposed here. With respect to the original contour detection algorithm, the level set function ?(t) defining the evolving contour S _t  = {x; ?(t, x) = 0} is iteratively computed from a nonlinear parabolic boundary value problem that is well posed in the space of functions with bounded variations. We provide a robust mathematical justification of the proposed level-set model.     

三维有限差分波束传输法:用于Ti∶ LiNbO_3方向耦合器的模拟

提出求解三维傍轴近似波动方程的交替方向隐式差分格式，并用它模拟Ｔｉ∶ＬｉＮｂＯ３方向耦合器     

Numerical methods for a fixed domain formulation of the glacier profile problem with alternative boundary conditions

In this paper we develop a set of numerical techniques for the simulation of the profile evolution of a valley glacier in the framework of isothermal shallow ice approximation models. The different mathematical formulations are given in terms of a highly nonlinear parabolic equation. A first nonlinearity comes from the free boundary problem associated with the unknown basal extension of the glacier region. This feature is treated using a fixed domain complementarity formulation which is solved numerically by a duality method. The nonlinear diffusive term is explicitly treated in the time marching scheme. A convection dominated problem arises, so a characteristic scheme is proposed for the time discretization, while piecewise linear finite elements are used for the spatial discretization. The presence of infinite slopes in polar regimes motivates an alternative formulation based on a prescribed flux boundary condition at the head of the glacier instead a homogeneous Dirichlet one. Finally, several numerical examples illustrate the performance of the proposed methods.     

高精度高分辨率迎风格式应用于不同速度范围内粘性流动

提出了一种适合于不同速度范围的高精度高分辨率的迎风有限差分格式，并基于此数植模型发展了适应于速度范围极宽的非定常粘性流动通用软件，不仅适用于超音速下捕捉强间断面，跨音速及高亚音速下捕捉弱间断面和滑移面，还可以精确地模拟低速情况下的粘性流动。此软件可分别用于研究内流和外流的流动特性以及预估其粘性损失。     

A NEW NONCONFORMING MIXED FINITE ELEMENT SCHEME FOR THE STATIONARY NAVIER-STOKES EQUATIONS

In this article, a new stable nonconforming mixed finite element scheme is proposed for the stationary Navier-Stokes equations, in which a new low order CrouzeixRaviart type nonconforming rectangular element is taken for approximating space for the velocity and the piecewise constant element for the pressure. The optimal order error estimates for the approximation of both the velocity and the pressure in L2-norm are established, as well as one in broken H1-norm for the velocity. Numerical experiments are given which are consistent with our theoretical analysis.     

Analysis of an Implicit Finite Difference Scheme for Time Fractional Diffusion Equation

Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.     

一类强色散DGH方程的多辛普雷斯曼格式

DGH方程作为一类重要的非线性水波方程有着许多广泛的应用前景.基于Hamilton系统的多辛理论研究了一类强色散DGH方程的数值解法,利用多辛普雷斯曼方法构造了一种典型的半隐式的多辛格式.分析了该格式的局部能量和动量守恒律误差,并给出了数值算例.数值算例结果表明该多辛离散格式具有较好的长时间数值稳定性.