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

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

动量矩方程的剪切积分解析表示式

本文采用两层涡粘性模型和速度剖面的指数律,导出了可压缩动量矩方程的剪切积分解析表示式,从而大大简化了边界层计算。和其它理论结果相比,本结果更接近于实验数据。     

Navier—Stokes方程粘性项的有限体积离散法研究

Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程粘性项的有限体积离散法研究封建湖西北工业大学应用数学系［摘要］：三维Ｎ－Ｓ方程是描述真实气体流动的非线性偏微分方程组，如何快速准确地计算其粘性项至关重要。这不仅因为粘性项的计算量非常巨大，而且粘性项对急剧变化的流动梯度的...     

高超音速化学非平衡层流粘性激波层钝头细长体绕流数值计算

本文研究了高超音速化学非平衡钝头细长体绕流的计算方法,给出了耦合求解连续方程和法向动量方程的具体做法,以及在高空情况下考虑滑移边界条件后的数学处理。计算结果表明,对钝头细长体,可以将计算推进到三十倍头部半径的后身区域。     

火箭水下发射复杂流场的近似数值模拟

本文对火箭水下直接点发射过程中，喷管内高温燃气的推进和气囊的形成、发展等流动演变过程提出了一类简化的模型，并进行了非定常粘性 可压缩流Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值模拟。计算结果表明，简化模型的合理性，对水下火箭发射的流动过程分析具有一定的参考价值和一种应用价值。     

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

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

再入钝头体近尾流动计算

本文采用联立求解粘性剪切层和有旋元粘区的方法计算了再入钝头体层流近尾流场，有旋无粘区采用有旋特征线法计算，粘性剪切层采用改进的流管法计算，这两个区域间的基本流线上的压力，温度，速度等由两个区域迭代匹配确定，粘性剪切层的下边界区分流线采用的实验确定，回避了回流区复杂计算，计算结果与有关试验和文献一致，本文的目的是提供远尾流计算所需的颈部初剖面。     

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

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

用高精度紧致差分格式分块耦合求解二维粘性不可压缩复杂流场

本文提出了张量形式的粘性不可压缩流场方程的高精度分块耦合求解方法。在分块算法中用重叠一层的方法来实现子块间的数据交换，同时克服了块剖分所造成的数值奇性。利用本文所提出的算法成功地模拟了椭圆柱绕流问题的时间发展过程，如尾流，旋涡的产生，脱落过程以及Karman涡街的形成等。