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

유동장(flow field)은 공기 유동의 공간적·시간적 분포를 기술하는 수학적 구조이며, 공기역학 해석의 대상이자 출발점이다. 유동장을 정량적으로 기술하기 위해서는 밀도, 압력, 속도, 온도, 엔트로피 등의 기본 변수와 이들의 시공간 분포가 정의되어야 하며, 이들 변수는 물리학적 의미와 측정 가능성 모두에서 엄밀히 규정된다. 본 절에서는 유동장의 수학적 표현과 기본 변수의 정의, 라그랑주 기술과 오일러 기술의 구분, 스칼라장과 벡터장의 개념, 물질 도함수, 정상 유동과 비정상 유동의 분류, 유동 표현을 위한 개념적 도구를 체계적으로 정리한다.

1. 유동장의 수학적 표현

유동장은 공간과 시간에 대한 변수의 함수로 표현된다. 공간의 한 점 $\mathbf{r} = (x, y, z)$ 와 시간 $t$ 에서 유동의 상태는 속도 벡터 $\mathbf{V}(\mathbf{r}, t)$ , 압력 $p(\mathbf{r}, t)$ , 밀도 $\rho(\mathbf{r}, t)$ , 온도 $T(\mathbf{r}, t)$ 등으로 기술된다. 이 변수들은 연속적인 함수이며, 적절한 미분 가능성을 만족한다는 연속체 가정 하에서 Navier-Stokes 방정식의 해로서 결정된다. 유동장의 해석은 이 연속 함수들의 시공간 분포를 찾는 문제이며, 해석적 방법과 수치적 방법이 모두 사용된다. 비압축성 유동에서는 밀도가 상수이므로, 유동장의 기본 변수는 속도와 압력으로 축소된다.

2. 라그랑주 기술과 오일러 기술

유동의 운동을 기술하는 두 가지 근본적 관점이 있다. 라그랑주 기술(Lagrangian description)은 개별 유체 입자(fluid particle)의 시간에 따른 궤적과 상태를 추적하는 방법이며, 고전역학의 입자 역학과 유사하다. 각 유체 입자는 초기 위치 $\mathbf{r}_0$ 으로 구분되며, 시간 $t$ 에서의 위치는 $\mathbf{r}(\mathbf{r}_0, t)$ 로 표현된다. 오일러 기술(Eulerian description)은 공간의 고정된 점에서 시간에 따라 통과하는 유체의 상태를 관찰하는 방법이며, 공간과 시간을 독립 변수로 사용한다. 공기역학에서는 오일러 기술이 일반적으로 사용되며, Navier-Stokes 방정식은 오일러 기술의 수학적 표현이다. 라그랑주 기술은 입자 추적, 오염물질 분산, 자유 표면 유동 등 특정 응용에서 유용하다.

3. 스칼라장과 벡터장

유동장의 변수는 수학적 성격에 따라 스칼라장과 벡터장으로 구분된다. 스칼라장(scalar field)은 공간과 시간의 각 점에 단일 숫자 값을 할당하는 장이며, 압력, 밀도, 온도, 엔트로피 등이 이에 해당한다. 벡터장(vector field)은 각 점에 벡터를 할당하는 장이며, 속도 벡터가 가장 중요한 예이다. 소용돌이 벡터(vorticity), 운동량 벡터, 힘 벡터 등도 벡터장이다. 벡터장의 공간적 구조는 발산(divergence), 회전(curl), 구배(gradient) 등의 미분 연산자를 통해 특성화된다. 텐서장(tensor field)은 응력 텐서와 변형률 텐서 등에서 나타나며, 각 점에 텐서 값을 할당한다.

4. 속도 벡터와 그 성분

속도 벡터 $\mathbf{V} = (u, v, w)$ 는 유체 입자의 운동을 기술하는 기본 벡터장이다. 직교 좌표계에서 속도 성분은 각각 $x$ , $y$ , $z$ 방향의 성분이며, 원통 좌표계와 구면 좌표계에서는 해당 좌표의 방향 성분으로 표현된다. 속도의 크기 $V = |\mathbf{V}| = \sqrt{u^2 + v^2 + w^2}$ 는 유동의 속력이며, 방향은 $\mathbf{V}/V$ 로 주어진다. 속도 벡터의 시공간 분포는 유동장 해석의 중심이며, 유선, 유적선, 유맥선 등의 기하학적 도구로 가시화된다.

5. 압력과 정적/정체/동적 압력

압력은 유체가 단위 면적당 물체 표면 또는 다른 유체 요소에 가하는 수직 방향 힘이며, 스칼라 물리량이다. 정적 압력(static pressure)은 유체와 함께 움직이는 관찰자가 측정하는 압력이며, 유체의 열역학적 상태를 특성화한다. 정체 압력(stagnation pressure)은 유체가 등엔트로피 과정을 거쳐 감속되어 정지할 때 도달하는 압력이며, 정적 압력과 운동 에너지 성분의 합으로 표현된다. 동적 압력(dynamic pressure)은 유체의 운동 에너지를 압력 단위로 표현한 값이며 다음과 같이 정의된다.

$q = \dfrac{1}{2} \rho V^2$

동적 압력은 양력과 항력의 정량화에 기준으로 사용되며, 양력 계수와 항력 계수의 정의에 필수적이다. 비압축성 유동에서는 Bernoulli 정리에 의해 정적 압력과 동적 압력의 합이 정체 압력이 된다.

21.4.6 물질 도함수와 가속도

오일러 기술에서 유체 입자의 가속도를 계산할 때에는 시간 미분이 특별한 형태를 가진다. 물질 도함수(material derivative) 또는 전체 도함수는 한 유체 입자를 따라가면서 관측한 변화율이며, 다음과 같이 정의된다.

$\dfrac{D}{Dt} = \dfrac{\partial}{\partial t} + \mathbf{V} \cdot \nabla$

여기서 첫 항은 특정 점에서의 시간에 따른 변화(국소 변화), 두 번째 항은 유체 입자가 공간을 이동하면서 겪는 변화(대류 변화)이다. 물질 도함수는 유체 입자의 가속도를 계산하는 데 사용되며, Navier-Stokes 방정식의 좌변을 구성한다. 물질 도함수의 개념은 오일러와 라그랑주 기술을 연결하는 핵심적 수학적 도구이다.

6. 정상 유동과 비정상 유동

유동은 시간 의존성에 따라 정상 유동(steady flow)과 비정상 유동(unsteady flow)으로 구분된다. 정상 유동은 공간의 고정된 점에서 유동의 모든 변수가 시간에 따라 변하지 않는 유동이며, 수학적으로 $\partial/\partial t = 0$ 을 만족한다. 비정상 유동은 시간에 따라 변하는 유동이며, 과도 상태, 주기적 유동, 난류 등이 이에 해당한다. 정상 유동의 가정은 해석을 크게 단순화하며, 많은 실용적 공기역학 문제가 정상 유동 가정 하에서 해석된다. 그러나 비정상 효과는 호버링에서 전진 비행으로의 천이, 돌풍 응답, 맥동 유동, 블레이드 통과 등에서 중요하며, 비정상 공기역학의 연구 주제가 된다.

7. 균일 유동과 비균일 유동

유동의 공간적 분포에 따라 균일 유동(uniform flow)과 비균일 유동(non-uniform flow)으로도 구분된다. 균일 유동은 공간의 모든 점에서 유동의 변수가 동일한 유동이며, 해석의 이상화된 경우이다. 자유 흐름(free stream)은 물체로부터 충분히 떨어진 영역에서 균일 유동으로 근사되며, 풍동 시험과 공력 해석의 기준 조건으로 사용된다. 비균일 유동은 공간에 따라 변수가 변화하는 일반적 유동이며, 실제 유동은 거의 모두 비균일 유동이다. 균일 자유 흐름에서 비균일 국소 유동으로의 변화가 물체 주변의 공기역학적 상호작용을 결정한다.

8. 1차원·2차원·3차원 유동

유동의 공간 차원에 따라 1차원, 2차원, 3차원 유동으로 분류된다. 1차원 유동은 한 좌표만의 함수로 기술되는 유동이며, 파이프 유동이나 노즐 유동의 단순화된 해석에 사용된다. 2차원 유동은 두 좌표의 함수이며, 평면 유동 또는 축대칭 유동의 해석에 사용된다. 무한 스팬의 날개 주변 유동은 2차원 유동으로 이상화되며, 익형의 공력 특성 해석의 기반이 된다. 3차원 유동은 세 좌표 모두의 함수이며, 실제 유한 날개와 동체 주변의 유동이 이에 해당한다. 차원의 감소는 해석을 단순화하지만, 실제 유동의 복잡한 현상을 정확히 포착하려면 3차원 해석이 필요하다.

9. 점성 유동과 비점성 유동

유동은 점성의 고려 여부에 따라 점성 유동(viscous flow)과 비점성 유동(inviscid flow)으로 구분된다. 비점성 유동은 점성을 무시한 이상 유체의 유동이며, Euler 방정식에 의해 기술된다. 경계층 외부의 자유 흐름 영역에서는 점성 효과가 무시할 만하여 비점성 근사가 유효하다. 점성 유동은 점성 효과를 포함한 실제 유체의 유동이며, Navier-Stokes 방정식으로 기술된다. 점성은 경계층, 박리, 후류, 마찰 항력, 층류-난류 천이 등의 핵심 현상을 지배한다. Prandtl의 경계층 이론은 비점성 유동과 점성 유동을 접합하여 실제 유동을 해석하는 기법을 제공하였다.

10. 회전 유동과 비회전 유동

유동은 소용돌이의 존재에 따라 회전 유동(rotational flow)과 비회전 유동(irrotational flow)으로 구분된다. 소용돌이 벡터는 속도의 회전(curl)으로 정의되며,

$\boldsymbol{\omega} = \nabla \times \mathbf{V}$

비회전 유동은 $\boldsymbol{\omega} = 0$ 을 만족하는 유동이며, 경계층 외부의 비점성 유동이 이에 해당한다. 비회전 유동은 포텐셜 유동(potential flow)으로도 불리며, 속도가 스칼라 포텐셜의 구배로 표현된다. 포텐셜 유동 이론은 해석적 해를 구할 수 있는 단순한 유동 모델을 제공하며, 얇은 익형 이론과 리프팅 라인 이론의 기반이 된다. 회전 유동은 점성 효과, 경계층, 후류, 소용돌이가 지배적인 영역의 유동이다.

21.4.12 압축성 유동과 비압축성 유동

유동은 밀도 변화의 고려 여부에 따라 압축성 유동(compressible flow)과 비압축성 유동(incompressible flow)으로 구분된다. 비압축성 유동은 밀도가 상수로 간주되는 유동이며, 저속 유동( $M < 0.3$ )에서 좋은 근사이다. 압축성 유동은 밀도 변화가 중요한 유동이며, 고속 유동과 충격파, 팽창파가 나타나는 유동을 포함한다. 대부분의 드론과 소형 무인기는 비압축성 영역에서 동작하며, 해석이 크게 단순화된다. 고정익 UAV 중 일부는 아음속 압축성 효과를 고려해야 할 수 있으며, 프로펠러 블레이드 끝단에서는 국소적으로 압축성 효과가 나타날 수 있다.

21.4.13 유동 해석의 통합적 관점

유동장의 기본 변수와 분류 체계는 유동 해석의 출발점이다. 실제 해석에서는 유동의 성질에 따라 적절한 가정과 방법이 선택되며, 이는 해석의 정확도와 효율을 결정한다. 예를 들어 저속 비압축성 정상 비점성 비회전 유동은 포텐셜 유동 이론으로 해석되며, 해석적 해가 가능한 가장 단순한 경우이다. 점성 효과가 중요한 경우 경계층 이론이 적용되며, 난류가 지배적이면 난류 모델이 도입된다. 고속 유동에서는 압축성이 고려되고, 비정상 현상에서는 시간에 따른 해석이 수행된다. 이러한 분류와 가정의 계층 구조를 이해하는 것은 공기역학 해석의 핵심적 출발점이다.

21.4.14 요약과 후속 연결

유동장의 기본 변수는 속도, 압력, 밀도, 온도 등이며, 이들은 시공간의 함수로 표현된다. 라그랑주 기술과 오일러 기술은 유동 운동의 두 가지 기본 관점을 제공하며, 오일러 기술이 공기역학 해석의 주류이다. 스칼라장과 벡터장, 물질 도함수, 정상/비정상 유동, 점성/비점성, 회전/비회전, 압축성/비압축성 등의 분류는 유동 해석의 체계를 구성한다. 이러한 기본 개념과 분류는 이후 다룰 지배 방정식과 해석 기법의 전제이며, 로봇 비행체의 공기역학 해석에서 반복적으로 사용된다. 다음 절에서는 비압축성 유동과 압축성 유동의 구분을 더 상세히 다루어, 저속 드론과 고속 비행체의 해석 접근법의 차이를 체계적으로 이해하는 관점을 제공한다.

출처

Anderson, J. D., Fundamentals of Aerodynamics, 6th ed., McGraw-Hill, 2016.
Batchelor, G. K., An Introduction to Fluid Dynamics, Cambridge University Press, 2000.
Panton, R. L., Incompressible Flow, 4th ed., Wiley, 2013.
Currie, I. G., Fundamental Mechanics of Fluids, 4th ed., CRC Press, 2012.
Aris, R., Vectors, Tensors, and the Basic Equations of Fluid Mechanics, Dover, 1989.
Fox, R. W., McDonald, A. T., and Pritchard, P. J., Fluid Mechanics, 9th ed., Wiley, 2016.

버전

v1.0