initial condition and coding

Question

Basheer 2015-7-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/231754-initial-condition-and-coding

回答： Basheer 2015-8-1

Hello guys, I have the 1D channel as shown in the pic. I use finite volume method to code and I depend on number of cells (m) in my coding. I would like to add the water depth initial condition based on the distance (x) not (m), so the water depth is h=0.75 when x<15.5 else it is zero.

in the same time I wanna keep the using of the definition of h(i),as vector, my previous work without of interact of x ,of other cases,is:

 h=zeros(1,m);
 for i=1:m;
     if (i>m/2)
         h(i)=0.000001;
     else
         h(i)=1;
     end
 end

do u have any suggestions? thanks in advance.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Basheer 2015-8-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/231754-initial-condition-and-coding#answer_187932

I did it, thx

x=0; for i=1:m; if x<=25.5; Zb(i)=0; elseif x>25.5 && x<=28.5; Zb(i)=Zb(i-1)+dx*0.1333; elseif x>25.5 && x<=31.5; Zb(i)=Zb(i-1)+dx*(-0.1333); elseif x>31.5; Zb(i)=0; end x=x+dx; end