ملخص تقني: قواعد ضغط تكيفية موجهة بالبواقي للقواعد المقلصة الأبعاد في إسقاط FLIP

بيان المشكلة

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

المنهجية

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

1. بناء القاعدة التكيفية

يتم بناء فضاء حل الضغط من خلال الجمع بين:

• قاعدة سطحية خشنة عالمية ($U$) تتعامل مع معظم النطاق.
• قواعد هوية دقيقة محلية ($U_a$) تُدرج بشكل انتقائي في النظام المُقلص.

تُعرف القاعدة الكلية بأنها $B = \hat{\Phi}[U \ U_a]$، حيث تتضمن $\hat{\Phi}$ شروط ديريكليه الحدية للسطح الحر. ويُمثل حقل الضغط كـ $p = Bq$، حيث تحتوي $q$ على معاملات لكل من المتغيرات الخشنة والتكيفية الدقيقة.

2. معايير الاختيار (حقل الأهمية)

يتم اختيار الخلايا للتحسين التكيفي بناءً على حقل أهمية مركب $I_i$، والذي يجمع بين ثلاثة مؤشرات:

• متنبئ البواقي الزمني ($R_i$): الإشارة الأساسية هي باقي الشبكة الدقيقة لحل الضغط المُقلص من الخطوة الزمنية السابقة. لضمان التماسك الزمني وتجنب الوميض بين الإطارات، يتم تحديث هذا الباقي باستخدام مرشح زمني أسي: $\bar{r}^n = \alpha_r |r^n| + (1-\alpha_r)\bar{r}^{n-1}$. وهذا يسمح ببناء القاعدة قبل الحل الحالي، باستخدام بيانات الخطأ التاريخية.
• انحناء السطح ($K_i$): مقدار انحناء السطح الحر، مُطبع ضمن نطاق ضيق. يعمل هذا على اكتشاف الهياكل عالية الانحناء التي تشكلت حديثاً (مثل الرذاذ) والتي قد لا تمتلك تاريخاً للبواقي.
• تغير السرعة ($G_i$): معيار فروبينيوس لتدرج السرعة المحلي، والذي يحدد المناطق ذات التباين الكينماتيكي القوي.

يُحسب سجل الأهمية كـ $I_i = w_r R_i + w_\kappa K_i + w_g G_i + S_i$، حيث $S_i$ هو حد أولوية السطح الذي يضمن حداً أدنى من مستوى التحسين عند السطح الحر.

3. إدارة الميزانية

لمنع نمو النظام المُقلص بشكل مفرط، يتم فرض ميزانية درجات حرية (DOF) صارمة. يتم فرز المرشحين الذين يتجاوزون عتبة $\tau_I$ حسب الأهمية. تختار الخوارزمية "البذور" عالية الأهمية وتوسع دعمها إلى الخلايا المجاورة (متصلة بـ 4 في 2D، و6 في 3D) حتى يصل عدد درجات الحرية التكيفية إلى حد محدد مسبقاً (على سبيل المثال، 3.5% من خلايا السائل أو سقف صلب قدره 120 درجة حرية إضافية).

4. مسار التنفيذ

1. إعادة بناء الحالة: نقل سرعات الجسيمات إلى الشبكة وإعادة بناء مستوى الطقم (level set).
2. تقييم المؤشرات: حساب المتنبئ الباقي (من الخطوة السابقة)، والانحناء، وتدرجات السرعة لتحديد حقل الأهمية.
3. اختيار القاعدة: اختيار الخلايا التكيفية تحت ميزانية درجات الحرية وبناء القاعدة المعززة $B$.
4. الحل المُقلص: حل النظام المُقلص $B^T A B q = B^T b$ باستخدام طريقة التدرج المترافق المسبقة الشرط (preconditioned conjugate gradient).
5. الاستطالة والتصحيح: استطالة الحل إلى الشبكة الدقيقة وتطبيق عدد ثابت (3) من تكرارات ج Jacobi لتقليل التباعد في النطاق الدقيق.
6. تحديث البواقي: حساب باقي الشبكة الدقيقة الجديد (قبل ج Jacobi) لتحديث المتنبئ الزمني للخطوة التالية.

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

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

النتائج التجريبية

تم تقييم الطريقة عبر 11 مشهداً متنوعاً في 2D (مثل انهيارات السدود، الأمواج المتصادمة، النوافير) وتم توسيعها إلى 3D.

• الدقة: مقارنة بالتقليص المنتظم للأبعاد (بخطوة 4)، قللت الطريقة المقترحة متوسط الخطأ النسبي لتقريب الضغط من 0.127 إلى 0.049 (أي بنسبة تقليل ~61%) في 2D. وفي 3D، انخفض الخطأ من 0.134 إلى 0.057.
• التباعد: انخفض التباعد بعد الإسقاط بشكل كبير. في 2D، انخفض متوسط جذر متوسط المربعات (RMS) للتباعد من $1.82 \times 10^{-3}$ (التقليص المنتظم) إلى $6.7 \times 10^{-4}$ (الطريقة المقترحة)، مما قلص الفجوة نحو الحل كامل الدقة.
• الكفاءة: تضيف الطة عدداً صغيراً فقط من درجات الحرية التكيفية (~84 في 2D، ~108 في 3D) إلى النظام الخشن. وبالتالي، زاد وقت حل الضغط بنسبة ~25% فقط عن طريقة التقليص المنتظم (1.58 مللي ثانية مقابل 1.26 مللي ثانية في 2D) ولكنه ظل أسرع بنسبة ~73% من الإسقاط كامل الدقة (5.84 مللي ثانية).
• دراسة الاستئصال (Ablation): أدى إزالة حد البواقي إلى زيادة الخطأ إلى 0.081، مما يؤكد أن الباقي هو الإشارة المهيمنة. وأدى إزالة حدود الانحناء/السرعة إلى زيادة الخطأ إلى 0.061، مما يوضح قيمتها في اكتشاف الهياكل الجديدة. حقق الجمع الكامل للمؤشرات أفضل النتائج.
• المتنبئ الزمني: استخدام متنبئ باقي مُفلتر قلل من معدل دوران الخلايا التكيفية إلى 17.9% (مقارنة بـ 29.7% للبواقي غير المفلترة من الإطار السابق) مع الحفاظ على دقة مماثلة لنهج "التحسين في نفس الإطار + إعادة الحل" الأكثر تكلفة.

الأهمية والادعاءات

يدعي البحث أن الطريقة المقترحة تقدم مقايضة مثالية بين الدقة والتكلفة. فهي تعالج قصور التقليص المنتظم للأبعاد، الذي يفشل في التقاط تغيرات الضغط المحلية تحت ظروف التقليص الشديد، دون تحمل التكلفة الحسابية العالية للإسقاط كامل الدقة أو دورات التحسين التكرارية.

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