이상 광섬유 솔리톤 전파를 위한 유한 차분·가우스‑시델 수치 해법

본 논문은 광통신 분야에서 솔리톤을 이용한 장거리 전송 기술의 핵심 이론과 수치 해법을 종합적으로 다룬다. 서론에서는 1970년대 Hasegawa와 Tappert가 제안한 비선형 파동 보상 개념을 시작으로, 1980‑2000년대에 걸친 실험적 진전과 현재 상용화 단계까지의 흐름을 개괄한다. 특히, 솔리톤 전파를 기술하는 기본 방정식으로 비선형 슈뢰딩거 연립 방정식(식 1)을 채택하고, 이 방정식이 광섬유의 군속 분산(α), 두 번째 고조파 발생률(β), 파동 간 속도 차이(δ), 그리고 정상·비정상 분산을 구분하는 매개변수 r에 의해 어떻게 정의되는지를 상세히 설명한다. 이후, 이상 광섬유(δ=0) 상황에서 해석적 솔리톤 해(식 2‑3)를 제시하고, 이 해가 파라미터 α, β, r에 따라 파동의 폭, 진폭, 전파 속도를 어떻게 결정하는지를 수식적으로 분석한다. 이러한 해석적 해는 수치 모델의 검증 기준으로 활용된다. 수치 모델링 파트에서는 유한 차분법을 적용해 시간(ξ) 방향은 전진 차분, 공간(s) 방향은 중앙 차분으로 근사한다. 이산화 과정에서 복소 변수 a₁, a₂를 실수·허수 성분으로 분리해 2 × 2 블록 행렬 형태의 선형 시스템(식 4)을 도출한다. 시스템은 암시적이므로 각 시간 단계마다 선형 방정식 집합을 풀어야 하며, 이를 위해 완화 가우스‑시델 반복법을 채택한다. 완화 인자 ω=0.1은 실험적으로 최적화된 값으로, 수렴 속도와 안정성 사이의 균형을 맞춘다. 반복 종료 기준은 상대 오차가 6 × 10⁻⁴ 이하가 될 때까지이며, 이는 실제 시뮬레이션에서 10⁴ ~ 10⁵ 회 반복으로 달성된다. 구현은 Fortran 90 기반으로 작성되었으며, 코드 흐름도(그림 2)를 통해 입력 파라미터 설정, 초기 조건(해석적 솔리톤), 경계 조건(충분히 큰 L, T), 이산화 단계, 가우스‑시델 반복, 결과 저장 순서를 명확히 제시한다. 실험 환경은 AMD Athlon 64X2 듀얼코어(2.29 GHz)와 1 GB RAM을 갖춘 PC이며, 격자 간격 Δξ=Δs=0.001을 사용해 L=50, T=10 구간을 5 × 10⁴개의 점으로 discretize하였다. 수치 결과에서는 기본 파와 제2 고조파 모두에 대해 해석적 솔리톤과 거의 일치하는 파형을 얻었다. 최대 절대 오차는 각각 3 × 10⁻⁵, 3 × 10⁻⁷ 수준이며, 이는 그래프(그림 3‑4)에서 눈에 띄게 확인된다. 전체 연산 시간은 약 1 × 10⁵ 초(≈28 시간)로, 현재 구현의 계산 효율성을 보여준다. 또한, δ≠0인 경우에도 동일한 프레임워크를 적용할 수 있음을 언급하며, 파라미터 δ가 증가하면 시스템의 강인성이 감소하고 수렴 속도가 느려질 수 있음을 지적한다. 향후 연구로는 다중‑그리드, 적응형 격자, 그리고 고차 차분 스킴을 도입해 계산 비용을 절감하고 비이상 섬유의 복잡한 물리 현상을 더 정밀히 모델링하는 방안을 제시한다. 결론에서는 제안된 유한 차분·가우스‑시델 수치 절차가 이상 광섬유 내 솔리톤 전파를 정확히 재현함을 확인하고, 이를 기반으로 비이상 섬유, 다중 파장 전송, 그리고 실시간 시스템 시뮬레이션 등 다양한 응용 분야에 확장 가능함을 강조한다.