[논문 리뷰] Numerical investigations of traveling singular sources problems via moving mesh method

Zhicheng Hu, Keiwei Liang|arXiv (Cornell University)|2012. 03. 08.

Fluid Dynamics and Vibration Analysis참고 문헌 22인용 수 3

한 줄 요약

이 논문은 여러 개의 이동하는 특이 소스가 서로 다른 속도로 이동하는 1차원 PDE를 해결하기 위해 도메인 분할과 결합한 이동 메쉬 방법을 제안한다. 각 소스 주변의 도메인을 하위 영역으로 나누고, 각 하위 영역에서 독립적으로 국소 메쉬 방정정식을 풀음으로써 해의 도함수 점프를 계산하지 않아도 되어, 이차 수준의 공간 수렴성과 효율적인 병렬 계산이 가능해진다. 이 방법은 다양한 속도로 이동하는 소스에 대해 폭발 현상을 정확히 포착할 수 있으며, 이는 이전에 수치적으로 연구되지 않은 케이스이다.

ABSTRACT

This paper studies the numerical solution of traveling singular sources problems. In such problems, a big challenge is the sources move with different speeds, which are described by some ordinary differential equations. A predictor-corrector algorithm is presented to simulate the position of singular sources. Then a moving mesh method in conjunction with domain decomposition is derived for the underlying PDE. According to the positions of the sources, the whole domain is splitted into several subdomains, where moving mesh equations are solved respectively. On the resulting mesh, the computation of jump $[\dot{u}]$ is avoided and the discretization of the underlying PDE is reduced into only two cases. In addition, the new method has a desired second-order of the spatial convergence. Numerical examples are presented to illustrate the convergence rates and the efficiency of the method. Blow-up phenomenon is also investigated for various motions of the sources.

연구 동기 및 목표

• 시간에 따라 변화하는 특이 소스가 서로 다른 속도로 이동하는 PDE를 수치적으로 해결하는 데 도전하는 것.
• 해의 도함수 점프 [u_x]로 인해 표준 유한차분/유한요소 방법이 이동 소스 위치를 통과할 때 실패하는 문제를 해결하는 것.
• 소스 위치에서 [u_x] 점프 항을 명시적으로 계산하지 않는 안정적인, 이차 수렴성을 확보한 수치적 스킴을 개발하는 것.
• 도메인 분할과 하위 영역에서의 국소 메쉬 균일 분포를 통해 효율적인 병렬 계산을 가능하게 하는 것.
• 비균일 및 주기적 속도를 가진 소스의 운동 패턴에 따라 해의 폭발 행동을 연구하는 것.

제안 방법

• q개의 특이 소스의 시간에 따라 변화하는 위치 α_i(t)를 제어하는 ODE dα_i/dt = ψ_i(t, α_i(t), u)에 따라 예측-수정 알고리즘을 사용하여 추적한다.
• 계산 도메인은 q개의 소스로 분리된 q+1개의 하위 영역으로 나뉘며, 각 하위 영역에서 독립적으로 국소 이동 메쉬 방정식을 푼다.
• 메쉬는 균일 분포 원칙에 기반한 MMPDEs(MMPDE4/5/6 등)를 통해 생성되며, 소스 주변에 메쉬 점이 집중되도록 한다.
• 물리적 PDE는 유도된 적응 메쉬 위에서 유한차분 스킴을 사용하여 이산화하며, 도메인 분할 덕분에 두 종류의 다른 스킴만 필요로 한다.
• 각 하위 영역을 독립적으로 다룸으로써 [u_x] 점프를 직접 계산하지 않아도 되어 이산화 과정이 단순화된다.
• 특수한 시간 스텝 구현을 통해 전체 스킴의 이차 수준 공간 수렴성을 보장한다.

실험 결과

연구 질문

• RQ1 도메인 분할을 통한 이동 메쉬 방법이 다수의 특이 소스가 비균일한 속도로 이동하는 PDE를 효과적으로 다룰 수 있는가?
• RQ2 이동 소스 위치에서 해의 도함수 점프 [u_x]를 수치 이산화 과정에서 어떻게 피할 수 있는가?
• RQ3 제안된 방법의 공간 수렴 속도는 얼마이며, 이차 정확도를 유지할 수 있는가?
• RQ4 소스가 서로 다른 속도로 이동할 경우 폭발 현상을 방법이 정확히 포착할 수 있는가?
• RQ5 이 방법은 다양한 운동 패턴을 가진 다수의 소스(q > 2)로 확장 가능한가?

주요 결과

• 수치 실험을 통해 특정 시간 스텝 전략 하에서 제안된 방법이 이차 수준 공간 수렴성을 확보함을 확인하였다.
• 속도 k=2로 이동하는 단일 소스의 경우, t ≈ 2.0397에 위치 x ≈ 4.0794에서 폭발이 발생하며, u_max ≈ 3.16×10^6에 도달한다.
• 동일한 속도 k=2로 이동하는 두 개의 소스의 경우, 두 번째 소스에서 t ≈ 1.6896에 폭발이 발생하며, u_max ≈ 3.16×10^6에 도달한다.
• A=π인 대칭적 주기적 운동의 경우, 두 소스에서 동시에 t ≈ 2.4969에 폭발이 발생하며, u_max ≈ 3.16×10^6에 도달한다.
• 다른 속도로 이동하는 소스에 대해 폭발 현상을 성공적으로 시뮬레이션하였으며, 이는 이전에 수치적으로 연구되지 않은 케이스이다.
• 국소 흡수 경계 조건을 사용함으로써 관측 도메인을 작게 유지해도 해의 정확도에 영향을 주지 않으며, 수치 결과로 이를 입증하였다.