22.3 베르누이 정리와 압력 분포 (Bernoulli’s Theorem and Pressure Distribution)

1. 베르누이 정리의 유도와 적용 조건

베르누이 정리(Bernoulli’s theorem)는 비압축성 비점성 정상 유동에서 유선(streamline)을 따라 전압(total pressure)이 보존됨을 진술하는 에너지 보존의 법칙이다:

$p + \frac{1}{2}\rho V^2 + \rho gz = \text{const.} \quad (\text{유선을 따라})$

여기서 $p$ 는 정적 압력(static pressure), $\frac{1}{2}\rho V^2$ 는 동압(dynamic pressure), $\rho gz$ 는 위치 에너지 항이다. 공기역학에서 중력 항( $\rho gz$ )은 밀도가 낮은 공기의 경우 무시할 수 있을 정도로 작으므로, 다음과 같이 단순화된다:

$p + \frac{1}{2}\rho V^2 = p_0 = \text{const.}$

여기서 $p_0$ 는 전압(total pressure) 또는 정체 압력(stagnation pressure)이다.

베르누이 정리의 적용 조건:

비압축성 유동( $M < 0.3$ )
비점성 유동(경계층 외부)
정상 유동(시간 불변)
동일 유선 위의 두 점 사이에 적용

이 조건은 경계층 외부의 자유 유동(free-stream)에서 잘 만족되며, 익형 주위의 비점성 유동 해석에 널리 적용된다. 경계층 내부에서는 점성 효과로 인해 전압이 보존되지 않으므로 베르누이 정리를 직접 적용할 수 없다.

2. 속도와 압력의 역관계

베르누이 정리의 핵심적 함의는 유선을 따라 속도가 증가하면 압력이 감소하고, 속도가 감소하면 압력이 증가한다는 것이다. 이 속도-압력 역관계(velocity-pressure inverse relationship)는 양력 발생의 기본 원리를 이해하는 데 핵심적이다.

자유류 조건( $p_\infty$ , $V_\infty$ )을 기준으로 임의 점에서의 압력은:

$p = p_\infty + \frac{1}{2}\rho(V_\infty^2 - V^2)$

$V > V_\infty$ 이면 $p < p_\infty$ (흡입 압력)
$V < V_\infty$ 이면 $p > p_\infty$ (양압)
$V = 0$ (정체점)이면 $p = p_0 = p_\infty + \frac{1}{2}\rho V_\infty^2$ (최대 압력)

22.3.3 익형 표면의 압력 분포

양의 받음각을 갖는 익형 주위의 유동에서 표면 압력 분포는 다음의 특징적 양상을 나타낸다:

정체점: 앞전(leading edge) 부근에서 유동 속도가 영이 되는 정체점이 존재한다. 양의 받음각에서 정체점은 앞전의 하면 쪽으로 이동한다. 정체점에서의 압력은 전압 $p_0$ 이며, 압력 계수는 $C_p = 1$ 이다.

상면(upper surface): 정체점으로부터 상면을 따라 유동이 급격히 가속되어 속도가 자유류보다 빠르게 된다. 이에 따라 상면의 압력은 자유류 압력보다 낮아진다( $C_p < 0$ ). 최소 압력점(흡입 피크, suction peak)에서 속도가 최대이고 압력이 최소이다.

하면(lower surface): 유동이 상면보다 상대적으로 느리게 흐르며, 압력이 자유류 압력에 가깝거나 약간 높다( $C_p \geq 0$ ). 큰 받음각에서는 하면의 앞전 부근에서 국소적 가속이 나타나기도 한다.

뒷전(trailing edge): 상면과 하면의 압력이 수렴하며, 쿠타 조건에 의해 두 값이 거의 동일해진다.

22.3.4 압력 계수에 의한 무차원 표현

압력 분포를 무차원 형태로 표현하기 위해 압력 계수(pressure coefficient, $C_p$ )가 사용된다:

$C_p = \frac{p - p_\infty}{q_\infty} = 1 - \left(\frac{V}{V_\infty}\right)^2$

$C_p$ 분포도는 시위 방향 위치 $x/c$ 를 가로축으로, $C_p$ 를 세로축(반전)으로 나타내는 것이 관례이다. 세로축 반전에 의해 흡입 영역(음의 $C_p$ )이 도표의 상부에 위치하여, 상면의 흡입 압력 분포를 직관적으로 파악할 수 있다.

$C_p$ 분포의 면적(상면과 하면 곡선이 둘러싸는 영역)은 법선력 계수에 비례하며, 이는 양력 계수와 직접 관련된다:

$C_n = \frac{1}{c}\int_0^c (C_{p,l} - C_{p,u})\,dx$

여기서 $C_{p,l}$ 과 $C_{p,u}$ 는 각각 하면과 상면의 압력 계수이다.

22.3.5 압력 분포와 양력의 관계

양력은 표면 압력 분포를 양력 방향으로 적분하여 산출된다. 상면의 흡입 압력(음의 $C_p$ )은 양력의 약 60~70%를 기여하고, 하면의 양압(양의 $C_p$ )은 나머지 30~40%를 기여한다. 따라서 양력 생성에서 상면의 흡입이 하면의 양압보다 더 지배적인 역할을 수행한다.

이 사실은 양력 발생의 설명에서 “상면의 저압“이 “하면의 고압“보다 더 큰 기여를 한다는 점에서 중요하며, 상면 유동의 가속과 이에 따른 압력 감소가 양력의 핵심 메커니즘임을 확인시켜 준다.

22.3.6 순압력 구배와 역압력 구배

베르누이 정리에 의해 유동의 가감속은 압력 구배(pressure gradient)를 형성한다:

순압력 구배(favorable pressure gradient, $dp/dx < 0$ ): 유동이 가속되는 영역에서 나타나며, 경계층의 안정성을 향상시키고 박리를 억제한다. 앞전으로부터 최소 압력점까지의 구간이 이에 해당한다.

역압력 구배(adverse pressure gradient, $dp/dx > 0$ ): 유동이 감속되는 영역에서 나타나며, 경계층의 성장을 촉진하고 박리를 유발할 수 있다. 최소 압력점으로부터 뒷전까지의 구간이 이에 해당한다.

역압력 구배의 강도는 익형의 형상과 받음각에 의해 결정되며, 과도한 역압력 구배는 경계층 박리(boundary layer separation)를 통해 양력 손실과 항력 증가를 초래한다.

22.3.7 받음각 변화에 따른 압력 분포의 변화

받음각의 증가에 따라 익형의 압력 분포는 다음과 같이 변화한다:

정체점이 하면 쪽으로 이동한다.
상면의 흡입 피크 강도가 증가하고 위치가 앞전 쪽으로 이동한다.
상면의 역압력 구배가 강화된다.
상면과 하면의 $C_p$ 차이가 증가하여 양력이 증가한다.
실속 받음각에 도달하면 상면의 과도한 역압력 구배에 의해 대규모 박리가 발생하고, 흡입 피크가 붕괴하여 양력이 감소한다.

22.3.8 비정상 베르누이 방정식

비정상 유동(unsteady flow)에서는 정상 베르누이 방정식에 시간 변화 항이 추가된다:

$\frac{\partial \phi}{\partial t} + \frac{p}{\rho} + \frac{V^2}{2} = \text{const.}$

여기서 $\phi$ 는 속도 포텐셜이다. 비정상 항 $\partial\phi/\partial t$ 는 비정상 유동에서의 추가적인 압력 성분(부가 질량 효과, added mass effect)을 나타내며, 돌풍 응답이나 동적 실속의 해석에서 중요하다.

3. 압축성 유동에서의 압력 관계

마하 수가 0.3을 초과하는 압축성 유동에서는 비압축성 베르누이 방정식 대신 등엔트로피 유동 관계식이 적용된다:

$\frac{p_0}{p} = \left(1 + \frac{\gamma-1}{2}M^2\right)^{\frac{\gamma}{\gamma-1}}$

아음속 압축성 유동에서 압력 분포에 대한 프란틀-글라우어트 보정(Prandtl-Glauert correction)은 비압축성 해의 압력 계수를 $1/\sqrt{1-M_\infty^2}$ 로 나누어 압축성 효과를 근사적으로 반영한다.

22.3.10 로봇 공학에서의 압력 분포 활용

비행 로봇의 설계에서 익형 표면의 압력 분포 해석은 다음의 목적에 활용된다:

양력 예측: 압력 분포의 적분으로부터 양력 계수를 산출한다.
실속 예측: 역압력 구배의 강도 분석을 통해 박리와 실속의 발생 가능성을 평가한다.
구조 설계: 표면 압력 분포는 날개와 블레이드의 구조적 하중 분포를 결정한다.
익형 설계: 목표 압력 분포를 지정하고 이를 실현하는 익형 형상을 역설계(inverse design)한다.

XFOIL과 같은 점성-비점성 연성 해석 코드는 비점성 유동 해석(패널 방법)에 의한 압력 분포 계산과 경계층 해석을 결합하여, 점성 효과를 포함한 보다 정확한 압력 분포를 제공한다 (Drela, 1989).

참고 문헌

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.

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