分享
分享赚钱 收藏 举报 版权申诉 / 9

类型二维超声速流动的数值解——普朗特迈耶稀疏波.doc

  • 上传人:精品资料
  • 文档编号:8856136
  • 上传时间:2019-07-14
  • 格式:DOC
  • 页数:9
  • 大小:538.10KB
  • 配套讲稿:

    如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。

    特殊限制:

    部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。

    关 键  词:
    二维超声速流动的数值解——普朗特迈耶稀疏波.doc
    资源描述:

    1、二维超声速流动的数值解普朗特迈耶稀疏波主程序数值计算%二维超声速流动的数值解普朗特迈耶稀疏波clc;clear;%常量设置CFL=0.5;Cy=0.2;r=1.4;R=287;%初始物理条件Ma0=2;P0=1.01*105;M0=1.23;T0=286.1;%流场几何外形描述sita=5*pi/180;H=40;L=65;%网格节点描述i=6501;j=4001;x=linspace(0,L,i);point1=length(find(x=0.00001elsexm=xm(2) xm(2) xm(2);endn=n+1;endMa(1,k+1)=xm(2);%下边界处物理量修正值P(1,k+

    2、1)=P(1,k+1)*(1+(r-1)/2*Ma_cal2)/(1+(r-1)/2*Ma(1,k+1)2)(r/(r-1);T(1,k+1)=T(1,k+1)*(1+(r-1)/2*Ma_cal2)/(1+(r-1)/2*Ma(1,k+1)2);M(1,k+1)=P(1,k+1)/R/T(1,k+1);if kpoint1V(1,k+1)=0;elseV(1,k+1)=-U(1,k+1)*tan(sita);endF1(1)=M(1,k+1)*U(1,k+1);F2(1)=M(1,k+1)*U(1,k+1)2+P(1,k+1);F3(1)=M(1,k+1)*U(1,k+1)*V(1,k+1)

    3、;F4(1)=r/(r-1)*P(1,k+1)*U(1,k+1)+M(1,k+1)*U(1,k+1)*(U(1,k+1)2+V(1,k+1)2)/2;kend计算结果后处理close all;%横向速度point=5*floor(1:20:130);figurefor i=1:length(point)subplot(1,length(point),i) plot(U(:,point(i),y(:,point(i),-,linewidth,2)ylabel(y/cm)xlabel(u/(m/s)title(strcat(x=,num2str(point(i)axis(660 760 -15 4

    4、0)% axis offend%纵向速度point=5*floor(1:20:130);figurefor i=1:length(point)subplot(1,length(point),i) plot(V(:,point(i),y(:,point(i),-,linewidth,2)ylabel(y/cm)xlabel(v/(m/s)title(strcat(x=,num2str(point(i)axis(-250 50 -15 40)% axis offend%马赫数point=5*floor(1:20:130);figurefor i=1:length(point)subplot(1,l

    5、ength(point),i) plot(Ma(:,point(i),y(:,point(i),-,linewidth,2)ylabel(y/cm)xlabel(Ma)title(strcat(x=,num2str(point(i)axis(1 3 -15 40)% axis offend% 密度point=5*floor(1:20:130);figurefor i=1:length(point)subplot(1,length(point),i) plot(M(:,point(i),y(:,point(i),-,linewidth,2)ylabel(y/cm)xlabel(M/(kg/m3)

    6、title(strcat(x=,num2str(point(i)axis(0.5 1.5 -15 40)% axis offend% 压力point=5*floor(1:20:130);figurefor i=1:length(point)subplot(1,length(point),i) plot(P(:,point(i),y(:,point(i),-,linewidth,2)ylabel(y/cm)xlabel(P/(Pa)title(strcat(x=,num2str(point(i)axis(0.1*105 2*105 -15 40)% axis offend%温度point=5*floor(1:20:130);figurefor i=1:length(point)subplot(1,length(point),i) plot(T(:,point(i),y(:,point(i),-,linewidth,2)ylabel(y/cm)xlabel(T/(k)title(strcat(x=,num2str(point(i)axis(200 300 -15 40)% axis offEnd计算结果(粘性 0.2+角度 5 度)

    展开阅读全文
    提示  道客多多所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    关于本文
    本文标题:二维超声速流动的数值解——普朗特迈耶稀疏波.doc
    链接地址:https://www.docduoduo.com/p-8856136.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    道客多多用户QQ群:832276834  微博官方号:道客多多官方   知乎号:道客多多

    Copyright© 2025 道客多多 docduoduo.com 网站版权所有世界地图

    经营许可证编号:粤ICP备2021046453号    营业执照商标

    1.png 2.png 3.png 4.png 5.png 6.png 7.png 8.png 9.png 10.png



    收起
    展开