激波冲击V形界面重气体导致的壁面与旋涡作用及其对湍流混合的影响

基于多组分混合物质量分数模型,采用色散最小耗散可控的高分辨率有限体积方法,数值模拟了弱激波冲击V形空气/SF_6界面后,界面不稳定性生成的旋涡与固体壁面作用问题.激波冲击V形界面之后,因斜压效应诱导涡量沉积在界面附近,形成沿界面规则排列的多个涡对结构.旋涡的诱导作用使界面不断变形和卷起,同时旋涡之间不断发生相互并对,诱导更多更小尺度的旋涡产生.旋涡诱导作用的叠加效应,使界面尖端处的初始涡对向上下壁面发展.随后,涡结构开始与壁面发生复杂的相互作用.旋涡与壁面作用后沿壁面加速,使得物质界面沿壁面伸展,随后,旋涡从壁面回弹,并诱导二次旋涡产生.旋涡与壁面相互作用的过程,能够明显加剧物质混合.本文从物质混合的角度研究了该过程的机理,分析了旋涡与壁面作用对物质混合的影响.     

求解可压缩流的高精度非结构网格WENO有限体积法

提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析.     

基于特征理论的二维无粘Lagrange流体力学有限体积法

提出一个求解二维无粘Lagrange流体力学方程的中心型有限体积方法.采用特征理论求解网格节点处的速度及压力,并利用这些物理量更新节点位置及计算网格界面通量.方法适用于结构网格与非结构网格.典型数值实验的结果表明,格式具有较好的收敛性、对称性、能量守恒性及鲁棒性,且能自然地求解多物质流动问题.     

基于特征理论的二维可压缩流动的二阶拉氏算法

提出一种求解二维拉氏可压缩流体力学方程的中心型二阶精度有限体积方法.利用特征理论构造网格节点处的局部近似演化算子,算子用来求解网格节点处的速度及压力,利用这些物理量更新节点位置及计算网格界面通量.通过结合一定的重构方案,该方法达到时、空二阶精度,并且形式简单、计算量小,适用于结构网格与非结构网格.典型数值实验表明,本文格式具有良好的收敛性、对称性及鲁棒性,且能自然地求解多物质流动问题.     

数值求解不可压缩流动的投影方法研究进展

在过去的20多年中, 投影方法通过速度和压力的解耦计算, 获得了比全耦合方法更高的计算效率, 这个显著优点使之得以广泛应用. 目前, 在计算非定常不可压缩流动的原始变量形式的数值方法中, 投影方法得到了越来越广泛的应用. 本文根据投影方法的构造思路,将众多的投影方法分成了3类, 即: Helmholtz-Hodge分解类投影方法、算子分裂类投影方法和局部连续投影方法, 并详细的介绍了3类投影方法的发展历程和求解步骤. 从投影方法的求解过程不难发现, 通过速度和压力的解耦计算, 提高了投影方法的计算效率, 但同时也给投影方法的时间精度分析带来了困难, 并长期成为大家争论的焦点. 普遍认为, 速度的时间精度比较容易达到高阶, 但是压力一般来说只有一阶精度. 但通过对3类投影的对比分析后, 我们认为, 局部连续投影方法将有助于澄清目前投影方法存在的相关争议, 并使得发展高阶精度的投影方法在理论上和技术上成为可能.      

基于样条逼近有限体积法的通量分裂算法

研究了任意曲线坐标系中求解双曲型守恒律的高精度、无波动样条逼近有限体积方法,比较了三种不同的通量分裂技术在这种方法中的应用。数值实验表明,在不同的通量分裂技术下,该方法对流场中的激波和接触间断都有很高的分辨率。     

庆祝中国力学学会成立50周年大会暨中国力学学会学术大会'2007会议简介

1 会议概况 值此中国力学学会50华诞之际,为回顾中国力学的光荣历程、展现我国力学的优秀成果、弘扬中国力学的优秀传统、展望新世纪力学学科的发展趋向,中国力学学会成立50周年暨中国力学学会学术大会'2007(Chinese Conference of Theoretical and Applied Mechanics-2007,CCTAM'2007)于2007年8月20～22日在北京召开.大会由中国力学学会主办,39个单位协办.     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或"红斑"现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.