[논문 리뷰] Quasi-geostrophic model of the instabilities of the Stewartson layer

Nathanaël Schaeffer, Philippe Cardin|arXiv (Cornell University)|2003. 12. 22.

Geomagnetism and Paleomagnetism Studies참고 문헌 24인용 수 6

한 줄 요약

이 논문은 질량 수율을 기여하는 비수렴 흐름을 고려하여, 회전하는 축대칭 용기 내에서 스타우어트슨 shear 계층의 불안정성 문제를 연구하기 위해 준지구학적 2차원 모델을 개발한다. 임계 로시오 수의 점근적 스케일링 법칙을 유도하여, 평평한 용기의 경우 킬빈-헬름홀츠 유사 메커니즘에 기반한 $E^{3/4}$ 의존성과 기울은 용기의 경우 로시오 웨이브 불안정성에 기반한 $\beta E^{1/2}$ 의존성을 보여주며, 저익먼 수에서의 수치 시뮬레이션과 실험 결과로 검증된다.

ABSTRACT

We study the destabilization of a shear layer, produced by differential rotation of a rotating axisymmetric container. For small forcing, this produces a shear layer, which has been studied by Stewartson and is almost invariant along the rotation axis. When the forcing increases, instabilities develop. To study the asymptotic regime (very low Ekman number $E$), we develop a quasi-geostrophic two-dimensional model, whose main original feature is to handle the mass conservation correctly, resulting in a divergent two-dimensional flow, and valid for any container provided that the top and bottom have finite slopes. We use it to derive scalings and asymptotic laws by a simple linear theory, extending the previous analyses to large slopes (as in a sphere), for which we find different scaling laws. For a flat container, the critical Rossby number for the onset of instability evolves as $E^{3/4}$ and may be understood as a Kelvin-Helmoltz shear instability. For a sloping container, the instability is a Rossby wave with a critical Rossby number proportional to $βE^{1/2}$, where $β$ is related to the slope. We also investigate the asymmetry between positive and negative differential rotation and propose corrections for finite Ekman and Rossby numbers. Implemented in a numerical code, our model allows us to study the onset over a broad range of parameters, determining the threshold but also other features such as the spatial structure. We also present a few experimental results, validating our model and showing its limits.

연구 동기 및 목표

• 유한한 기울기 경계를 가진 고속으로 회전하는 용기 내에서 스타우어트슨 shear 계층의 선형 불안정성을 분석하기 위해.
• 기존의 점근적 스케일링 법칙을 평평한 기하구조를 초월하여 구형 및 기울인 용기로 확장하기 위해.
• 질량 수율을 정확히 처리하기 위해 비수렴 2차원 흐름을 포함하는 준지구학 모델을 개발하기 위해.
• 다양한 익먼 수와 용기 기하구조에서 불안정성 발생의 임계 로시오 수를 결정하기 위해.
• 특히 지구물리학적 및 행성 유동에 관련된 저익먼 수 영역에서 수치 시뮬레이션과 실험 결과와의 일치를 검증하기 위해.

제안 방법

• 회전하는 shear 계층에서 질량 수율을 정확히 구현하기 위해 비수렴 흐름을 포함하는 준지구학적 2차원 모델을 수립한다.
• 평평한 및 기울인 용기 기하구조에 대해 선형 안정성 이론을 사용하여 점근적 스케일링 법칙을 도출한다.
• 익먼 마찰과 코리olis 효과를 포함하는 수정된 압력 방정식을 도입하며, 용기 기울기를 나타내는 $\beta$ 매개변수를 포함한다.
• 스펙트럼 방법을 사용하여 적응형 반경 방향 격자 해상도를 적용하고, 슬립 경계 조건 하에서 스트림함수 $\psi$ 의 고유값 문제를 해결한다.
• 익먼 수 범위가 넓은 영역($10^{-5}$ 에서 $10^{-10}$)에서 수치 시뮬레이션을 수행하여 불안정성 임계 조건을 규명한다.
• 수중 탱크 실험에서의 실험 데이터와 모델 예측을 비교한다. 이는 구형 및 디스크형 기하구조에서 수행된다.

실험 결과

연구 질문

• RQ1 스타우어트슨 계층에서 임계 로시오 수의 스케일링 법칙은 용기 기하구조에 따라 어떻게 달라지며, 특히 기울인 경계와 평평한 경계의 경우에 어떻게 다를까?
• RQ2 준지구학 모델에서 비수렴 흐름의 역할은 무엇이며, 비수렴 근사보다 불안정성 임계 조건을 더 잘 포괄하기 위해 어떻게 개선되었는가?
• RQ3 왜 평평한 용기에서는 임계 로시오 수가 $E^{3/4}$ 에 비례하고 기울인 용기에서는 $\beta E^{1/2}$ 에 비례하는가? 이 차이를 이끌어내는 물리적 메커니즘은 무엇인가?
• RQ4 유한한 익먼 수와 로시오 수는 불안정성 임계 조건에서 양의 및 음의 미세한 회전 차이에 어떤 영향을 미치며, 비대칭성을 어떻게 변화시키는가?
• RQ5 매우 낮은 익먼 수($E < 10^{-5}$) 영역에서 준지구학 모델이 행성 내부 유동에 관련된 실험 및 수치 결과를 어느 정도 정확하게 재현할 수 있는가?

주요 결과

• 평평한 용기의 경우, 임계 로시오 수는 $\mbox{{Ro}}_c \sim E^{3/4}$ 로 스케일링되며, 이는 킬빈-헬름홀츠 shear 불안정성 메커니즘과 일치한다.
• 기울인 용기의 경우, 임계 로시오 수는 $\mbox{{Ro}}_c \sim \beta E^{1/2}$ 로 스케일링되며, 이는 용기의 기울기로 인해 유도되는 로시오 웨이브 불안정성에 기인한다.
• 모델은 익먼 수의 넓은 범위($10^{-5}$ 에서 $10^{-10}$)에서 불안정성 임계 조건을 정확히 예측하며, 수치 시뮬레이션과 실험 결과 모두와 일치한다.
• 구형 케이지에서의 실험 데이터는 모델과 양호한 일치를 보이며, 익먼 수 $E = 10^{-5}$ 에서 임계 로시오 수는 $4.0 \times 10^{-2}$ 에서부터 $E = 10^{-10}$ 에서는 $5.4 \times 10^{-5}$ 로 변화한다.
• 모델은 양의 및 음의 미세한 회전 간에 뚜렷한 비대칭성을 드러내며, 음의 회전에서 더 강한 불안정성이 발생함을 보여주며, 이는 유한한 익먼 수와 로시오 수 효과를 포함함으로써 보정된다.
• 경계 조건 민감도 분석 결과, 슬립 조건과 익먼 펌프 효과가 $\mbox{{Ro}}_c$ 에 강한 영향을 미치며, 수평적 비수렴과 부피 점성도는 작지만 무시할 수 없는 영향을 미친다.