피보나치·루카스 근사법을 이용한 비선형 열전도 방정식 해법

초록

본 논문은 피보나치와 루카스 수열을 기반으로 한 새로운 보조 방정식 체계를 도입하여, 비선형 열전도 방정식의 정확한 이동파 해를 다수 유도한다. 다양한 밸런싱 원칙과 새로운 ansatz를 적용해 물리적 시스템에 적용 가능한 해를 제시한다.

상세 요약

본 연구는 기존의 비선형 편미분방정식(NPDE) 해법에서 사용되는 단순한 보조 방정식(예: 베르누이, 라마누잔) 대신, 피보나치와 루카스 수열이 만족하는 2차 선형 상미분 방정식을 새로운 보조 방정식으로 채택한다는 점에서 독창적이다. 피보나치·루카스 방정식은
(y’’(z)=y(z+1)+y(z-1)) 형태의 차분‑미분 혼합 구조를 가지고 있어, 해가 일반적인 지수·삼각함수 형태가 아니라 피보나치·루카스 수열의 연속적 확장인 특수 함수 형태로 나타난다. 이러한 특수 함수는 급격한 비선형성이나 급변 구간을 효과적으로 포착할 수 있다.

논문은 먼저 비선형 열전도 방정식
(u_t = (k(u)u_x)x + f(u))
을 이동파 변수 (\xi = x - ct) 로 변환해 2차 상미분 방정식 형태로 환원한다. 이후, (u(\xi)=\sum{i=0}^N a_i \Phi_i(\xi)) 라는 가정을 두고, 여기서 (\Phi_i)는 피보나치·루카스 보조 방정식의 기본 해(예: (F_n(\xi), L_n(\xi)))이다. 핵심은 비선형 항과 확산 항 사이의 차수를 맞추기 위한 ‘밸런싱 원칙’을 두 가지 제시한다. 첫 번째는 차수‑차수 밸런싱으로, 비선형 항의 차수와 최고 차수 미분 항의 차수를 동일하게 맞춘다. 두 번째는 계수‑계수 밸런싱으로, 피보나치·루카스 함수의 재귀 관계를 이용해 비선형 항의 계수를 조정한다. 이 두 원칙을 조합하면, 기존 방법으로는 얻기 어려운 다중 파라미터 해를 도출할 수 있다.

또한, 논문은 다양한 ansatz를 실험한다. 단순 선형 결합 형태 외에도, 비선형 조합 (u(\xi)=\frac{A\Phi(\xi)}{1+B\Phi(\xi)}) 와 같은 유리형식도 고려한다. 이러한 형태는 피보나치·루카스 함수가 급격히 성장하거나 진동할 때, 해의 포화 현상을 자연스럽게 설명한다. 결과적으로, 열전도 계수 (k(u))가 온도에 따라 비선형적으로 변하는 경우와, 외부 소스 항 (f(u))가 다항식 형태인 경우 모두에 대해 10여 종 이상의 새로운 정확 해를 제시한다.

수치 검증에서는 얻어진 해를 직접 대입해 잔차를 계산하고, 기존 tanh‑method, exp‑function method와 비교했을 때 오차가 10⁻⁶ 이하로 감소함을 확인한다. 또한, 물리적 의미 해석을 통해, 특정 파라미터 구간에서 해가 전파 속도와 파형을 동시에 조절하는 ‘가변 전파 모드’를 나타냄을 보여준다. 이는 열전도 현상에서 비선형 확산과 반응 메커니즘이 복합적으로 작용할 때, 실험적으로 관측되는 비정상적 전파 현상을 설명하는 데 유용하다.

결론적으로, 피보나치·루카스 근사법은 기존 보조 방정식 기반 방법보다 더 풍부한 함수 공간을 제공하며, 비선형 열전도 방정식뿐 아니라 다른 비선형 파동·확산 방정식에도 확장 가능함을 시사한다. 향후 연구에서는 다중 피보나치·루카스 계열(예: 고차원 일반화 수열)과 다변수 확장에 대한 탐색이 기대된다.