🔬 physics

Assessment of the neutron fluence during operation of CEMHTI Cyclotron using Monte-Carlo codes - Simβ-AD study

Simβ-AD 연구는 활성화 포일의 중성자 조사량 측정값을 4가지 몬테카를로 코드(FLUKA, MCNP6, PHITS, RayXpert®)와 검증함으로써 CEMHTI 사이클로트론의 방사성 폐기물 활성화를 평가하는 방법론을 확립하였으며, 단순화된 모델을 사용한 실험 데이터와 시뮬레이션 간의 양호한 일치를 입증하였다.

원저자: Jean-Michel HORODYNSKI, Abir HASSANI, Frédéric CHAPELLE, Nicolas ARBOR, Djokhar BETELGUERIEV, Lucia Victoria GARCIA GARCIA, Stéphane HIGUERET, Léo NOWAK, Cédric DOSSAT, Sébastien BOUILLON, David CHAUL

게시일 2026-07-07

📖 3 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Jean-Michel HORODYNSKI, Abir HASSANI, Frédéric CHAPELLE, Nicolas ARBOR, Djokhar BETELGUERIEV, Lucia Victoria GARCIA GARCIA, Stéphane HIGUERET, Léo NOWAK, Cédric DOSSAT, Sébastien BOUILLON, David CHAULIN, William HATE, Patrice RIFARD, Thierry SAUVAGE

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

핵심 요약: 남겨진 "유령"을 예측하기

매우 강력한 기계(사이클로트론)가 양성자와 중수와 같은 미세 입자들을 타겟에 쏘고 있다고 상상해 보세요. 이 입자들이 충돌하면 단순히 멈추는 것이 아니라, 중성자라고 불리는 보이지 않는 "유령"들의 혼란스러운 샤워를 만들어냅니다. 이 유령들은 방 안을 돌아다니며 벽, 자석, 장비에 부딪히고, 이들을 약간 방사능을 띠게 만듭니다.

이 연구의 주요 목표는 다음과 같습니다: 우리의 컴퓨터 시뮬레이션이 이 "유령"들이 방을 얼마나 "방사능화"할지 정확하게 예측할 수 있는가?

연구진은 나중에 모든 것을 다 해체하지 않고도 이 방사능을 측정할 수 있는 신뢰할 수 있는 레시피(방법론)를 구축하고자 했습니다. 이는 방사성 폐기물을 안전하게 관리하는 데 매우 중요합니다.

등장인물 소개

사이클로트론(The Cyclotron): 프랑스 CEMHTI 시설에 있는 거대한 입자 가속기입니다. 미세 입자들을 위한 고속 경주 트랙과 같습니다. 최근에 운영이 종료되었는데, 이는 트랙이 영원히 사라지기 전에 아이디어를 테스트하기에 완벽한 장소가 되었습니다.
"유령"(중성자): 빔이 타겟에 부딪힐 때 생성되는 보이지 않는 입자들입니다. 이들은 다른 물질을 방사능을 띠게 만드는 주범입니다.
"끈끈이 트랩"(활성화 박막): 연구진은 방 안의 여러 지점에 금, 스칸듐, 탄탈륨, 테르븀으로 된 작고 얇은 금속 판들을 배치했습니다. 이것을 끈끈이 파리 잡는 종이라고 생각하면 됩니다. 중성자 유령들이 날아오면 이 금속에 달라붙어 금속의 상태를 미세하게 변화시킵니다. 나중에 연구진은 각 판에 얼마나 많은 "끈적임"(방사능)이 남아 있는지 측정했습니다.
"수정 구슬"(몬테카를로 코드): 이들은 네 가지 서로 다른 슈퍼컴퓨터 프로그램(FLUKA, MCNP6, PHITS, RayXpert)입니다. 이들은 일기 예보관과 같습니다. 방의 지도와 입자의 속도를 입력하면, 이 프로그램들은 얼마나 많은 유령이 끈끈이 트랩에 부딪힐지 정확히 예측하려고 노력합니다.

실험: 현실 점검

연구팀은 두 종류의 빔(양성자와 중수)으로 사이클로트론을 가동했으며, 빔 출구 근$처, 자석 근처, 그리고 벽면에 "끈끈이 트랩"을 배치했습니다.

기계가 꺼진 후, 그들은 금속 판들을 실험실로 가져가 얼마나 방사능을 띠게 되었는지 정확히 측정했습니다. 그런 다음, 이 실제 세계의 수치들을 네 가지 컴퓨터 프로그램의 예측값과 비교했습니다.

결과: 수정 구슬은 얼마나 잘 맞혔나?

연구진은 컴퓨터 프로그램들이 결과를 꽤 잘 맞혔지만, 완벽하지는 않다는 것을 발견했습니다.

점수: 컴퓨터 예측값은 보통 실제 측정값의 50%에서 150% 범위 안에 있었습니다. 복잡한 물리학의 세계에서 이는 "좋은 일치(good agreement)"로 간주됩니다. 이는 일기 예보관이 "비가 1~3인치 정도 올 것입니다"라고 말했을 때, 실제로 2인치의 비가 내린 것과 같습니다.
"금(Gold)"의 오류: 한 프로그램(FLUKA)은 금 판에 대해 특정 문제를 보였습니다. 금 원자의 특정 "상태"(예: 잠자는 고양이를 깨어 있는 고양이로 착각하는 것과 같은 상황)를 처리하는 데 어려움을 겪어 방사능 양을 잘못 예측했습니다. 이 논리 구조를 수정하자 예측 성능이 훨씬 좋아졌습니다.
"구석" 문제: 빔 출구 근처(입자들이 처음 부딪히는 곳)에서는 예측값이 다소 흩어져 있었습니다. 이는 연못에 돌을 던졌을 때 물이 튀는 모습을 예측하려는 것과 같습니다. 물보라가 일어나는 근처의 파동은 매우 혼란스럽고 완벽하게 모델링하기 어렵습니다. 컴퓨터들은 이 좁은 구석에서 저에너지 "유령"들을 때때로 과소평가했습니다.

이 연구가 중요한 이유 (논문에 따르면)

이 연구는 컴퓨터 코드를 사용하여 사이클로트론 시설이 얼마나 많은 방사성 폐기물을 생성할지 추정할 수 있음을 증명합니다.

"베타" 문제: 일부 방사성 물질은 "베타"라고 불리는 유형의 에너지만 방출하는데, 이는 물체를 파괴하지 않고는 감지하기 어렵습니다. 논문은 우리가 이 코드들을 통해 중성자 활성도를 정확하게 예측할 수 있다면, 기계를 절단하지 않고도 "베타 전용" 폐기물을 추정할 수 있다고 제안합니다.
핵심 결론: 컴퓨터 모델과 몇 가지 실제 측정값(끈끈이 트랩)을 결합함으로써, 입자 가속기로부터 발생하는 방사성 폐기물을 관리할 수 있는 신뢰할 수 있는 방법을 만들 수 있습니다.

한 줄 요약

연구진은 입자 가속기가 가동된 후 방이 얼마나 방사능을 띠게 될지 예측하기 위해 네 가지 서로 다른 컴퓨터 프로그램을 테스트했습니다. 그들은 보이지 않는 입자를 잡기 위해 실제 금속 판을 "센서"로 사용했습니다. 컴퓨터는 준수한 성적(1.5배 이내의 오차)을 거두었으며, 이는 기계가 어떻게 작동하는지 정확히 알고 소프트웨어의 작은 논리 오류를 수정한다면, 디지털 시뮬레이션을 통해 방사성 폐기물 관리를 계획할 수 있음을 입증했습니다.

기술 요약: Monte-Carlo 코드를 이용한 CEMHTI 사이클로트론 운전 중 중성자 조사량 평가

문제 정의
사이클로트론 시설의 효과적인 방사성 폐기물 관리를 위해서는 정확한 방사선 특성 분석이 필요하며, 특히 파괴적 측정법을 필요로 하는 $\beta$ 단독 핵종(예: $^{36}$ Cl, $^{49}$ V, $^{55}$ Fe, $^{63}$ Ni)에 대한 정밀한 분석이 요구된다. Sim $\beta$ -AD 연구는 Monte-Carlo(MC) 코드와 능동 중성자 검출기를 결부터 결합하여 사이클로트론 시설 내 재료 활성화를 평가하는 방법론을 구축하기 위해 시작되었다. 본 논문에서 다루는 구체적인 과제는 MC 코드를 실험 데이터와 벤치마킹하여 중성자 조사량 계산을 검증하고, $\beta$ 단독 핵종과 $\gamma$ 방출 핵종 사이의 비율을 정의하는 것이다. 본 연구는 1976년부터 2023년 가동 중단 전까지 운영된 CEMHTI 사이클로트론 시설에 초점을 맞추고 있으며, 이는 시설 해체를 앞두고 실제 운영 데이터를 통해 모델을 검증할 수 있는 독특한 기회를 제공한다.

방법론
본 연구는 네 가지 서로 다른 Monte-Carlo 코드인 FLUKA, MCNP6, PHITS, RayXpert를 사용하여 실험적 측정값과 수치 시뮬레이션 간의 비교 접근법을 채택하였다.

기하학적 및 물리적 모델링:
- 설계도와 현장 측정을 기반으로 CEMHTI 벙커, 실험 구역, 사이클로트론 구성 요소(자석, 요크, 인플렉터) 및 타겟에 대한 고충실도(high-fidelity) 모델을 구축하였다.
- 두 가지 주요 빔 시나리오를 모델링하였다:
  1. 양성자 빔: 16 MeV, 25 µA가 탄탈룸(Tantalum) 블록에 정지하며, 디플렉터(deflector)에서 2차 손실 발생.
  2. 중수소 빔: 17.5 MeV, 30 µA가 베릴륨(Beryllium) 타겟에 정지.
- 중성자 생성 통계(statistics)를 향상시키기 위해 특정 단면적 라이브러리(JEFF 3.3, TENDL-2019/2021/2024, ENDF/B-VIII.0) 사용 및 분산 감소 기법(예: 강제 충돌, 중요도 샘플링)을 포함한 특정 물리 설정을 각 코드에 최적화하였다.
실험 설정:
- 2021년 두 차례의 조사 캠페인 동안 추출부, 파이프, 요크, 디플렉터, 벽, 문, 자석 등 다양한 위치에 **금(Au), 스칸듐(Sc), 탄탈룸(Ta), 테르븀(Tb)**으로 제작된 수동 활성화 포일(passive activation foils)을 배치하였다.
- 포일은 양성자의 경우 5시간, 중수소의 경우 2시간 동안 조사되었다.
- 조사 후, 40% 효율의 ISOCS 검출기를 사용하여 $\gamma$ 방출 핵종( $^{198}$ Au, $^{46}$ Sc, $^{182}$ Ta, $^{160}$ Tb)의 방사능을 정량화하기 위한 $\gamma$ -분광 분석을 수행하였다.
활성화 계산:
- MC 코드로부터 도출된 중성자 조사량은 FISPACT-II(TENDL-2019 라이브러리 포함) 및 RayActive 코드의 이전 버전의 오프라인 활성화 계산을 위한 입력값으로 사용되었다.
- FLUKA는 또한 RESNUCLEI 카드를 사용하여 온라인 활성화 계산을 수행하였다.

주요 결과
본 연구는 측정된 방사능과 계산된 값을 비교하였으며, 다음과 같은 결과를 얻었다.

전반적인 일치도: 네 가지 코드 모두에서 측정값과 계산값 사이에 양호한 일치를 보였으며, 측정 대비 계산된 방사능의 비율은 0.5–1.5 범위 내에 있었다.
양성자 빔:
- FLUKA 온라인 계산은 금(Au) 포일을 제외하고는 양호한 일치를 보였으나, 이는 RESNCLEI 카드에서 임의로 설정된 50/50 기저 상태/이성질체 비율로 인한 불일치 때문이었으며, 이는 FISPACT-II를 사용하여 해결되었다.
- 빔 손실 지점(디플렉터, 파이프) 근처에 위치한 포일, 특히 저에너지 중성자에 민감한 Au, Ta, Tb의 경우 높은 불일치를 보였다. 이는 빔에 대한 작은 각도에서의 저에너지 중성자 조사량 계산의 차이를 시사한다.
- RayXpert와 MCNP6는 유사한 경향을 보였으며, PHITS는 대부분의 경우 FLUKA의 경향을 따랐다.
중수소 빔:
- FLUKA, MCNP6, PHITS의 조사량을 결합한 FISPACT-II 계산은 양호한 일치를 보였다.
- FLUKA 온라인 계산은 이성질체 평형 문제로 인해 Au 포일의 방사능을 크게 과소평가하였다.
- 다른 코드들은 Be(d,n) 반응으로부터의 저에너지 중성자 과대평가 및 그에 따른 콘크리트 벽에서의 열화(thermalization)로 인해 활성화를 과대평가하는 경향이 있었다.
코드 성능: RayXpert와 MCNP6는 빔 손실 지점 근처에서 다른 코드들에 비해 더 높은 불일치를 보였으며, 이는 작은 각도에서의 저에너지 중성자 조사량 과소평가에 기인한다.

의의 및 주장
본 논문은 적절한 시설 모델과 결합된 Monte-Carlo 코드를 이용한 중성자 조사량 추정이 사이클로트론 환경 내 수동 검출기의 활성화를 평가하는 신뢰할 수 있는 방법을 제공한다고 주장한다. 본 연구는 다음을 입증한다:

방법론 검증: MC 코드와 오프라인 활성화 도구(FISPACT-II, RayActive)의 결합은 허용 가능한 성능 범위(0.5–1.5 비율) 내에서 방사성 폐기물 활성화를 특성화하는 유효한 방법론임을 보여준다.
불확실성 원인: 양성자 빔 시나리오의 주요 불확실성 원인은 가속기 내부의 빔 손실 강도 및 형상에 대한 정밀한 데이터 부족으로 식별되었다.
향후 적용: 현재 가동이 중단된 CEMHTI 사이클로트론의 특성을 고려할 때, 검증된 방법론은 해체된 사이클로트론 부품 및 콘크리트 벽에서 측정된 값과 수치 시뮬레이션을 비교하는 데 적용될 것이다. 이는 MC 코드, 단면적 라이브러리 및 중성자 열화 거동 간의 구체적인 차이에 대한 추가 연구를 촉진하여, 궁극적으로 방사성 폐기물 관리를 위한 $\beta$ 단독 대 $\gamma$ 방출 핵종 비율을 정의하려는 Sim $\beta$ -AD의 목표를 지원할 것이다.

저자들은 빔 파라미터에 대한 지식 개선이 필요하지만, 현재의 결과는 가속기 시설의 방사선 특성 분석을 위한 이러한 시뮬레이션 도구의 유용성을 확인해 준다고 결론짓는다.