Evolution of passive scalar mixing layers in stratified and unstratified homogeneous turbulence

고해상도 와류 시뮬레이션은 성층화된 균질 난류에서 수동 스칼라 혼합이 횡방향에서는 비성층화 경우와 유사하게 거동하지만, 성층화가 대규모 교란을 억제함으로써 수직 혼합을 심각하게 저해하여 스칼라 플럭스에 대한 특정 모델링 접근법을 informing 하는 고유한 성장 한계와 변동 강도를 초래함을 보여줍니다.

핵심

잉크 한 방울이 물 한 컵을 통해 퍼지는 모습을 상상해 보세요. 물을 부드럽게 저으면 잉크는 고르게 퍼집니다. 이것이 과학자들이 유체 내 혼합을 일반적으로 생각하는 방식입니다: 난류는 거대한 숟가락처럼 작용하여 모든 것을 균일해질 때까지 함께 저어줍니다.

하지만 그 물이 층을 이루고 있다면 어떻게 될까요? 바닥의 물은 무겁고 짜고, 위의 물은 가볍고 싱겁다고 상상해 보세요. 이를 성층화라고 합니다. 실제 세계에서는 깊은 바닷물이 더 밀도가 높고, 대기는 위로 올라갈수록 얇아지고 가벼워지는 곳에서 발생합니다.

이 논문은 고도의 컴퓨터 실험으로, 이러한 층화가 난류 유체 내의 "오염"(예: 오염물질이나 연기) 이 퍼지는 방식을 어떻게 변화시키는지를 묻습니다.

여기서 그들이 발견한 내용을 간단한 개념으로 나누어 설명합니다.

설정: 두 가지 유형의 "오염"

연구자들은 처음에 혼란스럽게 소용돌이치는 (난류) 가상 유체를 만들었습니다. 그런 다음 두 가지 다른 "오염"(수동 스칼라) 을 도입하여 어떻게 행동하는지 관찰했습니다:

1. 수평 오염: 물에 떠 있는 평평한 종이처럼 옆으로 퍼진 잉크 층.
2. 수직 오염: 위아래로 퍼진 잉크 층, 즉 수직으로 세워진 색의 벽.

그들은 두 가지 시뮬레이션을 실행했습니다: 하나는 정상적인 물 (비성층화) 에서, 다른 하나는 "층화된" 물 (성층화) 에서입니다.

주요 발견: "위 - 아래"와 "옆 - 옆"의 차이

1. 옆 - 옆 오염 (횡단 층)

무슨 일이 일어났는가: 오염이 옆으로 퍼졌을 때, 층화는 그 확산을 막지 않았습니다. 오히려 층화된 물에서 정상적인 물보다 약간 더 빠르게 퍼졌습니다.
비유: 복도에서 사람들이 뛰는 상황을 상상해 보세요. 바닥이 완벽하게 평평하다면 (비성층화), 그들은 모든 방향으로 달립니다. 바닥에 보이지 않는 완만한 경사가 있다면 (성층화), 그들은 여전히 옆으로 잘 달립니다. 오히려 조금 더 활기차게 움직일 수도 있습니다. 두 경우 모두 "잉크"는 광범위하게 퍼집니다.
주의할 점: 전체적인 확산은 비슷했지만, 층화된 물의 "잉크"는 더 "뾰족"했습니다. 매끄러운 경사 대신 날카롭고 거친 가장자리를 가졌습니다. 더 "간헐적"이었으며,这意味着 순수한 잉크 주머니와 순수한 물 주머니가 존재하고 매끄러운 중간 지대는 적었습니다.

2. 위 - 아래 오염 (수직 층)

무슨 일이 일어났는가: 여기서 마법 (그리고 제한) 이 일어났습니다. 정상적인 물에서는 수직 오염이 옆으로 퍼지는 것처럼 위아래로 쉽게 퍼졌습니다. 하지만 층화된 물에서는 확산이 거의 완전히 멈췄습니다.
비유: 숟가락으로 두꺼운 밀크셰이크를 저으려 한다고 상상해 보세요. 숟가락을 위아래로 움직이려 하면 셰이크의 층들이 저항합니다. 무거운 것은 바닥에, 가벼운 것은 위에 머물고 싶어 합니다. "저어주는" 운동이 눌리게 됩니다.
결과: 수직 오염은 아주 처음에 약간 커졌지만, 그 후 "천장"에 부딪혔습니다. 유체의 안정된 층들이 뚜껑처럼 작용하여 난류가 수직으로 혼합되는 것을 막았기 때문에 더 이상 넓어질 수 없었습니다. 유체는 여전히 옆으로 소용돌이칠 수 있었지만, 위아래로 혼합할 수는 없었습니다.

왜 이것이 중요한가? ("무엇" 뒤에 있는 "왜")

연구자들은 수직 방향에서 유체가 용수철처럼 행동한다고 발견했습니다. 난류가 무거운 층을 위로 밀거나 가벼운 층을 아래로 밀려고 하면, 중력이 그것을 다시 끌어당깁니다. 이로 인해 "저어주는" 운동이 멈춥니다.

그러나 유체는 여전히 옆으로 소용돌이칠 수 있습니다. 따라서 난류의 "수직 길이"는 특정 크기로 고정됩니다 (중력과 유체 층의 강도에 의해 결정됨), 그리고 오염은 그 크기 이상으로 커질 수 없습니다.

예측을 위한 "레시피"

이 논문은 슈퍼컴퓨터 없이 이러한 오염이 어떻게 퍼질지 예측하기 위한 간단한 수학적 "레시피"를 만들려고 시도했습니다.

• 오염의 모양을 알고 있다면: 옆으로 퍼지는 속도를 예측하기 위해 간단한 한 숫자 공식을 사용할 수 있습니다. 매우 잘 작동합니다.
• 오염의 모양을 모른다면: 모양을 추측해야 합니다 (매끄러운 곡선처럼 보인다고 가정). 이렇게 하면 두 숫자 공식이 필요합니다. 이는 오염이 소용돌이치는 유체와 리듬을 맞추어 안정화된 이후에 매우 잘 작동합니다.

결론

• 옆으로의 혼합: 안정된 층들 (깊은 바다나 상부 대기처럼) 은 옆으로의 혼합을 막지 않습니다. 오히려 더 강렬하고 거칠게 만들 수도 있습니다.
• 위아래로의 혼합: 안정된 층들은 브레이크처럼 작용합니다. 유체가 위아래로 혼합되는 것을 거의 완전히 막습니다.
• "저어주기"와 "혼합"의 구분: 유체는 여전히 옆으로 "저어줄"(움직일) 수 있지만, 층들이 서로 바뀌는 것을 저항하기 때문에 위아래로 "혼합"(섞일) 수는 없습니다.

저자들은 그들의 실험이 특정 유형의 유체 특성 (프란틀 수 0.7) 을 사용했다고 지적합니다. 유체가 "더 두꺼우"거나 다른 특성 (프란틀 수 > 1) 을 가진다면, 혼합이 자체 부력을 생성하는 "역전" 효과로 인해 결과가 달라질 수 있다고 경고합니다. 하지만 그들이 테스트한 조건에서는 "옆으로는 자유롭고, 위아래는 차단된다"는 규칙이 유효합니다.

기술 요약

기술적 요약: 성층 및 비성층 균일 난류에서의 수동 스칼라 혼합층의 진화

문제 제기
상부 대류권, 성층권, 해양 혼합층과 같은 안정적으로 성층화된 환경에서 입자 및 스칼라의 확산은 대기 및 해양 과학에서 중요한 과제이다. 안정적 성층화가 밀도장과 난류 감쇠에 미치는 영향은 잘 문서화되어 있지만, 에어로졸, 오염물질, 미세플라스틱과 같은 수동 스칼라의 확산은 상대적으로 덜 이해되고 있다. 구체적으로, 안정화 밀도 구배의 영향 하에 감쇠하는 균일 등방성 난류 (HIT) 에 도입된 수동 스칼라 혼합층 (PSML) 이 어떻게 진화하는지는 명확하지 않다. 핵심 질문은 부력 규모와 이방성의 도입이 비성층 HIT 에서 관찰되는 스칼라 통계의 자기유사적 진화를 어떻게 변화시키는지에 관한 것이다.

방법론
저자들은 비성층 및 안정적으로 성층화된 균일 난류의 감쇠를 조사하기 위해 고해상도 대와류 시뮬레이션 (LES) 을 사용하였다. 시뮬레이션은 비수평압 근사 (non-hydrostatic Boussinesq approximation) 를 적용한 비발산 나비어 - 스토크스 방정식을 사용하였다.

• 유동 구성: 두 가지 시뮬레이션이 수행되었다. 첫 번째는 비성층 HIT 이었다. 두 번째는 초기 조건이 동일하지만 균일한 안정화 밀도 구배 ($g \neq 0$) 를 포함하여, 약 한 번의 부력 주기 후 유동이 강하게 성층화된 난류로 전환되도록 하였다.
• 수동 스칼라: 두 가지 수동 스칼라 혼합층이 초기화되었다. 하나는 수직 방향 ($\xi_z$) 에 평균 구배를, 다른 하나는 횡방향 ($\xi_y$) 에 평균 구배를 가진다. 이들은 혼합 분율 ($0 \le \xi \le 1$) 로 모델링되었다.
• 수치적 세부 사항: 자기유사성에서 벗어나게 하는 확산 효과를 최소화하기 위해 시뮬레이션은 동일한 계수를 가진 초점성 (hyperviscosity) 과 초확산 (hyperdiffusivity) 을 사용하였다. 분자 프란틀 수는 $Pr = 0.7$로 설정되었다. 레이놀즈 수는 높았다 (비성층 경우$Re_h > 15,000$및$Re_\lambda > 450$), 관성 - 대류 및 관성 - 확산 범위의 존재를 보장하였다. 해상도는 부력 길이 규모 ($L_b/\Delta \gg 1$) 를 해결하기에 충분하였다.
• 분석: 연구는 스칼라 프로파일 폭, 분산, 고차 모멘트 (왜도, 첨도), 플럭스, 그리고 혼합 효율의 진화를 추적하였다. 또한 평균 스칼라 프로파일이 알려진 경우와 가정해야 하는 경우 두 조건 하에서 횡방향 스칼라 플럭스에 대한 와동확산 계수 모델을 평가하였다.

주요 결과

1. 횡방향 혼합 ($\xi_y$): 중력에 수직인 방향에서 성층화된 경우의 수동 스칼라 진화는 비성층 경우와 광범위하게 유사하지만, 약간 더 빠르게 확산된다.

   • 스칼라 요동은 성층화된 경우 더 강렬하며, 피크 분산은 약 15% 더 높다.
   • 난류/비난류 경계면은 성층화된 경우 더 간헐적이다.
   • 스칼라장이 속도장과 준평형에 도달하면, 성층화 및 비성층 경우 모두 정규화된 플럭스 프로파일이 붕괴되어 성층화에도 불구하고 유사한 근본적인 수송 메커니즘을 나타낸다.

2. 수직 혼합 ($\xi_z$): 성층화는 수직 혼합을 극적으로 억제한다.

   • 수직 층 폭 ($\delta_z$) 은 초기에 성장하지만 약 한 번의 부력 주기 ($t^* \approx 3$) 후 급격히 확장을 멈춘다.
   • 이 성장 정지는 부력력에 의한 대규모 수직 교반의 억제로 귀결된다.
   • 층의 최종 폭은 수평 속도의 수직 적분 길이 규모 ($L_v$) 에 비례하며, 이는 수직 프루드 수를 1 차수 ($Fr_v \sim O(1)$) 로 유지하도록 제한된다.
   • 평균 프로파일은 분자 확산으로 인해 단조 증가를 보이지만, 플럭스 분석은 내부 파와 유사한 거동과 일치하는 스칼라 가닥을 원래 쪽으로 다시 쓸어가는 "역플럭스 (reverse flux)" 현상을 보여준다.

3. 혼합 효율: 성층화된 경우의 누적 혼합 (스칼라 분산의 소산) 은 후기 시간에 비성층 경우를 결국 초과한다. 이는 혼합이 교반을 따르며, 성층화된 유동에서의 수평 교반율이 수평 길이 규모에 비례하고, 이는 비성층 경우보다 더 크게 성장하기 때문이다.

4. 플럭스 모델링:

   • 알려진 프로파일: 평균 스칼라 프로파일이 알려져 있다면, 단일 상수 와동확산 계수 모델 ($F_y = -c_1 \delta_y v' \partial \xi / \partial y$) 은 $c_1 \approx 0.4$로 성층화 및 비성층 경우 모두 플럭스를 효과적으로 예측한다.
   • 가정된 프로파일: 프로파일 모양 (오차 함수로 가정됨) 과 스칼라 길이 규모를 속도 통계에서 추정해야 한다면, 두 상수 모델이 필요하다. 이 모델은 스칼라가 속도장과 준평형에 도달하여 스칼라 길이 규모를 운동 에너지 소산율에서 스케일링할 수 있는 후기 시간에 잘 작동한다.

의의 및 주장
본 논문은 이방성 확산에 대한 이론적 예측과 수치적 관찰 사이의 간극을 연결하는 감쇠 성층 난류에서의 수동 스칼라 혼합에 대한 상세한 특성을 제공한다고 주장한다.

• 이방성: 본 연구는 성층화가 혼합에 있어 극명한 이분법을 만든다는 것을 확인한다. 즉, 수평 혼합은 활발하게 유지되거나 심지어 강화되는 반면, 수직 혼합은 유동이 강하게 성층화되면 효과적으로 억제된다.
• 모델링 유용성: 저자들은 스칼라가 속도장과 준평형에 있다면 표준 와동확산 계수 접근법이 성층화된 유동에서도 유효함을 입증한다. 그들은 프로파일 모양을 가정하는 것보다 알려진 프로파일을 가정하는 것이 모델링 요구 사항을 크게 단순화한다고 강조한다.
• 한계 및 향후 작업: 저자들은 겸손하게 결과가 프란틀 수 0.7 에 구체적임을 지적한다. 최근 문헌은 더 높은 프란틀 수 ($Pr > 1$) 가 작은 규모에서 역부력 플럭스를 유발하여 유동 감쇠율과 결과적으로 수동 스칼라 혼합율을 변경할 수 있다고 명시적으로 밝히고 있다. 따라서 현재 발견은 추가 조사 없이는 고프란틀 수 시나리오에 직접 적용되지 않는다.

본 연구는 성층화가 혼합층의 수직 역학을 근본적으로 변화시키지만, 횡방향 혼합은 간헐성과 요동 강도가 증가했음에도 불구하고 비성층 난류의 많은 특성을 유지한다고 결론 내린다.