用隐式无网格算法求解层流N-S方程

来源 :第十二届全国计算流体力学会议 | 被引量 : 0次 | 上传用户:johndovl1
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文通过粘性布点离散,研究求解N-S方程的隐式无网格算法,并应用于翼型的流场计算.采用无网格算法,计算区域用点云离散代替通常的网格划分;计算点上的空间导数,用当地点云上引入的二次极小曲面逼近.求解的N-S方程用隐式时间后差离散,结合用Roe的近似Riemann解确定无粘通量,并用LU-SGS算法分布计算,数值模拟了翼型亚、跨、超层流绕流.
其他文献
在过去二、三十年里,CFD得到了飞速发展.求解NS方程、复杂外形计算及湍流模拟是主要的目标.当前和未来发展的方向是什么?本文介绍一些CFD的新近进展情况.CFD(计算流体力学)与CSD(计算结构力学)相结合在模拟结构干扰问题方面会发挥巨大的作用,CFD软件工程化、标准化、可视化和通用化是其发展的必然趋势.CFD在描述真实气体效应方面、热喷流及其噪音、燃烧流动结构、壁面流、非平衡流等方面也会得到飞速
本文通过组合采用已有的数值技术,发展了一套可用于大攻角分离流动的数值模拟方法,通过不同算例的验算和与实验数据的比较,证明该方法具有较好的工程实用性.
运用人工压缩方法求解不可压缩二维Navier-Stokes方程,数值模拟了鱼类转弯机动运动中的尾鳍摆动的简化模型,探讨了其中引起转弯的机理.结果表明,鱼尾鳍一侧的来回摆动产生的流体动力力矩使得鱼体能够完成转弯角度.前摆阶段提供转弯对应的力矩,停留阶段保持冲量矩,后摆阶段抑制了持续的转弯;摆幅和停顿时间是控制转弯角度的两个重要因素.
本文首先详细介绍了有限体积的动能BGK方法:单元界面宏观物理量Vanleer校限制器进行TVD插值方法、采用线性展开和质量、动量和能量三大守恒定律构造粒子热运动速度分布函数,通过对分布函数取矩,得到单元面上的通量,对所得到的流体力学半离散控制方程,采用隐式一阶迎风得到全离散方程,隐式离散方程采用了LU-SGS进行求解.在此基础上,求解了典型的高超声速绕流流动,得到的壁面压力和热流,同实验测量非常吻
在单相流动CFD数值模拟中,对计算域出口进行流量守恒的修正是保证程序稳定和计算结果合理的有效手段.本文结合作者实际的模拟工作,探讨了在复杂的湍流两相流动,湍流两相反应流动中出口流量修正问题,给出了基于双流体模型的出流边界条件,为复杂反应流动的CFD模拟提供了有益的经验.
昆虫前飞需要满足这样的气动力约束条件:拍翼产生的升力与体重相平衡,推力与身体阻力相平衡.本文基于理论模化方法来研究拍翼运动的气动特性,不仅证明了昆虫前飞消耗的气动功率与前飞速度间的"J"型曲线关系,还通过对气动力随时间的变化曲线分析得出:升力主要由下拍行程产生,推力主要由上挥行程产生,而功率消耗也集中在下拍行程.另外,本文还在同一前飞速度(J=0.6)条件下,研究了不同拍翼幅值对气动性能的影响,结
阐明万亿次高性能计算机的体系结构特点和并行计算技术在CFD领域的应用情况,给出航空气动力计算软件WiseMan的并行算法、并行实现以及适应于高性能计算机发展而即将开展的相关技术改进.
建立了一个城市小区尺度模式.该模式为三维非静力高分辨,取k-ε湍流闭合方案.除考虑空气动力学作用外,引入了作为城市特征的街区建筑物布局及其高度、朝向和对短波辐射的遮蔽以及不同地表利用类型等特征影响,用强迫—恢复法计算地面温度,同时加入了污染物平流扩散方程.模式经专门布置的实测资料验证并用其模拟了北京一个小区中的气象条件和污染物扩散特征.结果表明:模式能够较合理的反映由于建筑物遮蔽造成的温度差异以及
发展了预处理时间导数项的双时间推进的非定常N-S方程解算器.在有限体积法的框架下对N-S方程进行积分并采用显式的Runge-Kutta方法以主隐式的LU-SGS方法进行时间推进.计算了低马赫数、低雷诺数下翼型以及机翼的定常和非定常空气动力特性.数值模拟了低马赫数、低雷诺数下做俯仰、沉浮运动的翼型绕流流场.数值结果显示,预处理可压缩N-S方程的方法可以模拟不可压缩的定常以及非定常流动以及雷诺数对低速
本文将加速收敛技术—预处理方法应用到双时间步法的框架之下,通过内层时间的预处理提高低速流动条件下解的精度和求解效率;同时,利用外层时间来逼近真实的物理非定常过程.在本文中,详细介绍了双时间步、预处理方法的实现过程.并用这种方法计算了无粘Bump流动和绕二维圆柱的非定常流动,数值结果表明该算法显著提高了低速时流场求解的效率及求解精度,对所模拟低速的非定常流动而言,所需代价只是传统时间推进方法的几分之