21.10 유선, 유적선, 유맥선의 해석 (Streamlines, Pathlines, and Streaklines)

1. 유동 가시화 곡선의 개념과 분류

유동장(flow field)의 구조를 이해하고 해석하기 위해서는 유체 입자의 운동을 시각적으로 표현하는 방법이 필요하다. 유체역학에서는 유동의 기하학적 특성을 기술하기 위해 세 가지 기본적인 곡선, 즉 유선(streamline), 유적선(pathline), 유맥선(streakline)을 정의한다. 이 세 곡선은 유동장의 서로 다른 측면을 기술하며, 정상 유동(steady flow)에서는 세 곡선이 일치하지만, 비정상 유동(unsteady flow)에서는 각각 상이한 형태를 나타낸다.

이들 곡선은 실험적 유동 가시화(flow visualization)와 전산 유체역학(CFD) 해석에서 유동장의 구조를 파악하고, 박리 영역(separation region), 와류(vortex), 정체점(stagnation point) 등 주요 유동 현상을 식별하는 데 핵심적인 도구로 활용된다.

2. 유선의 정의와 수학적 표현

유선(streamline)은 특정 순간에 유동장 내의 각 점에서 속도 벡터(velocity vector)에 접하는 곡선이다. 즉, 유선 위의 모든 점에서 유선의 접선 방향이 해당 점에서의 속도 벡터 방향과 일치한다. 유선은 특정 시각 $t$ 에서의 속도장의 순간적 스냅샷(snapshot)을 나타내는 것으로서, 시간에 대한 정보가 아닌 공간적 속도 분포에 대한 정보를 제공한다.

2차원 유동에서 속도장이 $\mathbf{V} = u(x, y, t)\,\mathbf{i} + v(x, y, t)\,\mathbf{j}$ 로 주어질 때, 유선의 미분 방정식은 다음과 같다:

$\frac{dx}{u} = \frac{dy}{v}$

이를 벡터 형태로 표현하면, 유선 위의 미소 변위 벡터 $d\mathbf{r}$ 과 속도 벡터 $\mathbf{V}$ 의 외적이 영이 되는 조건이다:

$d\mathbf{r} \times \mathbf{V} = 0$

3차원 유동에서는 이 조건이 다음의 연립 미분 방정식으로 확장된다:

$\frac{dx}{u} = \frac{dy}{v} = \frac{dz}{w}$

여기서 $u$ , $v$ , $w$ 는 각각 $x$ , $y$ , $z$ 방향의 속도 성분이다.

유선의 중요한 성질은 정상 유동에서 유체 입자가 유선을 횡단(cross)할 수 없다는 점이다. 이는 유선이 속도 벡터에 항상 접하므로, 유선에 수직한 방향의 속도 성분이 존재하지 않기 때문이다. 다만, 비정상 유동에서는 유선 자체가 시간에 따라 변하므로 유체 입자가 특정 시각의 유선을 횡단하는 것처럼 보일 수 있다.

21.10.3 유선 함수

2차원 비압축성 유동에서 유선 함수(stream function, $\psi$ )를 도입하면 유선을 보다 체계적으로 기술할 수 있다. 유선 함수는 다음의 관계를 만족하도록 정의된다:

$u = \frac{\partial \psi}{\partial y}, \quad v = -\frac{\partial \psi}{\partial x}$

이 정의에 의해 연속 방정식(continuity equation) $\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$ 이 자동으로 만족된다. 유선은 $\psi = \text{const.}$ 인 등고선(contour line)으로 표현되며, 두 유선 사이의 유선 함수 차이 $\Delta\psi$ 는 두 유선 사이를 단위 시간에 통과하는 단위 깊이당 체적 유량(volume flow rate)과 같다:

$\Delta Q = \psi_2 - \psi_1$

유선 함수의 값이 증가하는 방향은 속도 벡터의 좌측에 해당하며, 유선 간의 간격이 좁을수록 유동 속도가 빠르다. 이러한 특성은 유동장의 정성적 해석에 유용한 정보를 제공한다.

비회전 유동(irrotational flow)에서는 유선 함수가 라플라스 방정식(Laplace equation) $\nabla^2 \psi = 0$ 을 만족하며, 이는 속도 포텐셜(velocity potential) $\phi$ 의 라플라스 방정식 $\nabla^2 \phi = 0$ 과 함께 포텐셜 유동 이론(potential flow theory)의 기초를 형성한다 (Batchelor, 1967).

21.10.4 유적선의 정의와 수학적 표현

유적선(pathline)은 하나의 유체 입자가 시간에 따라 이동하는 실제 궤적(trajectory)이다. 유적선은 라그랑주 관점(Lagrangian perspective)에서 유동을 기술하는 것으로, 특정 유체 입자를 추적하여 그 위치를 시간의 함수로 기록한 것이다.

유적선의 수학적 표현은 유체 입자의 운동 방정식으로부터 직접 유도된다. 유체 입자의 위치 벡터를 $\mathbf{r}(t)$ 라 하면:

$\frac{d\mathbf{r}}{dt} = \mathbf{V}(\mathbf{r}, t)$

성분별로 표현하면:

$\frac{dx}{dt} = u(x, y, z, t), \quad \frac{dy}{dt} = v(x, y, z, t), \quad \frac{dz}{dt} = w(x, y, z, t)$

초기 조건 $\mathbf{r}(t_0) = \mathbf{r}_0$ 이 주어지면, 위의 상미분 방정식 체계를 적분하여 유적선을 결정할 수 있다. 정상 유동에서 속도장이 시간에 무관하므로( $\mathbf{V} = \mathbf{V}(\mathbf{r})$ ), 유적선의 미분 방정식이 유선의 미분 방정식과 동일해지며, 따라서 유적선과 유선이 일치한다.

실험적으로 유적선은 유동 내에 투입된 개별 추적 입자(tracer particle)의 궤적을 장시간 노출(long exposure) 촬영하여 획득할 수 있다. 이 방법은 입자 추적 유속계(Particle Tracking Velocimetry, PTV)의 기본 원리를 구성한다.

21.10.5 유맥선의 정의와 수학적 표현

유맥선(streakline)은 유동장 내의 특정 고정점을 연속적으로 통과한 모든 유체 입자들이 특정 관측 시각 $t$ 에 형성하는 곡선이다. 즉, 시각 $t$ 이전의 모든 시점에서 해당 고정점을 통과한 입자들의 현재 위치를 연결한 곡선이 유맥선이다.

유맥선의 수학적 구성은 다음과 같다. 고정점 $\mathbf{r}_0$ 를 시각 $\tau$ ( $\tau \leq t$ )에 통과한 유체 입자의 관측 시각 $t$ 에서의 위치를 $\mathbf{r}(t; \tau, \mathbf{r}_0)$ 라 하면, 유맥선은 매개변수 $\tau$ 를 $-\infty < \tau \leq t$ 범위에서 변화시킬 때 $\mathbf{r}(t; \tau, \mathbf{r}_0)$ 가 그리는 곡선이다. 각 $\tau$ 에 대해 $\mathbf{r}(t; \tau, \mathbf{r}_0)$ 는 유적선의 방정식을 초기 조건 $\mathbf{r}(\tau) = \mathbf{r}_0$ 으로 시각 $t$ 까지 적분하여 구한다.

실험적으로 유맥선은 유동 내의 특정 지점에서 염료(dye)나 연기(smoke)를 연속적으로 주입할 때 형성되는 가시적 흔적으로 직접 관찰된다. 풍동 실험에서 연기선(smoke wire) 기법이나 수조(water channel) 실험에서 염료 주입법은 유맥선을 가시화하는 대표적 기법이다.

21.10.6 세 곡선의 상호 관계

정상 유동에서는 속도장이 시간에 무관하므로( $\partial \mathbf{V}/\partial t = 0$ ), 유선, 유적선, 유맥선이 모두 동일한 곡선을 형성한다. 이는 다음과 같이 설명된다:

정상 유동에서 유선은 시간에 따라 변하지 않으므로, 유체 입자는 항상 동일한 유선을 따라 이동한다. 따라서 유적선이 유선과 일치한다.
특정 고정점을 통과한 모든 유체 입자는 동일한 유선을 따라 이동하므로, 이들의 현재 위치를 연결한 유맥선 역시 유선과 일치한다.

비정상 유동에서는 속도장이 시간에 따라 변하므로, 세 곡선은 일반적으로 서로 다른 형태를 나타낸다. 비정상 유동에서 각 곡선이 담고 있는 정보의 차이를 정리하면 다음과 같다:

유선: 특정 순간의 속도장 구조를 나타내는 순간적(instantaneous) 곡선이다.
유적선: 특정 유체 입자의 시간 이력(time history)을 나타내는 누적적(cumulative) 곡선이다.
유맥선: 특정 공간점을 통과한 유체 입자 집합의 현재 분포를 나타내는 집합적(collective) 곡선이다.

비정상 유동에서 세 곡선이 일치하지 않는다는 사실은 유동 가시화 실험의 해석에서 특히 중요하다. 실험에서 관찰되는 가시적 패턴이 유맥선에 해당하는 경우가 많으므로, 이를 유선으로 잘못 해석하면 유동장의 순간적 구조에 대한 오류가 발생할 수 있다 (Hama, 1962).

21.10.7 시간선과 기타 유동 곡선

유선, 유적선, 유맥선 외에도 유동장을 기술하는 추가적인 곡선이 존재한다. 시간선(timeline)은 특정 시각에 유동장 내의 일직선상에 배치된 유체 입자들이 이후 시간에 형성하는 곡선이다. 시간선은 유체 요소(fluid element)의 변형과 속도 분포의 공간적 변화를 가시적으로 나타내며, 특히 속도 분포(velocity profile)의 형태를 직관적으로 파악하는 데 유용하다.

와선(vortex line)은 유동장 내 각 점에서 와도 벡터(vorticity vector) $\boldsymbol{\omega} = \nabla \times \mathbf{V}$ 에 접하는 곡선으로 정의된다. 와선의 미분 방정식은 유선과 유사한 형태를 가진다:

$\frac{dx}{\omega_x} = \frac{dy}{\omega_y} = \frac{dz}{\omega_z}$

와선은 유동장의 회전 구조를 기술하며, 와관(vortex tube)과 함께 와류 역학(vortex dynamics)의 기하학적 기초를 형성한다.

3. 유선 해석의 수치적 방법

전산 유체역학에서 유선, 유적선, 유맥선의 수치적 계산은 속도장 데이터에 대한 수치 적분(numerical integration)을 통해 수행된다. 유선의 경우, 특정 시각에서의 속도장을 이용하여 유선 방정식의 공간 적분을 수행하며, 유적선의 경우에는 시간에 따른 속도장의 변화를 반영하여 시간 적분을 수행한다.

수치 적분에는 다양한 방법이 적용되며, 대표적으로 오일러 방법(Euler method), 2차 및 4차 룽게-쿠타 방법(Runge-Kutta method)이 사용된다. 4차 룽게-쿠타 방법은 정확도와 계산 효율의 균형이 우수하여 유동 곡선의 수치 계산에 널리 채택된다.

유맥선의 수치적 계산은 유선이나 유적선보다 복잡하다. 유맥선을 구성하려면 고정점을 서로 다른 시각에 통과한 다수의 유체 입자에 대해 각각의 유적선을 계산한 후, 관측 시각에서의 각 입자 위치를 연결하여야 한다. 이러한 계산 과정은 비정상 유동의 후처리(post-processing)에서 중요한 가시화 기법으로 활용된다.

4. 유동 가시화 기법과 실험적 적용

유동 가시화(flow visualization)는 유동장의 구조를 시각적으로 파악하기 위한 실험 기법의 총칭이다. 주요 유동 가시화 방법과 이들이 생성하는 유동 곡선의 대응 관계는 다음과 같다:

연기선(smoke wire) 기법: 가열된 금속선에서 발생하는 연기를 이용하여 유동 패턴을 가시화한다. 연속 방출 시 유맥선을, 순간 방출(pulse) 시 시간선에 근사적인 패턴을 생성한다.
염료 주입법(dye injection): 수조 실험에서 특정 지점에 염료를 주입하여 유맥선을 관찰한다.
입자 영상 유속계(Particle Image Velocimetry, PIV): 유동 내에 산포된 추적 입자의 변위를 단시간 노출 영상 쌍으로부터 측정하여 순간 속도장을 획득한다. 획득된 속도장으로부터 유선을 계산할 수 있다.
입자 추적 유속계(Particle Tracking Velocimetry, PTV): 개별 추적 입자의 궤적을 추적하여 유적선을 직접 측정한다.
쉴리렌 기법(Schlieren method): 밀도 구배에 의한 굴절률 변화를 이용하여 압축성 유동의 구조(충격파, 팽창파 등)를 가시화한다. 이 기법은 직접적으로 유동 곡선을 나타내지는 않으나, 유동장의 밀도 분포 구조를 파악하는 데 활용된다.

5. 로봇 공학에서의 유동 곡선 활용

비행 로봇의 공기역학적 설계와 해석에서 유동 곡선의 분석은 필수적인 과정이다. CFD 해석에서 유선 분포를 통해 익형 주위의 유동 구조, 박리 영역의 위치와 크기, 후류(wake)의 형태를 파악할 수 있으며, 이는 공력 성능의 정성적 평가에 활용된다.

멀티로터(multirotor) 시스템에서는 로터 하류(downwash)의 유선 분포가 각 로터 간의 공기역학적 간섭 특성을 분석하는 데 사용된다. 로터에 의해 유도된 하강 기류(induced downwash)의 유동 구조는 유선 해석을 통해 시각적으로 파악할 수 있으며, 이를 통해 로터 배치 최적화에 필요한 정보를 획득한다.

비정상 유동 환경에서 운용되는 비행 로봇의 경우, 돌풍(gust)이나 난류(turbulence) 조건에서의 유동장은 본질적으로 비정상적이므로, 유선과 유적선, 유맥선의 차이를 정확히 인식하는 것이 해석 결과의 올바른 해석에 중요하다. 특히 비정상 유동의 CFD 해석 결과를 가시화할 때, 순간 유선과 유적선을 구분하여 해석하여야 물리적으로 의미 있는 결론을 도출할 수 있다 (Tropea et al., 2007).

소형 비행 로봇의 저레이놀즈 수 유동 환경에서는 층류 박리 거품(laminar separation bubble)과 와류 구조가 유선 분포에 뚜렷이 나타나며, 이를 통해 박리 지점, 천이(transition) 위치, 재부착(reattachment) 지점 등의 핵심적인 유동 특성을 식별할 수 있다.

참고 문헌

Anderson, J. D. (2017). Fundamentals of Aerodynamics (6th ed.). McGraw-Hill Education.
Batchelor, G. K. (1967). An Introduction to Fluid Dynamics. Cambridge University Press.
Hama, F. R. (1962). Streaklines in a perturbed shear flow. Physics of Fluids, 5(6), 644–650.
Tropea, C., Yarin, A. L., & Foss, J. F. (Eds.). (2007). Springer Handbook of Experimental Fluid Mechanics. Springer.

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