HarmonicBalance.jl: A Julia suite for nonlinear dynamics using harmonic balance

기술 요약: "HarmonicBalance.jl: 조화 균형법을 이용한 비선형 역학용 Julia 수트"

문제 정의
비선형 상미분 방정식(ODE)은 유체 역학, 로보틱스, 광기계 공학, 진동자 네트워크에 이르기까지 다양한 물리적 시스템의 시간 진화를 지배합니다. 선형 시스템은 일반적으로 고유한 정상 상태(steady state)로 수렴하지만, 비선형 시스템은 다수의 정상 상태 해를 가질 수 있으며, 이는 히스테리시스(hysteresis), 자발적 대칭성 깨짐, 노이즈 유도 전환(noise-induced switching)과 같은 현상으로 이어집니다. 이러한 시스템을 분석하는 데 있어 중요한 과제는 표준 수치 ODE 솔버(초기값 문제)가 실행당 단 하나의 정상 상태만을 찾는다는 점입니다. 따라서 초기 조건 공간을 완벽하게 탐색하려면 무한한 샘플링이 필요하며, 이는 계산적으로 불가능합니다. 또한, 구동되는 시스템을 회전 좌표계(rotating frame)로 변환하면 구동 항을 정적으로 만들 수 있지만, 결과로 발생하는 결합된 비선형 대수 방정식의 모든 근을 찾는 것은 여전히 어렵습니다. 왜냐하면 표준 근 찾기 알고리즘(예: 뉴턴-랩슨 방법)은 초기 추측값 근처의 해만을 찾아내기 때문입니다.

방법론
저자들은 조화 구동되는 비선형 시스템의 모든 정상 상태 해를 찾기 위해 설계된 Julia 패키지인 HarmonicBalance.jl을 소개합니다. 이 방법론은 세 가지 핵심 구성 요소를 통합합니다:

1. 조화 균형 안사츠(Harmonic Balance Ansatz): 이 패키지는 충분히 긴 시간이 흐른 후, 시스템의 응답이 서서히 변화하는 진폭을 가진 유한한 조화 항들의 합으로 근사될 수 있다고 가정합니다. 조화 시간 의존성을 갖는 $N$개의 2차 ODE 시스템에 대해, 해 $x_i(t)$는 다음과 같이 전개됩니다:
   $$x_i(t) = \sum_{j=1}^{M} [u_{i,j}(T)\cos(\omega_{i,j}t) + v_{i,j}(T)\sin(\omega_{i,j}t)]$$
   여기서 $T$는 거친 규모(coarse-grained)의 시간 척도입니다. 이 안사츠를 원래의 ODE에 대입하고 빠른 진동(진폭의 2차 시간 미분 항)을 제거함으로써, 문제는 $u$와 $v$에 대한 일련의 자율 1차 ODE 세트로 축소됩니다. 정상 상태는 이 변수들의 시간 미분을 0으로 설정함으로써 발견됩니다. 이는 결합된 다항 대수 방정식 시스템이 됩니다.

2. 호모토피 연속법(Homotopy Continuation): 결과로 얻은 다항 시스템을 풀기 위해, 본 패키지는 HomotopyContinuation.jl 라이브러리를 통해 호모토피 연속법을 채택합니다. 반복적인 근 찾기 방법과 달리, 호모토피 연속법은 해석적으로 풀 수 있는 시스템(알려진 근의 개수를 가진)을 목표 시스템으로 변형시킵니다. 이 접근 방식은 모든 복소수 근을 발견하는 것을 보장합니다. 패키지는 계산 비용을 관리하기 위해 2단계 알고리즘을 사용합니다: 일반적인 시스템으로부터 경로를 추적하여 특이 경로(singular paths)를 걸러내는 "웜업(warm-up)" 단계와, 남은 비특이 경로를 물리적 해로 추적하는 "파라미터 호모토피(parameter homotopy)" 단계가 있습니다.

3. 기호 및 수치 적분: 패키지는 조화 방정식과 야코비안(Jacobian) 행렬을 자동으로 유도하기 위해 기호 조작을 위한 Symbolics.jl을 사용합니다. 또한 시간 의존적 시뮬레이션을 위해 DifferentialEquations.jl과 인터페이스하여, 사용자가 정상 상태 결과를 검증하고 단열 파라미터 스윕(adiabatic parameter sweeps, 히스테리시스)을 연구할 수 있도록 합니다.

주요 기여

• 포괄적인 해 탐색: 본 패키지의 주요 기여는 초기 조건에 의존하는 단일 해가 아니라, 조화 구동되는 비선형 시스템에 대한 모든 정상 상태 해(실수 및 복소수 해)를 수치적으로 결정할 수 있다는 점입니다.
• 통합 프레임워크: 패키지는 기호 대수, 조화 균형, 호모토피 연속법을 단일 오픈 소스 Julia 환경으로 통합하여 높은 성능과 유연성을 제공합니다.
• 안정성 및 변동 분석: 야코비안 행렬의 고윳값을 기반으로 해를 실수 또는 복소수로 자동 분류하고, 안정 또는 불안정 상태를 판별하는 도구를 제공합니다. 또한 안정적인 고정점에 대한 선형 응답 스펙트럼 계산을 지원합니다.
• 시각화: 수트에는 1D 및 2D 위상도(phase diagrams), 해의 분기(solution branches), 분기 구조(bifurcation structures)를 시각화하기 위한 그래픽 사용자 인터페이스와 플로팅 루틴이 포함되어 있어 복잡한 해 토폴로지 탐색을 용이하게 합니다.

결과 및 예시
논문은 다음과 같은 예시를 통해 패키지의 능력을 입증합니다:

• 구동되는 더핑 진동자(Driven Duffing Oscillator): 패키지는 비선형성에 의해 생성되는 상조화 응답($3\omega_d$)을 포함하여 단일 더링 진동자의 다중 정상 상태 해를 성공적으로 식별합니다. 이는 표준 공명 곡선과 해가 하나, 둘 또는 셋인 영역을 보여주는 위상도를 재현합니다.
• 결합된 더링 공진기(Coupled Duffing Resonators): 두 개의 결합된 진동자 시스템에 대해, 패키지는 파라미터 공간 전체에서 해의 지형을 매핑하여 복잡한 토폴로지를 드러냅니다. 이는 호모토피를 통해 발견된 "완전한" 해 집합과, 단열 시간 스윕을 통해 접근 가능한 제한된 부분 집합 사이의 차이를 대조하며, 히스테리시스와 분기 간의 점프를 보여줍니다.
• 매개변수 구동 진동자(Parametrically Driven Oscillator): 패키지는 매개변수 구동 더링 진동자(마티외 방정식)를 분석하여, 단일 및 이중 안정 해 영역을 식별하고 안정 상태의 수가 변하는 피치포크 분기(pitchfork bifurcation)를 특징짓습니다.
• 성능 스케일링: $N$개의 결합된 더링 진동자 체인에 대한 벤치마크 결과, 베주 경계(Bézout bound, 이론적 최대 해의 수)는 지수적으로($3^{2N}$) 증가하지만, 물리적인 실수 해의 수는 이보다 훨씬 적음을 보여줍니다. 패키지는 단일 CPU 코어에서 몇 분 안에 $N=5$(수백 개의 경로 포함)인 시스템을 해결할 수 있어, 경로 추적 전략의 효율성을 입증합니다.

의의
저자들은 HarmonicBalance.jl을 이론적 분석과 실험적 관찰 사이의 간극을 메우는 도구로 정의합니다. 전체 해의 지형을 탐색할 수 있게 함으로써, 이 패키지는 표준 시간 진화 시뮬레이션에서는 놓치기 쉬운 히스테리시스 및 대칭성 깨짐과 같은 실험적으로 관찰 가능한 현상을 식별할 수 있게 합니다. 이 소프트웨어는 임의의 운동 방정식을 특수한 GUI 제약 없이 정의할 수 있는 "ab-initio" 접근 방식을 제공하여 학술 커뮤니티가 쉽게 접근할 수 있도록 설계되었습니다. 저자들은 이 도구가 구조 역학, 비선형 광학, 광기계 공학, 이징 머신(Ising machines)을 포함한 광범위한 분야에 적용 가능하다고 언급했습니다. 현재 버전은 결정론적 역학과 정상 상태에 집중하고 있으나, 모듈형 설계는 리미트 사이클(limit cycle) 탐지 및 고차 평균화 방법과 같은 향후 확장을 지원하도록 의도되었습니다.