收藏本站
《重庆大学》 2012年
收藏 | 手机打开
二维码
手机客户端打开本文

谱元方法求解不可压缩流体流动及流动线性稳定性分析

梅欢  
【摘要】:本文的主要内容是利用谱元法数值模拟流体力学中的相关问题。总体思路如下:首先开展对谱元法基本算法的研究,分别建立了直角坐标系和极坐标系下的谱元法,并通过一些具有解析解的数值算例验证方法的精度和程序的有效性;其次结合时间分裂法,求解不同坐标系下的非定常不可压缩流动,在求解精度上与其它数值方法进行了比较;最后对方腔顶盖驱动流进行了线性稳定性分析,方腔采用有限长模型,研究展向(z方向)长度变化对流动失稳的影响,通过能量分析,探索流动失稳的物理机制。 针对上述研究思路,本文具体开展了如下工作,并完成了相应的程序开发:(1)推导了直角坐标系下的Chebyshev和Legendre谱元方法,并提出了极坐标系下的求解Poisson-型方程Legendre谱元方法和Fourier-Legendre谱元方法。在极坐标系下的谱元方法中,在单元径向上,变量采用Legendre多项式展开,极点所在单元的径向采用Gauss-Radau积分点,其它单元的径向采用Gauss-Lobatto积分点;在单元周向上,变量分别采用Legendre多项式和Fourier多项式展开。最后求解了多个具有解析解的Dirichlet或Neumann边界条件下的Poisson-型方程(Helmholtz方程和Poisson方程),用于验证本文方法的精度和有效性。(2)结合谱元方法和时间分裂法,求解了一系列非定常不可压缩流体流动,其中时间分裂法离散了Navier-Stokes方程中的时间项,从而得到相应的Poisson-型方程,进而通过谱元法进行求解。在直角坐标系下,分别求解了方腔顶盖驱动流、自然对流和表面张力流等流动,计算结果与文献基准解或有限体积法数值结果进行了对比,从而验证方法的精度和可行性。在极坐标系下,利用Fourier-Legendre谱元方法,以圆盘驱动流作为算例,研究加速坩埚旋转技术(ACRT)对于流体浓度均匀化的影响;首先,采用数值同位素模型,分别研究加速坩埚旋转和匀速坩埚旋转对于同位素流体浓度均匀化的作用,其中同位素流体的浓度变化仅由对流决定,以浓度的标准差判断均匀化程度,最后比较了谱元法和有限差分法的数值精度,对比了两种方法下的数值扩散的大小;其次,研究高温溶液晶体生长中加速坩埚旋转技术对于流体浓度均匀化的作用,流体浓度的变化由对流和扩散决定,本文选择了六种典型的加速坩埚旋转模式,并对每种旋转模式分别施加了不同的旋转时间周期,通过对浓度的标准差曲线变化找出最优的加速坩埚旋转模式和旋转时间周期,并以此指导实验。(3)结合谱元法和线性稳定性分析理论,分别研究方腔(即x-y平面)施加无滑移边界条件和滑移边界条件的有限长模型的顶盖驱动流失稳特性;其基本思想是在二维基态解上施加三维小扰动,去掉扰动的高阶非线性部分,将小扰动写成正则模形式,从而得到扰动的控制方程组,通过谱元离散方程组,并将方程组表达成积分弱形式,最后可以得到相应的广义特征值问题,利用Arpack程序包求解此特征值问题,对于不同的Re数(Reynold数)和波数,根据特征值的实部大小来寻找流动失稳的临界值;对临界流动进行扰动能量分析,探索不同展向长度下的流动失稳物理机制。 通过对不同问题的计算和相应数值结果进行分析,结果表明: ①本文的谱元方法能够以较少节点获得高精度数值解;在极坐标系下的谱元方法中,极点所在单元的径向采用Gauss-Radau积分点,能够成功地避免r=0处的1/r坐标奇异性;另外,本文谱元法可以通过区域分解技术避免极点附近节点的聚集,从而使得其在求解利用显示时间格式离散的非定常问题时,减轻对时间步长的限制。 ②利用谱元法和时间分裂法,成功地求解直角坐标系下的不可压缩流体流动;首先通过求解具有解析解的二维非定常Burgers方程,验证了本文方法的可行性和高精度性;其次将方腔顶盖驱动流、自然对流和表面张力流的数值结果同文献基准解或有限体积法数值结果进行对比,发现各种结果之间吻合得都非常好,从而说明本文方法能够正确地用于流体流动的求解,并为后续谱元法在极坐标系下的流动求解和流动的稳定性分析打下基础。 ③基于Fourier-Legendre谱元方法和有限差分法的数值结果比较表明:一阶迎风差分格式存在严重的数值假扩散现象,数值误差很大,增加节点数能稍微地改善结果;二阶迎风差分格式的数值扩散较严重,存在较大的误差,增加节点数时可以有效地改善数值结果;谱元方法求解的数值结果显示,谱元法存在极小的数值扩散,在匀速旋转时数值解几乎与理论解一致,在加速旋转时,标准差做非常均匀的周期性变化,不随时间变化发生幅值的偏移,与数学模型的物理意义十分吻合,说明谱元法是一种以较少节点获得高精度解的数值方法,具有很好的稳定性和收敛性。 ④利用Fourier-Legendre谱元方法深入地研究了晶体生长中加速坩埚旋转技术对于浓度均匀化的作用,数值结果显示:浓度均匀化最优的加速坩埚旋转模式是具有双向旋转的对称梯形模式,最优的无量纲时间周期为T=0.1;在本模型中溶液的完全混合是由扩散和对流共同决定的,而对流是从时间尺度上加快了溶液总体混合的过程,即改变了局部浓度梯度,从而加速扩散。 ⑤通过对方腔x-y平面上的边界施加无滑移边界条件的有限长模型的顶盖驱动流失稳特性分析发现:首先,对于立方体方腔流动,失稳的临界Re数和波数为818.41和3,流动的失稳属于静态失稳,失稳后的流动在壁面附近比方腔中心要明显得多,通过能量分析发现最危险的区域位于上流线固壁附近,这与展向施加周期性边界条件的无限长模型的数值是一致的,流动的失稳机制为与静态TGL(Taylor-Goertler-like)模式相关的离心失稳。其次研究了Λ为整数的几个算例,根据结果提出了两个判断不同Λ下流动失稳临界参数的预测,以Λ为小数且1Λ2的算例数值结果验证了上述预测的正确性;本文结果同直接数值模拟的数值结果进行了对比,两种思路得到的临界波数吻合得很好,仅Λ=1.8时存在不符,可能是数值方法的精度导致了这个差异,另外,不同Λ下的流动失稳均为静态失稳,且改变Λ并不改变流动失稳的物理机制。最后,当Λ=2π时,扰动的控制方程同基于无限长模型的扰动方程在形式上是一样的,而本文Λ=2π时的有限长模型失稳临界值同无限长模型的数值结果也符合地很好,且当波数k≤11流动为振荡失稳,即Hopf分叉,而当波数k11时流动为静态失稳,与文献结果十分一致;进一步验证了本文有限长模型的线性稳定性分析结果的正确性,能够用于解释流动的失稳机制。 ⑥通过对方腔x-y平面上的边界施加滑移边界条件的有限长模型的顶盖驱动流失稳特性分析发现:首先,立方体方腔流动的失稳临界Re为337.02,不到无滑移边界下的有限长模型和无限长模型临界Re的一半,说明了无滑移边界条件能够起到稳定流体流动的作用;流动的失稳属于静态失稳;通过直接数值模拟验证了临界波数k=2的正确性;基于能量分析发现,正的总能传递速率在方腔壁面附近达到顶峰,导致了流动的失稳,这是流动失稳的物理机制。由于角涡的消失,即缺失了形成TGL涡的主涡与下流线角涡之间的分界面,因此本模型下的流动失稳与无滑移边界下的有限长模型失稳是不同的。其次,展向长度的变化对于临界雷诺数的影响较小,且不同Λ下流动失稳均属于静态失稳;同样,由于施加了滑移边界条件,流动相对于无滑移边界下的流动更加地不稳定,因此,进一步说明无滑移边界条件能够稳定流体流动。
【学位授予单位】:重庆大学
【学位级别】:博士
【学位授予年份】:2012
【分类号】:O35

手机知网App
【相似文献】
中国期刊全文数据库 前10条
1 薛淑艳;张会强;王兵;王希麟;;空间发展超音速混合层入口扰动的最不稳定波方法[J];清华大学学报(自然科学版);2011年05期
2 彭岚;张伟;李友荣;孟海泳;;微重力条件下Marangoni数对分离结晶过程的影响[J];人工晶体学报;2011年03期
3 ;[J];;年期
4 ;[J];;年期
5 ;[J];;年期
6 ;[J];;年期
7 ;[J];;年期
8 ;[J];;年期
9 ;[J];;年期
10 ;[J];;年期
中国重要会议论文全文数据库 前10条
1 姚朝晖;郝鹏飞;何枫;;摆动射流机理研究[A];第七届全国实验流体力学学术会议论文集[C];2007年
2 程永光;;用Lattice Boltzmann方法模拟二维方腔流[A];第十三届全国水动力学研讨会文集[C];1999年
3 张靖煊;聂超群;陈静宜;;叶顶微喷气提高旋转畸变条件下轴流压气机失稳裕度的机理分析[A];大型飞机关键技术高层论坛暨中国航空学会2007年学术年会论文集[C];2007年
4 马东军;孙亮;秦丰华;孙德军;尹协远;;对称驱动方腔定常流动的多解现象研究[A];非定常空气动力学研讨会论文选集[C];2007年
5 孙智利;颜大椿;张汉勋;;自然对流边界层的转捩[A];中国力学学会学术大会'2005论文摘要集(上)[C];2005年
6 万德成;;用FEM-Level Set方法数值模拟Rayleigh-Taylor不稳定性问题[A];第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会论文集[C];2009年
7 曹洪建;查晶晶;万德成;;内置障碍物的方腔驱动流的数值模拟[A];现代数学和力学(MMM-XI):第十一届全国现代数学和力学学术会议论文集[C];2009年
8 刘兆存;肖勇;;二维平行壁面剪切流内部结构初探[A];第五届全国水动力学学术会议暨第十五届全国水动力学研讨会文集[C];2001年
9 梁晓波;张艺;张建伟;李世琼;程云君;;铸态Ti-22Al-25Nb合金热变形机理研究[A];动力与能源用高温结构材料——第十一届中国高温合金年会论文集[C];2007年
10 马汉东;;高速剪切流动及其光学传输效应研究[A];第五届全国青年计算物理学术交流会论文摘要[C];2008年
中国博士学位论文全文数据库 前8条
1 梅欢;谱元方法求解不可压缩流体流动及流动线性稳定性分析[D];重庆大学;2012年
2 张靖煊;畸变条件下轴流压气机叶顶间隙流对流动失稳影响的非定常机制[D];中国科学院研究生院(工程热物理研究所);2007年
3 祁志国;微尺度气体流动与传热的直接Monte Carlo方法模拟[D];中国科学院研究生院(工程热物理研究所);2007年
4 童志庭;轴流压气机中叶尖泄漏涡、失速先兆、叶尖微喷气非定常关联性的实验研究[D];中国科学院研究生院(工程热物理研究所);2006年
5 孙学卫;剪切稀化流体分层流动的线性稳定性研究和数值模拟[D];清华大学;2011年
6 徐国亮;三维边界层流动失稳与Bypass转捩模式研究[D];清华大学;2011年
7 莫东鸣;环形双层液池内热毛细对流的线性稳定性分析[D];重庆大学;2012年
8 朴明日;沿波纹壁面下落薄膜流动的动力学特性研究[D];上海大学;2011年
中国硕士学位论文全文数据库 前10条
1 李琳;压气机流动失稳DSP在线控制方法的研究[D];中国科学院研究生院(工程热物理研究所);2008年
2 徐建冬;时-空守恒元解元数值方法研究[D];南京航空航天大学;2009年
3 王沛;压气机叶顶泄漏流动的控制策略及其扩稳机理研究[D];中国科学院研究生院(工程热物理研究所);2007年
4 王永红;叶轮机械三维非定常流场数值方法研究及其并行策略[D];南京航空航天大学;2007年
5 商奇伟;环形浅液池内中等Pr数流体热毛细—浮力对流[D];重庆大学;2006年
6 张文杰;环形双层液体内热对流过程的二维数值模拟[D];重庆大学;2010年
7 黄婷婷;时空守恒元解元方法及其在气动噪声计算中的应用[D];南京航空航天大学;2012年
8 李碧仙;旋转管式微滤膜器内流型及分离过程实验研究[D];四川大学;2005年
9 王国斌;微尺度下油膜缝隙流动中的壁面滑移特性研究[D];北京工业大学;2011年
10 万钊;微可压缩模型预处理求解方法研究[D];中国空气动力研究与发展中心;2010年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026