순환연료 원자로 확률점동학 모델 완전 혼합 근사

초록

본 논문은 순환연료 원자로(CFR)에서 지연중성자 전구체(DNP)의 흐름을 완전 혼합 가정 하에 모델링하고, 이를 확률점동학(점-동역학) 형태로 전환한다. 두 개의 완전 혼합 체적(코어와 외코어)로 구성된 시스템을 이산 사건 과정과 Itô 확률미분방정식(SDE)으로 기술하고, 아날로그 몬테카를로(AMC)와 반암시적 Milstein SDE 해석기를 구현한다. 벤치마크 결과는 평균값이 결정론적 해와 일치함을 보이며, SDE 해석기가 특정 상황에서 DNP 분산을 과소평가한다는 한계를 제시한다. 또한 전구체 이동에 따른 반응성 손실 추정식이 편향됨을 분석한다.

상세 분석

이 연구는 기존의 정적 연료 원자로에 적용된 확률점동학 모델을 순환연료 시스템에 확장하는 데 초점을 맞춘다. 핵심 아이디어는 코어와 외코어 두 구역을 각각 완전 혼합으로 가정하고, 전구체의 체류시간 τ_c, τ_e 로 흐름을 묘사함으로써 지연항을 없애는 것이다. 이렇게 하면 전통적인 지연 미분식(2)의 메모리 효과를 제거하고, 현재 상태만으로 미래를 기술할 수 있는 마코프 과정으로 전환한다.

이산 사건 모델(Table 1)에서는 중성자 손실(포획·누출), 즉각적 및 지연 중성자 발생, 외부 소스, 전구체 붕괴·전이 등을 각각 확률적 사건으로 정의하고, 각 사건의 발생률을 N·γ(t), N·ϕ, S(t) 등으로 표현한다. 여기서 ϕ와 γ(t)는 점-동역학 파라미터 β, Λ, ρ(t) 로부터 유도된다(식 6‑7). 이러한 정의는 Poisson 과정에 기반한 확률적 전이 행렬을 구성하게 하며, 연속시간 한계에서 Itô SDE 로 수렴한다.

SDE 형태(식 8‑10)는 중성자 수 N_t 에는 확산항 D(t)·|N_t|·dW_t 가 포함되지만, 전구체 C_c, C_e 에는 확산항이 없다는 점이 핵심적인 한계로 지적된다. 저자들은 D(t) 를 ν_p 의 분산과 손실률 γ(t) 로 정의하고, 절대값을 취해 음수값 발생을 방지한다. 그러나 전구체의 확산항을 생략함으로써 DNP 변동성이 실제보다 작게 예측되는 현상이 관찰된다. 이는 특히 전구체 수가 적은 그룹(예: 그룹 1)에서 두드러진다.

수치 구현에서는 AMC 엔진(MARS)와 반암시적 Milstein 스키마를 사용한다. MARS는 사건 기반 시뮬레이션으로, 각 사건을 개별적으로 샘플링해 정확한 통계량을 제공한다. 반면 Milstein 해석기는 시간 스텝을 비균등하게 조정하고, 행렬 연산을 사전 계산해 효율성을 높인다. 두 방법 모두 평균값은 결정론적 해와 완벽히 일치했지만, 분산에서는 차이가 나타났다. 저자들은 이 차이가 SDE 모델에서 DNP 노이즈가 누락된 것에 기인한다고 추정한다.

또한 전구체 이동에 따른 반응성 손실 ρ₀(식 18)를 추정하는 실험을 수행했으며, 추정값 ˆρ₀ 가 Jensen 부등식에 의해 결정론적 값보다 음의 편향을 보임을 확인했다. 이는 N·C_c 비율이 비선형 함수이기 때문에 평균값을 직접 대입하면 과대평가가 발생한다는 통계적 해석과 일치한다. 편향은 1/N₀² 수준으로 수렴하지만, 분산이 과소평가될 경우 편향이 더 크게 나타날 수 있다.

전체적으로 이 논문은 순환연료 원자로의 저인구 동역학을 다루는 최초의 확률점동학 프레임워크를 제시하고, 모델링 가정(완전 혼합, 전구체 노이즈 무시)의 영향을 정량적으로 평가한다. 향후 연구에서는 전구체 확산항을 재도입하고, 실제 혼합 정도를 반영한 다중 체적 모델을 구축함으로써 현재 모델의 한계를 보완할 필요가 있다.