22.12 삼차원 날개의 양력 분포와 리프팅 라인 이론

1. 차원 익형에서 3차원 날개로의 확장

2차원 익형 이론은 무한 스팬(infinite span)을 가정한 해석이다. 실제 날개는 유한한 스팬을 가지므로, 날개 끝(wingtip)에서의 유동 현상이 전체 양력 분포에 근본적인 영향을 미친다. 유한 날개에서는 날개 끝에서 하면의 고압 영역에서 상면의 저압 영역으로 유동이 회전하여 날개 끝 와류(tip vortex)가 형성되며, 이로 인하여 날개 뒤쪽으로 와류 시트(vortex sheet)가 방출된다.

이 후류 와류 시스템은 날개 위치에서 하향세류(downwash)를 유도하여 유효 받음각(effective angle of attack)을 감소시킨다. 그 결과, 유한 날개의 양력은 동일한 기하학적 받음각에서 2차원 익형보다 낮아지며, 유도 항력(induced drag)이라는 추가적인 항력 성분이 발생한다 (Anderson, 2017).

2. 리프팅 라인 이론의 기본 개념

2.1 프란틀의 리프팅 라인 모형

리프팅 라인 이론(lifting line theory)은 Ludwig Prandtl이 1918~1919년에 제시한 유한 날개의 공력 해석 이론으로, 3차원 날개 공력학의 초석이 된다. 이 이론에서 날개는 시위 방향으로 수축되어 1/4 시위점(quarter-chord)을 연결하는 단일 속박 와류선(bound vortex line)으로 대체된다.

핵심 가정은 다음과 같다.

1. 고종횡비 날개: 날개의 종횡비(aspect ratio, AR)가 충분히 커서 날개 단면에서의 유동이 국소적으로 2차원 유동에 가까운 것으로 간주할 수 있다.
2. 비압축성, 비점성 유동: 포텐셜 유동을 가정한다.
3. 소교란 가정: 유동 교란이 자유류에 비하여 작다.

2.2 속박 와류와 자유 와류

헬름홀츠 와류 정리(Helmholtz vortex theorem)에 의하면, 와류선은 유체 내에서 자유롭게 끝날 수 없다. 따라서 속박 와류의 순환 강도 \Gamma(y)가 스팬 방향 위치 y에 따라 변하면, 순환의 변화분에 해당하는 자유 와류(free vortex)가 날개에서 후류로 방출된다.

위치 y에서 미소 구간 dy에 걸쳐 방출되는 자유 와류의 강도는 다음과 같다.

d\gamma = -\frac{d\Gamma}{dy} dy

이 자유 와류들은 날개 후방으로 무한히 뻗어나가며, 전체적으로 후류 와류 시트(trailing vortex sheet)를 형성한다.

3. 유도 하향세류

후류 와류 시트가 날개의 스팬 방향 위치 y_0에서 유도하는 하향세류 속도(downwash velocity) w(y_0)는 비오-사바르 법칙(Biot-Savart law)에 의하여 다음과 같이 계산된다.

w(y_0) = -\frac{1}{4\pi}\int_{-b/2}^{b/2} \frac{(d\Gamma/dy)}{y_0 - y} dy

여기서 b는 날개 스팬이며, 적분은 코시 주치(Cauchy principal value)의 의미에서 수행된다. 하향세류에 의한 유도 받음각(induced angle of attack)은 다음과 같다.

\alpha_i(y_0) = -\frac{w(y_0)}{U_\infty} = \frac{1}{4\pi U_\infty}\int_{-b/2}^{b/2} \frac{(d\Gamma/dy)}{y_0 - y} dy

4. 프란틀의 기본 적분 방정식

스팬 방향 위치 y_0에서의 유효 받음각은 기하학적 받음각에서 유도 받음각을 뺀 값이다.

\alpha_{\text{eff}}(y_0) = \alpha(y_0) - \alpha_i(y_0)

리프팅 라인 이론의 핵심 가정에 의하여, 각 스팬 위치에서의 국소 양력은 2차원 양력 기울기를 사용하여 다음과 같이 표현된다.

\Gamma(y_0) = \frac{1}{2} c(y_0) a_0 U_\infty \left[\alpha_{\text{eff}}(y_0) - \alpha_{L=0}(y_0)\right]

여기서 c(y_0)는 국소 시위 길이, a_0는 2차원 양력 기울기(일반적으로 2\pi), \alpha_{L=0}(y_0)는 국소 영양력 받음각이다.

위의 두 식을 결합하면 프란틀의 기본 적분 방정식(Prandtl’s fundamental equation)을 얻는다.

\Gamma(y_0) = \frac{1}{2} c(y_0) a_0 U_\infty \left[\alpha(y_0) - \alpha_{L=0}(y_0) - \frac{1}{4\pi U_\infty}\int_{-b/2}^{b/2} \frac{(d\Gamma/dy)}{y_0 - y} dy\right]

이 식은 미지 함수 \Gamma(y)에 대한 적분-대수 방정식이다.

5. 글라우어트 급수 해법

5.1 좌표 변환

기본 적분 방정식을 풀기 위하여 스팬 좌표를 각도 변수로 변환한다.

y = -\frac{b}{2}\cos\theta, \quad 0 \leq \theta \leq \pi

이 변환에 의하여 y = -b/2(좌측 날개 끝)는 \theta = 0에, y = b/2(우측 날개 끝)는 \theta = \pi에, y = 0(날개 중심)은 \theta = \pi/2에 대응한다.

순환 분포를 다음의 정현 급수(sine series)로 전개한다.

\Gamma(\theta) = 2bU_\infty \sum_{n=1}^{N} A_n \sin(n\theta)

\sin(n\theta)의 사용은 날개 끝(\theta = 0, \pi)에서 순환이 자연스럽게 영이 되도록 보장한다.

5.2 급수 계수의 결정

급수 전개를 기본 적분 방정식에 대입하면, 글라우어트의 적분 공식에 의하여 하향세류 적분이 해석적으로 평가된다.

\frac{1}{4\pi}\int_0^\pi \frac{\sum n A_n \cos(n\theta)}{\cos\theta - \cos\theta_0}\sin\theta \, d\theta = \sum_{n=1}^{N} n A_n \frac{\sin(n\theta_0)}{\sin\theta_0}

이를 대입하면 각 스팬 위치 \theta_0에서 다음의 대수 방정식을 얻는다.

\sum_{n=1}^{N} A_n \sin(n\theta_0)\left[\frac{4b}{a_0 c(\theta_0)} + \frac{n}{\sin\theta_0}\right] = \alpha(\theta_0) - \alpha_{L=0}(\theta_0)

N개의 스팬 위치 \theta_0를 선택하면 N개의 미지수 A_n에 대한 N원 연립 일차 방정식이 구성되며, 이를 풀어 순환 분포를 결정한다.

6. 타원형 양력 분포

6.1 타원 분포의 특성

순환 분포가 타원형(elliptical distribution)인 경우는 리프팅 라인 이론의 가장 중요한 특수 해이다. 타원 분포는 다음과 같다.

\Gamma(y) = \Gamma_0\sqrt{1 - \left(\frac{2y}{b}\right)^2}

여기서 \Gamma_0는 날개 중심(y = 0)에서의 최대 순환이다. 급수 표현에서는 A_1만 비영이고 나머지 A_n = 0 (n \geq 2)인 경우에 해당한다.

\Gamma(\theta) = 2bU_\infty A_1 \sin\theta

6.2 균일 하향세류

타원형 양력 분포의 가장 중요한 특성은 유도 하향세류가 스팬 전체에 걸쳐 균일(uniform)하다는 것이다.

w = -\frac{\Gamma_0}{2b} = \text{const}

유도 받음각도 스팬 전체에 걸쳐 일정하다.

\alpha_i = \frac{\Gamma_0}{2bU_\infty} = \frac{C_L}{\pi AR}

여기서 C_L은 전체 날개의 양력 계수, AR = b^2/S는 종횡비이다.

6.3 양력 계수

타원 분포를 가지는 날개의 양력 계수는 다음과 같다.

C_L = \pi AR \cdot A_1

양력 기울기는 다음과 같이 유도된다.

a = \frac{dC_L}{d\alpha} = \frac{a_0}{1 + \dfrac{a_0}{\pi AR}}

이 결과는 유한 날개의 양력 기울기가 2차원 양력 기울기 a_0보다 항상 작음을 보여준다.

6.4 타원 분포의 실현 조건

타원형 양력 분포는 다음의 조건에서 실현된다.

1. 타원형 평면형(elliptical planform): 시위 분포가 c(y) = c_0\sqrt{1 - (2y/b)^2}인 경우
2. 기하학적 비틀림이 없는 경우: 받음각과 영양력 받음각이 스팬 방향으로 일정한 경우

Supermarine Spitfire의 날개가 타원형 평면형의 대표적인 역사적 사례이다.

7. 일반적인 양력 분포

7.1 비타원 분포와 유도 항력

실제 날개는 타원형 평면형이 아닌 사다리꼴(trapezoidal), 직사각형(rectangular), 또는 테이퍼형(tapered) 평면형을 가진다. 이 경우 순환 분포에서 고차 계수 A_n (n \geq 2)이 비영이 되어 양력 분포가 타원에서 벗어난다.

비타원 분포에서의 유도 항력 계수는 다음과 같다.

C_{D_i} = \frac{C_L^2}{\pi AR}(1 + \delta)

여기서 \delta = \sum_{n=2}^{N} n(A_n/A_1)^2 \geq 0이다. 타원 분포에서 \delta = 0이므로, 타원 분포가 주어진 양력에 대하여 유도 항력을 최소화하는 최적 분포임이 증명된다.

7.2 양력 분포에 영향을 미치는 요인

스팬 방향 양력 분포의 형태는 다음의 설계 매개변수에 의하여 결정된다.

7.3 테이퍼비의 효과

사다리꼴 날개에서 테이퍼비 \lambda = c_{\text{tip}}/c_{\text{root}}는 양력 분포의 형태를 결정하는 핵심 매개변수이다. 직사각형 날개(\lambda = 1)에서는 양력이 날개 중심에 과도하게 집중되어 비효율적이며, 적절한 테이퍼(\lambda \approx 0.4)를 적용하면 근사 타원 분포를 달성할 수 있다. 과도한 테이퍼(\lambda < 0.2)에서는 날개 끝의 국소 양력 계수가 과대하여 날개 끝 실속이 발생할 위험이 있다.

8. 수치 해법

리프팅 라인 이론의 기본 적분 방정식은 글라우어트 급수 외에도 다양한 수치적 방법으로 풀 수 있다.

8.1 다수점 배치법(Multhopp’s Method)

Multhopp(1938)가 제안한 방법은 스팬 방향으로 등간격 \theta 격자점을 배치하고, 코시 주치 적분을 이산화하여 연립 방정식을 구성한다. 이 방법은 수치적 정확도가 높고 구현이 비교적 용이하여 실무에서 널리 사용된다.

8.2 와류 격자법(Vortex Lattice Method)

와류 격자법은 날개를 시위 방향과 스팬 방향으로 분할된 패널로 이산화하고, 각 패널에 말굽 와류(horseshoe vortex)를 배치하여 유동을 모형화한다. 이 방법은 리프팅 라인 이론의 확장으로서, 후퇴각, 테이퍼, 비틀림 등 복잡한 날개 형상을 용이하게 처리할 수 있다 (Katz & Plotkin, 2001).

9. 드론 설계에의 적용

리프팅 라인 이론은 고정익 드론의 주익 설계에서 다음과 같이 활용된다.

1. 초기 날개 형상 설정: 요구 양력에 대하여 최적의 스팬, 종횡비, 테이퍼비를 결정한다.
2. 비틀림 최적화: 양력 분포를 타원에 가깝게 조정하기 위한 기하학적 또는 공력적 비틀림을 설계한다.
3. 구조 하중 산정: 스팬 방향 양력 분포로부터 날개 구조에 작용하는 굽힘 모멘트와 전단력을 계산한다.
4. 유도 항력 추정: 설계 양력 계수에서의 유도 항력을 예측하여 비행 성능을 평가한다.

참고 문헌

• Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
• Glauert, H. (1926). The Elements of Aerofoil and Airscrew Theory. Cambridge University Press.
• Katz, J., & Plotkin, A. (2001). Low-Speed Aerodynamics (2nd ed.). Cambridge University Press.
• Prandtl, L. (1919). Tragflügeltheorie. Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, 451–477.

v 0.1.0