A priori Assessment of Tensor-Network Encoding for Isotropic Turbulent Flows
تُثبت هذه الورقة أن نموذج "حالة ناتج ضرب المصفوفات" (MPS) في شبكة التنسور يمكنه ضغط بيانات تدفق الاضطراب المتماثل من المحاكاة العددية المباشرة بفعالية، محققاً دقة إعادة بناء تصل إلى 99.8% باستخدام 5-15% فقط من الذاكرة الأصلية مع الحفاظ بدقة على الخصائص الإحصائية الرئيسية مثل الطاقة الحركية ومعدلات التبديد.
المؤلفون الأصليون:Massen Esmaeili, Hirad Alipanah, Robert Pinkston, Peyman Givi, Daniel Livescu, Andrew J. Daley, Dieter Jaksch, Juan José Mendoza-Arenas
المؤلفون الأصليون: Massen Esmaeili, Hirad Alipanah, Robert Pinkston, Peyman Givi, Daniel Livescu, Andrew J. Daley, Dieter Jaksch, Juan José Mendoza-Arenas
الاضطراب هو الحركة الفوضوية والدوامية للسوائل التي تحيط بنا، من الدخان المتصاعد من سيجارة إلى الهواء المندفع فوق جناح طائرة. إن التنبؤ بكيفية سلوك هذه التدفقات يعد أحد أصعب التحديات في الفيزياء لأن الحركة تحدث على مستويات عديدة ومختلفة في آن واحد. تحتوي الدوامات الكبيرة معظم الطاقة، بينما تكون الدوامات الأصغر هي المكان الذي تتبدد فيه تلك الطاقة في النهاية لتتحول إلى حرارة. وللمحاكاة بدقة على جهاز كمبيوتر، يجب على العلماء تتبع كل دوامة بمفردها، من الضخمة منها وصولاً إلى المجهرية. هذا المتطلب يخلق مشكلة بيانات ذات أبعاد هائلة؛ إذ يمكن لمحاكاة حاسوبية واحدة عالية الجودة لمثل هذا التدفق أن تولد قدراً من المعلومات كفيلاً بملء ذاكرة أقوى الحواسيب الفائقة في العالم، مما يجعل من المستحيل تقريباً تخزين النتائج أو تحليلها أو مشاركتها. ولعقود من الزمن، سعى الباحثون لإيجاد طرق لضغط هذه البيانات دون فقدان الفيزياء الأساسية التي تحكم التدفق.
لقد اختبر فريق من الباحثين الآن نهجاً جديداً لهذه المشكلة، مستعيرين أداة صُممت في الأصل لمجال علمي مختلف تماماً: الفيزياء الكمومية. في العالم الكمومي، يستخدم العلماء إطاراً رياضياً يسمى "شبكات الموتر" (tensor networks) لوصف الأنظمة ذات الأجزاء المتعددة المتفاعلة، مما يجد طريقة لتمثيل الحالات المعقدة دون الحاجة إلى سرد كل الاحتمالات الممكنة. طبق الباحثون هذا الإطار نفسه على دراسة السوائل المضطربة؛ حيث أخذوا مجموعتين من البيانات التفصيلية المولدة حاسوبياً — إحداهما تمثل الحركة الدوامية لسائل، والأخرى تمثل مادة خاملة، مثل الصبغة، يتم خلطها في ذلك السائل — وحاولوا ضغطهما باستخدام نوع محدد من شبكات الموتر يُعرف باسم "حالة ناتج المصفوفة" (matrix product state). كان الهدف هو معرفة ما إذا كان بإمكانهم تقليص حجم البيانات إلى كسر من حجمها الأصلي مع الاستمرار في الحفاظ على الإحصائيات الدقيقة للتدفق، مثل سرعة حركة السائل أو سرعة اختلاط الصبغة.
عملت العملية من خلال تقسيم الشبكة ثلاثية الأبعاد الضخمة للبيانات إلى سلسلة من القطع الأصغر والمتصلة ببعضها البعض. ثم قام الباحثون بتحليل الروابط بين هذه القطع لتحديد الأجزاء التي كانت ضرورية حقاً وتلك التي كانت زائدة عن الحاجة. ومن خلال التخلص من الروابط الأقل أهمية، أنشأوا نسخة مضغوطة أصغر بكيداً من التدفق الأصلي. وعندما قاموا بتوسيع هذه البيانات المضغودة لتعود إلى حجمها الكامل لمقارنتها بالأصل، كانت النتائج دقيقة بشكل مذهل. بالنسبة لسرعة السائل، احتفظت النسخة المضغوطة بنسبة 99.8% من المعلومات الأصلية مع استخدام 5% فقط من مساحة الذاكرة. وبالنسبة لسكالار الخلط (mixing scalar)، وصلت النسخة المضغوطة إلى نفس مستوى الدقة باستخدام 15% من الذاكرة. وهذا يعني أن مجموعة بيانات قد تتطلب عادةً كمية هائلة من التخزين يمكن تقليصها إلى حجم ملف صغير دون فقدان القدرة على رؤية سلوك التدفق.
لم يكتفِ الباحثون بالتحقق مما إذا كانت الصور تبدو متشابهة فحسب؛ بل اختبروا الفيزياء الكامنة في البيانات بصرامة. فقد فحصوا طاقة التدفق، ومعدل تبددها، والأنماط الإحصائية لكيفية حركة السائل عند سرعات مختلفة. وحتى عند مستويات الضغط العالية، تم استرداد إجمالي طاقة التدفق بخطأ يقل عن 0.2%. كما كان متوسط معدل فقدان الطاقة بسبب الاحتكاك قريباً جداً من القيمة الأصلية، حيث ظل ضمن حدود 10% تقريباً من القيمة الحقيقية. ويعد هذا إنجازاً كبيراً لأن معدل فقدان الطاقة يعتمد على أدق التفاصيل وأكثرها فوضوية في التدفق، وهي التفاصيل التي عادة ما تكون أول ما يختفي عند ضغط البيانات. وجد الباحثون أن البيانات المضغوطة نجحت في الحفاظ على نبضات الحركة المعقدة والمتقطعة التي تميز الاضطراب، مما أثبت أن الطريقة لم تكتفِ بتنعيم الحواف الخشنة بل حافظت على الطبيعة المتعرجة والوعرة للتدفق.
ومع ذلك، كشفت الدراسة أيضاً عن المواضع التي تكمن فيها حدود هذا الضغط. فبينما تم الحفاظ على الهياكل الكبيرة للتدفق بشكل رائع، تطلبت الأحداث الأكثر تطرفاً ودقة في التفاصيل أعلى مستوى من التفاصيل ليتم التقاطها بدقة. فعندما استخدم الباحثون ضغطاً أكثر حدة أدى إلى إزالة المزيد من البيانات، بدأت الإحصائيات المتعلقة بالمقاييس الأصغر، مثل التدرجات الحادة حيث يتغير اتجاه السائل بسرعة، تظهر أخطاء أكبر. يشير هذا إلى أنه بينما تعد الطريقة فعالة للغاية لتخزين وتحليل الجزء الأكبر من التدفق، فإن التقاط اللحظات الأكثر ندرة وكثافة في الاضطراب يتطلب الاحتفاظ بمزيد من البيانات الأصلية. ووجد الباحثون أن الطريقة التي رتبوا بها البيانات قبل ضغطها كانت مهمة للغاية؛ فمن خلال دمج المعلومات من اتجاهات مختلفة، تمكنوا من التقاط الارتباطات بين الحركات الدوامية بشكل أكثر فعالية مما لو قاموا بتجميعها بشكل منفصل.
يُظهر هذا العمل أن شبكات الموتر، وهي أداة من ميكانيكا الكم، يمكن أن تعمل كوسيلة جديدة قوية للتعامل مع مجموعات البيانات الضخمة الناتجة عن ديناميكا السوائل. إنها توفر مساراً للمضي قدماً في تخزين وتحليل التدفقات المضطربة المعقدة دون الحاجة إلى سعة التخزين الهائلة المطلوبة حالياً. وتشير النتائج إلى إمكانية استخدام هذه الطريقة لإنشاء نماذج مختزلة الرتبة تكون مدمجة ومخلصة للفيزياء في آن واحد، مما يسمح للعلماء بدراسة التدفقات المعقدة بكفاءة أكبر. وبينما أُجريت الاختبارات الحالية على تدفقات مثالية ومنتظمة، فإن نجاح الطريقة يوفر أساساً قوياً لاستكشاف استخدامها في سيناريوهات أكثر تعقيداً وواقعية، مثل التدفقات التي تحتوي على جدران، أو تفاعلات كيميائية، أو أطوار متعددة. وتؤكد الدراسة أن الفيزياء الأساسية للاضطراب يمكن تشفيرها في كمية صغيرة بشكل مفاجئ من المعلومات، مما يفتح الباب أمام عصر جديد من ديناميكا السوائل الحاسوبية حيث تكون البيانات قابلة للإدارة، بينما يظل العلم دقيقاً.
بيان المشكلة تعد المحاكاة العددية المباشرة (DNS) للتدفقات المضطربة مكلفة حوسبياً بسبب النمو الأسي في متطلبات الذاكرة والمعالجة مع زيادة رقم رينولدز. إن حل كامل طيف المقاييس، من الدوامات الحاملة للطاقة إلى المقاييس التبديدية، يتطلب شبكات دقيقة جداً وخطوات زمنية صغيرة، مما يتجاوز غالباً قدرات الحواسيب الفائقة الحالية. وبينما توجد تقنيات النمذجة ذات الرتبة المنخفضة (ROM)، إلا أنها تعتمد عادةً على التفكيك النمطي (مثل POD) أو طرق قائمة على الإسقاط التي قد لا تستطيع الحفاظ بشكل كامل على الارتباطات متعددة المقاييس المتأصلة في الاضطراب المتناحي. وبالتزامن مع ذلك، أثبتت طرق شبكة الموتر (TN)، وتحديداً حالات المصفوفة المترابطة (MPS)، فاعلية عالية في ضغط البيانات عالية الأبعاد في فيزياء الأنوية الكثيرة الأجسام الكمومية والمجالات الأخرى، ومع ذلك، لا يزال تطبيقها لتوصيف التدفقات المضطربة المتطورة بالكامل غير مستكشف بشكل كبير. هناك نقص في التقييمات المنهجية "القبلية" التي تحدد كيفية تأثير الضغط القائم على MPS على إحصاءات الاضطراب عبر المقاييس المكانية.
المنهجية تستخدم هذه الدراسة نموذج حالة المصفوفة المترابطة (MPS) كإطار عمل للنمذجة ذات الرتبة المنخفضة لبناء تمثيلات مختزلة لبيانات التدفق المضطرب متناحي الخواص. تتضمن المنهجية ما يلي:
ترميز البيانات: تم استخدام مجموعتي بيانات DNS: مجال سرعة هيدروديناميكي غير قابل للانضغاط ثلاثي الأبعاد (Reλ≈433) ومجال سكولار (كمية قياسية) فيكي (Fickian) محفوظ يتم نقله بواسطة تدفق مشابه (Reλ≈180). كلا المجالين معرفان على شبكة 10243.
إعادة تشكيل الموتر: تم إعادة تشكيل المجالات ثلاثية الأبعاد إلى موترات غير مهيكلة من الرتبة N=log2(n) باستخدام ترميز قائم على الكيوبت (d=2)، حيث يقابل كل مؤشر فيزيائي طولاً مكانياً محدداً.
ترتيب المؤشرات: لتحسين التقاط الارتباطات، تمت إعادة ترتيب مؤشرات الموتر المكانية. تقارن الدراسة ترتيبات "المقسم" (split)، و"الدرجي" (staircase)، و"المتداخل" (interleaved). تم اختيار الترتيب المتداخل (x1,y1,z1,x2,y2,z2,…) لأنه يضع المؤشرات ذات الأطوال المتقاربة من الاتجاهات المكانية المختلفة جنباً إلى جنب، مما يؤدي إلى تركيز الارتباطات السائدة بين الموترات ضمن مؤشرات الروابط قصيرة المدى.
الضغط عبر SVD: يتم تحليل الموتر إلى سلسلة من موترات المواقع المتصلة بمؤشرات روابط عبر تسلسل من تفكيك القيم المفردة (SVD). يتم التخلص من القيم المفردة الصغيرة بناءً على عتبة قطع محددة من قبل المستخدم (10−C، حيث C∈{2,3,4,5}). يعد هذا الاختزال أمثل في سياق معيار L2 (فروبينيوس)، حيث يحد من الخطأ بمجموع مربعات القيم المفردة المستبعدة.
إعادة البناء والتحقق: يتم دمج MPS المختزل مرة أخرى إلى الشبكة الكاملة 10243. وتتم مقارنة المجالات المعاد بناؤها مع بيانات DNS الأصلية باستخدام مقايسات عدم الوفاء (infidelity)، ونسب الضغط، ومجموعة شاملة من المقاييس الإحصائية.
المساهمات الرئيسية
تقييم منهجي قبلي: تقدم الورقة تقييماً صارماً لترميز MPS على بيانات DNS عالية الدقة، متجاوزة الدراسات المحدودة السابقة لفحص مجموعة واسعة من الإحصاءات من الرتب المنخفضة والعالية.
تحسين ترتيب المؤشرات: تثبت أن الترتيب المتداخل يتفوق بشكل كبير على الترتيب المقسم في الاضطراب المتناحي من خلال مواءمة بنية سلسلة MPS مع بنية الارتباط الفيزيائي للتدفق (حيث ترتبط المقاييس ذات الحجم المتماثل عبر x,y,z إحصائياً).
تحديد المقايضات: ترسم الدراسة بشكل منهجي المقايضة بين نسبة الضغط (تقليل التخزين) والدقة عبر مختلف الكميات الإحصائية، بما في ذلك أطياف الطاقة، ودوال الكثافة الاحتمالية (PDFs)، والدوال الهيكلية، والمقاييس المعتمدة على التدرج.
النتائج
كفاءة الضغط: يحقق تمثيل MPS تخفيضات كبيرة في التخزين. بالنسبة لمجال السرعة، يحقق القطع عند C=4 نسبة ضغط (CR) تبلغ حوالي 22 (الاحتفاظ بـ ~5% من الذاكرة الأصلية) مع عدم وفاء بنسبة 0.18%. بالنسبة للمجال السكولار، يحقق القطع عند C=4 نسبة ضغط تبلغ حوالي 6.5 (الاحتفاظ بـ ~15% من الذاكرة الأصلية) مع عدم وفاء بنسبة 0.18%. أما عند أعلى إعداد للدقة للمجال السكولار (C=5)، تنخفض نسبة الضغط إلى حوالي 1.7، مع الاحتفاظ بنحو 59% من الذاكرة الأصلية، وعدم وفاء بنسبة 0.02%. يتطلب المجال السكولار عموماً استخدام ذاكرة أعلى من مجال السرعة لتحقيق دقة مماثلة بسبب تدرجاته الأكثر حدة.
الدقة الإحصائية:
الإحصاءات منخفضة الرتبة: يتم استعادة إجمالي الطاقة الحركية وطاقة السكولار بأخطاء نسبية أقل من 0.2% عند C=4.
الإحصاءات المعتمدة على التدرج: تظهر الكميات التي تعتمد على المشتقات المكانية، مثل التبدد المتوسط وتبدد السكولار، حساسية أعلى للاختزال. عند C=4، تبلغ أخطاء التبدد المتوسط حوالي 4% للسرعة و10% للسكولار. تنخفض هذه الأخطاء إلى ∼2.6% و∼1.8% على التوالي عند C=5.
الميزات الطيفية والهيكلية: تتبع أطياف الطاقة المعاد بناؤها مقياس كولموغوروف k−5/3 في النطاق العطالي بدقة حتى النطاق التبديلي المعتمد على القطع. كما يتم الحفاظ جيداً على توزيع "قطرة الدموع" المميز لثوابت Q−R (ثوابت موتر تدرج السرعة) وقانون كولموغوروف 4/5 للدوال الهيكلية من الدرجة الثالثة عند C=4 وC=5.
دوال الكثافة الاحتمالية (PDFs): يلتقط MPS التوزيع القريب من التوزيع الطبيعي لمكونات السرعة والتوزيعات ذات الذيول الثقيلة لتدرجات السرعة والتبدد. ومع ذلك، فإن الأحداث المتطرفة (ذيول الـ PDFs) يتم تمثيلها بأقل من قيمتها الحقيقية عند القطوع المنخفضة (C=3) وتتطلب C=4 أو C=5 لإعادة إنتاج دقيق.
تحليل أطياف شميدت: تكشف الدراسة أن طيف شميدت يظهر بنية هرمية مع عدم تماثل؛ حيث تتطلب أبعاد الروابط الأكبر عند طرف "المقياس الدقيق" لسلسلة MPS التقاط الارتباطات واسعة المدى للتقلبات صغيرة المقياس، بينما يتم التقاط الهياكل كبيرة المقياس بعدد أقل من القيم المفردة.
الأهمية والادعاءات يدعي المؤلفون أن هذه النتائج تدعم ملاءمة MPS للتحليل المختزل القابل للتوسع لمجموعات البيانات المعقدة للاضطراب. يوضح العمل أن طرق TN يمكنها الحفاظ ليس فقط على الهياكل الحاملة للطاقة في المقاييس الكبيرة، بل أيضاً على المقاييس الأدق المرتبطة بالخلط والاضطراب (intermittency)، بشرما كانت عتبة القطع منخفضة بما يكفي. تسمح الطبيعة التكيفية لبناء MPS — من خلال الاحتفاظ بالمزيد من القيم المفردة حيث تكون الارتباطات بين الموترات أقوى — بتوازن في الحفاظ على فيزياء التدفق عبر المقاييس.
تخلص الورقة إلى أنه بينما يعد MPS فعالاً للضغط والتحليل القبلي، فإن تطبيقه على المحاكاة المتطورة زمنياً والتدفقات الأكثر تعقيداً (مثل التدفقات غير المتناحية، أو التدفقات المرتبطة بالجدران، أو التدفقات التفاعلية) يتطلب مزيداً من الاستكشاف. وتقترح أن العمل المستقبلي يجب أن يبحث في بنيات TN أخرى (مثل شبكات الموتر الشجرية - Tree Tensor Networks، أو PEPS) التي قد تلتقط بشكل أفضل هياكل الارتباط المحددة للاضطراب ثلاثي الأبعاد، واستكشاف التوازي بين اختزال TN ونمذجة المقياس الفرعي (subgrid-scale modeling) في محاكاة LES.