21.3 공기의 물리적 성질과 상태 방정식

21.3 공기의 물리적 성질과 상태 방정식

1. 공기의 조성과 열역학적 상태

공기는 공기역학 해석의 대상 유체이며, 그 물리적 성질에 대한 정확한 정의는 유동 현상의 정량적 기술을 위한 필수 전제이다. 건조 공기는 체적 기준으로 질소 약 78.08%, 산소 약 20.95%, 아르곤 약 0.93%, 이산화탄소 약 0.04%로 구성된 혼합 기체이며, 이로부터 산출되는 평균 분자량은 M_{\text{air}} \approx 28.9647\ \mathrm{g/mol}로 주어진다. 실제 대기에는 수증기가 포함되어 있으나, 공기역학의 일반적 해석에서는 건조 공기를 이상 기체로 근사하는 것이 표준적 관례이다. 이 근사는 대기권의 낮은 고도와 중간 속도 영역에서 충분한 정밀도를 제공하며, 실제 기체 효과가 유효해지는 극단적 조건에서만 보정이 요구된다. 따라서 비행 로봇의 해석 영역 대부분에서는 이상 기체 가정이 이론 전개의 출발점으로 채택된다.

공기는 열역학적으로 완전 기체(perfect gas)에 가까우므로 상태 방정식

p = \rho R T, \qquad R = R_{u}/M_{\text{air}} \approx 287.05\ \mathrm{J/(kg\cdot K)}

이 높은 정확도로 성립한다. 여기서 R_{u} = 8.3145\ \mathrm{J/(mol\cdot K)}는 일반 기체 상수이다. 이 관계는 온도·압력·밀도 중 두 변수가 주어지면 나머지 하나를 결정할 수 있게 하며, 따라서 공기역학의 지배 방정식을 닫힌 체계로 구성하는 핵심 연결 관계로 작용한다. 또한 공기는 비열비 \gamma = c_{p}/c_{v} \approx 1.4를 가지며, 정적 비열 c_{v} \approx 717\ \mathrm{J/(kg\cdot K)}, 정압 비열 c_{p} \approx 1\,005\ \mathrm{J/(kg\cdot K)}으로 근사된다. 이 값들은 적당히 넓은 온도 범위에서 일정한 것으로 취급할 수 있으며, 단열 관계 p V^{\gamma} = \text{const}과 같은 공기역학의 중요한 열역학적 관계에 직접적으로 투영된다.

2. 점성, 열전도, 음속의 정량적 기술

공기의 점성은 유동장 내부의 운동량 확산을 기술하는 핵심 물성이며, 공기역학에서 Reynolds 수의 규모를 결정한다. 해수면 표준 조건에서 동점성이 아닌 절대 점성 계수는 \mu_{0} \approx 1.81 \times 10^{-5}\ \mathrm{Pa\cdot s}로 주어지고, 그 온도 의존성은 Sutherland 법칙

\mu(T) = \mu_{0}\left(\frac{T}{T_{0}}\right)^{3/2}\frac{T_{0} + S}{T + S}, \qquad T_{0} = 273.15\ \mathrm{K},\ \ S = 110.4\ \mathrm{K}

로 근사된다. 동점성 계수는 \nu = \mu/\rho로 정의되며, 공기의 경우 \nu \approx 1.48 \times 10^{-5}\ \mathrm{m^{2}/s}이다. 이 값은 경계층의 두께 \delta \sim \sqrt{\nu x/V}와 Reynolds 수 \mathrm{Re} = \rho V L / \mu의 규모를 결정하므로, 점성의 온도 의존성은 저·고고도 운용 환경 사이의 공력 성능 차이를 설명하는 물성적 근거가 된다. 실제 대기의 점성이 표준값에서 수 퍼센트 벗어날 수 있음을 고려하면, 정밀 해석에서는 Sutherland 법칙 또는 더 정교한 다항 근사를 온도에 대한 함수로 사용해야 한다.

열전도 계수 k는 열유속 \mathbf{q} = -k\nabla T를 통하여 열 전달을 기술하며, 공기의 경우 k \approx 0.0257\ \mathrm{W/(m\cdot K)}의 값을 가진다. 운동량 확산과 열확산의 상대적 크기는 Prandtl 수 \mathrm{Pr} = \mu c_{p}/k \approx 0.71로 표현되며, 공기의 Prandtl 수가 거의 1에 가깝다는 점은 속도와 온도 경계층의 두께가 유사하게 전개된다는 중요한 특성을 함축한다. 공기 중에서의 음속은 완전 기체 근사 아래 a = \sqrt{\gamma R T}로 주어지며, 해수면 표준 조건에서 a_{0} \approx 340.29\ \mathrm{m/s}이다. 음속은 Mach 수 \mathrm{Ma} = V/a를 정의하여 압축성의 정도를 판정하고, 비행 로봇의 일반적 작동 영역인 \mathrm{Ma} < 0.3에서는 밀도 변화가 작아 비압축성 가정이 정당화된다.

표 21.3.1은 비행 로봇 해석에서 자주 사용되는 공기의 기본 물성을 정리한다.

물성기호해수면 표준값단위
밀도\rho1.225kg/m³
압력p101 325Pa
온도T288.15K
절대 점성\mu1.81 \times 10^{-5}Pa·s
동점성\nu1.48 \times 10^{-5}m²/s
열전도k0.0257W/(m·K)
정압 비열c_{p}1 005J/(kg·K)
비열비\gamma1.4
음속a340.29m/s

3. 저속 근사와 압축성의 도입

비행 로봇의 일반적 작동 영역에서는 저속 비압축성 근사가 널리 사용된다. 이 근사는 Mach 수가 작을 때 밀도 변화 \Delta \rho / \rho \sim \mathrm{Ma}^{2}가 미소하다는 섭동적 논의에 기반하며, \mathrm{Ma} < 0.3에서 밀도 변화가 5% 이내로 유지됨을 근거로 한다. 이러한 조건에서는 연속 방정식이 \nabla\cdot\mathbf{u} = 0으로 단순화되고, 운동량 방정식은 \rho\, D\mathbf{u}/Dt = -\nabla p + \mu \nabla^{2}\mathbf{u} + \rho \mathbf{g}의 비압축성 Navier–Stokes 방정식으로 축약된다. 그 결과 에너지 방정식이 유동장에 대한 역결합을 상실하여 해석이 상당히 단순해지며, 베르누이 정리와 같은 간결한 관계가 유효한 범위를 확보한다. 이러한 단순화는 저속 영역의 설계와 해석에 있어 이론과 수치 도구 모두에서 결정적 이점을 제공한다.

그러나 프로펠러 팁 속도가 음속에 근접하는 고RPM 소형 로터나, 고고도에서 음속이 감소하는 상황에서는 압축성이 무시할 수 없는 수준으로 작용할 수 있다. 이 경우 밀도의 공간적 변화를 유지하는 압축성 지배 방정식이 사용되며, Karman–Tsien 보정이나 Prandtl–Glauert 보정 C_{p} = C_{p,0}/\sqrt{1 - \mathrm{Ma}^{2}}과 같은 근사 보정이 저속 결과를 준(準)압축성 영역으로 확장하는 도구로 활용된다. 충격파가 발생하는 초·초음속 영역은 비행 로봇의 표준 운용 영역을 벗어나지만, 회전익의 국부 전이음속 영역에서는 여전히 실무적 관심사가 된다. 따라서 실제 해석에서는 Mach 수의 크기에 따라 비압축성 근사, 준압축성 보정, 완전 압축성 해석 사이의 선택이 체계적으로 이루어져야 한다.

4. 공기 성질의 로봇공학적 적용

공기의 물리적 성질은 비행 로봇의 성능과 제어 모델을 결정하는 매개변수로 전면적으로 활용된다. 예를 들어 양력 L = \tfrac{1}{2}\rho V^{2} S C_{L}과 추력 T = C_{T}\rho \pi R^{2}(\Omega R)^{2}은 모두 밀도에 선형 의존하므로, 고도 증가로 인한 밀도 감소는 동일한 형태의 기체에서 직접적인 성능 저하를 초래한다. 온도와 압력이 별도로 계측되는 소형 드론의 경우, 상태 방정식 p = \rho R T을 이용하여 밀도를 실시간으로 재구성하고 이를 제어기의 공력 게인 스케줄링에 반영하는 방식이 일반화되어 있다. 또한 Reynolds 수의 변화는 익형 선택과 프로펠러 설계에 직접 개입하므로, 점성과 밀도의 고도 의존성은 저고도·고고도 임무 사이의 기체 설계 차이를 설명하는 핵심적 근거로 작동한다.

기체의 비열비와 음속 또한 로봇공학 실무에서 간접적이지만 필수적인 입력으로 쓰인다. Mach 수 계측은 피토 튜브의 정체 압력과 음속의 온도 의존성을 결합하여 수행되며, 자율 비행체의 대기 데이터 시스템은 이러한 관계를 내부적으로 사용하여 지시 대기 속도, 실제 공기 속도, 마하 수를 분리하여 출력한다. 고정익 UAV의 비행 포락선은 임계 마하 수로 상한이 결정될 수 있으며, 회전익 로봇에서는 팁 마하 수가 소음과 공탄성 문제의 임계 지표로 기능한다. 이러한 사례들은 공기의 물리적 성질이 단순히 열역학적 상수로만 기능하는 것이 아니라, 설계·계측·제어·임무 계획의 각 층에서 정량적으로 개입하는 공학적 매개변수임을 보여 준다.

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.
  • White, F. M., Viscous Fluid Flow, 3rd ed., McGraw-Hill, 2006.
  • Incropera, F. P., DeWitt, D. P., Bergman, T. L., and Lavine, A. S., Fundamentals of Heat and Mass Transfer, 7th ed., Wiley, 2011.
  • U.S. Standard Atmosphere, 1976, NOAA/NASA/USAF, Washington, D.C., 1976.
  • Moran, M. J., Shapiro, H. N., Boettner, D. D., and Bailey, M. B., Fundamentals of Engineering Thermodynamics, 8th ed., Wiley, 2014.

6. 버전

v1.0