一类二阶守恒单调重映算法

在大变形流体力学问题的数值模拟中，经常会涉及到计算网格的重分。基于不同网格的物理量传递便是所谓的重映技术。重映算法是任意拉格朗日一欧拉方法的重要组成部分，本文描述了一类适用于任意网格的二阶守恒单调重映算法。该算法分为网格内物理量的多项式重构、近似积分计算、物理量的单调修正三个部分。本文采用了最小二乘法的思想构造了网格内物理量的梯度，并且通过对物理量的单调修正，保证了算法的精度和单调性。该算法计算公共表面上的网格间质量的改变，所以可以在任意网格上使用。通过一系列数值算例验证了该算法，并说明本文给出的算法是有效可行的。     

气囊充气过程流固耦合数值模拟

基于LS-DYNA有限元分析软件,采用任意拉格朗日-欧拉(Arbitrary Lagrangian-Eulerian,ALE)流固耦合方法模拟了气囊充气过程,获得了充气过程中,气囊外形和流场之间的动态关系,分析了囊内流场的动态变化情况,并和气囊分析的控制体积(Control volume, CV)法进行了对比.结果表明,CV方法尽管简单,但对于高速下的气囊充气过程,无法获得囊内瞬态的流场变化及气囊表面的展开细节,确立了ALE方法在气囊充气仿真的可行性;通过囊面结构应力的分析,获得了气囊展开时的危险位置.该仿真结果对于协助气囊设计,减少试验次数具有重要意义.     

基于SNGR方法的二维空腔噪声数值模拟

基于SNGR方法,采用有限差分法求解带源项的线化欧拉方程。数值方法采用色散关系保持(DRP)格式。低耗散色散龙格-库塔格式显式时间推进,采用无反射远场边界条件。对二维空腔,长深比为9∶1,在亚声速、跨声速和超声速多个马赫数下进行噪声预测,并与实验结果对比,两者吻合较好。     

随机微分系统的耗散性

提出了有关Ito型随机微分系统耗散性理论的新概念:按模耗散、按模等度耗散和按模一致耗散,并利用Lyapunov方法,借助于Ito微分公式沿着Ito型随机微分系统的解对所构造的Lyapunov函数求导数,给出了Ito型随机微分系统有关按模耗散理论的一些代数判据,获得了与确定性常微分系统耗散性理论相对应的结论,最后的算例证明了该方法的有效性和可行性。     

一种新型颤振激励系统的数值计算研究

针对一种新型颤振激励系统的特点,采用非结构旋转运动网格技术,通过数值求解非定常欧拉方程,计算了该激励系统工作时的非定常气动力特性,并与准定常计算结果进行了比较。研究结果表明:采用非定常方法所得结果与准定常计算结果有明显的差别;气动力的变化与激励系统中圆柱的旋转方向有关,旋转角速度越大,气动力变化的幅值越小;气动力的幅值与圆柱缝道宽度也有一定关系。     

带油箱结构的机身框段坠撞仿真分析

 主要研究了任意拉格朗日/欧拉耦合方法计算带油箱的机身框段的坠撞过程。首先,建立了带油箱结构的机身框段坠撞分析模型,包括机身框段结构模型和欧拉流体模型。采用水代替油箱内部燃油,考虑了不同装水量对机身框段耐撞性的影响。分析了坠撞过程中的液体晃动和泼溅与机身坠撞响应的影响,给出了装水量和机身框段最大垂向压缩位移、最大过载和能量吸收等参数之间的关系。通过仿真分析,揭示了冲击载荷作用下油箱内部燃油量对机身框段各个部分结构的损伤破坏的影响。研究指出,在垂直撞击环境下,机身加强框和蒙皮是主要的吸能结构。在油箱内装水量多的情况下,油箱结构的吸能作用不能忽略。给出了应急着陆或可生存坠撞事故发生前,飞行员所应采取的紧急措施。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法（Ⅱ）

建立了二维振荡机翼的非定常反问题的有限差分解法,并通过若干算例对求解方法进行了验证.首先根据基于欧拉方程的非定常气动反问题数学模型确定了具体求解方案,接着重点讨论了壁面非定常反问题边界条件和远场非定常无反射边界条件的处理方法,编写了相应的计算程序,并以绕1/4弦长点做小幅俯仰周期振动的NACA-65-0012翼型为初始机翼,在不同目标压强的条件下进行了气动反问题计算得到新翼型,它们的压强分布与预设目标相吻合,这表明本文的反问题求解方案与计算方法都是可行的.     

有限体积法求解任意回转面叶栅叶型反问题设计的欧拉方程

以有限体积法为基础,探讨了用θｍ坐标系下积分型Ｅｕｌｅｒ方程求解任意回转面上叶栅反问题的方法,应用这一技术设计叶型,能用较少的机时使型面上气流压力值满足预先给定的值。这种叶型设计方法考虑了任意回转面上旋转角速度、半径变化、流片厚度变化等多种影响因素。文中还详细讨论了其基本方程、数值格式和求解过程,给出了根据给定的压力分布进行静子叶栅和转动叶栅设计的验证算例,结果表明该叶型设计方法是有效的。     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.