에너지 보존을 위한 두 유한차분 스킴과 비선형 초전도 결합 진동계 연구

본 논문은 반감쇠 사인‑고든 방정식으로 모델링되는 1차원 반무한 조셉슨 접합 배열을 대상으로, 에너지 보존과 소산을 정확히 재현하는 두 가지 유한차분 스킴을 개발하고 그 수치적 특성을 체계적으로 분석한다. 1. **서론 및 물리적 배경** 사인‑고든 방정식은 장거리 조셉슨 접합, 비선형 파동 전파, 강체 펜듈럼 등 다양한 물리 현상을 기술한다. 특히, 내부·외부 감쇠(β, γ), 질량(m), 전류 파라미터(J, I) 등을 포함한 변형 형태는 비선형 초전도 초전송(supratransmission) 현상을 설명하는 데 핵심이다. 기존 연구는 보존형 케이스에 집중했으며, 감쇠가 포함된 경우는 아직 충분히 다루어지지 않았다. 2. **수학적 모델** 연속 모델은 N개의 조셉슨 접합이 서로 인접한 차분 연산 Δ₂ₓ와 시간 2차 미분으로 연결된 시스템(식 1)으로 제시된다. 경계 조건은 입력 전압 φ(t)=A sin(Ωt)와 출력 전류 I=ẋ_N/R 로 구성된다. 선형화 시 파동 전파 관계 ω²(k)=m²+1+2c²(1−cos k) 로부터 금지대역 Ω<√(m²+1) 를 도출한다. 3. **에너지 식** 보존형 경우 각 사이트의 해밀토니안 Hₙ(식 4)와 전체 에너지 E(식 5)를 정의하고, 감쇠가 존재할 때 총 에너지의 시간 미분을 식 6으로 전개한다. 이 식은 경계 입력, 내부 감쇠, 외부 감쇠 항을 명시적으로 포함한다. 무한 배열에 대한 극한을 취하면 Corollary 2와 같이 경계 전류만이 에너지 흐름을 결정한다는 중요한 결과를 얻는다. 4. **유한차분 스킴** - **명시적 스킴(식 9)**: 시간 중심 차분 δ²ₜ와 공간 중심 차분 δ²ₓ를 사용하고, 평균 연산자 ¯δ와 중앙 차분 δₜ를 결합해 2차 정확도를 확보한다. 비선형 포텐셜 V(u)의 차분은 (V(uₖ₊₁)−V(uₖ₋₁))/(uₖ₊₁−uₖ₋₁) 형태로 처리한다. 경계 조건은 φₖ를 이용해 첫 번째 노드에 적용한다. - **암시적 스킴(식 13)**: δ²ₓ를 4점 스템으로 확장해 β·Δₜδ²ₓ 항을 O(Δt²) 정확도로 근사한다. 이는 감쇠가 큰 경우에도 수치적 안정성을 크게 향상시킨다. 두 스킴 모두 뉴턴-크루트 반복과 삼중 대각 행렬 해법(Crout)으로 비선형 연립 방정식을 해결한다. 5. **에너지 보존 및 소산 검증** 각 스킴에 대해 이산 해밀토니안 Hₖₙ(식 11, 14)과 총 에너지 Eₖ(식 12, 15)를 정의하고, 차분 형태의 에너지 변화율을 식 17(명시적) 및 식 23(암시적)으로 도출한다. 두 식은 연속식(6)과 완전히 일치하며, 내부·외부 감쇠 항이 음의 기여를, 경계 입력이 양의 기여를 하는 것을 명시한다. 따라서 제안된 스킴은 “에너지 보존(또는 정확한 소산)” 특성을 갖는다. 6. **안정성 분석** 선형화된 시스템에 대해 von Neumann 분석을 수행하여 Δt와 Δx 사이의 CFL 조건을 도출한다. 명시적 스킴은 Δt ≤ 2/√(c²·(2−cos k)+m²+1) 의 제한을 받으며, 이는 감쇠가 클수록 더욱 엄격해진다. 암시적 스킴은 β, γ에 관계없이 무조건적인 안정성을 보이며, 실제 시뮬레이션에서도 Δt를 비교적 크게 잡아도 수렴한다. 7. **수치 실험** 파라미터 β, γ, m, J, I 를 다양하게 변동시키며 초전도 초전송 임계 진폭 Aₜₕ를 측정한다. 주요 결과는 다음과 같다. - 내부·외부 감쇠가 증가하면 Aₜₕ가 크게 상승, 즉 초전송이 억제된다. - 질량 m이 실수이면 금지대역이 넓어져 Ω가 작을 때도 초전송이 어려워진다. - 전류 J와 출력 전류 I가 비대칭일 경우, 에너지 흐름이 비대칭적으로 변형되어 초전송 임계가 비선형적으로 변조된다. - 제시된 스킴을 사용한 장기 시뮬레이션에서도 에너지 소산이 연속식과 일치함을 확인하였다. 8. **결론 및 향후 연구** 두 유한차분 스킴은 해와 에너지, 에너지 소산률 모두에서 2차 정확도를 유지하면서, 명시적 스킴은 조건부 안정, 암시적 스킴은 무조건적 안정이라는 장점을 제공한다. 이는 비선형 초전도 시스템의 장기 동역학, 초전송 임계값 탐색, 그리고 파라미터 스위프와 같은 고비용 계산에 매우 유용하다. 향후 연구에서는 2차원·3차원 배열, 비동질성(불균일한 파라미터) 및 외부 잡음 효과를 포함한 확장 모델에 적용할 계획이다.