21.8 베르누이 정리와 공기역학적 적용

1. 유도와 수학적 형태

베르누이 정리는 비점성·비압축성·정상·보존력 조건에서 Euler 방정식을 유선에 따라 적분하여 얻어지는 관계이며, 공기역학에서 가장 널리 사용되는 간결한 보존 법칙이다. 유선 $s$ 에 대한 방향 단위 벡터 $\hat{\mathbf{s}}$ 의 방향으로 Euler 방정식을 내적하면

$\rho V \frac{dV}{ds} = -\frac{dp}{ds} - \rho g \frac{dz}{ds}$

가 성립하며, 양변을 $s$ 에 대하여 적분하면 유선을 따라

$p + \tfrac{1}{2}\rho V^{2} + \rho g z = \text{const}$

의 관계가 얻어진다. 이 식의 각 항은 단위 체적당 정압, 동압, 중력 퍼텐셜 에너지에 대응하며, 에너지 보존의 유체역학적 재서술로 해석할 수 있다. 비회전 유동의 경우 이 상수가 유동장 전반에 걸쳐 동일하게 유지되며, 회전이 존재하는 경우에는 유선마다 상수값이 다르게 설정된다. 이러한 구분은 실제 적용에서 유동장의 회전성 여부를 먼저 판정해야 하는 이유를 제공한다.

베르누이 정리는 여러 변형 형태로 사용되며, 각각은 적용 대상에 대한 물리적 강조점을 달리한다. 압축성을 무시할 수 있는 경우에는 앞서 제시한 고전 형태가 기본이 되고, 고도 변화가 미미한 공기역학 해석에서는 $p + \tfrac{1}{2}\rho V^{2} = p_{0}$ 의 형태가 일상적으로 사용되며, 정체 압력 $p_{0}$ 은 속도가 0인 유선 종점에서의 압력으로 정의된다. 압축성이 유의한 경우에는 에너지 방정식을 결합하여

$h + \tfrac{1}{2} V^{2} = h_{0}$

의 형태가 성립하며, 등엔트로피 관계 $p/\rho^{\gamma} = \text{const}$ 를 이용하여 $(p/\rho)(\gamma/(\gamma-1)) + \tfrac{1}{2}V^{2} = \text{const}$ 와 같은 확장된 압축성 형태로 재구성된다. 이러한 변형들은 동일한 물리 원리의 서로 다른 표현이며, 해석 조건에 따라 가장 적합한 형태가 선택된다.

2. 적용 조건과 유효성의 한계

베르누이 정리의 엄밀한 유도는 비점성·비압축성·정상·보존력·유선 적분이라는 전제에 의존하므로, 이러한 조건이 성립하지 않는 영역에서는 직접 적용될 수 없다. 점성이 지배적인 벽 근방의 경계층 내부에서는 유체 입자 사이의 전단 응력이 에너지를 소산시키므로 베르누이 상수가 유선을 따라 감소하며, 이러한 손실은 Darcy–Weisbach 형태의 압력 강하나 경계층 변위 두께를 통한 보정으로 다루어진다. 유동의 분리와 재부착이 존재하거나, 강한 와류가 포함되거나, 유동이 주기적으로 변동하는 영역에서도 베르누이 정리는 근사적 도구로만 사용되며, 이때의 오차는 해석의 신뢰성을 심각하게 훼손할 수 있다. 따라서 실제 해석에서는 베르누이 관계의 적용이 가능한 영역을 먼저 식별한 뒤, 경계층과 분리 영역을 별도의 모델로 보완하는 이원 해석이 표준적이다.

압축성의 영향이 본격적으로 작용하는 고Mach 영역에서도 고전 베르누이 정리는 유효성을 상실한다. 정체점에서의 등엔트로피 관계

$\frac{p_{0}}{p} = \left(1 + \frac{\gamma - 1}{2}\mathrm{Ma}^{2}\right)^{\gamma/(\gamma - 1)}$

은 비행체 주위의 압축성 정체 압력 해석에 사용되며, 이 관계는 저Mach 근사에서 $p_{0} \approx p + \tfrac{1}{2}\rho V^{2}$ 로 회복된다. Mach 수가 0.3을 초과하는 순간 보정 없이 사용되는 고전 베르누이 관계는 수 퍼센트 수준의 오차를 유발하며, 천음속 이상에서는 충격파와 총 압력 손실이 발생하여 해석의 성격 자체가 달라진다. 이러한 한계를 명시적으로 인식하고 적용 범위를 엄밀히 제한하는 것이 베르누이 정리를 실무에서 안전하게 사용하기 위한 전제 조건이다.

3. 공기역학 해석에서의 응용

베르누이 정리는 공기역학의 다양한 실무 해석에서 핵심적 도구로 기능한다. 대기 데이터 계측에서는 피토 튜브의 정체 압력과 정적 압력 차이로부터 동압 $q = \tfrac{1}{2}\rho V^{2}$ 을 산출하여 공기속도를 $V = \sqrt{2(p_{0} - p_{\infty})/\rho}$ 로 계산한다. 풍동 수축부와 자유 단면의 속도 관계를 결정하는 데에도 연속 방정식과 결합된 베르누이 관계가 사용되며, 노즐과 디퓨저의 1차원 해석, 벤추리 유량계의 보정 관계 역시 동일한 원리에 기반한다. 익형 주위의 압력 분포 해석에서도 비점성 포텐셜 해가 주어지면 속도장의 크기로부터 표면 압력 계수 $C_{p} = 1 - (V/V_{\infty})^{2}$ 을 얻을 수 있으며, 이 결과는 경계층 해석의 외곽 조건을 제공하는 이중 구조 해석의 근간이 된다.

특히 양력 생성의 정량적 설명에서 베르누이 정리는 표면 위·아래의 속도 차이에 기인한 압력 차이를 직접 관계 짓는 수단을 제공한다. 표면 위쪽에서 속도가 증가하면 정압이 낮아지고, 아래쪽에서 속도가 상대적으로 낮으면 정압이 높아지며, 이 압력 차이의 적분이 양력으로 구현된다. 다만 이 서술은 현대적 관점에서 Kutta–Joukowski 정리와 순환 이론으로 보완되어야 완전해지며, 베르누이 관계는 양력의 크기에 대한 직접적 적분 도구가 아닌 압력 분포 해석의 매개 관계로 이해되는 것이 엄밀하다. 이러한 해석적 위상은 교육적 맥락에서 단순화된 설명이 때로 인과적 오해를 낳을 수 있다는 사실을 상기시키며, 정교한 해석을 위해서는 Euler 해석과 순환 이론과 결합된 서술이 요구된다.

4. 로봇공학적 활용과 해석 전략

비행 로봇의 해석과 계측에서 베르누이 정리는 일상적으로 사용된다. 대기 데이터 시스템은 피토–정압 튜브 구조를 통하여 동압을 계측하고, 온도 센서를 결합하여 밀도 $\rho = p/(RT)$ 를 산출한 뒤 실제 공기속도 $V = \sqrt{2 q / \rho}$ 를 계산한다. 이 산출값은 자세 제어, 포락선 보호, 수평 바람장 추정에 입력되며, 돌풍 경감과 자율 이착륙 알고리즘의 기초 정보로 기능한다. 또한 환기 덕트, 냉각 유로, 프로펠러 상류·하류의 압력 분석에 베르누이 관계가 활용되어 기계적 설계와 열 관리의 기본 지침을 제공한다. 이러한 활용은 단순한 교과서 공식이 아니라, 실시간 비행 제어와 장기 성능 유지에 관여하는 공학적 장치로서의 베르누이 정리의 역할을 보여 준다.

해석 전략의 관점에서는 베르누이 정리를 전반적 해석의 출발점으로 활용하되, 점성 경계층과 분리 영역에 대한 보완 모델을 반드시 결합하는 이원 구조가 표준이다. 설계 초기 단계에서는 포텐셜 유동 해석과 베르누이 관계로부터 익형 표면 압력 분포와 양력 계수의 초기 추정을 얻고, 이후 경계층 해석을 통하여 표면 마찰과 분리 위치를 평가한다. CFD가 활용되는 상세 설계 단계에서도 베르누이 관계는 후처리의 물리적 검증 수단으로 쓰이며, 계산된 속도장과 압력장 사이의 일관성을 확인하는 간접적 진단 도구로 기능한다. 이러한 다층적 활용은 베르누이 정리가 단순한 간이식을 넘어, 공기역학 해석 전반의 내적 일관성을 유지하는 지표로 작동함을 잘 보여 준다.

5. 출처

Anderson, J. D., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2016.
Kundu, P. K., Cohen, I. M., and Dowling, D. R., Fluid Mechanics, 6th ed., Academic Press, 2015.
Batchelor, G. K., An Introduction to Fluid Dynamics, Cambridge University Press, 2000.
White, F. M., Fluid Mechanics, 8th ed., McGraw-Hill, 2016.
Katz, J., and Plotkin, A., Low-Speed Aerodynamics, 2nd ed., Cambridge University Press, 2001.
Liepmann, H. W., and Roshko, A., Elements of Gasdynamics, Dover, 2001.

6. 버전

v1.0