对流扩散问题的交替方向特征有限元方法

1 引言 在对流扩散方程的数值方法研究中,近年来由Douglas等人提出一种特征线修正法,并广泛应用于油藏模拟问题,核废料污染问题,半导体器件瞬态问题等领域.采用这种方法处理对流为主扩散问题时可以在不降低计算精度的情况下提高计算效率.实际问题一般是多维的,无论是用有限元或是差分法进行数值解都要解高阶的代数方程组,计算是相当复杂的.因此,研究如何对多维问题进行降维处理的数值方法无疑有着很重要的理论和实际意义.由算子的近似分解理论导出的交替方向迭代法,即可以把多维问题化为一维问题迭代求解,具有存贮量少,计算效率高等优点.本文对矩形区域上的二维对流扩散方程     

非线性对流扩散问题的差分-流线扩散法

1．引言流线扩散法（简称SD方法）是由Huzhes和Brooks在1980年前后提出的一种数值求解对流占优扩散问题的新型有限元算法．随后，Johnson和N8vert将SD方法推广到发展型对流扩散问题（［1］，［2］，［3］）．熟知，对于对流扩散问题，标准有限元法虽具有高阶精度，但常产生数值振荡；古典人工粘性Galerkin法更具有较好的稳定性，但仅具有一阶精度．而（SD方法兼具良好的数值稳定性和高阶精度，因此得到了越来越多的重视，对于发展型对流扩散问题，传统的SD方法均采用时空有限元．这样做，虽然可使时间和空间方向上的精度很好的协调起…     

非线性对流扩散方程的逆风型差分格式

§1 对流扩散方程可以用来描述水中和大气中污染物质的分布、流体流动和流体中的传热等,以上均有对流扩散的特征。数值求解对流扩散方程很重要,近年来工作不少,如[1,2]。我们考虑二维非线性对流扩散方程的初边值问题     

对流扩散方程的局部坐标有限分析法

本文给出了一种在不规则四边形网格上求解二维对流扩散方程的局部坐标有限分析法.文中还证明了这样导出的非恒定对流扩散方程的解在一定条件下是L_∞稳定的.数值例子证实计算是有效的.     

非线性对流扩散问题的交替方向差分流线扩散法

1引言Peaceman,Douglas等人于1955年提出了差分格式的交替方向法。随后,Douglas,Dupont于1972年又提出了有限元格式的交替方向法[1]。其基本思想是:对两个或三个空间变量的二阶抛物型和双曲型问题,将交替方向法与Galerkin方法相结合,通过算子分裂技术,把高维问题转化为一系列低维问题,交替地沿各空间变量的方向求解。[2]、[3]和[4]给出了对更一般扩散问题(带对流项的抛物方程)的数值求解和误差分析。     

对流扩散方程带罚函数项的有限体积格式及其分析

1引言 有限体积方法[l]一l’]作为守恒型的离散技术，被广泛应用于工程计算领域.文【2，3} 基于分片常数和分片常向量函数空间，对二维驻定对流扩散方程提出了一类非协调混合 有限体积(Covolume)格式，证明了格式具有。(hl/2)收敛精度.但该格式要求对偶剖分 比较规则，即采用重     

对流扩散问题的Crank-Nicolson差分-流线扩散法

1　引　言Streamline- Diffusion method (SD方法 )是近年来 Hughes和 Brooks提出的一种求解定常的对流占优和对流扩散问题的人工粘性有限元方法[1 ] ,[2 ] ,它具有标准有限元方法的高阶精度特点和人工粘性 Galerkin方法的稳定性特点 ,因此越来越受到人们的重视 .现在 ,SD方法已被推广到 Euler方程和 Navier- Stokes方程等发展型对流扩散问题[3 ] [4] ,但是常常采用时空有限元 [3 ] [5] ,这样能把时间和空间的精度很好地统一起来 ,却增大了数值计算的复杂性 ,基于此 [6 ]对非线性的对流占优扩散问题提出一种 Finite Difference- Strea…     

对流扩散方程的四阶紧凑迎风差分格式

§1.引言 流动和传热传质的基本方程均是对流扩散型的.对流扩散方程的高阶紧凑差分格式,作为提高计算可靠性和节省计算量的一条有效途径,已引起相当的重视.作为该领域的一大进展,新近由Dennis推出的对流扩散方程四阶紧凑格式,在二维情形下呈九点式且勿须引入中间变量,只涉及对流扩散量本身,能在较粗网格下获取较为准确的数值结果.从本质上说,该格式系指数型四阶紧凑格式的多项式型翻版.它与指数型紧凑格     

二维非线性对流-扩散方程的特征-差分解法

1.引言 近年来,求解对流-扩散方程的修正特征线法日渐被人们重视,已有不少人讨论了这一方法,如[1]-[5]。其中[1]和[2]是该领域的奠基性工作。[6]讨论了对流-扩散方程的特征-差分方法,改善并推广了[1]中某些重要结果。本文将讨论二维非线性对流-扩     

非定常线性对流扩散问题的算子分裂半显式有限元分析

在现代科学及工程领域中,存在着许多同时伴有物质传输和动力耗散两种过程的物理系统.在数学上,它们常归结为对流占优的对流扩散方程或以这种方程占主导的方程组.这类方程具有殆双曲性质,其解函数常呈现局部大梯度变化,使得传统的求解抛物问题的数值方法常常     

Solution of a stochastic Darcy equation by polynomial chaos expansion

This paper deals with solving a boundary value problem for the Darcy equation with a random hydraulic conductivity field.We use an approach based on polynomial chaos expansion in a probability space of input data.We use a probabilistic collocation method to calculate the coefficients of the polynomial chaos expansion. The computational complexity of this algorithm is determined by the order of the polynomial chaos expansion and the number of terms in the Karhunen–Loève expansion. We calculate various Eulerian and Lagrangian statistical characteristics of the flow by the conventional Monte Carlo and probabilistic collocation methods. Our calculations show a significant advantage of the probabilistic collocation method over the directMonte Carlo algorithm.     

双比例延迟微分方程配置法的可达阶

近几十年来,随着延迟微分方程广泛应用于变最测试、信号传递、机械化工、生命科学及经济管理系统等实际中,人们对此类方程的数值计算要求越来越高．     

A hybrid collocation method for a nonlinear Volterra integral equation with weakly singular kernel

This work is concerned with the numerical solution of a nonlinear weakly singular Volterra integral equation. Owing to the singular behavior of the solution near the origin, the global convergence order of product integration and collocation methods is not optimal. In order to recover the optimal orders a hybrid collocation method is used which combines a non-polynomial approximation on the first subinterval followed by piecewise polynomial collocation on a graded mesh. Some numerical examples are presented which illustrate the theoretical results and the performance of the method. A comparison is made with the standard graded collocation method.     

Piecewise spectral collocation method for system of Volterra integral equations

The main purpose of this paper is to investigate the piecewise spectral collocation method for system of Volterra integral equations. The provided convergence analysis shows that the presented method performs better than global spectral collocation method and piecewise polynomial collocation method. Numerical experiments are carried out to confirm these theoretical results.     

Analysis of a moving collocation method for one-dimensional partial differential equations

A moving collocation method has been shown to be very efficient for the adaptive solution of second- and fourth-order time-dependent partial differential equations and forms the basis for the two robust codes MOVCOL and MOVCOL4.In this paper,the relations between the method and the traditional collocation and finite volume methods are investigated.It is shown that the moving collocation method inherits desirable properties of both methods: the ease of implementation and high-order convergence of the traditional collocation method and the mass conservation of the finite volume method.Convergence of the method in the maximum norm is proven for general linear two-point boundary value problems.Numerical results are given to demonstrate the convergence order of the method.     

多孔介质中两相不可压缩不易混溶渗流问题的特征配置法

多孔介质中两相不可压缩不易混溶渗流问题是非线性偏微分方程的耦合系统,其中压力方程是椭圆的用配置法逼近,而饱和度方程是对流占优的抛物方程,用特征配置法来逼近,并且证明了数值解的存在唯一性,最后得到了最优的误差估计.     

Indirect methods of collocation: Trefftz‐Herrera collocation

A nonstandard collocation method (TH‐collocation) is presented, where collocation is used to construct specialized weighting functions instead of the solution itself, as it is usual, so that in this sense it is an indirect method. TH‐collocation is shown to be as accurate as standard collocation, but computationally far more efficient. The present article is the first of a series devoted to explore thoroughly collocation methods. The following classification of collocation methods is introduced: direct‐nonoverlapping; indirect‐nonoverlapping; direct‐overlapping; and indirect‐overlapping. Most of the effort reported in the literature has gone to direct‐nonoverlapping methods. The procedure presented in this article falls into the indirect‐nonoverlapping category and it is based on Trefftz‐Herrera formulation. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 709–738, 1999     

Linear barycentric rational collocation method for solving heat conduction equation

The linear barycentric rational collocation method for solving heat conduction equation is presented. The matrix form of discrete heat conduction equation by collocation method is also obtained. With the help of convergence rate of the barycentric interpolation, the convergence rate of linear barycentric rational collocation method for solving heat conduction equation is proved. At last, several numerical examples are provided to validate the theoretical analysis.     

Numerical solution of calcium-mediated dendritic branch model

The standard algorithms for spatial discretizations of calcium-mediated dendritic branch models via finite difference methods are quite accurate, but they are also extremely slow. To improve computational efficiency we apply spatial discretization using a spectral collocation method. Simulations using the spectral collocation method are compared to the finite difference approach using a model for calcium-mediated restructuring with spine pruning. We find that the spectral collocation method is about fifteen times more efficient to achieve similar accuracy than the finite difference approach even though spectral collocation requires more steps.     

A collocation method for solving singular integro-differential equations

This paper describes a collocation method for numerically solving Cauchy-type linear singular integro-differential equations. The numerical method is based on the transformation of the integro-differential equation into an integral equation, and then applying a collocation method to solve the latter. The collocation points are chosen as the Chebyshev nodes. Uniform convergence of the resulting method is then discussed. Numerical examples are presented and solved by the numerical techniques.