[논문 리뷰] Nonlinear Geometric Optics method based multi-scale numerical schemes for highly-oscillatory transport equations

Nicolas Crouseilles, Shi Jin외 1명|arXiv (Cornell University)|2016. 05. 31.

Advanced Numerical Methods in Computational Mathematics참고 문헌 23인용 수 3

한 줄 요약

이 논문은 비선형 기하광학 앵색을 통해 진동하는 위상 요소를 독립 변수로 통합함으로써 고주파 진동 운반 방정식에 대한 새로운 점근적 보존 수치적 방법을 제안한다. 이 방법은 고주파 진동을 해상하지 않아도 파장 주기 동안 균일한 일계 정밀도를 확보하며, 소수 파arameter ε에 의존하지 않는 메쉬와 시간 스텝을 사용하여 안정적이고 정확한 시뮬레이션을 가능하게 한다. 다양한 척도에서의 수치 실험을 통해 검증되었다.

ABSTRACT

We introduce a new numerical strategy to solve a class of oscillatory transport PDE models which is able to captureaccurately the solutions without numerically resolving the high frequency oscillations {\em in both space and time}.Such PDE models arise in semiclassical modeling of quantum dynamics with band-crossings, and otherhighly oscillatory waves. Our first main idea is to use the nonlinear geometric optics ansatz, which builds theoscillatory phase into an independent variable. We then choose suitable initial data, based on the Chapman-Enskog expansion, for the new model. For a scalar model, we prove that so constructed model will have certain smoothness, and consequently, for a first order approximation scheme we prove uniform error estimates independent of the (possibly small) wave length. The method is extended to systems arising from a semiclassical model for surface hopping, a non-adiabatic quantum dynamic phenomenon. Numerous numerical examples demonstrate that the method has the desired properties.

연구 동기 및 목표

• 소파장 ε를 해상하지 않아도 고주파 운반 방정식의 해를 정확히 포착할 수 있는 수치적 방법을 개발하기 위해.
• 고주파 진동에서 Δx, Δt = O(ε) 가 필요로 하는 기존 방법의 계산 불가능성을 해결하기 위해.
• 기하광학 원리를 통합함으로써 점근적 보존(AP) 방법을 공간적·시간적 진동 문제로 확장하기 위해.
• Chapman-Enskog 전개를 기반으로 한 정밀하게 구성된 초기 조건을 통해 ε에 독립적인 균일한 오차 추정을 확보하기 위해.
• 비단조화 양자 전이를 포함하는 반도체 표면 점프 모델에서 발생하는 시스템에 적용하기 위해.

제안 방법

• 진동하는 위상을 독립 변수 τ로 변환하기 위해 비선형 기하광학(NGO) 앵색을 사용함으로써 고주파 진동을 물리적 공간과 시간에서 분리한다.
• 위상에 의존하는 계수를 갖는 수정된 운반 방정식을 따르는 진폭 함수 V에 대한 새로운 방정식 시스템을 유도한다.
• 원래 문제의 진동적 구조와 일치하는 스무스함과 일致성을 확보하기 위해 Chapman-Enskog 유형의 전개를 이용해 초기 자료를 구성한다.
• 공간 이산화에는 시간 분할 스펙트럼 방법을 사용하고, 운반 및 완화 항은 정확한 시간 적분을 수행하며, τ에서의 빠른 진동은 의사 스펙트럼 방법으로 처리한다.
• 시스템의 경우, 표면 점프를 모델링하는 2×2 시스템으로 확장되며, (x, p, τ) 변수로의 변환과 함께 τ에서의 정확한 적분이 수행된다.
• ε를 해상하지 않기 위해 굵은 메쉬(Δx ≫ ε)와 고정된 시간 스텝을 사용하지만, 위상 변수 변환과 적절한 초기 자료를 통해 정확도를 유지한다.

실험 결과

연구 질문

• RQ1 소파장 ε를 해상하지 않아도 운반 방정식의 고주파 해를 포착할 수 있는 수치적 방법을 설계할 수 있는가?
• RQ2 비선형 기하광학 앵색이 ε → 0 근처에서 수치 해의 균일한 수렴성을 가능하게 하는가?
• RQ3 Chapman-Enskog 전개를 사용해 변환된 시스템에서 스무스함과 정확도를 유지하는 초기 자료를 구성할 수 있는가?
• RQ4 비단조화 전이가 포함된 시스템, 예를 들어 표면 점프 모델에서 이 방법은 어떻게 작동하는가?
• RQ5 표준 방법이 계산적으로 금기되는 수준인 ε = 1/256일 때에도 정확도와 효율성을 유지할 수 있는가?

주요 결과

• 스칼라 모델에서 이 방법은 ε에 독립적인 파장 주기 동안 균일한 일계 오차 추정을 달성하며, Δx나 Δt가 ε에 비례해야 한다는 요구 조건 없이 수렴성을 입증한다.
• ε = 1일 때, 새로운 방법은 N_x = 32, N_p = 64, N_τ = 8를 사용하여 참값(Δt = 0.05, N_x = 256)과 동일한 정확도를 달성했으며, CPU 시간은 15초에서 1분으로 감소했다.
• ε = 1/32일 때, 새로운 방법은 동일한 굵은 메쉬(N_x = 32)를 유지하면서 정확도를 유지는 반면, 참값 방법은 Δt = 0.02와 N_x = 512 가 필요했고, CPU 시간은 75초로 증가했다.
• ε = 1/256일 때, 새로운 방법은 동일한 파rameter(Δt = 0.05, N_x = 32)로 정확한 결과를 생성했지만, 참값 방법은 Δt = 0.0005와 N_x = 4096 가 필요했고, CPU 시간은 1420초로 증가했다.
• 수치 결과는 메쉬가 공간 진동보다 더 굵어도 점별 해와 밀도를 정확히 포착함을 보여주었다.
• 이 방법은 비단조화 전이를 포함한 복잡한 시스템, 예를 들어 2×2 표면 점프 모델에서도 효과적으로 작동하여 다양한 ε 척도에서 강건성과 확장성을 입증했다.