طريقة بولتزمان الشبكية الكمومية: ترميز مصفوفة الكثافة للتدفقات ذات الجدران
السياق
تعد الحوسبة الكمومية واعدة بتوفير هائل في الذاكرة لمشاكل الديناميكا المائية، لكن اللاخطية لمعادلات نافير-ستوكس وتعقيد الشروط الحدودية أعاقت التقدم لفترة طويلة. تستبدل طريقة بولتزمان الشبكية التقليدية (LBM) النقل اللاخطي بخطوة انتشار خطية للجسيمات، مما يجعلها جذابة للتنفيذ الكمومي. ومع ذلك، فإن الإصدارات الكمومية الحالية (QLBM) إما تقتصر على حدود دورية أو تتطلب موارد تنمو بشكل أسي. تسد ورقة سو ووانغ ويان هذه الفجوة من خلال اقتراح مخطط بحدود ذات هندسة تعسفية وشروط دخول/خروج. يعتمد العمل على نهج المجموعة لوانغ وآخرين، لكنه يترجمه إلى لغة مصفوفات الكثافة، مما يفتح الطريق أمام النمذجة العملية.
المنهجية
يستخدم المؤلفون ترميز الكثافة: تُكتب احتمالات جميع تكوينات السرعات على جميع عقد الشبكة في قطر مصفوفة الكثافة. وهذا يسمح بتنفيذ المصفوفات العشوائية (مثل خطوة النقل أو التصادم) عبر مؤثرات كراوس، والتي يمكن دائمًا تضمينها في تطور وحدوي مع كيوبتات مساعدة. بالنسبة للجدران، يقدمون خطوة "تبادل المكونات" قبل النقل: إذا كان الجسيم متجهًا نحو الجدار، يتم تبادل مكونه مع المكون المنعكس في نقطة الحدود المجاورة، ثم يعيده النقل العادي تلقائيًا. يتم تعيين شروط الدخول والخروج عبر عملية SWAP متحكم بها مع سجل مساعد مُعد بتوزيع السرعات المطلوب. لا يتطلب المخطط تعديل خطوة النقل ويحافظ على التعقيد الكلي من رتبة O(2^q n^2) بوابة لكل خطوة.
النتائج
تم اختبار الطريقة كميًا على ثلاثة أنواع من التدفقات. بالنسبة لتدفق بويزيلي المخمد في قناة بجدران ثابتة، يتناقص الخطأ النسبي L2 مع زيادة دقة الشبكة (من 32 إلى 256 عقدة) وينمو خطيًا مع الزمن، مما يؤكد صحة شرط الالتصاق الحدودي. في التدفق خلف العتبة العكسية عند أرقام رينولدز 17.8 و35.5 و71، يختلف طول منطقة إعادة التدوير عن نتائج OpenFOAM بنسبة 17.0% و13.0% و10.3% على التوالي. بالنسبة للتدفق حول أسطوانة عند Re=6.6 و43، بلغت الأخطاء النسبية 17.34% و12.53%. أخيرًا، تُظهر محاكاة التدفق حول الحروف "PKU" على شبكة 2048×512 عند Re=64.8 التقاط الدوامات في الفجوات بين الحروف، مما يؤكد قابلية التطبيق على الهندسة المعقدة.
الأهمية
يُظهر العمل لأول مرة أن طريقة بولتزمان الشبكية الكمومية يمكنها مراعاة الجدران الصلبة ذات الشكل التعسفي بشكل صحيح دون تغيير خطوة النقل. وهذا يميزها عن الأساليب المجاورة التي إما تتجاهل الحدود أو تتطلب تعديل المخطط بأكمله. استخدام مصفوفات الكثافة بدلاً من الحالات النقية يزيل القيد الصارم للوحدوية ويفتح الطريق أمام فئة أوسع من العمليات غير الوحدوية في الديناميكا المائية الكمومية. عمليًا، هذا يقرب QLBM من المشاكل الهندسية الحقيقية حيث تكون الحدود ذات الهندسة المعقدة حتمية.
التطور المستقبلي
تبقى الخطوة H، التي تستعيد استقلالية مكونات السرعة بعد التصادمات، هي عنق الزجاجة الرئيسي: فهي تتطلب تكرار المخطط بأكمله عدة مرات وقد تلغي التسارع الكمومي. يخطط المؤلفون للبحث عن تطبيقات مبسطة للخطوة H أو صيغ بديلة بدونها. كما يجب توسيع الطريقة لتشمل الحدود المتحركة والقابلة للتشوه، الأمر الذي سيتطلب طرقًا جديدة لتحديد الشروط الحدودية في مصفوفة الكثافة.
التأثير
قد تؤثر الطريقة على الديناميكا المائية الحسابية، خاصة في مشاكل التدفق حول الأجسام ذات الأشكال المعقدة، والميكروفلويديك، والديناميكا الهوائية. كما تحفز تطوير خوارزميات كمومية للعمليات غير الوحدوية في مجالات أخرى حيث يلزم العمل مع الأنظمة الكمومية المفتوحة.
الخطوات التالية
تشمل الخطوات التالية تحسين الخطوة H لتقليل التكاليف الإضافية واختبارها على أجهزة كمومية حقيقية مع مراعاة الضوضاء. ومن المهم أيضًا دراسة دقة الشروط الحدودية عند أرقام رينولدز أعلى وفي مشاكل ثلاثية الأبعاد.
المسائل المفتوحة الرئيسية
في الديناميكا المائية، تظل مسألة مراعاة الشروط الحدودية بدقة وكفاءة في الخوارزميات الكمومية مفتوحة، خاصة للتدفقات المضطربة ذات أرقام رينولدز العالية. يخطو العمل خطوة نحو الحل، لكنه لا يحل مشكلة الحدود المتحركة ولا يزيل التعقيد الأساسي للخطوة H الذي يحد من التسارع العملي.
لمحاكاة التدفق حول الحروف "PKU" (اختصار جامعة بكين)، تطلب الأمر مخططًا كموميًا بشبكة 2048×512 عقدة و23 كيوبتًا، وهو ما يعادل شبكة تقليدية تضم أكثر من مليون نقطة.
تتردد فكرة المحاكاة الكمومية للسوائل مع قصة ستانيسلاف ليم "المنيع"، حيث تشكل الروبوتات النانوية عقلًا جماعيًا يشبه السائل. هنا، يحاكي الكمبيوتر الكمومي سلوك العديد من الجسيمات التي تخضع لقواعد محلية بسيطة، مما يذكرنا بالخصائص الناشئة في الخيال العلمي.
أرقام رئيسية
- الخطأ النسبي للعتبة العكسية: 10.3% عند Re=71
- الخطأ النسبي للأسطوانة: 12.53% عند Re=43
- عدد الكيوبتات لشبكة 128×128: 23
- عدد العقد لحروف PKU: 2048×512
- رقم رينولدز لحروف PKU: 64.8