[논문 리뷰] Hodograph Method and Numerical Solution of the Two Hyperbolic Quasilinear Equations. Part III. Two-Beam Reduction of the Dense Soliton Gas Equations

E. V. Shiryaeva, M. Yu. Zhukov|arXiv (Cornell University)|2015. 12. 21.

Nonlinear Waves and Solitons인용 수 4

한 줄 요약

이 논문은 보존 법칙에 기반한 호도그래프 방법을 사용하여 조밀한 솔리톤 가스 방정식의 두 빔 감소(보른-인펠트 방정식과 동치)에 대해 분석적이고 수치적인 해를 제시한다. 히퍼볼릭 비선형 편미분방정식을 선형 편미분방정식에 대한 명시적 리만-그린 함수를 통해 상미분방정식의 초기값 문제로 변환함으로써, 인위적 점성 없이 정확하고 근사하지 않은 수치적 통합이 가능해지며, 주기적 및 파동 패킷 초기 조건에 대해 높은 정밀도로 해를 도출할 수 있다.

ABSTRACT

The paper presents the solutions for the two-beam reduction of the dense soliton gas equations (or Born-Infeld equation) obtained by analytical and numerical methods. The method proposed by the authors is used. This method allows to reduce the Cauchy problem for two hyperbolic quasilinear PDEs to the Cauchy problem for ODEs. In some respect, this method is analogous to the method of characteristics for two hyperbolic equations. The method is effectively applicable in all cases when the explicit expression for the Riemann-Green function for some linear second order PDE, resulting from the use of the hodograph method for the original equations, is known. The numerical results for the two-beam reduction of the dense soliton gas equations, and the shallow water equations (omitting in the previous papers) are presented. For computing we use the different initial data (periodic, wave packet).

연구 동기 및 목표

• 두 히퍼볼릭 비선형 편미분방정식의 초기값 문제를 해결하기 위한 수치적 방법을 개발한다. 이 방정식들은 두 빔 감소된 조밀한 솔리톤 가스 방정식에서 유도된다.
• 보존 법칙에 기반한 호도그래프 방법을 적용하여 원래의 편미분방정식을 근사하지 않고도 상미분방정식의 초기값 문제로 감소시킨다.
• 기존에 생략된 얕은 수면 방정식을 포함한 주기적 및 파동 패킷 초기 조건에 대해 보른-인펠트 방정식에 대해 이 방법의 효과성을 입증한다.
• 다중값 및 불연속 해를 다룰 수 있는 프레임워크를 제공하며, 인위적 점성 또는 수치적 확산을 도입하지 않는다.

제안 방법

• 호도그래프 방법은 보존 법칙과 관련 선형 편미분방정식에 대한 명시적 리만-그린 함수를 활용하여 원래의 히퍼볼릭 비선형 편미분방정식을 상미분방정식 시스템으로 변환한다.
• 해는 리만 불변량 $ R^1, R^2 $ 을 통해 암시적으로 표현되며, 시간 $ t(a,b) $ 와 위치 $ x(a,b) $ 는 보존 법칙 밀도 $ \varphi $ 의 적분을 통해 정의된다.
• 이 방법은 라그랑주 변수 $ a, b $ 를 사용하며, $ R^1(x,t) = R^1_0(b) $, $ R^2(x,t) = R^2_0(a) $ 를 만족하고, 매 시간 등면 $ t = t_* $ 에서 매개변수 $ \tau $ 에 대한 상미분방정식을 통해 해를 재구성한다.
• 상미분방정식 시스템 $ a(\tau), b(\tau), X(\tau) $ 의 초기값 문제를 $ t_a, t_b $ 의 편미분을 포함한 상미분방정식으로 해결하며, 이들은 리만-그린 함수와 초기 조건에서 계산된다.
• 계산 효율성을 향상시키기 위해 $ \theta $ 를 사용한 매개변수 간소화 기법이 도입되며, 이로 인해 상미분방정식이 $ \frac{da}{d\theta} = R^2_0(a) - R^1_0(a) $ 등의 형태로 단순화된다.
• 유한차분, 유한체적, 리만 솔버 근사와 같은 방법을 회피하여, 해의 정확도가 오직 상미분방정식 통합 정밀도에만 의존하도록 보장한다.

실험 결과

연구 질문

• RQ1 보존 법칙에 기반한 호도그래프 방법이 조밀한 솔리톤 가스 방정식의 두 빔 감소에 효과적으로 적용될 수 있는가?
• RQ2 두 히퍼볼릭 비선형 편미분방정식의 초기값 문제를 수치적 근사 없이 상미분방정식 시스템으로 감소시킬 수 있는가?
• RQ3 보른-인펠트 방정식에 대해 주기적 및 파동 패킷 초기 조건에 적용했을 때 이 방법의 수치적 성질은 어떠한가?
• RQ4 이 방법은 인위적 점성 또는 수치적 확산을 도입하지 않고도 해의 구조를 유지하는가?
• RQ5 이 방법은 히퍼볼릭 시스템에서 발생하는 다중값 및 불연속 해를 다룰 수 있는가?

주요 결과

• 호도그래프 방법은 조밀한 솔리톤 가스 방정식의 두 빔 감소를 상미분방정식의 초기값 문제로 성공적으로 감소시켜 정확한 수치 통합을 가능하게 한다.
• 원래의 편미분방정식을 근사하지 않고 상미분방정식만 수치적으로 풀기 때문에, 인위적 점성 없이 고정밀도 해를 도출할 수 있다.
• 주기적 및 파동 패킷 초기 조건에 대한 수치 결과는 이 방법의 강건성과 복잡한 파동 구조에 대한 적용 가능성에 대한 증거를 제공한다.
• 등면에서의 해는 매개변수 $ \tau $ 에 대한 상미분방정식을 통해 계산되며, $ t_a $ 와 $ t_b $ 는 리만-그린 함수와 초기 조건에서 유도된다.
• 매개변수 변환 기법을 통해 $ \theta $ 를 도입함으로써 상미분방정식이 단순화되어 계산 효율성이 향상되면서도 정밀도는 그대로 유지된다.
• 이 방법은 리만-그린 함수가 명시적으로 존재하는 모든 방정식에 대해 광범위하게 적용 가능하며, 얕은 수면 방정식 및 보른-인펠트 방정식을 포함한다.