Stochastic compressible Navier-Stokes equations under location uncertainty and their approximations for ocean modelling

본 논문은 위치 불확실성 하의 확률론적 압축성 나비에-스톡스 방정식에 대한 이론적 및 수치적 연구를 제시하며, 이는 잠재 에너지에서 상당한 압축 효과를 드러내고 아격자 수직 혼합 모델의 에너지 일관성을 향상시키는 Boussinesq 근사를 통해 해양 모델링에 적용됨을 보여줍니다.

핵심

바다를 거대하고 격렬하게 끓어오르는 국물 한 그릇으로 상상해 보세요. 오랫동안 과학자들은 이 국물의 움직임을 완벽한 결정론적 규칙을 사용하여 그 '레시피'를 작성하려고 노력해 왔습니다. 그들은 모든 물 분자의 정확한 위치와 속도를 안다면 국물이 다음에 어디로 이동할지 정확히 예측할 수 있다고 가정했습니다.

그러나 바다는 그보다 너무 복잡합니다. 직접 보거나 측정하기에는 너무 작은 무수한 혼란스러운 소용돌이 (난류) 가 가득 차 있습니다. 모든 소용돌이를 하나씩 시뮬레이션하려는 시도는 해변의 모래 알갱이 하나하나를 세어 조수를 예측하려는 것과 같습니다. 불가능할 뿐만 아니라 계산 비용도 너무 많이 듭니다.

이 논문은 새로운 레시피 작성 방식을 제안합니다. 모든 모래 알갱이를 추적하려 하기보다, 저자들은 불확실성에 기반한 '수정 계수 (fudge factor)'를 추가할 것을 제안합니다. 이를 위치 불확실성 (Location Uncertainty, LU) 프레임워크라고 부릅니다.

다음은 핵심 아이디어를 간단한 개념으로 분해한 것입니다:

1. '취한 사람의 보행' 비유

혼잡한 시장을 걷고 있다고 상상해 보세요. 당신은 명확한 목적지 ('대규모' 흐름) 를 가지고 있지만, 당신을 부딪히는 사람들은 당신을 무작위 방향으로 약간 빗나가게 밀어냅니다.

• 옛 방식: 당신을 부딪히는 모든 사람의 정확한 경로를 계산하려 합니다.
• 새 방식 (LU): 당신이 밀려날 것이라고 받아들입니다. 당신의 움직임을 매끄러운 보행 더하기 무작위적이고 떨리는 '브라운 운동' (취한 사람의 보행과 유사) 으로 모델링합니다. 당신은 밀림이 당신을 정확히 어디로 밀어낼지 알지 못하지만, 밀림의 통계적 특성 (얼마나 강한지, 그리고 어떻게 상관관계가 있는지) 은 알고 있습니다.

2. '압축 가능한' 국물

대부분의 해양 모델은 물이 '비압축성'이라고 가정합니다. 즉, 찌그러뜨릴 수 없는 고체 젤리 덩어리와 같다는 것입니다. 하지만 실제로는 물이 약간 찌그러질 수 있으며, 특히 압력 변화나 온도 변화 시에는 더욱 그렇습니다.

• 저자들은 압축 가능한 물 (찌그러질 수 있는 물) 의 완전하고 복잡한 물리학으로 시작합니다.
• 그런 다음 이 복잡한 시스템에 그들의 '무작위 밀림' 수학을 적용합니다.
• 결과: 그들은 기존 방정식과 유사하지만 추가 항이 포함된 새로운 방정식 세트를 유도해 냅니다. 이 추가 항들은 무작위 밀림이 행하는 '일'을 나타냅니다. 군중이 당신을 밀 때 전달되는 에너지라고 생각하세요. 그것은 단순한 무작위 밀기가 아니라, 실제로 당신의 속도와 몸의 온도를 변화시킵니다.

3. 혼합물 속의 '숨겨진 열'

이 논문은 온도와 대류 (따뜻한 물은 상승하고 차가운 물은 가라앉음) 에 중점을 둡니다.

• 문제: 표준 모델에서 차가운 물이 가라앉을 때, 종종 혼합층 (바다의 상부 부분) 의 바닥에서 갑자기 멈춥니다. 실제로는 이러한 차가운 물의 '플룸 (plume)'이 창처럼 더 깊고 따뜻한 물 속으로 관통하는 경우가 많습니다. 이를 **관통 대류 (penetrative convection)**라고 합니다.
• 발견: 저자들이 새로운 확률론적 모델을 실행했을 때, '무작위 밀림' 항들이 자연스럽게 이 관통 효과를 재현한다는 것을 발견했습니다.
• 비유: 군중 (바다) 이 무거운 상자 (차가운 물 플룸) 를 이동시키려 한다고 상상해 보세요. 표준 모델은 상자를 멈추게 하는 단단한 벽처럼 작용합니다. 반면 새로운 모델은 혼란스러운 군중처럼 작용합니다. 무작위적인 밀치기는 상자가 벽을 넘어 예상보다 더 깊게 들어갈 수 있도록 충분한 추가 운동량을 제공합니다.

4. 에너지를 측정하는 두 가지 방법

저자들은 시스템의 에너지를 측정하는 방식에 대해 흥미로운 점을 발견했습니다:

• 내부 에너지 (열기): 열만 볼 때, '찌그러짐' (압축) 효과는 미미하여 중요하지 않았습니다. 이는 기존의 더 단순한 모델과 일치했습니다.
• 위치 에너지 (높이): 하지만 중력장 내에서 물이 얼마나 높은지와 관련된 에너지를 볼 때, '찌그러짐' 효과는 매우 중요해졌습니다.
• 교훈: 공이 튀는 것을 측정하는 것과 같습니다. 공이 바닥에 닿을 때 얼마나 뜨거워지는지만 측정한다면, 튀는 것은 중요하지 않은 것처럼 보입니다. 하지만 공이 얼마나 높이 튀는지 측정한다면, 그 충격은 엄청납니다. 저자들은 그들의 모델에 있는 무작위 압력 항이 숨겨진 스프링처럼 작용하여 에너지 예산에서 물이 얼마나 높이 '튀는지'에 영향을 미친다는 것을 발견했습니다.

5. '드리프트 (Drift)'와 '확산 (Diffusion)'

수학은 '수정 계수'로 작용하는 두 가지 구체적인 새로운 항을 생성합니다:

• 드리프트 (이토 - 스토크스 드리프트): 무작위 밀림이 완벽하게 균일하지 않기 때문에 발생하는 체계적인 밀림입니다. 이는 바위 (난류) 가 특정 패턴으로 배열되어 있기 때문에 약간 다르게 흐르는 강류와 같습니다.
• 확산: 무작위 밀림에 의해 발생하는 퍼져 나가는 효과입니다.

성과 요약

저자들은 바다의 messy하고 혼란스러운 현실과 이를 예측하는 데 사용하는 깔끔한 수학적 모델 사이의 다리를 성공적으로 구축했습니다.

• 그들은 가능한 가장 복잡한 물리학 (압축성, 무작위성, 열역학) 으로 시작했습니다.
• 그것을 '표준' 해양 관점 (부시네스크 근사) 으로 단순화했을 때, 새로운 방정식이 여전히 작동하며 실제로 차가운 물이 얼마나 깊게 가라앉는지 예측을 개선한다는 것을 보였습니다.
• 그들은 정답을 얻기 위해 모든 작은 소용돌이를 시뮬레이션할 필요가 없으며, 물이 어디로 가려는지에 대한 불확실성을 수학적으로 고려하기만 하면 된다는 것을 증명했습니다.

간단히 말해, 그들은 '모든 모래 알갱이를 세는' 불가능한 임무를 '모래가 흩어지는 경향을 고려하는' 더 현명한 전략으로 대체했으며, 이 접근법이 이전 방법들보다 바다의 깊고 관통하는 해류를 훨씬 잘 포착한다는 것을 발견했습니다.

기술 요약

기술 요약: 해양 모델링을 위한 확률적 압축성 나비어-스토크스 방정식及其 근사

문제 제기
해양 모델링은 기후 시뮬레이션 (예: 대서양 경순환) 에 중요한 하위 격자 규모 과정, 특히 수직 혼합과 침투성 대류를 표현하는 데 상당한 도전에 직면해 있습니다. 전통적인 접근법은 종종 비압축성을 가정하고 압축성 효과를 무시하는 Boussinesq 근사에 의존합니다. 많은 맥락에서 타당하지만, 이 근사는 온도와 운동량 플럭스 간의 에너지 일관성을 흐리게 하며, 압축 및 팽창과 관련된 음파 전파나 특정 열역학적 일 항을 포착하지 못합니다. 더 나아가, 해양 역학을 위한 기존 확률적 모델은 대부분 비압축성 유동에 대해 유도되어, 압축성 열역학적으로 활성인 유체를 위한 이론적 프레임워크에 공백이 존재합니다. 저자들은 위치 불확실성 (LU) 을 고려한 압축성 유동에 대한 엄밀한 확률적 형식을 개발하여, 하위 격자 규모 수직 모델의 에너지 일관성을 개선하기 위해 필요한 열역학적 항을 명시적으로 유도하고자 합니다.

방법론
본 논문은 유체 입자의 변위를 시간 미분 가능한 해결된 속도와 브라운 운동으로 모델링된 고도로 요동치는 시간 비상관 성분으로 분해하는 위치 불확실성 (LU) 프레임워크를 활용합니다.

1. 이론적 유도:

   • 저자들은 질량, 종 질량, 운동량, 총 에너지의 기본 보존 법칙에서 출발합니다.
   • 확률적 소스 항을 포함하도록 확장된 확률적 레이놀즈 수송 정리 (SRTT) 를 적용합니다. 이는 보존형과 비보존형 (기본 변수형) 모두에서 일반적인 확률적 압축성 나비어-스토크스 방정식 체계로 이어집니다.
   • 결정적으로, 이 유도는 확률적 변위에 작용하는 힘 (압력과 점성 응력) 의 일과 열 플럭스를 명시적으로 고려합니다. 그 결과, 노이즈 유발 확산, Itô-Stokes 표류, 확률적 미분 곱에서 비롯된 2 차 공변항을 포함한 특정 확률적 항이 도출됩니다.
   • 이후 이 체계는 Boussinesq 및 정수압 근사를 거칩니다. 저자들은 항의 크기 순서를 신중하게 추적하여 열역학적 효과, 특히 위치 에너지를 통해 표현되는 압축의 일을 유지하는 확률적 Boussinesq 체계를 유도합니다.

2. 수치 적용:

   • 이론적 프레임워크를 검증하기 위해 저자들은 CROCO 모델 (비정수압, 비 Boussinesq 기능) 을 사용하여 성층화된 해양의 온도 구동 자유 대류 사건에 대한 와류 대류 시뮬레이션 (LES) 을 수행합니다.
   • 해결된 규모와 하위 격자 규모 요동을 분리하기 위해 LES 데이터에 3 시간 비상관 시간 척도를 사용하여 시간 필터링 기법을 적용합니다.
   • 하위 격자 규모 잔차는 확률적 노이즈 (브라운 운동) 와 그 분산 텐서를 구성하는 데 사용됩니다.
   • 이후 진단된 통계량에 의해 강제된 오프라인 1 차원 확률적 모델을 통해 수평 평균 온도 프로파일의 진화를 시뮬레이션합니다. 이를 통해 확률적 모델의 출력과 기준 LES 결과 간의 직접적인 비교가 가능해집니다.

주요 기여

• 일반적 확률적 압축성 형식: 본 논문은 종 수송과 전체 에너지 예산을 포함하는 확률적 압축성 나비어-스토크스 방정식의 완전한 유도를 제공합니다. 이는 주로 비압축성 유동에 국한되었던 이전 LU 작업을 확장한 것입니다.
• 열역학적 일관성: 이 유도는 압축 및 팽창의 일을 명시적으로 포함합니다. 저자들은 이러한 항이 내부 에너지로 표현될 때 (Boussinesq 근사와 일관되게) 무시할 수 있지만, 위치 에너지로 표현될 때 중요해짐을 보여줍니다. 이는 혼합층 내에서 운동 에너지와 위치 에너지를 결합하는 새로운 확률적 압력 항을 드러냅니다.
• 열역학을 포함한 확률적 Boussinesq 모델: 저자들은 이러한 열역학적 효과를 통합한 확률적 Boussinesq 모델을 유도하여, 이전의 확률적 등적 모델보다 물리적으로 일관된 프레임워크를 제공합니다.
• 침투성 대류 진단: 본 연구는 확률적 항, 특히 Itô-Stokes 표류와 확률적 확산이 혼합층 하부에서의 수직 온도 플럭스를 재현할 수 있음을 확인합니다. 이는 표준 매개변수화가 종종 침투성 대류 효과 (혼합층 하부의 온난화) 를 포착하지 못하는 영역입니다.

결과

• 에너지 예산: 수치 실험에서 온도 변화의 주된 균형은 Boussinesq 근사와 일관되게 내부 에너지 수송에 의해 지배됩니다. 그러나 위치 에너지를 분석할 때, 확률적 압력과 관련된 2 차 공변항 ($A_u$) 과 표류 일 ($D_\sigma$) 같은 확률적 항은 상당한 기여를 보입니다.
• 혼합층 역학: 확률적 모델은 침투성 대류와 관련된 혼합층 하부의 온난화를 성공적으로 재현합니다. 이는 명시적인 확산 폐쇄 가정에 의존하지 않고 수평 변동성과 관성 과정에 의해 구동되는 비국소적 수직 역학을 포착할 수 있는 모델의 능력에 기인합니다.
• 한계: 혼합층 자체 내에서 확률적 모델은 기준 LES 에 비해 제한된 기술을 보입니다. 저자들은 이 영역의 난류장의 상대적으로 낮은 변동성이 현재 진단 설정에서 확률적 항의 기여를 무시할 수 있게 만든다고 지적합니다. 또한 모델은 표면 근처에서 어려움을 겪는데, 이는 아마도 확률적 표현이 LES 에 의해 해결된 특정 소규모 요동을 포착하지 못하기 때문일 것입니다.

의의 및 주장
저자들은 이 작업이 확률적 설정 내에서 해양 역학에 대한 새로운 관점을 여는 구조화된 모델링 프레임워크를 제공한다고 주장합니다.

• 물리적 정당성: 현상학적 접근법과 달리, 확률적 방정식의 추가 항은 보존 법칙과 확률 미적분학에서 엄밀하게 유도되어 그 형태와 역할에 대한 명확한 물리적 정당성을 제공합니다 (예: 압축성 LES 의 바로피크날 일과 유사).
• 에너지 일관성: 이 프레임워크는 역학과 열역학 간의 결합에 관한 물리적 불확실성과 근사를 명시적으로 고려함으로써, 해양 일반 순환 모델 (OGCM) 에서 사용되는 하위 격자 규모 수직 모델의 에너지 일관성을 개선할 수 있는 경로를 제시합니다.
• 침투성 대류: 확률적 항이 혼합층 하부에서 누락된 레이놀즈 응력을 정성적으로 재현할 수 있다는 점은, 완전한 비정수압 솔버의 계산 비용 없이 침투성 대류 효과를 표현할 수 있는 가능성을 시사합니다.
• 향후 방향: 본 논문은 현재 진단적 접근법이 유망하지만, 향후 작업은 노이즈 특성 (예: TKE 기반 모델을 통한) 의 예보적 모델링과 혼합층 및 표면 근처의 역학을 더 잘 포착하기 위한 무작위 표면 압력 강제력의 포함에 초점을 맞춰야 한다고 겸손하게 제안합니다.