Adaptive Multiphysics Coupling for Hyperbolic Systems

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

يقدم هذا البحث طريقة ذكية للقيام بذلك بالضبط باستخدام كود حاسوبي يسمى Trixi.jl. فكر في الأمر كأنه "مايسترو" بارع لفرقة موسيقية حيث يعزف بعض الموسيقيين موسيقى "الهيفي ميتال" (فيزياء معقدة) بينما يعزف آخرون موسيقى "اللالاي" أو المقطوعات الهادئة (فيزياء بسيطة). لقد بنى المؤلفون، S. Candelaresi وE. Faulhaber وM. Schlottke-Lakemper، نظاماً يسمح لهؤلاء "الموسيقيين" المختلفين بالتواصل فيما بينهم عند الحدود حيث تلتقي عوالمهم. ويسمون هذا "الاقتران متعدد الفيزياء التكيفي" (adaptive multiphysics coupling). إن الخدعة السحرية هي أن النظام لا يسمح لهم بالتواصل فحسب؛ بل يمكنه أيضاً تغيير القواعد أثناء العمل. فإذا هبت عاصفة مغناطيسية فجأة نحو الحي الهادئ، يقوم النظام فوراً بترقية ذلك المنطقة إلى النموذج عالي القدرة. وإذا رحلت العاصفة، يقوم بخفض المستوى في تلك المنطقة للعودة إلى النموذج البسيط. يتيح هذا للعلماء تشغيل عمليات محاكاة أسرع بكثير من ذي قبل، مما يوفر وقتاً هائلاً من وقت الحاسوب مع الحصول على نتائج دقيقة.

المشكلة: المقاس الواحد لا يناسب الجميع

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

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

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

الحل: مترجم عالمي وكائن متقلب الشكل

حل الفريق في جامعة أوغسبورغ وجامعة كولونيا هذه المشكلة عبر إنشاء نظام يعمل كمترجم عالمي وكائن متقلب الشكل.

المترجم العالمي (دوال الاقتران)
تخيل جارين، أحدهما يتحدث الفرنسية والآخر يتحدث اليابانية. يريدان تبادل التفاح. يقول الجار الفرنسي: "هذه تفاحة"، لكن الجار الياباني لا يفهم كلمة "تفاحة". إنهما بحاجة إلى مترجم. في هذا البحث، "المترجم" هو مجموعة من الدوال التي يحددها المستخدم. تأخذ هذه الدوال المتغيرات من جانب واحد (مثل الكثافة والسرعة) وتحولها إلى المتغيرات التي يفهمها الجانب الآخر.

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

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

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

كيف يعمل الأمر في الممارسة العملية

لإثبات نجاح ذلك، أجرى المؤلفون عدة اختبارات باستخدام الكود الخاص بهم، Trixi.jl.

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

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

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

النتائج: السرعة دون تضحية

السؤال الكبير هو: هل يوفر هذا الوقت؟ الإجابة هي نعم قاطعة.

في اختبارهم، قارنوا النظام التكيفي مقابل تشغيل نموذج MHD المعقد للـنطاق بأكمله.

• تشغيل النموذج المعقد الكامل: استغرق 133.5 ميلي ثانية لكل خطوة زمنية.
• التشغيل التكيفي: استغرق 44.9 ميلي ثانية لكل خطوة زمنية.

هذا يعني أن الطريقة التكيفية كانت أسرع بنحو 3 مرات لكل خطوة واحدة من خطوات المحاكاة. وفي كامل عملية التشغيل، انتهت المحاكاة التكيفية في 756.4 ثانية، بينما استغرق نموذج الـ MHD الكامل 1896 ثانية. وهذا يعني تسريعاً إجمالياً قدره 2.5 مرة في الوقت الكلي.

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

لماذا يهم هذا الأمر

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

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