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

1. 질량 보존의 적분 형태와 미분 형태

연속 방정식은 질량 보존 법칙을 유체의 연속체 역학에 맞게 기술한 관계이며, 운동량·에너지 보존과 함께 공기역학의 세 기본 지배 방정식을 이룬다. 임의의 공간에 고정된 제어 체적 $\mathcal{V}$ 와 그 경계면 $\partial \mathcal{V}$ 에 대하여, 제어 체적 내부의 총 질량 변화는 경계를 통한 질량 플럭스의 총합과 같아야 한다. 이를 수식으로 표현하면

$\frac{d}{dt}\int_{\mathcal{V}} \rho \, dV + \oint_{\partial \mathcal{V}} \rho \mathbf{u}\cdot \mathbf{n}\, dA = 0$

이 된다. 여기서 $\mathbf{n}$ 은 외향 단위 법선 벡터이며, 첫 항은 제어 체적 내부의 질량 축적률, 두 번째 항은 경계를 통과하는 질량의 순 유출률을 의미한다. 이 적분 형태는 제어 체적 해석의 기본 도구로서, 덕트·노즐·엔진 흡·배기·프로펠러 디스크와 같은 공학적 대상에 대하여 직접적이고 물리적으로 명확한 해석을 제공한다.

적분 형태에 Gauss의 발산 정리를 적용하면 경계 적분이 체적 적분으로 전환되며, 제어 체적의 임의성을 이용하여 피적분 함수가 0이어야 한다는 조건으로부터 미분 형태의 연속 방정식

$\frac{\partial \rho}{\partial t} + \nabla\cdot(\rho \mathbf{u}) = 0$

이 도출된다. 이 미분 형태는 유동장의 국소적 질량 보존을 표현하며, Navier–Stokes 방정식과 에너지 방정식과 결합하여 공기역학의 완전한 지배 방정식 체계를 구성한다. 물질 도함수를 이용하여 정리하면 $D\rho/Dt + \rho \nabla\cdot\mathbf{u} = 0$ 의 형태가 되며, 이는 유체 입자의 밀도 변화가 속도장의 발산과 반대 부호로 비례함을 보여 준다. 즉 유체 입자가 팽창하면 $\nabla\cdot\mathbf{u} > 0$ 이 되어 밀도가 감소하고, 수축하면 밀도가 증가한다는 직관적 의미를 제공한다.

2. 비압축성·정상·좌표계별 단순화

비압축성 유동에서는 밀도가 유체 입자를 따라 일정하게 유지되어 $D\rho/Dt = 0$ 이 성립하며, 미분 형태의 연속 방정식은

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

으로 축약된다. 이 조건은 속도장이 발산 자유임을 의미하고, 2차원 유동에서는 유동 함수 $\psi$ 의 존재를 허용하여 $u = \partial\psi/\partial y$ , $v = -\partial\psi/\partial x$ 로 표현할 수 있게 한다. 정상 유동에서는 추가로 $\partial \rho/\partial t = 0$ 이 성립하여 압축성 유동의 연속 방정식이 $\nabla\cdot(\rho\mathbf{u}) = 0$ 으로 줄어든다. 또한 두 단면 $A_{1}$ , $A_{2}$ 를 가지는 유관(stream tube)에 대한 적분 형태는 $\rho_{1} V_{1} A_{1} = \rho_{2} V_{2} A_{2}$ 의 질량 유량 불변 조건을 제공하며, 비압축성 극한에서는 $V_{1} A_{1} = V_{2} A_{2}$ 의 체적 유량 불변 형태로 단순화된다. 이러한 단순화는 풍동 수축부, 환기 덕트, 노즐, 프로펠러 유관의 해석에서 기초적이면서도 강력한 도구로 활용된다.

연속 방정식은 문제의 기하에 맞추어 다양한 좌표계로 표현된다. 원통 좌표계 $(r, \theta, z)$ 에서 비압축성 형태는

$\frac{1}{r}\frac{\partial (r u_{r})}{\partial r} + \frac{1}{r}\frac{\partial u_{\theta}}{\partial \theta} + \frac{\partial u_{z}}{\partial z} = 0$

로 주어지며, 구면 좌표계 $(r, \theta, \varphi)$ 에서는

$\frac{1}{r^{2}}\frac{\partial (r^{2} u_{r})}{\partial r} + \frac{1}{r\sin\theta}\frac{\partial (u_{\theta}\sin\theta)}{\partial \theta} + \frac{1}{r\sin\theta}\frac{\partial u_{\varphi}}{\partial \varphi} = 0$

의 형태로 표현된다. 원통 좌표계는 회전 대칭성을 가진 프로펠러·로터·제트 해석에 적합하며, 구면 좌표계는 흡입구 근방의 반경 방향 유입이나 포인트 소스 근사 해석에서 유용하다. 이러한 좌표계 선택은 수학적 편의의 문제에 그치지 않고, 기하학적 대칭성을 직접 반영하여 방정식의 차원을 효과적으로 축소함으로써 해석의 깊이를 확대하는 수단으로 기능한다.

3. 공기역학 해석에서의 적용 사례

연속 방정식은 공기역학의 다양한 해석 절차에서 핵심적 제약 조건으로 작동한다. 풍동의 수축부에서는 상류 대단면 $A_{1}$ 으로부터 시험 단면 $A_{2}$ 로의 속도 증가가 질량 유량 불변 $V_{2} = (A_{1}/A_{2}) V_{1}$ 로부터 예측되며, 덕트 프로펠러의 성능 해석에서도 유입·유출 단면의 속도 분포가 이 관계를 통하여 연결된다. 또한 프로펠러의 운동량 이론에서는 디스크 상류와 하류의 유관을 통과하는 질량 유량이 일정해야 한다는 조건으로부터 유도 속도 $v_{i}$ 와 추력 $T$ 사이의 관계가 확립된다. 이러한 해석은 해석적 폐형 해에서 수치 모사에 이르기까지 일관된 형태로 사용되며, 연속 방정식이 단순한 이론적 장치가 아니라 구체적 설계 식을 유도하는 공학적 근간임을 보여 준다.

CFD 수치 해석에서도 연속 방정식은 동일한 중심성을 유지한다. 비압축성 Navier–Stokes 솔버에서는 연속 방정식이 속도장의 발산 자유 조건을 부과하므로, 속도와 압력을 동시에 일관적으로 갱신하기 위하여 SIMPLE, PISO, fractional step, projection method와 같은 투영 기법이 사용된다. 이들 기법의 공통된 구조는 예측된 속도장에 대한 발산을 제거하도록 압력에 대한 Poisson 방정식을 풀어 속도장을 보정하는 것이며, 이 절차 전체는 연속 방정식의 대수적 표현이 수치 해의 일관성을 지배한다는 원리에 기반한다. 압축성 Navier–Stokes 해석에서는 밀도 기반 방법이 사용되어, 연속 방정식이 $\rho$ 의 시간 전진 방정식으로 직접 활용된다. 이러한 구현의 차이는 연속 방정식이 수치 기법 선택의 핵심 기준이 됨을 보여 준다.

4. 로봇공학적 함의

비행 로봇의 해석에서 연속 방정식은 여러 층위에서 사용된다. 프로펠러 기반 추진계의 예비 설계에서는 디스크 양쪽의 질량 유량이 같아야 한다는 조건과 운동량 보존을 결합하여 유도 속도 $v_{i} = \sqrt{T/(2\rho A_{\mathrm{disk}})}$ 을 얻으며, 이는 호버링 추력 계산과 전력 소요 추정의 출발점이 된다. 또한 덕트가 있는 프로펠러와 리프팅 덕트 프로펠러의 해석에서는 덕트 내부의 질량 유량이 단면을 통하여 일정하게 유지되므로, 덕트 형상 설계는 연속 방정식의 제약 아래에서 이루어진다. 흡기·배기가 있는 이상 압축 작동 기체 시스템에서도 유사한 제약이 부과되며, 비행 로봇의 열 관리·환기 설계에서 기본적 제어 변수로 작동한다. 이러한 사례는 연속 방정식이 추상적인 편미분 관계가 아니라 구체적인 공학적 선택을 규정하는 기초 방정식임을 보여 준다.

비행 로봇의 제어·추정 분야에서도 연속 방정식의 의의가 유지된다. 실시간 공력 모델은 비압축성 가정에서 출발하여 속도장의 발산 자유 조건을 자연스럽게 만족시키며, 이로 인하여 양력 및 항력 계수의 스케일링이 단순한 밀도 의존 형태로 표현될 수 있다. 한편 밀도와 속도의 변화가 유의한 영역에서는 압축성 연속 방정식이 요구되며, 이는 대기 데이터 시스템과의 결합을 통하여 온도·압력의 실시간 변화를 공력 모델에 반영하는 절차로 구체화된다. 이러한 통합적 사용은 연속 방정식이 비행 로봇의 설계·해석·운용 전반에서 지속적으로 참조되는 물리적·수학적 제약임을 분명히 해 준다.

5. 출처

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.
Batchelor, G. K., An Introduction to Fluid Dynamics, Cambridge University Press, 2000.
White, F. M., Viscous Fluid Flow, 3rd ed., McGraw-Hill, 2006.
Ferziger, J. H., Peric, M., and Street, R. L., Computational Methods for Fluid Dynamics, 4th ed., Springer, 2020.
Leishman, J. G., Principles of Helicopter Aerodynamics, 2nd ed., Cambridge University Press, 2006.

6. 버전

v1.0