[논문 리뷰] Analytic solutions for the Burgers equation with source terms

Gino I. Montecinos|arXiv (Cornell University)|2015. 03. 29.

Nonlinear Waves and Solitons참고 문헌 3인용 수 11

한 줄 요약

이 논문은 강성과 비선형 형태를 포함한 소스 항이 있는 버거스 방정식을 특성 기반 분석 방법으로 해결하는 방법을 제시한다. 편미분방정식(PDE)을 두 개의 해석 가능한 상미분방정식(ODE)으로 변환함으로써—해결 가능한 등가 ODE(특성 곡선 沿해 진화를 지배하는)와 특성 ODE(공간-시간에서 곡선을 정의하는)—소스 항의 역수에 원시 함수가 존재할 경우 정확한 해를 도출할 수 있다. 주요 기여는 수치적 방법의 검증을 위해 비선형성과 강성을 반영한 정확한 해를 생성하는 체계적인 절차를 제공하는 것이다. 특히 특성 곡선을 구하기 위해 이분법을 사용한 수치적 근사가 필요한 경우에도 적용 가능하다.

ABSTRACT

Analytic solutions for Burgers equations with source terms, possibly stiff, represent an important element to assess numerical schemes. Here we present a procedure, based on the characteristic technique to obtain analytic solutions for these equations with smooth initial conditions.

연구 동기 및 목표

• 강성과 비선형 형태를 포함한 소스 항이 있는 버거스 방정식의 정확한 해석적 해를 유도하기 위한 체계적인 방법을 개발하는 것.
• 비선형성과 강성을 반영한 정확한 해를 제공함으로써 수치적 방법의 엄밀한 검증을 가능하게 하는 것.
• 특성 기반 ODE 시스템이 해석 가능한 조건을 설정하는 것—특히 소스 항과 그 역수에 원시 함수가 존재하는지에 중점을 두고.
• 특성 곡선의 은둔적 의존성 문제를 해결하기 위해 초기 조건에 따라 발생하는 대수방정식을 해결하기 위한 이분법 도입

제안 방법

• 이 방법은 편미분방정식(PDE)을 특성 곡선을 따라 간소화하여 두 개의 상미분방정식(ODE) 시스템으로 변환한다.
• 특성 곡선을 따라 PDE를 유도한 등가 ODE는 $ \frac{d}{dt}\hat{q}(t) = s(\hat{q}(t)) $로 표현되며, 초기 조건은 $ \hat{q}(0) = h_0(y) $이다.
• 등가 ODE의 해석 가능성은 $ s(q)^{-1} $가 원시 함수 $ G(q) $를 갖는 경우 보장되며, 이 경우 해는 $ \mathcal{E}(t, h_0(y)) = G^{-1}(t + G(h_0(y))) $로 표현된다.
• 특성 ODE는 $ \frac{d}{dt}x(t) = \mathcal{E}(t, h_0(y)) $로 표현되며, $ \mathcal{E} $를 적분함으로써 $ x = y + \mathcal{F}(t, h_0(y)) $의 형태로 해석 가능해지며, 여기서 $ \mathcal{F} $는 $ \mathcal{E} $의 원시 함수이다.
• 전체 해는 $ q(x,t) = \mathcal{E}(t, h_0(y)) $로 암시적으로 정의되며, $ y $ 는 $ x = y + \mathcal{F}(t, h_0(y)) $를 만족한다. 이 경우 해석적 표현이 어려운 비표현적 형태에서는 수치적 근사가 필요하다.
• 비표현적 형태의 $ y $ 에 대해서는 특성 곡선 정의에서 유도된 대수방정식을 해결하기 위해 이분법을 적용한다.

실험 결과

연구 질문

• RQ1 강성 또는 비선형 소스 항이 있는 버거스 방정식에 대해 특성 기반 접근을 통해 해석적 해를 도출할 수 있는가?
• RQ2 등가 ODE가 해석 가능한 조건은 소스 항에 대해 어떤가?
• RQ3 특성 곡선 방정식이 은둔적일 경우 해는 초기 조건에 따라 어떻게 표현할 수 있는가?
• RQ4 특성 곡선 수식에서 발생하는 은둔적 대수방정식을 해결하는 데 효과적인 수치 방법은 무엇인가?

주요 결과

• 등가 ODE는 소스 항의 역수 $ s(q)^{-1} $가 원시 함수 $ G(q) $를 갖는 경우 해석 가능하며, 이 경우 해는 $ \mathcal{E}(t, h_0(y)) = G^{-1}(t + G(h_0(y))) $로 표현된다.
• 특성 ODE는 등가 ODE의 해 $ \mathcal{E}(t, h_0(y)) $가 원시 함수 $ \mathcal{F}(t, h_0(y)) $를 갖는 경우 해석 가능하며, 이 경우 곡선은 $ x = y + \mathcal{F}(t, h_0(y)) $로 표현된다.
• 선형 소스 항 $ s(q) = \beta q $ 의 경우 해는 $ \hat{q}(t) = \frac{h_0(y)}{1 - \beta t h_0(y)} $이며, 특성 곡선은 $ x = y - \frac{\ln(1 - \beta t h_0(y))}{\beta} $로 표현된다.
• 이차 소스 항 $ s(q) = \beta q^2 $ 의 경우 해는 $ \hat{q}(t) = \frac{h_0(y)}{1 - \beta t h_0(y)} $로 동일하며, 특성 곡선은 명시적으로 해석 가능하다.
• 비선형 소스 항 $ s(q) = e^{\beta q} $ 의 경우 해는 $ \hat{q}(t) = -\frac{1}{\beta} \ln(e^{-\beta h_0(y)} - \beta t) $이며, 특성 곡선은 복잡한 원시 함수를 통해 표현된다.
• 특성 곡선 방정식이 $ y $ 에 대해 명시적으로 해석되지 않을 경우, 이분법을 사용하여 $ y(x,t) $ 를 수치적으로 구하고, 이로써 임의의 $ (x,t) $ 에서 해를 평가할 수 있다.