입자 영상 유속계

유체역학 전산유체역학 난류 마지막 수정: 2026-09-02 04:37:18

1. 개요[편집]

입자 영상 유속계 Particle Image Velocimetry (PIV)
재는 것	면(또는 체적) 전체의 순간 속도장
원리	추적입자 두 장 촬영 → 심문창별 상호상관 → 변위 ÷ $\Delta t$
광원	이중 펄스 레이저 시트 (대개 Nd:YAG 532 nm)
부화소 정밀도	~0.05 ~ 0.1 px — 가우시안 3점 피팅
설계 국룰	「1/4 규칙」 · 심문창당 입자상 10개 이상
천적	피크 로킹 · 면외 손실 · 강한 속도 구배

유동장 한 점을 아주 정확히 재는 것과, 유동장 전체를 한꺼번에 재는 것 중 무엇이 더 쓸모 있는가. 지난 30년의 답은 후자였다.

입자 영상 유속계(Particle Image Velocimetry, PIV)는 유동에 뿌린 미세 추적입자를 얇은 레이저 시트로 밝혀 아주 짧은 간격으로 두 번 촬영하고, 두 영상을 작은 창 단위로 상호상관해 얻은 국소 변위를 시간간격으로 나눠 순간 속도장을 복원하는 비접촉 전역 계측법이다. 유동에 물리적 탐침을 꽂지 않고, 한 점이 아니라 면 전체를, 시간 평균이 아니라 순간 값을 준다는 세 가지 성질이 이 기법을 실험 유체역학의 표준으로 만들었다.

의미는 계측 자체보다 그것이 열어 준 비교에 있다. 피토관과 열선유속계는 한 점의 시계열을 준다. 그런데 난류의 핵심 대상은 와류·전단층·재순환 같은 공간 구조이고, 전산유체역학 해석이 뱉는 것도 한 점이 아니라 장(field)이다. PIV가 나오기 전까지 해석 결과와 실험을 맞대는 일은 몇 개의 점과 몇 개의 적분량으로 하는 협상이었고, PIV 이후에는 장과 장을 직접 겹쳐 보는 일이 되었다. 검증 및 확인의 무기고에서 이 기법이 차지하는 자리가 그것이다.

2. 계측 사슬[편집]

2.1. 시딩 — 실은 여기서 승부가 난다[편집]

입자가 유동을 따라가지 못하면 그 뒤의 모든 계산은 정교한 헛수고다. 지름 $d_p$ , 밀도 $\rho_p$ 인 구형 입자의 스토크스 응답시간은

\tau_p = \frac{\rho_p\, d_p^{2}}{18\,\mu}

이고, 유동의 특성시간 $\tau_f$ 와의 비 $\mathrm{St} = \tau_p/\tau_f$ 가 스토크스 수다. $\mathrm{St}\ll1$ 이어야 입자가 유동을 충실히 따라간다. 난류에서 $\tau_f$ 는 콜모고로프 스케일의 시간척도로 잡는 것이 보수적인 기준이다.

여기서 정면충돌이 일어난다. 입자를 키우면 따라가지 못하고, 줄이면 보이지 않는다. 산란광 세기는 미 산란 영역에서 대략 지름의 제곱에 비례하므로(미 산란), 지름을 절반으로 줄이면 밝기가 1/4로 떨어진다. 그래서 실무는 늘 “따라가는 한도 안에서 최대한 큰 입자”에 수렴한다. 기체 유동에서는 1~2 μm 오일 미스트(DEHS)나 연기, 액체에서는 10 μm 안팎의 폴리아미드·유리 중공구가 관례다. 부력에 따른 침강도 같이 봐야 한다.

2.2. 조명과 촬영[편집]

이중 공진기 Nd:YAG 레이저가 나노초 펄스를 $\Delta t$ 간격으로 두 번 쏜다. 펄스가 나노초라는 점이 중요하다 — 그 시간 동안 입자는 사실상 정지해 있으므로 상이 흐르지 않고, 노출 시간이 아니라 펄스 간격이 시간분해능을 정한다. 시트 두께는 대개 0.5~2 mm이며, 이 두께가 뒤에 나올 면외 손실을 지배한다.

카메라는 프레임 스트래들링(frame straddling)으로 동작한다. 첫 프레임의 노출 끝자락에 첫 펄스를, 두 번째 프레임 시작 직후에 두 번째 펄스를 넣어, 카메라의 프레임률과 무관하게 마이크로초 단위의 $\Delta t$ 를 만든다. PIV 시스템에서 가장 비싼 부품이 종종 레이저가 아니라 이 타이밍을 맞추는 동기화 장치인 이유다.

기록되는 입자상의 지름도 설계 변수다. 기하학적 상과 회절 한계의 합성으로

d_\tau \simeq \sqrt{\,(M d_p)^2 + \bigl(2.44\,\lambda\,f_\#\,(M+1)\bigr)^2\,}

정도가 되는데( $M$ 은 배율, $f_\#$ 는 조리개값), 목표는 2 px 안팎이다. 이유는 아래 피크 로킹 절에서 나온다. 입자상이 너무 작으면 조리개를 조여 회절을 키우거나 살짝 초점을 흐리는, 처음 들으면 어이없는 처방이 실제로 쓰인다.

3. 상호상관 — 심문창이라는 단위[편집]

영상 전체를 한꺼번에 다루지 않는다. 한 변 $W$ (전통적으로 32 또는 16 px)의 심문창(interrogation window)으로 잘라, 창 안 입자들이 한 덩어리로 같은 변위만큼 이동했다고 가정하고 그 변위 하나를 찾는다. 두 시각의 창 $I_1, I_2$ 에 대해

R(m,n) \;=\; \sum_{i}\sum_{j} I_1(i,j)\, I_2(i+m,\, j+n)

를 모든 이동량에 대해 계산하고 봉우리 위치를 변위로 삼는다. 직접 계산은 $O(W^4)$ 지만, 상관 정리에 따라 고속 푸리에 변환으로 $R = \mathcal F^{-1}\{\overline{\mathcal F(I_1)}\,\mathcal F(I_2)\}$ 를 쓰면 $O(W^2\log W)$ 로 떨어진다. 1990년대에 PIV가 연구실 밖으로 나간 것은 이 한 줄 덕이 크다.

FFT를 쓰면 상관이 순환(circular)이 되어 창 밖으로 밀려난 것이 반대편에서 되돌아오는 접힘이 생기므로 영채움(zero padding)을 한다. 또 하나 더 실질적인 편향이 있다 — 변위가 클수록 두 창이 겹치는 면적이 줄어 상관값 자체가 작아지고, 그 결과 봉우리가 원점 쪽으로 끌린다. 겹침 면적(삼각 가중치)으로 나눠 정규화하는 것이 표준 처방이다.

봉우리를 무너뜨리는 요인은 셋이다. 면내 손실(입자가 창 밖으로 나감), 면외 손실(입자가 레이저 시트 밖으로 나감), 입자상 밀도 부족. 셋 다 신호 봉우리를 낮추고 잡음 봉우리를 상대적으로 키워, 어느 순간 상관이 엉뚱한 곳을 가리킨다.

창을 겹쳐서(50 % 또는 75 % 중첩) 격자를 촘촘하게 뽑는 것이 관례인데, 중첩으로 늘어난 벡터는 독립이 아니다. 공간분해능은 어디까지나 창 크기가 정하고 중첩은 보간에 가깝다. 이 구분을 흐리면 “벡터가 많으니 분해능이 높다”는 착각으로 이어지고, 난류 통계에서 소용돌이 크기를 과대평가하게 된다.

진리값 유동장(램-오신 소용돌이 + 균일 이송)으로 만든 합성 추적입자 이미지쌍만 보고 16×16 심문창 직접 상호상관과 3점 가우시안 부화소 피팅으로 벡터장을 복원하고, 복원 변위의 소수부를 함께 쌓는다. 입자 지름 0.8 px 에서는 소수부의 87.8%가 정수 ±0.05 안에 몰리지만(균일이면 10%) 2.2 px 에서 9.9%로 풀린다 — 피크 로킹. 화소값은 화소 중심의 점 표본이다.

4. 부화소와 피크 로킹[편집]

정수 픽셀 단위 변위로는 쓸 만한 정밀도가 안 나온다. 가우시안 모양의 입자상들을 상관하면 봉우리도 가우시안이므로, 봉우리 화소와 그 양옆 두 점에 가우시안을 맞춘다. 로그를 취하면 포물선 피팅과 같아져 닫힌 식이 나온다.

\varepsilon \;=\; \frac{\ln R_{-1} - \ln R_{+1}}{2\ln R_{-1} - 4\ln R_{0} + 2\ln R_{+1}}

이 세 점 피팅 하나로 변위 정밀도가 약 0.05~0.1 px 까지 내려간다. 32 px 창에서 8 px를 움직인다면 상대오차 1 % 수준이라는 뜻이다.

그런데 여기에 PIV 최대의 체계 오차가 숨어 있다. 피크 로킹(peak locking)은 추정된 변위가 정수 픽셀 값 쪽으로 몰리는 편향이다. 입자상이 1 px보다 작으면 봉우리가 사실상 한 화소에만 얹혀 세 점 피팅이 형태 정보를 얻지 못하고, 추정치가 가장 가까운 정수로 빨려 들어간다. 앞에서 입자상 지름을 2 px 근처로 맞추라고 한 이유가 이것이다.

진단은 간단하다. 변위의 소수부 히스토그램을 그려 보라. 평평해야 정상이고, 0과 1 근처에 뾰족한 봉우리가 서면 피크 로킹이다. 이것이 무서운 이유는 평균 속도장은 멀쩡해 보이는데 변동 통계가 오염되기 때문이다. 난류강도가 작은 영역에서 변위 변동이 통째로 정수로 붙잡히면 분산이 인공적으로 줄고, 레이놀즈 응력이 실제보다 작게 나온다. 그리고 그 결과를 난류 모델링 검증에 쓰면, 모델이 틀린 게 아니라 계측이 틀린 것을 가지고 몇 달을 싸우게 된다.¹

5. 설계 규칙 — 「1/4 규칙」과 그 감옥[편집]

킨과 에이드리언이 몬테카를로 모사로 정리한 설계 지침이 오래도록 국룰이었다. 요지는 넷이다.

심문창 하나당 입자상이 10개 이상 들어와야 상관 봉우리가 안정적으로 선다.
면내 변위 ≤ 심문창 한 변의 1/4. 32 px 창이면 8 px까지.
면외 변위 ≤ 레이저 시트 두께의 1/4.
창 안의 속도 변화가 입자상 지름을 넘지 않을 것. 넘으면 봉우리가 변위 방향으로 늘어져 뭉개진다.

2번과 3번을 묶어 1/4 규칙이라 부른다. 문제는 이 규칙이 만드는 감옥이다. 최대 변위가 $0.25W$ 로 묶이고 최소 분해 변위가 부화소 정밀도 $\sim0.1\,\text{px}$ 라면, 한 장에서 표현 가능한 속도 동적 범위는

\mathrm{DVR} \;\approx\; \frac{0.25\,W}{0.1\ \text{px}} \;\xrightarrow{\;W=32\;} \; \approx 80

에 불과하다. 재순환 영역과 주류가 함께 있는 유동에서는 한쪽을 포기해야 한다는 뜻이고, $\Delta t$ 를 고르는 일이 늘 고통스러운 이유가 여기 있다. 창을 키우면 동적 범위는 늘지만 공간분해능이 죽는다. 정면돌파는 불가능하고, 아래의 반복법이 우회로다.

6. 다중격자와 창 변형 — 규칙을 우회하는 법[편집]

1990년대 말 이후의 표준은 반복적 다중격자 + 창 변형이다.

큰 창(예: 64 px)으로 한 번 상관해 거친 변위장을 얻는다.
그 변위만큼 두 번째 창을 미리 옮겨 놓고(window shifting) 다시 상관한다. 창이 이미 따라가 있으므로 잔여 변위가 거의 0이고, 면내 손실이 사라진다. 그러면 1/4 규칙의 2번이 최종 창 크기가 아니라 잔여 변위에만 걸리게 되어, 창을 16 px, 8 px까지 줄일 수 있다.
여기에 국소 속도 구배까지 반영해 창을 변형시킨다(window deformation). 전단층 안에서 사각형 창을 평행사변형으로 기울여 잡는 것이라, 규칙 4번의 봉우리 뭉개짐을 직접 공략한다. 영상을 부화소로 워핑해야 하므로 B-스플라인 같은 고차 보간이 필요하다.
격자를 세분하며 2~4회 반복한다.

효과는 극적이다. 공간분해능이 몇 배 좋아지고 동적 범위도 함께 오른다. 대신 반복이 불안정해질 수 있다 — 예측 변위장에 창 크기에 가까운 파장의 잡음이 실리면 다음 반복에서 증폭되어 격자무늬로 발산한다. 그래서 반복 사이에 예측장을 저역통과 필터로 매끄럽게 만드는 처리가 사실상 필수이고, 이 평활화의 세기가 “분해능을 얼마나 가져갈 것인가”를 정하는 숨은 손잡이가 된다.²

7. 벡터 검증 — 틀린 벡터를 골라내기[편집]

상관이 잡음 봉우리를 골랐을 때 나오는 가짜 벡터(spurious vector)는 크기와 방향이 이웃과 전혀 다르므로 눈에 잘 띈다. 자동 검출의 표준 도구는 둘이다.

봉우리 비(peak ratio) $Q = R_1/R_2$ — 최고 봉우리와 두 번째 봉우리의 비다. 승부가 애매하면 낮아지므로 대략 $Q < 1.2\!\sim\!1.5$ 를 버린다. 스테레오 비전의 유일성 비율과 완전히 같은 발상.
정규화 중앙값 검사 — 이웃 8개의 중앙값에서 얼마나 벗어났는지를, 그 이웃들의 중앙값 절대편차로 나눠 무차원화한 잔차를 본다. 여기에 부화소 잡음 수준에 해당하는 작은 상수를 더해 두어 유동이 균일한 곳에서 분모가 0으로 죽는 것을 막는다. 평균 대신 중앙값을 쓰는 이유는 로버스트 통계 그대로다 — 이웃 안에 이미 가짜 벡터가 하나 섞여 있어도 판정이 흔들리지 않아야 한다.

버린 자리는 이웃 보간으로 메우고, 가짜 벡터 비율 자체를 품질 지표로 보고한다. 관례적으로 5 %를 넘으면 실험 설정을 다시 본다. 대치된 벡터가 20 %인 예쁜 벡터장은 예쁜 보간 결과일 뿐이라는 사실을, 컬러 컨투어를 뽑고 나면 잊기 쉽다.

8. PIV·PTV·LDV[편집]

	PIV	PTV	LDV
관점	오일러 — 격자 위 속도장	라그랑주 — 입자별 궤적	오일러 — 한 점
공간 범위	면 또는 체적	면 또는 체적	점 하나(트래버스 필요)
시간분해능	레이저·카메라가 제한(고속형은 수 kHz)	같음	매우 높음(입자 통과마다)
시딩 농도	창당 10개 이상, 조밀	개별 식별이 되도록 희박	희박해도 됨
공간분해능	심문창 크기가 하한	입자 간격	측정체적 크기
잘하는 것	순간 구조, 통계, CFD 대조	궤적·가속도·라그랑주 통계	벽 근처 고정밀 시계열, 스펙트럼

PTV(입자 추적 유속계)는 개별 입자를 프레임마다 식별해 따라간다. 상관 평균이 없으니 심문창이 강제하는 공간분해능 하한이 사라지고 가속도 같은 라그랑주 양을 직접 얻지만, 입자가 조밀해지면 대응 문제가 폭발한다. 최근에는 궤적의 시간 연속성을 예측자로 써서 조밀한 시딩에서도 개별 추적을 해내는 계열이 자리 잡아, PIV와 PTV의 경계가 흐려지고 있다.

LDV(레이저 도플러 유속계)는 원리부터 다르다. 두 레이저 빔을 교차시켜 만든 간격 $d_f$ 의 간섭 무늬를 입자가 통과할 때 산란광이 $f_D = u_\perp/d_f$ 로 깜빡이는데, 이 주파수를 재면 속도가 나온다. 영상도 상관도 없이 주파수 계측이므로 교정이 기하학만으로 끝나고(무늬 간격은 파장과 교차각으로 결정된다) 정확도가 매우 높다. 대신 한 점씩이다. 벽 근처 마찰속도나 난류 스펙트럼처럼 한 점의 시계열이 필요한 곳에서 LDV는 여전히 대체 불가다.

9. 3차원으로, 그리고 시뮬레이션과의 접점[편집]

스테레오 PIV — 카메라 두 대로 같은 면을 다른 각도에서 찍어 면외 속도 성분까지 얻는다(2D3C). 비스듬히 보므로 상면·렌즈면·물체면을 기울여 초점을 맞추는 샤임플러그 배치가 들어가고, 두 시점을 묶는 것은 스테레오 비전의 기하 그대로다. 덤으로 면외 성분을 알게 되어 면외 손실을 직접 진단할 수 있다.
단층 PIV(tomographic PIV) — 시트 대신 두꺼운 체적을 밝히고 카메라 3~6대로 찍어 입자 분포를 체적으로 재구성한 뒤 3차원 상관을 건다. 재구성은 전형적인 부정문제이고 대수적 반복법(MART 계열)이 쓰인다. 얻는 것은 3차원 3성분 속도장, 따라서 와도 전체와 속도구배 텐서 전체다.
압력장 복원 — 시간분해 PIV로 속도장의 시간미분을 얻으면 운동량 방정식에서 압력 구배를 계산하고 적분해 압력장을 만들 수 있다. 비접촉으로 압력을 재는 셈이라 매력적이지만, 미분과 적분을 데이터에 연달아 거는 구조라 잡음에 극도로 민감하다.
데이터 축소와 융합 — PIV 스냅숏 행렬에 적합직교분해나 동적 모드 분해를 걸어 지배 구조를 뽑는 것이 표준 후처리가 됐다. 시야 밖이나 가려진 영역의 결측을 저계수 구조로 메우는 갭피 POD는 행렬 완성의 문제 그대로이며, 계측장을 해석에 사전정보로 밀어 넣는 자료동화형 접근도 활발하다.
검증 데이터로서의 무게 — 직접수치모사나 대와류 모사 결과와 PIV를 겹칠 때는 필터 폭을 맞춰야 한다. PIV의 심문창은 물리적으로 공간 저역통과 필터이므로, 필터링하지 않은 DNS 장과 직접 비교하면 계측이 작은 스케일을 놓친 것이 마치 해석 오차처럼 보인다. 계측기의 전달함수를 해석 쪽에 걸어 주는 것까지가 공정한 비교다.³

10. 관련 문서[편집]

소수부 히스토그램을 그려 보는 데 30초가 걸린다. 그 30초를 아끼고 논문을 쓴 뒤 심사에서 “peak locking을 확인했는가”라는 한 줄 코멘트를 받는 것이 이 바닥의 통과의례다. 참고로 이 편향은 평균장에는 잘 드러나지 않아서, 평균 속도 프로파일이 문헌과 잘 맞는다는 사실은 아무런 면죄부가 되지 못한다. ↩
“반복하면 좋아진다”는 직관이 배신하는 대표적 사례다. 반복이 예측장을 되먹임으로 쓰는 구조이므로, 평활화가 부족하면 창 크기 근처 파장이 지수적으로 자라 벡터장에 체스판 무늬가 뜬다. 상용 소프트웨어가 반복 횟수와 평활 세기를 사용자에게 다 열어 주지 않는 데에는 나름의 이유가 있다. ↩
그래서 PIV와 해석이 안 맞을 때 던져야 할 첫 질문은 “모델이 틀렸나”가 아니라 “우리가 같은 스케일을 보고 있나”다. 심문창 32 px가 물리적으로 몇 mm인지, 그것이 그 유동의 테일러 스케일 대비 얼마인지를 계산해 보면 논쟁의 절반이 그 자리에서 끝난다. ↩