21.16 항력 계수와 항력 분석 (Drag Coefficient and Drag Analysis)

1. 항력 계수의 정의

항력 계수(drag coefficient, $C_D$ )는 물체에 작용하는 항력(drag)을 무차원화한 공기역학적 계수로서, 다음과 같이 정��된다:

$C_D = \frac{D}{q_\infty S} = \frac{D}{\frac{1}{2}\rho_\infty V_\infty^2 S}$

여기서 $D$ 는 항력, $q_\infty$ �� 자유류 동압(dynamic pressure), $S$ 는 기준 면적(reference area)이다. 2차원 익형에서는 단위 스팬당 항력 계수 $C_d$ 를 사용하며, 기준 길이로 시위 길이 $c$ 를 취한다:

$C_d = \frac{D'}{q_\infty c}$

기준 면적의 선택은 물체의 유형에 따라 상이하다. 날개와 익형에서는 평면형 면적(planform area)을 사용하며, 동체와 같은 회전체에서는 전면 투영 면적(frontal area) 또는 습윤 면적(wetted area)이 사용되기도 한다. 항력 계수 값을 비교할 때는 기준 면적의 정의를 반드시 확인��여야 한다.

2. 항력 극곡선

항력 극곡선(drag polar)은 양력 계수( $C_L$ )와 항력 계수( $C_D$ )의 관계를 나타내는 곡선으로서, 항공기 및 비행체의 공력 성능을 종합적으로 표현하는 기본적 도표이다. 가장 널리 사용되는 포물선 근사(parabolic approximation)는 다음과 ��다:

$C_D = C_{D,0} + K C_L^2$

여기�� $C_{D,0}$ 는 영양력 항력 계수(zero-lift drag coefficient), $K$ 는 양력 유발 항력 인자(lift-induced drag factor)로서 $K = 1/(\pi e AR)$ 이다. $e$ 는 오스왈드 효율 인자(Oswald efficiency factor), $AR$ 은 종횡비(aspect ratio)이다.

오스왈드 효율 인자는 타원 양력 분포로부터의 편차와 양력에 따른 프로파일 항력 증가를 모두 포함하는 경험적 인자로서, 일반적인 항공기에서 $e \approx 0.7 \sim 0.85$ 의 값을 갖는다. $e = 1$ 은 타원 양력 분포를 갖는 이상적 날개에 해당한다.

보다 정밀한 항력 극곡선 모델에서는 최소 항력이 영양력 조건이 아닌 특정 양력 계수 $C_{L,\text{min}D}$ 에서 발생하는 경우를 반영한다:

$C_D = C_{D,\text{min}} + K(C_L - C_{L,\text{min}D})^2$

3. 2차원 익형의 항력 계수 특성

2차원 ��형의 항력 계수(프로파일 ��력 계수, $C_d$ )는 마찰 항력과 형상 항력의 합으로 구성된다. 유도 항력은 3차원 효과이므로 2차원 익형에는 존재하지 않는다.

층류 익형(laminar airfoil)의 항력 곡선에서 특징적인 현상은 항력 버킷(drag bucket)의 존재이다. 항력 버킷은 설계 양력 계수 근처의 좁은 양력 계수 범위에서 항력 계수가 현저히 낮게 유지되는 구간이다. 이 구간에서는 상면과 하면 모두에서 순압력 구배(favorable pressure gradient)가 유지되어 층류 경계층이 넓은 범위에서 보존되며, 마찰 항력이 최소화된다. 항력 버킷의 범위를 벗어나면 한쪽 표면에서 역압력 구배에 의한 조기 천이가 발생하여 항력이 급증한다 (Abbott & von Doenhoff, 1959).

저레이놀즈 수( $Re < 5 \times 10^5$ ) 영역에서는 층류 박리 거품(laminar separation bubble)의 존재가 프로파일 항력에 현저한 영향을 미친다. 층류 박리 거품에 의한 추가적인 압력 항력은 저레이놀즈 수 익형의 항력 계수를 고레이놀즈 수 대비 현저히 증가시��다.

4. 마찰 항력 계수의 해석

평판(flat plate)에 대한 마찰 항력 계수는 경계층 이론으로부터 해석적 또는 반경험적으로 구할 수 있으며, 이는 익형 및 물체의 마찰 항력 추정에서 기준값으로 활용된다.

층류 경계층: 블라지우스 해(Blasius solution)에 의한 평판 전체의 층류 마찰 항력 계수는:

$C_f = \frac{1.328}{\sqrt{Re_L}}$

여기서 $Re_L = \rho_\infty V_\infty L / \mu$ 는 평판 길이 $L$ 에 기반�� 레이놀즈 ��이다.

난류 경��층: 슈리히팅(Schlichting)의 경험식에 의한 평판 전체의 난류 마찰 항력 계수는:

$C_f = \frac{0.455}{(\log_{10} Re_L)^{2.58}}$

혼합 경계층(mixed boundary layer): 실제 유동에서는 앞전으로부터 일정 거리까지 층류 경계층이 유지된 후 난류로 천이하는 경우가 대부분이다. 천이 위치 $x_\text{tr}$ 에서의 레이놀즈 ��를 $Re_\text{tr}$ 이라 하면, 혼합 경계층의 마찰 항력 계수는:

$C_f = \frac{0.455}{(\log_{10} Re_L)^{2.58}} - \frac{A}{Re_L}$

여기서 상수 $A$ 는 천이 레이놀즈 수에 의존하며, $Re_\text{tr} = 5 \times 10^5$ 일 때 $A \approx 1700$ 이다 (Schlichting & Gersten, 2017).

21.16.5 형상 항력 계수와 등가 평판 면적

항공기 전체의 항력 분석에서 등가 평판 면적(equivalent flat plate area, $f$ )은 각 구성 요소의 항력 기여를 직관적으로 비교하기 위한 지표이다:

$f = C_D S$

등가 평판 면적은 해당 물체와 동일한 항력을 발생시키는 평판의 면적으로서, 항력의 절대적 크기를 단일 스칼라 값으로 나타낸다.

항공기 전체의 영양력 항력 계수는 각 구성 요소의 항력 기여를 합산하여 추정할 수 있다:

$C_{D,0} = \frac{1}{S_\text{ref}} \sum_i C_{D,i} S_i + \Delta C_{D,\text{interference}} + \Delta C_{D,\text{misc}}$

여기서 $C_{D,i}$ 와 $S_i$ 는 각 구성 요소( $i$ )�� 항력 계수와 기준 면적, $\Delta C_{D,\text{interference}}$ 는 간섭 항력 증분, $\Delta C_{D,\text{misc}}$ 는 기타 항력 증분(리벳, 틈새, 안테나 등)이다.

21.16.6 둔두 물체�� 항력 계수

유선형이 아닌 둔두 물체(bluff body)�� 항력 계수는 유동 박리의 특성에 의해 결정되며, 형상 항력이 지배적이다. 대표적인 둔두 물체의 항력 계수(전면 투영 면적 기준)�� 다음과 같다:

무한 원통(2차원, $Re > 10^4$ ): $C_D \approx 1.2$
구(sphere, $Re > 10^4$ ): $C_D \approx 0.4 \sim 0.5$
평판(유동에 수직, 무한 스팬): $C_D \approx 2.0$
정사각형 단면(2차원): $C_D \approx 2.1$

둔두 물체의 항력 계수는 레이놀즈 수에 대해 비단조적(non-monotonic) 의존성을 보일 수 있다. 특히 구와 원통에서 임계 레이놀즈 수(critical Reynolds number) 부근에서 항력 계수가 급격히 감소하는 항력 위기(drag crisis) 현상이 관찰된다. ��는 경계층의 층류-난류 천이에 의해 박리점이 후방으로 이동하고 후류(wake)의 폭이 축소되기 때문이다 (Anderson, 2017).

21.16.7 항력의 실험적 측정

항력의 실험적 측정은 풍동(wind tunnel)에서 다음의 방법으로 수행된다:

힘 측정법(force balance method): 풍동 저울(balance)을 이용하여 물체에 작용하는 항력을 직접 측정한다. 3분력 또는 6분력 저울이 사용되며, 양력, 항력, 모멘트를 동시에 측정할 수 있다. 이 방법은 전체 항력을 직접 제공하지만, 항력의 성분별 분해에는 추가적 분석이 필요하다.

후류 측량법(wake survey method): 물체 후방의 후류에서 전압(total pressure)과 정압(static pressure)의 분포를 측정하고, 운동량 적분(momentum integral)을 통해 항력을 산출한다. 이 방법은 2차원 익형의 프로파일 항력 측정에 널리 활용되며, 후류에서의 운동량 결손(momentum deficit)으로부터 항력을 정밀하게 결정할 수 있��. Jones(1936)가 제시한 후류 측량 공식은 다음과 같다:

$C_d = \frac{2}{c} \int_{-\infty}^{\infty} \frac{u}{V_\infty}\left(1 - \frac{u}{V_\infty}\right) dy$

여기서 $u$ 는 후류에서의 유동 방향 속도이다.

표면 압력 적분법: 물체 표면의 압력 분포를 측정하고 적분하여 압력 항력을 산출한다. 마찰 항력은 이 방법으로 직접 측정되지 않으므로, 전체 항력에서 압력 항력을 차감하여 간접적으로 구한다.

5. 항력 분석의 수치적 방법

전산 유체역학(CFD)에서 항력 계수의 정확한 예측은 도전적인 과제이며, 수치 해석의 정밀도가 중요하다. 항력은 양력에 비해 절대적 크기가 작으므로, 수치 오차의 상대적 영향이 크다.

RANS(Reynolds-Averaged Navier-Stokes) 방법: 난류 모델(turbulence model)을 사용하여 시간 평균된 유동장을 해석한다. 마찰 항력과 형상 항력의 예측에 널리 활용되지만, 난류 모델의 선택에 따라 예측 정확도가 달라진다. Spalart-Allmaras 모델, $k$ - $\omega$ SST 모델 등이 공기역학 해석에서 표준적으로 사용된다.

원장 방법(far-field method)과 근장 방법(near-field method): 항력의 수치적 산출에는 물체 표면에서의 압력과 전단 응력을 적분하는 근장 방법과, 물체로부터 먼 거리의 유동장 변수를 이용하는 원장 방법이 있다. 원장 방법은 항력을 물리적 성분(점성 항력, 유도 항력, 파항력 ��)으로 분해할 수 있어 항력 저감 설계에 유용한 정보를 ��공한다.

점성-비점성 연성 방법: XFOIL 등의 점��-비점성 연성 해석 코드는 2차�� 익형의 프로파일 항력을 효율적으로 예측하며, 특히 저레��놀즈 수 영역에서의 층류 박리 거품과 천이 현상을 적절히 모델링한다 (Drela, 1989).

6. 항력 분할 분석

항공기의 항력 저감 설계를 위해서는 전체 항력을 물리적 발생 원천별로 분할하여 각 성분의 기여도를 정량적으로 파악하는 것이 필요하다. 항력 분할(drag breakdown) 분석의 일반적 결과는 다음과 같다:

아음속 고정익 항공기의 순항 조건에서:

마찰 항력: 전체 항력의 약 40~60%
형상 항력: 전체 항력의 약 5~15%
유도 항��: 전체 항력의 약 30~40%
간섭 항력 및 기타: 전체 항력의 약 5~10%

이 비율은 항공기의 형태, 종횡비, 순항 양력 계수에 따라 상당히 변화한다. 소형 UAV에서는 동체 비율이 크고 날개 종횡비가 상대적으로 낮은 경우가 많아, 기생 항력의 비율이 유인 항공기보다 높은 경향�� 있다.

7. 로봇 공학에서의 항력 분석 실무

비행 로봇의 설계에서 항력 분석은 요구 동력(required power)과 에너지 소비의 예측에 직접 연결되며, 비행 지속 시간과 항속 거리의 결정에 핵심적 역할을 한다.

고정익 UAV에서 수평 정상 비행 시 요구 추력은 항력과 같으므로, 요구 동력은:

$P_\text{req} = DV = \frac{1}{2}\rho V^3 S C_D$

최소 요구 동력 속도와 최대 항속 거리 속도는 항력 극곡선으로부터 결정되며, 이는 비행 로봇의 임무 계획(mission planning)에서 기본적인 성능 매개변수이다.

멀티로터 시스템에서는 호버링 시 로터 블레이드의 프로파일 항력이 요구 동력의 주요 성분 중 하나이며, 전진 비행 시에는 기체의 기생 항력에 의한 추가 요구 추력이 비행 효율을 저하시킨다. 따라서 멀티로터의 전진 비행 효율 향상을 위해서는 기체 형상의 유선형화와 기생 항력의 저감이 중요하다 (Leishman, 2006).

참고 문헌

Abbott, I. H., & von Doenhoff, A. E. (1959). Theory of Wing Sections: Including a Summary of Airfoil Data. Dover Publications.
Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
Drela, M. (1989). XFOIL: An analysis and design system for low Reynolds number airfoils. In Low Reynolds Number Aerodynamics, Lecture Notes in Engineering, Vol. 54, Springer.
Jones, B. M. (1936). The measurement of profile drag by the pitot-traverse method. ARC R&M No. 1688.
Leishman, J. G. (2006). Principles of Helicopter Aerodynamics (2nd ed.). Cambridge University Press.
Schlichting, H., & Gersten, K. (2017). Boundary-Layer Theory (9th ed.). Springer.

버전: 2026-04-16-v1.0