21.6 연속 방정식의 공기역학적 적용

연속 방정식(continuity equation)은 질량 보존 법칙을 유체역학의 관점에서 수학적으로 기술한 것이며, 공기역학의 세 가지 기본 지배 방정식 중 하나이다. 질량은 어떠한 물리 과정에서도 생성되거나 소멸되지 않으므로, 임의의 제어 체적에서 질량의 시간에 따른 변화는 체적 경계를 통해 유입되는 질량과 유출되는 질량의 차이에 의해서만 결정된다. 이 보존 법칙은 적분 형태와 미분 형태로 표현되며, 다양한 공기역학 응용에서 필수적으로 사용된다. 본 절에서는 연속 방정식의 유도, 적분 형태와 미분 형태의 관계, 비압축성과 압축성 유동에서의 단순화, 원통 좌표계와 구면 좌표계에서의 표현, 그리고 실제 공기역학 문제에서의 적용을 체계적으로 정리한다.

1. 질량 보존의 기본 원리

질량 보존은 물리학의 가장 근본적인 법칙 중 하나이며, 고립된 시스템에서 질량의 총합은 시간에 따라 변하지 않는다는 명제이다. 열역학적 개방 시스템에서는 시스템의 경계를 통한 질량의 유출입이 허용되며, 이 경우에도 시스템의 내부 질량 변화율은 유입과 유출의 차이와 같다. 유체역학에서는 이 원리를 임의로 선택된 공간 영역(제어 체적)에 적용하며, 유체의 유입과 유출이 연속적이고 점 단위로 발생한다는 특성을 반영한다. 이 원리는 미시적 입자 수준에서도 성립하며, 연속체 가정 하에서 거시적 장 방정식으로 확장된다.

2. 적분 형태의 연속 방정식

임의의 공간 고정 제어 체적 $V$ 에 대해 질량 보존을 적용하면, 적분 형태의 연속 방정식이 다음과 같이 유도된다.

$\dfrac{\partial}{\partial t} \int_V \rho \, dV + \oint_S \rho \mathbf{V} \cdot \mathbf{n} \, dS = 0$

여기서 첫 항은 제어 체적 내부의 질량의 시간 변화율, 둘째 항은 제어 체적 경계 $S$ 를 통한 질량 유출률이며, $\mathbf{n}$ 은 경계의 외향 단위 법선 벡터이다. 이 적분 형태는 제어 체적의 임의성과 경계면에 대한 명시적 표현으로 인해 실무적 계산에서 직관적이며, 덕트 유동, 노즐 유동, 프로펠러 후류 등의 해석에 자주 사용된다. 정상 유동의 경우 첫 항이 0이 되어 유입 유량과 유출 유량이 같다는 단순한 관계로 축약된다.

21.6.3 미분 형태의 연속 방정식

적분 형태의 연속 방정식에 발산 정리(divergence theorem)를 적용하면 미분 형태가 유도된다.

$\dfrac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{V}) = 0$

이 방정식은 공간의 각 점에서 성립해야 하며, 유동장의 점 단위 해석에 사용된다. 첫 항은 특정 점에서의 밀도의 시간 변화율이고, 둘째 항은 질량 유량 벡터 $\rho \mathbf{V}$ 의 발산이다. 미분 형태는 편미분 방정식의 형태를 가지며, 운동량 방정식과 에너지 방정식과 결합하여 완전한 유체역학 시스템을 구성한다. 물질 도함수를 사용하여 다음과 같이 재배치할 수도 있다.

$\dfrac{D\rho}{Dt} + \rho \nabla \cdot \mathbf{V} = 0$

이 형태는 유체 입자를 따라가면서 관측한 밀도 변화율과 속도의 발산의 관계를 명확히 보여준다.

21.6.4 비압축성 유동의 연속 방정식

비압축성 유동에서는 밀도가 상수이므로 $D\rho/Dt = 0$ 이며, 연속 방정식은 다음과 같이 단순화된다.

$\nabla \cdot \mathbf{V} = 0$

이는 속도장의 발산이 0인 조건이며, 유동의 체적이 보존됨을 의미한다. 비압축성 유동의 연속 방정식은 운동량 방정식과 독립적으로 제약 조건으로 작용하며, 압력은 이 조건을 만족시키기 위한 Lagrange 승수의 역할을 한다. 직교 좌표계에서는 다음과 같이 전개된다.

$\dfrac{\partial u}{\partial x} + \dfrac{\partial v}{\partial y} + \dfrac{\partial w}{\partial z} = 0$

이 단순한 형태는 저속 공기역학의 해석에서 반복적으로 사용된다.

21.6.5 1차원 유동의 연속 방정식

1차원 정상 유동에서 연속 방정식은 단면적과 속도의 관계로 축약되며, 단면적 $A$ 를 가진 덕트를 따라가는 유동에 대해 다음과 같이 표현된다.

$\rho A V = \text{const}$

이를 질량 유량 $\dot{m} = \rho A V$ 로 표현하면, 정상 1차원 유동에서 질량 유량이 일정하다는 결론이 도출된다. 비압축성 경우에는 더 단순히 $A V = \text{const}$ 가 되며, 단면적이 감소하면 속도가 증가하고 단면적이 증가하면 속도가 감소하는 관계가 나타난다. 이는 노즐, 디퓨저, 벤추리 미터, 풍동의 시험 단면 등 여러 실무 응용의 이해에 기초를 제공한다.

3. 압축성 1차원 유동의 면적-속도 관계

압축성 유동에서는 밀도의 변화를 고려해야 하며, 1차원 정상 유동의 연속 방정식을 미분하면 다음과 같다.

$\dfrac{d\rho}{\rho} + \dfrac{dA}{A} + \dfrac{dV}{V} = 0$

Mach 수에 따라 면적 변화와 속도 변화의 관계는 다음과 같이 나타난다.

$\dfrac{dA}{A} = (M^2 - 1) \dfrac{dV}{V}$

$M < 1$ 의 아음속 영역에서는 면적이 감소하면 속도가 증가하고, $M > 1$ 의 초음속 영역에서는 면적이 증가해야 속도가 증가한다. 이는 라발 노즐(Laval nozzle)의 수렴-발산 형상과 초음속 풍동의 원리를 설명하는 기본 관계이다.

4. 원통 좌표계와 구면 좌표계에서의 연속 방정식

유동의 기하학에 따라 원통 좌표계 또는 구면 좌표계에서의 연속 방정식이 유용하다. 원통 좌표계 $(r, \theta, z)$ 에서 비압축성 유동의 연속 방정식은 다음과 같다.

$\dfrac{1}{r}\dfrac{\partial (r u_r)}{\partial r} + \dfrac{1}{r}\dfrac{\partial u_\theta}{\partial \theta} + \dfrac{\partial u_z}{\partial z} = 0$

이 형태는 축대칭 유동, 회전 유동, 파이프 유동 등의 해석에 사용된다. 구면 좌표계에서는 구형 물체 주변의 유동과 점 소스·싱크 모델에 활용된다. 좌표계의 선택은 문제의 대칭성과 경계 조건에 따라 결정되며, 적절한 좌표계의 선택은 해석의 효율을 크게 향상시킨다.

21.6.8 유동 함수와 연속 방정식

2차원 비압축성 유동에서는 유동 함수(stream function) $\psi$ 가 정의될 수 있으며, 속도 성분은 다음과 같이 표현된다.

$u = \dfrac{\partial \psi}{\partial y}, \quad v = -\dfrac{\partial \psi}{\partial x}$

이 정의는 연속 방정식을 자동으로 만족시키며, 이 표현을 통해 연속 방정식의 제약이 해소된다. 유동 함수는 2차원 유동의 해석을 크게 단순화하며, 유선의 가시화와 유량 계산에도 유용하다. 유동 함수의 차이는 두 유선 사이의 단위 폭당 유량과 같으며, 물리적 직관을 제공한다. 3차원 비압축성 유동에서는 벡터 포텐셜이 유사한 역할을 한다.

5. 정상 유동과 비정상 유동에서의 적용

정상 유동에서 연속 방정식은 시간 편미분 항이 사라지며, $\nabla \cdot (\rho \mathbf{V}) = 0$ 이 된다. 비정상 유동에서는 시간 편미분 항이 포함되어야 하며, 밀도 파의 전파, 맥동 유동, 과도 현상 등의 해석에서 중요한 역할을 한다. 공기역학에서는 대부분의 해석이 정상 유동 가정 하에서 수행되지만, 돌풍 응답, 블레이드 통과 주파수, 비정상 공력 하중의 해석에서는 비정상 효과가 필수적으로 고려된다.

6. 제어 체적 해석의 응용

적분 형태의 연속 방정식은 제어 체적 해석(control volume analysis)의 기본 도구이며, 유동 시스템의 통합적 성능을 평가하는 데 사용된다. 프로펠러 후류의 질량 유량 계산, 덕트 유동의 유량 보존, 풍동의 시험 단면 설계 등이 대표적 응용이다. 프로펠러 운동량 이론에서 연속 방정식은 후류의 속도 분포를 결정하는 기본 관계 중 하나이며, 추력과 효율의 계산에 기여한다. 제어 체적의 선택은 문제의 특성에 따라 유연하게 이루어지며, 적절한 선택은 해석의 단순성과 직관성을 높인다.

7. 수치 해석에서의 연속 방정식

전산 유체역학(CFD)에서 연속 방정식은 이산화된 형태로 수치적으로 해결된다. 유한 체적법에서는 각 격자 셀에 대한 질량 보존이 엄밀히 유지되며, 유한 차분법에서는 연속 방정식이 이산화 점에서 만족되어야 한다. 비압축성 유동의 수치 해석에서는 압력-속도 결합 알고리즘(SIMPLE, PISO 등)이 연속 방정식의 제약을 만족시키기 위해 사용된다. 압축성 유동의 수치 해석에서는 밀도가 주 변수로 계산되며, 연속 방정식과 운동량·에너지 방정식이 연립하여 해결된다. 수치 해석의 정확도와 안정성은 연속 방정식의 이산화 방식에 크게 의존한다.

8. 로봇 공학에서의 응용

로봇 비행체의 공기역학 해석에서 연속 방정식은 다양하게 활용된다. 첫째, 프로펠러와 로터의 운동량 이론에서 연속 방정식은 후류의 속도 분포와 질량 유량을 결정한다. 둘째, 덕트 프로펠러와 슈라우드 설계에서 단면적 변화에 따른 속도 변화의 분석에 사용된다. 셋째, 드론의 전산 유체역학 해석에서 연속 방정식은 지배 방정식의 필수 구성 요소이며, 이산화된 형태로 해결된다. 넷째, 풍동 시험 단면의 설계와 흐름 균일성 평가, 공력 측정 장비의 교정 등에도 연속 방정식이 기본 원리로 사용된다. 다섯째, 대기 환경 시뮬레이션에서 바람장의 연속성 조건이 수치 해석의 기본 제약이 된다.

9. 요약과 후속 연결

연속 방정식은 질량 보존의 유체역학적 표현이며, 적분 형태와 미분 형태로 기술된다. 비압축성 유동에서는 속도장의 발산이 0이 되는 단순한 형태로 축약되며, 압축성 유동에서는 밀도 변화가 포함되는 완전한 형태가 사용된다. 1차원, 2차원, 3차원 유동에서 각각의 적용 형태가 있으며, 유동 함수, 면적-속도 관계 등의 파생 개념이 유용하게 활용된다. 연속 방정식은 운동량 방정식 및 에너지 방정식과 결합하여 완전한 공기역학의 지배 시스템을 구성하며, 이론, 수치, 실험의 모든 접근에서 기본 도구로 사용된다. 다음 절에서는 Euler 방정식과 이상 유체의 흐름을 다루어, 비점성 유동의 운동량 보존 법칙과 그 해석적 기법을 체계적으로 이해하는 관점을 제공한다.

10. 출처

Anderson, J. D., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2016.
Kundu, P. K., Cohen, I. M., and Dowling, D. R., Fluid Mechanics, 6th ed., Academic Press, 2015.
White, F. M., Fluid Mechanics, 8th ed., McGraw-Hill, 2016.
Panton, R. L., Incompressible Flow, 4th ed., Wiley, 2013.
Batchelor, G. K., An Introduction to Fluid Dynamics, Cambridge University Press, 2000.
Versteeg, H. K., and Malalasekera, W., An Introduction to Computational Fluid Dynamics: The Finite Volume Method, 2nd ed., Pearson, 2007.

11. 버전

v1.0