Two pathways to diapycnal mixing in strongly stratified flows with no initial vertical shear

이 연구는 선형 이론과 직접 수치 시뮬레이션을 결합하여, 초기 수직 전단이 없는 강한 성층 흐름에서 수평 전단 불안정성이 두 가지 뚜렷한 경로—수직 전전의 직접적인 발생 또는 기둥 모양 와류로의 비선형적 진화 중 하나를 통해—를 거쳐 필연적으로 연직 확산을 유도하며, 이 두 경로 모두 궁극적으로 소규모 켈빈-헬름홀츠 불안정성을 촉발하지만 서로 다른 연직 규모를 자극하기 때문에 서로 다른 혼합 효율을 나타낸다는 것을 밝혀냈다.

핵심

바다나 대기를 거대한 층이 있는 케이크라고 상상해 보세요. 이 층들은 서로 다른 밀도를 가진 유체(마치 서로 다른 맛의 케이크 층처럼)로 이루어져 있으며, 쉽게 섞이지 않습니다. 보통 과학자들은 이 층들을 옆으로 밀면서 동시에 위아래로도 밀 때(수직 전단) 어떤 일이 일어나는지 연구합니다. 하지만 이 논문은 다른 질문을 던집니다. 만약 매우 강력하게 층이 나뉜 환경에서, 위아래 움직임 없이 오직 옆으로만 밀기만 한다면 어떤 일이 벌어질까요?

연구진은 비록 처음에 수평적인 흐름만 있더라도, 자연은 혼돈(난류)과 층의 혼합을 만들어내는 두 가지 뚜렷한 "레시피" 또는 경로를 가지고 있다는 것을 발견했습니다. 자연이 어떤 레시피를 선택하느냐는 전적으로 실험에 어떤 "씨앗"을 심느냐, 즉 실험 시작 시점에 유체에 아주 작은 자극을 어떻게 주느냐에 달려 있습니다.

다음은 두 가지 경로를 쉬운 비유를 통해 설명한 내용입니다.

설정: 잔잔한 강

넓고 잔잔한 강이 수평으로 흐르고 있다고 상상해 보세요. 물은 층층이 쌓인 팬케이크처럼 층이 나뉘어 있습니다(강한 성층). 처음에는 흐름이 매끄럽고 2차원적입니다(오직 좌우로만 움직이며 위아래로는 움직이지 않습니다).

경로 1: "붐비는 방" 효과 (직접적인 경로)

시작: 강 전체에 아주 작고 무작위적인 자극을 동시에 줍니다(마치 공중에 꽃가루 한 줌을 뿌리는 것처럼).
과정:

1. 물결: 층 때문에 유체는 단순히 좌우로 출렁이는 것이 아니라, 즉시 다양한 크기로 위아래로 출렁이기 시작합니다. 이것은 마치 방 안의 사람들이 모두 동시에 움직이려 하여 혼란스러운 다방향적 요동을 일으키는 것과 같습니다.
2. 전단: 이 물결들은 매우 빠르게 강력한 수직 흐름(전단)을 만들어냅니다.
3. 붕괴: 이 수직 흐름들이 너무 강해져서 끊어지며, 작은 소용돌이(마치 작은 소용돌이 같은 것)를 생성합니다. 이것이 바로 바람이 불 때 물 위에서 보이는 부서지는 파도와 같은 "켈빈-헬름홀츠(Kelvin-Helmholtz)" 불안정성입니다.
   결과: 혼합이 효율적으로 일어납니다. 에너지가 다양한 크기의 물결에 분산되어 있기 때문에 "마찰"(점성 소산)이 낮아져 혼합 과정이 상대적으로 효율적입니다.

경로 2: "싱크로나이즈드 댄스" 효과 (간접적인 경로)

시작: 강에 매우 구체적이고 조직적인 자극을 줍니다(마치 지휘자가 지휘봉을 휘둘러 모든 사람이 특정 패턴에 맞춰 움직이게 하는 것처럼).
과정:

1. 와동: 무작위적인 물결 대신, 유체는 수직으로 길게 늘어선 소용돌이 기둥(강 속에 서 있는 거대한 토네이도 같은 것)으로 스스로를 조직화합니다. 한동안 흐름은 이 거대한 소용돌이 기둥들만 있는 완벽한 2차원 상태를 유지합니다.
2. 흔들림: 결국 이 거대한 기둥들이 불안정해집니다. 이들은 매우 특정한 고주파 방식으로 흔들리기 시작합니다. 연구진은 이를 "쌍곡선 불안정성(hyperbolic instability)"이라고 부릅니다. 이는 마치 팽이가 쓰러지기 직전에 격렬하게 흔들리는 모습과 같습니다.
3. 붕괴: 이 격렬한 흔들림은 매우 얇고 날카로운 수직 전단 층을 만듭니다. 이 얇은 층들은 경로 1에서와 마찬가지로 작은 소용돌이로 부서집니다.
   결과: 혼합은 일어나지만, 덜 효율적입니다. 왜일까요? 이 경로는 극도로 얇고 날카로운 층을 만들기 때문입니다. 이 작은 층들을 만들고 깨뜨리는 데는 많은 에너지(마찰)가 필요합니다. 이것은 마치 두꺼운 치즈를 둔한 칼로 자르는 것(경로 1)과 면도날로 자르는 것(경로 2)의 차이와 같습니다. 면도날은 훨씬 더 날카롭고 에너지를 많이 소모하는 절단을 만들어냅니다.

핵심 결론

이 논문은 수직 전단(위아래 움직임)이 처음에 반드시 존재할 필요는 없다는 것을 증명합니다. 유체가 충분히 두껍다면(높은 레이놀즈 수), 수직 전단은 수평 전단의 필연적인 부산물입니다.

• 무작위 노이즈로 시작하면: 경로 1(직접적이고 효율적인 혼합)이 나타납니다.
• 특정한 패턴으로 시작하면: 경로 2(간접적이고 덜 효율적인 혼합)가 나타납니다.

연구진은 강력한 컴퓨터 시뮬레이션을 사용하여 이 두 경로가 실재하며 서로 구별된다는 것, 그리고 처음에 선택한 "레시피"가 열로 낭비되는 에너지와 실제로 층을 섞는 데 사용되는 에너지의 양을 결정한다는 것을 보여주었습니다.

요약하자면: 완벽하게 층이 나뉜 유체라 할지라도, 수평적인 밀기(push)는 결국 수직적인 혼돈을 만들어냅니다. 하지만 그 밀기를 어떻게 시작하느냐에 따라 그 혼돈의 모습은 달라지며, 층을 섞는 효율성 또한 달라집니다.

기술 요약

문제 정의
대규모의 강한 성층화된 수평 흐름이 3차원(3D) 섭동으로 전이되는 안정성은 지구물리학 및 천체물리학 시스템에서 에너지의 하향 스케일 전이와 이성분 확산(diapycnel) 수송을 결정하는 핵심 기제이다. 수직 전단이 존재하는 성층화 흐름은 잘 연구되어 왔으나, 수직 전단이 없는 수평 전단 흐름에 대해서는 상대적으로 연구가 부족했다. 핵심적인 질문은 다음과 같다: 초기에 2차원(2D)이며 강하게 성층화된($Fr \ll 1$) 수평 전단 흐름이 어떻게 난류로 전이되고 수직 전단을 생성하는가? 기존 문헌은 두 가지 뚜렷한 경로를 제시한다: 하나는 3D 일차 모드의 직접적인 성장으로 인해 케엘빈-헬름홀츠(Kelvin-Helmholtz, KH) 불안정성을 유도하는 경로이고, 다른 하나는 기둥형 와류(columnar vortices)로의 비선형 진화 후 이차 불안정성으로 이어지는 경로이다. 그러나 이러한 경로들이 초기 조건(예: 하이퍼볼릭 탄젠트 vs 사인파 프로파일)의 차이로 인해 서로 배타적인 것인지, 아니면 문제 설정(초기값 문제 vs 체적 강제 방식)의 차이인지, 혹은 동일한 물리적 시스템 내에서 이용 가능한 대안적 경로인지는 불분명했다. 또한, 선택된 경로가 이성분 혼합 효율에 미치는 영향이 체계적으로 정량화되지 않았다.

방법론
저자들은 삼중 주기 도메인에서 이러한 전이를 조사하기 위해 선형 안정성 이론과 직접 수치 시뮬레이션(DNS)을 결합하였다. 연구는 선형 성층화된 유체 내에서 일정한 부력 주파수 $N$을 갖는 수직 불변, 수평 평행 콜모고로프 흐름($\bar{u}(y) = \sin(y)\mathbf{e}_x$)에 초점을 맞춘다. 시뮬레이션은 높은 레이놀즈 수($Re = 10,000$)와 페클레 수($Pe = 10,000$), 그리고 낮은 프루드 수($Fr = 0.15$및$Fr = 0.1$)에서 수행되었으며, 그 결과 부력 레이놀즈 수($Re_b = Fr^2 Re$)는 각각 225와 100이 되었다.

초기 조건이 경로 선택에 미치는 영향을 격리하기 위해, 두 가지 DNS 세트가 수행되었다:

1. DNS1 (경로 P1): 격자 규모의 백색 잡음(white noise)으로 초기화하여, 일차 수평 전단 불안정성의 모든 불안정한 모드가 성장할 수 있도록 하였다.
2. DNS2 (경로 P2): 격자 규모의 잡음과 더불어 일차 불안정성의 가장 빠르게 성장하는 2D($k_z=0$) 고유 모드에 대한 특정 "시드(seed)"를 추가하여 초기화하였다.

추가로, 저자들은 DNS2에서 나타나는 기둥형 흐름의 "동결된(frozen-in)" 해석적 모델에 대한 선형 안정성 분석을 수행하였다. 이 모델은 비선형적으로 진화한 2D 흐름을 테일러-그린 와류 배열과 사인파 흐름이 중첩된 형태로 근사하며, 이를 통해 이차 3D 섭동의 고유 모드와 성장률을 계산할 수 있다.

주요 기여 및 결과

1. 단일 설정 내에서의 두 경로 확인:
   본 연구는 초기 섭동 스펙트럼에 의해서만 구별되는, 동일한 모델 설정 내에서 두 가지 경로가 모두 발생할 수 있음을 입증하였다:

• 경로 P1 (직접 3D 성장): 백색 잡음으로 초기화될 경우, 일차 수평 전단 불안정성의 $k_z \neq 0$ 고유 모드들이 직접 성장한다. 이 모드들은 광범위한 스케일에서 수직 전단을 생성하며, $Re_b$가 충분히 클 때 이후 이차적인 소규모 KH 불안정성에 의해 불안정해진다. 이 경로는 체적 강제된 콜모고로프 흐름에서의 관측 결과와 일치한다.
• 경로 P2 (기둥형 진화): $k_z=0$ 모드가 시딩될 경우, 이 모드가 흐름을 지배하여 배경 흐름을 장기간 지속되는 시간 의존적 2D 기둥형(와류형) 흐름으로 몰아넣는다. 이 2D 상태는 이후 이차 3D 모드에 의해 불안정해진다. 저자들은 이 이차 모드들을 와류 배열의 쌍곡점(hyperbolic stagnation points)에서 발생하는 **쌍곡형 불안정성(hyperbolic instabilities)**으로 식별하였다. 이 불안정성은 높은 수직 파수($k_z \sim Fr^{-1}$)를 가지며 수직 전단의 출현을 유도하고, 이는 결국 삼차 KH 불안정성을 촉발한다.

2. 경로 P2의 이차 불안정성 특성 규명:
   동결된 기둥형 흐름에 대한 선형 안정성 분석을 통해, 저자들은 이차 불안정성이 고-$k_z$ 쌍곡형 불안정성임을 보여주었다. 이 불안정성의 성장률은 쌍곡점의 국소 변형률($\Delta_h$)에 비례하며, 큰 진폭에서는 성층화 강도와 거의 무관하게 $\lambda \approx 0.46 \Delta_h$로 나타난다. 이 불안정성은 고립된 와류 쌍에서 관찰되는 전통적인 "지그재그(zigzag)" 불안정성과 구별되는, 쌍곡점 주변의 수직 속도장에 집중된 "주름(ripples)" 형태로 나타난다.

3. 혼합 효율에 미치는 영향:
   연구는 두 경로에 대한 이성분 혼합 및 혼합 효율($\eta$)을 정량화하였다.

• 부력 소산 ($\chi$): 최대 부력 변동 소산율은 두 경로에서 유사하게 나타나는데, 이는 두 경우 모두 소규모 운동($k_z \gtrsim Fr^{-1}$)에 의해 제어되기 때문이다.
• 점성 소산 ($\epsilon$): 경로 P2는 경로 P1보다 현저히 높은 점성 소산을 보인다. 이는 P2의 쌍곡형 불안정성이 좁고 높은-$k_z$ 스펙트럼을 자극하여 얇은 수직 전단층에서 강렬한 소산을 일으키기 때문이다. 반면, P1은 넓은 수직 스케일 스펙트럼을 자극하여 소산을 더 넓게 분산시킨다.
• 혼합 효율: 결과적으로, 최대 혼합 효율은 경로 P2($\eta_{max} \approx 0.165$)보다 경로 P1($\eta_{max} \approx 0.25$)에서 유의미하게 높다. 이러한 차이는 P2가 혼합되는 부력량에 비해 더 많은 운동 에너지를 소산시키기 때문에 발생한다.

의의 및 주장
본 논문은 굴곡점(inflection point)이 존재하는 충분히 큰 $Re_b$ 조건의 강한 성층화된 흐름에서, 수직 전단이 만드는 소규모 KH 불안정성의 출현은 수평 전단 불안정성의 필연적인 부산물이라고 주장한다. 두 경로 사이의 선택은 초기 조건, 즉 $k_z=0$ 모드의 초기 진폭과 $k_z \neq 0$ 모드 간의 비율에 의해 결정된다.

저자들은 기존 문헌(특히 하이퍼볼릭 탄젠트 전단층 연구)에서 경로 P2가 빈번하게 나타나는 현상이, 계산을 가속화하기 위해 $k_z=0$ 시드를 흔히 추가하여 직접적인 3D 모드의 성장을 억제하는 초기화 관행의 결과일 수 있다고 시사한다. 반대로, 일부 맥락에서 P1이 나타나지 않는 것은 성층화가 약해 $k_z \neq 0$ 모드의 불안정한 범위가 제한적이기 때문일 수 있다.

궁극적으로, 이 연구는 수직 전단이 미리 존재하지 않더라도 수평 전단 불안정성이 성층화된 난류와 이성분 혼합을 생성하는 핵심 요소임을 확립한다. 이는 난류의 경로(직접 3D 성장 vs 2D-to-3D 전이)가 에너지 소산 스펙트럼과 그에 따른 혼합 효율을 근본적으로 변화시키며, 이는 기후 및 지구물리학 모델의 난류 매개변수화 시 반드시 고려되어야 할 요소임을 강조한다. 저자들은 본 결과가 이상적인 흐름에서는 견고하지만, 실제 지구물리학적 환경에서의 회전 및 주변 파동장과의 상호작용은 여전히 미해결 과제로 남아 있다고 언급하였다.