1. 개요
오일러 방정식(Euler's equation)은 레온하르트 오일러에 의해 만들어진 방정식이다. 변분법과 유체역학에서 지칭하는 오일러 방정식은 서로 다른 식이므로 유의. 또한 오일러 등식과도 다르다.2. 변분법에서의 오일러 방정식
변분법에서 범함수의 최소, 최대를 찾는 방법으로 개발된 방정식이다.
2.1. 정의
변수를 하나만 가지는 함수 [math(y\left(x\right))]와 그 도함수 [math(y^\prime \left(x\right))]를 변수로 가지는 어떤 범함수(functional) [math( J )] 가 있다고 하자.| [math(\displaystyle J = \int_{ x_1 }^{ x_2 }{ f \left\{ y \left( x \right) , y^\prime \left( x \right) ; x \right\} } dx )] |
이러한 범함수의 극값(극대 또는 극소)를 주는 함수 [math( f )]가 만족시키는 미분방정식은 다음과 같다.
[math(\displaystyle \frac{ \partial f }{ \partial y } - \frac{ d }{ dx }\left( \frac{ \partial f }{ \partial y' } \right) = 0 )]
이 미분방정식을 오일러 방정식이라고 한다.
2.2. 유도
우선 일반적인 경우를 증명하기 전에 이해를 돕기 위해서 다음 예시를 먼저 살펴보자.두 점을 잇는 가장 짧은 곡선은 직선임을 보여라.
먼저, 두 점을 잇는 선을 어떤 현악기의 줄이라고 생각해 보자. 그 줄의 길이를 [math( L )]이라 하자.
그 다음 줄을 손으로 1mm 정도 잡아당겨 보자. 그럼 아래 그림처럼 줄이 휘어질 것이다. (그림은 과장해서 그렸다)
이렇게 휘어진 줄과 원래 직선과의 차이를 [math( \eta (x) )]라고 하자. (줄을 함수의 그래프처럼 생각하면 된다. x축이 바이올린의 줄과 평행한 좌표계에서 처음 상태를 [math( f(x) )]라고 하고 휘어진 그래프를 [math( g(x) )]라고 하면 [math( \eta (x) = g(x) -f(x) )]인 것이다.) 그런데, 줄의 끝은 고정되어 있다. 즉 차이가 없는 것이다. 따라서 [math( \eta (x_1) = \eta(x_2) = 0 )]이다. (줄의 끝의 좌표가 [math( x_1, x_2 )]이다)
그럼 여기서 줄을 더 잡아당겨서 5mm 정도 잡아당기면 어떻게 될까? 차이가 5배로 커질 것이다. 즉 차이는 [math( 5 \eta (x) )]가 된다.
마찬가지로, 2mm를 잡아당기면 차이는 대략 [math( 2 \eta (x) )]이고, 10mm를 잡아당기면 [math( 10 \eta (x) )]이고, 아니면 반대 방향으로 3mm를 잡아당기면 차이는 [math( -3 \eta (x) )]가 될 것이다. 즉 [math( \alpha )] mm 잡아당기면 차이는 [math( \alpha \eta (x) )]가 된다.
이때, 줄을 함수의 그래프 [math( y = f(x) )]라고 생각해 보자. (이제부터 간단하게 [math( y(x) )]로 쓰자.) 먼저, 악기를 가만히 놔두면 당연히 줄은 직선(최단거리) 모양이 된다. 이것을 [math( y(0,x) )]라고 하자. 그리고 [math( \alpha )] mm만큼 줄을 잡아당겼을 때의 모양을 아래 그림처럼 [math( y( \alpha, x))]라고 하자.
즉 줄을 5mm 잡아당긴 함수는 [math( y(5,x) )]이고, 반대쪽으로 3mm 잡아당기면 [math( y(-3,x) )]이다. 그런데 아까 줄을 [math( \alpha )] mm 잡아당기면 원래 상태와의 차이는 [math( \alpha \eta (x) )]가 된다고 했다. 따라서 다음과 같이 쓸 수 있다.
|
[math( \displaystyle y(\alpha ,x)=y(0,x)+\alpha \eta (x) )] (잡아당긴 줄은, 원래 줄에다가 잡아당김을 더한 것이다.) |
즉, 그래프 [math( y )]는 이제 [math( \alpha )]에 따라 변하는 함수가 된 것이다. 예를 들어, 줄을 튕겨서 진동시키는 것은 [math( \alpha )]를 진동시키는 것이라고 생각하면 된다. 그러면 당연히 줄의 길이 [math( L )]도 [math( \alpha )]에 대한 함수가 된다! 예를 들어, 원래 길이가 50cm였는데 줄을 잡아당기면 51cm가 될 수도 있는 것이다. 그러면 줄의 길이 [math( L )]를 잡아당긴 거리 [math( \alpha )]에 따라 그래프를 그릴 수 있다.
그런데 줄이 가장 짧을 때는 언제일까? 당연히 하나도 안 잡아당겼을 때, 즉 [math( \alpha=0 )]일 때이다. 그리고 미적분을 배우면 알겠지만 [math( \alpha = 0 )]에서 줄의 길이가 극소가 된다는 것은 미분계수가 0이라는 것을 의미한다.
| \displaystyle \left[ {\frac{ \partial L }{ \partial \alpha }} \right] _ {\alpha = 0} = 0 </math> |
결국 위 식을 풀면 최단거리는 직선일 때임을 증명할 수 있다. 이제 계산을 해 보자.
두 점 사이의 곡선의 길이 공식은 다음과 같다.
| \displaystyle L = \int_{x_1}^{x_2}{\sqrt{1+{y'}^2}dx} </math> |
이때 [math( f = \sqrt{1+{y'}^2} )]라고 하면 [math( L = \int_{x_1}^{x_2}{f dx} )]라고 쓸 수 있다. 적분의 위끝 아래끝이 상수이므로 적분 기호 안에서 미분해도 된다. 이것을 [math( \alpha )]로 미분하면
| \displaystyle \frac{ \partial L }{ \partial \alpha } = \frac{ \partial }{ \partial \alpha } \int_{ x_1 }^{ x_2 }{ f dx } = \int_{ x_1 }^{ x_2 }{ \frac{ \partial f }{ \partial \alpha} dx } </math> |
이 되는데, 그래프가 직선이라면 위 식은 0이 되어야 한다. 여기서 연쇄 법칙(합성함수의 미분법)을 쓰면, [math( \displaystyle \frac{ \partial f }{ \partial \alpha} = \frac{ \partial f }{ \partial y'} \frac{ \partial y' }{ \partial \alpha} )]이다. 이를 계산하면
| \displaystyle \frac{ \partial L }{ \partial \alpha } = \int_{ x_1 }^{ x_2 }{ \frac{ \partial f }{ \partial y' } \frac{ \partial y' }{ \partial \alpha } dx} = \int_{ x_1 }^{ x_2 }{ \frac{ y' }{ \sqrt {1+{y'}^2} } \frac{ \partial y' }{ \partial \alpha } dx} = 0 </math> |
한편, 위에서 [math( y(\alpha ,x)=y(0,x) + \alpha \eta (x) )]라는 관계식을 보자. 이를 [math( x )]로 미분해서 도함수를 구하면 [math( y'(\alpha ,x) = y'(0,x) + \alpha \eta ' (x) )]이다. 그런데 이것을 다시 [math( \alpha )]로 편미분하면 [math( \displaystyle \frac{ \partial y' }{ \partial \alpha } = \eta ' (x) )]이다. ([math( \alpha )]를 제외하고는 다 상수로 보고 미분한다.) 이를 대입하면
| \displaystyle \frac{ \partial L }{ \partial \alpha } = \int_{ x_1 }^{ x_2 }{ \frac{ y' }{ \sqrt {1+{y'}^2} } \eta ' (x) dx} = 0 </math> |
그런데 이 적분은 부분적분이 되는 꼴이다. [math( \eta ' (x) )]를 적분, [math( \frac{ y' }{ \sqrt {1+{y'}^2}} )]를 미분할 함수로 두면 된다. 이를 계산하면 다음 식을 얻는다.
| \displaystyle \frac{ \partial L }{ \partial \alpha } = \left[ \frac{ y' }{ \sqrt {1+{y'}^2}} \eta (x) \right]_{x_1}^{x_2} - \int_{ x_1 }^{ x_2}{ \frac{d}{dx} \left( \frac{ y' }{ \sqrt {1+{y'}^2}} \right) \eta (x) dx } = 0</math> |
그런데 아까 줄의 끝이 고정되어 있다고 했다. 즉 [math( \eta (x_1) = \eta (x_2) = 0 )]이므로, 첫 항은 사라진다. 그러면 뒤에 있는 적분이 남는다.
| \displaystyle \int_{ x_1 }^{ x_2}{ \left[ \frac{d}{dx} \left( \frac{ y' }{ \sqrt {1+{y'}^2}} \right) \right] \eta (x) dx } = 0</math> |
여기서 잘 생각해 보자. [math( \eta (x) )]는 잡아당긴 줄과 원래 줄의 차이었다. 그런데 줄을 한 손가락으로 당길 수도 있는 것이고, 두 손가락으로 당길 수도 있는 것이고, 활로 그을 수도 있고, 줄 맨 끝을 튕길 수도 있고, 정중앙에 손가락을 대고 튕겨서 하모닉스 소리를 낼 수도 있는 것이다. 이게 무슨 말이냐면, [math( \eta (x) )]가 도대체 어떻게 생겼는지 알 수가 없다는 것이다. 즉 어떤 모양의 [math( \eta (x) )]를 가져와도 무조건 저 식이 0이 되도록 해야 한다. 항등식의 성질을 생각해 보면, 그렇게 하려면 대괄호 안이 0이 되는 수밖에 없다는 것을 알아챌 수 있을 것이다.[1] 따라서
| [math( \displaystyle \frac{d}{dx} \left( \frac{ y' }{ \sqrt {1+{y'}^2}} \right) = 0 )] |
이어야 한다. 양변을 [math( x )]로 적분하면
| \displaystyle \frac{ y' }{ \sqrt {1+{y'}^2}} = C</math> |
이고, [math( C )]는 적분상수이다. 이를 [math( y' )]에 대하여 정리하면
| \displaystyle y' = \sqrt{ \frac{ C^2 }{ 1 - C^2 }} = \text{상수} = a</math> |
따라서 기울기는 [math( a )]로 일정한 직선이다. 즉 [math( y = ax+b )]를 얻는다!
여기까지가 오일러 방정식을 유도하는 과정을 예로 들어 본 것이다. 이제 일반적인 범함수에서 오일러 방정식을 유도해 보자.
2.2.1. 오일러 미분방정식
한 지점 [math( (x_1,y_1) )]에서 다른 지점 [math( (x_2,y_2) )]까지 적분으로 정의된 범함수 [math( J = \int_{ x_1 }^{ x_2 }{ f \left\{ y(x), y'(x); x \right\} }dx )]를 생각해 보자. 우리의 목표는 이러한 [math( J )]가 극값(극대 또는 극소)가 되는 [math( y(x) )]를 찾는 것이고, 이때의 [math( y(x) )]를 [math( y(0,x) )]라고 하자. 그러면 가능한 모든 [math( y(x) )]를 아래 그림과 같이 다음과 같은 꼴로 쓸 수 있다.| \displaystyle y(\alpha ,x)=y(0,x)+\alpha \eta (x) </math> |
이때 [math( \eta (x) )]는 미분가능하고 [math( \eta (x_1) = \eta (x_2) = 0 )]을 만족하는 임의의 함수이고, [math( \alpha )]는 임의의 작은 실수이다. 위 식의 의미는 [math( y(x) )]가 정답으로부터 [math( \alpha \eta (x) )]만큼 살짝 벗어나 있다는 것이다. 이제 [math( J )]는 다음과 같이 [math( \alpha )]에 대한 함수가 되었다.
| \displaystyle J(\alpha) = \int_{ x_1 }^{ x_2}{ f \left\{ y(\alpha ,x), y'(\alpha ,x); x \right\} }dx </math> |
그런데 [math( \alpha = 0 )]를 대입하면 [math( y(\alpha ,x) )]는 [math( y(0,x) )], 즉 극값이 된다. [math( \alpha = 0 )]에서 극값을 가진다는 것을 수학적으로 표현하면 다음과 같다.
| \displaystyle \left[ {\frac{ \partial J }{ \partial \alpha }} \right] _ {\alpha = 0} = 0 </math> |
(위 식은 간단하게 [math( \delta J = 0 )]로 쓰기도 한다. 변분 참고)
이제 편미분을 계산해 보자. [math( J )]를 [math( \alpha )]로 편미분하면 다음과 같다.
| \displaystyle \frac{ \partial J }{ \partial \alpha } = \frac{ \partial }{ \partial \alpha } \int_{ x_1 }^{ x_2 }{ f \left\{ y, y'; x \right\} }dx </math> |
그런데 적분 기호의 위끝과 아래끝은 상수이므로, 적분 기호 안에서 미분할 수 있다. 다변수 함수 [math( f )]에 연쇄 법칙을 이용하면 위 식은 다음과 같이 쓸 수 있다.
| \displaystyle \frac{ \partial J }{ \partial \alpha } = \int_{ x_1 }^{ x_2 }{ \left\{ \frac{ \partial f }{ \partial y } \frac{ \partial y }{ \partial \alpha } + \frac{ \partial f }{ \partial y' } \frac{ \partial y' }{ \partial \alpha } + \frac{ \partial f }{ \partial x } \frac{ \partial x }{ \partial \alpha }\right\} }dx </math> |
그런데 [math( y(\alpha ,x)=y(0,x)+\alpha \eta (x) )]이고, 양변을 [math( x )]로 미분하면 [math( y' (\alpha ,x)=y'(0,x) + \alpha \eta ' (x) )]가 된다. 그리고 당연히 [math( \frac{\partial}{\partial \alpha} x = 0)]이다. [math( \alpha )]로 편미분하는 것은 [math( x )]를 상수 취급 하기 때문이다. 하지만 [math( y )] 및 [math( y' )]는 [math( \alpha )]에 대한 함수라는 것에 주목하자.[2] 따라서 각 변수를 [math( \alpha )]로 편미분하면, 다음을 얻는다.
| \displaystyle \frac{ \partial y }{ \partial \alpha } = \eta (x), \frac{ \partial y' }{ \partial \alpha } = \eta ' (x), \frac{\partial x}{\partial \alpha} = 0 </math> |
이를 위 식에 대입하면,
| \displaystyle \frac{ \partial J }{ \partial \alpha } = \int_{ x_1 }^{ x_2 }{ \left\{ \frac{ \partial f }{ \partial y } \eta \left( x \right) + \frac{ \partial f }{ \partial y' } \eta ' (x) \right\} }dx </math> |
여기서 두 번째 항에 [math( \eta ' (x) )]를 적분, [math( \frac{ \partial f }{ \partial y' } )]를 미분할 함수로 두고 부분적분법을 적용하면 다음과 같다.
| \displaystyle \frac{ \partial J }{ \partial \alpha } = \int_{ x_1 }^{ x_2 }{ \frac{\partial f}{\partial y} \eta \left( x \right) dx } + \left[ \frac{\partial f}{\partial y'} \eta \left( x \right) \right]_{x_1}^{x_2} - \int_{ x_1 }^{ x_2}{ \frac{d}{dx} \left( \frac{\partial f}{\partial y'} \right) \eta \left( x \right) dx } </math> |
위에서[math( \eta \left( x \right) )]를 정의할 때 [math( \eta (x_1) = \eta (x_2) = 0 )]였으므로 가운데 항은 날아간다. 이제 [math( \eta (x) )]로 묶어서 정리하면
| \displaystyle \frac{ \partial J }{ \partial \alpha } = \int_{ x_1 }^{ x_2 }{ \left[ \frac{\partial f}{\partial y} - \frac{d}{dx} \frac{\partial f}{\partial y'} \right] \eta \left( x \right) dx } </math> |
이 되고, 이것이 0이 되어야 한다. 그런데 이 식은 임의의 [math( \eta (x) )]에 대하여 항상 0이 되어야 하므로, 결국 대괄호 안 전체가 0이 되는 방법밖에 없다. (사실 엄밀하게 이것은 변분법의 기본정리에 의한 것으로, 자세한 증명은 항목에 있다.) 따라서 다음의 오일러 방정식을 얻는다.
| \displaystyle \frac{ \partial f }{ \partial y } - \frac{ d }{ dx }\left( \frac{ \partial f }{ \partial y' } \right) = 0 </math> |
2.2.2. 예제
위에서 했던 평면에서 두 점 [math( (x_1,y_1) )]과 [math( (x_2,y_2) )] 사이의 최단경로는 직선임을 증명해 보자.
두 점 사이를 잇는 곡선의 길이 공식은 다음과 같다.
| \displaystyle L = \int_{x_1}^{x_2}{\sqrt{1+{y'}^2}dx} </math> |
이때 [math( f = \sqrt{1+{y'}^2} )]이고, 목표는 [math( L )]를 최소화시키는 것이다.
오일러 방정식 [math( \frac{ \partial f }{ \partial y } - \frac{ d }{ dx }\left( \frac{ \partial f }{ \partial y' } \right) = 0 )]에 [math( f = \sqrt{1+{y'}^2} )]를 대입하면 바로 풀린다. 여기서 [math( f )]는 [math( y' )]만의 함수이므로 당연히 [math( \frac{ \partial f }{ \partial y } = 0 )]이고, 오른쪽 항만 계산하면
| \displaystyle \frac{d}{dx} \left( \frac{ y' }{ \sqrt {1+{y'}^2}} \right) = 0 </math> |
이다. 여기서부터는 위랑 똑같이 하면 된다. 미분 안은 상수이고, 따라서 [math( y' )]는 상수이다. 기울기가 상수인 함수는 직선이다.
2.3. 벨트라미 항등식
오일러 방정식을 풀 때, 위의 예제처럼 [math( f )]가 [math( y', x )]만의 함수이고, [math( y )]와는 독립일 경우 첫째 항 [math( \displaystyle \frac{\partial f}{\partial y} )]가 사라지기 때문에 풀기 쉽다. 한편, [math( f )]가 [math( y, y')]만의 함수이고 [math( x )]는 들어가지 않을 때, 즉 [math( f(y,y') )]일 때도 쉽게 변형해서 푸는 방법이 있는데, 이를 벨트라미 항등식(Beltrami identity)라고 한다. [math( \displaystyle \frac{\partial f}{\partial x} = 0 )]일 때, 오일러 방정식은 다음 방정식과 동치이다.[math( \displaystyle f - y' \frac{\partial f}{\partial y' } = \mathsf{constant})]
(벨트라미 항등식)
(벨트라미 항등식)
참고로, 좌변은 [math( f )]를 [math( y' )]에 대해 르장드르 변환한 것이다.
- [ 증명 보기 · 숨기기 ]
[math( f )]는 [math( y,y' )]의 함수이다. 따라서 연쇄 법칙에 의해
\displaystyle \frac{df}{dx} = \frac{\partial f}{\partial y} \frac{dy}{dx} + \frac{\partial f}{\partial y'} \frac{dy'}{dx} = y' \frac{\partial f}{\partial y} + y'' \frac{\partial f}{\partial y'} </math>
한편, [math( \displaystyle y' \frac{\partial f}{\partial y'} )]를 [math( x )]로 미분하면 곱의 미분법에 의해
\displaystyle \frac{d}{dx} \left( y' \frac{\partial f}{\partial y'} \right) = y'' \frac{\partial f}{\partial y'} + y' \frac{d}{dx} \frac{\partial f}{\partial y'} </math>
위 두 식에서 [math(y'' \frac{\partial f}{\partial y'})]를 소거하면
\displaystyle \frac{d}{dx} \left( y' \frac{\partial f}{\partial y'} \right) = \frac{df}{dx} - y' \frac{\partial f}{\partial y} + y' \frac{d}{dx} \frac{\partial f}{\partial y'} </math>
여기서 마지막 두 항은 [math( y' )]으로 묶을 수 있다. 정리하면
\displaystyle \frac{d}{dx} \left( f - y' \frac{\partial f}{\partial y'} \right) - y' \left( \frac{\partial f}{\partial y} - \frac{d}{dx} \frac{\partial f}{\partial y'} \right) = 0 </math>
그런데 마지막 괄호 안에 있는 식은 오일러 방정식이랑 똑같은 모양이다! 따라서 괄호 안은 0이 되어서 사라진다. 남은 항을 [math(x)]로 적분하면 증명이 끝난다.
[math( \displaystyle f - y' \frac{\partial f}{\partial y' } = \mathsf{constant})]
참고로 이 식의 [math(f)]에다가 라그랑지언 [math(\mathscr L)]을 대입하면 이는 해밀토니언 [math(\mathcal H)]의 정의가 된다! 따라서 이는 해밀토니언이 보존된다는 결과를 의미한다.
2.3.1. 예제
[math( (x_1,y_1) )]과 [math( (x_2,y_2) )]를 이은 곡선을 [math( x )]축을 중심으로 회전한 회전체의 겉넓이가 최소가 되는 경로를 찾아보자. 회전체의 겉넓이 [math( S = \int_{x_1}^{x_2}{2 \pi y \sqrt{1 + {y'}^2} dx} )]이다. 그런데 적분변수를 [math( dy )]로 치환하면| \displaystyle S = \int_{x_1}^{x_2}{2 \pi y \sqrt{1 + {y'}^2} dx} = \int_{y_1}^{y_2}{2 \pi y \sqrt{1 + \frac{1}{{x'}^2}}x'dy} = \int_{y_1}^{y_2}{2 \pi y \sqrt{1 + {x'}^2}dy} </math> |
| \displaystyle \frac{ d }{ dy }\left( y \frac{x'}{\sqrt{1 + {x'}^2}} \right) = 0 </math> |
이다. 양변을 [math( y )]로 적분하면
| \displaystyle \left( y \frac{x'}{\sqrt{1 + {x'}^2}} \right) = a </math> |
이고, [math( a )]는 적분상수이다. 이를 [math( x' )]에 대해서 풀면
| \displaystyle x' = {dx \over dy} = \frac{a}{\sqrt{y^2 - a^2}} </math> |
| \displaystyle x = a \cosh^{-1}{{y \over a}} + b </math> |
| \displaystyle y = a \cosh \frac{x-b}{a} </math> |
2.4. 응용
물리적으로는 라그랑주 역학의 핵심인 최소작용의 원리를 수학적으로 표현한 것이다. 대표적으로 [math(x)]가 시간, [math(y)]가 거리, [math(f)]가 라그랑지언이라고 하면 라그랑지안 역학의 기본꼴이 되며, 에너지의 표현식을 알면 힘의 작용점 분석 없이도 운동방정식을 얻을 수 있다.3. 유체역학에서의 오일러 방정식
|
유체역학 Fluid Mechanics |
||
|
{{{#!wiki style="margin:0 -10px -5px; min-height:2em; word-break:keep-all" {{{#!folding [ 펼치기 · 접기 ] {{{#!wiki style="margin:-6px -1px -11px" |
<colbgcolor=#0D98BA><colcolor=#fff> 유체와 힘 | <colbgcolor=#fff,#1f2023> 유체 · 뉴턴 유체 · 비뉴턴 유체(멱법칙 유체 · 오스트발트-드 웰 관계식 · 허쉘-버클리 유체 · 리-아이링 이론) · 압력 · 부력 ( 아르키메데스의 원리) · 항력 ( 수직항력 · 스토크스 법칙) · 응력( 응용) · 양력 · 표면장력 · 상 · 밀도 · 기체 법칙 ( 이상 기체 법칙) · 달랑베르의 역설 |
| 유체동역학 | 유동 ( 압축성 · 탄성 · 점성/ 점성계수) · 난류 및 층류 · 레이놀즈 수송 정리 (체적 검사) | |
| 무차원수 | 마하 수 · 레이놀즈 수 · 프란틀 수 · 레일리 수 · 그라스호프 수 · 슈미트 수 · 네버러 수 · 프루드 수 | |
| 방정식 | 나비에-스토크스 방정식 · 연속 방정식 · 오일러 방정식 · 구성 방정식 · 베르누이 방정식 · 파스칼의 원리 · 브라운 운동 방정식 · 하겐-푸아죄유 법칙 · 글래드스톤-데일 방정식 | |
| 응용 및 현상 | 날씨 · 모세관 현상 · 마그누스 효과 · 케이 효과 · 카르만 효과 · 사이펀의 원리 · 대류 현상 · 슬립 스트림 · 최대동압점 · 스탈링 방정식 · 벤추리 효과 · 레인-엠든 방정식 · 엠든-찬드라세카르 방정식 · 라이덴프로스트 효과 | |
| 유체역학 연구 | 전산유체역학( CFD) · 풍동 실험 · 차원분석 | }}}}}}}}} |
유체역학에서 비점성 유체, 즉 전단 응력(shearing stress)이 0인 유체의 운동량의 변화를 묘사하는 물질 (시간) 도함수(material derivative 또는 material time derivative)에 관한 미분 방정식이다.
| \displaystyle \rho \frac{D \mathbf{u}}{Dt}=-\nabla P+\rho \mathbf{g}</math> |
3.1. 유도 1
밀도가 [math(\displaystyle \rho)]인 유체가 중력장 [math(\displaystyle \mathbf{g})] 내에서 유체 사이의 압력이 [math(\displaystyle P)]로 주어진 상황을 생각해보자. 어떤 폐곡면 [math(\mathbf{S})]로 둘러싸인 영역 [math(\mathbf{V})] 내에 있는 유체가 갖는 운동량의 변화는| \displaystyle \frac{\partial}{\partial t} \int_\mathbf{V} \rho\mathbf{u} d^3 r = - \oint_\mathbf{S} \mathbf{u}(\rho\mathbf{u} \cdot d\mathbf{a}) - \oint_\mathbf{S} P d\mathbf{a} + \int_\mathbf{V}\rho\mathbf{g} d^3 r</math> |
와 같은 방정식으로 표현된다.
여기서 좌변은 영역 [math(\mathbf{V})] 내에 있는 유체의 단위 시간당 운동량의 변화량을 뜻하고 우변의 첫째항, 둘째항, 셋째항은 각각 단위 시간당 경계면 [math(\mathbf{S})]를 통해 유입되는 유체가 갖는 운동량, 압력에 의해 영역 [math(\mathbf{V})] 내의 유체가 받는 힘, 영역 [math(\mathbf{V})] 내에 있는 유체가 받는 중력을 의미한다.
위 방정식에 발산 정리를 적용하면
| \displaystyle \frac{\partial}{\partial t}(\rho\mathbf{u}) = -\mathbf{e}_{i} \nabla \cdot [u_{i} (\rho \mathbf{u})] - \nabla P + \rho \mathbf{g}</math> |
와 같이 정리된다.
이 방정식의 우변의 첫째항을 이항하여 정리하면 다음과 같이 오일러 방정식(Euler equation)이 유도된다. (참고: 질량에 대한 연속 방정식)
| \displaystyle \rho \frac{D \mathbf{u}}{Dt}=-\nabla P+\rho \mathbf{g}</math> |
여기서 [math(\displaystyle\frac{D}{Dt} \equiv \frac{\partial}{\partial t}+(\mathbf{u}\cdot \nabla)=\frac{\partial}{\partial t}+u_{i} \frac{\partial}{\partial x_{i}})]는 물질 도함수(material derivative)라 불린다.
homentropic flow(단위 질량당 엔트로피가 공간에 걸쳐 균일)로 가정하면 [math(w)]를 단위 질량당 엔탈피(enthalpy)라 할 때 열역학에 따르면 [math(\displaystyle \frac{dP}{\rho} = dw)]라 표현할 수 있는 것이 알려져 있다. 정리하면 다음과 같은 결과를 얻는다.
| \displaystyle \frac{ \partial \mathbf{u} }{ \partial t } + \left(\mathbf{u} \cdot \nabla\right) \mathbf{u} = -\nabla w + \mathbf{g} </math>[3] |
한편, 이 식이 다루는 대상이 '비점성'이라고 했는데, 이제 여기 좌변에다가 점성항 [math(- \nu \nabla^2 \mathbf{u})]를 얹어주면 바로 나비에-스톡스 방정식이 된다.
3.2. 유도 2
길이, 너비, 높이가 각각 [math(\text{d}x,\text{d}y,\text{d}z)]인 유체의 미소부피를 고려해보자. 점성과 외력이 없다면 이 미소부피가 가속하는 이유는 각 면들에 가해지는 압력들의 차이 때문이다. 예를들어, [math(x)]방향으로 작용하는 압력 차이는 [math(\displaystyle \frac{\partial P}{\partial x} \text{d}x)]이다.
힘 = 압력 x 면적이므로 이 압력 차이 때문에 나타나는 힘은 [math(\displaystyle \frac{\partial P}{\partial x} \text{d}x\text{d}y\text{d}z=\frac{\partial P}{\partial x} \text{d}V)]이다.
이게 양수라면 이 미소부피는 음의 [math(x)]방향으로 가속하므로 [math(\displaystyle \text{d}F_x = -\frac{\partial P}{\partial x} \text{d}V)]가 성립한다.
| \displaystyle \text{d} \textbf{F} = -\nabla P \, \text{d}V</math> |
이제 뉴턴의 제 2법칙을 적용하자. [math(\displaystyle \textbf{F} =m \textbf{a} \Rightarrow -\nabla P \, \text{d}V=\text{d}m \, \textbf{a}=\rho \, \text{d}V \textbf{a} \Rightarrow -\nabla P=\rho \, \textbf{a}\Rightarrow -\nabla P=\rho \, \frac{\text{D} \textbf{u}}{\text{D}t})].
여기서 [math(\displaystyle \frac{\text{D} \textbf{u}}{\text{D}t})]는 "material derivative" 또는 "total derivative"라고 하는 항이다. 단, [math(\displaystyle \textbf{u})]는 유체의 속도를 나타내는, [math(\displaystyle x,y,z,t)]의 벡터함수이다.
[math(\displaystyle \frac{\text{D} \textbf{u}}{\text{D}t}=\frac{\partial \textbf{u}}{\partial t}+\frac{\partial \textbf{u}}{\partial x} \, \frac{\text{d}x}{\text{d}t}+\frac{\partial \textbf{u}}{\partial y} \, \frac{\text{d}y}{\text{d}t}+\frac{\partial \textbf{u}}{\partial z} \, \frac{\text{d}z}{\text{d}t}=\frac{\partial \textbf{u}}{\partial t}+\Bigg(\frac{\partial}{\partial x} \, u_x+\frac{\partial}{\partial y} \, u_y+\frac{\partial}{\partial z} \, u_z\Bigg)\textbf{u}=\frac{\partial \textbf{u}}{\partial t}+(\textbf{u}\cdot \nabla)\textbf{u})]
이제 양쪽을 [math(\displaystyle \rho)]로 나누고 중력 같은 외력에 의한 가속도 [math(\displaystyle \textbf{g})]를 더해주면 된다. [math(\displaystyle \frac{\nabla P}{\rho}=\nabla w)]를 사용하면 익숙한 오일러 방정식 완성.
| \displaystyle \frac{\partial \textbf{u}}{\partial t}+(\textbf{u}\cdot \nabla)\textbf{u}=-\nabla w+\textbf{g}</math> |
보너스로, 점성을 고려하고 싶다면 우항에 [math(\nu \nabla^2 \textbf{u})]를 더해줘야 한다. [math(\nabla^2)]은 라플라시안 연산자며, [math(\nu)]는 뉴턴의 점성법칙에 나오는 점성상수[4].
점성으로 인한 가속도가 왜 [math(\nu \nabla^2 \textbf{u})]일까? 라플라시안 항목에서도 서술되어 있듯이, 여기서 [math(\nabla^2 \textbf{u})]의 물리적 의미는 어떤 미소부피의 속도와 이 미소부피 근방의 평균적인 유체 속도의 차이다. 이 속도 차이가 클수록 이 미소부피는 점성에 의해 받는 힘이 늘어나며, 이 미소부피의 속도가 주위 유체들로 (또는 주위 유체의 속도가 미소부피로) '전달'된다[5].
이렇게 [math(\nu \nabla^2 \textbf{u})]를 우변에 더해주고 좌변으로 옮기면 나오는게 그 유명한 (비압축성) 나비에-스토크스 방정식. 더 자세한 정보는 나비에-스토크스 방정식 문서에.
| \displaystyle \frac{\partial \textbf{u}}{\partial t}+(\textbf{u}\cdot \nabla)\textbf{u}-\nu \nabla^2 \textbf{u}=-\nabla w+\textbf{g}</math> |
3.3. 의미
뉴턴법칙 중 제2법칙을 유체의 운동에 적용한 것에 지나지 않는다.이 방정식은 (위치에 따른) 압력의 차이와 유체에 가해지는 중력이 유체의 운동량의 (시간에 따른) 변화를 유도함을 표현한다.
3.4. 운동량에 대한 연속 방정식
오일러 방정식을 텐서를 이용하여 정리한 방정식이다. 여기서 [math( \displaystyle \mathbf{\Pi} )]는 운동량 선속 텐서(momentum flux tensor)라 불리는 2차 텐서이다.| [math( \displaystyle \frac{\partial}{\partial t}(\rho \mathbf{u}) + \nabla \cdot \mathbf{\Pi} = \rho \mathbf{g} )] |
3.4.1. 유도
오일러 방정식에서[math(\displaystyle \rho \left [ \frac{\partial \mathbf{u}}{\partial t} + \left ( \mathbf{u} \cdot \nabla \right ) \mathbf{u} \right ] = \frac{\partial}{\partial t} \left ( \rho \mathbf{u} \right ) - \mathbf{u} \frac{\partial \rho}{\partial t} + \left ( \rho \mathbf{u} \cdot \nabla \right ) \mathbf{u} = \frac{\partial}{\partial t} \left ( \rho \mathbf{u} \right ) + \mathbf{u} \left [ \nabla \cdot \left ( \rho \mathbf{u} \right ) \right ] + \left ( \rho \mathbf{u} \cdot \nabla \right ) \mathbf{u} = \frac{\partial}{\partial t} \left ( \rho \mathbf{u} \right ) + \mathbf{e}_{i} \nabla \cdot \left ( \rho u_{i} \mathbf{u} \right ) )]
[math(\displaystyle \nabla P = \mathbf{e}_{i} \frac{\partial P}{\partial x_{i}} = \mathbf{e}_{i} \frac{\partial}{\partial x_{j}}(P \delta_{ij}) = \mathbf{e}_{i} \nabla \cdot\left [ \mathbf{e}_{j}(P \delta_{ij}) \right ] )]
라 할 수 있으므로 이를 대입하여 정리하면 다음과 같이 정리된다. (참고: 질량에 대한 연속 방정식)
| \displaystyle \frac{\partial}{\partial t} (\rho \mathbf{u}) + \mathbf{e}_{i} \nabla \cdot\left [ \mathbf{e}_{j}(P \delta_{ij} + \rho u_{i} u_{j}) \right ] = \rho \mathbf{g} </math> |
[math(\displaystyle \Pi_{ij} \equiv P \delta_{ij} + \rho u_{i} u_{j} )]라 정의하면
| [math( \displaystyle \frac{\partial}{\partial t}(\rho \mathbf{u}) + \nabla \cdot \mathbf{\Pi} = \rho \mathbf{g} )] |
과 같이 운동량에 대한 연속 방정식을 얻는다.
4. 관련 문서
[1]
이것의 엄밀한 증명은
변분법의 기본정리를 참고하라.
[2]
이것은 상당히 중요한 사실이다. 지금 [math( \alpha )]의 변화를 보고 있는데, [math( x )]는 [math( \alpha )]와 독립이고 [math( y, y' )]만 [math( \alpha )]의 함수이다. 이것은 극값을 따질 때 [math( y,y' )]의 변화만 보겠다는 것이다. 이를 통해
라그랑주 역학에서
액션이 시간의 변화와 관계 없다는 사실을 알 수 있다.
[3]
표기법은 나비에-스톡스 방정식 항목의 식과 동일하다
[4]
점성이 일정한 유체를 '뉴턴 유체'라고 부른다. 비뉴턴 유체에는 이 형태의 나비에 스톡스 방정식이 안 먹힌다.
[5]
뜨거운 곳에서 차가운 곳으로 흐르는 열 전달과 매우 비슷하다. 사실
열 방정식도 라플라시안으로 서술된다. 여기선 이 점성항을 직관적이게 설명했는데, 좀 더 수학적으로 엄밀하게 구하려면
텐서 개념이 필요하다.