An entropy-stable oscillation-eliminating dgsem for the euler equations on curvilinear meshes

1. 연구 배경 및 문제 제기 (Problem)

• 오일러 방정식의 복잡성: 압축성 오일러 방정식은 비선형 쌍곡형 보존 법칙으로, 충격파 (shock waves), 접촉 불연속면, 희박파 등 다양한 구조를 가집니다.
• 수치적 난제:
  • 고차 정확도 vs. 진동: 고차 정확도 (high-order accuracy) 를 가진 방법 (예: DG) 은 매끄러운 영역에서는 정확하지만, 강한 불연속면 근처에서는 깁스 현상 (Gibbs phenomenon) 으로 인해 비물리적 진동이 발생합니다.
  • 엔트로피 안정성: 물리적으로 타당한 해를 얻기 위해서는 이산 수준에서도 엔트로피 부등식 (entropy inequality) 을 만족해야 합니다.
  • 곡선 격자의 어려움: 복잡한 형상을 모델링하기 위해 곡선 격자를 사용할 경우, 기하학적 항 (metric terms) 의 이산화가 엔트로피 안정성과 자유류 보존 (free-stream preservation) 을 해칠 수 있어 추가적인 주의가 필요합니다.
  • 기존 OEDG 의 한계: 기존 진동 제거 DG (OEDG) 방법은 직교 격자 (Cartesian meshes) 에 국한되었으며, 곡선 격자로의 확장이 어렵고, 모든 셀에 적용될 경우 계산 비용이 매우 높았습니다.

2. 제안된 방법론 (Methodology)

이 논문은 엔트로피 안정 DGSEM과 수정된 OEDG를 결합한 새로운 프레임워크를 제시합니다.

A. 엔트로피 안정 DGSEM (Entropy-Stable DGSEM)

• SBP (Summation-By-Parts) 성질: 가우스 - 르장드 - 로바토 (GLL) 점을 기반으로 한 노드형 DGSEM 을 사용하여 SBP 성질을 만족시킵니다.
• 곡선 격자 처리: 곡선 좌표 변환에 따른 메트릭 항 (metric terms) 을 이산적으로 처리할 때, **이산 메트릭 항등식 (discrete metric identities)**을 만족하도록 설계하여 엔트로피 안정성을 보장합니다.
• 엔트로피 보존 체적 플럭스: 체적 항 (volume terms) 에서 점별 플럭스 대신 **엔트로피 보존 2 점 플럭스 (entropy-conservative two-point flux, Chandrashekar flux)**를 사용하여 이산 엔트로피 보존 법칙을 유도합니다.
• 인터페이스 플럭스: 요소 간 결합에는 엔트로피 안정성을 가진 수치 플럭스 (예: 로컬 랙 - 프리드리히스 플럭스) 를 사용하여 전역 엔트로피 부등식을 만족시킵니다.

B. 수정된 진동 제거 DG (Modified OEDG)

• 사후 처리 (Post-processing): 시간 적분 단계 후, 진동이 발생하는 영역에서 DG 해에 대해 OEDG 연산자를 적용합니다.
• 투영 연산자 기반 곡선 격자 확장: 기존 OEDG 가 직교 격자에서의 국소 직교 모드 기저에 의존했던 것과 달리, 본 논문은 **투영 연산자 (projection operators)**를 사용하여 일반 곡선 격자에서도 체계적으로 적용 가능하도록 재구성했습니다.
• 충격 지시자 (Shock Indicator) 의 발견:
  • OEDG 과정에 사용되는 **0 차 감쇠 계수 (zero-order damping coefficient)**가 자연스럽게 효과적인 충격 지시자 역할을 한다는 것을 발견했습니다.
  • 이를 통해 **고통받는 셀 (troubled cells)**만 선별적으로 OEDG 를 적용할 수 있게 되어 계산 비용을 크게 절감했습니다.
• 보존성 유지: OEDG 연산자는 0 차 모드 (셀 평균) 를 변경하지 않으므로, 질량, 운동량, 에너지의 보존성을 유지합니다.

3. 주요 기여 (Key Contributions)

1. 곡선 격자를 위한 엔트로피 안정 OEDG 프레임워크 개발: 복잡한 기하학적 형상에서도 엔트로피 안정성과 고차 정확도를 동시에 보장하는 최초의 DGSEM 기반 OEDG 방법 중 하나입니다.
2. 계산 효율성 극대화: OEDG 감쇠 계수를 충격 지시자로 활용하여, 진동 제어 메커니즘을 불연속면이 있는 영역에만 국소화함으로써 계산 비용을 획기적으로 줄였습니다.
3. 일반 격자 확장: 투영 연산자를 도입하여 OEDG 절차를 직교 격자에서 곡선 격자로 확장하는 체계적인 방법을 제시했습니다.
4. 강건한 수치 실험: 다양한 난제 (충격파, 접촉 불연속, 복잡한 와류 등) 에 대한 검증으로 방법론의 유효성을 입증했습니다.

4. 수치 실험 결과 (Results)

논문에서는 다음과 같은 테스트 케이스를 통해 제안된 방법의 성능을 검증했습니다.

• 등엔트로피 오일러 와류 (Isentropic Euler Vortex):
  • 곡선 격자에서 정확한 해를 가진 문제를 수행하여 $(N+1)$차 수렴률을 확인했습니다.
  • 전체 엔트로피가 시간에 따라 감소하거나 일정하게 유지됨을 확인하여 엔트로피 안정성을 입증했습니다.
• 2 차원 리만 문제 (Riemann Problems):
  • 복잡한 파동 상호작용을 가진 리만 문제 (Case 12, 13) 를 해결하여 충격파와 접촉 불연속면을 정밀하게 포착함을 보였습니다.
  • 충격 지시자가 진동이 발생하는 영역에만 정확히 활성화됨을 확인했습니다.
• 더블 매치 반사 (Double Mach Reflection):
  • 강한 충격파와 켈빈 - 헬름홀츠 (Kelvin-Helmholtz) 불안정성 등 미세한 유동 구조를 잘 해상도했습니다.
  • **충격 지시자 임계값 (C)**을 조절하여 계산 비용을 최적화할 수 있음을 보였습니다 (OEDG 적용 시 시간 단계 제한이 완화되어 전체 계산 시간이 단축됨).
• 원통 주변의 초음속 유동 (Supersonic Flow Around Cylinder):
  • 복잡한 곡선 격자에서 bow shock, fishtail shock, 와류 구조 등을 성공적으로 시뮬레이션했습니다.
  • 감쇠 계수 (scaling factor $s$) 조절을 통해 수치 확산과 미세 구조 해상도 사이의 균형을 조절할 수 있음을 보였습니다.
• 다중 원통 내 폭발 (Explosion in Multiple Cylinders):
  • 복잡한 장애물이 있는 곡선 격자에서 충격파의 반사, 회절, 와류 생성 등 복잡한 유동 현상을 강건하게 처리했습니다.

5. 의의 및 결론 (Significance)

이 연구는 고차 정확도, 엔트로피 안정성, 진동 제어라는 세 가지 핵심 요구사항을 일반 곡선 격자 환경에서 동시에 충족하는 수치 기법을 제시했습니다.

• 실용성: 복잡한 기하학적 형상을 가진 실제 공학 문제 (항공기, 터빈 등) 에 적용 가능한 강력한 도구로 평가됩니다.
• 효율성: 불필요한 영역에서의 진동 제거 연산을 배제함으로써 고차 방법의 계산 부담을 크게 줄였습니다.
• 확장성: 제안된 프레임워크는 나비에 - 스토크스 방정식이나 자기유체역학 (MHD) 등 더 복잡한 시스템으로 확장될 수 있는 잠재력을 가지고 있습니다.

결론적으로, 이 논문은 복잡한 기하학적 구조에서의 압축성 유동 해석을 위한 강건하고 정확한 고차 수치 해법의 새로운 기준을 제시했다는 점에서 중요한 의의를 가집니다.