for문을 통해 총 300개의 속도장 데이터를 읽어들이고, 각각의 파일을
% 읽음과 동시에 상기 선언된 변수(sum_u, sum_v) 위로 누적 합계를 덮어쓰면서 300
% 프레임의 속도 벡터 총합을 구한다.
for i = 1:9
filename = strcat('2조000',int2str(i),'.vec');
[x y u v null] = textread(filename,'%f%f%f%f%f','headerlines',7);
sum_u = sum_u + u;
Streamline를 보면 거의 Re가 매우 작은 상황과 비슷한 모양을 나타내고 있지만, 순간속도장의Streamline을 보면 위아래로 요동치는 Streamline을 볼 수 있다. 이를 Average 시켰기 때문에 Laminar flow와 비슷한 모양으로 나타난 것으로 생각 할 수 있다.
6. 3번에서 구한 Re의 값으로 이 유동이 Laminar인지 Turbulen
한다. 여기서 적절한 velocity는 Time-averaging된 속도장에서 자유 유동 속도(Free stream velocity)에 가장 가까운 값을 선택하면 된다. 자유 유동 속도는 무한한 유동장 내에서 아무런 장애를 받지 않는 유속을 말한다. 즉, 1장에서 구해진 속도장에서 실린더의 영향을 가장 적게 받는 유속을 선정하여야 한다.
for count=2:1:300
filename=strcat(int2str(count), '.vec'); %파일 이름은 0001->1로 바꾸어주었다.
[x y u v empty] = textread(filename,'%f %f %f %f %f','headerlines',7);
if flag0==0
Sum_u=zeros(length(u),1);
Sum_v=zeros(length(v),1);
flag0=1;
end
Sum_u=Sum_u+u;
Sum_v=Sum_v+v;
end
%평균 속도를 구한다.
Avg_u=Sum_u./299
6차 근사 다항식이기 때문에 hue값이 120을 넘어갈 경우 그에 대응하는 온도 값은 의 자리수를 갖는 값으로 추출이 된다. 이는 매우 이상한 결과이기에 두 번째 그래프에서는 이러한 부분의 값은 모두 0으로 대체하여 그렸고 이 경우 첫 번째 방법으로 구한 그래프보다는 조금 더 신뢰성이 높다고 할 수 있