Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

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

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

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

الطريقة الأولى، المعروفة باسم "نهج الجدول" (table approach)، تقوم بحساب مسبق لقائمة من الخصائص لكل حالة ممكنة لكتلة التربة وتبحث عن الإجابة بناءً على متوسط درجة حرارة الكتلة. أما الطريقة الثانية، وهي "المخطط العقدي" (nodal scheme)، فتقرر حالة التربة عند زوايا الكتل وتقوم بمتوسط تلك القرارات. أما الطريقة الثالثة، وتسمى "تنعيم h" (h-smoothing)، فتقوم بحساب الجزء الدقيق من الكتلة الذي ذاب بناءً على تدرج حراري خطي، أي أنها تقيس فعليًا شريحة الكتلة التي تحولت إلى ماء. أجرى الباحثون هذه المحاكاة في أبعاد ثنائية وثلاثية الأبعاد، مقارنين النتائج بحل "مرجعي" عالي التفصيل يتتبع كل طبقة من طبقات التربة.

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

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

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