19.8 오일러 방정식과 이상 유체 운동

1. 개요

오일러 방정식은 비점성 유체의 운동을 기술하는 운동량 보존 법칙의 미분 형태이며, 유체역학의 고전적 방정식 중 하나이다. 오일러가 1757년에 제시한 이 방정식은 점성을 무시할 수 있는 이상 유체 운동의 해석에서 중심적 역할을 하며, 고전 공기역학의 여러 결과들의 기반이 된다. 본 절은 오일러 방정식의 유도, 벡터 형태와 성분 형태, 경계 조건, 비회전 유동과 포텐셜 이론의 연결, 와도와 순환의 개념, Kelvin의 순환 정리, 이상 유체 운동의 고전적 해, 로봇 공학적 응용을 학술적으로 정리한다.

2. 오일러 방정식의 유도

오일러 방정식은 유체 요소에 작용하는 힘과 가속도의 관계로부터 유도된다. 점성을 무시한 유체에서 응력 텐서는 압력과 관련된 등방성 부분만을 가지며, 운동량 방정식은 다음과 같이 표현된다.

\rho\,\frac{D\mathbf{u}}{Dt} = -\nabla p + \rho\,\mathbf{b}

여기서 D/Dt는 물질 도함수, \mathbf{b}는 단위 질량당 체적력이다. 중력장에서 \mathbf{b} = -g\hat{\mathbf{k}}이다. 이 식은 질량 보존의 연속 방정식과 결합되어 비점성 유체 운동의 지배 방정식을 구성한다.

물질 도함수의 전개

물질 도함수는 유체 입자를 따라가는 시간 변화율을 나타내며, 오일러 기술에서 다음과 같이 표현된다.

\frac{D\mathbf{u}}{Dt} = \frac{\partial\mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u}

여기서 첫째 항은 국소 가속도, 둘째 항은 대류 가속도이다. 오일러 방정식을 명시적으로 쓰면 다음과 같다.

\frac{\partial\mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u} = -\frac{1}{\rho}\nabla p + \mathbf{b}

이 형태는 유체역학의 수치 해석에서 직접 사용되며, 비정상과 정상 유동의 해석에 모두 적용된다.

성분 형태와 좌표계

오일러 방정식은 직교 좌표계, 원통 좌표계, 구 좌표계 등 다양한 좌표계에서 표현될 수 있다. 직교 좌표계에서의 세 성분 방정식은 다음과 같다.

\frac{\partial u}{\partial t} + u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} + w\frac{\partial u}{\partial z} = -\frac{1}{\rho}\frac{\partial p}{\partial x} + b_x

y와 z 성분도 유사한 형태를 가진다. 특정 문제의 대칭성에 따라 적절한 좌표계를 선택하면 해석이 크게 단순화된다. 축대칭 유동에서는 원통 좌표계가, 구 대칭 유동에서는 구 좌표계가 자연스럽다.

3. 경계 조건

오일러 방정식의 해를 결정하기 위해서는 적절한 경계 조건이 필요하다. 비점성 유체에서는 고체 표면에서의 비관통 조건(no-penetration)만 적용되며, 이는 유체 속도의 법선 성분이 고체 표면의 속도와 같아야 함을 의미한다.

\mathbf{u}\cdot\mathbf{n} = \mathbf{u}_{\text{wall}}\cdot\mathbf{n}

이 조건은 점성 유동의 점착 조건(no-slip)과 대비되며, 접선 속도에 대한 제약이 없음을 나타낸다. 이러한 경계 조건의 단순성은 오일러 해석이 실제 유동의 접선 속도를 정확히 예측하지 못하는 원인이 된다.

비회전 유동과 포텐셜 이론

오일러 방정식의 중요한 특수 경우는 비회전(irrotational) 유동이며, 이 경우 와도 \boldsymbol\omega = \nabla\times\mathbf{u} = 0이 성립한다. 비회전 유동은 속도 포텐셜 \phi로 기술되며, \mathbf{u} = \nabla\phi의 관계가 성립한다. 이를 연속 방정식에 대입하면 비압축성 비회전 유동의 경우 라플라스 방정식이 유도된다.

\nabla^2\phi = 0

이는 포텐셜 유동(potential flow) 이론의 기초이며, 복소 해석과 조화 함수 이론의 도구를 통해 엄밀한 해가 도출된다. 균일 유동, 소스와 싱크, 더블릿, 와류와 같은 기본 해의 중첩으로 복잡한 유동장을 구성할 수 있다.

4. 와도와 순환

와도는 유동의 국소적 회전 강도를 기술하는 벡터장이며, \boldsymbol\omega = \nabla\times\mathbf{u}로 정의된다. 순환(circulation)은 폐곡선을 따라의 속도의 선적분이며, 다음과 같이 정의된다.

\Gamma = \oint_C \mathbf{u}\cdot d\mathbf{l}

스토크스 정리에 의해 순환은 폐곡선으로 둘러싸인 면적의 와도 플럭스와 같다. 와도와 순환은 이상 유체 운동의 해석에서 핵심적 개념이며, 양력과 항력의 해석에 직접 연결된다.

Kelvin의 순환 정리

Kelvin의 순환 정리는 이상 유체의 비점성 비압축성 운동에서 물질 곡선의 순환이 시간에 따라 일정하게 보존됨을 기술한다.

\frac{D\Gamma}{Dt} = 0

이 정리는 이상 유체의 와도가 유체 입자와 함께 이동함을 의미하며, 와류 관(vortex tube)의 보존과 헬름홀츠 정리의 기반이 된다. 이러한 결과는 날개 주위의 순환과 양력의 해석에 직접 적용된다.

5. 이상 유체 운동의 고전적 해

오일러 방정식의 고전적 해는 여러 기본 유동 구성에 대해 엄밀히 도출되었다. 균일 유동은 일정한 속도 벡터로 기술되며, 자유 공간에서의 기본 유동이다. 소스와 싱크는 반경 방향 유동을 기술하며, 더블릿은 무한소로 가까운 소스와 싱크의 조합이다. 실린더 주위의 이상 유동은 균일 유동과 더블릿의 중첩으로 표현되며, 순환을 추가하면 양력이 발생하는 유동이 된다. 구 주위의 이상 유동도 유사한 방법으로 도출된다. 이러한 고전적 해는 공기역학의 기초를 제공한다.

6. D’Alembert의 역설

이상 유체 가정 하에서 균일 유동에 놓인 임의 형상의 물체에 작용하는 항력은 영이라는 결과가 엄밀히 도출된다. 이를 D’Alembert의 역설이라 한다. 이 결과는 실제 관측과 명백히 모순되며, 점성과 경계층의 중요성을 보여주는 고전적 예이다. 역설의 해결은 프란틀의 경계층 이론을 통해 이루어졌으며, 경계층의 얇은 영역에서 점성이 지배적으로 작용하여 실제 항력이 발생함을 설명한다. 이러한 해석은 유체역학의 역사적 발전에서 중요한 이정표이다.

7. 이상 유체 모델의 유효성과 한계

이상 유체 모델은 점성 효과가 무시될 수 있는 영역에서는 정확한 해를 제공하지만, 경계층과 후류(wake) 영역에서는 실제와 크게 다른 결과를 제공한다. 이상 유체 해석은 양력과 압력 분포의 대략적 예측에는 유용하지만, 항력과 실제 유동의 세부 구조를 기술하지 못한다. 실무적 공기역학에서는 이상 유체 해석과 경계층 해석을 결합하여 유동의 완전한 기술을 제공한다. 현대의 수치 해석에서는 나비에-스토크스 방정식의 직접 해석이 이러한 한계를 극복한다.

8. 로봇 공학적 응용

오일러 방정식과 이상 유체 이론은 드론의 날개와 프로펠러의 초기 설계, 수중 로봇 주위의 유동장의 대략적 해석, 고속 유동의 이론적 기준 해의 도출에 사용된다. 특히 포텐셜 유동 이론은 복잡한 형상 주위의 유동의 1차 근사 해를 제공하며, 패널 방법(panel method)을 통해 실용적 수치 해법으로 발전되었다. 이러한 방법은 CFD 해석보다 훨씬 빠른 계산을 제공하여 초기 설계 단계에서 유용하다.

9. 본 절의 의의

본 절은 오일러 방정식과 이상 유체 운동의 수학적 기초와 고전적 해, 로봇 공학적 응용을 체계적으로 정리한다. 이 방정식은 비점성 유체 운동의 표준 모델이며, 포텐셜 유동 이론과 함께 유체역학의 고전적 결과들의 출발점을 제공한다. 본 절의 내용은 이후 절에서 다룰 나비에-스토크스 방정식의 이해를 위한 중요한 기반이 된다.

10. 학습 권장사항

독자는 간단한 2차원 포텐셜 유동(소스, 싱크, 더블릿, 와류)의 중첩으로 실린더 주위의 유동장을 구성해 보고, 순환의 유무에 따른 양력의 발생을 관찰할 것을 권장한다. 또한 Kelvin의 순환 정리를 증명하고, 그 결과가 날개 주위의 양력 해석에 어떻게 활용되는지를 이해하는 것도 유익하다. D’Alembert의 역설을 직접 도출해 보는 것은 이상 유체 모델의 한계를 체감하는 좋은 방법이다.

11. 참고 문헌

• Euler, L. (1757). Principes généraux du mouvement des fluides. Mémoires de l’Académie Royale des Sciences et Belles-Lettres de Berlin, 11, 274–315.
• Batchelor, G. K. (2000). An Introduction to Fluid Dynamics. Cambridge University Press.
• Lamb, H. (1932). Hydrodynamics (6th ed.). Cambridge University Press.
• Milne-Thomson, L. M. (1968). Theoretical Hydrodynamics (5th ed.). Macmillan.
• Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill.
• Katz, J., & Plotkin, A. (2001). Low-Speed Aerodynamics (2nd ed.). Cambridge University Press.
• Kundu, P. K., Cohen, I. M., & Dowling, D. R. (2015). Fluid Mechanics (6th ed.). Academic Press.

version: 1.0