收藏本站
《湘潭大学》 2007年
收藏 | 手机打开
二维码
手机客户端打开本文

多介质可压缩大变形流体及辐射热传导数值模拟研究

张瑗  
【摘要】: 本文致力于研究惯性约束聚变(ICF)内爆压缩过程数值模拟中需要解决的两个关键问题,一是探索内爆过程中遇到的多介质可压缩大变形流体及流体界面不稳定性的高精度Euler数值模拟方法,二是寻求内爆过程数值模拟所涉及的三温辐射热传导方程的更为有效的数值解法。 本文首次系统地探索了高阶加权本质上无振荡格式(weighted essentially nonoscillatoryscheme,简记为WENO)用于ICF内爆压缩过程数值模拟的实际可能性。对于ICF内爆过程所遇到的各种流体不稳定性问题及多介质可压缩大变形流体问题,本文将一律使用经典的高阶WENO方法或者针对具体疑难问题设计基于高阶WENO的新的计算方法进行计算和测试,并与低阶方法进行比较,最终目标是试图在国内外率先研制一个基于高阶WENO方法的高精度ICF内爆程序,为我国国防和现代化建设服务。这既是本文的主旨,也是本文的主要创新之一。本文主要工作如下: (1)提出了求解辐射热传导方程的一类无网格方法(见本文第二章)。这项工作首次将无网格方法成功地应用于求解ICF数值模拟中经常遇到的二维三温方程。理论分析和大量数值试验表明,本文提出的方法在不规则网格上的计算精度以及对于网格不规则程度的适应性均明显优于目前在国内应用最广的九点差分格式。由此可见,本文提出的无网格方法不仅对于求解几何边界复杂的非线性抛物型方程有重要实用价值,而且在ICF数值模拟研究领域具有广阔应用前景。 (2)针对多介质可压缩大变形流体,提出了一类基于准守恒γ-模型或准守恒体积分数模型的高阶混合型有限体WENO方法,简记为FV-WENO-MT,而且在求解控制方程的同时,我们用水平集方法(或者Lagrange方法)为多介质问题提供了一张清晰的物质界面(见本文第七章)。这是国内外此前无人涉足的一项创造性成果。这项工作不仅解决了流体混合型方法界面不十分清晰的问题,同时突破了二阶精度方法的局限,为各种高密度比、强激波、大变形、具有多个界面且界面拓扑结构变化十分复杂的高维问题的高精度计算开辟了新的途径。 (3)用高阶有限差分WENO格式(简记为FD-WENO)模拟了ICF内爆过程中经常遇到的各种高密度比Rayleigh-Taylor不稳定性(简称RT不稳定性),激光烧蚀RT不稳定性及高马赫数Richtmyer-Meshkov不稳定性(简称RM不稳定性),均获得了令人满意的数值结果,而且通过进行定性及定量比较,表明高阶FD-WENO格式的确明显优于已有的低阶方法,如MUSCL和PPM等(见本文第三章和第四章)。由此证明高阶FD-WENO格式用于内爆压缩过程数值模拟不仅是可行的,而且是十分可取的。 (4)在综述多介质可压缩流体界面捕捉方法及界面追踪方法的基础上,对基于高阶FD-WENO格式的流体体积方法(volume of fluid,简记为VOF),基于高阶FD-WENO格式的水平集方法(level set)及界面跟踪方法(front tracking)进行了数值测试、比较和理论分析,表明后二种方法效果很好,但第一种方法欠佳(见本文第六章)。从而得出了一个有指导意义的结论:在使用高阶FD-WENO格式进行ICF数值模拟时,应尽量避免使用VOF方法捕捉界面。 (5)为了避免辐射流体力学方程组的表达形式过于复杂和降低求解难度,其中的三温能量方程通常采用内能(而不是总能量)作未知函数,这就导致整个方程组不是守恒型的,以致各种经典的求解双曲守恒律组的数值格式,如高阶FD-WENO格式等在这里不能直接使用。为了解决这个矛盾,我们于本文第五章专门设计了两类用于求解以内能作未知函数的Euler方程组的高阶FD-WENO方法,数值试验表明计算效果较好,这两类新的方法明显优于文献【29】中所介绍的方法。
【学位授予单位】:湘潭大学
【学位级别】:博士
【学位授予年份】:2007
【分类号】:O241.8

手机知网App
【参考文献】
中国期刊全文数据库 前3条
1 柏劲松,陈森华,李平,张展冀;多介质可压缩流体动力学界面捕捉方法[J];爆炸与冲击;2004年01期
2 张学莹;赵宁;朱君;;多流体界面不稳定性的守恒和非守恒高精度数值模拟方法[J];爆炸与冲击;2006年01期
3 符尚武,沈隆钧,黄书科;激光驱动内爆二维数值模拟方法及其程序(Lared-I)的研究[J];高技术通讯;1998年07期
【共引文献】
中国期刊全文数据库 前10条
1 张飞,汪继文;一维溃坝问题的间断Galerkin方法[J];安徽大学学报(自然科学版);2005年02期
2 邓义华;;一类非线性微分方程单支θ-方法的稳定性[J];安徽大学学报(自然科学版);2006年04期
3 王保国;郭延虎;刘秋生;沈孟育;;HIGH-ORDER ACCURATE AND HIGH-RESOLUTION UPWIND FINITE VOLUME SCHEME FOR SOLVING EULER/REYNOLDS-AVERAGED NAVIER-STOKES EQUATIONS[J];Acta Mechanica Sinica;1998年01期
4 王汝权;申义庆;;SOME WEIGHT-TYPE HIGH-RESOLUTION DIFFERENCE SCHEMES AND THEIR APPLICATIONS[J];Acta Mechanica Sinica;1999年04期
5 姜宗林;;ON DISPERSION-CONTROLLED PRINCIPLES FOR NON-OSCILLATORY SHOCK-CAPTURING SCHEMES[J];Acta Mechanica Sinica;2004年01期
6 ;A robust algorithm for moving interface of multi-material fluids based on Riemann solutions[J];Acta Mechanica Sinica;2006年06期
7 ;Non-linear characteristics of Rayleigh-Taylor instable perturbations[J];Acta Mechanica Sinica;2008年02期
8 ;Numerical investigation on the flowfield of “swallowtail” cavity for supersonic mixing enhancement[J];Acta Mechanica Sinica;2009年01期
9 姚兴军;张凤阳;王正东;章文俊;;基于VOF和CSF方法的变截面毛细流动数值建模[J];半导体技术;2011年02期
10 涂侯杰,恽寿榕,赵衡阳;破爆型串联战斗部第一级爆炸对第二级影响的研究[J];兵工学报;1994年03期
中国硕士学位论文全文数据库 前10条
1 侯建军;高温矿井热环境数值模拟及热害控制技术研究[D];河南理工大学;2010年
2 李蒙;特殊构型靶对激光质子加速特性的影响[D];中国工程物理研究院;2010年
3 任健;辐射热传导方程九点格式的研究与应用[D];中国工程物理研究院;2010年
4 吕信;激光能量沉积光路追踪法及其并行化[D];中国工程物理研究院;2009年
5 李连举;轴流式通风机气动噪声分布特性的数值研究[D];山东科技大学;2010年
6 罗富强;扁管管壳式换热器数值模拟研究[D];哈尔滨工程大学;2010年
7 侯鹏;基于双曲型方程RKDG法在流场及声场中的数值模拟研究[D];哈尔滨工程大学;2010年
8 高慧颖;基于CFD的大型船舶球艏线型设计比较研究[D];哈尔滨工程大学;2010年
9 郭涛;偏流板性能的分析与计算[D];哈尔滨工程大学;2010年
10 刘丽莉;电磁场驱动下的液体流动实验与模拟研究[D];辽宁工程技术大学;2009年
【二级参考文献】
中国期刊全文数据库 前1条
1 符尚武,沈隆钧,黄书科;激光驱动内爆二维数值模拟方法及其程序(Lared-I)的研究[J];高技术通讯;1998年07期
【相似文献】
中国期刊全文数据库 前10条
1 娄路亮,曾攀;影响无网格方法求解精度的因素分析[J];计算力学学报;2003年03期
2 曾清红,卢德唐;耦合径向基函数与多项式基函数的无网格方法[J];计算物理;2005年01期
3 向勇,冯斌,马驰,赵春茁,王成程,陈骥;原型装置真空靶室真空形变分析研究[J];真空与低温;2005年03期
4 张占文;李波;唐永建;王朝阳;漆小波;陈素芬;刘一杨;翟世明;;玻璃微球壁厚和预充气工艺对气体渗透性能的影响[J];强激光与粒子束;2006年06期
5 谢琴;刘更;刘天祥;;再生核粒子方法及其应用[J];机械科学与技术;2006年07期
6 温万治;董师舜;杭旭登;胡晓燕;;二维MLSPH无网格方法[J];计算物理;2008年02期
7 陈光南,常铁强,张钧,张兴宏,裴文兵,尤希文;激光与靶非平衡耦合的数值模拟[J];计算物理;1998年04期
8 史宝军,袁明武,李君;基于核重构思想的最小二乘配点型无网格方法[J];力学学报;2003年06期
9 曾清红,卢德唐;含有启动压力梯度的渗流问题及其无网格解法[J];计算力学学报;2005年04期
10 柴立群;徐建程;许乔;;加窗后波前功率谱密度的计算值修正[J];强激光与粒子束;2005年12期
中国重要会议论文全文数据库 前10条
1 吴兴;李志鹏;王昌生;;流固耦合问题的非边界适应方法及无网格方法的研究进展[A];第二十三届全国水动力学研讨会暨第十届全国水动力学学术会议文集[C];2011年
2 刘寿先;陈光华;黄文斌;袁树云;王德田;李泽仁;彭其先;;差分干涉 CT 技术在界面不稳定性研究中应用[A];中国光学学会2011年学术大会摘要集[C];2011年
3 陈勇;黄建明;陆明万;;地震波传播计算中的无网格方法[A];1999年中国地球物理学会年刊——中国地球物理学会第十五届年会论文集[C];1999年
4 代晓青;曾竟成;肖加余;江大志;邢素丽;;水平集方法模拟反应性树脂体系一维流动界面[A];复合材料:创新与可持续发展(下册)[C];2010年
5 敬闰宇;孙婧;李梦龙;;使用改进的水平集方法对数据进行分类[A];中国化学会第28届学术年会第14分会场摘要集[C];2012年
6 倪国喜;王瑞利;林忠;;无网格方法中的粒子初始均匀分布[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
7 王丽丽;;界面不稳定性和湍流混合现象的数值模拟研究[A];中国工程物理研究院科技年报(2010年版)[C];2011年
8 刘高联;李孝伟;;基于当地笛卡儿架构的无网格方法[A];中国力学学会学术大会'2005论文摘要集(下)[C];2005年
9 雷雨冰;袁亚雄;;无网格数值方法研究及其在流体力学中的应用[A];加入WTO和中国科技与可持续发展——挑战与机遇、责任和对策(下册)[C];2002年
10 侯中喜;梁剑寒;王承尧;;高阶WENO格式的分区并行计算[A];第十届全国计算流体力学会议论文集[C];2000年
中国重要报纸全文数据库 前4条
1 记者 侯志明 冯昌勇;我国“氢弹之父”的晚年生活[N];新华每日电讯;2003年
2 特约记者 苏扩善、记者 范炬炜;我国“863”胜利推进15载[N];解放军报;2001年
3 本报记者 李鹏;贺贤土:秘密研制核武器[N];北京科技报;2009年
4 组稿 王毅俊;激光 让生活更美好[N];上海科技报;2010年
中国博士学位论文全文数据库 前10条
1 张瑗;多介质可压缩大变形流体及辐射热传导数值模拟研究[D];湘潭大学;2007年
2 李坤;基于无网格方法的声学问题数值模拟研究[D];华中科技大学;2011年
3 任红萍;插值型无网格方法研究[D];上海大学;2010年
4 刘箐;水平集方法在预混合湍流V型火焰模拟和形状复原问题上的应用[D];浙江大学;2010年
5 SIDIBE Kalilou(西迪贝 卡里鲁);接触碰撞问题的非线性无网格方法及在金属成型分析中的应用[D];湖南大学;2005年
6 李九红;复变量无网格方法及其应用研究[D];西安理工大学;2004年
7 Mohammed Aman Ullah(阿曼);用守恒型跟踪算法数值模拟流体中的介质界面[D];上海大学;2011年
8 林颖;基于水平集方法的图像分割关键技术研究[D];哈尔滨工程大学;2010年
9 葛全文;ICF中流体不稳定性的数值模拟研究与Hamilton-Jacobi方程的运动网格方法[D];中国工程物理研究院;2003年
10 吴俊峰;收缩几何中惯性约束聚变流体不稳定性的理论和数值模拟研究[D];中国工程物理研究院北京研究生部;2003年
中国硕士学位论文全文数据库 前10条
1 屈小妹;PPMLR格式用于惯性约束聚变中流体不稳定性数值模拟研究[D];湘潭大学;2007年
2 王晓光;弹性及弹塑性问题的无网格局部边界元法[D];浙江工业大学;2003年
3 丁慧芬;对流扩散方程的径向基函数法[D];辽宁师范大学;2008年
4 黄娟;改进广义移动最小二乘近似的无网格方法[D];苏州大学;2007年
5 马修彦;弹性问题无网格局部边界元法及其应用研究[D];浙江工业大学;2003年
6 苏芳;基于局部多顶式近似空间的单位分解方法及其误差分析[D];湘潭大学;2005年
7 孙阳光;自然边界元的无网格方法和拟小波方法研究[D];华中科技大学;2005年
8 郑兴;光滑质点流体动力学(SPH)算法研究[D];哈尔滨工程大学;2005年
9 郑志成;基于Level Set方法的散乱数据曲面重建[D];大连理工大学;2005年
10 张小华;无网格方法在计算流体力学中的应用研究[D];西北工业大学;2006年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62791813
  • 010-62985026