论文二维不可压缩层流模型数值解和稳定性是关于对写作层流论文范文与课题研究的大学硕士、相关本科毕业论文紊流与层流的区别论文开题报告范文和相关文献综述及职称论文参考文献资料下载有帮助。
摘 要层流流动稳定性的研究主要在于数学模型的建立以及求数值解.以Poiseuille流动为例,运用谱方法对二维不可压缩层流模型OrrSommerfeld方程进行了展开和数值计算,得到了相应的层流稳定性数值条件．计算结果表明,谱方法具备较高的数值精度和较少的计算时间．
关键词OrrSommerfeld方程；稳定性；数值计算；谱方法
中图分类号O357.1 文献标识码A 文章编号10002537（2012）05002004
层流流体稳定性问题的研究长期以来一直倍受数学家的关注．当雷诺数在某个特定范围内逐步增大时,流体将经历一个由层流向紊流转换的特殊过程［1］,即由一个稳定的、亚临界的层流流体转化成为一个超临界的紊流流体．为此,通过对流体雷诺数和波长的计算来预测层流的稳定性问题变得尤为重要．目前,人们通过对层流流体加入微小外部扰动来预测层流流体稳定性的方法已被广泛地接受［2］．如果在此过程中,外部扰动能够逐渐减弱直至消失,那么流体将继续保持稳定．反之,如果这些外部的扰动随着时间的推移逐渐增大,那么流体将经历一个由层流向紊流流动转化的过渡过程．
层流流体的稳定性问题可归结为OrrSommerfeld方程求解问题,该方程的求解在数学界引起了广泛的兴趣．Davey［3］通过有限元的方法来求解OrrSommerfeld方程,取得了良好的结果．然而由于有限元方法本身的局限性,运算精度会受到很大的限制．如在计算误差（Δx）p渐进减小的过程中,运用有限元方法只能得到一些有限的p值,同时计算时间过长,不适合广泛运用于求解此类问题．切比雪夫多项式在逼近理论中有重要的应用．第一类切比雪夫多项式的根（被称为切比雪夫节点）可以用于多项式插值．相应的插值多项式能最大限度地降低龙格现象,并且提供多项式在连续函数的最佳一致逼近．Orszag［4］运用切比雪夫多项式以及QR矩阵特征值算法所计算的结果有着较高的计算精度以及较少的计算时间．
当前,谱方法已成为计算流体力学领域的三大主流数值方法之一,并受到广大学者的关注．它是上世纪70年 展起来的一种数值求解偏微分方程的方法,该方法在理论上具有“无穷阶”收敛性,可采用快速算法,现已被广泛用于气象、物理、力学等诸多领域,成为继差分法和有限元法之后又一种重要的数值方法［57］．由于谱方法的“无穷阶”收敛性且具备高的计算精度,相对于有限元方法和切比雪夫多项式的方法,它能够对误差形成指数型收敛,提高了计算精度并缩短了计算时间,为求解OrrSommerfeld方程提供了新的思路.
本文运用谱方法的数值分析,以Poiseuille流动为例对OrrSommerfeld方程进行了展开和数值计算,并利用计算结果对层流流体的稳定性进行了分析讨论,并验证了谱方法的高精度性和快速收敛性．
图1Poiseuille流动稳定性计算图求解方程（12）的特征值c,需要给定α和Re值．如果ci<0,那么层流流动将保持稳定,相反地,如果ci>0,那么层流流体就会向紊流进行转化,ci等于0表示临界稳定状态．图1表示了ci随Re和α值变化的曲线（具体计算程序见附录1）．
图中分别显示了ci等于0.0,ci等于0.04,ci等于006的等值线．从图中可以看出由等值线ci等于00所包围的中间区域为不稳定区,也就是说当Re和α的值在这个范围内时,层流流体是不稳定的,极易向紊流进行转化．等值线ci等于0.0 被视为稳定区域．同时可以得到临界的雷诺数值为Recrit等于5773,也就是说如果当流体的雷诺数小于5773时,无论α取值多少,流体始终保持层流流动．同时在临界雷诺数处,首先产生不稳定流体的αc等于1.02．对于本文得到的临界数值可以和其他文献所得到的值进行比对,Thomas［9］运用有限元的方法得到Recrit≈5780,αc等于1.026．Orszag［4］通过切比雪夫多项式的方法得到Recrit≈5772,αc等于1.021．由此验证了通过谱方法求解OrrSommerfeldequation有着很高的计算精度以及较少的计算时间,为解决此类方程提供了新的思路．
参考文献：
［1］SCHOBEIRIMT.Fluidmechanicsforengineers［M］.Berlin：SpringerVerlag,2010.
［2］D EYA,DRAZINPG.Thestabilityofpoiseuilleflowinthepipe［J］.JFluidMech,1969,36（02）：209218.
［3］D EYA,NGUYENHPF.Finiteamplitudestabilityofpipeflow［J］.JFluidMech,1971,45（04）：701720.
［4］ORSZAGSA.Accuratesolutionoftheorrsommerfeldstabilityequation［J］.JFluidMech,1971,50（04）：689703.
［5］尚新民.谱方法的数值分析［M］.北京：科学出版社,2000.
［6］CANUTOC,HUSSAINIMY,QUARTERONIA,etal.Spectralmethodsinfluiddynamics［M］.Berlin：SpringerVerlag,1993.
［7］BALTENSPERGERR,TRUMMERMR.Spectraldifferencingwithatwist［J］.SiamJSciComput,2002,24（5）：14651487.
［8］DONW,SOLOMONOFFS.Accuracyandspeedincomputingthechebyshevcollocationderivative［J］.SiamJSciComput,1995,16（6）：12531268.
［9］THOMASLH.Thestabilityofplanepoiseuilleflow［J］.PhysRev,1953,91（4）：780783.
总结:这篇层流论文范文为免费优秀学术论文范文,可用于相关写作参考。
参考文献：
1、 有界域上一维双极量子力学模型解的经典极限 Abstract：We consider a onedimensional bipolar isentropic quantum hydrody。
2、 一种电能计量设备运维作业动态优化模型算法 摘 要：针对电网企业传统模式下运维作业所存在的任务分配不科学、路径规划不合理和无法实时响应动态需求等问题，提出了一种新的动态优化模型及求解算法 。
3、 动力学神经网络模型的构建与其稳定性 摘 要： 针对静态人工神经网络具有在反映系统动态行为时网络结构复杂、不能很好地反映系统动态性能的缺点，提出一种由带有积分器和可调反馈系数的神经元。
4、 求解矩阵方程XAHX—2AQ正定解数值方法 摘 要：非线性矩阵方程，其中，是正定矩阵，在控制理论、动态规划、随机筛选和统计、高斯互反过程、纳米研究以及铁路交通的动力分析中等领域有着广泛的应。
5、 基于AHP网络广告三维决策模型 网络广告是利用数字技术制作和表示的，在互联网上传播、发布的广告。自从1997年3月中国出现第一条商业性网络广告以来，网络广告得到了飞速发展，据艾。
6、 互联网多维价值网络模型 【摘 要】 以互联网产业化为代表的新技术革命给传统经济带来严峻考验，许多企业纷纷创新以求生存及发展。美的集团从2011年起两次自主战略转型，创新。