가변 밀도 현의 진동 모드: 정밀 근사와 비정규 해법

초록

본 논문은 양끝이 고정된 가변 밀도 현의 진동 고유모드를 세 가지 정리 기반의 반복 알고리즘으로 근사한다. 제시된 방법은 밀도 분포와 무관하게 적용 가능하며, 구체적인 예제에서 WKB 해와 콜로케이션 수치 해와 비교해 높은 정확도를 보인다. 또한 2차 섭동 이론으로 얻은 에너지의 비대칭적 차수가 고차원 경우와 다른 형태임을 확인한다.

상세 분석

이 연구는 비균일한 밀도를 가진 현의 고유진동 문제를 해석적으로 다루기 위해 세 가지 핵심 정리를 도입한다. 첫 번째 정리는 변분 원리를 이용해 현의 고유값을 Rayleigh 비율 형태로 표현하고, 임의의 시험함수에 대해 상한·하한을 제공한다는 점에서 기존의 Sturm‑Liouville 이론을 일반화한다. 두 번째 정리는 이러한 시험함수를 이용해 고유함수를 점진적으로 개선하는 반복 공식(Iterative Improvement Scheme)을 제시한다. 구체적으로, 현재 근사함수 ψ_n(x)를 입력으로 하여 새로운 함수 ψ_{n+1}(x)=∫_0^x G(x,ξ) ρ(ξ) ψ_n(ξ)dξ 형태의 적분 연산을 적용하면, 고유값과 고유함수가 수렴한다는 것이 증명된다. 여기서 G는 고정된 경계조건을 만족하는 Green 함수이며, ρ(x)는 현의 밀도 함수이다. 세 번째 정리는 위 반복 과정이 L^2 노름에서 선형 수렴성을 갖는다는 것을 보이며, 수렴 속도는 밀도 변동의 스무스함과 초기 시험함수의 선택에 의존한다는 점을 강조한다.

논문은 이러한 정리를 바탕으로 두 가지 대표적인 밀도 프로파일—(i) 선형적으로 증가하는 밀도 ρ(x)=ρ_0(1+αx)와 (ii) 급격히 변하는 단계형 밀도—에 대해 구체적인 계산을 수행한다. 각 경우에 대해 1차·2차 섭동 전개와 반복 알고리즘을 동시에 적용해 고유주파수 ω_n을 구한다. 섭동 전개에서는 ρ(x)−⟨ρ⟩를 작은 파라미터 ε로 두고, 고유값 λ_n=ω_n^2에 대해 λ_n=λ_n^{(0)}+ε λ_n^{(1)}+ε^2 λ_n^{(2)}+… 형태로 전개한다. 여기서 λ_n^{(0)}는 균일 밀도 현의 고유값이며, λ_n^{(1)}, λ_n^{(2)}는 각각 밀도 변동의 1차·2차 보정항이다. 특히 2차 보정항 λ_n^{(2)}는 밀도 함수의 제곱 평균과 고유함수의 공간적 구조가 복합적으로 작용함을 보여준다.

WKB(윅켈) 근사는 고주파(큰 n) 영역에서 λ_n≈π^2 n^2 /∫_0^L ρ^{-1/2}(x)dx 라는 비정규적 형태를 제시한다. 저자들은 섭동 이론으로부터 얻은 2차 보정이 이 WKB 결과와 정확히 일치함을 증명함으로써, 두 방법이 동일한 비대칭적 차수를 공유한다는 점을 확인한다. 흥미롭게도, 2차 섭동 결과는 ρ(x) 의 평균값이 아닌, ρ^{-1/2}(x) 의 평균에 의존한다는 점에서 2차원·3차원 파동 방정식에서 나타나는 ρ^{−1} 평균과는 다른 스케일링을 보인다. 이는 차원에 따라 밀도에 대한 효과적인 가중치가 달라진다는 물리적 의미를 내포한다.

수치 검증을 위해 저자들은 콜로케이션 방법(특히 Chebyshev‑Gauss‑Lobatto 격자 기반)으로 직접 고유값 문제를 풀었다. 비교 결과, 반복 정리 기반 근사는 n=110 범위에서 상대 오차가 10^{-5} 이하로, WKB 근사보다 현저히 정확했으며, 특히 저주파 모드에서 큰 차이를 보였다. 또한, 초기 시험함수를 단순한 사인 함수 sin(π n x/L) 로 선택했음에도 불구하고, 56 회의 반복만으로 수렴이 이루어졌다. 이는 제시된 정리들의 강인함과 일반성을 뒷받침한다.

마지막으로, 저자들은 이론적 확장을 논의한다. 현재의 정리와 알고리즘은 1차원 Sturm‑Liouville 형태의 문제에 국한되지 않고, 변형된 질량·강성 분포를 갖는 막대, 원통형 파이프, 심지어 비선형 탄성체에도 적용 가능하다는 가능성을 제시한다. 특히, 고차원(2D, 3D) 파동 방정식에서 밀도와 강성의 비정규적 결합이 나타날 때, 섭동 전개의 차수가 현재 1D 결과와 어떻게 달라지는지를 탐구하는 것이 향후 연구 과제로 제시된다.