[논문 리뷰] The Riccati Differential Equation and a Diffusion-Type Equation

Erwin Suazo, Sergeĭ K. Suslov외 1명|ArXiv.org|2008. 07. 28.

Differential Equations and Numerical MethodsMathematics인용 수 19

한 줄 요약

이 논문은 시간에 의존하는 계수를 가진 일차원 확산형 미분방정식의 코시 초기값 문제에 대한 명시적 해석적 해를 제시한다. 문제를 리카티 방정식과 그 관련된 이阶 선형 특성 방정식으로 환원함으로써, 기본 함수(열핵)에 대한 닫힌 표현식을 초등 함수와 적분을 이용해 유도한다. 이는 동차 및 비동차 케이스에 대해 두하멜 원리를 통해 정확한 해를 가능하게 한다.

ABSTRACT

We construct an explicit solution of the Cauchy initial value problem for certain diffusion-type equations with variable coefficients on the entire real line. The corresponding Green function (heat kernel) is given in terms of elementary functions and certain integrals involving a characteristic function, which should be found as an analytic or numerical solution of the second order linear differential equation with time-dependent coefficients. Some special and limiting cases are outlined. Solution of the corresponding non-homogeneous equation is also found.

연구 동기 및 목표

• 실수선 위에서 변수 계수를 가진 일반적인 확산형 방정식의 코시 초기값 문제에 대한 명시적 해를 구성하는 것.
• 기본 함수(열핵)가 초등 함수와 적분으로 표현될 수 있는 조건을 규명하는 것.
• 확산 방정식의 해와 특성 함수로부터 유도된 리카티 방정식의 해 사이의 관계를 설정하는 것.
• 두하멜 원리를 적용하여 비동차 방정식으로 이 방법을 확장하는 것.
• hyperbolic 및 삼각함수 형태의 시간에 의존하는 계수를 가진 특수 케이스에 대한 명시적 해를 제공하는 것.

제안 방법

• 해는 $ u = A(t) e^{S(x,y,t)} $ 형태의 가설을 통해 구성되며, $ A(t) = 1 / \sqrt{2\pi\mu(t)} $ 이고, $ S $ 는 시간에 의존하는 계수를 가진 $ x $ 와 $ y $ 에 대한 이차형식이다.
• 가설을 PDE에 대입함으로써 유도된 ODE 시스템을 통해 $ S $ 의 계수를 결정하며, 이는 $ \alpha(t) $ 에 대한 리카티 방정식으로 이어진다.
• 리카티 방정식은 특성 함수 $ \mu(t) $ 에 대한 이계 선형 ODE로 변환되며, $ \mu'' - \tau(t)\mu' - 4\sigma(t)\mu = 0 $ 이다. 여기서 $ \tau $ 와 $ \sigma $ 는 시간에 의존하는 계수들로 정의된다.
• 기본 함수는 $ K(x,y,t) = \frac{1}{\sqrt{2\pi\mu(t)}} \exp(S(x,y,t)) $ 로 명시적으로 주어지며, 여기서 $ S $ 는 ODE 시스템의 해에 의존한다.
• 해는 초기 조건 $ \mu(0) = 0 $, $ \mu'(0) = 2a(0) \neq 0 $ 를 만족시키는 특성 방정식 $ \mu'' - \tau(t)\mu' - 4\sigma(t)\mu = 0 $ 의 해를 구함으로써 유도된다. 이는 초기 행동이 정확하게 반영되도록 보장한다.
• 비동차 방정식의 경우, 두하멜 원리가 적용되어 해가 동차 해의 연산자와 강제항을 포함하는 시간 적분 형태로 표현된다.

실험 결과

연구 질문

• RQ1 시간에 의존하는 계수를 가진 확산형 PDE의 기본 해가 초등 함수와 적분을 이용해 닫힌 형태로 표현될 수 있는 조건은 무엇인가?
• RQ2 시간에 의존하는 계수를 가진 확산 방정식의 해는 시간에 의존하는 계수를 가진 리카티 미분방정식의 해와 어떻게 관련이 있는가?
• RQ3 특성 함수 $ \mu(t) $ 는 어떤 역할을 하는가? 그리고 그 함수는 기본 함수의 구조를 어떻게 결정하는가?
• RQ4 두하멜 원리는 시간에 의존하는 계수를 가진 비동차 확산 방정식의 해를 유도하는 데 어떻게 적용될 수 있는가?
• RQ5 hyperbolic 및 삼각함수 형태의 시간에 의존하는 계수를 가진 특정 함수에 대해 열핵의 명시적 형태는 무엇인가?

주요 결과

• 확산 방정식의 기본 함수는 $ K(x,y,t) = \frac{1}{\sqrt{2\pi\mu(t)}} \exp(\alpha x^2 + \beta xy + \gamma y^2 + \delta x + \varepsilon y + \kappa) $ 로 명시적으로 구성되며, 계수들은 ODE 시스템의 해를 통해 결정된다.
• 계수 $ \alpha(t) $ 는 리카티 방정식 $ \frac{d\alpha}{dt} + b(t) - 2c(t)\alpha - 4a(t)\alpha^2 = 0 $ 를 만족하며, 이는 해법의 핵심이다.
• 특성 함수 $ \mu(t) $ 는 이계 선형 ODE $ \mu'' - \tau(t)\mu' - 4\sigma(t)\mu = 0 $ 를 만족하며, 여기서 $ \tau(t) = \frac{a'}{a} + 2c - 4d $, $ \sigma(t) = ab + cd - d^2 + \frac{d}{2}\left(\frac{a'}{a} - \frac{d'}{d}\right) $ 이다.
• 특수 케이스 $ a(t) = \sinh^2 t $, $ b(t) = \cosh^2 t $, $ c(t) = \frac{1}{2}\sinh 2t $ 에서 특성 방정식은 $ \mu'' - 2\tanh t\, \mu' + 2\mu = 0 $ 가 되며, 해는 $ \mu_1(t) = \cos t \sinh t + \sin t \cosh t $ 로 주어지며, $ 0 < t < T_1 \approx 0.93755 $ 에서 유효하다.
• 이 경우 열핵은 수렴 조건이 $ \gamma(t) $ 의 부호와 연결된 형태로 명시적으로 주어지며, 적절한 초기 자료에 대해 적분 가능성을 보장한다.
• 유사한 해는 $ a(t) = \cos^2 t $, $ b(t) = \sin^2 t $, $ c(t) = -\frac{1}{2}\sin 2t $ 인 경우에도 유도되며, 열핵의 형태는 동일하지만 $ x $ 와 $ y $ 가 서로 바뀌어 있으며, $ 0 < t < T_2 \approx 2.347 $ 에서 유효하다.