收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

基于DLR型k-ε模型的AMG方法对渐扩管路内紊流数值仿真

蒋光彪  
【摘要】:对紊流的数值仿真是流动、传热传质中最基本的课题.何永森等研究者通过一系列数值实验发现,用DLR型k—ε紊流模型·BFC法(边界拟合曲线坐标变换法),能够对总扩散角为8~0、扩散度为4的锥形渐扩管路内完全发展的不可压粘性紊流场较精确地数值仿真(又称数值模拟,计算机仿真)。多重网格方法(Multi-Grid Method,简称MGM)也称为多格子方法或多层网格法,是求解偏微分方程(组)大规模离散化方程最有效的方法,它一般可以分为几何多重网格方法(Geometric Multigrid Method,简称CMC)和代数多重网格方法(Algebraic multigrid method,简称AMG)。近年来,由于大型计算机的迅速发展和功能的日趋完善,从而使得多重网格方法做为最优的算法将理论用于现实。由于实际应用问题的错综复杂性,以及数值商业软件对“即插即用”型求解器的要求,使得几何多重网格方法的应用变得越来越困难,而代数多重网格方法的高效性和稳健性(robustness“鲁棒性”)使之成为了当今多重网格方法的研究热点. 本文针对边界拟合曲线坐标贴体网格系统下代数多重网格方法在紊流数值预测中的应用问题,在何永森教授工作的基础上,详细探讨了代数多重网格方法在数值预测紊流中的实施过程,并围绕如何提高锥形渐扩管路内紊流数值仿真效率与精度展开研究,进行了一系列数值实验及与物理实验结果的比较,得到了一些有意义的成果,丰富和充实了代数多重网格方法,拓宽了代数多重网格方法在一些领域中的应用。具有理论和工程应用价值.主要内容和成果包括: (一)原来的程序软件判断数值收敛与否是根据数值计算结果的情况分次沿时间步推进求解并以人工方式进行.本文对这一情况进行了改进,通过添加了一个判断收敛的子程序,变人工方式为自适应控制收敛方式.数值实验表明,用这种方式控制迭代次数,当迭代结束时,各种流动参数的误差(前后两时间步的差)曲线处于平稳并且很小,计算结果与实验结果较好符合。 (二)将AMG方法引入到紊流数值预测的有限差分计算领域,详细探讨了代数多重网格方法在紊流数值预测中的实施过程.研究了不同有限差分格式下,紊流模型离散得到的大规模代数系统对应系数矩阵的整合及边界条件的嵌入方法,讨论了系数矩阵的“三元组”压缩存储方式,利于节约内存及便于利用AMG方法求解。编制了代数多重网格方法数值求解紊流模型离散得到的大规模代数离散系统的程序接口;将代数多重网格方法与DLR、DHR型k—ε紊流模型.BFC法结合,应用编制的程序DLRAMG和DHRAMG对锥形渐扩管路内紊流进行了数值预测,将数值计算结果与物理实验结果进行了比较,多重网格方法的计算结果与实验结果符合较好。并且与POINT-SOR方法相比,可以节约近三分之一的CPU时间,提高了数值预测效率。 (三)针对原来十三点格式所用的模板节点多,方程组系数矩阵的带宽大,非零元多,数值求解费时问题,设计了一种基于五节点模板的新五点差分格式,并将其与AMG方法结合,进一步提高了紊流数值预测的效率.数值实验结果表明,在AMG方法求解的条件下,新五点差分格式比原十三点差分格式可以节约近三分之一的CPU机时。 (四)本文提出了一种可以同时实现控制网格正交性和任意控制边界网格间距的一种BFC网格生成的新方法。该方法可对生成的网格边界间距大小任意控制,同时生成的BFC网格还具有边界及内部较好的正交性.应用实例的计算结果表明,该方法能够对复杂边界的单连通域或多连通域生成较理想的BFC网格。 (五)另外,研究了k-ε紊流模型对数值预测正弦波壁流动的应用。基于有限体积法结合非正交同位网格系统,压力与速度耦合采用SIMPLE方法,对该种流动进行了数值预测,得到了与实验结果符合较好的计算结果。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 谢伟,陈材侃,程尔升;用于曲线坐标生成的多重网格与有限分析法[J];华中理工大学学报;1995年05期
2 刘之行,封卫兵;三维区域上的多重网格算法[J];西安交通大学学报;1998年12期
3 卫加宁,王伟沧,皮新明,章社生;一类虚拟边界预条件多重网格并行算法[J];武汉理工大学学报(交通科学与工程版);2001年01期
4 葛永斌,吴文权,田振夫;二维波动方程的加权平均隐格式及多重网格算法[J];上海理工大学学报;2002年03期
5 欧阳洁;求解跨音速势流的自适应并行多重网格算法[J];计算机工程与科学;1995年04期
6 葛永斌,吴文权,卢曦;解二维扩散方程的高精度多重网格方法[J];工程热物理学报;2002年S1期
7 刘星,卞恩荣,童谷生;非结构网格上二维欧拉方程的多重网格解[J];华东交通大学学报;1999年03期
8 葛永斌,吴文权,卢曦;求解二维扩散方程的加权平均隐式多重网格方法[J];内蒙古大学学报(自然科学版);2003年05期
9 田小波,吴庆举,曾融生;弹性波场数值模拟的隐式差分多重网格算法[J];地球物理学报;2004年01期
10 冯树常,彭爱和;一个在球面半规则网格上的求解Helmholtz方程的多重网格法[J];气象科技;1994年02期
11 徐正权;多重网格区域分裂分布式计算[J];数值计算与计算机应用;1996年01期
12 纪欢;王宁宁;;二维抛物方程初边值问题的拟多重网格预处理迭代法[J];科技信息;2011年05期
13 李桦,王承尧,王磊;用多重网格TVD方法求解三维高超音速粘性流场[J];国防科技大学学报;1998年02期
14 石钟慈,许学军;Wilson非协调元的V循环多重网格法[J];中国科学A辑;1999年10期
15 石钟慈,许学军;非对称不定问题的非协调多重网格法[J];计算数学;1999年04期
16 刘学强,伍贻兆,夏健;多重网格法在非结构网格中的应用[J];计算物理;2002年04期
17 曾晓清,吴雄华,王云飞;有限体积法与多重网格法用于提高激波捕捉质量及加速收敛[J];应用数学与计算数学学报;1995年01期
18 曾晓清,吴雄华,王云飞;多重网格法与有限体积法在流体力学中的应用[J];同济大学学报(自然科学版);1995年05期
19 黎明,张政;计算流体力学与传热中的多重网格块修正算法[J];北京化工大学学报;1999年03期
20 葛永斌,吴文权,卢曦;基于二维泊松方程六阶紧致格式的多重网格方法[J];上海理工大学学报;2002年04期
中国重要会议论文全文数据库 前10条
1 周磊;谭伟伟;;聚合多重网格的体网格可视化技术[A];探索 创新 交流(第4集)——第四届中国航空学会青年科技论坛文集[C];2010年
2 袁礼;;三维不可压Navier-Stokes方程的多重网格求解[A];第十届全国计算流体力学会议论文集[C];2000年
3 王青平;白武明;王洪亮;;多重网格在二维泊松方程有限元分析中的应用[A];中国科学院地质与地球物理研究所第11届(2011年度)学术年会论文集(上)[C];2012年
4 王斐;赵杰煜;;基于多重网格的水平集图像分割方法[A];第十四届全国图象图形学学术会议论文集[C];2008年
5 刘刚;刘伟;牟斌;肖中云;;涡流发生器绕流模拟的数值方法研究[A];庆祝中国力学学会成立50周年暨中国力学学会学术大会’2007论文摘要集(下)[C];2007年
6 葛永斌;田振夫;吴文权;;三维非定常对流扩散方程的高精度全隐式多重网格方法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
7 朱培烨;;基于非结构网格的欧拉方程聚合多重网格法[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
8 李笑郁;毛士艺;;基于相位信息加权的合成多重网格干涉SAR相位展开算法[A];中国体视学学会图像分析专业、中国体视学学会仿真与虚拟现实专业、中国航空学会信号与信息处理专业第一届联合学术会议论文集[C];2000年
9 吴爱弟;牟永光;杨文杰;曹思远;符力耘;;微分方程的小波解法[A];1995年中国地球物理学会第十一届学术年会论文集[C];1995年
10 陈明;柯涛;丁大志;樊振宏;陈如山;;电磁散射问题中的等级基二重网格法[A];2007年全国微波毫米波会议论文集(上册)[C];2007年
中国博士学位论文全文数据库 前10条
1 刘学强;基于混合网格和多重网格上的N-S方程求解及应用研究[D];南京航空航天大学;2002年
2 阳述林;高维中子输运方程的离散格式与并行算法研究[D];中国工程物理研究院;2003年
3 冯国峰;波动方程反问题的多尺度—信赖域反演方法[D];哈尔滨工业大学;2006年
4 谭勤学;p型多重网格间断有限元及其在燃机冷却系统中的应用[D];清华大学;2013年
5 张挺;基于多点地质统计的多孔介质重构方法及实现[D];中国科学技术大学;2009年
6 李壮;探地雷达反问题的同伦算法研究[D];哈尔滨工业大学;2007年
7 付朝江;集群MPI环境下有限元结构分析并行计算研究[D];上海大学;2006年
8 葛长江;基于交错、同位和非正交网格求解流场的半隐式快速算法研究[D];吉林大学;2010年
9 李文军;块代数多重网格算法的研究及其实现技术在在油藏数值模拟软件中的应用[D];中国科学院软件研究所;2001年
10 牛海峰;流体及非线性最优控制问题的有限元方法:状态受限与超收敛分析[D];山东大学;2011年
中国硕士学位论文全文数据库 前10条
1 李林;一维抛物型方程拟多重网格预处理迭代法及数值试验[D];渤海大学;2013年
2 杨磊;基于PC机群环境MPI的多重网格并行算法研究[D];武汉理工大学;2004年
3 刘鑫;基于PDE的几种图像处理方法及其数值解法[D];大连理工大学;2005年
4 杨艳南;二维抛物型方程初边值问题拟多重网格预处理迭代法[D];渤海大学;2013年
5 赵振宇;求解第一类算子方程的多重网格算法[D];河北工业大学;2004年
6 陈飞;基于多重网格的图像恢复算法[D];福州大学;2006年
7 万建军;椭圆边值问题非协调有限元多重网格法分析[D];郑州大学;2004年
8 陆华才;基于摩擦学理论的缸套激光表面微造型参数优化的研究[D];江苏大学;2005年
9 张杨;图形处理器并行计算应用研究[D];西南交通大学;2006年
10 王军华;多尺度计算方法的初探及应用[D];西北工业大学;2006年
中国重要报纸全文数据库 前1条
1 张伟;管理出“妙方” 治理全覆盖[N];泸州日报;2009年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978