对流扩散方程引言

  • 格式:doc
  • 大小:351.50 KB
  • 文档页数:6

下载文档原格式

  / 7
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

对流扩散方程的定解问题是指物质输运与分子扩散的物理过程和黏性流体流动的数学模型,它可以用来描述河流污染、大气污染、核污染中污染物质的分布,流体的流动和流体中热传导等众多物理现象。关于对流扩散方程的求解很也备受关注,因此寻找一种稳定实用的数值方法有着重要的理论与实际意义。

求解对流扩散方程的数值方法有多种,尤其是对流占优扩散方程,这些方法有迎风有限元法,有限体积法,特征有限体积法,特征有限差分法和特征有限元法,广义差分法,流线扩散法,以及这些方法与传统方法相结合的方法如迎风广义差分法,迎风有限体积法有限体积——有限元法等这些方法数值求解效果较好,及有效的避免了数值震荡,有减少了数值扩散,但是一般计算量偏大

近年,许多研究者进行了更加深入的研究,文献提出了对流扩散方程的特征混合元法,再次基础上,陈掌引入了特征间断混合元方法,还有一些学者将特征线和有限体积法相结合,提出了特征有限体积元方法(非线性和半线性),于此同时迎风有限元也得到较大的发展,胡建伟等人研究了对流扩散问题的Galerkin 部分迎风有限元方法和非线性对流扩散问题的迎风有限元,之后又有人对求解发展型对流扩散问题的迎风有限元法进行了理论分析

有限差分法和有限元是求解偏微分方程的常用数值方法,一般情况下考虑对流占优的扩散方程,当对流项其主导作用时,其解函数具有大梯度的过渡层和边界层,导致数值计算困难,采用一般的有限元或有限体积方法虽然具有形式上的高精度,不能解决数值震荡的问题,虽然我们不能简单的将对流占优扩散方程看做对流方程,但由于次方程中含有一阶不对称的导数,对流扩散方程仍会表现出“对流效应”,从而采用迎风格式逼近,尽量反应次迎风特点,此格式简单,克服了锋线前沿的数值震荡,计算结果稳定,之前的迎风格式只能达到一阶精度,我们采用高精度的广义迎风格式,此格式是守恒的,精度高,稳定性好,具有单调性,并且是特征线法的近似,有效的避免了锋线前沿的数值震荡。

有限体积是求解偏微分方程的新的离散技术,日益受到重视。有限体积与有限差分、有限元法最大的区别及优点在于有限体积将求解区域内的计算转化到控制体积边界上进行计算,而后二者均是直接(或间接)在域内计算,故有限体积有着明显的物理涵义,在很大程度上减少计算工作量又能满足计算精度要求,加快收敛速度。由于此方法讲散度的积分化为子域边界积分后子啊离散,数值解满足离散守恒,而且可以采用非结构网格,所以在计算物理特别是计算流体力学领域上有限体积有广阔的前景。

间断Galerkin(DG)方法是在1973年,Reed和Hill在求解种子迁移问题时,针对一阶双曲问题的物理特点提出的。之后C.Johnson,G.R.Richter等人对双曲问题的DG方法做了进一步的研究,并且得到了该机的误差分析结果,由于这种方法具有沿流线从“上游”到“下游”逐层逐单元计算的显示求解的特点,并且可以进行并行计算,所以被广泛应用于各类方程的求解。最近Douglas等人在{25}中处理二阶椭圆问题时,得到DG方法的有限元空间不需要满足任何连续性条件,因此空间构造简单,具有较好的局部性和并行性。DG发展的一个重要方面是对对流占优扩散方程的应用。G.R.Richter等在1992年提出利用DG方法求解定长对流扩散问题

近年DG方法有了新的发展,其中YeXiu提出间断体积元方法备受人们关注,2004年,她将有限体积法与DG相结合,提出了椭圆问题的间断有限体积法,此方法解除了逼近函数在跨越边界上连续的限制,之后更多的研究者应用到Stokes问题,抛物问题,双曲问题,并得到了较好的结果,该方法不但继承了有限体积的高精度计算简单及保持物理间局部守恒等优点,而且有限元空间无需满足任何连续性要求,空间构造简单,有较好的局部和并行性。

当对流扩散方程中的对流项占主导地位时,方程具有双曲方程的特点,这是由于对流扩散方程中的非对称的对流项所引起的迎风效应使对流扩散方程的数值求解更困难,用传统的中心差分法和标准的有限元求解会差生数值的震荡,从而使数值模拟失真,为了克服这一困难,早在20世纪50年代,就有人提出了迎风思想,由于使用迎风技巧可以有效的消除数值解不稳定性,因此吸引了众多学者的关注,从1977年,Tabata等人就针对对流扩散方程提出了三角形网格上的迎风格式{42,38},并进行了深入的研究,梁栋基于广义差分法,提出并分析了一类建立在三角网格上的广义迎风差分格式,袁益让2001年就多层渗流方程组合系统提出并分析了迎风分数步长差分方法,以上均是讨论的线性对流扩散问题,胡建伟等通过引入质量集中算子,构造并分析了一类基于三角网格的质量集中型的部分有限元方法处理线性和非线性对流扩散问

题,近年,陈传军和杨青等人研究了半导体器件问题的迎风有限体积元方法,进一步推动了迎风方法的发展。(张

针对对流扩散方程的特点,尽管迎风差分格式能够消除对流项非对称效应引起的震荡现象,但它却“过量”的反应了解得情况,导致数值解得扩散。考虑到迎风效应是由对流扩散方程中的不对称的对流项引起的[1,2],本文直接考虑守恒的对流扩散方程,在前人工作的基础上,本文利用间断有限体积法处理扩散项,迎风间断有限体积处理对流项。通过数值试验本文的精度及收敛性是比较好的,并且解没有出现数值震荡和扩散现象。迎风间断有限体积元方法的优点主要是计算简单,保持物理间局部守恒,有限元空间无需满足任何连续性条件,空间构造简单,高精度,可以进行并行计算,节省工作量,有效的消除数值弥散和非物理震荡现象,具有较高的稳定性等。

这里我们以环境中污染物的扩散为例介绍一下方程中各部分的物理意义:u表示污染物的浓度,b代表水流的速度,D代表污染物扩散的能力,飞为源汇项,为了避免边界处理遇到的困难,我们假设问题是-周期的。

对于对流扩散方程的对流部分沿着迎来流方向取边界的值进行离散,得到较小的时间截断误差,当采用较大时间步长时,锋线前沿可以避免数值弥散现象。扩散部分采用间断有限体积元方法离散,保持局部质量守恒,当采用较大时间步长时,整个逼近格式不能保持质量守恒,因此,我们提出了守恒迎风间断有限体积元方法,并得到了离散解得最优L2误差估计。