건물 전체 온도를 온도계 하나로 잰다면 — 점동특성 방정식이 하는 일

큰 건물 하나를 상상해 보세요. 볕이 잘 드는 남쪽 사무실은 후덥지근하고, 그늘진 북쪽 창고는 서늘하고, 기계실 옆방은 늘 조금 더 뜨겁습니다. 층마다, 방마다 온도가 다 다르죠. 그런데 건물 관리실에는 온도계가 딱 하나 걸려 있고, 관리인은 그 숫자 하나만 보고 "지금 건물 온도가 오르고 있다" 또는 "떨어지고 있다"고 판단합니다. 방마다 정확히 몇 도인지는 모르지만, "건물 전체가 지금 뜨거워지는 중인가, 식는 중인가" 정도는 그 온도계 하나로 충분히 판단할 수 있는 거죠.

원자로 안에서도 똑같은 상황이 벌어집니다. 노심 안에는 위치마다 중성자가 얼마나 활발히 움직이는지를 나타내는 값(중성자속)이 다 다릅니다. 중심부는 활발하고 가장자리는 상대적으로 잠잠하죠. 이 복잡한 공간별 분포를 하나하나 정확히 계산하는 방법(앞서 다룬 확산이론이 바로 이 역할을 합니다)도 있지만, 원자로가 "지금 출력이 늘어나는 중인가, 줄어드는 중인가, 얼마나 빠르게?"라는 질문에 답하기 위해서라면 꼭 그 복잡한 공간 분포를 다 알 필요는 없습니다. 노심 전체를 대표하는 숫자 하나, 즉 "노심 전체 출력"이 시간에 따라 어떻게 움직이는지만 봐도 충분한 경우가 훨씬 많거든요. 오늘 이야기할 점동특성 방정식(point kinetics equations)이 바로 이 "대표 온도계 하나로 건물 전체를 판단하는" 방식의 원자로 버전입니다.

온도계 하나, 혹은 GDP 숫자 하나로 충분한 이유

"방마다 온도가 다 다른데, 온도계 하나로 뭉뚱그리는 게 말이 되나?" 싶으실 수도 있습니다. 그런데 우리는 일상에서 이런 식의 요약을 아주 자연스럽게 받아들입니다. 한 나라의 경제 상태를 이야기할 때도, 그 나라 안 수천만 가구의 소득과 지출을 낱낱이 따지는 대신 국내총생산(GDP)이라는 숫자 하나로 "경제가 성장하는 중인지, 위축되는 중인지"를 가늠하잖아요. GDP 하나로는 어느 동네가 잘살고 어느 동네가 어려운지는 알 수 없지만, "나라 경제가 전반적으로 지금 어느 방향으로 가고 있는가"라는 큰 흐름을 보는 데는 아주 요긴한 숫자입니다.

원자로의 점동특성 방정식도 정확히 이 발상입니다. 노심 구석구석의 중성자속 분포(어느 자리가 더 활발한지)라는 "동네별 사정"은 과감히 접어두고, 노심 전체를 대표하는 딱 하나의 값(전체 중성자 개체수, 또는 이에 비례하는 전체 출력)이 시간이 지나면서 어떻게 변하는지만 계산합니다. 그래서 이름도 "점(point) 동특성"입니다 — 노심이라는 넓은 공간을 계산상 하나의 점으로 뭉뚱그려 놓고, 그 점이 시간에 따라 오르내리는 모습만 쫓는다는 뜻이죠.

실제로 무엇을 뭉개고, 무엇을 남기는가

여기서 조금 더 정확히 짚어볼 지점이 있습니다. 원자로 안의 진짜 중성자속은 위치(공간)와 시간 두 가지에 동시에 의존하는 값입니다. 이걸 있는 그대로 풀어내려면 공간과 시간이 함께 얽힌 복잡한 방정식(공간-시간 동특성)을 매 순간 노심 전체에 대해 풀어야 하는데, 계산량이 어마어마합니다.

점동특성 방정식은 여기서 대담한 가정을 하나 겁니다. "중성자속의 공간적인 모양(어디가 상대적으로 더 활발하고 덜 활발한지의 패턴)은 시간이 지나도 거의 그대로 유지되고, 오직 그 모양 전체의 '크기'만 시간에 따라 커졌다 작아졌다 한다"는 가정입니다. 건물 비유로 치면, "남쪽 사무실이 북쪽 창고보다 늘 더 따뜻하다"는 상대적인 패턴 자체는 그대로 유지된 채, 건물 전체가 통째로 조금 더워졌다 식었다만 반복한다고 가정하는 셈이죠. 이 가정 덕분에 복잡한 공간-시간 방정식이, 시간에 대한 몇 개의 간단한 방정식으로 확 줄어듭니다.

점동특성 방정식, 실제로 이렇게 생겼습니다

말로 풀어낸 이야기를 수식 하나로 정리하면 이렇습니다 (선행핵을 대표 하나로 단순화한 형태).

dn/dt = [(ρ − β) / Λ] × n + λ × C
dC/dt = (β / Λ) × n − λ × C

기호가 낯설어 보이지만 하나씩 풀면 어렵지 않습니다.

• n: 노심 전체를 대표하는 중성자 개체수, 즉 이 방정식이 실제로 알고 싶어 하는 "노심 전체 출력"에 해당하는 값입니다. 앞서 말한 "대표 온도계 숫자"가 바로 이 n이에요.
• ρ(반응도): 원자로가 "정확히 균형 잡힌 상태(임계)"에서 얼마나 벗어나 있는지를 나타내는 값입니다. ρ가 0보다 크면 출력이 늘어나는 방향, 0보다 작으면 줄어드는 방향, 정확히 0이면 지금 이 순간 출력이 그대로 유지된다는 뜻입니다. 흔히 "반응도를 넣는다·뺀다"는 표현을 쓰는데, 이는 실제 출력 자체를 직접 바꾼다는 뜻이 아니라 "출력이 앞으로 변화할 방향과 속도를 조절한다"는 뜻에 가깝습니다.
• β(지발중성자 비율): 핵분열에서 나오는 중성자 대부분은 거의 즉시 튀어나오지만(즉발중성자), 극히 일부는 핵분열 부산물이 나중에 붕괴하면서 시간차를 두고 뒤늦게 튀어나옵니다(지발중성자). β는 전체 중성자 중 이 "뒤늦게 나오는 몫"의 비율로, 흔히 교재에서 예시로 드는 대표값은 1% 미만(우라늄-235 열중성자 핵분열 기준 약 0.65% 안팎) 수준입니다.
• Λ(중성자 세대시간): 한 세대의 중성자가 태어나 다음 세대의 핵분열을 일으키기까지 걸리는 아주 짧은 시간입니다. 열중성자로 기준으로 보통 십만분의 1초보다도 짧은 수준으로 소개됩니다.
• λ와 C: λ는 지발중성자를 뒤늦게 내놓는 "선행핵"이 붕괴하는 속도, C는 그 선행핵이 지금 노심 안에 얼마나 쌓여 있는지를 나타내는 값입니다. 앞으로 튀어나올 지발중성자의 "예비 재고"라고 생각하면 됩니다.

이 지발중성자 항(λ, C)이 왜 중요하냐면, 만약 이 항이 없고 즉발중성자만으로 원자로가 반응한다면 출력 변화 속도가 오직 아주 짧은 Λ에만 좌우되어 사람이 손쓸 틈도 없이 순식간에 변해버립니다. 지발중성자가 섞여 있는 덕분에 원자로의 실제 반응 속도는 사람과 제어계통이 충분히 따라잡을 수 있는 수준으로 늦춰집니다. 점동특성 방정식에서 이 항은 곁다리 보정이 아니라, 원자로가 통제 가능한 이유 그 자체를 담고 있는 핵심 부분입니다.

이 방정식이 계산해주는 것 — 시간에 따라 출력이 그리는 곡선

결국 점동특성 방정식을 풀면 무엇을 얻을까요? 바로 시간이 지남에 따라 노심 전체 출력(n)이 어떤 모양의 곡선을 그리며 오르내리는지, 그 궤적입니다. 반응도 ρ를 어떻게 넣어주느냐(제어봉을 얼마나 뽑는지, 온도가 어떻게 변하는지 등)에 따라 이 곡선의 기울기와 모양이 달라지고, 지발중성자 항 덕분에 이 곡선은 급격한 계단이 아니라 사람이 대응할 수 있는 정도의 매끄러운 변화로 나타납니다.

확산이론과는 무엇이 다른가 — 뭉개는 대상 자체가 다르다

앞서 다룬 확산이론 이야기와 헷갈리기 쉬운 지점이 있어 짚고 넘어가겠습니다. 확산이론은 "중성자 한 마리 한 마리의 무작위한 궤적을 다 추적하는 대신, 노심 공간 곳곳의 중성자속 분포를 매끄러운 함수로 계산해내는" 근사였습니다. 즉 확산이론은 "노심 어디에 중성자가 얼마나 있는가"라는 공간 정보를 여전히 정교하게 풀어냅니다. 다만 그 계산 방식을 개별 입자 추적 대신 연속적인 분포 계산으로 가볍게 만든 것뿐이죠.

점동특성은 이것과는 전혀 다른 축의 단순화입니다. 확산이론이 애써 계산해준 "노심 어디가 얼마나 활발한가"라는 공간 정보 자체를 아예 통째로 무시합니다(또는 그 모양이 고정되어 있다고 가정합니다). 대신 "노심 전체가 지금 얼마나 활발한가"라는 시간 축 하나에만 집중하죠. 비유하자면 확산이론이 "건물 안 모든 방의 온도 지도를 정밀하게 그려내는 일"이라면, 점동특성은 "그 지도는 신경 쓰지 않고, 관리실 온도계 숫자 하나가 시간에 따라 어떻게 변하는지만 기록하는 일"에 해당합니다. 같은 "복잡한 것을 근사로 단순화한다"는 공학적 태도를 공유하지만, 무엇을 포기하고 무엇을 남기는지가 서로 다른 두 가지 도구인 셈입니다.

그래서 어디에 쓰이는가 — 실용적 공학 도구로서의 쓸모와 한계

점동특성 방정식은 원자로 운전원 훈련용 시뮬레이터, 초기 단계의 안전해석, 제어계통 설계 검토처럼 "노심 전체 출력이 시간에 따라 대략 어떻게 반응하는가"를 빠르게 파악해야 하는 실무 곳곳에서 기본 도구로 쓰입니다. 노심 전체에 고르게 영향을 주는 변화(예: 냉각재 온도가 노심 전반에 걸쳐 서서히 변하는 경우)라면, 공간적인 모양이 크게 흐트러지지 않으므로 이 근사가 잘 들어맞습니다.

다만 한계도 분명합니다. 노심의 아주 좁은 한 자리에서만 반응도가 급격하게 튀는 것과 같이, 공간 분포 자체가 실제로 크게 일그러지는 상황에서는 "공간 모양이 고정되어 있다"는 점동특성의 전제 자체가 깨집니다. 이럴 때는 공간과 시간을 함께 고려하는 더 정교한 3차원 동특성 해석으로 넘어가야 합니다. 즉 점동특성 방정식은 모든 상황에 만능으로 쓰이는 도구가 아니라, "노심 전체가 거의 같은 리듬으로 함께 움직이는" 많은 실용적 상황에서 계산 부담을 극적으로 줄여주는 출발점입니다.

자주 헷갈리는 지점들

몇 가지 오해를 짚어보겠습니다. 먼저 "점(point)"이라는 이름 때문에 원자로가 실제로 점처럼 작다는 뜻으로 오해하기 쉬운데, 전혀 아닙니다. 계산 모델 안에서 노심 전체를 대표하는 하나의 값으로 뭉뚱그렸다는 수학적 근사를 가리키는 이름일 뿐, 원자로의 실제 물리적 크기와는 아무 상관이 없습니다.

또, 반응도(ρ)를 원자로의 출력 그 자체로 착각하기도 쉽습니다. 반응도는 "지금 균형 상태로부터 얼마나 벗어나 있는가"를 나타내는 값이지, 출력의 절대적인 크기(몇 메가와트인지)를 뜻하지 않습니다. 반응도가 커졌다고 출력이 이미 높다는 뜻이 아니라, 출력이 앞으로 늘어나는 방향으로 움직이고 있다는 뜻입니다.

마지막으로 점동특성 방정식을 확산이론의 "시간 버전"쯤으로 여기기 쉬운데, 둘은 서로의 연장선이 아니라 애초에 무엇을 뭉갤지 자체가 다른 별개의 근사입니다. 확산이론은 공간 정보를 포기하지 않는 도구이고, 점동특성은 그 공간 정보 자체를 통째로 접어두는 도구입니다. 둘은 경쟁 관계가 아니라, 노심 설계처럼 공간 분포가 중요한 질문에는 확산이론(또는 더 정밀한 수송이론)을, 노심 전체 출력의 시간 변화가 궁금한 질문에는 점동특성을 각각 골라 쓰는, 서로 다른 목적의 도구입니다.

오늘의 한 줄 정리

• 점동특성 방정식은 노심 안 위치별로 다른 복잡한 중성자속 분포(공간 의존성)를 무시하고, 원자로 전체를 대표하는 하나의 값(전체 출력)이 시간에 따라 어떻게 변하는지만 계산하는 근사적 공학 도구다.
• 건물 전체 온도를 온도계 하나로 판단하거나, 나라 경제를 GDP 하나로 요약하는 것과 같은 발상이다.
• 확산이론이 "공간 분포 자체를 정교하게 계산하는" 도구라면, 점동특성은 "그 공간 분포를 아예 통째로 뭉개고 시간 변화만 보는" 도구로, 서로 다른 축의 단순화다.
• 반응도(ρ), 지발중성자 비율(β), 중성자 세대시간(Λ), 선행핵 농도(C) 같은 몇 개의 값만으로 노심 전체 출력이 시간에 따라 그리는 곡선을 계산할 수 있으며, 특히 지발중성자 항 덕분에 그 변화가 사람이 대응 가능한 속도로 느려진다.
• 다만 노심 안 특정 자리만 국소적으로 크게 요동치는 상황에서는 이 근사의 전제가 깨져, 공간까지 함께 고려하는 더 정교한 해석이 필요해진다.

노심 안에서 중성자가 "어디에 얼마나" 있는지를 살펴봤던 지금까지의 이야기에 이어, 이제 원자로가 시간에 따라 실제로 어떻게 반응하고 움직이는지를 들여다보는 여정이 이렇게 시작됩니다. 반응도라는 조작 언어, 그리고 그 반응도가 노심 전체 출력을 어떻게 밀고 당기는지에 관한 이야기가 앞으로도 여러 각도에서 계속 이어질 예정입니다.