飞行器流动仿真讲稿第7章-拟一维喷管流动的数值模拟.pptVIP

飞行器流动仿真讲稿第7章-拟一维喷管流动的数值模拟.ppt

  1. 1、本文档共104页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
飞行器流动仿真讲稿第7章-拟一维喷管流动的数值模拟

程序编写举例 U1m(n)=U1(n) U2m(n)=U2(n) U3m(n)=U3(n) pm(n)=p(n) C dU1dtm=-(U2m(i)-U2m(i-1))/deltx cJ2=g7*pm(i)*(a(i)-a(i-1))/deltx dU2dtm=-(g5*(U2m(i)**2/U1m(i)-U2m(i-1)**2/U1m(i-1))+ + g6*(U3m(i)-U3m(i-1)))/deltx+cJ2 dU3dtm=-(g0*(U2m(i)*U3m(i)/U1m(i)-U2m(i-1)*U3m(i-1)/U1m(i-1)) + -g4*(U2m(i)**3/U1m(i)**2-U2m(i-1)**3/U1m(i-1)**2))/deltx C dU1dtav=(dU1dt+dU1dtm)/2.0 !计算时间导数的算术平均值 dU2dtav=(dU2dt+dU2dtm)/2.0 dU3dtav=(dU3dt+dU3dtm)/2.0 程序编写举例 if (icase.eq.2) then U1(i)=U1(i)+dU1dtav*dt !计算校正步变量值 U2(i)=U2(i)+dU2dtav*dt U3(i)=U3(i)+dU3dtav*dt else dU1m=U1m(i+1)-2.0*U1m(i)+U1m(i-1) dU2m=U2m(i+1)-2.0*U2m(i)+U2m(i-1) dU3m=U3m(i+1)-2.0*U3m(i)+U3m(i-1) dpm=abs(pm(i+1)-2.0*pm(i)+pm(i-1))/(pm(i+1)+2.0*pm(i)+pm(i-1)) s1m=Cx*dpm*dU1m s2m=Cx*dpm*dU2m s3m=Cx*dpm*dU3m C U1(i)=U1(i)+dU1dtav*dt+s1m !计算校正步变量值(含人工黏性) U2(i)=U2(i)+dU2dtav*dt+s2m U3(i)=U3(i)+dU3dtav*dt+s3m end if 程序编写举例 r(i)=U1(i)/a(i) u(i)=U2(i)/U1(i) t(i)=g1*U3(i)/U1(i)-g4*u(i)**2 p(i)=r(i)*t(i) am(i)=u(i)/sqrt(t(i)) fm(i)=r(i)*a(i)*u(i) if (i.eq.it) write(12,22) j,r(i),t(i),p(i),am(i) 400 continue C进口边界条件i=1 r(1)=1.0 !物理边界条件 t(1)=1.0 U1(1)=r(1)*a(1) U2(1)=2.0*U2(2)-U2(3) !一阶外推 u(1)=U2(1)/U1(1) U3(1)=U1(1)*(g2*t(1)+g3*u(1)**2) p(1)=r(1)*t(1) am(1)=u(1)/sqrt(t(1)) fm(1)=r(1)*a(1)*

文档评论(0)

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

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

1亿VIP精品文档

相关文档