(精品word)拟一维喷管流动的数值解法(MATLAB)代码.doc资源-CSDN文库

需积分: 9 200 浏览量 2022-10-23 23:31:41 上传评论收藏 19KB DOC 举报

资源详情

资源评论

资源推荐

(精品 word)拟一维喷管流动的数值解法(MATLAB)代码

拟一维喷管流动的数值解法（MATLAB）代码

数值计算代码

%拟一维喷管流动的数值解

%亚声速-超声速，非守恒形式

function main()

clear；

clc;

r=1。4; ％绝热指数

N=1001； %时间步长

i=31；％网格数目

L=3；％喷管长度

C=0.5; %柯朗数

dx=L/(i-1）; %空间步长

dt（N）=0； %时间步长

x=linspace(0,L，i)；％网格点横坐标

A=1+2。2*（x—1。5).^2； %喷管面积

%赋值

M(N,i)=0；

T（N,i）=0;

V（N，i)=0；

%初始条件

M(1，：）=1—0.3146＊x；

T（1,:)=1—0.2314＊x;

V（1，：）=（0.1+1。09*x）。*(1—0。2314＊x)。^0。5；

％按时间步长推进

for k=1:N—1

％预估偏导数

M_t(1:i—1)=—V(k，1：i—1）.＊(M（k，2：i）—M(k，1：i—1）)/dx-M(k，1：i-1)。＊（V（k，

2:i）-V(k,1：i—1））/dx-M(k，1:i-1）。*V(k，1：i—1).*log（A（2：i）./A（1：i-1)）/dx;

V_t(1:i-1）=—V（k，1：i—1）。＊(V（k,2：i）-V（k，1:i—1))/dx-1/r.＊((T(k,2:i)—T（k,1：

i-1））/dx+T（k,1:i—1）。/M（k，1:i-1).*（M（k,2：i）-M(k,1：i-1))/dx)；

T_t(1：i-1）=—V(k，1：i—1）.＊(T(k，2:i）-T（k,1：i—1））/dx-(r-1).*T（k,1:i—1）。＊

（（V(k，2：i）-V（k,1:i-1)）/dx+V(k，1:i-1).*log(A（2:i)。/A(1:i—1））/dx）;

％求取内部网格点处最小时间步长

t=C*dx。/（V(k，2：i-1)+sqrt（T(k,2：i—1）)）;

dt(k）=min(t);

%预估值

M1（1：i—1）=M(k，1：i—1）+M_t(1:i-1)＊dt（k);

V1（1:i—1)=V(k，1：i—1）+V_t(1：i—1）*dt（k）；

T1（1:i—1）=T（k,1：i-1）+T_t(1：i—1）＊dt（k)；

%校正偏导数

M_t_1 （ 2 ： i-1 ） =-V1 （ 2 ： i—1). ＊（ M1(2 ： i-1 ） -M1(1:i-2) ） ./dx—M1(2:i-1 ） .*(V1(2 ：

i—1)-V1(1:i-2））./dx—M1（2：i-1)。*V1（2:i-1).＊log（A（2:i—1)。/A(1：i—2）)。/dx；

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余2页未读，立即下载

评论0

内容反馈

智慧安全方案

粉丝: 3650
资源: 59万+

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈

feedback-tip