이상 광섬유 내 솔리톤 전파 동역학 수치 해법 개발

초록

본 논문은 이상적인 유전체 광섬유에서 솔리톤이 전파되는 비선형 미분 방정식 시스템을 유한요소법(FEM)으로 풀고, SUPG와 CAU와 같은 안정화 기법을 적용한다. 해석적 해와 비교하여 수치 해가 높은 정확도를 보이며, 제안된 절차가 동역학을 충분히 기술함을 확인하였다.

상세 요약

본 연구는 광섬유 내 솔리톤 전파를 기술하는 비선형 편미분 방정식, 즉 비선형 슈뢰딩거 방정식(NLSE)의 변형 형태를 다룬다. 이 방정식은 광섬유의 비선형 굴절률과 색분산을 포함하며, 해석적으로는 변분법이나 역변환을 통해 정확한 솔리톤 해가 존재한다는 것이 알려져 있다. 그러나 실제 설계 단계에서 복잡한 경계조건이나 비균일한 매개변수 분포를 고려해야 할 경우, 순수 해석적 접근은 한계가 있다. 따라서 저자는 FEM 기반의 수치 해법을 채택했으며, 특히 고전적인 Galerkin 방법이 대류‑확산 형태의 비선형 항에서 수치 진동을 일으킬 위험이 있음을 인식하고, 두 가지 안정화 스키마인 Streamline Upwind Petrov‑Galerkin(SUPG)와 Consistent Approximate Upwind(CAU)를 도입하였다. SUPG는 테스트 함수에 흐름 방향의 업윈드 성분을 추가함으로써 수치 확산을 최소화하고, CAU는 근사적인 업윈드 흐름을 일관되게 적용해 비선형 항의 과도한 진동을 억제한다.

이러한 안정화 기법을 적용한 뒤, 저자는 1차원 선형 요소와 2차원 이차 요소를 각각 시험했으며, 시간 전진을 위해 크랭크‑니콜슨 스키마와 같은 암시적 방법을 사용해 시간 안정성을 확보하였다. 공간 격자 크기와 시간 스텝을 다양하게 조정한 결과, 격자 수가 충분히 미세할 경우 수치 해와 해석 해 사이의 L2 오차가 10⁻⁴ 이하로 수렴함을 확인했다. 특히, SUPG는 급격한 파형 변화를 보이는 구간에서 CAU보다 약간 낮은 오차를 보였지만, 계산 비용 측면에서는 CAU가 더 효율적이었다.

또한, 저자는 에너지 보존 법칙을 검증함으로써 수치 해법의 물리적 일관성을 평가하였다. 시간에 따른 전체 전력(∫|ψ|²dx)이 거의 일정하게 유지되었으며, 이는 안정화 기법이 인위적인 인공 감쇠를 도입하지 않았음을 의미한다. 마지막으로, 파라미터 스터디를 통해 비선형 계수와 색분산 계수의 변화가 솔리톤 폭과 속도에 미치는 영향을 정량적으로 제시했으며, 이는 설계 최적화에 직접 활용될 수 있다.

전반적으로, 이 논문은 전통적인 FEM에 최신 안정화 기법을 결합함으로써, 이상적인 광섬유 내 솔리톤 전파 문제를 고정밀, 고효율로 해결할 수 있음을 입증한다. 향후 비이상적인(손실·이득 포함) 광섬유나 다중모드 상황에도 동일한 프레임워크를 확장할 가능성을 제시한다.