반응형


유체역학의 핵심 지배방정식 : 나비어-스토크스 방정식 (Naview stokes equation) 소개  공학이야기 

2008/10/08 10:10

복사http://blog.naver.com/eshssy/130036019164


클레이 수학 연구소가 밀레니엄 난제로 지목한 7개 문제중 하나에 들어가는 나이어스토크스 관계식은 사실 수십년전부터 기계공학을 비롯한 각종 공학에서는 유용한 지배방정식으로 이미 활용되고 있었다. 단, 수학적 일반해를 구해서 이용한 것이 아니라 수치해석적으로 접근하여 공학적으로 이용한것에 불과하다. 이 방정식의 일반해를 구하는 것은 앞으로도 오랫동안 영구미제로 남을수도 있을 것이다. 비선형 방정식은 아주 오래전부터 지독하게도 어려운 방정식으로 정평이 나 있었기도 하다... 이미 공학을 접었으나 여전히 수학교육이라는 다분히 이공학적 마인드와 관계가 있는 분야에 종사하고 있는 관계로 과거의 악몽이었었던 나비어스토크스 관계식은 여전히 나의 가슴에 영원한 관심사로 남아있기도 하다. 수학교육과 의학의 도상에서 잠시 공학의 기억을 되살려보고자한다.

 

 

 

 

   이동현상은 크게 세가지의 전달현상을 다룬다.  운동량과 열 그리고 물질전달이 이 세가지이다.  이들의 각각 이동현상은 전달매체가 다를 뿐만 아니라 물리적인 이동과정도 매우 상이하다.  그러나 이들의 현상을 설명하는 수학적 체계는 매우 단순하리만큼 일치성을 보이고 있다는 것은 놀라운 점이다.

   여기서 내가  다루고자 하는 분야는 바로 운동량 전달이다.  이름대로 운동량의 전달현상을 목적으로 하기에 당연히 운동량의 변화를 설명할 수 있는 관계식으로부터 현상을 이해하여야 할 것이다.  에서 유체 흐름의 운동관계식을 적분형으로 다루어 볼수도 있으나, 적분형 식은 여러 현상을 적용하기에 한계가 있기에 이를 더욱 유용한 미분형으로 바꾸어 일반식을 만들어 놓는 것이 여러 가지 이유로 필요할 것이다.    따라서 여기서는 운동량 관계식을 미분형으로 유도하는 과정을 설명하고, 이에 대한 여러 관점을 함께 논하고자 한다.

 

일단 나비어 스토크스 방정식을 이해하기 위한 기본 개념을 짚고 넘어가자.

 

 

 

 

 

 

   질량보존의 법칙은 자연의 기본 원리로서 잘 알려진 법칙 중 하나이다.  즉 어떠한 물리적 과정을 거치게 될 때, 과정 전과 후의 총 질량이 변하지 않다는 원리이다.  이 법칙을 유체 흐름에 적용시키면 매우 유용하게 사용할 수 있는 관계식을 얻을 수 있게 된다.  우리가 관심을 가지는 대상부피에 대한 질량보존의 개략적 표현은 다음과 같이 나타낼 수 있다.

 

 

(1)

 

이 관계를 잘 살펴보면 매우 당연한 표현이다.  즉, 질량의 유출과 유입의 차이가 바로 대상부피 내의 질량 축적량을 의미하는 것이다.  이 관계를 수식으로 표현하기 위하여 다음의 간단한 모형을 생각해 볼 수 있다.

 

 

   이 모형은 어느 특정 시간에 형성되는 유선을 표현하였다.  dA는 대상부피 표면의 미세면적이고, v는 유선방향의 속도벡터, n는 dA의 법선벡터이다.  당연히 θ는 v와 n의 차이각이다.  따라서 dA로 통과되는 질량의 유출속도는 가 된다.  여기서 를 질량플럭스(mass flux)라고 하며 단위시간당 및 단위면적을 통과하는 질량을 의미한다.  이러한 질량의 유출속도를 벡터 표현으로 바꾸면 다음과 같이 나타낼 수 있다.

 

 

미세면적의 질량 유출속도 = 

(2)

 

이 값을 전체 대상표면에 대하여 적분을 하면 총 대상면적을 통과하는 질량의 순유출량을 알 수 있다.  이 값은 윗 식을 적분하면 얻게 된다.

 

 

(3)

 

이 값이 양수이면 질량의 순유출이 일어나며, 반대로 음수이면 순유입이 된다.  또한 0이면 대상부피 내의 질량은 항상 정량을 유지하게 된다.

 

   한편, 질량 축적속도는 결국 시간에 대한 밀도변화이므로 다음과 같이 표현할 수 있다.

 

 

(4)

 

결과적으로 질량보존의 최종 관계식은 식(3)과 (4)를 식(1)에 적용시켜서 다음의 적분형 물질수지식을 얻게 된다.

 

(5)

 

 

 

 

 

 

 

   유체흐름의 선운동량 식은 질량보존의 관계식과 매우 비슷하게 유도될 수 있다.  아래 그림에서 운동량의 유출속도는 질량유출속에 속도를 곱하는 것으로 생각할 수 있다.  자세한 관계식 유도는 질량보존의 관계식을 참고하도록 하고, 여기서는 최종 결과식만을 제시하고자 한다.

 

 

 

(1)

 

 

  1.  유체흐름의 운동량 관계식 모형

 

   유체 운동량의 적분형 관계식을 다음과 같이 유도됨을 전 절에서 다루었다.

 

 

(1)

 

   이 관계식을 풀어서 설명하면 다음과 같은 개념관계로 도시할 수 있다.

 

 

(2)

 

   이 관계식을 쉽게 설명하기 위하여 식(1)을 대상부피로 나눈 다음 극한값을 취하면 다음과 같이 쓸 수 있다.

 

 

(3)

 

이제 이 식을 하나하나의 항으로 분류하여 풀어보도록 하자.

 

 

 

 

  2.  외부 힘의 합

 

   대상부피에 작용하는 외부 힘의 영향으로는 수직응력과 전단응력, 그리고 중력과 외부로부터 받는 압력의 영향을 적용할 수 있다.  먼저 x 방향으로 작용하는 힘만을 고려한다면 다음과 같이 쓸 수 있다.

 

 

 

(4)

 

식(3)의 형식을 맞추기 위하여 식(4)를 미소부피()로 나누어 극한을 취하고, 동일한 방법으로 y 방향과 z 방향을 적용하면 다음과 같이 종합할 수 있다.

 

 

 

 

(5)

 

 

 

 

  3.  대상부피를 통과하는 선운동량의 유출속도

 

 

 

   위 그림은 대상 미소부피를 통과하는 운동량 flux를 표현한 것이다.  이 그림을 참조하여 식(2)와 (3)의 선운동량의 순 유출속도항을 플어서 쓰면 다음과 같이 표현할 수 있다.

 

 

 

                

 

 

                  

 

(6)

 

마지막 식은 연속방정식을 적용시킨 결과이다.

 

 

 

 

  4.  대상부피내의 운동량 변화 속도

 

 

(7)

 

 

 

 

  5.  종합된 결과와 Navier-Stokes 방정식

 

   지금까지 얻어진 결과들을 종합하면 다음과 같이 정리할 수 있다.

 

 

 

 

(8)

 

식(8)의 결과식을 식(3)에 대입하고 성분별로 분리하여 쓰면 다음과 같다.

 

 

(9)

 

여기서 D/Dt는 전미분을 의미한다.  식(9)는 응력-변형률 관계에 상관없이 모든 유체에 적용할 수 있으며, Stokes의 점도법칙으로부터 다음과 같이 나타낼 수 있다.

 

 

 

 

(10)

 

이렇게 유도된 식(10)을 Navier-Stokes 방정식이라고 부른다.  물론 이 식은 직교좌표의 경우이다.

 

 

 

 

  6.  ρ와 μ가 일정한 비압축성 흐름에서의 Navier-Stokes 방정식

 

   일반적인 경우는 이러한 구속조건을 가지는 경우가 많이 발생한다.  비압축성 흐름의 경우는 이므로, 이 경우에는 다음과 같이 전 좌표의 상황을 하나의 식으로 표현할 수 있다.

 

 

(11)

 

   특히 더욱 특수한 경우로서 비점성인 경우 즉, 점도가 '0'인 경우의 식을 Euler의 식이라고 불리며, 다음과 같이 나타낸다.

 

 

(12)

 

 

 

 

  7.  역학적 에너지관계식의 표현

 

   Navier-Stokes 방정식을 조금 관점을 달리 하여, 흐르는 유체상에서 에너지 관계성이 어떠한지에 대하여 알아보고자 한다.  식(9)를 벡터형식으로 통합하여 다음과 같이 나타낼 수 있다.

 

 

(13)

 

식(13)을 에너지 rate형식으로 나타내기 위하여 양 변에 지역속도를 곱하면 다음과 같이 쓸 수 있다.

 

 

(14)

 

식(14)를 각각 의미를 가지는 항으로 풀어서 정리하면 다음처럼 종합할 수 있다.

 

 

   

:  net rate of input of kinetic energy by virtue of bulk flow

(15)

: rate of increase

in kinetic energy

per unit volume

   

:  rate of work done by pressure od surroundings on volume element

   

:  rate of reversible conversion to internal energy

 

   

:  rate of work done by viscous force on volume element

 

   

:  rate of irreversible conversion to internal energy

 

   

:  rate of work done by gravity force on volume element

 

여기서 항은 Newtonian 흐름에서 항상 (+)값을 가지며, 다음과 같은 의미를 지닌다.

 

 

 

(16)

 








http://canoas27.egloos.com/2072656

 


 






반응형

+ Recent posts