안정한 성층 대기 흐름을 위한 에너지와 플럭스 예산 난류 폐쇄 모델과 내부 중력파의 역할

초록

본 논문은 기존의 에너지·플럭스·예산(EFB) 난류 폐쇄 모델을 확장하여 대규모 내부 중력파(IGW)가 난류에 미치는 영향을 정량화한다. IGW는 난류 운동에너지(TKE), 난류 퍼텐셜 에너지(TPE) 및 온도 플럭스에 추가 생산·소산 항을 제공한다. 결과적으로 무파 상태에서 일정한 0.2‑0.25였던 플럭스 리차드슨 수가 파에너지에 따라 크게 변동하고, 균일 성층에서는 1을 초과하며, 비균일 성층에서는 음수가 될 수 있다(역온도 플럭스). IGW는 난류 이방성을 감소시키고 TPE 비중을 높이며, 파의 전파 방향에 따라 총 수직 모멘텀 플럭스를 강화하거나 약화시킨다. 모델 예측은 관측·실험·DNS·LES 자료와 일치한다.

상세 분석

이 연구는 Zilitinkevich 등(2007, 2008)이 제시한 에너지·플럭스·예산(EFB) 폐쇄 모델을 기반으로, 대규모 내부 중력파(IGW)가 안정적인 성층 대기 난류에 미치는 복합적인 영향을 체계적으로 도입하였다. 기존 EFB 모델은 난류 운동에너지(TKE)와 난류 퍼텐셜 에너지(TPE)의 상호 전환, 그리고 평균 흐름에 대한 모멘텀·열 플럭스의 예측에 초점을 맞추었으나, 파동에 의한 추가적인 에너지·플럭스 전달 메커니즘을 고려하지 않았다. 본 논문은 IGW가 제공하는 수직 모멘텀 플럭스와 온도 플럭스의 추가 생산항을 각각 τ_w와 F_θ,w 로 정의하고, 이들을 TKE·TPE 방정식에 삽입함으로써 파-난류 상호작용을 정량화한다.

주요 결과는 다음과 같다. 첫째, 플럭스 리차드슨 수(R_f)라는 무차원 지표가 파에너지 수준에 따라 크게 변동한다. 균일한 성층에서는 IGW 에너지가 증가할수록 R_f는 0.2‑0.25에서 1을 초과하는 값까지 상승한다. 이는 파가 제공하는 추가적인 퍼텐셜 에너지 생산이 열 플럭스를 강화시켜, 전통적인 ‘임계’ R_f 값을 초월함을 의미한다. 반면, 성층이 수직으로 비균일(즉, 위쪽으로 강해지는 경우)일 때는 파가 강한 성층 쪽으로 전파하면서 온도 플럭스가 역방향(반경사)으로 전환될 수 있다. 이때 R_f는 0을 지나 음수까지 내려가며, 이는 난류가 평균 온도 구배와 반대 방향으로 열을 수송한다는 비직관적인 현상을 설명한다.

둘째, IGW는 난류의 이방성을 완화한다. 파가 제공하는 수직 진동 성분이 기존의 수평 주도 난류 구조에 섞이면서, 수직 속도 성분의 분산이 증가하고, 결과적으로 난류 길이 스케일의 수직·수평 비율이 1에 가까워진다. 이는 기존 EFB 모델이 예측한 ‘수평 얇은 난류’와는 대조적인 특성이다.

셋째, TPE가 전체 난류 에너지(TTE = TKE + TPE)에서 차지하는 비중이 상승한다. 파에 의해 유도된 퍼텐셜 에너지 생산이 TKE보다 상대적으로 크게 기여하기 때문에, 특히 강한 안정성(큰 브루넬 수) 구역에서 TPE 비율이 30 % 이상으로 증가한다. 이는 대기 경계층 상부와 같은 고도에서 파동‑난류 상호작용이 에너지 저장 메커니즘을 재편함을 시사한다.

넷째, 파의 전파 방향에 따라 총 수직 모멘텀 플럭스(τ_total = τ_turb + τ_w)의 부호와 크기가 달라진다. 파가 하향 전파(하강 파)일 경우 τ_w가 τ_turb와 같은 부호를 가져 전체 모멘텀 전송을 강화한다. 반대로 상승 파는 τ_w가 τ_turb와 반대 부호가 되어 전송을 약화시킨다. 이러한 효과는 관측된 고도별 풍속 프로파일의 비대칭성을 설명하는 데 유용하다.

마지막으로, 제안된 모델은 다양한 실험·관측 자료와 비교했을 때 높은 일치성을 보인다. 대기 현장 측정, 회전식 수조 실험, 직접수치시뮬레이션(DNS) 및 대규모 난류 시뮬레이션(LES)에서 얻은 R_f, 난류 이방성, TPE 비중, 모멘텀 플럭스 변화를 모두 재현한다. 따라서 IGW를 포함한 확장 EFB 모델은 안정적인 성층 대기 흐름을 이해하고 예측하는 데 강력한 도구가 된다.