21.3 공기의 물리적 성질과 상태 방정식
공기는 공기역학 해석의 대상 유체이며, 그 물리적 성질에 대한 정확한 이해는 유동 현상의 정량적 기술에 필수적이다. 공기는 질소, 산소, 아르곤, 이산화탄소 등 여러 성분 기체의 혼합물이지만, 공기역학의 일반적 해석에서는 균질한 이상 기체 또는 완전 기체로 근사된다. 공기의 주요 물리적 성질로는 밀도, 압력, 온도, 점성, 열전도도, 비열, 음속 등이 있으며, 이들은 상태 방정식과 열역학 관계식에 의해 서로 연관된다. 본 절에서는 공기의 조성과 기본 성질, 상태 방정식, 점성과 열전도, 비열비와 음속, 저속 근사와 압축성의 고려, 그리고 공기역학 해석에서의 응용을 체계적으로 정리한다.
1. 공기의 조성과 평균 분자량
건조 공기의 조성은 부피 기준으로 약 78% 질소(N_2), 21% 산소(O_2), 0.93% 아르곤(Ar), 0.04% 이산화탄소(CO_2), 그리고 미량의 네온, 헬륨, 메탄 등으로 구성된다. 이러한 조성은 지구 대기의 주요 부분에서 거의 일정하게 유지되며, 이를 균질 공기(homogeneous air)라 한다. 공기의 평균 분자량은 약 28.9647 g/mol이며, 이 값은 공기의 기체 상수 R_{\text{air}} = 287.058 J/(kg·K)의 계산에 사용된다. 습윤 공기는 수증기를 포함하며, 수증기의 분자량이 건조 공기보다 작으므로 습윤 공기의 밀도는 건조 공기보다 낮다. 고정밀 해석에서는 수증기의 영향을 반영한 습윤 공기 모델이 사용된다.
2. 이상 기체 상태 방정식
공기는 대기 조건에서 이상 기체(ideal gas)로 잘 근사되며, 상태 방정식은 다음과 같이 표현된다.
p = \rho R T
여기서 p는 압력, \rho는 밀도, R은 기체 상수, T는 절대 온도이다. 이 방정식은 공기의 세 가지 기본 상태 변수를 연결하며, 공기역학 해석에서 반복적으로 사용된다. 이상 기체 가정은 분자 간 힘이 무시할 만하고 분자의 크기가 분자 간 거리에 비해 매우 작다는 조건에서 성립하며, 대기 조건에서의 공기에 대해 높은 정확도를 제공한다. 극저온 또는 극고압 조건에서는 실제 기체 효과가 나타나며, 반데르 발스 방정식 등의 수정된 상태 방정식이 사용되지만, 일반적 공기역학에서는 이상 기체 가정이 충분하다.
21.3.3 공기의 비열과 비열비
공기의 비열은 정적 비열 c_v와 정압 비열 c_p로 구분된다. 정압 비열은 일정 압력 하에서 온도를 1 K 상승시키는 데 필요한 단위 질량당 열량이며, 공기의 경우 약 1005 J/(kg·K)이다. 정적 비열은 일정 부피 하에서의 대응 값이며, 공기의 경우 약 718 J/(kg·K)이다. 두 비열의 비율은 다음과 같이 정의되며, 공기역학과 열역학에서 핵심적 파라미터이다.
\gamma = \dfrac{c_p}{c_v}
공기의 경우 \gamma \approx 1.4이며, 이는 이원자 분자 이상 기체의 특징이다. 비열비는 음속, 단열 과정, 충격파, 압축성 효과의 정량화에 필수적으로 사용된다.
3. 공기의 점성과 Sutherland 공식
공기의 동점성 계수(dynamic viscosity) \mu는 온도에 따라 변화하며, Sutherland 공식으로 잘 근사된다.
\mu(T) = \mu_0 \dfrac{T_0 + S}{T + S} \left(\dfrac{T}{T_0}\right)^{3/2}
여기서 \mu_0은 기준 온도 T_0 = 273 K에서의 점성 1.716 \times 10^{-5} kg/(m·s), S = 110.4 K는 Sutherland 상수이다. 동점성 계수를 밀도로 나눈 운동 점성 계수(kinematic viscosity) \nu = \mu / \rho도 유체역학에서 자주 사용되며, 해수면 표준 조건에서 약 1.461 \times 10^{-5} m²/s이다. 점성은 공기와 물체 표면 사이의 전단 응력을 생성하며, 마찰 항력, 경계층, 층류-난류 천이의 물리적 기반이다.
21.3.5 공기의 열전도도
공기의 열전도도(thermal conductivity) k는 공기 내부의 열 전달 능력을 나타내며, 온도에 따라 변화한다. 해수면 표준 조건에서 공기의 열전도도는 약 0.0257 W/(m·K)이다. 열전도도는 열 경계층과 온도 분포의 해석에 중요하며, 고속 비행과 공력 가열 해석에서 필수적으로 사용된다. 점성 계수와 열전도도는 분자 운동론의 관점에서 비슷한 메커니즘으로 연관되며, Prandtl 수 Pr = \mu c_p / k는 두 확산 현상의 상대적 크기를 나타낸다. 공기의 Prandtl 수는 약 0.71이며, 이는 공기의 속도 경계층과 열 경계층이 비슷한 두께를 가짐을 의미한다.
21.3.6 공기 중 음속과 Mach 수
공기 중의 음속은 압력 섭동이 전파되는 속도이며, 이상 기체에 대해 다음과 같이 표현된다.
a = \sqrt{\gamma R T}
해수면 표준 조건(15°C)에서 음속은 약 340.3 m/s이며, 온도의 제곱근에 비례한다. 음속은 비행체 속도와의 비율인 Mach 수 M = V / a의 정의에 사용되며, 압축성 효과의 중요도를 결정한다. 일반적으로 M < 0.3에서는 밀도 변화가 5% 이내로 유지되어 비압축성 근사가 타당하며, 0.3 < M < 0.8에서는 아음속 압축성 효과가, 0.8 < M < 1.2에서는 천음속 효과가, M > 1에서는 초음속 효과가 중요해진다. 대부분의 드론과 소형 무인기는 저속 영역에서 동작하므로, 공기역학 해석은 비압축성 가정 하에서 수행된다.
4. 연속체 가정과 Knudsen 수
공기역학 해석은 유체를 연속체로 간주하는 가정에 기반한다. 이 가정이 유효한 조건은 분자의 평균 자유 행로(mean free path)가 특징 길이에 비해 충분히 작을 때이며, 이 비율은 Knudsen 수 Kn = \lambda / L로 정량화된다. Kn < 0.01이면 연속체 유체역학이 유효하며, 이는 대기 조건에서의 일반적 공기역학에 해당한다. 매우 작은 스케일(MEMS)이나 매우 높은 고도(희박 대기)에서는 Knudsen 수가 증가하여 분자 운동의 영향이 중요해지며, 분자 역학적 접근이 필요해진다. 일반적 드론과 무인기의 공기역학은 연속체 가정 하에서 문제없이 다루어진다.
5. 이상 기체와 완전 기체의 구분
공기역학에서 “이상 기체“와 “완전 기체“의 용어는 때로 구분되어 사용된다. 이상 기체는 상태 방정식 pv = RT를 만족하는 기체를 의미하며, 이는 분자 간 상호작용이 무시될 수 있는 조건에서 성립한다. 완전 기체(calorically perfect gas)는 이상 기체이면서 비열 c_p와 c_v가 온도에 무관하게 일정한 기체를 의미한다. 공기는 저온 영역에서는 완전 기체로 잘 근사되지만, 매우 높은 온도(예: 초음속 비행에 의한 공력 가열)에서는 비열이 온도에 의존하므로 열적 불완전 기체(thermally perfect gas)로 취급된다. 일반적 드론과 무인기의 공기역학에서는 완전 기체 가정이 충분한 정확도를 제공한다.
6. 열역학적 관계식
공기역학과 밀접하게 연관된 열역학 관계식으로는 이상 기체의 내부 에너지와 엔탈피의 표현이 있다. 이상 기체의 내부 에너지와 엔탈피는 온도만의 함수이며 다음과 같이 표현된다.
du = c_v dT, \quad dh = c_p dT
단열 가역 과정(등엔트로피 과정)에서는 다음 관계가 성립한다.
\dfrac{p}{\rho^{\gamma}} = \text{const}, \quad T \rho^{1-\gamma} = \text{const}, \quad T p^{(1-\gamma)/\gamma} = \text{const}
이러한 등엔트로피 관계식은 압축성 유동의 해석과 노즐 유동, 충격파 이전의 유동 등에서 반복적으로 사용된다.
7. 수증기와 습윤 공기의 취급
실제 대기에는 수증기가 포함되며, 수증기의 분자량이 건조 공기의 분자량보다 작으므로(약 18 g/mol 대 29 g/mol), 습윤 공기의 평균 분자량은 건조 공기보다 작고 밀도도 낮다. 습윤 공기의 밀도는 건조 공기의 부분 압력과 수증기의 부분 압력을 분리하여 계산된다. Dalton의 분압 법칙에 따라 전체 압력은 각 성분의 분압의 합이며, 이상 기체 가정 하에서 각 성분은 독립적 상태 방정식을 만족한다. 고정밀 공력 해석에서는 습윤 공기의 효과가 고려되며, 특히 열대 지역과 해상 운용 조건에서 중요하다.
8. 저속 공기역학에서의 성질 단순화
대부분의 드론과 소형 무인기는 저속 영역(M < 0.3)에서 운용되며, 이 조건에서는 공기의 밀도 변화가 무시할 만하다. 따라서 비압축성 가정(\rho = \text{const})이 적용되며, 상태 방정식과 에너지 방정식의 역할이 축소된다. 저속 공기역학에서는 연속 방정식과 운동량 방정식(Navier-Stokes 방정식)이 중심이 되며, 에너지 방정식은 독립적으로 해결될 수 있다. 온도와 밀도는 근사적으로 일정하게 취급되며, 공기의 점성만이 유동의 주요 성질로 남는다. 이러한 단순화는 저속 비행체의 해석을 크게 용이하게 한다.
9. 로봇 비행체의 해석에서의 공기 성질
로봇 비행체의 공기역학 해석에서 공기의 성질은 다음과 같이 활용된다. 첫째, 밀도는 양력과 항력의 계산에 직접 사용되며, 고도와 기상 조건에 따른 밀도 변화가 성능에 영향을 준다. 둘째, 점성은 레이놀즈 수의 계산에 사용되어 유동의 층류·난류 특성과 마찰 항력을 결정한다. 셋째, 음속은 Mach 수의 계산에 사용되며, 일반 드론의 경우 비압축성 가정이 유효하다. 넷째, 비열비는 압축성 유동과 로켓 추진 해석에 필요하지만, 일반 드론 해석에서는 간접적으로만 사용된다. 다섯째, 열전도도는 공력 가열과 온도 경계층이 중요한 특수 응용(예: 재돌입체)에서만 고려된다. 일반적 드론 설계에서는 표준 대기 조건의 공기 성질을 기반으로 해석이 수행되며, 필요에 따라 실측 대기 조건으로 보정된다.
10. 요약과 후속 연결
공기는 이상 기체로 잘 근사되며, 그 물리적 성질은 밀도, 압력, 온도, 점성, 열전도도, 비열, 음속 등으로 특성화된다. 이상 기체 상태 방정식은 세 가지 상태 변수를 연결하는 기본 관계식이며, Sutherland 공식은 점성의 온도 의존성을, a = \sqrt{\gamma R T}는 음속의 표현을 제공한다. 이러한 기본 성질은 공기역학 해석의 정량적 기반이며, 레이놀즈 수와 Mach 수 등의 무차원 매개변수를 통해 유동의 특성을 결정한다. 저속 드론과 무인기의 경우 비압축성 가정이 유효하며, 공기의 밀도와 점성이 주된 고려 대상이다. 다음 절에서는 유동장의 기본 변수와 정의를 다루어, 공기의 성질이 유동 해석의 수학적 기술로 연결되는 과정을 체계적으로 제시한다.
11. 출처
- Anderson, J. D., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2016.
- White, F. M., Viscous Fluid Flow, 3rd ed., McGraw-Hill, 2006.
- Schlichting, H., and Gersten, K., Boundary Layer Theory, 9th ed., Springer, 2017.
- Liepmann, H. W., and Roshko, A., Elements of Gasdynamics, Dover, 2001.
- Çengel, Y. A., and Boles, M. A., Thermodynamics: An Engineering Approach, 8th ed., McGraw-Hill, 2015.
- Shapiro, A. H., The Dynamics and Thermodynamics of Compressible Fluid Flow, Vol. 1, Ronald Press, 1953.
12. 버전
v1.0