Dissolution of carbonate stones caused by CO2 pollutant: an erosion model

إليك ملخص تقني مفصل للورقة البحثية بعنوان "تحلل الأحجار الكربوناتية الناتج عن ملوث ثاني أكسيد الكربون: نموذج تآكل".

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

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

• السياق: تتعرض المباني والمنحوتات القديمة لعمليات تجوية بمعدلات متسارعة بسبب زيادة مستويات ثاني أكسيد الكربون في الغلاف الجوي والأمطار الحمضية (التي يقل أسها pH < 5).
• الآلية: يذوب ثاني أكسيد الكربون ($CO_2$) في الماء لتكوين حمض الكربونيك ($H_2CO_3$)، والذي يتغلغل في المصفوفة المسامية للحجر. يتفاعل هذا الحمض مع كربونات الكالسيوم ($CaCO_3$)، مما يسبب الذوبان وما يتبعه من زيادة في المسامية حتى ينهار الهيكل المادي (التآكل).
• قصور النماذج الحالية: توفر الأدوات الإحصائية الحالية (مثل مؤشر ليبفرت للتآكل) تقديرات متوسطة بناءً على الهطول وترسب الملوثات، لكنها تفشل في مراعاة:
  • الخصائص النوعية للمادة (التجانس الكيميائي، المسامية).
    _ التطور الزمني للتدهور.
  • الظواهر المناخية الدقيقة المحلية.
  • الديناميكيات المحددة لتغلغل ثاني أكسيد الكربون.

2. المنهجية

يقترح المؤلفون نظام تفاعل ونقل غير خطي محكوم بتدفق دارسي في الأوساط المسامية لنمذجة عملية التآكل.

أ. النموذج الرياضي

يربط النظام بين ثلاثة متغيرات أساسية:

1. تشبع الماء ($\theta$): الجزء الحجمي للماء في المسام.
2. تركيز حمض الكربونيك ($c_a$): تركيز الحمض المذاب.
3. المسامية ($n$): المسامية المتطورة.

تتضمن المعادلات الحاكمة (النظام 3) ما يلي:

• نقل الرطوبة: معادلة انتشار غير خطية تعتمد على معادلة ريتشاردز، باستخدام دالة $B(\theta/n)$ لربط الضغط الشعري ومحتوى الرطوبة.
• نقل الحمض: معادلة حمل-انتشار-تفاعل تصف حركة حمض الكربونيك واستهلاكه عبر التفاعل الكيميائي.
• تطور المسامية: معادلة تفاضلية عادية حيث تزداد المسامية مع تفاعل الحمض مع مصفوفة الكربونات ($CaCO_3 + H_2CO_3 \rightarrow Ca^{2+} + 2HCO_3^-$).
• الشروط الحدية: مقترنة بشروط نيومان/ديريكليه التي تمثل التبادل بين الرطوبة والحمض بين الهواء المحيط وسطح الحجر.

دوال الامتصاص: حدد المؤلفون صيغتين لعلاقة الضغط الشعري/المحتوى المائي:

1. المتماثلة ($B'$): نموذج أبسط بـ 3 معاملات.
2. غير المتماثلة ($B'_{kP}$): نموذج أكثر تعقيداً يفصل بين النفاذية والضغط الشعري، مما يسمح بوصف فيزيائي أغنى (تمت معايرته مقابل البيانات التجريبية).

ب. الخوارزمية العددية

للتعامل مع الحدود المتحركة الناتجة عن فقدان المادة (التآكل)، طور المؤلفون مخططاً عددياً هجيناً:

• تقسيم النطاق: طريقة الفروق المحدودة (FDM) على شبكة ديكارتية منتظمة.
• طريقة المستوى المجموع (Level-Set Method): تُستخدم لتتبع الواجهة المتحركة (جبهة التآكل). يتم تضمين النطاق الفيزيائي في شبكة أكبر، وتحدد دالة المستوى $\phi(x) = n(x) - n_{max}$ الحدود.
• طريقة النقاط الشبحية (Ghost Point Method): للتعامل مع الحدود التي لا تتوافق مع الشبكة، يتم إسقاط "نقاط شبحية" على الحدود الفيزيائية على طول المتجه العمودي. يتم فرض الشروط الحدية عبر استكمال لاغرانج (ثنائي التربيع) على هذه النقاط الشبحية.
• التكامل الزمني: تُستخدم طريقة أويلر الضمنية للسماح بخطوات زمنية أكبر وتجنب قيود الاستقرار.

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

1. إطار رياضي مبتكر: تقديم نموذج تفاعل ونقل مقترن مصمم خصيصاً لذوبان ثاني أكسيد الكربون في الأحجار الكربوناتية، متجاوزاً المتوسطات الإحصائية إلى الفيزياء الميكانيكية.
2. مخطط عددي قوي: تطوير خوارزمية (Level-Set Ghost-Point Finite Difference) قادرة على محاكة هندسات ثنائية الأبعاد معقدة وتتبع جبهة التآكل دون الحاجة لإعادة التقسيم (remeshing).
3. معايرة المعاملات: تمت معايرة معاملات النموذج (ثوابت التفاعل، معاملات الانتشار، عتبات المسامية) مقابل بيانات أدبية تتعلق بمعدلات تراجع شواهد القبور في البيئات الحضرية والضواحي.
4. المرونة: يسمح الإطار بدمج البيانات البيئية في الوقت الفعلي (درجة الحرارة، الرطوبة، مستويات ثاني أكسيد الكربون) لمحاكاة ظروف موقع محدد.

4. النتائج

أجرى المؤلفون اختبارات عددية في بُعد واحد (1D) وبُعدين (2D) للتحقق من صحة النموذج.

• محاكاة البعد الواحد (1D):
  • أظهرت أن كلاً من دوال الامتصاص المتماثلة وغير المتماثلة تعطي نتائج متطابقة تقريباً بالنسبة للرطوبة ومواقع جبهة التآكل في الظروف القياسية.
  • يصل النظام إلى حالة مستقرة للرطوبة وتركيز الحمض، بينما تزداد المسامية حتى يحدث فقدان للمادة.
  • معدل التآكل: توقع النموذج معدل تراجع يبلغ حوالي $2.36 \times 10^{-6}$ سم/ساعة (حوالي 2 ملم/100 عام) للظروف في الضواحي، وهو ما يتماشى مع القيم الأدبية.
• محاكاة البعدين (2D):
  • أظهرت قدرة النموذج على التعامل مع هندسات معقدة (مثل عينة حجرية مستطيلة).
  • صورت تطور جبهة التآكل على مدار عام واحد، مما أظهر تراجعاً واقعياً للسطح.
• السيناريو الكارثي:
  • محاكاة حالة "الغمر" (التشبع الكامل).
  • أظهرت النتائج تغلغلاً سريعاً للماء ولكن مع تراكم محدود للحمض بسبب الاستهلاك عبر التفاعل، مما أدى إلى منحدر تآكل أعلى بكثير ($8.25 \times 10^{-4}$ سم/ساعة).
• التقارب: أظهر المخطط العددي رتبة دقة عالمية قدرها 1، وهو ما يتوافق مع تقسيم الوقت والمكان من الدرجة الأولى المستخدم.

5. الأهمية والعمل المستقبلي

• حفظ التراث الثقافي: يوفر النموذج أداة تنبؤية لربط مستويات التلوث بتدهور المواد، مما يمكن من وضع استراتيجيات حفظ خاصة بالمواقع (مثل منتزه أوشيا أنتيكا الأثري).
• القابلية للتوسع: على الرغم من اختباره حالياً في 1D و 2D، إلا أن الخوارزمية مصممة ليتم توسيعها بسهولة إلى 3D.
• الاتجاهات المستقبلية:
  • معايرة المعاملات ($K_w, K_a, K_n, n_{max}$) باستخدام بيانات الاستشعار في الموقع.
  • التوسع في عمليات المحاكة ثلاثية الأبعاد (3D).
  • دمج تأثيرات الجاذبية (الخاصية الشعرية الصاعدة) ودورات التجمد والذوبان.
  • نمذجة التفاعل بين التآكل وتكون القشور (الترسيب).
  • رفع رتبة المخطط العددي (على سبيل المثال، Crank-Nicolson) لتحقيق دقة أعلى.

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