数值模拟运动激波对物体的冲击效应.pdfVIP

数值模拟运动激波对物体的冲击效应.pdf

  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文档。上传文档
查看更多
数值模拟运动激波对物体的冲击效应 王惠玲 乐嘉陵 (中甯空气动力研究与发展中心,621000) 摘要本文采用显式二阶迎风TVD格式数值模拟激波绕流流场,计算T--维、兰维情况下运动激波对物 体的冲击效应.二维冲击波效应计算了激波绕两个方块、长方块运动的流场,三维冲击波效应计算了激 波绕~个方柱运动的流场.由于网格生成采用以二维、三维Laplace方程为变换式的微分方法,且采用分 块网格计算.所以网格能适应物面形状.计算结果显示了包括激波、接触问断和涡在内的详细的流场结 构. 关键词:TVD格式、数值模拟、激波、冲击效应 一、引言 通过数值模拟的方法研究激波沿管道的运动和运动激波对各类物体的冲击效应在工 程上十分重要。由于在激波绕流流场中存在间断面、描性干扰、分离涡等复杂的物理现象, 对于这样复杂的流场的数值模拟。差分法具有更多的优越性。流场中存在激波时,一般采 用两种方法进行数值计算,即激波装配法和激波捕捉法.激波装配法将激波作为边界米处 理,计算精度及分辨率都比较高,但必须事先知道激波的位置。对于复杂的激波绕流流场 更有效的是激波捕捉法。最初人们使用一阶格式捕捉激波.分辨率较低,而二阶格式在激 波附近的非物理振荡比较大。八十年代初,Harten首先提出了保单调格式即TVD格式 Variation (Total dimini$1fing)111,12]。使激波捕提法有了大发展,提高了计算精度,基 本上解决了高阶精度差分解在激波附近的非物理振动问题。以后TVD格式有了进~步发展 pI|H】-I”.形成了迎风型和对称型两种TVD格式。本文采用二阶显示迎风TVD格式数值求解 层流N—s方程,较好地捕捉到流场中的激波。国内外应川TVD格式解决运动激波绕流问题 因受计算机条件及经费的限制一般物形比较简单且多计算二维绕流问题。本文计算了二维、 三维激波效应,二维计算得到了详细的流场结果,三维计算结果比较简单。本文网格生成 \、、采用微分方法,并采用分块技术,使网格更好地适应物面,且物形复杂时程序仍适用。 、)、 二、二维冲击波效应的数值模拟 l控制方程采用无量纲二维N-S方程。 。 ‘ 户b 口。 上 L 无量纲量分别为:万2去·西=iu,可=言,i=主,歹=专,芦=寿0口。 彤o‘ 3 1 一:—;,=:乓,Re =——丁, 2—下t 为气体常数。月、仲吊裂· p口o‘0 C10。 R 无量纲状态方程为:=p;亍 p 2计算格式 采用文献[6]中的预估矫正法,具体步骤为 w吒一砜At(丘一一厂n 00~芳k;“t;刁 (1) +尝(E,j*t,t-E^u)+尝(彤”一彤”~,) U1.t‘1’=△%,I‘1’+q,I“ (2) +割尝(鹳“Ⅲ。∞川.。)+茜(掣”一甜”卅一,)} (3) ’ 一 u,《1=△E。t仁’+%.t。’ (4) E,、F。为数值通量.采用

文档评论(0)

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

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

1亿VIP精品文档

相关文档