[논문 리뷰] Commentary on `The Three-Dimensional Current and Surface Wave Equations' by George Mellor

Fabrice Ardhuin, Alastair D. Jenkins외 1명|ArXiv.org|2005. 04. 14.

Ocean Waves and Remote Sensing참고 문헌 64인용 수 14

한 줄 요약

이 논문은 파장이 얕은 곳으로 올라오는 파면에서 메일러(2003)의 3차원 수류 및 표면파 방정식에 기인한 기본적인 모순을 규명한다. 이는 삼각좌표계에서 파동에 의해 유도되는 운동량 유속의 부적절한 근사화 때문이며, 이로 인해 수치적 보정을 제안하지만, 다른 좌표계—예를 들어 일반화된 라그랑주 평균(GLM) 또는 준유동속도 기반 좌표계—는 파동과 평균류의 운동량을 철저히 분리함으로써 문제를 완전히 피할 수 있다. 이는 수치적 파동 힘 평가 없이도 일관되고 해석적으로 다룰 수 있는 방정식을 가능하게 한다.

ABSTRACT

The lowest order sigma-transformed momentum equation given by Mellor (J. Phys. Oceangr. 2003) takes into account a phase-averaged wave forcing based on Airy wave theory. This equation is shown to be generally inconsistent due to inadequate approximations of the wave motion. Indeed the evaluation of the vertical flux of momentum requires an estimation of the pressure $p$ and coordinate transformation function $s$ to first order in parameters that define the large scale evolution of the wave field, such as the bottom slope. Unfortunately there is no analytical expression for $p$ and $s$ at that order. A numerical correction method is thus proposed and verified. Alternative coordinate transforms that allow a separation of wave and mean flow momenta do not suffer from this inconsistency nor require a numerical estimation of the wave forcing. Indeed, the problematic vertical flux is part of the wave momentum flux, thus distinct from the mean flow momentum flux, and not directly relevant to the mean flow evolution.

연구 동기 및 목표

• 메일러(2003)의 삼각좌표 운동량 방정식이 얕은 곳에서 파동이 올라오는 조건에 적용되었을 때 발생하는 모순을 규명하고 진단하는 것.
• 메일러의 수식에서 파동에 기인한 운동량 유속 항이 파동 압력과 좌표 변환에 대해 O(ε₂) 수준의 추정이 필요하며, 이는 해석적 해가 존재하지 않는다는 것을 보여주는 것.
• 모순의 근본 원인이 원래 수식에서 파동 운동량 유속이 평균류 운동량 유속와 혼동되어 있기 때문임을 보여주는 것.
• 메일러의 방정식에 대한 파동 힘 항의 수치 보정 방법을 제안하는 것.
• 특히 GLM 기반 또는 준유동속도 방정식과 같은 대체 수식을 주장함으로써 파동과 평균류 운동량을 본질적으로 분리함으로써 수치적 파동 힘 평가가 필요 없이도 일관된 방정식을 얻을 수 있음을 설명하는 것.

제안 방법

• 작은 파고 기울기(ε₁ ≪ 1)와 느린 해저 기울기(ε₂ ≪ 1) 조건 하에서 메일러(2003) 수식의 최저차수 삼각좌표 변환 운동량 방정식을 해석적으로 분석하는 것.
• 수직 운동량 유속 p∂s/∂x가 파동 압력 p와 좌표 함수 s에 의존하며, 이 둘 다 O(ε₂) 수준에서 해석적 표현이 없음을 규명하는 것.
• 메일러의 방정식에서 일관성을 회복하기 위해 O(ε₂) 수준에서 파동 힘 항을 추정하는 데 수치적 보정 방법을 제안하는 것.
• 앤드류스와 마킨타이어(1978)의 일반화된 라그랑주 평균(GLM) 방정식과 메일러의 수식을 비교하여 GLM가 파동과 평균류 운동량 유속을 별도로 처리함으로써 모순을 피할 수 있음을 보여주는 것.
• 파동 운동량을 평균류 진화에서 분리하기 위해 ū = U − u_S 라는 준유동속도 개념을 도입하여, 명시적 파동 힘 항 없이도 일관된 방정식을 가능하게 하는 것.
• 라플라스 방정식의 수치적 해를 사용하여 이론적 논거를 검증하고, 마그네 등(2006)에서 사용한 다섯 모드 근사화조차도 정확한 파동 힘 평가에 부족함을 보여주는 것.

실험 결과

연구 질문

• RQ1 에너지 소산이 없는 조건에서 메일러(2003)의 삼각좌표 운동량 방정식이 얕은 곳에서 파동이 올라오는 경우에 왜 일관성이 없나요?
• RQ2 메일러 수식에서 파동에 기인한 운동량 유속 항이 왜 O(ε₂) 수준에서 수치적 평가가 필요하며, 이는 실용적인 해양 모델링에 있어 왜 문제인가요?
• RQ3 메일러의 방정식에서 발생하는 모순은 수치적 보정을 통해 해결될 수 있으며, 그 방법의 효과는 어느 정도인가요?
• RQ4 GLM 기반 또는 준유동속도 수식은 왜 메일러의 방정식에서 발생하는 파동 힘 모순을 피할 수 있나요?
• RQ5 표면를 따라가는 좌표계인 메일러의 방식은 어떤 조건에서 물질 표면을 올바르게 표현하지 못하며, 그 결과 운동량 유속 평가에 어떤 영향을 미치나요?

주요 결과

• 메일러(2003)의 최저차수 삼각좌표 변환 운동량 방정식은 파동에 기인한 압력 p와 좌표 함수 s가 O(ε₂) 수준에서 해석적 표현이 없기 때문에 일관성이 없다.
• Eulerian 운동량 유속 항 uw를 대체하는 파동 운동량 유속 항 S_x3 = p̃^ξ ∂s̃/∂x 는 파동 운동량을 포함하고 있으며, 평균류 진화와 직접적인 관련이 없어 개념적이고 계산적 혼란을 야기한다.
• O(ε₂) 수준에서 파동 힘 항을 추정하기 위한 수치적 보정 방법을 제안하지만, 이는 일반적으로 운영용 해양 모델에 비현실적으로 높은 계산 비용을 수반한다.
• 일반화된 라그랑주 평균(GLM) 기반 또는 준유동속도 ū = U − u_S 기반의 대체 수식은 파동과 평균류 운동량을 철저히 분리함으로써 모순을 피하고, 명시적 파동 힘 평가가 필요 없게 된다.
• GLM 기반 준유동속도 방정식은 해저 기울기와 유속 기울기 존재하에서도 최저차수에서 일관되고 해석적으로 다룰 수 있으며, 파동 운동량 유속 항이 평균류 운동량 방정식에 나타나지 않기 때문이다.
• 수치 실험 결과는 마그네 등(2006)에서 사용한 다섯 모드 근사화조차도 정확한 파동 힘 평가에 부족함을 보이며, 더 높은 해상도 또는 대체 수식이 필요함을 시사한다.