CFD研究及CFD软件应用及平台开发.doc

  1. 1、本文档共11页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
CFD研究与CFD软件应用和平台开发 高 智 中国科学院力学研究所 摘要: CFD软件已得到相当普遍的应用,不少单位正在开发CFD专业平台。在CFD软件应用和CFD平台开发中进一步发挥CFD研究的作用,处理好三者之间的正确关系,形成三者之间互有借鉴,彼此促进的良性循环,这是一个重要的课题,课题的解决需要CFD研究和应用工作者的共同关心和努力。本文讨论这方面的某些问题。 引言,讨论本话题的起因 气动中心牵头的《十·五》联合公关重大项目“实用的气动数值计算系统开发和实用算法研究”的宗旨是改进和提高已有CFD平台性能和档次,一个核心措施是把国内有关单位(三家和三家外)的好算法,好模型,好的网格生成等吸收进CFD平台中。我提出的数值摄动方法(包括摄动有限差分和有限体积方法)、离散流体力学算法也在被吸收之列。参加此项目二年多,对“CFD研究与软件应用和平台开发”的关系问题,与有关专家有过一些讨论,这个关系问题很重要,因此在这里讲讲,目的是抛砖引玉。 CFD软件应用中的某些问题 CFD软件已有很多,如CFX,Flow3D,Fluent,Phoenics,Star-CD,Fine等。“国内设计部门普遍对气动计算软件既缺乏信任,又抱怨不好用。而西方先进国家的情况并非如此,美AIAA1998.6月__发“Guide for Verification and Validation of CFD Simulations, AIAA G-077-1998)同时AIAA J.出版了“Credible CFD Simulations”,标志着CFD软件得到设计部门的使用和认可。 应用CFD软件要消化并吃透它,但不能局限于应用,发现问题,用CFD研究解决问题,进而反馈到平台开发中,例如气动力计算已得到认可,但摩阻和热流计算很成问题。对此,我们从算法,边界处理和网格生成作些讨论。 实用算法:有限体积方法 以模型方程为例 一阶迎风格式,什么都能算,数值耗散与物理耗散之比为 二阶中心格式,精度、效率、简单、适用性好等,综合考虑为首选, 网格雷诺数大产生数值振荡。 Jameson正型格式 这些格式中积分近似均取中点规则为二阶精度,重构近似为直接差分。 多基点高精度格式:包括积分近似和重构近似高精度。边界处理复杂,代数方程复杂,难构造,内点精度高,边界点精度低。 Derivatives Derivatives Compact method Multi-node method Direct difference and PFD method Space点 X-△X X X+△X 数值计算与物理分析相结合的一些问题 4.1.摩阻、热流计算和边界处理:算好摩阻、热流要求 ? ; 边界用单边差分近似、精度为一阶和二阶时,数值耗散与物理耗散之比约为0.50和0.167. 因此的要求不够准确。 高Re数NS方程计算的若干物理问题研讨 §1 引言 高Re数流动计算的历史概貌 理想无粘流方程 ?D’Alembert佯谬 边界层理论(无粘、边界层渐进匹配求解理论,解决了D’Alembert佯谬,把粘性流研究从佯谬笼罩的领域引入科学领域) 扩散抛物化Navier-Stokes(DPNS)方程理论(无粘-边界层匹配求解失效,空间推进求解,三维?二维) 完全NS方程计算 LISTNUM a ,,, ,, ,, , ,, 讨论: 算不了NS时,“2”和“3”可认为是“不得已而为之”。虽有历史作用,且DPNS已与Euler,NS一起作为流体力学基本守恒方程在教材和CFD软件中得到使用。NS计算已成为主流。NS方程是很准确的,但已有NS方程对摩阻和热传导率的计算问题不少,不同计算结果之间彼此相差很大,需要从物理和数值方面加以探讨。 §2 完全NS方程计算的某些问题 2.1 现况与评估 张涵信等的估计:NS方程组(1.1)与m阶精度差分格式等价的修正方程是 LISTNUM b 令, LISTNUM b LISTNUM b NS计算要求 LISTNUM b 讨论: 以?x=?y=10-2,m=2为例,最大ReL?104,若ReL104,为Euler方程计算;边界层法向网格局部加密,为DPNS计算 提高m,取m=3,4,最大ReL?106,108 在固壁邻域,m?2;有限体积方法 总之按(2.1)—(2.4) ,只要提高m即可实现高Re数NS计算 张涵信等“网格与高精度差分计算问题”,力学学报1999,31(4):398-405 2.2 非均匀流动(剪切流)NS计算对网格间距,精度的要求 均匀场分析 剪切流分析(高智 “强粘性层流理论及在粘性流计算中的应用” 空

文档评论(0)

189****6140 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档