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

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

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

베르누이 정리는 비점성·비압축성·정상·보존력 조건에서 Euler 방정식을 유선을 따라 적분하여 얻어지는 에너지 보존의 유체역학적 표현이다. 유선 s의 방향 단위 벡터를 따라 운동량 방정식을 내적하고 적분하면

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

의 관계가 성립하며, 각 항은 단위 체적당 정압, 동압, 중력 퍼텐셜 에너지에 대응한다. 비회전 유동에서는 이 상수가 유동장 전반에서 동일하게 유지되므로, 멀리 떨어진 자유류 조건과 표면의 국소 조건을 직접 연결하는 도구가 된다. 회전 성분이 존재하는 유동에서는 유선마다 상수값이 달라지며, 이 경우 베르누이 관계는 개별 유선을 따라서만 적용된다. 이러한 적용 조건의 명확한 인식은 실제 해석에서 오용을 방지하는 필수적 전제이다.

저Mach 영역에서 공기의 밀도 변화가 미소하다는 가정 아래 비압축 형태가 유효하지만, Mach 수가 상승하면 압축성의 영향이 유의해지며 확장된 형태가 필요하다. 정상·비점성·단열 유동에서는 엔탈피와 운동 에너지의 합이 유선을 따라 보존되므로 h + \tfrac{1}{2}V^{2} = h_{0}이 성립하며, 등엔트로피 관계 p/\rho^{\gamma} = \text{const}와 결합되어

\frac{\gamma}{\gamma - 1}\frac{p}{\rho} + \tfrac{1}{2}V^{2} = \text{const}

의 압축성 베르누이 형태로 표현된다. 이러한 확장은 회전익의 팁 영역과 같이 국소 Mach 수가 일정 수준을 넘는 경우의 압력 예측에 사용되며, 비행 로봇의 일반적 저Mach 영역에서는 고전 형태가 여전히 주된 도구로 기능한다.

2. 표면 압력 분포와 압력 계수

익형 주위의 유동장은 자유류 상태로부터 표면의 속도 분포에 의해 변형되며, 이 속도 분포와 압력 분포 사이의 관계를 정량적으로 표현한 것이 압력 계수(pressure coefficient) C_{p}이다. 기준 자유류 정적 압력 p_{\infty}과 동압 q_{\infty} = \tfrac{1}{2}\rho_{\infty} V_{\infty}^{2}을 사용하여

C_{p} = \frac{p - p_{\infty}}{q_{\infty}}

로 정의되며, 비점성·비압축 조건에서는 베르누이 관계로부터 C_{p} = 1 - (V/V_{\infty})^{2}의 간결한 형태가 성립한다. 정체점에서는 국소 속도가 0이므로 C_{p,\text{stag}} = 1이며, 멀리 떨어진 자유류와 동일한 속도로 회복되는 영역에서는 C_{p} = 0이 된다. 이러한 정의는 서로 다른 크기와 속도 조건의 유동을 공통된 기준 위에서 비교할 수 있게 하며, 공기역학 해석과 실험의 공용 언어로 기능한다.

익형 표면의 C_{p} 분포는 상면과 하면에서 서로 다른 경향을 보인다. 상면에서는 속도가 증가하여 C_{p}가 음의 값을 가지고, 하면에서는 상대적으로 속도가 감소하여 C_{p}가 양의 값을 가지는 것이 일반적이다. 이 두 곡선의 차이를 현 방향으로 적분한 결과가 단면 양력 계수 C_{l}과 직접 연결되며, 실험적으로 얻어진 표면 압력 분포와 수치 해석의 결과를 직접 비교할 수 있는 정량적 지표를 제공한다. 또한 표면 압력의 극값과 분리점의 위치, 분리 버블의 흔적과 같은 세부 구조를 C_{p} 분포로부터 식별할 수 있으므로, 이는 공력 설계와 해석의 진단 도구로 폭넓게 사용된다.

3. 압축성 보정과 임계 압력 계수

저속 비압축 해석에서 얻어진 C_{p}를 아음속 압축성 영역으로 확장하기 위해 Prandtl–Glauert 보정이 사용된다. 이 보정은

C_{p} = \frac{C_{p,0}}{\sqrt{1 - \mathrm{Ma}_{\infty}^{2}}}

의 형태로 주어지며, C_{p,0}는 동일한 기하와 받음각에 대한 비압축성 해이다. 이 관계는 1차 섭동 해석의 결과로서 \mathrm{Ma}_{\infty}가 0.7 이하의 영역에서 실용적 정확도를 가진다. 보다 정밀한 보정으로는 Karman–Tsien 방법이 있으며, 이는 C_{p,0}의 크기 의존성을 추가로 반영한다. 천음속에 근접할수록 보정 오차가 증가하므로, 이 영역에서는 전체 압축성 해석 또는 천음속 CFD 해석이 요구된다. 이러한 보정의 체계적 사용은 저속 해석 결과를 중속 영역으로 확장하는 실용적 수단이 된다.

압축성 효과가 중요해지는 임계 조건은 임계 압력 계수 C_{p,\text{cr}}로 표현되며, 이는 표면 국소 Mach 수가 1에 도달하는 자유류 Mach 수에 대응하는 C_{p}의 임계값이다. 완전 등엔트로피 관계로부터 임계 압력 계수는

C_{p,\text{cr}} = \frac{2}{\gamma \mathrm{Ma}_{\infty}^{2}}\left\{\left[\frac{2 + (\gamma - 1)\mathrm{Ma}_{\infty}^{2}}{\gamma + 1}\right]^{\gamma/(\gamma - 1)} - 1\right\}

로 주어지며, 익형 표면의 최소 C_{p}가 이 값과 같아지는 자유류 Mach 수가 해당 익형의 임계 Mach 수 \mathrm{Ma}_{\text{cr}}이다. 기체 설계에서는 이 임계 Mach 수가 비행 포락선의 상한으로 부과되며, 파동 항력의 발생과 비선형 공력 변화를 피하기 위한 실용적 지표로 작동한다.

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

비행 로봇의 공력 해석에서 베르누이 정리와 C_{p}는 여러 층위에서 활용된다. 대기 데이터 계측에서는 피토 튜브의 정체 압력과 정적 압력 차이로부터 공기속도를 V = \sqrt{2(p_{0} - p_{\infty})/\rho}로 산출하며, 이 관계는 제어와 포락선 보호에 필요한 속도 정보의 원천이 된다. 익형과 날개의 설계에서는 표면 C_{p} 분포가 양력 기여 영역을 시각화하고, 분리 버블의 존재 여부를 진단하는 도구로 사용되어, 수치 해석과 풍동 시험의 결과를 교차 검증하는 지표로 기능한다. 덕트 프로펠러와 냉각 유로의 설계에서도 내부 유동의 속도·압력 관계가 동일한 관계식으로 기술되며, 기계적 설계의 공력적 근거를 제공한다.

해석 전략 관점에서 베르누이 정리는 비점성 이상 이론의 출발점이지만, 실제 설계에는 점성 경계층 효과와의 결합이 반드시 함께 수행되어야 한다. 설계 초기에는 포텐셜 유동 해석과 C_{p} 분포로부터 받음각에 대한 양력 곡선의 초기 추정을 산출하고, 이후 경계층 해석을 통해 천이 위치와 분리 거동을 평가한다. CFD 해석의 후처리에서도 C_{p} 분포는 결과 검증의 핵심 지표로 사용되며, 실험과의 비교를 통해 난류 모델과 격자 해상도의 적절성을 판단한다. 이와 같이 베르누이 정리와 압력 분포는 공기역학 해석의 고전적 도구에 머물지 않고, 현대 비행 로봇 공학의 설계·시험·운용 전반을 관통하는 공용 해석 언어로 지속적으로 기능한다.

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.
  • 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.
  • Drela, M., Flight Vehicle Aerodynamics, MIT Press, 2014.
  • Abbott, I. H., and von Doenhoff, A. E., Theory of Wing Sections, Dover, 1959.

6. 버전

v1.0