基于多层增量未知元方法的一类三维对流扩散方程的研究

对于一类一般形式的三维对流扩散方程, 运用有限差分方法, 在增量未知元方法(IU)下, 可以得到一个IU型正定但非对称的线性方程组.其系数矩阵条件数要远远优于不用IU方法的情形^[1]. 考虑到IU方法的这一优点, 作者在文中将IU方法与几种经典的迭代方法相结合, 来求解上述系统. 作者从理论上对该系统的IU型系数矩阵条件数进行了估计, 并通过数值试验验证了这几种IU型迭代方法的有效性.     

Local projection stabilisation for higher order discretisations of convection-diffusion problems on Shishkin meshes

We consider a singularly perturbed convection-diffusion equation on the unit square where the solution of the problem exhibits exponential boundary layers. In order to stabilise the discretisation, two techniques are combined: Shishkin meshes are used and the local projection method is applied. For arbitrary r≥2, the standard Q _r -element is enriched by just six additional functions leading to an element which contains the P _r+1. In the local projection norm, the difference between the solution of the stabilised discrete problem and an interpolant of the exact solution is of order uniformly in ε. Furthermore, it is shown that the method converges uniformly in ε of order in the global energy norm.      

A Petrov－Galerkin Method Based on Symmetrization for Parabolic Convection－Diffusion Problems

ＡＰｅｔｒｏｖ－Ｇａｌｅｒｋｉｎ　ＭｅｔｈｏｄＢａｓｅｄｏｎＳｙｍｍｅｔｒｉｚａｔｉｏｎｆｏｒＰａｒａｂｏｌｉｃＣｏｎｖｅｃｔｉｏｎ－ＤｉｆｆｕｓｉｏｎＰｒｏｂｌｅｍｓＷｕＷｅｉ（吴微）（ＤｅｐａｒｔｍｅｎｔｏｆＭａｔｈｅｍａｔｉｃｓＪｉｌｉｎＵｎｉ...     

Asymptotic behavior of solution to nonlinear evolution equations with damping

In this paper, we consider the global existence and the asymptotic behavior of solutions to the Cauchy problem for the following nonlinear evolution equations with ellipticity and dissipative effects

(E)     

流动注射分析中试样分散的对流-扩散理论

用两种方法求解适合于单一平直管道且无化学反应的流动注射分析(FIA)体系的对流-扩散方程,得到了描述塞状进样在管道中分散后的浓度分布公式,从而得到F曲线。讨论了反应管和采样管的内径和长度、流速和扩散系数对峰形和分散度的影响;采用数值法解出半峰宽与扩散系数的关系,提出了两种测定扩散系数的方法,为FIA体系的设计和指示物质、操作条件的选择提供了理论依据。     

对流扩散方程的涡区分离解法

对流扩散方程是流体计算中一个基本方程,常用的数值方法导至解一个高阶的代数方程组,要求较大的存贮量和较长的计算时间。本文提出一种涡区分离解法,它利用对流扩散方程的迎风性质,把涡区从对流支配区分离出来,仅在各个涡区建立代数方程组并求解。而在对流支配区,则充分利用其抛物性,只需采用显式格式进行计算。由于在各涡区建立的这些方程组阶数和带宽都较小,因此要求存贮量较小,计算速度较快。对于雷诺数较大,涡区范围较小的问题,该方法特别有效。     

A numerical method for mass conservative coupling between fluid flow and solute transport

We present a new coupled discretization approach for species transport in an incompressible fluid. The Navier-Stokes equations for the flow are discretized by the divergence-free Scott-Vogelius element on barycentrically refined meshes guaranteeing LBB stability. The convection-diffusion equation for species transport is discretized by the Voronoi finite volume method. In accordance to the continuous setting, due to the exact integration of the normal component of the flow through the Voronoi surfaces, the species concentration fulfills discrete global and local maximum principles. Besides of the numerical scheme itself, we present important aspects of its implementation. Further, for the case of homogeneous Dirichlet boundary conditions, we give a convergence proof for the coupled scheme. We report results of the application of the scheme to the interpretation of limiting current measurements in an electrochemical flow cell with cylindrical shape.     

A characterisation of oscillations in the discrete two-dimensional convection-diffusion equation

It is well known that discrete solutions to the convection-diffusion equation contain nonphysical oscillations when boundary layers are present but not resolved by the discretisation. However, except for one-dimensional problems, there is little analysis of this phenomenon. In this paper, we present an analysis of the two-dimensional problem with constant flow aligned with the grid, based on a Fourier decomposition of the discrete solution. For Galerkin bilinear finite element discretisations, we derive closed form expressions for the Fourier coefficients, showing them to be weighted sums of certain functions which are oscillatory when the mesh Péclet number is large. The oscillatory functions are determined as solutions to a set of three-term recurrences, and the weights are determined by the boundary conditions. These expressions are then used to characterise the oscillations of the discrete solution in terms of the mesh Péclet number and boundary conditions of the problem.