테일러 급수로 푸는 원자로 점동역학 방정식

초록

본 논문은 테일러 급수 전개를 이용해 원자로 점동역학 방정식의 수치 해법을 제시한다. 1차 전개만으로도 단계, 램프, 사인파, 지그재그 형태의 반응도 입력에 대해 기존 고급 방법들과 동등한 정확도를 얻는다. 알고리즘은 구현이 간단하고 투명하며, 다양한 초기 조건과 입력 데이터를 통해 검증되었다.

상세 요약

원자로 점동역학 방정식은 중성자 밀도와 전구체 농도의 시간 의존성을 기술하는 강하게 비선형인 미분방정식이다. 전통적으로는 정밀도가 높은 런지-쿠타, 선형 다항 근사, 혹은 고차 차분법 등이 사용되어 왔지만, 구현 복잡도와 연산 비용이 문제점으로 지적되어 왔다. 본 연구는 테일러 급수를 시간 단계마다 1차까지 전개함으로써, 각 단계에서 현재 시점의 미분값만을 이용해 다음 시점의 값을 직접 계산한다는 매우 단순한 접근법을 제안한다. 핵심 아이디어는 중성자 밀도 n(t)와 전구체 농도 C_i(t) 를 각각
n(t+Δt) ≈ n(t) + Δt·dn/dt|_t,
C_i(t+Δt) ≈ C_i(t) + Δt·dC_i/dt|_t
와 같이 1차 테일러 전개로 근사하고, dn/dt, dC_i/dt는 점동역학 방정식 자체에서 바로 얻는다. 이때 시간 간격 Δt는 충분히 작게 잡아야 하지만, 실험 결과는 Δt = 0.001 s 정도에서도 높은 정확도를 유지한다는 것을 보여준다.

알고리즘의 장점은 첫째, 미분식 자체가 선형 형태이므로 전구체 수가 늘어나도 연산량이 O(N) 수준에 머문다. 둘째, 각 단계에서 새로운 비선형 방정식을 푸는 과정이 없으므로 수치적 안정성이 뛰어나며, 특히 급격한 반응도 변화(스텝, 램프)에도 오버슈팅 없이 정확히 추적한다. 셋째, 구현 코드가 직관적이어서 교육용이나 빠른 프로토타이핑에 적합하다.

비교 실험에서는 고전적인 프루드-헨리, 고차 루프-카우프만, 그리고 고정점 반복법과 결과를 대조하였다. 스텝 반응도(ρ = 0.003)에서는 중성자 밀도 피크와 감쇠율이 모두 0.1 % 이하의 오차로 일치했으며, 램프 반응도(ρ = 0.001·t)에서도 시간 연속성이 유지되었다. 사인파(ρ = 0.0005·sin ωt)와 지그재그(ρ = ±0.001 교대로) 입력에서는 주파수 스펙트럼 분석 결과, 테일러 1차 해법이 고차 방법과 거의 동일한 위상과 진폭을 보였다.

한계점으로는 Δt가 지나치게 크게 설정될 경우 1차 전개에 의한 truncation error가 누적되어 정확도가 떨어진다. 또한, 급격한 재활용(재활용 계수 β가 큰 경우)에서는 2차 이상 전개가 필요할 수 있다. 그러나 이러한 경우에도 테일러 전개를 2차까지 확장하면 기존 방법과 동일한 수준의 정확도를 손쉽게 확보할 수 있다.

결론적으로, 테일러 급수 기반 1차 전개는 복잡한 수치 기법을 대체할 수 있는 충분히 강력하고 효율적인 도구이며, 특히 실시간 제어 시스템이나 교육용 시뮬레이터에 적합한 해법으로 평가된다.