On the closed solution of a problem coupling fluid infiltration with a hydration reaction

ملخص تقني: حل مغلق لمسألة اقتران ترشيح السوائل والتميؤ

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

يفترض النظام الفيزيائي ما يلي:

1. وسط مسامي متجانس بمسامية $\phi$.
2. جبهة ترشح حادة تفصل بين المناطق المشبعة والجافة (افتراض "جرين-أمبت" Green-Ampt).
3. تفاعل تميؤ من الدرجة الأولى يستهلك الماء المتسرب.
4. تدفق "دارسي" (Darcy flow) داخل المنطقة المشبعة ($0 < x < s(t)$).

تجمع المعادلات الحاكمة بين قانون دارسي، وحفظ الكتلة، وحركية التفاعل. يحقق الضغط $p(x,t)$ معادلة من نوع "بواسون" حيث يمثل حد المصدر استهلاك الماء بفعل التفاعل. ويُغلق النظام بشرط من نوع "ستيفان" عند الجبهة المتحركة $s(t)$، يربط تدرج الضغط بسرعة الجبهة، مع شروط حدودية تثبت الضغط عند الأصل وتجعله صفراً عند الجبهة.

المنهجية
اشتق المؤلفون صيغة لا بُعدية للمسألة، حيث استُخدم مقياس طول $\ell^*$ ومقياس زمن $t^*$ بناءً على النفاذية، واللزوجة، ورأس الضغط، ومعدل التفاعل. يتضمن النظام لا البعدي إيجاد موضع الجبهة الحرة $s(t)$، ومجال الضغط $p(x,t)$، ومتغير تقدم التفاعل $\Gamma(x,t)$.

للحصول على حل مغلق الصيغة، استخدم المؤلفون نهجاً تحليلياً يرتكز على دالة مساعدة $J(t)$، تُعرف بأنها تكامل متغير تقدم التفاعل عبر النطاق المشبع. ومن خلال دمج المعادلات الحاكمة وتطبيق نظرية "فوبيني" (Fubini's theorem)، اشتقوا معادلة تفاضلية عادية (ODE) خطية لـ $J(t)$. ويسمح حل هذه المعادلة بالتحديد الصريح لموضع الجبهة الحرة $s(t)$. وبمجرد معرفة $s(t)$، يتم إعادة بناء متغير التفاعل $\Gamma$ ومجال الضغط $p$ باستخدام معكوس الدالة الرتيبة الصارمة $s(t)$ وتكامل معادلة تدرج الضغط.

المساهمات والنتائج الرئيسية
تتمثل المساهمة الأساسية في اشتقاق حل تحليلي دقيق ومغلق الصيغة للجبهة الحرة $s(t)$ والمجالات المرتبطة بها، وهو حل صالح لجميع الأوقات $t \geq 0$. الحل فريد ومتزايد بشكل صارم.

يكشف التحليل عن أنماط ديناميكية متميزة لجبهة الترشح بناءً على العلاقة بين المسامية $\phi$ ومعلمة استهلاك الماء لا بُعدية $Y$ (التي تمثل الحجم الإجمالي للماء المستهلك لكل وحدة حجم من الوسط):

1. سلوك الوقت المبكر: بالنسبة لـ $t \to 0$، يُظهر الحل عالمياً اعتماداً على الجذر التربيعي للزمن ($s(t) \sim \sqrt{t/\phi}$)، بغض النظر عن قوة تفاعل التميؤ. وهذا يطابق سلوك الترشح الكلاسيكي في غياب التفاعل، حيث تتقدم الجبهة بسرعة أكبر مما يمكن للتفاعل استهلاكه من ماء.
2. سلوك الوقت المتأخر: يعتمد السلوك التقاربي عندما $t \to \infty$ على إشارة $\phi + Y$:
   • في حالة $\phi + Y > 0$: تتقدم الجبهة بسلوك يعتمد على الجذر التربيعي للزمن ($s(t) \sim \sqrt{t/(\phi+Y)}$)، ولكن بمعامل معدل يتحدد بكل من المسامية واستهلاك التفاعل.
   • في حالة $\phi + Y = 0$: تتقدم الجبهة خطياً مع الزمن ($s(t) \sim t$).
   • في حالة $\phi + Y < 0$: يظهر تقدم الجبهة نمواً أسياً.

توضح الورقة صراحةً وجود "انتقال" (crossover) من نظام الجذر التربيعي في الوقت المبكر إلى النظام في الوقت المتأخر، مستردةً تحليلياً السلوكيات التي لوحظت سابقاً عددياً فقط في الدراسات [2, 9].

الأهمية
تزعم الورقة أن هذا الحل المغلق يوفر أداة تحليلية دقيقة لقياس أنظمة النقل والتفاعل في باطن الأرض. وتحديداً:

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

يؤكد المؤلفون أنه بينما تفرض الفيزياء أن $\phi > 0$ و $Y > 0$ (استهلاك الماء)، فإن التحليل الرياضي يسري على نطاق أوسع من المعلمات، بما في ذلك الحالات التي يكون فيها $Y < 0$ (إطلاق الماء)، رغم الإشارة إلى أنها سيناريوهات متميزة فيزيائياً.