23. Solving the Neutron Diffusion Equation, and Criticality Relations

MIT에서의 방사선 문제 해결

수업 개요

• 화요일에 MIT에서 사용할 가장 큰 방정식을 개발하고, 목요일에는 이를 간단한 형태로 줄였으며, 오늘은 다양한 원자로 문제를 해결하는 방법을 보여줄 것임.

• AP 1000 원자로와 같은 두 그룹 원자로 문제를 다루며, 중성자를 빠른 그룹과 열 그룹으로 나누는 방식에 대해 설명할 예정임.

방정식 및 그래프

• 중성자의 수를 설명하는 방정식이 있으며, 이 방정식을 통해 플럭스(Flux) 대 x의 그래프 형태를 이해하려고 함.

• 오늘 수업에서는 이 방정식을 실제로 풀어보며 어떻게 사용하는지에 대한 실습을 진행할 것임.

단순화 과정

• NIN 반응과 포토 분열을 생략하여 더 간단한 형태의 방정식으로 시작할 것임.

• 라플라시안 연산자는 차원과 좌표계에 따라 다른 형태를 가지며, 1차원 카르테시안에서는 간단하지만 원통형에서는 복잡함.

무한 슬래브 케이스

• 무한 슬래브 케이스부터 시작하여 해석적으로 쉽게 풀 수 있는 방법을 모색할 것임.

• 라플라시안을 분리하고 플럭스 단위로 나누어 간단한 형태로 정리함.

평균화된 단면적

• 평균화된 단면적을 사용하여 에너지 범위 전체에서의 교차 섹션을 평균내는 과정을 설명함.

• 카르테시안 공간에서 두 번째 미분이 자기 자신과 같아지는 함수가 무엇인지 질문하며 지수 함수와 사인/코사인 함수에 대해 논의함.

플럭스 프로파일 도출

• 플럭스 솔루션이 특정 형식을 가져야 한다고 가정하며, 대칭성을 고려해 사인 항목은 제거됨.

• 플럭스가 원자로 가장자리에 도달하지 않도록 설정하며, 이는 중성자가 계속해서 흐르고 있음을 나타냄.

확장 길이 및 확산 상수

각종 물리적 성질과 재료의 상관관계

코사인 함수와 원자 질량

• 각도 코사인은 대략적으로 산란되는 물체의 평균 원자 질량의 2/3배에 해당한다. 이 d는 물리적 비유에서 도입된 것으로, 교차 섹션 및 재료 특성에서 유래한 표현이 있다.

재료와 기하학적 조건

• 모든 정보는 Janis 라이브러리 또는 사용 중인 교차 섹션 라이브러리에서 확인할 수 있으며, 실제로 사용하는 숫자 밀도와 결합된다. 이제 상황이 매우 현실적이고 물리적으로 다가오고 있다.

방정식 변형 및 기하학적 요소

• 코사인 bx를 대입하여 방정식을 다시 작성하면, 최종적으로 b² = a의 여러 재료 특성과 같은 간단한 형태로 정리된다. 여기서 기하학에 대한 정보는 포함되지 않는다.

k 효과적인 값 계산

• 만약 반응기의 기하학과 재료를 알면, 최종 미지수인 k 효과적인 값을 해결할 수 있는 조건을 설정했다. bg가 무엇인지 알아보아야 한다.

임계 상태 이해하기

• k(임계성)는 중성자의 이득과 손실 비율이다. fission에 의한 유일한 이득 메커니즘과 흡수 또는 누출에 의한 손실 메커니즘이 존재한다. 이 조건이 만족되면 반응기는 임계 상태에 있다고 할 수 있다.

온도 변화가 k 효과적인 값에 미치는 영향

흡수 증가 시나리오

• 흡수를 증가시키면 k 효과적인 값은 감소해야 한다. 이는 분모가 커지기 때문에 비율이 작아진다.

온도 상승 시 결과 분석

• 대부분의 교차 섹션은 온도가 상승하면 감소하는 경향이 있다. 이는 Doppler broadening 현상 때문이며, 온도가 오르면 밀도가 줄어들기 때문이다.

확산 상수 변화 예측

반응로의 효과적인 k 값은 어떻게 변할까?

온도 피드백의 개념

• 반응로에서 k effective는 재료의 상대적 양에 따라 달라지며, 이는 긍정적 또는 부정적 온도 피드백을 초래할 수 있다.

• 온도가 상승하면 k effective가 하락해야 하지만, 특정 재료를 사용할 경우 온도가 상승하면서 k effective가 증가하는 상황이 발생할 수 있으며, 이는 설계상 금지된다.

반응로 크기 변화의 영향

• 반응로의 크기를 늘리면 어떤 요소들이 증가하거나 감소하는지를 탐구한다.

• 더 많은 반응로를 추가하면 k effective는 증가하며, 이는 직관적으로 이해할 수 있는 패턴이다.

전력과 비판성 상태

• 전력은 중성자의 운동 에너지에서 나오지만, 비판성과는 반드시 연결되지 않는다.

• 비판적인 상태에서도 전력이 매우 적거나 아예 없는 경우가 있을 수 있으며, 이러한 시스템은 중성자 물리학을 테스트하는 데 유용하다.

동적 평형과 제어봉 조작

• k effective가 1이 아닌 경우 반응로는 동적 평형 상태에 있지 않으며, 이때 전력이 변화하기 시작한다.

• 제어봉을 철회했을 때 전력이 즉시 상승하고 이후에는 느려지는 현상이 나타나며, 이는 지연 중성자 때문이라는 점이 강조된다.

두 그룹 균형 방정식으로 일반화하기

• 단일 그룹 균형 방정식을 통해 얻은 통찰력을 바탕으로 두 그룹 균형 방정식으로 일반화하려고 한다.

시그마 열 및 빠른 중성자

중성자의 기여

• 시그마 열과 시그마 빠른 중성자가 모두 핵분열에 기여할 수 있음.

• 모든 변수를 빠른 부분과 열 부분으로 나누어 표현해야 함.

• 1 eV 이하에서 태어나는 중성자는 없다고 가정함.

손실 메커니즘

• 빠른 그룹의 중성자는 흡수로 인해 떠날 수 있음.

• 스캐터링이 또 다른 손실 메커니즘임을 언급함.

스캐터링과 누출

스캐터링 비율

• 모든 스캐터링 사건이 빠른 그룹을 떠나게 하지는 않음.

• 수소의 경우, 스캐터링 확률이 일정하다고 설명함.

누출

• 중성자가 원자로를 떠날 수 있는 누출 현상도 존재함.

열 그룹의 손익

열 그룹으로의 진입

• 열 에너지 그룹으로 들어오는 유일한 출처는 빠른 그룹에서의 스캐터링임.

• 빠른 그룹에서 나가는 모든 중성자는 열 그룹으로 들어감.

손실 메커니즘

• 열 그룹에서도 누출이 발생할 수 있으며, 이는 별도의 확산 계수로 표현됨.

• 흡수 또한 중요한 손실 메커니즘임.

k 효과와 임계 상태

k 효과 정의

• k 효과는 원자로의 기하학적 구조와 재료에 따라 결정됨.

중성자 생산과 손실의 이해

중성자의 효과적인 수와 손실

• 교수는 k effective가 중성자의 총 원천 아래에 위치한다고 설명하며, 이는 모든 중성자의 생산을 설명하고 나머지는 손실이라고 언급함.

• 실험적으로 중성자가 1에서 10 MeV 사이에서 생성된다는 것을 알고 있으며, 문제를 확대하여 사고 실험을 제안함.

열적 및 빠른 중성자 스펙트럼

• 청중이 "스펙트럼?"이라고 질문하자, 교수는 열적 및 빠른 중성자를 구분하는 Chi 변수를 추가해야 한다고 강조함.

• 교수는 일부 중성자가 열적으로 태어날 경우의 모델링 방법을 설명하며, 실제로는 드물지만 시험 문제로 적합하다고 언급함.

방정식 정리 및 플럭스 대체

• 두 개의 방정식과 세 개의 미지수가 있는 상황에서 하나의 플럭스로 모든 항을 표현할 수 있다고 설명함.

• 상단 방정식에서 phi thermal을 다른 변수들로 고립시키기 위해 재배열하는 과정을 보여줌.

최종 결과 도출

• phi thermal에 대한 표현식을 얻은 후 이를 상단 방정식에 대입하여 최종 답변에 가까워짐.

• 모든 항이 빠른 플럭스에 대한 표현으로 남아있으며, 이를 통해 k effective를 물질 특성과 기하학으로만 정의할 수 있음을 강조함.

비율과 설계 조건

• k effective의 비율이 손실과 이득 간의 관계를 나타내며, 이는 핵 반응기의 설계와 관련됨.

• 교수는 이득과 손실 비율이 어떻게 작용하는지를 다시 한번 강조하며 질문을 받음.

질문 및 응답 시간

• 청중이 빠른 방정식에서 산란 항목이 누락되었는지 질문하자, 교수가 이를 확인하고 수정함.

수업 내용 요약

중성자 수송 방정식의 이해

• 교수는 중성자 수송 방정식에 대해 설명하며, 학생들이 이 방정식을 기억하고 그 의미를 설명할 수 있기를 원한다고 언급함.

• 학생들은 중성자 수송 방정식의 용어와 데이터 획득 방법, 그리고 간소화 단계의 정당성을 설명할 수 있어야 함.

• 교수는 물리적 이유로 확산 근사가 적용되는 지역과 그 한계를 질문하며, 학생들은 제어봉이나 연료 근처에서 확산 근사가 부적합하다고 답함.

• 확산은 장거리 안정 상태 솔루션을 설명하지만, 급격한 변화가 있는 곳에서는 실패한다는 점을 강조함.

• 교수는 중성자 확산 방정식을 단순화하여 비판 조건을 도출하는 문제를 제시할 것이라고 말함.

엔지니어링과 물질 특성

• 교수는 이 강좌가 핵공학 및 이온화 방사선 입문이라는 점을 상기시키며, 공학적인 측면이 중요하다고 강조함.

• 학생들은 다양한 기하학적 상황에서 중성자 수송 방정식이 어떻게 적용될 수 있는지를 이해해야 함.

• 라플라시안 연산자를 사용하여 1차원 무한 반응기의 경우를 다루고 있으며, 원통 좌표계에서도 유사하게 적용 가능하다고 설명함.

• 원통형 기하학에서 발생하는 베셀 함수에 대해 언급하며, 이러한 함수들이 사인 및 코사인과 유사하게 행동한다고 설명함.

• 효과적인 k 값에 대한 직관적인 이해가 필요하며, 평형 상태에서 벗어나는 물리적 상황에 대한 질문도 제기됨.

그래프와 플럭스

• 원통형 그래프에서 플럭스의 형태에 대해 질문이 이어지고, 대칭성과 관련된 내용을 교수는 설명함.

• r 및 z 방향으로 플럭스 형태가 첫 번째 베셀 함수로 나타날 것이라고 언급하며 다차원 반응기의 해법을 간단히 소개함.

• r과 z 부분이 분리 가능한 경우 해결이 더 쉬워진다는 점을 강조하고 Cartesian 좌표계만 다룰 것이라고 말함.