扩散界面浸入边界法结合壁面模型在湍流模拟中的应用

提出了一种模拟高雷诺数湍流的扩散界面浸入边界法（IBM）。该方法采用壁面模型来减少壁面附近的网格量。为了实施壁面模型和边界条件，在物面外部设置了2个辅助层（一系列拉格朗日点）：外侧的参考层用于实施壁面模型来得到壁面剪应力，内侧的强制层用于实施离壁的边界条件。在实施壁面模型时，将动量方程沿物面法向积分，从而把物面切向的速度修正量与由壁面模型得到的壁面剪应力联系起来，而速度的法向分量则按二次曲线分布来近似重构。在实施边界条件时，为了严格满足无穿透条件，应用了基于隐式速度修正的IBM。最后，利用绕平板湍流和绕NACA0012翼型湍流来对方法的可行性进行了验证。     

不可压缩Navier－Stokes方程的基本解法

利用不同阶次的Ｈｅｌｍｈｏｔｚ算子方程的基本解，建立了时间离散不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的完全边界积分表示式。对俯仰谐振的ＮＡＣＡ００１２翼型所计算得的升力和阻力的滞后现象，与实验结果相一致     

基于ANCF的松弛绳索动力学建模与仿真

传统绝对节点坐标方法(ANCF)绳索模型是基于梁单元建立的,其不能反映纤维绳索的不可抗压和松弛特性。考虑纤维绳索初始松弛余量,给出了绳索非线性轴向应力-应变关系。绳索在松弛状态下轴向力接近于零,在张紧状态下表现出线弹性特性。在此基础上,采用绝对节点坐标方法推导了松弛绳索动力学模型。通过静力学和动力学仿真,将传统绳索模型和松弛绳索模型进行对比,结果表明:在重力及不同集中载荷作用下,松弛绳索相比于传统绳索都具有更大的变形;分析移除集中载荷后绳索的动力学响应可知,传统绳索在振动过程中始终处于张紧状态,且绳索上各点振动同步,而松弛绳索会在张紧和松弛状态之间不断转换,绳索上各点的振动存在相位差,能够更好地反映绳索在松弛状态下的动力学特性。     

快速反设计翼型的涡源交替分布迭代法

本文在高阶板块法任意翼型分析程序的基础上,发展了给定表面压力分布设计翼型的反问题程序。根据表面涡分布对控制弯度有效而源分布对控制厚度分布有效的原理,提出了在反设计过程中交替采用改变涡分布和改变源分布,进行迭代。克服了一般算法单纯采用改变一种分布进行迭代的缺陷。算例结果表明,该法迭代次数少,结果好。     

NUMERICAL SIMULATION OF VORTEX BREAKDOWN ON A DELTA WING AT LOW SPEED

ＮＵＭＥＲＩＣＡＬ　ＳＩＭＵＬＡＴＩＯＮ　ＯＦ　ＶＯＲＴＥＸ　ＢＲＥＡＫＤＯＷＮ　ＯＮ　Ａ　ＤＥＬＴＡ　ＷＩＮＧ　ＡＴ　ＬＯＷ　ＳＰＥＥＤＺｈｕＺｉｑｉａｎｇ；ＪｉａＪｉａｎｂｏ（ＩｎｓｔｉｔｕｔｅｏｆＦｌｕｉｄＭｅｃｈａｎｉｃｓ，ＢｅｉｊｉｎｇＵｎ...     

后台阶流中自由旋涡的稳定性

最近的实验表明用侧壁小孔抽吸可在后台阶流中捕捉到静止的分离旋涡。笔者将用理论证明,二维后台阶流中存在自由旋涡的驻定态位置,并且有中性的驻定态稳定性。还进一步证明,用侧壁小孔抽吸,可使自由涡达到真正的驻定态稳定     

Numerical simulation of airfoil flows with a turbulence model

Turbulent flows over AS240 and NACA4412 airfoil were simulated numerically using a two-equation turbulence model named k-ξ model. The predictions of velocity profiles and the pressure coefficient of airfoil AS240 at 8°/19° attack angle and NACA4412 at 13.87° attack angle were calculated. The results were compared with those using k-ε and k-ω models,as well as experimental data. It indicates that the new k-ξ model offers more realistic prediction than the other two models. The main finding shows that the new k-ξ model is good at predicting separated flows around airfoils, and it captures the flow feature of pressure-induced separation adequately. All calculations are implemented as per openFOAM 1.7.1(open source field operation and manipulation).      

不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

不可压湍流反应流直接模拟和RANS反应模型检验

用谱方法对三维不可压槽道湍流反应流动进行了直接数值模拟,得到了温度和质量分数耦合的瞬态数据库.结果显示温度与质量分数脉动在近壁区都有条带结构.推导了雷诺平均方法中湍流二阶矩反应模型中所求关联量的精确输运方程,发现耗散项是化学反应率系数-质量分数关联量封闭的关键,耗散项需要考虑化学反应影响.基于数据库的统计结果,对关联量模型方程中各项进行了先验研究,发现产生项和耗散项的贡献最大,扩散项和反应项的贡献较小,化学反应对各项大小和分布形状有明显影响.在算例中,直接模拟统计得到的化学反应率系数-质量分数关联和用代数二阶矩模型的模拟值很接近,说明ASOM模型具有一定的合理性.      

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。