Chapter 21. 공기역학 기초 (Fundamentals of Aerodynamics)

1. 유동장의 기술과 대기의 물리적 성질

공기역학은 공기 중에서 운동하는 물체 주변의 유동장을 기술하고, 그로부터 발생하는 힘과 모멘트를 정량적으로 해석하는 학문이다. 유동장은 공간과 시간의 함수로 정의되는 속도장 \mathbf{u}(\mathbf{x}, t), 압력장 p(\mathbf{x}, t), 밀도장 \rho(\mathbf{x}, t), 온도장 T(\mathbf{x}, t)의 집합으로 표현되며, 대기를 연속체 가정 아래에서 다루는 전제에서 출발한다. 대류권 대기는 근사적으로 이상 기체 상태 방정식 p = \rho R T로 기술되며, 해수면 표준 조건에서 공기의 밀도는 약 \rho = 1.225\ \mathrm{kg/m^3}, 점성 계수는 \mu \approx 1.81 \times 10^{-5}\ \mathrm{Pa\cdot s}, 동점성 계수는 \nu = \mu/\rho \approx 1.48 \times 10^{-5}\ \mathrm{m^2/s}로 주어진다. 이러한 물성값은 온도와 고도에 따라 변화하며, 점성 계수의 온도 의존성은 Sutherland 법칙 \mu(T) = \mu_0 (T/T_0)^{3/2} (T_0 + S) / (T + S)로 근사된다.

유동장을 기술하는 관점은 공간의 한 점에서 유동 변수의 시간 변화를 추적하는 오일러 기술(Eulerian description)과 특정 유체 입자를 따라가며 상태를 추적하는 라그랑주 기술(Lagrangian description)로 구분되며, 두 기술의 연결은 물질 미분 D/Dt = \partial/\partial t + \mathbf{u}\cdot\nabla로 이루어진다. 공기역학의 수학적 전개는 대부분 오일러 기술을 채택하고, 유선(streamline), 유적선(pathline), 와선(streakline) 등의 기하학적 개념을 이용하여 유동 구조를 시각화한다. 또한 압축성의 유무는 마하 수 \mathrm{Ma} = V/a에 의하여 판정되며, 비행 로봇이 주로 운용되는 \mathrm{Ma} < 0.3 영역에서는 밀도 변화가 5% 미만으로 유지되어 비압축성 가정이 타당하게 적용된다. 이러한 전제 위에서 이후의 지배 방정식이 간결한 형태로 축약된다.

2. 유체 지배 방정식과 베르누이 정리

유체의 연속체 역학은 질량, 운동량, 에너지 보존 법칙에 의하여 지배되며, 이를 미분 형태로 표현한 것이 공기역학의 지배 방정식이다. 질량 보존은 연속 방정식 \partial \rho / \partial t + \nabla\cdot(\rho \mathbf{u}) = 0으로 표현되고, 비압축성 유동에서는 이를 \nabla\cdot\mathbf{u} = 0으로 단순화할 수 있다. 운동량 보존은 Navier–Stokes 방정식

\rho \frac{D\mathbf{u}}{Dt} = -\nabla p + \mu \nabla^{2}\mathbf{u} + \rho \mathbf{g}

로 기술되며, 점성을 무시할 수 있는 영역에서는 Euler 방정식 \rho \, D\mathbf{u}/Dt = -\nabla p + \rho\mathbf{g}로 축약된다. 에너지 보존은 전체 엔탈피 h_0 = h + \tfrac{1}{2}|\mathbf{u}|^2의 수송 방정식으로 표현되며, 단열·비점성·정상 유동의 조건에서는 h_0가 유선을 따라 일정하게 유지된다. 이러한 보존 법칙은 유동 전반에 걸쳐 성립하며, 공기역학의 모든 정량적 해석은 이 미분 방정식들의 적절한 이상화와 경계 조건 부여로부터 시작된다.

이상 유동의 정상·비압축·비점성 조건에서는 운동량 방정식을 유선을 따라 적분하여 베르누이 정리

p + \tfrac{1}{2}\rho V^{2} + \rho g h = \text{const}

를 얻는다. 이 관계는 공기역학에서 가장 널리 사용되는 근사적 도구로서, 익형 주위의 압력 분포 추정, 피토 튜브를 이용한 속도 측정, 프로펠러 디스크 이론의 유도 등에 활용된다. 다만 베르누이 관계는 점성의 영향이 무시되지 않는 경계층 내부, 분리된 유동, 강한 회전 유동에서는 직접 적용될 수 없으며, 이러한 제한은 이후의 경계층 이론과 점성 유동 해석을 통하여 체계적으로 보완된다. 또한 압축성이 본격적으로 작용하는 고속 영역에서는 엔탈피가 일정한 조건 아래 확장된 형태의 베르누이 관계 h + \tfrac{1}{2}V^2 = \text{const}가 사용되며, 이는 비행 로봇의 일반적 작동 영역보다는 고속 항공기 해석에서 주로 요구된다.

3. 경계층, 무차원 매개변수, 해석 방법론

실제 유동에서는 점성에 의하여 물체 표면 근처에 얇은 경계층(boundary layer)이 형성되며, 경계층 내부에서는 속도가 벽면의 점착 조건에 따른 0으로부터 자유류 속도까지 점진적으로 변화한다. Prandtl의 경계층 방정식은 벽 근방의 운동량 방정식을 두께 방향으로 축약한 형태로,

u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = -\frac{1}{\rho}\frac{\partial p}{\partial x} + \nu \frac{\partial^{2} u}{\partial y^{2}}

로 표현된다. 경계층은 Reynolds 수 \mathrm{Re} = \rho V L / \mu에 의하여 층류와 난류로 구분되며, 비행 로봇의 익현 기반 \mathrm{Re}는 대체로 10^{4}에서 10^{6} 사이에 분포한다. 이 영역은 층류 분리와 천이가 공존하는 대역이며, 저레이놀즈 영역의 익형 성능이 전통적 항공공학의 고레이놀즈 영역과 뚜렷이 구별되는 원인을 제공한다. 경계층의 박리와 재부착, 천이 위치, 두께 성장률은 표면 압력 분포와 벽 전단 응력을 결정하여 양력과 항력의 발생 양상에 직접적인 영향을 미친다.

공기역학의 해석에는 다양한 무차원 매개변수가 사용되며, 대표적으로 앞서 언급한 Reynolds 수와 Mach 수 외에 Strouhal 수 \mathrm{St} = fL/V, Froude 수 \mathrm{Fr} = V/\sqrt{gL}, Prandtl 수 \mathrm{Pr} = \nu/\alpha 등이 유동 현상의 상사성과 지배 기구를 식별하는 지표로 활용된다. 이러한 무차원화는 동일한 물리 법칙 아래 서로 다른 규모의 유동을 상호 비교할 수 있도록 하며, 풍동 시험과 실기 운용 사이의 축적비 상사성을 확보하는 이론적 기반을 제공한다. 표 21.1은 대표적 무차원 수의 정의와 물리적 의미를 정리한다.

해석 방법론은 이론적 접근, 수치적 접근, 실험적 접근의 삼각 구도로 구성된다. 이론적 접근은 포텐셜 유동 이론, 등각 사상 기법, 섭동 해석 등 폐형 해를 구성하는 방법을 포괄하며, 얇은 익형 이론과 리프팅 라인 이론과 같은 고전적 결과를 산출한다. 수치적 접근은 유한 차분, 유한 체적, 유한 요소, 격자 볼츠만, 패널법 등을 포함하는 전산 유체역학(CFD)으로 대표되며, Reynolds 평균 Navier–Stokes(RANS), 대와류 모사(LES), 직접 수치 모사(DNS) 등 서로 다른 수준의 난류 모델링이 적용된다. 실험적 접근은 풍동 시험을 중심으로 하며, PIV, LDV, 핫 와이어, 표면 압력 측정, 힘 균형 등 다양한 계측 기법이 통합적으로 활용된다. 비행 로봇의 공기역학 해석은 이러한 세 접근을 조합하여, 예비 설계 단계에서는 이론과 저차 수치 모델을, 상세 설계 단계에서는 CFD와 풍동 시험을 병행하는 방식이 표준으로 자리 잡아 왔다.

4. 출처

• 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.
• Schlichting, H., and Gersten, K., Boundary-Layer Theory, 9th ed., Springer, 2017.
• White, F. M., Viscous Fluid Flow, 3rd ed., McGraw-Hill, 2006.
• Katz, J., and Plotkin, A., Low-Speed Aerodynamics, 2nd ed., Cambridge University Press, 2001.
• Pope, S. B., Turbulent Flows, Cambridge University Press, 2000.

5. 버전

v1.0