2025-11-13T23:07:10.449786

Spatiotemporal statistics of the dissipation rate at the boundary of a turbulent flow using Diffusing-Wave Spectroscopy

Francisco, Lambret, AumaÃ®tre

We use Diffusing Wave Spectroscopy (DWS) to perform the first direct space- and time-resolved measurement of the dissipation rate~$Îµ$ at the boundary of a turbulent flow. We have shown in a previous publication that this technique provides maps of the dissipation rate of Newtonian fluids~\cite{Francisco}. Here, we apply the technique at the boundary of a turbulent flow generated in a square box by an impeller stirring the fluids. Although the measurement is made on a small region near the boundary, we show that the dissipation remains proportional to the injected power and follows the turbulent scaling $Îµ\propto \mathrm{Re}^3$, with Re being the Reynolds number ranging from $1.5 \times 10^4$ to $6 \times 10^5$. With this flow, there is no need for logarithmic corrections to reproduce the dissipation near the flat boundary. In addition, our setup allows us to measure the spatio-temporal fluctuations of the dissipation near the boundary. These fluctuations are quite large (the relative fluctuations are about 50\%) and are well described by a log-normal distribution, as expected for the dissipation rate in the bulk of homogeneous and isotropic turbulence (HIT) but Power Density Spectra (PDS) do not correspond to those expected for HIT \cite{Li07,Graham16,K62}

academic

확산파 분광법을 이용한 난류 흐름 경계에서의 소산율의 시공간 통계

기본 정보

논문 ID: 2511.03462
제목: Spatiotemporal statistics of the dissipation rate at the boundary of a turbulent flow using Diffusing-Wave Spectroscopy
저자: Enzo Francisco, Julien Lambret, Sébastien Aumâıtre
소속: Service de Physique de l'État Condensé, Université Paris-Saclay, CNRS, CEA, France
분류: physics.flu-dyn (유체역학)
제출일: 2025년 11월 5일
논문 링크: https://arxiv.org/abs/2511.03462

초록

본 연구는 확산파 분광법(DWS)을 이용하여 난류 경계에서 소산율 ε의 직접적인 시공간 분해 측정을 처음으로 달성했습니다. 연구팀은 정사각형 용기에서 임펠러 교반을 통해 난류를 생성하고 경계 근처에서 측정을 수행했습니다. 측정 영역이 작음에도 불구하고, 결과는 소산율이 주입 전력에 정비례하며 난류 스케일링 관계 ε∝Re³(레이놀즈 수 범위: 1.5×10⁴~6×10⁵)를 따름을 보여줍니다. 관로 흐름과 달리, 이 흐름은 평탄한 경계 근처에서 로그 보정이 필요하지 않습니다. 또한 경계 근처의 소산 시공간 변동이 매우 크며(상대 변동 약 50%), 로그정규분포를 따르지만, 전력 스펙트럼 밀도는 균질 등방성 난류(HIT)의 예상과 일치하지 않습니다.

연구 배경 및 동기

연구 문제

난류 소산율은 유체역학의 핵심 물리량이지만, 특히 경계 근처에서의 직접 측정은 극히 어렵습니다. 본 연구는 난류 경계에서 소산율의 시공간 분포를 직접 측정할 수 있는 실험 기술 개발을 목표로 합니다.

문제의 중요성

공학적 응용 가치: 소산율은 차량의 에너지 소비와 관로 내 압력 손실을 결정합니다
지구물리학적 의의: 기후 모델에서 에너지 균형과 전달 효율에 영향을 미칩니다
기초 과학적 의의: 소산 구조는 난류 간헐성을 주도할 수 있습니다
경계 문제: 많은 실제 응용에서 표면 근처의 소산과 전단율에 관심을 가집니다(침수 물체의 항력, 지구물리 흐름의 침식 효율 등)

기존 방법의 한계

열선 기술: Taylor 동결 난류 가정이 필요하며, 다중 선 프로브는 흐름을 교란합니다
입자 영상 속도 측정(PIV):
- 정확도가 질의 창 격자 크기에 의해 제한됨
- 중간 레이놀즈 수 범위에만 적용 가능
- 이중 평면 스테레오 PIV 기술이 복잡함
직접 수치 모의(DNS): 계산 비용이 높고, 장시간 모의가 어려우며, 드문 사건을 과소평가할 수 있습니다

연구 동기

뉴턴 유체의 소산율은 다음과 같이 정의됩니다: $\epsilon = \sum_{i,j} \frac{\nu}{2}(\partial_j u_i + \partial_i u_j)^2$

이는 모든 방향에서 속도장의 공간 도함수, 즉 변형률 텐서의 노름을 측정해야 합니다. DWS 기술은 이 문제를 해결하기 위한 새로운 경로를 제공합니다.

핵심 기여

최초 달성: 난류 경계에서 소산율의 직접적인 시공간 분해 측정을 최초로 수행
스케일링 관계 검증: 경계 소산율이 주입 전력에 정비례하며 ε∝Re³ 스케일링 법칙을 따르고 로그 보정이 필요 없음을 증명
변동 특성 규명: 경계 소산 변동이 50%에 달하며 로그정규분포를 따름을 발견
HIT 대비 분석: 경계 소산의 통계 특성이 부분적으로 HIT 이론(PDF)과 일치하지만 전력 스펙트럼 특성은 다름을 발견
기술 확장: DWS 기술을 Taylor-Couette 흐름에서 완전 난류 흐름으로 성공적으로 적용
고 레이놀즈 수 측정: Re가 6×10⁵에 달하는 측정을 달성하여 기존 방법의 한계를 초과

방법 상세 설명

작업 정의

입력: 난류 경계에서의 산란광 강도 시계열
출력: 경계 근처 소산율 ε의 시공간 분포 맵
제약 조건:

입자 반경 r과 광파장 λ는 광 전송 평균 자유 경로 l*보다 훨씬 작음
기울기 특성 길이는 유체 내 광 경로보다 큼
시간 분해능은 흐름 진화 시간 척도보다 작아야 함

DWS 기술 원리

기본 방정식

상호 간섭 광이 작은 입자를 포함한 운동 유체를 조사하면 스페클 패턴을 생성합니다. 역산란광 강도의 자기상관함수는:

$g_2(\tau) = \frac{\langle I(t+\tau)I(t)\rangle}{\langle I(t+\tau)\rangle\langle I(t)\rangle} = \beta \exp\left(-\gamma\sqrt{6(\tau/\tau_o + \tau^2/\tau_V^2)}\right) + 1$

여기서:

β: 강도 대비도
γ: 광학 기하학과 경계 조건에만 의존하는 매개변수
τₒ = 1/(Dk²): 브라운 운동으로 인한 시간 척도
τᵥ: 유체 운동으로 인한 시간 척도

소산율 추출

τᵥ와 소산율의 관계:

$\tau_V = \sqrt{30}/\left(l^*k\sqrt{\langle\epsilon\rangle/\nu}\right)$

τᵥ ∝ λ/l* ≪ 1/√(⟨ε⟩/ν)이므로, ⟨ε⟩의 시간 진화를 포착할 수 있습니다.

실험 장치

흐름 시스템

용기: 60 cm 변의 정사각형 수조
임펠러: 반경 R=20 cm, 1.2 cm 두께 원판과 4개의 3.6 cm 높이 블레이드 포함
회전 속도 범위: Ω = 0.02-0.8 회전/초
레이놀즈 수 정의: Re = 2πUL/ν, 여기서 U = RΩ
레이놀즈 수 범위: 1.5×10⁴ ≤ Re ≤ 6×10⁵
유체: 170리터 탈이온수 + 1.7 kg TiO₂ 입자(질량 분율 1%)

示踪 입자 특성

평균 수력학적 반경: r = 200 nm
Stokes 수: Sk ≤ 0.1(수동 示踪 거동 보장)
침강 속도: < 10⁻⁷ μm/s(무시 가능)
전송 평균 자유 경로: l* = 80 μm

특성 길이 척도 비교(표 1 참조)

매개변수	수치 범위
파장 λ	0.532 μm
입자 반경 r	0.2 μm
전송 자유 경로 l*	80 μm
용기 크기 L	60 cm
임펠러 반경 R	20 cm
Taylor 미시 척도 Λₜ	5-0.8 mm
Kolmogorov 척도 ηₖ	440-30 μm
점성 저층 δᵥ	1mm-25 μm

광학 시스템

조명 시스템: 2W Nd:YAG 레이저(λ=532 nm), 20× 현미경 렌즈를 통한 빔 확장
두 가지 감지 방식:
- PMT + 상관기: 직경 12 cm 원판 영역의 평균 소산 감지
- 고속 카메라: 5.1×5.1 cm² 영역(128×128 픽셀)의 변동 분해
시간 분해능:
- 상관기 샘플링 시간: 1.28 μs
- 카메라 프레임 속도: 430,000 fps
- 상관함수 수렴 필요: 25,000 프레임(0.06 s)

측정 프로토콜

정지 유체에서 g₂(τ) 측정, Stokes-Einstein 공식으로 알려진 τₒ로부터 γ 도출
임펠러를 설정 회전 속도로 실행, 방정식(3)에서 τᵥ 결정
정지 측정으로 γ 값 검증
대부분의 측정은 약 2초 지속(카메라 RAM 제한)
Re=2×10⁵에서 90초 장시간 측정 수행(45개의 2초 세그먼트)

기술 혁신점

공간 분해 장점: 변형률 텐서 추정 정확도가 카메라 분해능에 의존하지 않으며, 동일 구성에서 큰 레이놀즈 수 범위 탐색 가능
비침투성: 광학 측정으로 흐름장을 교란하지 않음
직접 측정: Taylor 가정이나 복잡한 속도장 재구성 불필요
시공간 분해: 공간 분포와 시간 진화 정보를 동시에 획득
고 레이놀즈 수 적용성: 기존 방법의 레이놀즈 수 한계 돌파

실험 설정

유체 시스템 매개변수

작동 유체: 탈이온수(동점성 ν ≈ 1×10⁻⁶ m²/s)
示踪 입자: TiO₂ 나노입자
- 밀도: ρₛ > ρw(침강 무시 가능)
- 농도: 1% 질량 분율
- 현탁 안정성: 탈이온수로 보장
전력 측정: 토크 측정기로 토크 C 측정, 주입 전력 P = C·Ω 계산

광학 매개변수 검증

l 현장 측정*: 공간 분해 반사 기술
- 다른 250 mW Nd:YAG 레이저 사용
- CCD 카메라로 확산광 반점의 방사형 강도 감쇠 측정
- ∝(r/l*)⁻³ 법칙 준수
- 48시간 실험 기간 동안 l* 안정 유지(연속 교반 조건)

측정 영역

PMT 감지 영역: 직경 12 cm 원판, 깊이 약 5l* ≈ 400 μm
카메라 감지 영역: 5.1×5.1 cm² 정사각형, 깊이 약 5l* ≈ 400 μm
픽셀 크기: 약 400 μm(l*보다 크며 요구사항 충족)
측정 위치: 임펠러 높이의 용기 벽면

실험 조건 제어

교반 전략: 고저 회전 속도 교대로 입자 균일 혼합 유지
측정 시간:
- 일반 측정: 약 2초
- 장시간 측정: 90초(PDF 통계용)
거시 시간 척도: 1/(4Ω) = 0.3-12.5초

실험 결과

주요 결과

1. 소산율 스케일링 관계(그림 2)

핵심 발견:

경계 소산율 ⟨ε⟩은 주입 전력 P에 정비례
난류 영점 법칙 준수: ⟨ε⟩ ∝ P ∝ (ν³/L⁴)·Re³
로그 보정 불필요: 관로/채널 흐름과 다름

정량적 결과:

측정된 경계 소산율은 평균 주입 전력의 약 2.4배
이러한 과잉 소산은 임펠러 근처 경계 과잉 소산의 예상과 일치

물리적 해석:

정사각형 용기의 단일 임펠러가 생성하는 흐름은 고도로 불균일
소산 구조의 공간 분포가 불균일
임펠러 근처 경계에 과잉 소산 존재

2. 공간 변동 특성(그림 3)

관찰 결과:

상대 변동: σε/⟨ε⟩ ≈ 50%(Re와 무관)
표준편차 스케일링: σε ∝ (ν³/L⁴)·Re³
소산 구조 척도가 픽셀 크기보다 훨씬 큼

구조 특성:

소산 구조가 대규모 흐름에 의해 평류 수송됨
구조가 늘어나고 변형됨
주요 흐름 방향: 환형 흐름(y 방향) + 상승 극 성분

시간 진화:

PMT 측정(12 cm 직경 평균): 변동 약 30%
블레이드 통과 주파수와 관련된 특성 주파수 미관찰

3. 확률 밀도 함수(그림 4)

통계 분석:

90초 장시간 측정(Re=2×10⁵)
변수: (log(ε) - ⟨log(ε)⟩)/σlog(ε)

주요 발견:

PDF가 로그정규분포 준수
HIT 이론 예상과 일치
그러나 측정 위치는 난류 경계층, 비-HIT 영역

이론 대비:

HIT의 로그정규분포는 캐스케이드 과정 기반
경계층에는 복잡한 계층 구조 존재
실험 결과는 로그정규성이 경계층에서도 성립함을 시사

4. 전력 스펙트럼 밀도(그림 5)

실험 스펙트럼 특성:

소규모 변동이 Re 증가에 따라 강화됨
픽셀 분해능이 모든 관련 구조를 포착하기에 충분
kηₖ = 1/ηₖ 이전에 노이즈 수준에 도달(관성 영역에서 가정)

HIT 이론과의 대비:

정밀 Kolmogorov 이론 예측: $\langle\epsilon(x)\epsilon(x+r)\rangle \propto \langle\epsilon^2\rangle(L/r)^\mu$ 해당 전력 스펙트럼: k⁻⁽¹⁻μ⁾, μ ≈ 0.2-0.3
실험 관찰: 스펙트럼 감쇠가 k⁻⁰·⁸보다 가파름

이방성 분석:

우선 방향 식별 어려움
방사형 파수 k = √(k²ₓ + k²y) 사용

DNS 검증(Johns Hopkins 데이터베이스)

모의 매개변수

흐름 유형: 난류 채널 흐름
레이놀즈 수: Re ≈ 4×10⁴
격자: 2048×1536×512 노드
계산 영역: 8π × 3π × 2
경계 조건: x, y 방향 주기, z 방향 무슬립
관찰 창: 20×40 노드(실험 5.1×5.1 cm²에 해당)

DNS 결과 대비

소산 구조(그림 3c, 3d):

실험과 DNS가 유사한 구조 표시
DNS에서 이방성이 더 명확함(평균 흐름 방향 명확)

PDF 대비(그림 4):

DNS 변동도 로그정규분포 준수
실험 결과와 일치

전력 스펙트럼 대비(그림 5 삽입):

횡 방향 및 흐름 방향 전력 스펙트럼
두 방향 모두 감쇠가 HIT 예측 k⁻⁰·⁸보다 가파름
실험 관찰과 일치

중요한 결론: DNS 결과는 실험 발견이 측정 오류가 아니라 경계 난류의 진정한 특성임을 증명합니다.

실험 발견 요약

스케일링 법칙 검증: 경계 소산율이 ε∝Re³를 준수하며 로그 보정 불필요
국소 과잉 소산: 임펠러 높이에서 소산율이 평균값의 2.4배
대폭 변동: 상대 변동 50%, Re와 무관
로그정규성: PDF가 HIT 이론 예상 준수
스펙트럼 특성 차이: PDS 감쇠가 HIT 예측보다 빠름
DNS 일관성: 채널 흐름 DNS가 경계 특성 증명

방법	장점	한계
열선 기술12	고주파 응답	Taylor 가정 필요; 다중 선 프로브 흐름 교란
점 측정 PIV13	직접 속도 측정	단일점 측정; 중간 Re
이중 평면 스테레오 PIV14	평면 분해	기술 복잡; 격자 정확도 제한
DNS11,15	완전한 정보	계산 비용; 드문 사건 과소평가
DWS(본 논문)	직접 측정; 고 Re; 시공간 분해	경계 측정; 투명 유체 필요

결론 및 논의

주요 결론

기술 돌파: DWS가 난류 경계 소산율의 직접적인 시공간 분해 측정을 최초로 달성하며, Re가 6×10⁵에 달할 때도 우수한 공간 및 시간 분해능 유지
스케일링 관계: 경계 소산율이 주입 전력에 정비례하며 순수 거듭제곱 스케일링 ε∝Re³를 따르고, 관로 흐름의 로그 보정 불필요
국소 특성: 임펠러 높이에서 측정 가능한 과잉 소산 존재(2.4배)
변동 특성:
- 소산 구조 변동이 강함(상대 변동 50%)
- 로그정규분포 준수(HIT와 일치)
- 전력 스펙트럼 감쇠가 HIT 예측보다 빠름
보편성: 채널 흐름 DNS가 측정 특성이 난류 경계의 보편적 거동을 나타냄을 증명

물리적 의의

경계 소산 메커니즘

임펠러 구동 흐름의 경계 소산이 고전적 벽 난류와 다름
체 소산 구조가 경계층 소산을 주도
점성 저층의 통계 가중치가 Re에 따라 변함(93.6%→2.9%)

통계 특성의 이중성

로그정규성 유지: 캐스케이드 과정의 영향이 경계까지 확장
스펙트럼 특성 변화: 경계층 복잡 구조가 공간 상관성 변경

한계

측정 위치 제한:
- 경계 근처 얇은 층만 측정(깊이~400 μm)
- 체 소산 직접 측정 불가
광학 조건 요구:
- 투명 또는 반투명 유체 필요
- 적절한 示踪 입자 필요(r, λ ≪ l*)
시간 분해능 제약:
- 상관함수 수렴에 0.06초 필요
- 흐름 진화 시간 척도가 충분히 길어야 함
- 연구 가능한 흐름 유형 제한
공간 분해능:
- 픽셀 크기~400 μm, 최소 척도보다 큼
- l*=80 μm이 Re 범위에서 Kolmogorov 척도 통과
- 고 Re에서 실패 가능(l* ≤ ηₖ)
흐름 유형:
- 본 논문 흐름이 고도로 불균일
- HIT 및 고전적 벽 난류와 모두 다름
- 보편성 추가 검증 필요
통계 수렴:
- 일반 측정은 2초만
- PDF 분석은 단 한 번의 90초 측정 기반
- 통계 수렴성 개선 필요