自适应笛卡尔网格湍流计算研究.pdfVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第十四届全国计算流体力学会议论文集 自适应笛卡尔网格湍流计算研究 刘周,周伟江 (中国航天空气动力技术研究院,北京100074) 捅 要:复杂外形的粘性网格生成已经成为计算流体力学在工程应用中的主要瓶颈,高效而鲁棒的 网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡尔网格方法生成计算网格,投影方 法拟和壁面边界,将投影得到的柱形单元沿法向分层获得适于粘性计算的网格。流场解算中应用最小二乘 一方程湍流模型加入以获得湍流解。LU-SGS格式隐式时间推进,对网格重排序保证格式的稳定性和平衡 性.算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于 实际工程计算。 关键词:笛卡尔网格;sA湍流模型;Lu—SGS格式;网格重排序 O引言 现代飞行器的外形已经越来越复杂。在飞行器的前期设计阶段中,快速,可信的气动力计算对 飞行器的选型工作起到重要的作用,这对网格生成和流场求解提出了新的要求。实践经验表明,几 何处理与结构网格生成的时间总和约占整个CFD计算周期的8096甚至更多,而在选型阶段,网格生 成上的大量时间消耗是不合算的。因此,非结构网格越来越受到人们的重视。 非结构网格中的一大类是三角形(二维)和四面体(三维)网格。经过过去20余年的发展,此 类网格的生成技术已经相当的完善并在工程上得到了大量的应用n~。然而,由于网格本身的特点, 此类网格并不适于粘性的计算,常用的解决办法是在壁面通过拉伸得到棱柱形网格捕捉边界层,在 外围使用三角形或四面体网格幅’3。。 非结构网格中的另一大类是笛卡尔网格。它非常容易生成,与四面体非结构网格相比,计算花 销和存储需求较少,然而由于笛卡尔网格是从“空间到边界”的方法,故它的一个重大缺陷是不能生 成贴体网格,因此对边界的精确表示是一个成功的笛卡尔网格方法最重要的部分。常用的边界处理 方法有浸入边界方法(ImmersedBoundary)m1,嵌入边界方法(EmbeddedBoundary)M和投影方法‘删。 本文的研究采用投影方法处理边界,将笛卡尔流场单元和物面之间的空腔用基于最近距离的投 影方法填充为柱形网格,同时得到新的表面网格。将投影形成的柱形网格进一步划分为多层,自然 的得到粘性边界层网格。 隐式时间推进n到。应用重编号技术¨副对单元进行重编号以保证格式的稳定性和平衡性,并提高计算 效率。 1 网格生成方法 、 以一个类CAV外形为例细述网格生成方法。使用三角形来离散初始的表面模型。这种表面描述 大量的用在包括CFD在内的多种领域内,得到各种CAD包的广泛支持。如图l。 1.1笛卡尔网格生成 使用1个或多个矩形单元覆盖整个计算域。将那些与物面相交的单元称为切割单元。通过对切 422 刘周,周伟扛:自适应笛卡尔月格湍流*算Ⅻ究 割单元的一分为八加密来得到整个计算网格,如图2所示。 本文使用(1)最大加密次数;(2)单元的尺寸;(3)局部曲率自适应三种加密方式。对第三种 加密,是指笛卡尔网格单元包含的表面三角形序列中任意两三角形法向矢量夹角的余弦小于某个指 定值时,该单元应该加密,如F式, “ia(n,_)‘‘。《射,Viu,f≠, (1) 在加密时,除了切割单元需要加密外.与切割单元相邻的几层单元也应该加密,使得阿格层次 之佩的过渡更加光滑。在加密过程中,应避免相邻单元之间的层次大于l的情况出现。 氆唇衙 目1#!《Ⅲ=自*目 目2*{$月*∞gH8 网格生成过程中使用八叉树的数据结构,父.子单元之间的继承关系使得同格的遍历和单元邻 居的查找都非常方便。另个重要的数据结构是存放表面三角形数据的ADT数据结构”1,它对有限 尺寸目标的近似搜索非常有效。利用这两种数据结构,能够极大提高在单元切割判断时的效率。 计算中仅用剜流场单元,切割单元和内部单元应被屏蔽。切割单元在阿格生成过程中己被标定, 对剩下的单元.首先用射线法判定某个单元的类别,然后使用P

文档评论(0)

july77 + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档