Topological cell-openness index for porous materials

ملخص تقني: مؤشر الانفتاح الخلوي الطوبولوجي للمواد المسامية

بيان المشكلة
يعد توصيف انفتاح المواد المسامية أمراً حاسماً لفهم خصائص النقل، ومع ذلك تعتمد المعايير الحالية بشكل كبير على قياس الكثافة بالغاز (gas pycnometry). توفر هذه الطريقة التجريبية معلماً واحداً، وهو "كسر المسام المفتوحة" ($\phi_0$)، الذي يمثل الكسر الحجمي للمسام المتصلة بحدود العينة. ومع ذلك، فإن $\phi_0$ يعاني من قيود: فهو يتطلب معرفة بالكثافة الحقيقية للمادة، ويواجه صعوبات في العينات التي لا ترتبط فيها المسام المفتوحة بالسطح، ولا يمكنه تمييز الانفصالات الهيكلية الداخلية ضمن طور "مفتوح" واحد (على سبيل المثال، وجود مكونات متعددة تلامس الحدود ولكنها غير متصلة ببعضها البعض). ومع تزايد إمكانية الوصول إلى التصوير ثلاثي الأبعاد الدقيق (مثل التصوير المقطعي المحوسب للأشعة السينية الدقيقة X-ray micro-CT)، تبرز الحاجة إلى طرق مكملة تعتمد على الصور وتلتقط معلومات الاتصال الطوبولوجي بما يتجاوز مجرد الكسور الحجمية البسيطة.

المنهجية
يقترح المؤلفون طريقة تعتمد على الصور باستخدام تحليل البيانات الطوبولوجية (TDA)، وتحديداً التماثل المستمر (persistent homology)، لتقدير نسبة الخلايا المفتوحة والمغلقة في البنى المجهرية الثنائية (binary) ثنائية وثلاثية الأبعاد (2D/3D).

1. توليد البيانات والمصادر: تستخدم الدراسة بنى مجهرية ثنائية اصطناعية تم توليدها عبر عملية عشوائية تتضمن مساماً كروية/دائرية على شبكة، متصلة بقنوات ذات احتمالات اتصال ($p$) متفاوتة. تم إنشاء أربع عائلات من مجموعات البيانات: ثقوب متصلة جزئياً (2D/3D)، خلايا مغلقة تماماً (2D)، خلايا مفتوحة تماماً (2D)، وخلايا مختلطة ثلاثية الأبعاد لتحليل انتقال الحرارة. كما استُخدمت مجموعة بيانات ثنائية الأبعاد خارجية مكونة من 31 صورة مسامية ذات قيم نفاذية معروفة.
2. القياس الحجمي الرقمي ($\phi_0$): تم تنفيذ نظير رقمي لقياس الكثافة بالغاز من خلال تحديد المكونات المتصلة لطور المسام التي تلامس حدود المجال وحساب كسر حجمها بالنسبة لإجمالي حجم المسام.
3. التحليل الطوبولوجي: يحسب المؤلفون أعداد بيتي (Betti numbers) ($\beta_0, \beta_1, \beta_2$) على المجموعات الفرعية لـ "تحويل المسافة الموقعة" (Signed Distance Transform - SDT) للبنية الثنائية.
   • $\beta_0$: عدد المكونات المتصلة.
   • $\beta_1$: عدد الحلقات المستقلة (1D) أو الأنفاق.
   • $\beta_2$: عدد الفجوات المحصورة (في 3D).
     يتم تصفية الضوضاء من خلال الاحتفاظ فقط بالميزات التي تبلغ درجة استمراريتها (persistence) $\ge 1.5$.
4. مؤشر الانفتاح الخلوي ($\tau$): يتم تعريف مؤشر قياسي جديد $\tau \in [0, 1]$ بناءً على أعداد بيتي عند قيمة ترشيح محددة ($t=-1$):
   • في الحالة ثنائية الأبعاد (2D): $\tau = \frac{\beta_1}{\beta_0 + \beta_1}$
   • في الحالة ثلاثية الأبعاد (3D): $\tau = \frac{\beta_1 + \beta_2}{\beta_0 + \beta_1 + \beta_2}$
     يشير $\tau = 0$ إلى نظام من المسام المنفصلة والمغلقة؛ بينما يشير $\tau = 1$ إلى شبكة متصلة ومفتوحة بالكامل.
5. المحاكاة الفيزيائية: تم حساب الموصلية الحرارية الفعالة للهياكل ثلاثية الأبعاد عن طريق حل معادلة الانتشار المستقر مع موصلية متجانسة المناحي (piecewise isotropic). وتم أخذ بيانات النفاذية من مجموعة البيانات المرجعية الخارجية.

المساهمات والنتائج الرئيسية

• تعريف $\tau$: يقدم البحث $\tau$ كمكمل طوبولوجي لـ $\phi_0$. وبينما يوجد ارتباط وثيق بين $\phi_0$ و $\tau$، فإن $\tau$ يلتقط تفاصيل الاتصال التي يفتقدها $\phi_0$. وتحديداً، في الحالات التي يصل فيها $\phi_0$ إلى حالة التشبع بالقرب من 1 (مما يشير إلى انفتاح عالٍ)، يمكن لـ $\tau$ التمييز بين الشبكة المتصلة حقاً وبين تلك المكونة من عدة مكونات تلامس الحدود فردياً ولكنها غير متصلة ببعضها البعض.
• الارتباط مع الخصائص الفيزيائية:
  • النفاذية في 2D: في مجموعة بيانات مكونة من 31 صورة مسامية، أظهر $\tau$ علاقة مكافئة (parabolic) قوية مع $\log_{10}(\text{permeability})$ ($R^2 = 0.95$)، رغم أن $\tau$ تغير بأقل من 1% عبر مجموعة البيانات.
  • الموصلية الحرارية في 3D: عبر 250 هيكلاً ثلاثي الأبعاد اصطناعياً، أظهر $\tau$ ارتباطاً موجباً متوسطاً مع أثر مؤثر الموصلية الحرارية الفعالة ($\text{tr}(K)$). وفي جميع المجموعات الفرعية، كان $\tau$ متنبئاً أفضل لـ $\text{tr}(K)$ من $\phi_0$.
• تقدير المقياس الطولي عبر منحنيات بيتي: يوضح المؤلفون أن منحنيات بيتي ($\beta_i(t)$) المستمدة من ترشيحات SDT يمكنها استعادة الأطوال الهندسية المميزة:
  • الخلايا المغلقة: نجحت نقاط أقصى تدرج في $\beta_0$ و $\beta_1$ في تقدير نصف قطر المسام المتوسط ($r$)، ونصف عرض ما بين المسام ($d/2$)، ونصف سمك اللب الصلب ($h/2$) بدقة عالية ($R^2$ تصل إلى 0.95).
  • الخلايا المفتوحة: بينما تعمل هذه الطريقة بفعالية في الهياكل الدورية المثالية، كانت القدرة التنبؤية لعرض العنق ($w/2$) في مجموعات البيانات المفتوحة الواقعية والمشوبة بالضوضاء أضعف، ويرجع ذلك على الأرجية إلى تداخل المقاييس الطولية الذي يحجب البصمات الطوبولوجية.

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

يدعي البحث ما يلي:

1. يوفر $\tau$ معلومات هيكلية حقيقية فيما يتعلق بالاتصال الداخلي التي لا يستطيع $\phi_0$ تمييزها.
2. إن عدم التطابق بين $\tau$ و $\phi_0$ في الأنظمة شديدة الانفتاح ليس خطأً، بل هو ميزة تشير إلى وجود شبكات فرعية غير متصلة.
3. توفر منحنيات بيتي مساراً لتقدير أحجام الميزات المميزة (نصف قطر المسام، عرض العنق) مباشرة من البصمات الطوبولوجية، وإن كان ذلك أكثر قوة في أنظمة الخلايا المغلقة.
4. يرتبط $\tau$ بخصائص النقل (النفاذية والموصلية الحرارية)، مما يشير إلى أنه يرمز لبصمات هيكلية ذات معنى فيزيائي.

يخلص المؤلفون إلى أنه في حال توفر تصوير ثلاثي الأبعاد عالي الجودة، ينبغي تسجيل $\tau$ جنباً إلى جنب مع $\phi_0$ لتوفيد توصيف أكثر اكتمالاً للبنى المجهرية للمواد المسامية.