21.5 비압축성 유동과 압축성 유동의 구분

유동을 해석할 때 가장 근본적이고 결정적인 분류 중 하나는 비압축성 유동과 압축성 유동의 구분이다. 이 구분은 밀도의 변화 정도에 의해 결정되며, 해석에 사용되는 지배 방정식, 수학적 기법, 물리적 해석이 모두 이 분류에 따라 달라진다. 저속 드론과 소형 무인기는 일반적으로 비압축성 영역에서 운용되며, 해석이 크게 단순화된다. 그러나 고정익 UAV의 고속 운용이나 프로펠러 블레이드 끝단에서는 압축성 효과가 나타날 수 있으며, 이러한 영역을 정확히 식별하고 적절한 해석 방법을 선택하는 것이 중요하다. 본 절에서는 비압축성과 압축성의 구분 기준, Mach 수의 역할, 두 해석 접근법의 차이, 천음속과 초음속 영역의 특성, 로봇 비행체에서의 적용을 체계적으로 정리한다.

1. 압축성의 물리적 정의

유체의 압축성은 압력 변화에 대한 밀도 변화의 정도를 나타내는 물리적 성질이다. 정량적으로 등엔트로피 과정에서 정의되는 압축성 계수는 다음과 같이 표현된다.

\tau = \dfrac{1}{\rho} \left(\dfrac{\partial \rho}{\partial p}\right)_s

모든 실제 유체는 유한한 압축성 계수를 가지며, 엄격한 의미에서 비압축성 유체는 존재하지 않는다. 그러나 유동 해석의 관점에서 “비압축성“은 유동의 관심 영역에서 밀도 변화가 무시할 만하여 밀도를 상수로 취급할 수 있는 경우를 의미한다. 이는 유체 자체의 성질이 아니라 유동의 특성과 관련된 실용적 분류이다.

21.5.2 Mach 수와 압축성의 기준

압축성의 정도를 정량화하는 가장 일반적인 지표는 Mach 수이며, 다음과 같이 정의된다.

M = \dfrac{V}{a}

여기서 V는 유동 속도, a는 국소 음속이다. Mach 수가 낮을수록 밀도 변화가 작으며, 높을수록 압축성 효과가 두드러진다. 경험적으로 M < 0.3일 때 밀도 변화가 5% 이내로 유지되므로, 비압축성 가정이 공학적 관점에서 충분히 정확하다. 0.3 < M < 0.8은 아음속 압축성 영역이며, 밀도 변화가 중요해진다. 0.8 < M < 1.2는 천음속 영역이며, 충격파와 국부 초음속 영역이 공존한다. 1.2 < M < 5는 초음속 영역이며, M > 5는 극초음속 영역이다.

2. 비압축성 유동의 지배 방정식

비압축성 유동의 지배 방정식은 연속 방정식과 운동량 방정식(Navier-Stokes 방정식)이며, 밀도가 상수이므로 크게 단순화된다. 연속 방정식은 다음과 같이 속도 장의 발산이 0이 되는 조건으로 축약된다.

\nabla \cdot \mathbf{V} = 0

비압축성 Navier-Stokes 방정식은 다음과 같다.

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

밀도가 상수이므로 에너지 방정식은 독립적으로 해결될 수 있으며, 일반적인 저속 유동 해석에서는 에너지 방정식이 배제된다. 이러한 단순화는 해석의 효율과 직관성을 크게 향상시키며, 저속 공기역학의 이론적 틀의 기반이 된다.

3. 압축성 유동의 지배 방정식

압축성 유동의 지배 방정식은 연속 방정식, 운동량 방정식, 에너지 방정식, 상태 방정식으로 구성된 완전한 시스템이다. 밀도가 유동의 변수가 되므로, 에너지 방정식이 독립적으로 해결될 수 없고 다른 방정식들과 결합하여 해석되어야 한다. 연속 방정식은 다음과 같다.

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

이 방정식은 밀도와 속도의 곱에 대한 보존 법칙이며, 비압축성의 경우와 구조적으로 다르다. 압축성 유동의 해석은 비선형성과 방정식의 결합으로 인해 해석이 훨씬 복잡하며, 특수한 수치 기법과 이론적 도구가 필요하다.

21.5.5 비압축성 가정의 정당화

저속 유동에서 비압축성 가정이 정당화되는 이유는 Mach 수가 낮을 때 등엔트로피 관계에서 밀도 변화가 작기 때문이다. 등엔트로피 관계로부터 유도되는 밀도비의 근사식은 다음과 같다.

\dfrac{\rho}{\rho_0} = \left(1 + \dfrac{\gamma - 1}{2} M^2\right)^{-1/(\gamma-1)}

M = 0.3에서 \gamma = 1.4의 경우 이 비율은 약 0.956이며, 밀도 변화가 약 4.4%이다. 이 수준의 오차는 대부분의 공학 응용에서 허용 가능하다. 더 높은 Mach 수에서는 이 오차가 빠르게 증가하여 비압축성 가정의 정확도가 저하된다. 따라서 Mach 수 0.3이 비압축성 가정의 실용적 상한으로 널리 인정된다.

4. 아음속 압축성 효과와 Prandtl-Glauert 보정

아음속 영역(0.3 < M < 0.8)에서는 압축성 효과가 중요해지지만, 충격파는 아직 형성되지 않는다. 이 영역에서는 비압축성 결과를 압축성 조건으로 보정하는 Prandtl-Glauert 보정이 사용된다. 예를 들어 얇은 익형의 압력 계수는 다음과 같이 보정된다.

C_p = \dfrac{C_{p,0}}{\sqrt{1 - M_\infty^2}}

여기서 C_{p,0}은 비압축성 해석의 압력 계수, M_\infty는 자유 흐름 Mach 수이다. 이 보정은 아음속 영역에서 합리적인 정확도를 제공하며, 고속 고정익 UAV의 초기 설계에 사용될 수 있다. 그러나 Mach 수가 1에 근접할수록 이 보정의 유효성은 저하된다.

21.5.7 천음속 영역과 충격파

천음속 영역(0.8 < M < 1.2)에서는 자유 흐름이 아음속이지만 익형 표면의 일부에서 국부적으로 초음속 영역이 형성되며, 아음속 영역으로 복귀할 때 충격파가 발생한다. 이 충격파는 에너지 손실, 경계층 박리, 파동 항력(wave drag)을 유발하며, 천음속 영역의 항력이 급격히 증가하는 원인이다. 천음속 영역의 해석은 비선형 편미분 방정식의 해결을 요구하며, 해석적 방법이 제한적이고 수치 해석에 크게 의존한다. 대부분의 드론과 무인기는 이 영역에 진입하지 않지만, 일부 고성능 고정익 UAV와 프로펠러 블레이드 끝단에서 국부적으로 천음속 조건이 나타날 수 있다.

21.5.8 초음속과 극초음속 영역

초음속 영역(M > 1)에서는 모든 유동 교란이 Mach 원추 내에서만 전파되며, 충격파와 팽창파가 지배적 현상이 된다. 초음속 공기역학은 항공기, 미사일, 재돌입체 등의 설계에 적용되며, 특수한 이론과 수치 기법이 사용된다. 극초음속 영역(M > 5)에서는 공력 가열, 실제 기체 효과, 해리, 화학 반응 등이 추가적으로 고려되어야 하며, 가장 복잡한 영역이다. 일반적 로봇 비행체는 이러한 영역에서 운용되지 않지만, 로켓 추진 드론과 실험적 고속 비행체의 연구에서는 고려될 수 있다.

21.5.9 저속 드론과 비압축성 해석

대부분의 드론과 소형 무인기는 최고 속도가 수십 m/s 수준이며, 이는 Mach 수가 약 0.1 이하에 해당한다. 이 영역에서는 비압축성 가정이 극히 정확하며, 공기역학 해석은 크게 단순화된다. 양력과 항력은 동적 압력에 기반한 계수 형태로 표현되고, 밀도는 대기 조건에 따른 상수로 취급된다. 비압축성 Navier-Stokes 방정식과 Euler 방정식, 포텐셜 유동 이론이 해석의 기본 도구이며, 이러한 저속 공기역학의 틀은 소형 드론 설계의 근간이다.

21.5.10 프로펠러 끝단의 압축성

프로펠러와 로터 블레이드의 끝단에서는 회전 속도와 비행 속도의 합성으로 인해 국소적으로 높은 속도가 발생할 수 있다. 블레이드 끝단 Mach 수(tip Mach number)가 약 0.7을 초과하면 국부적으로 천음속 효과가 나타나며, 0.9 이상에서는 충격파가 형성되어 효율 저하와 소음 증가를 유발한다. 드론의 소형 프로펠러는 일반적으로 끝단 Mach 수가 낮게 유지되지만, 고성능 헬리콥터와 틸트로터의 경우 이 한계가 설계의 중요한 제약이다. 블레이드 끝단의 압축성 효과를 고려한 블레이드 형상 설계와 회전수 제한이 실무적 대응 방안이다.

21.5.11 성층권 무인기와 밀도 효과

성층권 무인기(HAPS, high-altitude platform stations)는 18 km 이상의 고도에서 장기 체공하며, 통신 중계와 지구 관측을 수행한다. 이 영역의 공기 밀도는 해수면의 약 7% 수준이며, 음속도 낮아 동일한 비행 속도에서 Mach 수가 증가한다. 밀도가 낮으므로 양력과 항력이 감소하고, 레이놀즈 수도 감소한다. 비압축성 가정은 여전히 유효할 수 있지만, 저밀도·저레이놀즈 수 영역의 고유한 공력 특성이 나타나며, 특수한 익형 설계와 넓은 날개 면적이 요구된다.

21.5.12 해석 방법의 선택 기준

유동 해석의 첫 단계는 유동의 성격을 평가하고 적절한 해석 방법을 선택하는 것이다. Mach 수의 평가는 압축성 고려 여부를 결정하며, 레이놀즈 수는 점성 효과와 층류·난류 특성을 결정한다. 이들 지표에 기반하여 해석 방법이 선택된다. 저 Mach 저 Reynolds 수 영역에서는 비압축성 포텐셜 유동 이론이나 저레이놀즈 수 RANS(Reynolds-averaged Navier-Stokes) 해석이 사용된다. 중간 Mach 수에서는 압축성 RANS 해석이 일반적이며, 높은 Reynolds 수에서는 난류 모델의 선택이 중요하다. 천음속과 초음속 영역에서는 특수한 수치 기법이 필요하며, LES(large eddy simulation)와 DNS(direct numerical simulation)가 고정밀 해석에 사용된다.

21.5.13 로봇 공학에서의 적용

로봇 비행체의 공기역학 해석에서 비압축성/압축성의 구분은 다음과 같이 반영된다. 첫째, 일반 드론과 멀티로터는 비압축성 가정 하에서 해석되며, 양력과 항력은 동적 압력과 무차원 계수로 표현된다. 둘째, 고속 고정익 UAV는 자유 흐름 Mach 수가 0.3을 초과할 수 있으며, 이 경우 Prandtl-Glauert 보정 또는 압축성 해석이 필요하다. 셋째, 헬리콥터와 틸트로터의 블레이드 끝단 압축성은 설계와 성능 해석에서 중요한 고려 사항이다. 넷째, 성층권 무인기와 고고도 운용 로봇은 저밀도 환경에서 특수한 공기역학적 특성을 가진다. 다섯째, 초음속 드론과 극한 환경 비행체는 압축성 효과가 설계의 중심이며, 전산 유체역학의 광범위한 활용이 필수이다.

21.5.14 요약과 후속 연결

비압축성 유동과 압축성 유동의 구분은 Mach 수에 기반하며, M < 0.3이 비압축성 가정의 실용적 상한이다. 비압축성 유동은 연속 방정식과 운동량 방정식으로 축약되어 해석이 크게 단순화되며, 저속 드론과 소형 무인기의 해석에 적합하다. 압축성 유동은 에너지 방정식과 상태 방정식을 포함한 완전한 시스템을 요구하며, 고속 비행체와 블레이드 끝단 등에서 중요하다. 천음속, 초음속, 극초음속 영역은 각각 고유한 특성과 해석 기법을 요구한다. 다음 절에서는 연속 방정식의 공기역학적 적용을 다루어, 질량 보존이 공기역학 해석의 기본 방정식으로 어떻게 활용되는가를 체계적으로 이해하는 관점을 제공한다.

출처

• Anderson, J. D., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2016.
• Anderson, J. D., Modern Compressible Flow: With Historical Perspective, 3rd ed., McGraw-Hill, 2003.
• Liepmann, H. W., and Roshko, A., Elements of Gasdynamics, Dover, 2001.
• Shapiro, A. H., The Dynamics and Thermodynamics of Compressible Fluid Flow, Vol. 1 and 2, Ronald Press, 1953.
• Ashley, H., and Landahl, M., Aerodynamics of Wings and Bodies, Dover, 1985.
• Hirsch, C., Numerical Computation of Internal and External Flows, 2nd ed., Butterworth-Heinemann, 2007.

버전

v1.0