Structure-Preserving High-Order Methods for the Compressible Euler Equations in Potential Temperature Formulation for Atmospheric Flows

기술 요약: 전위온도 정식화 기반 압축성 오일러 방정식을 위한 구조 보존 고차 방법론

문제 정의
본 논문은 대기 흐름에 특화된 압축성 오일러 방정식의 구조 보존 수치 방법을 개발하는 문제를 다룬다. 불연속 갈레르킨(Discontinuous Galerkin, DG) 방법은 고차 정확도와 병렬 확장성을 갖추고 있어 널리 사용되지만, 장기 시뮬레이션에서 안정성을 유지하기 위해 추가적인 안정화 기법(예: 필터, 인공 점성)을 요구하는 경우가 많으며, 이는 정확도를 저하시킬 수 있다. 대기 모델링의 특정 과제는 예측 변수의 선택이다. 총 에너지 정식화는 일반적인 전산유체역학(CFD)에서 표준적으로 사용되지만, 전위온도($\theta$) 정식화는 기상학적 이론 역할과 열역학 방정식을 단순한 수송 방정식으로 축소할 수 있는 능력 때문에 대기 모델에서 선호된다. 그러나 일반 곡선 좌표 격자 및 중력 소스 항이 존재하는 환경에서 전위온도 정식화를 사용하여 구조 보존 스킴(엔트로피, 총 에너지, 운동 에너지 보존)을 구축하는 것은 여전히 복 복잡한 작업이다. 또한, 기존의 잘 균형 잡힌(well-balanced) 스킴들은 곡선 격자 상에서 등온 상태와 일정한 전위온도 배경 상태 모두에 대해 정수압 평형을 유지하는 데 어려움을 겪는 경우가 많다.

방법론
저자들은 플럭스 차분(flux-differencing) 접근 방식을 사용하는 DGSEM(Discontinuous Galerkin Spectral Element Method) 기반의 프레임워크를 개발한다. 핵심 방법론은 다음과 같다:

1. 정식화: 연구는 두 가지 방정식 세트에 집중한다: 전통적인 총 에너지 정식화와 전위온도 정식화이다. 후자의 경우, 예측 변수는 밀도($\varrho$), 운동량($\varrho \mathbf{V}$), 그리고 전위온도($\varrho \theta$)이다.
2. 엔트로피 및 에너지 보존: 저자들은 특정 보존 특성을 만족하도록 설계된 새로운 이점(two-point) 수치 플럭스를 유도한다:
   • 엔트로피 보존 (EC): 열역학적 엔트로피의 보존을 보장한다.
   • 총 에너지 보존 (TEC): 총 에너지(내부 에너지 + 운동 에너지)의 보존을 보장한다.
   • 운동 에너지 보존 (KEP): 압력 구배가 없는 상태에서 운동 에너지를 보존한다.
   • 압력 평형 보존 (PEP): 일정한 압력 및 속도 상태를 유지한다.
     저자들은 전위온도 정식화의 경우, EC, TEC, KEP, PEP를 동시에 달성하기 위해서는 밀도 플럭스가 압력에 의존해야 함을 엄밀한 점프 조건 및 엔트로피 변수 분석을 통해 확인하였으며, 이는 필수적인 트레이드오프임을 입증하였다.
3. 소스 항 이산화: 중력 퍼텐셜 항($\phi$)을 처리하기 위해, 저자들은 소스 항을 비보존적 곱(non-conservative product)으로 취급한다. 저자들은 플럭스 차분 접근 방식을 채택하여, 스킴이 잘 균형 잡히도록 특정 평균값(등온 상태의 경우 로그 평균, 일정한 전위온도 상태의 경우 스톨라르스키 평균)을 사용하여 소스 항을 이산화한다.
4. 곡선 확장: 이 방법들을 3차원 곡선 격자로 확장한다. 저자들은 유도된 플럭스와 소스 항 이산화가 이산 메트릭 항등식(discrete metric identities)을 만족하는 한, 임의의 좌표 사상 하에서도 구조 보존 특성(TEC, KPEP, well-balancedness)을 유지함을 증명한다.
5. 구현: 본 방법론은 Julia 언어와 Trixi.jl 프레임워크를 사용하여 구현되었으며, SSPRK43 시간 적분법을 채택하였다.

주요 기여

• 전위온도를 위한 새로운 플럭스: 저자들은 전위온도 정식화에 특화된 새로운 수치 플럭스(TEC, EC, 그리고 결합된 ETEC 플럭스)를 유도하고 그 특성을 규명하였다. ETEC 플럭스가 EC, TEC, KEP, PEP를 만족하는 유일한 플럭스임을 증명하였으나, 이는 압력 의존적인 밀도 플럭스를 필요로 한다.
• 중력을 포함한 일반화된 KPEP 및 TEC: 저자들은 중력 퍼텐셜 항이 존재하는 상황에서 운동 및 위치 에너지 보존(KPEP) 특성을 일반화하였으며, 총 에너지 정식화와 전위온도 정식화 모두에 대한 총 에너지 보존 조건을 확립하였다.
• 일반화된 잘 균형 잡힌(Well-Balanced) 스킴: 곡선 격자 상에서 등온($T=$ const) 및 일정한 전위온도($\theta=$ const) 배경 상태 모두에 대해 정수압 평형을 보존할 수 있는 잘 균형 잡힌 스킴을 개발한 것이 중요한 기여이다. 기존 연구들은 종종 이러한 능력을 데카르트 격자나 특정 배경 상태로 제한하였다.
• 비보존적 소스 항 처리: 저자들은 DGSEM 프레임워크 내에서 중력 소스 항에 대한 비보존적 곱 이산화를 엄밀하게 적용하였으며, 이를 통해 원하는 구조적 특성을 보존하는 플럭스 차분 형태의 반-이산(semi-discrete) 방정식 재작성 방법을 제시하였다.

결과
저자들은 일련의 수치 실험을 통해 이론적 발견을 검증하였다:

• 보존 특성: 중력이 없는 상태에서, 새로운 플럭스는 매끄러운 흐름(밀도파) 및 난류 흐름(Taylor-Green vortex)에 대해 기계 정밀도(machine precision)까지 엔트로피와 총 에너지를 정확하게 보존함을 보여준다.
• 잘 균형 잡힌 특성 (Well-Balancedness): 뒤틀린 곡선 격자 위에서, 스킴은 등온 및 단열(일정한 $\theta$) 배경 상태 모두에 대해 정수압 평형을 유지한다. 오차는 500,000 타임스텝의 장기 적분 동안 기계 정밀도를 유지하며, 이는 비보존적 소스 항 이산화의 견고함을 확인시켜 준다.
• 수렴성: 관성 중력파에 대한 수렴성 분석 결과, 소스 항의 비보존적 이산화는 점별(point-wise) 이산화에 비해(특히 곡선 격자에서) 수직 속도 오차를 최대 4자리(orders of magnitude)까지 낮추는 것을 보여주었다. 스킴은 기대되는 고차 수렴율을 달로 달성한다.
• 대기 시나리오: 방법론을 고전적인 산악파 벤치마크(선형 정수압, 선형 비정수압, Schär 산악파)에 테스트하였다. 결과는 비보존적 정식화가 점별 소스 항 이산화에 비해 더 견고하며, 불필요한 진동이 적음을 보여준다. 전위온도 정식화는 이러한 시나리오에서 총 에너지 정식화와 대등한 성능을 보인다.
• 바록클리닉 불안정성 (Baroclinic Instability): 저자들은 큐브드 스피어(cubed-sphere) 격자 상에서 바록클리닉 불안정성 벤치마크(DCMIP 2016) 테스트를 시작함으로써, 복잡한 전 지구적 대기 흐름에 대한 방법론의 적용 가능성을 추가로 검증하였다.

의의
본 논문은 이러한 발전이 대기 모델링 커뮤니티에 큰 관심 대상이라고 주장한다. 전위온도 정식화에 특화된 구조 보존 방법을 제공함으로써, 본 연구는 총 에너지 변수로 전환하지 않고도 기존 대기 모델에 고차, 안정적, 정확한 수치 스킴을 통합할 수 있도록 돕는다. 복잡한 지형과 장기 적분이 표준인 실제 대기 시뮬레이션에서, 일반 곡선 격자 상에서 여러 배경 상태에 대해 정수압 평형을 보존할 수 있는 능력은 핵심적인 진전으로 강조된다. 이 연구는 전위온도 정식화의 이론적 이점과 구조 보존 수치 이산화의 실질적 요구 사항 사이의 간극을 메우는 역할을 한다.