[논문 리뷰] Fractional Calculus in Wave Propagation Problems

Francesco Mainardi|arXiv (Cornell University)|2012. 02. 01.

Fractional Differential Equations Solutions참고 문헌 37인용 수 8

한 줄 요약

이 논문은 선형 점탄성 매질에서 확산과 고전적 파동 운동 사이의 중간 파동 전파 현상을 기술하기 위해 분수적 미분을 프레임워크로 제안한다. 파동 방정식에서 정수계 수식의 시간 도함수를 분수계 도함수(0 < β ≤ 2)로 대체함으로써, 이 모델은 거듭제곱 법칙에 따르는 찌그러짐 행동을 포착하며, M-Wright 함수로 지배되는 해를 도출하여 확산적 행동과 진동적 행동 사이를 보간한다.

ABSTRACT

Fractional calculus, in allowing integrals and derivatives of any positive order (the term "fractional" kept only for historical reasons), can be considered a branch of mathematical physics which mainly deals with integro-differential equations, where integrals are of convolution form with weakly singular kernels of power law type. In recent decades fractional calculus has won more and more interest in applications in several fields of applied sciences. In this lecture we devote our attention to wave propagation problems in linear viscoelastic media. Our purpose is to outline the role of fractional calculus in providing simplest evolution processes which are intermediate between diffusion and wave propagation. The present treatment mainly reflects the research activity and style of the author in the related scientific areas during the last decades.

연구 동기 및 목표

• 확산과 파동 전파 사이의 중간 진화 과정을 기술하는 데 있어 분수적 미분의 역할을 조사하는 것.
• 거듭제곱 법칙에 따르는 찌그러짐 행동를 보이는 선형 점탄성 재료에서의 파동 전파를 분석하는 것. 이러한 재료는 점성 유체와 탄성 고체 사이의 행동을 보인다.
• 시간 도메인에서 분수계 수준 β(0 < β ≤ 2)의 도함수가 표준 확산 방정식과 파동 방정식을 일반화하는 적분미분 방정식을 도출하는 방법을 보여주는 것.
• 표준 초깃값 조건과의 호환성 덕분에 물리적 응용에서 초깃값 문제에 가장 적합한 도구로 캡투(Caputo) 분수 도함수를 설정하는 것.
• 분수 파동 방정식의 기본 해로서 M-Wright 함수의 해석적 및 점근적 성질을 탐색하는 것.

제안 방법

• 표준 파동 방정식 또는 확산 방정식에서 일阶 또는 이阶 시간 도함수 대신 캡투 분수 도함수(0 < β ≤ 2)를 사용하여 일반화된 진화 방정식을 수립한다.
• 정수계 도함수를 포함한 물리적 초깃값 조건과의 호환성을 확보하기 위해 분수 도함수의 캡투 정의를 사용한다.
• 분수 적분을 적용하여 분수 파동 방정식의 등가 적분미분 형태를 유도하고, 초깃값을 명시적으로 포함한다.
• 분수 미분 방정식에 라플라스 변환을 적용하여 주파수 도메인에서의 해 표현식을 유도한다.
• M-Wright 함수 Mν(x)를 사용하여 기본 해를 특성화하며, 여기서 ν = β/2 − 1/2이며, 해석적 점근적 분석과 수치적 플롯을 통해 그 행동을 분석한다.
• 선형 척도와 로그 척도를 사용하여 ν 값(0에서 1 사이)에 따른 M-Wright 함수의 행동을 비교하고, ν → 1⁻에 해당하는 경우 딜라크 델타 함수에 수렴하는 극한을 분석한다.

실험 결과

연구 질문

• RQ1 분수적 미분은 어떻게 확산과 고전적 파동 운동 사이의 중간 파동 전파 과정을 모델링할 수 있는가?
• RQ2 파동 전파 문제에서 1계 또는 2계 시간 도함수를 분수계 도함수 β(0 < β ≤ 2)로 대체할 경우의 물리적 의미는 무엇인가?
• RQ3 왜 점탄성 파동 모델링의 초깃값 문제에서는 리만-리우빌 도함수보다 캡투 도함수가 선호되는가?
• RQ4 M-Wright 함수 Mν(x)는 어떻게 분수 파동 방정식의 기본 해를 기술하며, 그 형태는 β에 따라 어떻게 변화하는가?
• RQ5 특히 ν → 1⁻일 때, M-Wright 함수의 원점 근처 및 무한대에서의 점근적 및 수치적 행동은 어떠한가?

주요 결과

• 캡투 도함수를 갖는 분수 파동 방정식은 β < 1일 경우 확산적 행동, β > 1일 경우 진동적 행동 사이를 보간하는 해를 도출한다.
• β = 1일 경우 방정식은 표준 확산 방정식으로 간소화되며, β = 2일 경우 고전적 파동 방정식으로 간소화된다.
• 분수 파동 방정식의 기본 해는 M-Wright 함수 Mν(x)로 표현되며, 여기서 ν = β/2 − 1/2이다. 이 함수는 ν = 0일 때 지수 감쇠에서 시작하여 ν = 1/2일 때 가우시안 형태로 변화하고, 최종적으로 ν = 1일 때는 딜라크 덴타 함수로 수렴한다.
• β → 2⁻(ν → 1⁻)일 때 M-Wright 함수는 x ≈ ±1 근처에 날카운 피크를 형성하며, 유한한 전파 속도를 갖는 파동 유사 전파 행동이 나타남을 나타낸다.
• ν가 1에 가까운 경우(예: ν = 0.99 또는 ν = 0.999)의 Mν(x) 수치 플롯은 x ≈ 1 근처에 고고도, 좁은 피크를 보이며, 파동 유사 행동으로의 전이를 확인한다.
• 캡투 도함수의 라플라스 변환은 초깃값의 물리적 형태를 유지하므로, 측정 가능한 初기 상태를 갖는 실제 점탄성 시스템을 모델링하는 데 이상적이다.