열전달 5-43.docx



열전달 HEAT AND MASS TRANSFER 4th Edition SI Units.

Fundamentals and Applications

-YUNUS A. CENGEL

-AFSHIN J. GHAJAR    

-유성연, 김경훈, 김병철, 김창녕, 이종붕, 조형희 공역

McGraw-Hill



문제 5-43


풀이: EES를 이용하여 주어진 연립 대수 방정식을 풀면 다음과 같다.
(a) Ctrl + N 또는 Equations Window 아이콘을 선택하여 수식을 적은 후 F2 키 또는 Solve 아이콘을 선택하여 대수 방정식을 풀 수 있다.


따라서 각각의 값은 아래와 같다.


매트랩 matlab의 solve 함수를 이용하여 주어진 연립 대수 방정식을 풀면 다음과 같다.
syms x y z
[x, y, z]=solve((3*x)-y+(3*z)==0, -x+(2*y)+z==3, (2*x)-y-z==2)


(b) Ctrl + N 또는 Equations Window 아이콘을 선택하여 수식을 적은 후 F2 키 또는 Solve 아이콘을 선택하여 대수 방정식을 풀 수 있다.


따라서 각각의 값은 다음과 같다.


매트랩 matlab의 solve 함수를 이용하여 주어진 연립 대수 방정식을 아래와 같이 풀면 다음과 같다.
syms x y z;
[x y z] = solve((4*x)-2*(y^2)+(0.5*z)==-2, (x^3)-y+z==11.964, x+y+z==3, 'Real', true, 'PrincipalValue', true)

x_1=2.532,x_2=2.364,x_3=-1.896


추가: 매트랩의 solve 함수를 이용한 풀이 중 위와 같이 추가 옵션 인수를 적용하지 않을 경우
파라미터 parameters와 조건 conditions로 다양한 해를 구할 수 있습니다.


Posted by 귀여운촌아
:

열역학 1-78.docx



열역학 Thermodynamics 5th Edition.

Fundamentals and Applications

-Yunus A. Cengel

-Michael A. Boles

-부준홍 김덕줄 김세웅 김수현 신세현 이교우 정우남 최경민 공역

McGraw-Hill



문제 1-78


공학용 소프트웨어를 이용하여 방정식의 해를 찾는다.

풀이: 가정 먼저 EES Limited Academic Version을 이용하여 방정식의 해를 구하며, 코드 내용은 아래와 같다.
((x^2)*y)-z=1
x-(3*y^0.5)+(x*z)=-2
x+y-z=2

새 Equations Window에서 위와 같이 코드 작성 후 solve 아이콘을 선택하거나 F2키를 이용하여 위의 연립 방정식을 계산한다.


따라서 EES를 이용한 연립방정식의 실근은 다음과 같다.
x=1, y=1, z=0

매트랩(matlab)의 solve() 함수를 이용한 연립방정식의 해는 다음과 같이 계산할 수 있다.
syms x y z real
[x,y,z]=solve(((x^2)*y)-z==1,x-(3*y^0.5)+(x*z)==-2,x+y-z==2)

매트랩(matlab)의 solve() 함수를 이용하면 경고 메시지와 함께 여러 해들이 출력되며, 긴 소수점을 가진 해들은 엄밀한 검증이 필요하다.


Posted by 귀여운촌아
:

열역학 1-77.docx



열역학 Thermodynamics 5th Edition.

Fundamentals and Applications

-Yunus A. Cengel

-Michael A. Boles

-부준홍 김덕줄 김세웅 김수현 신세현 이교우 정우남 최경민 공역

McGraw-Hill



문제 1-77


공학용 소프트웨어를 이용하여 방정식의 해를 찾는다.

풀이: 가정 먼저 EES Limited Academic Version을 이용하여 방정식의 해를 구하며, 코드 내용은 아래와 같다.
(2*x)-y+z=5
(3*x^2)+(2*y)=z+2
(x*y)+(2*z)=8

새 Equations Window에서 위와 같이 코드 작성 후 solve 아이콘을 선택하거나 F2키를 이용하여 위의 연립 방정식을 계산한다.


따라서 EES를 이용한 연립방정식의 실근은 다음과 같다.
x=1.141, y=0.8159, z=3.535

매트랩(matlab)의 solve() 함수를 이용한 연립방정식의 해는 다음과 같이 계산할 수 있다.
syms x y z
[x, y, z]=solve((2*x)-y+z==5,(3*x^2)+(2*y)==z+2,(x*y)+(2*z)==8,x,y,z,'Real',true)


solve() 함수에 실수근만 출력하기 위해 ‘Real’, true 옵션을 추가하여 계산하며, 출력된 결과가 복잡하므로 아래와 같이 출력되도록 한다.
syms x y z
 [x, y, z]=solve((2*x)-y+z==5,(3*x^2)+(2*y)==z+2,(x*y)+(2*z)==8,x,y,z,'Real',true);
double([x, y, z])


Posted by 귀여운촌아
:

열역학 1-76.docx



열역학 Thermodynamics 5th Edition.

Fundamentals and Applications

-Yunus A. Cengel

-Michael A. Boles

-부준홍 김덕줄 김세웅 김수현 신세현 이교우 정우남 최경민 공역

McGraw-Hill



문제 1-76


공학용 소프트웨어를 이용하여 방정식의 해를 찾는다.

풀이: 가정 먼저 EES Limited Academic Version을 이용하여 방정식의 해를 구하며, 코드 내용은 아래와 같다.
(x^3)-(y^2)=7.75
(3*x*y)+y=3.5

새 Equations Window에서 위와 같이 코드 작성 후 solve 아이콘을 선택하거나 F2키를 이용하여 위의 연립 방정식을 계산한다.


작성한 코드에 문제가 없다면 다음과 같이 계산 완료 메시지 확인 후 결과를 볼 수 있다.



EES를 이용한 연립방정식의 실근은 다음과 같다.



매트랩(matlab)의 solve() 함수를 이용한 연립방정식의 해는 다음과 같이 계산할 수 있다.
syms x y
[Sx, Sy]=solve((x^3)-(y^2)==7.75, (3*x*y)+y==3.5, x, y, 'Real', true)



solve() 함수에 실수근만 출력하기 위해 ‘Real’, true 옵션을 추가하여 계산한다.


Posted by 귀여운촌아
:

열역학 1-75.docx



열역학 Thermodynamics 5th Edition.

Fundamentals and Applications

-Yunus A. Cengel

-Michael A. Boles

-부준홍 김덕줄 김세웅 김수현 신세현 이교우 정우남 최경민 공역

McGraw-Hill



문제 1-75


공학용 소프트웨어를 이용하여 방정식의 해를 찾는다.

풀이: 가정 먼저 EES Limited Academic Version을 이용하여 방정식의 해를 구하며, 코드 내용은 아래와 같다.
(2*x^3)-(10*x^0.5)-(3*x)=-3


EES의 Calculate>Slove를 선택하거나 F2키, 또는 아래의 아이콘을 선택하여 방정식을 푼다.


아래와 같이 오류가 있을 경우


Options의 Variable Info를 선택하거나 F9키 또는 를 바로 눌러 변수 정보 창을 연다.


그 후 변수 x의 초기 추측값을 2로 수정하여 Calculate>Slove 또는 F2키, Slove 아이콘을 선택하여 방정식을 푼다.


계산된 해는 아래와 같다.


매트랩(Matlab)의 fzero() 함수를 이용하여 주어진 방정식의 해를 구한다.
fzero() 함수를 사용하기 전에 ezplot() 함수를 이용하여 방정식의 해가 정수 2 근처에 있음을 알 수 있다.
ezplot('(2*x^3)-(10*x^0.5)-(3*x)+3'), grid on

다음 fzero() 함수를 이용하여 방정식의 해를 구한다.
fzero('(2*x^3)-(10*x^0.5)-(3*x)+3',2)


아래와 같이 초기값 선택이 좋지 않으면 해를 제대로 계산되지 않을 수 있다.


매트랩(matlab)의 slove()를 이용하여 주어진 방정식의 양의 실근은 구한다.
syms s
fzero('(2*x^3)-(10*x^0.5)-(3*x)+3',[0.1 3])


slove() 함수를 이용하면 방정식의 모든 해를 계산할 수 있다.


Posted by 귀여운촌아
:


열전달 1-124.docx



열전달 HEAT AND MASS TRANSFER 4th Edition SI Units.

Fundamentals and Applications

-YUNUS A. CENGEL

-AFSHIN J. GHAJAR

-유성연, 김경훈, 김병철, 김창녕, 이종붕, 조형희 공역

McGraw-Hill



1-124


가정 : 정상상태이며 공기온도와 벽의 온도는 일정하다. 전열기는 전력을 모두 열에너지로 변환한다.

방사율과 전열기의 표면온도는 전열기의 표면에 대해 균일하고 일정한 값을 가진다.

전열기는 방의 벽에 완전히 둘러싸여있다. 기타 전도에 의한 열전달은 없는 것으로 가정한다.

 

풀이 : 전열기가 500W의 전력을 소모할 때 복사에 의한 열손실은 다음과 같다.


따라서 복사에 의한 열손실을 제외한 대류에 의한 열손실은 다음과 같다.


그러므로 대류열전달계수는 다음과 같다.


따라서 전열기가 (a) 700W를 소모하고 복사열전달을 무시하는 경우 표면온도는 대류에 의한 열손실로만 결정된다.


(b) 복사열전달을 고려하는 경우 표면온도는 다음과 같이 계산할 수 있다.


위의 연립 방정식을 EES를 이용하여 계산하면


이며, 복사 열전달에 의한 효과로 전열기의 표면온도가 대류 열전달만 고려한 경우보다 훨씬 낮은 것을 알 수 있다.

Posted by 귀여운촌아
:


열전달 1-123.docx


Q_loss.m



열전달 HEAT AND MASS TRANSFER 4th Edition SI Units.

Fundamentals and Applications

-YUNUS A. CENGEL

-AFSHIN J. GHAJAR

-유성연, 김경훈, 김병철, 김창녕, 이종붕, 조형희 공역

McGraw-Hill



1-123


가정 : 정상상태이며, 창문의 틈이나 열고 닫거나에 대한 열손실은 고려하지 않는다.


EES를 이용한 풀이 :

메뉴의 [File] [New] 또는 ctrl + N 을 이용하여 새로운 Equations Window를 연다.

Q=U*A*(Tin-Tout)

A=1.2*1.8

Tin=20

Tout=-8

위 식을 Equations Window에 적는다.

각각의 식의 구분은 줄바꿈 또는 세미콜론(;)으로 구분한다.

사용한 문자는 직접 계산하여 기입해도 되며 그림과 같이 문자에 대입하여도 된다.






인자 도표화를 사용하기 위해 메뉴의 Tables > New Parametric Table을 선택하거나 바로가기 아이콘을 선택한다.


No. of Runs를 설정하고 변수 두 개를 선택한다. 독립변수인 U-factor U와 종속변수인 열손실 Q를 각각 선택하고 Add 버튼을 클릭하여 Variables in Table로 설정한다.


독립변수 U에 적당한 값을 설정합니다.



Parametric Table은 다음과 같으며 그래프로 표현하기 위해

메뉴의 Plots > New Plot Window를 선택하거나 바로가기 아이콘을 클릭한다.



다음 Y-Axis에 종속변수 Q를 선택한다.



































Matlab을 이용한 풀이 :

크게 복잡한 문제는 아니지만 매트랩에 사용자 정의 함수를 이용해보겠습니다.

 

매트랩의 메뉴에 New Script 또는 ctrl + N 또는 New > Script (또는 Function)으로 Script Editor를 연다. 다음과 같이 스크립트를 작성한다.

function [ Q ] = Q_loss( U )

A = 1.2 * 1.8;

Tin = 20;

Tout = -8;

Q = U .* A * (Tin - Tout);

end


함수 정의가 모두 끝나고 입력값을 Command Window에서 바로 생성합니다.


U_range = [처음값 : 증가량 : 마지막값]

위와 같이 요령으로 적당한 값의 범위를 만들어 줍니다.

그 다음 미리 정의해둔 함수에 값을 전달하면 원하는 결과값을 얻을 수 있습니다.


마지막으로 plot을 이용하여 그래프로 나타냅니다.


검토

문제 1-122에서 U-factor1.25W/m2•K일 때와 6.25W/m2•K일 때 계산한 열손실 값이 같다.

Posted by 귀여운촌아
:

 

열전달 1-106.docx





열전달 HEAT AND MASS TRANSFER 4th Edition SI Units.

Fundamentals and Applications

-YUNUS A. CENGEL

-AFSHIN J. GHAJAR

-유성연, 김경훈, 김병철, 김창녕, 이종붕, 조형희 공역

McGraw-Hill




 

1-106

EES는 Engineering Equation Solver의 약자로써 컴퓨터 소프트웨어이다.

여러 공학적 문제 및 수학적 방정식을 이 소프트웨어를 이용하여 수치해적으로 풀 수 있도록 도와주는 도구이며,

EES를 제외하고 많은 공학/수학용 소프트웨어들이 있다.

 

EES를 이용한 풀이 :

 

메뉴의 [File] [New] 또는 ctrl + N 을 이용하여 새로운 Equations Window를 연다.

y=x^2

위 식을 Equations Window에 적는다.

인자 도표화(Parametric Table) 기능을 이용하기 위해서 위의 아이콘을 클릭하거나

메뉴의 Tables > New Parametric Table을 선택한다.

그 다음 나오는 팝업창에서 필요한 설정을 한다.


No. of Runs는 100, Table은 원하는 테이블명을 입력하고 Variables in equations에 있는 방정식의 변수들을 선택한 후 Add 버튼을 선택한다.


OK를 선택하여 인자 도표를 생성한다. 생성된 도표에 변수 x에 1부터 100까지 넣는다.


위와 같이 우클릭하여 Alter Values 또는 작은 역삼각형을 클릭한다.

(또는 직접 셀을 더블 클릭하여 일일이 값을 넣는다.)

위와 같이 First Value에 1, Last value를 선택하고 100을 입력한 뒤 Apply/OK를 선택한다.

x변수가 있는 column 행에 1부터 100까지 값이 지정되었다.


초록색의 화살표를 클릭하면 1부터 100까지 각각의 값을 제곱한 값을 계산할 수 있다.

다음 그래프를 그리기 위해

New Plot Window 아이콘을 클릭하거나 메뉴의 Plots > New Plot Window > X-Y Plot을 선택한다.

원하는 그래프 이름을 넣고, x 축과 y 축 변수를 클릭하여 선택한 뒤 OK를 선택한다.

Table 항목에 알맞은 Table을 선택했는지 확인한다.


Matlab을 이용한 풀이 :

Matlab을 이용하여 위의 문제를 해결할 때는 매우 다양한 방법이 있다.

하지만 지금의 풀이에서는 가장 간단하고 기본적인 방법만 사용하도록 한다.

Matlab을 실행한 후 Command Window에서 다음과 같이 명령문을 작성하고 실행한다.

위의 명령문으로 그래프까지 한번에 작성할 수 있다.

linspace() 함수를 이용하여 1부터 100까지 행렬을 생성하고 plot() 함수 안에서

power() 함수를 이용하여 각 행렬의 수를 제곱하고 이를 그래프로 그린다.


linspace() 함수는 주어진 구간에 최대한 등간격으로 주어진 개수만큼 행렬을 생성해 주는 함수이다.

linspace(x, y, n)과 같이 사용할 때 x부터 y 사이에 n개의 개수만큼 등간격 행렬을 생성해주는 함수이다.

이때 만들어지는 행렬은 x, y를 포함하며 n 값을 명시하지 않으면 100개의 원소

즉, 100개의 주어진 구간 내의 등간격 숫자를 생성한다. 또한 n=1이면 y 값을 돌려준다.

plot() 함수는 그래프를 그리는 그려주며 여러 옵션들이 있다.

power() 함수는 일반 거듭제곱이 아닌 배열(행렬)에 대한 거듭제곱으로

power(x, y)는 x의 각 요소에 y의 각 요소만큼 거듭제곱한다.

즉, x = [1 , 2, 3, 4, 5]이고, y = [9, 2, 3, 2, 4] 일 때 power(x, y)는 19, 22, 33, 42, 55

각각 계산하여 ans = [1, 4, 27, 16, 625]를 보여준다. 그러므로 x와 y는

하나의 스칼라 값을 제외하면 반드시 같은 크기의 행렬이어야 한다.

또한 power(x, y)는 x.^y와 같다.

정리하면 일반적인 거듭제곱은 a^b가 되지만 행렬의 경우 power(x, y) 또는 x.^y로 해야한다.


추가

위의 EES와 Matlab을 이용한 풀이만이 전부가 아니며 다양한 소프트웨어를 이용하여 문제를 해결할 수 있습니다.

또한 EES와 Matlab에서 여러가지 옵션을 적용할 수 있고 그 방법들은 쉽게 검색하여 알아볼 수 있습니다.

이후에 비슷한 유형의 문제가 나온다면 각각 다른 방법으로 시도해 보겠습니다.

Posted by 귀여운촌아
:


열전달 1-104.docx



열전달 HEAT AND MASS TRANSFER 4th Edition SI Units.

Fundamentals and Applications

-YUNUS A. CENGEL

-AFSHIN J. GHAJAR

-유성연, 김경훈, 김병철, 김창녕, 이종붕, 조형희 공역

McGraw-Hill



1-104

EES Engineering Equation Solver의 약자로써 컴퓨터 소프트웨어이다. 여러 공학적 문제 및 수학적 방정식을 이 소프트웨어를 이용하여 수치해적으로 풀 수 있도록 도와주는 도구이며, EES를 제외하고 많은 공학/수학용 소프트웨어들이 있다.

 

EES를 이용한 풀이 :

메뉴의 [File] [New] 또는 ctrl + N 을 이용하여 새로운 Equations Window를 연다.

(2*x)-y+z=5

(3*x^2)+(2*y)=z+2

(x*y)+(2*z)=8

위 식을 Equations Window에 적는다.

계산기 모양의 [Solve] 아이콘 또는 F2를 이용하여 근을 구한다.

연립 방정식의 경우 두 개의 식을 줄바꿈(엔터)를 이용하여 구분하거나 한 줄에 한번에 입력하는 경우에는 식과 식 사이에 세이콜론(;)으로 구분하면 된다.

x=1.141

y=0.8159

z=3.535 이다.

 

Matlab을 이용한 풀이 :

Matlab 실행 후 기본 Command Window에서 Matlab solve 함수를 이용하여 방정식의 해를 구한다. solve 함수를 이용하여 위의 방정식을 해결하기 위해서는 solve 함수에 3개의 방정식과 근을 구할 변수를 입력해야 한다. syms 키워드 또는 sym() 함수를 이용하여 3개의 변수 x, y, z를 선언한다. 그 후 solve() 함수를 이용하여 미지수를 구한다.

(연립 방정식의 해에 허수 값도 있으므로 실수로 선언한다.)

따라서 Command Window

>>syms x y z real

>> [x, y, z] = solve('(2*x)-y+z==5', '(3*x^2)+(2*y)==z+2', '(x*y)+(2*z)==8')

결과가 분수형으로 나올수 있으므로 적절히 바꾸어준다.

>> x=double(x);y=double(y);z=double(z);

결과값

x=1.140599…

y=0.815895…

z=3.534695… 이다.

 

검산

각각의 방정식에 결과값을 대입해서 계산해 본다.

 

Posted by 귀여운촌아
:


열전달 1-103.docx



열전달 HEAT AND MASS TRANSFER 4th Edition SI Units.

Fundamentals and Applications

-YUNUS A. CENGEL

-AFSHIN J. GHAJAR

-유성연, 김경훈, 김병철, 김창녕, 이종붕, 조형희 공역

McGraw-Hill



1-103

EES Engineering Equation Solver의 약자로써 컴퓨터 소프트웨어이다. 여러 공학적 문제 및 수학적 방정식을 이 소프트웨어를 이용하여 수치해적으로 풀 수 있도록 도와주는 도구이며, EES를 제외하고 많은 공학/수학용 소프트웨어들이 있다.

 

EES를 이용한 풀이 :

메뉴의 [File] [New] 또는 ctrl + N 을 이용하여 새로운 Equations Window를 연다.

위의 방정식을 그대로 적어넣는다.

계산기 모양의 [Solve] 아이콘 또는 F2를 이용하여 근을 구한다.

연립 방정식의 경우 두 개의 식을 줄바꿈(엔터)를 이용하여 구분하거나 한 줄에 한번에 입력하는 경우에는 식과 식 사이에 세이콜론(;)으로 구분하면 된다.

미지수 x, y의 값은 각각 2.215 0.6018 입니다.

 

Matlab을 이용한 풀이 :

Matlab 실행 후 기본 Command Window에서 Matlab solve 함수를 이용하여 방정식의 해를 구한다. solve 함수를 이용하여 위의 방정식을 해결하기 위해서는 solve 함수에 2개의 방정식과 근을 구할 변수를 입력해야 한다. syms 키워드 또는 sym() 함수를 이용하여 2개의 변수 x, y를 선언한다. 그 후 solve() 함수를 이용하여 미지수를 구한다.

(연립 방정식의 해에 허수도 값도 있으므로 실수로 선언한다.)

x=2.214650…

y=0.601783…

또는

으로 구할 수 있다.

 

검산

직접 대입해본다.

Posted by 귀여운촌아
: