공력해석

전산유체역학 유체역학 마지막 수정: 2026-07-06 04:37:09

1. 개요[편집]

새는 그냥 난다. 인간은 계수를 열두 개 뽑고도 왜 나는지 매번 헷갈린다.

공력해석(aerodynamic analysis)은 공기 속을 움직이는 물체에 작용하는 힘과 모멘트를 예측하고 분석하는 유체역학의 한 분야로, 날개·차체·건물 등에 걸리는 양력(lift)·항력(drag)·모멘트(moment)를 정량화하는 것을 목표로 한다. 풍동 실험과 전산유체역학이 양대 도구이며, 항공기 설계에서 자동차 연비, 풍하중 평가까지 응용이 광범위하다.

공력해석의 언어는 힘 그 자체가 아니라 무차원 계수다. 절대적인 힘은 속도·밀도·크기에 따라 제멋대로 변하지만, 이를 동압으로 정규화한 계수는 형상만 같으면 스케일이 달라도 비교 가능하기 때문이다. 이 무차원화 덕분에 손바닥만 한 풍동 모형의 데이터로 실물 항공기를 논할 수 있다. 무차원수의 축복이다.

2. 공력 3계수: 양력·항력·모멘트[편집]

물체에 걸리는 힘과 모멘트는 동압 $\frac{1}{2}\rho V^2$ 과 기준 면적 $S$ 로 나눠 무차원화한다.

C_L = \frac{L}{\tfrac{1}{2}\rho V^2 S}, \qquad C_D = \frac{D}{\tfrac{1}{2}\rho V^2 S}, \qquad C_M = \frac{M}{\tfrac{1}{2}\rho V^2 S c}

여기서 $L$ 은 흐름에 수직인 양력, $D$ 는 흐름 방향의 항력, $M$ 은 피칭 모멘트, $c$ 는 기준 시위길이(chord)다. 세 계수는 받음각(angle of attack) $\alpha$ 의 함수로 표현되며, 공력 특성의 전부라 해도 과언이 아니다.

양력계수 $C_L$ : 받음각에 대해 거의 선형으로 증가하다가, 실속각에서 급락한다. 얇은 날개 이론에 따르면 저받음각에서 $C_L \approx 2\pi\alpha$ 의 기울기를 갖는다.
항력계수 $C_D$ : 낮을수록 좋다. 연비와 최고속도를 직접 좌우한다.
모멘트계수 $C_M$ : 기체의 세로 안정성(pitch stability)을 결정한다. 무게중심과 공력중심의 관계에서 부호가 갈린다.

3. 항력의 두 얼굴: 압력항력 vs 마찰항력[편집]

항력은 발생 메커니즘에 따라 크게 둘로 쪼갤 수 있다.

마찰항력(friction drag): 물체 표면을 스치는 공기의 점성 전단응력이 누적된 것. 경계층 내부의 속도 구배가 벽면에 만드는 마찰이 원인이다. 표면이 매끈하고 흐름이 층류일수록 작다.
압력항력(pressure drag / form drag): 물체 앞뒤의 압력차가 만드는 항력. 특히 경계층이 물체 뒤에서 박리(separation)되어 후류(wake)가 넓어지면, 뒷면 압력이 회복되지 못해 압력항력이 폭증한다. 골프공에 딤플을 파는 이유가 바로 이 박리를 늦춰 압력항력을 줄이려는 것.¹

유선형(streamlined) 물체는 압력항력이 작아 마찰항력이 지배적이고, 뭉툭한(bluff) 물체는 박리 때문에 압력항력이 압도적이다. 원기둥 뒤에 생기는 카르만 와열이 대표적인 압력항력의 원흉. 둘의 비율은 형상과 레이놀즈수에 따라 극적으로 달라진다.

4. 유도항력과 극곡선[편집]

3차원 유한 날개에는 2차원 익형에는 없는 항력이 추가로 붙는다. 바로 유도항력(induced drag)이다. 날개 끝에서 아랫면의 고압 공기가 윗면으로 말려 올라가며 날개끝 와류(wingtip vortex)를 만드는데, 이 와류가 하강기류(downwash)를 유도해 유효 받음각을 깎고 양력 벡터를 뒤로 기울인다. 그 뒤로 기운 성분이 곧 유도항력이다.²

유도항력계수는 양력계수의 제곱에 비례하는 것이 핵심 특성이다.

C_{D,i} = \frac{C_L^2}{\pi e\, AR}

여기서 $AR$ 은 가로세로비(aspect ratio), $e$ 는 오스왈드 효율계수다. 글라이더의 날개가 길고 가느다란 이유가 바로 이것 — $AR$ 을 키워 유도항력을 줄이려는 것. 전체 항력은 양력과 무관한 유해항력(parasite drag, 마찰+압력)과 양력 제곱에 비례하는 유도항력의 합이다.

C_D = C_{D,0} + \frac{C_L^2}{\pi e\, AR}

이 관계를 $C_L$ - $C_D$ 평면에 그린 것이 극곡선(drag polar)이다. 항공기 성능 해석의 심장부로, 곡선 위 특정 점들이 최대 양항비(순항 효율 최적점), 최소 항력 속도 등을 알려준다. 원점에서 극곡선에 그은 접선의 접점이 최대 양항비 $(L/D)_{max}$ 지점이라는 것은 공력 엔지니어의 상식.

5. 실속: 양력이 무너지는 순간[편집]

받음각을 계속 키우면 어느 순간 경계층이 날개 윗면의 역압력구배(adverse pressure gradient)를 못 버티고 대규모로 박리한다. 이 순간 양력이 급락하고 항력이 폭증하는 현상이 실속(stall)이다. 실속각은 익형과 레이놀즈수에 따라 대략 12~18도 부근.

실속은 CFD로 예측하기 가장 까다로운 현상 중 하나다. 박리 지점과 후류의 비정상성을 제대로 잡으려면 난류 모델링이 정교해야 하는데, 흔히 쓰는 RANS 모델은 박리 예측에서 자주 삐끗한다. Spalart-Allmaras나 SST $k$ - $\omega$ 모델이 그나마 박리에 강하다고 알려져 있으나, 완전 실속 이후의 대규모 박리 영역에서는 LES나 DES 같은 비정상 해석이 필요해진다. “실속각 예측 오차 2도”가 논문 한 편이 되는 세계다.³

6. 압축성과 고속 공력[편집]

비행 속도가 음속에 가까워지면 공기를 더는 비압축성으로 못 본다. 마하수 $M = V/a$ 가 0.3을 넘으면 밀도 변화가 무시할 수 없어지고, 천음속(transonic) 영역에서는 날개 윗면에 국소 초음속 영역과 충격파가 생겨 조파항력(wave drag)이 추가된다. 이 항력 발산(drag divergence)이 여객기 순항 마하수를 0.85 언저리에 묶어두는 물리적 족쇄다. 고속 공력해석은 압축성 유동 지배 방정식과 충격파 포착에 강한 차분 도식을 요구한다.

7. 여담: 계수 하나에 밥줄이 걸린다[편집]

자동차 회사에서 $C_D$ 를 0.01 낮추면 부서 회식을 한다는 말이 있다. 항공기에서는 유해항력 카운트(drag count, $C_D$ 0.0001) 하나가 수억 원의 연료비로 환산된다. 그래서 공력 엔지니어는 소수점 넷째 자리를 두고 격자를 갈아넣고 난류 모델링 상수를 만지작거리며 밤을 새운다. 새는 그냥 나는데 말이다.

8. 관련 문서[편집]

딤플이 층류 경계층을 강제로 난류로 전이시킨다. 난류 경계층은 운동량이 커서 역압력구배를 더 오래 버티고, 박리점이 뒤로 밀린다. 그 결과 후류가 좁아져 압력항력이 준다. 골프공이 매끈했다면 비거리가 절반으로 줄었을 것이다. ↩
그래서 이론상 무한히 긴 날개(2차원 익형)에는 유도항력이 없다. 현실의 날개는 반드시 끝이 있으므로 유도항력을 피할 수 없다. 세상에 공짜 양력은 없다. ↩
최대양력계수 $C_{L,max}$ 와 실속각 예측은 CFD 검증 워크숍(AIAA High Lift Prediction Workshop)의 단골 주제다. 매번 참가 코드들의 결과가 실속각에서 사이좋게 흩어진다. 실속은 겸손을 가르치는 현상이다. ↩