🔬 physics

An elementary method to determine the critical mass of a sphere of fissile material based on a separation of neutron transport and nuclear reaction processes

تقدم هذه الورقة طريقة تعليمية مبسطة بحتة لحساب الكتلة الحرجة لكرية انشطارية من خلال فصل عمليات التفاعل النووي وعمليات النقل عبر حجج إحصائية حول أطوال مسارات النيوترونات، مما يعطي نتائج ضمن بضعة بالمائة من القيم المعروفة ومن نهج معادلة الانتشار دون الحاجة إلى معادلات تفاضلية معقدة.

المؤلفون الأصليون: Steven K. Lamoreaux

نُشر 2026-03-03

📖 5 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Steven K. Lamoreaux

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

تخيل أنك تحاول الحفاظ على اشتعال نار مخيم. لديك كومة من الخشب (الوقود)، وتحتاج إلى كمية محددة تماماً من الخشب وبالشكل الصحيح حتى لا تنطفئ النار، ولكن أيضاً حتى لا تنفجر بشكل خارج عن السيطرة.

هذه الورقة البحثية لعالم الفيزياء س. ك. لاموريو (S.K. Lamoreaux) تدور حول إيجاد منطقة "غولديلوكس" (المنطقة المثالية) للوقود النووي — وتحديداً، تحديد الحد الأدنى للحجم الذي يجب أن تكون عليه كرة من اليورانيوم أو البلوتونيوم لضمان استمرار التفاعل المتسلسل النووي. وهذا ما يسمى الكتلة الحرجة.

عادةً، يتطلب حساب ذلك حواسيب فائقة ورياضيات معقدة تبدو كأنها شفرات فضائية. أراد لاموريو أن يوضح أنه يمكنك فهم الفيزياء الجوهرية باستخدام المنطق البسيط، وحساب التفاضل والتكامل الخاص بالمرحلة الثانوية، وبعض التشبيهات الذكية.

إليك تفصيل طريقته باستخدام لغة يومية بسيطة:

1. الفكرة الكبرى: شد الحبل

تخيل النيوترونات (الجسيمات الصغيرة التي تبدأ التفاعل) كعدائين في سباق داخل قلعة ضخمة ومنفوخة (كرة الوقود).

لكي يظل التفاعل النووي مستقراً (حرجاً)، يجب أن يظل عدد العدائين داخل القلعة ثابتاً بمرور الوقت.

الهدف: في كل مرة يتم فيها "امتصاص" عداء بواسطة قطعة خشب (مما يسبب انشطاراً)، يجب أن يولد عدداً كافياً من العدائين الجدد ليعوض نفسه بالإضافة إلى أولئك الذين تم امتصاصهم.
المشكلة: العداءون يستمرون في السقوط خارج القلعة عبر الجدران (التسرب).
الحل: يجب أن تكون القلعة كبيرة بما يكفي لكي يرتد العداءون داخلها لفترة طويلة قبل أن يجدوا المخرج. هذا الارتداد يمنحهم فرصاً أكثر للاصطدام بقطعة خشب وإنتاج عدائين جدد.

2. الرقصة ذات الخطوتين

يقسم لاموريو المشكلة إلى رقصتين منفصلتين لا تتواصلان مع بعضهما البعض إلا في النهاية:

الرقصة (أ): التفاعل النووي ("المصنع")
هذا يتعلق بالكيمياء. إذا قطع نيوترون مسافة معينة عبر الوقود، فما هي احتمالية اصطدامه بذرة يورانيوم وانقسامها؟

يحسب لاموريو "المسافة السحرية" ( $\ell_c$ ). وهي متوسط المسافة التي يجب أن يقطعها النيوترون داخل المادة لضمان أنه، في المتوسط، ينتج ما يكفي من النيوترونات الجديدة لتعويض تلك المفقودة بسبب الامتصاص والتسرب.
تشبيه: تخيل أنك تسير في غابة. يجب أن تمشي مسافة كافية لتصطدم بما يكفي من الأشجار لتبني منزلاً جديداً لنفسك، وإلا ستنفد منك الأخشاب.

الرقصة (ب): المشية العشوائية ("القلعة المنفوخة")
هذا يتعلق بالفيزياء والهندسة. النيوترونات لا تسير في خطوط مستقيمة؛ بل ترتد عن الذرات مثل كرات البينبول في آلة "الباشينكو".

بسبب ارتدادها، تسلك مسارات ملتوية وطويلة للوصول إلى خارج الكرة، حتى لو كانت صغيرة.
يستخدم لاموريو قاعدة إحصائية بسيطة: إذا اتخذت خطوات عشوائية، فإن المسافة التي تنتهي عندها من نقطة البداية ترتبط بالجذر التربيعي لعدد الخطوات التي اتخذتها.
تشبيه: إذا كنت مخموراً وتحاول الخروج من غرفة، فلن تسير في خط مستقيم. ستتعثر يساراً، يميناً، للأمام، وللخلف. يستغرق الأمر منك وقتاً أطول للوصول إلى الباب مما لو كنت تسير في خط مستقيم. هذا "التسكع" يبقي النيوترونات بالداخل لفترة أطول، مما يسمح للتفاعل بالحدوث.

3. تجميع الأجزاء معاً

تربط الورقة بين هاتين الرقصتين:

الخطوة 1: حساب "المسافة السحرية" اللازمة لموازنة التفاعل النووي (الرقصة أ).
الخطوة 2: حساب حجم الكرة المطلوب بحيث تجبر عملية "التسكع" (الرقصة ب) النيوترونات على قطع تلك "المسافة السحرية" المحددة قبل أن تهرب.

النتيجة هي معادلة بسيطة تخبرك بـ نصف القطر الحرج. وبمجرد حصولك على نصف القطر، يمكنك بسهولة حساب الكتلة (مدى ثقل الكرة).

4. لماذا هذا مهم (ولماذا هو أمر مفاجئ)

يعترف المؤلف بأن هذه الطريقة "تعليمية" (pedagogical). هو لم يستخدم "معادلة الانتشار" الثقيلة والمعقدة التي يستخدمها مهندسو النووي عادةً، بل استخدم الإحصاء البسيط.

النتيجة الصادمة:
على الرغم من أن طريقته بسيطة، إلا أنها دقيقة بشكل مذهل.

بالنسبة لـ البلوتونيوم-239، توقعت رياضيات بسيطة منه كتلة حرجة قدرها 10.0 كجم.
محاكاة الحاسوب الفائق (MCNP) تقول إنها 10.2 كجم.
هذا فرق يقل عن 2%!

يثبت هذا أنك لست بحاجة إلى دكتوراه في الرياضيات المتقدمة لفهم الفيزياء وراء كيفية عمل القنبلة النووية أو لماذا يحتاج المفاعل إلى حجم معين. المفهوم الأساسي هو مجرد موازنة معدل الإنشاء مقابل معدل الهروب.

5. الآثار في العالم الحقيقي

الأمان والانتشار: تساعدنا هذه الطريقة في تقدير كمية اليورانيوم المخصب المطلوبة لصنع قنبلة بسرعة. وهي توضح أنه إذا كان لديك وقود "غير نقي" أو أقل نقاءً (مختلط بنظائر أخرى)، فإن الحجم المطلوب يصبح أكبر بكثير، مما يجعل تصنيع السلاح أمراً أصعب.
التحقق من التصميم: يمكن للمهندسين استخدام هذه الرياضيات "التقديرية السريعة" للتحقق بسرعة مما إذا كانت محاكاة الحاسوب المعقدة منطقية. إذا قال الحاسوب إن الكتلة الحرجة هي 500 كجم بينما تقول هذه الرياضيات البسيطة 10 كجم، فأنت تعلم أن هناك خطأ ما في المحاكاة.
المقياس الزمني: تقدر الورقة أيضاً مدى سرعة حدوث الانفجار النووي. إنه سريع للغاية — حوالي نصف ميكروثانية (0.0000005 ثانية). لهذا السبب يجب أن يكون توقيت آلية الانفجار مثالياً؛ لأنه إذا انفجرت القنبلة نفسها بعيداً عن بعضها ولو بجزء ضئيل جداً من الثانية، فإن التفاعل المتسلسل يتوقف قبل أن يطلق كامل طاقته.

الملخص

تذكرنا ورقة لاموريو بأن الطبيعة غالباً ما تتبع قواعد بسيطة. من خلال فصل "التفاعل الكيميائي" عن "هندسة الحركة"، أظهر أن الكتلة الحرجة لكرات المواد النووية هي مجرد مسأ matter تتعلق بضمان أن النيوترونات لديها الوقت الكافي للارتداد والقيام بمهمتها قبل أن تخرج من الباب. إنها طريقة جميلة وبسيطة لفهم واحدة من أكثر الظواهر تعقيداً في الفيزياء.

إليك ملخص تقني مفصل لورقة البحث بعنوان "طريقة أولية لتحديد الكتلة الحرجة لكرات من المواد الانشطارية بناءً على فصل عمليات النقل النووي والتفاعلات النووية" للكاتب س. ك. لامورو (S.K. Lamoreaux).

1. بيان المشكلة

تتناول الورقة التحدي المتمثل في حساب الكتلة الحرجة (والنصف قطر الحرج) لكرات من المواد الانشطارية (مثل $^{235}\text{U}$ و $^{239}\text{Pu}$ ). تقليديًا، تتطلب هذه المشكلة حل معادلة انتشار النيوترونات، وهي معادلة تفاضلية جزئية معقدة غالبًا ما تستلزم تقنيات رياضية متقدمة وبرمجيات عددية (مثل MCNP). يهدف المؤلف إلى تطوير طريقة تعليمية أولية تهدف إلى:

تجنب حل معادلة الانتشار.
ట استخدام الحسابات التفاضلية والحجج الإحصائية الأولية فقط.
فصل فيزياء التفاعل النووي (المقاطع العرضية، حصيلة الانشطار) عن ميكانيكا النقل الهندسي (المسار العشوائي، التسرب).
توفير تعبير مغلق للصيغة الرياضية للنصف قطر الحرج يظل دقيقًا (في حدود بضعة بالمائة) مقارنة بالمحاكاة المتطورة.

2. المنهجية

يشتق المؤلف نصف القطر الحرج من خلال موازنة إنتاج النيوترونات مقابل فقدان النيوترونات (الامتصاص والتسرب) تحت شرط ثبات عدد النيوترونات ( $N = \text{constant}$ ). تنقسم الاشتقاقات إلى جزأين متميزين:

أ. مكون التفاعل النووي (تحديد طول المسار الحرج $\ell_c$ )

يفترض المؤلف وجود غاز نيوتروني متجانس داخل الكرة.

المنهج الاحتمالي: احتمال انتقال النيوترون مسافة $\ell$ دون أن يُمتص أو يسبب انشطارًا هو $P(\ell) = e^{-n\sigma_0 \ell}$ ، حيث $n$ هي الكثافة العددية النووية و $\sigma_0 = \sigma_a + \sigma_f$ هو إجمالي المقطع العرضي للامتصاص.
معادلة التوازن: من أجل الوصول إلى الحالة الحرجة، يجب أن يكون عدد النيوترونات الناتجة عن الانشطار مساويًا لمجموع النيوترونات المفقودة بسبب الامتصاص غير الانشطاري والنيوترونات المتسربة خارج الكرة.
- النيوترونات المنتجة: $N_f = \nu N \frac{\sigma_f}{\sigma_0} (1 - e^{-n\sigma_0 \ell})$
- النيوترونات المفقودة (امتصاص + تسرب): $N_0 + N = N(1 - e^{-n\sigma_0 \ell}) + N$
شرط العتبة: يؤدي ضبط التغير الصافي إلى الصفر إلى معادلة لطول المسار الحرج $\ell_c$ :
$\ell_c = -\frac{1}{n\sigma_0} \ln \left( 1 - \frac{\sigma_0}{\nu \sigma_f - \sigma_0} \right)$
تعزل هذه الخطوة الخصائص النووية ( $\nu, \sigma_f, \sigma_a$ ) عن الهندسة.

ب. مكون النقل (ربط طول المسار بنصف القطر $R_c$ )

يصوّر المؤلف حركة النيوترونات كـ مسار عشوائي لربط طول المسار الخطي $\ell_c$ بنصف القطر الفيزيائي $R_c$ .

إحصاءات المسار العشوائي: يخضع النيوترون لـ $N_s$ من خطوات التشتت بمتوسط مسار حر $\ell_s = 1/(n\sigma_{tr})$ ، حيث $\sigma_{tr}$ هو المقطع العرضي لنقل التشتت.
التصحيح الهندسي: متوسط مربع المسافة من نقطة عشوائية داخل كرة إلى السطح هو $4R^2/5$ . ومع ذلك، فإن جذر متوسط مربع الإزاحة للمسار العشوائي يبالغ في تقدير الإزاحة الفعلية.
عامل التصحيح ( $\epsilon$ ): يتم إدخال عامل تصحيح عديم الأبعاد $\epsilon \approx 0.89$ لأخذ الهندسة والطبيعة الإحصائية للمسار العشوائي في الاعتبار.
العلاقة النهائية: من خلال الجمع بين خطوات المسار العشوائي ( $N_s \ell_s = \ell_c$ ) والقيد الهندسي ( $R^2 \approx N_s \epsilon^2 \ell_s^2$ ) نصل إلى صيغة نصف القطر الحرج:
$R_c = \epsilon \sqrt{\ell_c \ell_s} = \frac{\epsilon}{n\sqrt{\sigma_{tr}\sigma_0}} \left[ -\ln \left( 1 - \frac{\sigma_0}{\nu \sigma_f - \sigma_0} \right) \right]^{1/2}$

التقريبات الرئيسية

استقلال الطاقة: يتم التعامل مع المقاطع العرضية و $\nu$ كمتوسطات مستقلة عن الطاقة حول 1 ميجا إلكترون فولت (طاقة نيوترونات الانشطار النموذجية).
التجانس: يُفترض أن كثافة النيوترونات موحدة داخل الكرة.
مقطع النقل العرضي: يتم تقريب $\sigma_{tr}$ كـ $\sigma_e + \frac{1}{2}\sigma_i$ (المرن + نصف غير المرن)، مع مراعاة التشتت الموجه للأمام.

3. المساهمات الرئيسية

فصل الفيزياء عن الهندسة: نجحت الورقة في فصل حركية التفاعل النووي المعقدة عن مشكلة النقل المكاني، مما يسمح بفهم معياري للحرجة.
اشتقاق أولي: توفر حلاً تحليلياً مغلقاً دون الحاجة إلى حل معادلة الانتشار أو المعادلات المتسامية المرتبطة عادةً بالشروط الحدية.
معالجة شرط حدودي جديد: في الملحق (D)، يوضح المؤلف أن تطبيق شرط حدي محدد على معادلة الانتشار يعطي نتيجة مطابقة تقريبًا لطريقة المسار العشوائي الأولية، مما يؤكد صحة النهج المبسط.
عامل التصحيح: يعد اشتقاق العامل $\epsilon \approx 0.89$ (الذي يجمع بين العامل الهندسي $4/5$ وتصحيح المسار العشوائي الإحصائي) مساهمة تقنية مبتكرة تزيد من دقة التقديرات مقارنة بتقديرات المسار العشوائي البسيطة.

4. النتائج

تم اختبار الطريقة مقابل الكتل الحرجة المعروفة ومحاكاة MCNP6 عالية الدقة باستخدام بيانات المقاطع العرضية ENDF/B-VIII.0:

المادة	$R_c$ المحسوب (سم)	$m_c$ المحسوبة (كجم)	MCNP6 / الأدبيات $m_c$ (كجم)	الاتفاق
$^{239}\text{Pu}$ (1.9 MeV)	4.96	10.00	10.2	< 2%
$^{235}\text{U}$ (1.5 MeV)	8.37	46.5	46.4	< 0.5%
$^{235}\text{U}$ (93.7% enriched)	-	46.5	49.0 (المرجع [10])	جيد (المرجع [10] يستخدم تفاصيل تخص نسبة الإثراء)

تحليل الحساسية: تحلل الورقة تأثير شوائب $^{238}\text{U}$ ومستويات الإثراء. وتوضح أنه بالنسبة لمستويات الإثراء المنخفضة، يصبح تدهور طاقة النيوترونات (التهدئة) مهمًا، مما يتطلب تعديلات في المقطع العرضي للانشطار ( $\sigma_f$ ) لمطابقة البيانات التجريبية.
المقارنة مع أوبنهايمر-بيث: تمت مقارنة الصيغة المستمدة (المعادلة 11) بصيغة أوبنهايمر-بيث التاريخية. وتوفر طريقة المؤلف نتائج أكثر دقة (على سبيل المثال، 10.0 كجم مقابل 8.6 كجم للبلوتونيوم) مع الحفاظ على اشتقاق أكثر شفافية.

5. الأهمية

القيمة التعليمية: تعتبر هذه الطة مثالية لتدريس الفيزياء النووية ودراسات عدم الانتشار، حيث تشرح الآليات الفيزيائية للحرجة (توازن الإنتاج مقابل التسرب) دون اللجوء إلى "الصندوق الأسود" للمحللات العددية المعقدة.
دليل تصميم: تعمل كأداة سريعة للتحقق من رتبة المقدار (back-of-the-envelope) في تصميمات النيوترونات وللتحقق من نتائج أكواد المحاكاة واسعة النطاق مثل MCNP.
الرؤية الفيزيائية: يوضح التحليل أن تهدئة النيوترونات (الإبطاء) تلعب دورًا ثانويًا في التجمعات الحرجة السريعة لأن عمر النيوترون قصير جدًا؛ حيث يهيمن الطيف على نيوترونات الانشطار الجديدة عالية الطاقة.
التوسع: أظهرت الورقة إمكانية تكييف الإطار مع هندسات أخرى (مثل القضبان) والأنظمة دون الحرجة، بما في ذلك تقدير التضاعف النيوتروني في تجارب الانشطار المستحث بالنيوترينو.

في الختام، يقدم لامورو إطارًا تحليليًا مبسطًا وقويًا يحقق دقة عالية في التنبؤ بالكتل الحرجة، مما يسد الفجوة بين مفاهيم الفيزياء الأولية وحسابات الهندسة النووية المتطورة.