复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

类升力体外形高超声速粘性绕流拓扑结构分析

采用交替方向隐式分解的隐式NND格式求解全N-S方程模拟“类升力体”外形在高超声速下的大攻角流动,给出了“类升力体”外形表面极限流线随攻角变化的拓扑结构及40°攻角下垂直于体轴的横截面流线拓扑结构。结果表明:类升力体外形三维流场结构十分复杂,攻角从0°～40°变化时,背风面表面极限流线依次由不分离、开式分离向起始于鞍、结点组合的高阶奇点的分离方式转化,翼面横向分离亦随攻角增大而增大;垂直于体轴的横截面流动拓扑结构与文献[2]给出的理论分析一致。     

第四届亚洲计算流体力学会议学术交流情况

第四届亚洲计算流体力学会议于２０００年９月１９日至２１日在四川绵阳富乐山国际大酒店举办．出席会议的有１２６人,除我国大陆外,还有日本、韩国、印度、新加坡、美国、德国、英国、芬兰、瑞典等国家以及台湾、香港、澳门等地区的代表,远远超出了亚洲的范围．会议发表论文１０４篇,其中大会邀请报告９篇,专题邀请报告１２篇,其它论文８３篇．内容涉及计算流体力学的前沿方向和课题,重点领域的研究进展以及在网格技术、计算格式、并行计算、复杂流动及湍流等方面的新成果,也包含计算流体力学在航空航天及其它工业部门的应用研究．…     

二维粘性不可压缩流动的通用分离判据

本文研究了定常、非定常二维粘性不可压缩流动在分离点附近的流线性状,给出了判定分离的通用判据。并指出:M.R.S判据、奇点分离判据、零涡判据等仅适用于部分分离情况,不是普遍的分离判据。     

有关力学工作的一些反思

近年来,我国的一部分力学工作者产生了一些困惑,他们感觉力学学科受到的重视程度在下降,远不如上世纪五六十年代.     

求解Boltzmann模型方程的气体运动论大规模并行算法

研究各流域三维流动问题的Boltzmann模型方程计算方法,建立直接求解分子速度分布函数的气体运动论耦合迭代数值格式;基于变量依赖关系、数据通信与并行可扩展性分析,使用区域分解并行化方法,建立气体运动论数值算法并行方案,发展求解各流域三维绕流问题的气体运动论并行算法.拟定高低不同马赫数下来自不同流域的三维球体及返回舱绕流算例,进行高性能Fortran(HPF)大规模并行计算,将计算结果与有关实验数据、相关理论预测等进行比较分析,研究揭示不同流区复杂绕流现象及流动机理.研究表明,所发展的气体运动论并行算法具有很好的并行独立性,基本达到线性加速的并行效果,显示出良好的并行可扩展性.     

A mixed explicit-implicit antidiffusive method of navier-stokes equations for supersonic and hypersonic separated flows

The supersonic and hypersonic.laminar and turbulent vi-scous separated flows over two and three-dimensional com-pression corners are solved by the antidissipative explicit-implicit difference scheme presented in this paper.The re-sults obtained show that a high accuracy has been achieved bythis method and the time needed for computation is greatlyreduced.     

数值研究平板方舵激波-湍流边界层干扰

数值研究了平板方舵激波-湍流边界层干扰流场。模拟出了分离激波与弓型激波砬撞后形成的“λ”激波结构;消晰地显示了分离区中的旋涡结构,发现流场中会出现二次分离涡,并从理论上分析了流场对称面涡心形态与非定常的关系,得到了涡心为不稳定螺旋点或出现极限环是非定常流动特征的新结论。