流体热动力学方程组的混合有限元解法

来源 :西安交通大学 | 被引量 : 0次 | 上传用户:reza777
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
由于流体流动用非线性微分方程(比如Navier-Stokes方程,流体热动力学方程组)描述,并且它的非线性来自动力学而不是物理假设,因此给研究它带来了巨大的困难.在流体流动的过程中如果伴有热传导就要用流体热动力学方程组来描述,它是Navier-Stokes方程和能量方程耦合,由Boussinesq逼近得到的方程组.事实上,在流体运动的过程中,由于流体有黏性必然产生热量.在流体流动的过程中就伴随着速度、压力和温度相互影响与转化.所以,研究流体热动力学方程组就有重要的理论和实践意义,也更能反映流体流动的普遍规律.本文主要研究流体热动力学方程组的混合有限元方法,其研究内容如下:  一、本文给出了求解定常流体热动力学方程组的Newton迭代混合有限元算法,构造了三种Newton迭代解法和一种两重网格算法.第一种Newton迭代混合有限元算法要求给出速度和温度初值,由动量方程和守恒方程求出速度和压力,然后由能量方程求出温度,一直循环直到满足迭代终止条件为止.在第二种Newton迭代混合有限元算法要求给出速度初值,由能量方程求出温度,然后由动量方程和守恒方程求出速度和压力,一直循环直到满足迭代终止条件为止.这两种算法在求解的过程中,把速度、压力和温度进行解耦求解,大大的减小了计算规模,节省计算时间.若把能量方程同样利用Newton迭代格式,就得到了耦合Newton迭代混合有限元格式.该算法需要同时给出速度和温度初值,把动量方程、守恒方程和能量方程耦合求解.算法虽然需要计算时间较多,但是能求解较小的粘性系数,计算结果最好.在第二种格式的基础上,本文给出了求解流体热动力学方程组的一种两层网格算法.在粗网格上通过m次线性迭代求解一个非线性流体热动力学方程组,然后在细网格上求解一个线性化的方程,从而达到节省计算量得目的,提高计算效率.本文对以上算法进行了详尽的误差分析,分析了有限元离散误差和非线性迭代误差.数值算例求解了方腔内自然对流换热问题和环腔内自热对流换热问题,并把求解结果和经典文献结果比较,表明了算法的可靠性和理论分析的正确性.  二、本文给出了求解定常流体热动力学方程组的亏量校正(Defect Correction)混合有限元算法.亏量校正算法先求解原方程的粘性系数加上σh作为人工粘性系数得到的一个近似方程,然后在同一套网格上经过亏量步求解一个线性问题进行校正,能很好的求解小粘性系数问题.本文给出了它的稳定性分析和误差估计以及数值算例.数值试验表明亏量校正有限元方法能求解更小的粘性系数,同时也证明了本文的理论分析结果正确.  三、本文给出了求解非定常流体热动力学方程组亏量校正混合有限元方法.本文首先考虑了非定常流体热动力学方程组的半离散亏量校正混合有限元方法,给出它的误差估计,得到了许多的理论结果.由于半离散方法是全离散基础,本文也研究了时间导数项用Euler格式离散的全离散亏量校正混合有限元方法,给出了理论分析和数值试验.理论分析表明该算法是稳定的并且具有良好的数值精度.数值试验证明了算法的优越性和理论分析的正确性.
其他文献
受贿700多万韩桂芝一审被判死缓 2005年12月15日,北京市第一中级人民法院依法对黑龙江省政协原主席韩桂芝受贿案作出一审判决,以受贿罪判处韩桂芝死刑,缓期两年执行。北京市
自从Daubechies构造了一维紧支集正交小波基以来,出现了各种各样的新的小波基.一维小波基被成功地应用于图像处理、数值计算和统计等领域.但在许多实际问题如图像压缩中用到
该论文将多元统计学中的判别分析方法和数学图象处理技术相结合,在模式识别原理基础上建立了分散型模式识别的模型并在计算机上的实现,成为实用软件.论文分别从胶料图像的图
大型辛矩阵特征值的计算在子结构链的动(静)力分析中十分重要,在离散时间最优控制大系统分析及金融数学中都有重要应用。在保证Hamilton结构始终不变的原则下来求解Hamilton矩
群G的子群H称为ss-拟正规的,如果存在K
在文献[1]中给出了许多求解多目标规划的方法,其中绝大多数方法是通过将多目标 规划问题转化成一个或多个非线性或线性规划问题,并借助于传统的求解非线性或线性规划的方法进
该文是基于硕研期间的研究工作写成,既概括地论述了小波变换的基本理论,又包含了自己的体会和总结.针对小波分析与实际应用的密切关系,该文从基本应用的角度研究了:1.基于多
近十年来,随机线性规划又引起了研究人员与实际工作者的重视.其中具补偿的随机性规划在不确定性可表为离散的Scenario集时,为决策者提供了对具有不确定因素或噪声数据的序列
简化Markovitz的证券组合选择模型的计算是证券组合选择理论和应用研究的一个重 要研究课题.指数模型算法是其中的一种简化算法,该文研究基于多指数模型原EGP算法, 同时研究
近年来,关于集值分析理论和微分包含理论已经逐步完善,可控性和优化控制方面的研究工作也得到进一步发展。本文研究了一类发展微分包含的可控性问题,利用不动点定理给出了发展微