21.7 오일러 방정식과 이상 유체의 흐름

오일러 방정식(Euler equation)은 점성을 무시한 이상 유체에 대한 운동량 보존 법칙의 수학적 표현이며, 유체역학의 가장 고전적인 지배 방정식 중 하나이다. 이 방정식은 18세기 Leonhard Euler에 의해 도출되었으며, 유체 운동의 기본 구조를 명확히 제시한 최초의 체계적 이론이다. 오일러 방정식은 Navier-Stokes 방정식의 점성 항을 제외한 특수한 경우이며, 경계층 외부의 자유 흐름 영역이나 점성이 무시할 만한 조건에서 유동을 정확히 기술한다. 공기역학에서 오일러 방정식은 포텐셜 유동 이론, 베르누이 정리, 충격파 이론 등 여러 이론의 출발점이 되며, 수치 해석에서도 CFD의 기본 모델 중 하나로 사용된다. 본 절에서는 오일러 방정식의 유도, 수학적 구조, 보존 형태와 비보존 형태, 비압축성 및 압축성 형태, Bernoulli 정리와의 관계, 그리고 공기역학 응용을 체계적으로 정리한다.

1. 이상 유체의 정의와 가정

이상 유체(inviscid fluid)는 점성이 0인 가상의 유체이며, 전단 응력이 존재하지 않는다. 실제 유체는 모두 유한한 점성을 가지지만, 많은 공기역학 문제에서 점성의 영향이 국소적으로 제한되어 있거나(경계층 내부) 전반적으로 작은 경우 이상 유체 근사가 유용하다. 이상 유체 가정 하에서는 응력 텐서가 압력 항만을 포함하며, 운동량 방정식이 Navier-Stokes 방정식보다 크게 단순화된다. 이상 유체의 유동은 비점성 유동(inviscid flow)으로 불리며, 오일러 방정식으로 기술된다. 이 가정의 타당성은 유동의 특성과 해석의 목적에 따라 판단되어야 한다.

2. 오일러 방정식의 비보존 형태

뉴턴의 제2법칙을 유체 입자에 적용하고 점성을 무시하면, 오일러 방정식의 비보존 형태는 다음과 같이 표현된다.

\rho \dfrac{D\mathbf{V}}{Dt} = -\nabla p + \rho \mathbf{g}

좌변은 유체 입자의 가속도와 밀도의 곱이며, 우변은 압력 구배에 의한 힘과 중력 등의 체적력의 합이다. 물질 도함수를 전개하면 다음과 같다.

\rho \left(\dfrac{\partial \mathbf{V}}{\partial t} + (\mathbf{V} \cdot \nabla) \mathbf{V}\right) = -\nabla p + \rho \mathbf{g}

이 형태는 속도의 시간 변화와 공간적 변화(대류 가속도)를 명시적으로 보여주며, 유체 운동의 직관적 이해에 유용하다.

3. 오일러 방정식의 보존 형태

연속 방정식과 결합하여 오일러 방정식을 재배치하면, 보존 형태(conservation form)가 얻어진다.

\dfrac{\partial (\rho \mathbf{V})}{\partial t} + \nabla \cdot (\rho \mathbf{V} \otimes \mathbf{V}) + \nabla p = \rho \mathbf{g}

여기서 \mathbf{V} \otimes \mathbf{V}는 속도 벡터의 외적으로 형성되는 2차 텐서이다. 보존 형태는 질량 유량과 운동량 플럭스의 보존을 명시적으로 드러내며, 수치 해석에서의 충격파 처리와 질량·운동량 보존의 엄밀한 유지에 유리하다. 이 형태는 현대의 압축성 CFD 기법에서 표준적으로 사용된다.

21.7.4 비압축성 오일러 방정식

비압축성 유동에서 밀도가 상수이면, 오일러 방정식은 다음과 같이 단순화된다.

\dfrac{\partial \mathbf{V}}{\partial t} + (\mathbf{V} \cdot \nabla) \mathbf{V} = -\dfrac{1}{\rho} \nabla p + \mathbf{g}

이 비압축성 오일러 방정식은 연속 방정식 \nabla \cdot \mathbf{V} = 0과 함께 해석된다. 저속 공기역학의 해석에서는 이 형태가 가장 자주 사용되며, 포텐셜 유동 이론의 출발점이 된다. 점성을 무시하고 비압축성 가정을 적용한 이 형태는 경계층 외부의 자유 흐름, 익형 주변의 국소 유동(경계층 외), 단순화된 프로펠러 모델 등의 해석에 적합하다.

4. 압축성 오일러 방정식

압축성 유동에서는 밀도가 변수이므로, 연속 방정식, 운동량 방정식(오일러 방정식), 에너지 방정식, 상태 방정식이 결합된 완전한 시스템이 필요하다. 1차원 압축성 오일러 방정식의 보존 형태는 다음과 같이 표현된다.

\dfrac{\partial}{\partial t}\begin{pmatrix}\rho \\ \rho u \\ \rho E\end{pmatrix} + \dfrac{\partial}{\partial x}\begin{pmatrix}\rho u \\ \rho u^2 + p \\ (\rho E + p) u\end{pmatrix} = 0

여기서 E = e + u^2/2는 전체 에너지이다. 이 쌍곡선형(hyperbolic) 방정식 시스템은 특성 방법(method of characteristics)과 Riemann 해법으로 해석되며, 충격파와 팽창파의 전파를 기술한다. Euler 방정식은 충격파를 불연속으로 수용하는 대표적 예이며, Rankine-Hugoniot 조건이 불연속면에서의 보존 법칙을 규정한다.

21.7.6 Euler 방정식과 Navier-Stokes 방정식의 관계

Euler 방정식은 Navier-Stokes 방정식에서 점성 항을 제외한 특수한 경우이다. 점성이 중요한 영역(경계층 내부, 박리 영역, 후류)에서는 Navier-Stokes 방정식이 사용되어야 하며, 점성이 무시할 만한 영역(경계층 외부, 비점성 주류)에서는 Euler 방정식이 유효하다. Prandtl의 경계층 이론은 이 두 영역을 결합하는 고전적 접근이며, 전체 유동장을 비점성 주류와 점성 경계층으로 분할하여 해석한다. 현대의 CFD는 Navier-Stokes 방정식을 전 유동장에 대해 풀 수 있으나, 계산 효율을 위해 Euler 방정식이 여전히 중요한 도구로 사용된다.

21.7.7 Bernoulli 정리의 유도

Euler 방정식으로부터 Bernoulli 정리가 유도된다. 정상 비점성 비압축성 유동에서 유선(streamline)을 따라 운동량 방정식을 적분하면 다음 관계가 얻어진다.

p + \dfrac{1}{2} \rho V^2 + \rho g z = \text{const}

이는 Bernoulli 정리의 가장 일반적인 형태이며, 정압, 동압, 위치 압력의 합이 유선을 따라 일정함을 나타낸다. 비회전 유동의 경우에는 이 관계가 전체 유동장에서 성립하며, 유선에 한정되지 않는다. Bernoulli 정리는 양력 발생의 기본 이해, 피토관 기반 속도 측정, 벤추리 효과, 익형 주변의 압력 분포 등 여러 공기역학 응용의 기본 원리이다.

5. 비회전 유동과 속도 포텐셜

비회전 유동 \nabla \times \mathbf{V} = 0에서는 속도가 스칼라 포텐셜 \phi의 구배로 표현된다.

\mathbf{V} = \nabla \phi

이를 비압축성 연속 방정식 \nabla \cdot \mathbf{V} = 0에 대입하면 Laplace 방정식이 유도된다.

\nabla^2 \phi = 0

Laplace 방정식은 선형 편미분 방정식이며, 해석적 해법과 수치적 해법이 모두 잘 발달되어 있다. 포텐셜 유동 이론은 단순한 흐름 요소(균일 흐름, 점 소스·싱크, 이중자, 소용돌이)의 중첩으로 복잡한 유동을 구성하며, 익형 주변의 양력 해석, 물체 주변의 유동 가시화, 리프팅 라인 이론 등의 기반이 된다. 이 이론의 한계는 점성 효과와 경계층 박리를 포착할 수 없다는 점이다.

6. Bernoulli 정리의 압축성 확장

압축성 유동에서는 Bernoulli 정리가 수정되어야 한다. 등엔트로피 흐름에서 정상 비점성 유동에 대해 다음과 같은 관계가 유도된다.

\dfrac{\gamma}{\gamma-1} \dfrac{p}{\rho} + \dfrac{V^2}{2} = \text{const}

또는 엔탈피의 형태로 h + V^2/2 = \text{const}로 표현된다. 이 관계는 노즐 유동, 초음속 유동, 정체점 조건 등의 해석에 사용된다. 충격파가 존재하는 유동에서는 엔트로피가 증가하므로 이 관계가 충격파를 넘어서는 성립하지 않으며, Rankine-Hugoniot 조건이 대신 사용된다.

21.7.10 Euler 방정식의 수치 해석

Euler 방정식의 수치 해석은 CFD의 중요한 한 분야이다. 비점성 압축성 유동의 해석에는 유한 체적법, Godunov 기법, Roe 평균, Lax-Friedrichs 기법, 고차 정확도 기법(MUSCL, WENO 등)이 사용된다. 충격파의 정확한 포착, 인공적 진동의 억제, 엔트로피 조건의 만족이 수치 기법의 중요한 요구 사항이다. 현대의 압축성 Euler 솔버는 복잡한 기하학과 비정상 유동에 대해 고정밀 해를 제공하며, 고정익 UAV의 고속 비행, 프로펠러 블레이드 끝단의 천음속 해석, 충격파-경계층 상호작용 등의 연구에 활용된다. 비점성 해석은 점성 효과를 무시하므로 항력 평가에 한계가 있지만, 압축성 현상과 양력 계산에는 효율적이다.

21.7.11 Euler 방정식의 한계와 보완

Euler 방정식은 점성을 무시하므로, 점성이 지배적인 현상을 정확히 기술할 수 없다. 경계층, 박리, 후류, 층류-난류 천이, 마찰 항력 등은 Euler 방정식의 범위를 넘어서며, Navier-Stokes 방정식의 해석이 필요하다. 또한 비회전 가정 하의 포텐셜 유동 이론은 강한 비선형성과 비정상 현상에서 한계를 보인다. 실제 공기역학 해석에서는 Euler 방정식과 경계층 이론을 결합한 접근, 또는 Euler 방정식의 해를 Navier-Stokes 해석의 초기 조건으로 사용하는 접근 등이 사용된다. 이러한 조합적 접근은 계산 효율과 정확도의 균형을 제공한다.

21.7.12 로봇 공학에서의 응용

로봇 비행체의 공기역학 해석에서 Euler 방정식과 이상 유체 이론은 다양하게 활용된다. 첫째, 저속 드론의 익형 주변 양력 해석은 포텐셜 유동 이론(Euler 방정식의 특수한 경우)에 의해 수행되며, Bernoulli 정리가 압력 분포의 직관적 이해를 제공한다. 둘째, 프로펠러의 단순화된 모델(actuator disk, blade element theory)은 이상 유체 가정에 기반하며, 초기 설계와 성능 예측에 효율적이다. 셋째, 고정익 UAV의 고속 공력 해석에는 압축성 Euler 방정식의 수치 해석이 사용될 수 있다. 넷째, 풍동 시험 단면과 노즐 설계에서 1차원 Euler 이론이 기본 관계식을 제공한다. 다섯째, 공력 설계 최적화의 초기 단계에서는 Euler 기반 해석이 빠른 평가를 가능하게 하며, 후속의 Navier-Stokes 해석으로 정밀화된다.

21.7.13 요약과 후속 연결

Euler 방정식은 점성을 무시한 이상 유체의 운동량 보존 법칙이며, 유체역학의 고전적 지배 방정식이다. 비압축성과 압축성 형태, 보존 형태와 비보존 형태가 각각의 해석 맥락에서 사용되며, Bernoulli 정리와 포텐셜 유동 이론의 이론적 기반을 제공한다. Euler 방정식은 경계층 외부의 비점성 주류 해석과 충격파를 포함한 압축성 유동의 해석에 유용하며, CFD에서 중요한 모델로 사용된다. 점성이 지배적인 영역에서는 Navier-Stokes 방정식의 해석이 필요하며, 두 접근의 적절한 결합이 실무 해석의 표준이다. 다음 절에서는 Bernoulli 정리와 공기역학적 적용을 더 상세히 다루어, 양력 발생과 압력 분포의 이해에 필요한 기본 도구를 체계적으로 이해하는 관점을 제공한다.

출처

• Anderson, J. D., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2016.
• Landau, L. D., and Lifshitz, E. M., Fluid Mechanics, 2nd ed., Pergamon Press, 1987.
• Chorin, A. J., and Marsden, J. E., A Mathematical Introduction to Fluid Mechanics, 3rd ed., Springer, 1993.
• Saffman, P. G., Vortex Dynamics, Cambridge University Press, 1992.
• Toro, E. F., Riemann Solvers and Numerical Methods for Fluid Dynamics, 3rd ed., Springer, 2009.
• Leveque, R. J., Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, 2002.

버전

v1.0