ASP-Assisted Symbolic Regression: Uncovering Hidden Physics in Fluid Mechanics

ملخص تقني: الانحدار الرمزي المدعوم ببرمجة التقييد المنطقي (ASP) في ميكانيكا الموائع

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

المنهجية
تستخدم الأبحاث منهجية ذات شقين: نهج انحدار رمزي قياسي يليه إطار عمل هجين مبتكر يدمج برمجة التقييد المنطقي (ASP).

1. توليد البيانات والإعداد:

• تتضمن المشكلة الفيزيائية تدفقاً صفيحياً ثلاثي الأبعاد في مجرى مستطيل متماثل ($L=5, H=W=1$).
• المعادلات الحاكمة هي معادلات نافيير-ستوكس غير البعدية ومعادلة حفظ الكتلة.
• تم توليد الحلول العددية باستخدام طريقة الحجم المحدود (FVM) المنفذة في برنامج MATLAB على شبكة مكونة من $50 \times 25 \times 25$.
• يتم دفع التدفق بواسطة تدرج ضغط ثابت في اتجاه التيار ($c$)، والذي يعمل كمعلمة تحكم رئيسية. تم حساب رقم رينولدز الفعال ($Re$) لاحقاً من الحل العددي ليكون بمثابة معلمة تحكم قياسية لنماذج الانحدار الرمزي.
• تشمل مجموعات البيانات مجموعات التدريب ($Re \in \{34, 105, 174, 209, 279\}$) ومجموعات الاختبار ($Re \in \{70, 139, 244\}$).

2. مرحلة الانحدار الرمزي (SR):

• تُستخدم مكتبة PySR (الإصدار 0.19.4) لاكتشاف التعبيرات الرمزية للسرعة المحورية ($u$) والضغط ($p$) كدوال للإحداثيات المكانية ($X, Y, Z$) و $Re$.
• العمليات الأولية: تم تقييد فضاء البحث بالعمليات المستوحاة من المجال: العمليات الحسابية الأساسية ($+, -, *, /$) والحدود متعددة الحدود من الدرجة المنخفضة (التربيع، التكعيب).
• القيود: لضمان القابلية للتفسير، تم وضع حد أقصى لتعقيد التعبير عند 25 عملية. كما تمنع القيود المتداخلة الوظائف عالية الرتبة المكررة.
• التحسين: يتم تطوير النماذج عبر 1000 تكرار لتقليل متوسط مربع الخطأ (MSE)، مما يولد جبهة باريتو (Pareto front) من الحلول التي توازن بين الدقة والبساطة.

3. مرحلة الهجين (SR/ASP):

• لمعالجة قصور النماذج المعتمدة على البيانات فقط، قام المؤلفون بدمج الانحدار الرمزي (SR) مع برمجة التقييد المنطقي (ASP)، وهي نموذج برمجة تصريحية لإرضاء القيود.
• سير العمل: يقوم الانحدار الرمزي بتوليد المعادلات المرشحة. ثم يتم تحليل هذه المعادلات إلى حقائق (على سبيل المثال: التعقيد، الخسارة، وجود مصطلحات محددة مثل $Re$ أو $Y^2$).
• التصفية: يقوم برنامج ASP بترميز القيود الفيزيائية الخاصة بالمجال (مثل استبعاد المصطلحات غير الفيزيائية مثل $X^3$ أو $Y^3$، وضرورة تضمين $Re$، والالتزام بشروط عدم الانزلاق، ومتطلبات التماثل).
• الحل: يقوم برنامج حل ASP (وهو CLASP) بمعالجة هذه الحقائق والقيود لتصفية مخرجات الانحدار الرمزي، والاحتفاظ فقط بالمعادلات التي تكون دقيقة إحصائياً وصحيحة فيزيائياً.

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

• اشتقاق القوانين الفيزيائية الموجزة: نجحت الدراسة في اشتقاق معادلات رمزية صريحة مغلقة الشكل لسرعة التدفق ثلاثي الأبعاد والضغط، والتي تطابق الحلول العددية والاتجاهات التحليلية الراسخة بشكل وثيق.
• التكامل المبتكر بين SR و ASP: تقدم الورقة، حسب علم المؤلفين، أول جهد لدمج الانحدار الرمزي مع برمجة التقييد المنطقي. يستفيد هذا الإطار الهجين من القدرة التوليدية للانحدار الرمزي والاستنتاج التصريحي لـ ASP لضمان التزام النماذج المشتقة بالمبادئ الفيزيائية.
• التحقق من المتانة: تم التحقق من صحة النهج من خلال حقن الضجيج (إضافة ضجيج بنسبة 2% إلى البيانات) وتوسيع مجموعة العمليات (بما في ذلك الدوال المثلثية/الأسية)، مما يثبت أن العلاقات متعددة الحدود المكتشفة هي انعكاس حقيقي لفيزياء التدفق وليست مجرد نواتج عرضية لفضاء البحث.

النتائج

• نموذج السرعة: استنتج نموذج الانحدار الرمزي المعادلة $u_{symb} = Re(2.18 - 8.46 Y^2)(1 - 3.89 Z^2)$. تلتقط هذه التعبير ملف السرعة المكافئ المتوقع. وقد حقق نموذج متوسط الخطأ المطلق المعياري (NMAE) قدره 0.008% على مجموعة التدريب و 0.010% على مجموعة الاختبار.
• نموذج الضغط: معادلة الضغط المستنتجة هي $p_{symb} = Re(143.43 - 28.69 X)$، مما يعكس انخفاضاً خطياً في الضغط. حقق هذا النموذج NMAE في نطاق $10^{-8}\%$، مما يشير إلى توافق شبه مثالي مع البيانات العددية.
• المقارنة مع الأدبيات: تقترب معادلة السرعة المستنتجة بشكل وثيق من الحل التحليلي للمجاري المستطيلة الموجود في الأدبيات (على سبيل المثال، $u \propto (1-4y^2)(1-4z^2)$)، مع اختلافات طفيفة تُعزى إلى دقة البيانات العددية المحدودة.
• تصفية ASP: في مثال نموذجي، نجحت وحدة ASP في تصفية جبهة باريتو مكونة من 14 معادلة. اختارت أكثر المرشحين دقة وصحة فيزيائية (بما في ذلك معادلة السرعة النهائية) في حوالي 0.1 ثانية، وهو وقت ضئيل مقارنة بـ 15 دقيقة المطلوبة لمرحلة تدريب الانحدار الرمزي. نجحت عملية التصفية في استبعاد المعادلات التي تنتهك التماثل أو شروط الحدود.

الأهمية والادعاءات
تدعي الورقة أن منهجية SR/ASP المقترحة تثبت إمكانية الجمع بين نهج البيانات والنهج القائم على تمثيل المعرفة لتعزيز القابلية للتفسير، والموثوقية، والمواءمة الفيزيائية لنماذج الذكاء الاصطناعي في ديناميكا الموائع.

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

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