微尺度圆管内气体流量的实验测量

以氮气和氦气为工作介质,研究了微尺度(φ=17.6μm,17.9μm)圆管内压力驱动的气体流动,测量了在不同进出口压力下的流量,并与计算值进行了比较.结果表明,即使对于极低Mach数(M=0.003～0.03)流动,其可压缩性仍十分明显,它来源于微尺度效应.在本文的实验参数下,Kundsen数范围为3×10-3～7×10-3,表面滑移不明显,稀薄气体效应可以忽略.     

叶栅可压缩流场的数值分析

将一维ＡＵＳＭ＾＋格式与三阶精度ＭＵＳＣＬ格式融合,给出了其在任意风线坐标下的二维形式,采用有限体积法,将ＡＵＳＭ＾＋格式与Ｒｕｎｇｅ－Ｋｕｔｔａ格式结合,对叶栅可压缩流场进行了数值模拟,典型算例的计算结果与有关献结果相符很好,表明ＲＫ－ＡＵＳＭ＾＋混合格式对叶栅可压缩流场进行数值模拟是可行的。     

三维多组分可压缩驱动问题的分数步特征差分方法

本文提出三维多组分驱动问题的分数步长特征差分格式,并应用变分形式,能量方法,粗细网格配套,叁二次插值,高阶差分算子的分解和乘积交换性理论和技巧,得到最佳阶L^2误差估计,该方法已成功应用到油资源评估,强化采油数值模拟和海水入侵预测和防治的数值模拟中。     

高分辨率激波捕捉格式及其CFD应用

在曲线坐标下将高分辨率AUSMPW激波捕捉格式扩展到三维问题。为了提高计算精度,采用三阶MUSCL格式,将AUSMPW格式与隐械时间推进法LU－SGS方法结合,应用于可压缩无粘和有粘流动的数散模拟,计算结果与文献中计算结果和实验数据相符,并将该格式与CD、TVD、FVS和FD格式进行了比较。     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法,它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点,该方法被有效应用于充填流动过程的数值模拟.     

类橡胶材料变形直到破坏的显式本构模型

本文通过直接、显式的方法提出一个多轴可压缩应变能弹性势来模拟类橡胶材料受载荷直到软化破坏的变形行为．首先，我们提出一个多轴可压缩应变能函数；其次，通过特定的不变量，该多轴应变能函数在单轴拉伸，平面应变和等双轴拉伸三个基准实验的情况下，可以退化为各自的单轴形函数形式；再次，我们显式给出带有软化破坏特性的形函数；最后，模型结果和试验数据可以精确匹配，同时可以预测材料临近破坏以后，接下来的变形行为．     

圆柱面上微小后向台阶层流绕流场的数值研究

1前言钩体绕流是燃烧、传热领域基本的流动现象，具有流动分离、流态转被、旋涡脱落、涡致振动等复杂特性。降低钩体绕流阻力一直是人们关心的问题。自从Prandtl圆球绊线实验以来，Stuper，Roshko，Zukausks，Achenbach，Klebanoff等［‘一列从不同角度研究了钝体绕流阻力的构成因素及降低钝体绕流阻力的若干方法。一致的观点是：因绕流钝体必然伴随着流动分离，流动阻力可分为与粘性有关的摩擦阻力及与流动分离有关的压差阻力两部分，其中，压差阻力约占总阻力的90N95％，其大小与流动分离点的位置密切相关。若流动分离时边界层内为层流…     

不可压缩非牛顿流体非定常流动初边值问题的非线性半群方法

本应用非线性半群理论构造了不可压缩非牛顿流体非定常流动初边值问题的强解。     

一种捕捉可压缩流体多重交汇界面的改进型LS方法

 在普通Level Set（LS）方法基础上，给出了一种捕捉多重交汇界面的改进型LS方法，该方法能够大大减小多重交汇界面之间因光滑而产生的空穴区域，而不影响多重交汇以外区域的界面位置。同时对改进算法进行了测试，并与可压缩流体耦合，给出了多介质流体多重交汇界面的计算结果。     

水下爆炸过程的高精度数值计算

采用三阶精度的Parabolic Piecewise Method(PPM)方法和Volume Of Fluid(VOF)方法相结合，运用Lagrange-Remapping算法，数值计算爆轰产物、水、空气多种介质之间的相互作用过程，该方法可以用来处理界面两边高密度比可压缩流动以及强剪切滑移运动等问题。