آشنایی با چشمه‌های بحرانی و کارت KCODE در MCNP

در بخش‌های قبلی آموختیم که چگونه با کارت SDEF یک چشمه مشخص (مانند یک چشمه کبالت یا منبع نوترونی) را تعریف کنیم. اما در طراحی قلب یک رآکتور هسته‌ای یا بررسی ایمنی یک انبار سوخت، ما با “چشمه‌ای” سرکار داریم که خودش تابع چیدمان مواد است. در این حالت، هدف ما پیدا کردن ضریب تکثیر موثر ($k_{eff}$​) است. برای این منظور، MCNP از ساختار متفاوتی به نام KCODE استفاده می‌کند.

تفاوت مبنایی چشمه ثابت و محاسبات بحرانیت

در یک مسئله چشمه ثابت، ذرات از یک نقطه یا سطح مشخص گسیل می‌شوند و سفر آن‌ها با جذب یا فرار تمام می‌شود. اما در محاسبات بحرانیت (Criticality)، شبیه‌سازی به صورت “چرخه‌ای” (Cycle) انجام می‌شود. نوترون‌های حاصل از شکافت در یک نسل، چشمه نسل بعدی را می‌سازند. هدف این است که بدانیم آیا جمعیت نوترون‌ها در هر نسل ثابت می‌ماند ($k_{eff}=1$)، رشد می‌کند ($k_{eff}>1$) یا کاهش می‌یابد ($k_{eff}<1$).

کالبدشکافی کارت KCODE: فرمانده محاسبات رآکتور

کارت KCODE ابزار اصلی برای تعریف پارامترهای این فرآیند چرخه‌ای است. این کارت معمولاً دارای چهار پارامتر اصلی است:

1. NSRCK: تعداد نوترون‌ها در هر نسل (مثلاً ۵۰۰۰ ذره).
2. RKK: حدس اولیه برای مقدار $k_{eff}$​ (معمولاً عدد ۱.۰).
3. IKZ: تعداد چرخه‌های “صرف‌نظر شده” یا Settling Cycles.
4. KCT: تعداد کل چرخه‌هایی که شبیه‌سازی خواهد شد.

فرمول پایه که کد سعی در حل آن دارد، نسبت جمعیت نوترون‌ها در دو نسل متوالی است:

$k_{eff}=\frac{N_{n+1}}{N_n}$

کارت KSRC: محل تولد اولین نسل

در یک رآکتور، ما دقیقاً نمی‌دانیم نوترون‌ها در کجا متولد می‌شوند (این همان چیزی است که کد باید پیدا کند). اما برای شروع بازی، نیاز به یک “جرقه” اولیه داریم. کارت KSRC (Kernel Source) مختصات نقاطی را دریافت می‌کند که اولین نسل نوترون‌ها باید از آنجا سفر خود را آغاز کنند. این نقاط باید حتماً درون مواد شکافت‌پذیر (مانند میله‌های سوخت) قرار داشته باشند. پس از چند چرخه اولیه، توزیع چشمه به حالت تعادل فیزیکی رسیده و دیگر به نقاط اولیه کارت KSRC وابسته نخواهد بود.

اهمیت چرخه‌های صرف‌نظر شده (Inactive Cycles)

یکی از نکات کلیدی در کاهش واریانس و افزایش دقت در محاسبات بحرانیت، تعیین صحیح پارامتر IKZ است. در چرخه‌های ابتدایی، توزیع نوترون‌ها هنوز با هندسه واقعی قلب رآکتور تطبیق نیافته است. اگر این چرخه‌ها را در میانگین‌گیری نهایی دخالت دهیم، نتیجه $k_{eff}$​ دچار خطا می‌شود. MCNP این چرخه‌ها را اجرا می‌کند تا “توزیع چشمه” شکل بگیرد، اما نتایج آن‌ها را در آمار نهایی لحاظ نمی‌کند. این فرآیند باعث می‌شود که شبیه‌سازی از نظر آماری همگرا (Converged) شود.

کاربرد در تحلیل ایمنی و طراحی قلب

استفاده از کارت KCODE تنها به طراحی رآکتور محدود نمی‌شود. در مباحث “ایمنی بحرانیت” (Criticality Safety)، برای اطمینان از اینکه مخازن نگهداری سوخت یا پسماندهای هسته‌ای تحت هیچ شرایطی به حالت بحرانی نمی‌رسند، از این کارت استفاده می‌شود. دقت در تعریف صحیح ابعاد و چگالی مواد در کنار تنظیمات دقیق KCODE، تضمین‌کننده امنیت یک تأسیسات هسته‌ای است.