求解对流扩散方程和Navier-Stokes方程的积分方程法

发布时间：2018-03-13 10:15

  本文选题：对流扩散方程　切入点：Navier-Stokes方程　出处：《山东大学》2016年博士论文　论文类型：学位论文

【摘要】：对流扩散方程和Navier-Stokes方程是流体力学中的两类基本方程,它们在物理、化学和工程中有着广泛的应用。在简化的情形下通过求解这两类方程还能得到解析的结果,但当遇到复杂问题时就很难求得解析解,这时数值求解成为一种非常有效的手段,于是构造精确、稳定和高效的数值方法成为研究这两类问题的重要内容。近来,无网格方法,特别是基于积分方程理论提出的精度高和计算量少的算法正引起一些学者的兴趣。本文提出了一种积分方程法用于数值求解对流扩散方程和Navier-Stokes方程,通过和其它方法的比较,该方法展现出了一些优势。在研究稳态对流扩散问题时,首先引入一个关于格林函数的拉普拉斯方程,将格林函数展开成级数形式并代入拉普拉斯方程,可将拉普拉斯方程变为一代数方程组,求解该方程组可得到级数形式的格林函数。利用格林函数的性质,可将对流扩散方程转化为积分方程。此时,再用相同的正交多项式把所要求解的未知量写成级数形式,然后利用多项式的正交性质,可把积分方程化为一个代数方程组,求解该方程组即可得到对流扩散方程级数形式的近似解。最后,借助Chebyshev多项式和Fourier级数,应用积分方程法求解了非齐次边界条件的一维对流扩散问题和齐次边界条件的二维对流扩散问题。在与有限体积法、有限元法和迎风差分法的比较中积分方程法都表现出了很高的精度,尤其在处理对流占优的对流扩散问题时更是展示出了很好的稳定性。对于非稳态对流扩散方程,由于与稳态情形相比,方程中多了时间变量,因此,需要考虑时间变量的离散方式。这里选用Crank-Nicolson方法对方程的时间变量进行离散,这种方法既保证了离散格式的简洁,又使结果达到了较好的精度。对时间变量离散后得到的方程可以看做是稳态形式的方程,这时构造与此方程相关的关于格林函数的拉普拉斯方程,并利用格林函数的性质,就可以将离散后的对流扩散方程转化为关于空间变量的积分方程。通过把格林函数和未知变量展开成级数形式,可进一步将积分方程简化为一个代数方程组。求解该方程组就可以得到用级数形式的有限和表示的非稳态对流扩散问题的近似解。在数值实验部分,用两个一维非稳态对流扩散问题和四个二维非稳态对流扩散问题检验了该方法。在一维问题中,给出的左右边界条件一侧为第一类边界条件,另一侧为第二类边界条件。在二维问题中,给出的既有对流占优的问题,也有非常数对流速度的问题,并将计算结果与多种方法进行了比较。在与变分多尺度方法的比较中,两种方法都在网格数较小时就达到了很高的精度,但积分方程法的精度(尤其是计算效率)要明显好于变分多尺度方法。在与有限体积元法的比较中,当网格数较少时,积分方程法的优势并不明显,但随着网格数的增加和时间步长的减小,积分方程法误差减小的速度显著快于有限体积元法。在求解流体力学中的Navier-Stokes方程时,如何处理速度-压力的耦合是一棘手的问题。通过比较现有的各种处理该问题的方法,并考虑到积分方程法的特点,本文将采用投影法处理Navier-Stokes方程中的速度-压力耦合问题。为了将投影法离散后的方程转化为积分方程,文中根据离散方程的形式分别引入了格林函数满足的拉普拉斯方程。然后利用格林函数的性质,将离散后的Navier-Stokes方程转化为积分方程。由于所研究的Navier-Stokes方程中速度场满足的是一般性的边界条件,并且投影法中的中间变量满足的是Neumann边界条件,所以为了计算的方便,统一采用了Chebyshev多项式对格林函数和未知变量进行展开。应用Chebyshev多项式的性质可以将所要求解的方程离散化为代数方程组,求解这些方程组就可以得到相应方程的解。最后,本文给出了一个算例用于检验该方法。计算结果表明与分数步法相比,积分方程法具有很好的精确性和收敛性,并且对于相同节点数,所用的CPU时间要明显少于分数步法。文中还专门讨论了周期边界条件下的对流扩散方程和Navier-Stokes方程。在这种边界条件下,由于可以采用Fourier级数对格林函数和未知变量进行展开,所以计算格式相对简单。这里也是首先构造格林函数的拉普拉斯方程,再利用格林函数的性质将对流扩散方程或者Navier-Stokes方程转化为积分方程,并进一步利用级数的正交性,将积分方程简化为一个常微分方程组,最后应用TVDRunge-Kutta方法对该方程组进行数值求解。计算结果表明,在求解对流扩散问题时,积分方程法的精度以及收敛性要好于局部间断Galerkin方法。在求解具有不同Reynolds数的Navier-Stokes方程时,积分方程法和投影法相比也有很高的精度。尤其对于高Reynolds数的不可压缩流动问题,积分方程法的优势更加明显。
[Abstract]:......
【学位授予单位】：山东大学
【学位级别】：博士
【学位授予年份】：2016
【分类号】：O241.82

【相似文献】

相关期刊论文 前10条

1 樊明武;王敬国;;积分方程法计算静磁场[J];华北电力学院学报;1984年04期

2 刘兴业,郑建军;旋转变厚度圆盘的样条积分方程法[J];航空动力学报;1989年03期

3 刘宝平;多自由度系统非线性振动的积分方程法[J];数学学报;1978年01期

4 云天铨;线载荷积分方程法的位移和应力场的唯一性定理[J];应用数学和力学;1985年03期

5 王有成;圆板非线性分析的样条积分方程法[J];合肥工业大学学报;1985年02期

6 云天铨;线载荷积分方程法分析嵌在非均匀半空间的轴向受力桩[J];力学学报;1986年S1期

7 苏继超,吴礼义;积分方程法计算翼型的跨音速绕流[J];航空学报;1987年11期

8 黄烈德;航天数学中积分方程法的作用[J];中国空间科学技术;1987年03期

9 李志能 ,邹红武 ,沈庆垓;用积分方程法计算带极靴的磁透镜[J];浙江大学学报(自然科学版);1988年05期

10 季海波,王敏中,武际可;用积分方程法分析承受轴向载荷的桩[J];北京大学学报(自然科学版);1988年03期

相关会议论文 前5条

1 陈小凤;;电磁学数值计算方法简介[A];中国地球物理学会第22届年会论文集[C];2006年

2 陈桂波;汪宏年;姚敬金;;用积分方程法模拟各向异性地层三维电磁响应[A];中国地球物理学会第二十四届年会论文集[C];2008年

3 李帝铨;底青云;王妙月;;“地—电离层”模式有源电磁场三维积分方程法正演[A];中国地球物理学会第二十七届年会论文集[C];2011年

4 陆洋;;解测高一重力混合边值问题的椭球谐积分方程法[A];《大地测量与地球动力学进展》论文集[C];2004年

5 孙建国;;稳定电流场的电位反射函数[A];中国地球物理.2003——中国地球物理学会第十九届年会论文集[C];2003年

相关博士学位论文 前1条

1 魏涛;求解对流扩散方程和Navier-Stokes方程的积分方程法[D];山东大学;2016年

相关硕士学位论文 前7条

1 吴红利;混合边界裂缝散射问题的积分方程方法[D];东华理工大学;2014年

2 彭笠;基于积分方程法的电法测井响应数值计算研究[D];电子科技大学;2012年

3 王劲松;大地电磁测深积分方程法三维正演问题的研究[D];中国地质大学（北京）;2006年

4 彭飞;移动源激发下CSAMT三维正演研究[D];中国地质大学（北京）;2007年

5 郝爱芹;用积分方程法数值模拟球形导电区域的磁流体自激发电现象[D];山东大学;2014年

6 胡俊华;瞬变电磁积分方程法正演模拟研究[D];中国地质大学;2014年

7 秦三团;电磁散射的时域积分方程法研究[D];西安电子科技大学;2005年

，

本文编号：1605940