用于非定常粘性流动计算的无网格算法(英文)

把无网格算法发展用于求解涉及动边界的非定常粘性流动问题。在处理粘性流动的布点问题时,通过在远离物面的区域采用无粘流动计算时所用的各向同性点云,而在物面附近引入各向异性点云,从而在保证物面法线方向布点较密以准确模拟边界层的同时,有效控制了布点总量,减少了计算时间。对于动边界问题,本文在上述方法获得的初始布点的基础上,采用基于扰动衰减规律的点云移动技术,快速得到物面运动到其他位置时流场求解所需点云。在所获得的点云上,采用一种带权系数的二次极小曲面逼近方法来离散Navier-Stokes方程的空间导数,权系数的引入使得流场求解时对各向同性点云和各向异性点云无需分开考虑,而是采用统一的求解方法,从而简化了编程求解。用发展的无网格算法,结合Navier-Stokes方程求解双时间推进方法,并耦合Spalart-Allmaras湍流模型,先成功地模拟出NLR7301翼型后缘摆动诱发的非定常流,并通过与实验比较,验证了本方法.接着进行了NACA0012失速模拟,给出了动态失速过程,展示出用本文方法处理复杂动边界问题的效果。     

自适应笛卡尔网格超声速黏性流动数值模拟

复杂外形/流场的高质量网格的生成往往需要占用大量人力资源,而自适应笛卡尔网格方法能够自动化生成高质量网格,具有很好的工程实用价值和应用前景。基于笛卡尔网格方法,采用叉树数据结构进行数据的存储和访问,分别从几何特征和流场解特征出发进行网格的自适应加密和粗化,发展了一种二维情况下自动、高效的自适应笛卡尔网格生成方法。从浸入边界方法出发,结合虚拟镜像对称方法和曲率修正技术进行黏性物面边界条件的处理,同时建立了多值点问题的处理技术,发展了一种在笛卡尔网格下可有效模拟黏性物面边界条件的方法。针对自适应笛卡尔网格非均匀的特点,发展了悬挂网格的处理方法,并构建了适用于自适应笛卡尔网格的黏性数值求解器。通过典型算例的考核,验证了所发展的自适应笛卡尔网格生成技术和构建的数值求解器具有较高的精度和可靠性。     

有烧蚀气体产物影响的化学非平衡粘性流数值计算研究

本文是针对高超声速绕流钝锥体再入飞行器，研究有烧蚀气体产物影响的化学非平衡反应粘性流ＮＳ方程数值计算方法。通过带有化学源项的ＮＳ方程对流动作数值模拟，总的连续方程由单个组元守恒方程组代替。对１３个化学组元，２６个气相化学反应，以及２个气相－固相化学反应，研究烧烛碳－空气的化学反应系统，并与纯空气化学反应系统的计算结果进行比较。利用二阶迎风ＴＶＤ差分格式和激波捕捉法，计算球锥体头部和身部化学非平衡粘性流场     

FAST SOLUTION OF TURBULENT SEPARATION BY STRAWN''''S METHOD

ＦＡＳＴＳＯＬＵＴＩＯＮＯＦＴＵＲＢＵＬＥＮＴＳＥＰＡＲＡＴＩＯＮＢＹＳＴＲＡＷ＇ＳＭＥＴＨＯＤＹｉｎＪｕｎｆｅｉ（ＤｅｐａｒｔｍｅｎｔｏｆＰｏｗｅｒＥｎｇｉｎｅｅｒｉｎｇ，ＮＵＡＡ２９ＹｕｄａｏＳｔｒｅｅｔ，Ｎａｎｊｉｎｇ２１００１６，Ｐ．Ｒ．Ｃｈ...     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

小振幅振动圆柱在偏心圆域中诱导的二次定常旋涡流动

研究在充满粘性不可压缩静止流体的偏心圆柱域中当内柱在垂直于柱轴的平面内作小振幅简谐振动时诱导的二次定常旋涡流动。当Ｒｗ＝Ｕ＾２∞／ωγ为有限时量，Ｓｔｏｋｅｓ层次的二次定常整流流动必须用Ｎ－Ｓ方程来描述。     

地效飞行器三维地面效应的数值模拟

采用三维非结构贴体直角网格和N-S方程研究了三维地面效应的数值计算方法。用最小二乘法对N-S方程进行了MUSCL有限体积法的粘性求解,提出了带约束边界条件的最小二乘法的处理办法,修正了下垫面的数值边界条件,模拟了下垫面在低马赫数情况下对飞行器的地面效应作用。结果表明计算方法可以反映地面效应的基本原理和特性,模拟结果显示飞行器的地面效应不仅有正效应,而且可能有负效应。     

粘性不可压流动压力修正的有限元计算

将ＳＩＭＰＬＥ方法压力修正思想应用于有限元方法，对粘性不可压缩流动的ＳＵＰＧ有限元方程，建立将压力与速度分离，进行压力修正迭代的计算格式；对典型算例进行多个Ｒｅ数的层流，湍流计算，显示此方法具有良好的收敛性和精度。     

小振辐振荡圆柱诱导的二次旋涡流动

本文给出充满粘性不可压缩流体时同心圆柱间内柱作小振幅振荡所诱导的二次旋涡流动的解析解,对典型的柱径比给出容器外边界对二次旋涡流动流线、涡量分布和动能分布的影响。给出了等涡量线和等动能线分布的鞍点。采用数值求解非定常Stokes流动考察了二次旋涡流动的发展过程。二次定常旋涡流动的显示实验表明,解析解、数值计算和实验显示所示流谱基本一致。     

强激光与高超声速球锥流场干扰数值模拟研究

本文采用有限差分技术，数值求解带能量源项与空气11组分化学动力学模型的轴对称粘性流体力学方程组，采用类氢原子理论模型计算空气对激光的吸收系数，建立了模拟强激光与超声速钝体流场干扰的计算软件，得到了高能浙江与超声速球锥流场干扰数值模拟的结果。