21.4 유동장의 기본 변수와 정의

1. 유동장의 수학적 구성과 기본 변수

유동장(flow field)은 공기 유동의 공간적·시간적 분포를 기술하는 수학적 구조이며, 공기역학 해석의 대상이자 출발점이다. 유동장은 관심 영역 $\Omega \subset \mathbb{R}^{3}$ 의 모든 점 $\mathbf{x} \in \Omega$ 와 시간 $t \in \mathbb{R}$ 에 대하여, 유체의 상태를 기술하는 기본 변수를 함수적으로 부여하는 사상으로 정의된다. 즉 밀도장 $\rho(\mathbf{x}, t)$ , 압력장 $p(\mathbf{x}, t)$ , 속도장 $\mathbf{u}(\mathbf{x}, t) = (u, v, w)^{\top}$ , 온도장 $T(\mathbf{x}, t)$ , 엔트로피장 $s(\mathbf{x}, t)$ 와 같은 변수들이 유동장의 구성 요소를 이룬다. 이 변수들은 연속체 가정 아래에서 일반적으로 미분 가능한 함수로 간주되며, 충격파와 같은 불연속이 존재하는 영역에서는 Rankine–Hugoniot 관계에 의한 도약 조건이 별도로 부여된다.

유동장의 기본 변수는 열역학 관계에 의하여 서로 결합된다. 건조 공기를 완전 기체로 근사하면 상태 방정식 $p = \rho R T$ 가 성립하고, 엔탈피와 내부 에너지는 $h = c_{p} T$ , $e = c_{v} T$ 로 주어지며, 비열비는 $\gamma = c_{p}/c_{v}$ 로 정의된다. 이러한 열역학 관계는 유동장의 독립 변수를 두 개로 축소하며, 따라서 $p$ , $\rho$ , $T$ 가운데 두 변수와 속도장이 주어지면 나머지 변수가 결정된다. 실제 계측에서는 압력과 온도가 가장 직접적으로 측정되므로, 이들을 독립 변수로 선택하고 밀도는 상태 방정식을 통하여 산출하는 방식이 표준적으로 사용된다. 이러한 구조화는 유동장의 해석과 계측을 일관된 수학적 체계 안에서 결합할 수 있도록 한다.

2. 라그랑주 기술과 오일러 기술, 그리고 물질 도함수

유동장을 기술하는 관점은 관찰자의 위치와 운동 방식에 따라 라그랑주 기술(Lagrangian description)과 오일러 기술(Eulerian description)로 나뉜다. 라그랑주 기술은 초기 시각 $t_{0}$ 에 위치 $\mathbf{X}$ 에 있던 유체 입자가 시간이 지남에 따라 어떤 궤적 $\boldsymbol{\phi}(\mathbf{X}, t)$ 을 따라 이동하는지를 추적하며, 각 입자에 대하여 상태 변수의 시간 이력을 부여한다. 오일러 기술은 공간의 고정점 $\mathbf{x}$ 에서 유체 변수의 시간 변화를 관찰하며, 유동장을 공간 분포의 시간적 변화로 기술한다. 공기역학의 수학적 전개는 경계 조건 부여가 용이하고 계측 관점에 부합하는 오일러 기술을 주로 채택하며, 라그랑주 기술은 입자 추적, 후류 구조 분석, 공정 경로 해석과 같은 특정 응용에서 제한적으로 활용된다.

두 기술은 물질 도함수(material derivative)를 통하여 수학적으로 연결된다. 임의의 스칼라 장 $\phi(\mathbf{x}, t)$ 에 대한 물질 도함수는

$\frac{D\phi}{Dt} = \frac{\partial \phi}{\partial t} + \mathbf{u}\cdot \nabla \phi$

로 정의되며, 이는 유체 입자가 흐름을 따라 이동하는 동안 관측되는 $\phi$ 의 시간 변화율을 뜻한다. 여기서 우변의 첫 항은 공간의 고정점에서 관찰되는 시간 변화, 두 번째 항은 유체가 공간적 비균일성을 통과하면서 누적하는 변화를 각각 의미한다. 벡터 장 $\mathbf{u}$ 에 대한 물질 도함수는 $D\mathbf{u}/Dt = \partial \mathbf{u}/\partial t + (\mathbf{u}\cdot\nabla)\mathbf{u}$ 이며, 이는 가속도장을 오일러 변수로 표현하는 정의이다. 물질 도함수는 Navier–Stokes 방정식을 비롯한 지배 방정식의 표준적 서술 방식으로 사용되며, 라그랑주 관점의 입자적 직관을 오일러 기술 안에서 정량적으로 회복하는 수단으로 작동한다.

3. 정상·비정상 유동, 장의 기하학적 구조, 유동 분류

유동장의 시간 의존성은 정상 유동(steady flow)과 비정상 유동(unsteady flow)의 구분을 가능하게 한다. 정상 유동에서는 공간의 모든 점에서 기본 변수가 시간에 대하여 불변이므로 $\partial \phi/\partial t = 0$ 이 성립하며, 물질 도함수는 $D\phi/Dt = \mathbf{u}\cdot\nabla \phi$ 로 축약된다. 비정상 유동에서는 $\partial \phi/\partial t \neq 0$ 이 일반적이며, 회전익의 후류, 돌풍 하의 비행, 천이 경계층, 플러터와 같은 현상이 이에 해당한다. 비정상성은 Strouhal 수 $\mathrm{St} = fL/V$ 로 정량화되며, 이 값이 1 근처의 규모를 가질 때 비정상 효과가 유동장의 지배적 구조를 결정한다. 비행 로봇 해석에서는 호버링을 정상 유동으로, 돌풍 경감 제어와 기동 상태를 비정상 유동으로 각각 이상화하여 해석 전략을 구분한다.

유동장의 기하학적 구조는 유선(streamline), 유적선(pathline), 와선(streakline)과 같은 곡선에 의하여 시각화된다. 유선은 특정 시각의 속도장에 접하는 곡선으로 $d\mathbf{x}/ds = \mathbf{u}(\mathbf{x}, t_{0})$ 의 해로 정의되고, 유적선은 한 유체 입자의 실제 궤적 $\boldsymbol{\phi}(\mathbf{X}, t)$ 을 의미하며, 와선은 특정 공간점을 통과한 유체 입자들의 시간적 자취이다. 정상 유동에서는 세 곡선이 일치하지만 비정상 유동에서는 일반적으로 서로 다르며, 이러한 구분은 유동 시각화 실험과 계산 결과 해석에서 혼동을 방지하는 수단이 된다. 또한 유동장의 국소적 성질은 속도 구배 텐서 $\nabla \mathbf{u}$ 의 대칭부 변형률 텐서 $\mathbf{S} = \tfrac{1}{2}(\nabla\mathbf{u} + \nabla\mathbf{u}^{\top})$ 와 반대칭부 회전 텐서 $\boldsymbol{\Omega} = \tfrac{1}{2}(\nabla\mathbf{u} - \nabla\mathbf{u}^{\top})$ 로 분해되며, 이 분해는 유동의 변형과 회전 특성을 분리하여 해석하는 기하학적 기초를 제공한다.

유동장은 추가적인 무차원·물리적 기준에 따라 세분된다. 첫째, Reynolds 수에 따라 층류와 난류로 구분되고, 둘째, Mach 수에 따라 비압축성·아음속·천음속·초음속·극초음속으로 구분되며, 셋째, 와도 $\boldsymbol{\omega} = \nabla \times \mathbf{u}$ 의 유무에 따라 비회전 유동과 회전 유동으로 구분된다. 비회전·비점성·비압축의 조합은 포텐셜 유동(potential flow) 이론을 성립시켜 $\mathbf{u} = \nabla \phi$ , $\nabla^{2}\phi = 0$ 의 단순한 해석 구조를 제공하며, 얇은 익형 이론과 리프팅 라인 이론의 이론적 기반으로 기능한다. 이러한 분류는 단순한 명명법이 아니라, 서로 다른 근사와 해석 기법을 선택하는 공학적 판단의 기준을 제공한다.

4. 유동장 기술의 로봇공학적 활용

비행 로봇의 해석과 제어에서 유동장의 기본 변수는 상태 방정식·계측·제어 합성의 전 단계에 걸쳐 활용된다. 대기 데이터 시스템은 정압 $p_{\infty}$ 와 정체 압력 $p_{0}$ 의 차이를 이용하여 공기속도를 $V = \sqrt{2(p_{0} - p_{\infty})/\rho}$ 로 산출하고, 온도 센서와 압력 센서의 조합을 통하여 실시간 밀도 $\rho = p/(RT)$ 를 재구성한다. 이러한 계측값은 비행 제어기의 공력 게인 스케줄링, 추력 관리, 포락선 보호에 입력으로 사용되며, 유동장의 기본 변수가 단순한 해석적 추상이 아니라 운용 중의 정량적 제어 신호로 작동함을 보여 준다. 온보드 대기 데이터의 정확성은 공력 모델의 유효성과 직접적으로 결합되어 있다.

또한 유동장의 기하학적 구조는 비행 로봇의 궤적 계획과 군집 제어에서 실질적 의미를 갖는다. 선행기의 후류는 유동장에 시간 지연된 속도 섭동을 유발하며, 이는 후속 기체의 받음각과 받는 힘을 변동시켜 군집 안정성에 영향을 미친다. 도심 비행에서는 건물 사이의 가속·감속 통로와 와류 구조가 국부적으로 강한 속도 구배를 형성하여 자세 제어의 난제로 작용하고, 그 해석에는 $\nabla \mathbf{u}$ 기반의 변형률·회전 분해가 사용된다. 이와 같은 실무적 맥락은 유동장의 기본 변수와 정의가 공기역학의 순수 이론을 넘어, 비행 로봇의 설계와 임무 수행에서 지속적으로 참조되는 공학적 언어임을 명확히 드러낸다.

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., Viscous Fluid Flow, 3rd ed., McGraw-Hill, 2006.
Panton, R. L., Incompressible Flow, 4th ed., Wiley, 2013.
Katz, J., and Plotkin, A., Low-Speed Aerodynamics, 2nd ed., Cambridge University Press, 2001.

6. 버전

v1.0