[재료 수치 해석] 수치해석 프로젝트 레포트와 수행결과 C 프로그램 소스

소개글

재료 수치 해석 레포트와 C언어 소스입니다.
소스는 파일에 첨부해 두었습니다.
문제는 아래와 같으며 양적해법과 음적해법으로 두가지 방식으로 C 프로그램을 구현하였으며
각각의 결과 및 실행 그래프를 캡처해서 올려 놓았습니다.

(문제) 열전달 문제
두께 100cm의 붉은 벽돌로 된 벽이 있다. 내면온도 275℃, 외면온도 40℃라고 가정할 때,
k=0.25, c=0.76, 벽의 밀도는 3.2g/cm3이다.
따라서,
L=100cm
편미분방정식 계수(alpha) =k/밀도/c= 0.103
초기조건 u( x, 0 )=275℃
경계조건 u( 0, t )=40℃

목차

1. 문제

2. 프로그램

3. 실행결과

4. 결과 그래프
- 양적해법
- 음적해볍

5. 느낀점

6. 참고문헌

본문내용

열전달 문제>
두께 100cm의 붉은 벽돌로 된 벽이 있다. 내면온도 275℃, 외면온도 40℃라고 가정할 때,
k=0.25, c=0.76, 벽의 밀도는 3.2g/cm3이다.
따라서,
L=100cm
편미분방정식 계수(alpha) =k/밀도/c= 0.103
초기조건 u( x, 0 )=275℃
경계조건 u( 0, t )=40℃

<중 략>

error_max = 0.0;
temper[mesh+1] = temper[mesh-1];
for(i=1; i < mesh+1; i++){
temp_new =temper[i+1].start-2*temper[i].start + temper[i-1].start;
temp_old =temper[i+1].old-2*temper[i].old + temper[i-1].old;
temper[i].new = temper[i].old +
r*(rel_coeff*temp_new+(1-rel_coeff)*temp_old);
error = fabs(temper[i].start-temper[i].old);
if(error > error_max)

참고 자료

없음