⚛️ quantum physics

An HHL-Based Quantum-Classical Solver for the Incompressible Navier-Stokes Equations with Approximate QST

تقدم هذه الورقة حلاً هجيناً كمياً-كلاسيكياً يدمج خوارزمية "هارو-هاسيديم-لوي" (HHL) مع تقنية "تشيبيشيف" لتقصي الحالة الكمية التقريبي، وذلك لحل معادلات نافييه-ستوكس غير القابلة للانضغاط بكفاءة، حيث تم التحقق من نجاح هذا النهج من خلال محاكاة دقيقة لتدفقات تجويف الغطاء وتدفقات دوامة "تايلور-غرين" باستخدام إطار عمل "كيسكيت" (Qiskit) التابع لشركة IBM.

المؤلفون الأصليون: Moshe Inger, Steven Frankel

نُشر 2026-03-20

📖 4 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Moshe Inger, Steven Frankel

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

الصورة الكبيرة: تعاون بين الكم والتقليدي

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

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

المشكلة: على جهاز كمبيوتر عادي، يكون حل هذا اللغز مثل محاولة عد كل حبة رمل على الشاطئ واحدة تلو الأخرى. هذا يستهلك 90% من وقت الكمبيوتر ويبطئ كل شيء.

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

الشخصيات الرئيسية

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

كيف اختبروا ذلك: روتينات "الجمباز" المزدوجة

لإثبات أن فريقهم يعمل، اختبروه على تحديين مشهورين في ديناميكا السوائل:

1. تجربة التجويف مدفوع الغطاء (صندوق الجيلي)

تخيل صندوقاً مربعاً مليئاً بالماء. يتحرك الغطاء العلوي جهة اليمين، مما يسحب الماء معه، ويخلق دوامة ضخمة.

الاختبار: قاموا بمحاكاة ذلك على شبكة 16×16.
النتيجة: نجح الفريق الهجين في إعادة إنشاء الدوامة الدوارة. كانت إجابة الكمبيوتر الكمي دقيقة بنسبة 92% تقريباً مقارنة بالإجابة الكلاسيكية المثالية. كانت الأخطاء تتركز في الزوايا حيث يصبح الماء ساكناً جداً (بالقرب من الصفر)، وهو أمر يصعب قياسه بدقة.

2. دوامة تايلور-غرين (الدونات المثالية الدوارة)

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

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

العقبات التي واجهوها

على الرغم من نجاح التجربة، إلا أن الفريق واجه بعض العثرات في الطريق:

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

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

هذه الورقة هي إثبات لمفهوم. إنها تشبه أول رحلة للأخوين رايت: لم تكن تحمل ركاباً عبر المحيط، لكنها أثبتت أن الآلة يمكنها الطيران.

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

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

إليك ملخص تقني مفصل للورقة البحثية بعنوان "محلل هجين (كمي-كلاسيكي) يعتمد على خوارزمية HHL لمعادلات نافييه-ستوكس غير القابلة للانضغاط مع تصوير الحالة الكمية التقريبي (QST)" من إعداد موشيه إنغر وستيفن فرانكل.

1. بيان المشكلة

يعد الحل العددي لمعادلات نافييه-ستوكس غير القابلة للانضغاط حجر الزاوية في ديناميكا الموائع الحسابية (CFD). ومع ذلك، تتطلب طرق الإسقاط القياسية حل معادلة بواسون للضغط عند كل خطوة زمنية لفرض شرط عدم الانضغاط ( $\nabla \cdot u = 0$ ).

عنق الزجاجة: يؤدي تمثيل معادلة بواسون رقمياً إلى نظام ضخم من المعادلات الجبرية الخطية المتفرقة ( $Ax = b$ ). ويستهلك حل هذا النظام عادةً ما يصل إلى 90% من إجمالي وقت التنفيذ في محاكاة CFD الكلاسيكية.
الفرصة الكمية: توفر خوارزمية هارو-هسيديم-للويد (HHL) تسارعاً نظرياً أسياً لحل الأنظمة الخطية المتفرقة. ومع ذلك، هناك عقبتان رئيسيتان تمنعان تطبيقها العملي في CFD:
1. محدودية القراءة (Readout Limitation): تنتج خوارزمية HHL حالة كمية $|x\rangle$ تتناسب مع الحل. ويتطلب استخراج المتجه الكامل (متجه الحالة) عدداً أسياً من القياسات، مما يلغي ميزة التسريع.
2. تحضير الحالة (State Preparation): يمثل تحميل البيانات الكلاسيكية في حالة كمية (ترميز السعة) تحدياً مفتوحاً.
3. القيود الهجينة: بينما تقوم HHL بترميز $2^n$ من نقاط البيانات في $n$ من الكيوبتات، فإن الحاجة لنقل البيانات ذهاباً وإياباً بين المعالجات الكمية والكلاسيكية تحد حالياً من الميزة، حيث لا تزال هناك حاجة لـ $2^n$ بت كلاسيكي للمخرجات.

2. المنهجية

يقترح المؤلفون محللاً هجيناً (كمياً-كلاسيكياً) تكرارياً يدمج خوارزمية HHL في مخطط إسقاط كلاسيكي لمعادلات نافييه-ستوكس ثنائية الأبعاد غير القابلة للانضغاط.

أ. الإطار الكلاسيكي

المعادلات الحاكمة: يستخدم المحلل طريقة إسقاط الخطوات المنفصلة، حيث يتم حساب سرعة وسيطة $u^*$ بشكل صريح، يتبعها خطوة تصحيح الضغط لفرض عدم الانضغاط.
التمثيل الرقمي (Discretization): يتم تمثيل النطاق باستخدام مخطط فرق مركزي من الدرجة الثانية على شبكة ديكارتية ثنائية الأبعاد. يتم تحويل الشبكة ثنائية الأبعاد إلى متجه أحادي البعد لتكوين مصفوفة متفرقة، ثلاثية الأقطار الكتلية $A$ لمعادلة بواسون ( $\nabla^2 p = \rho/\Delta t \nabla \cdot u^*$ ).
التحويل المنحني (Curvilinear Mapping): للتعامل مع الطبقات الحدودية (تحديداً في تدفق التجويف ذي الغطاء المتحرك)، تُستخدم دالة تمدد زائدي لتجميع نقاط الشبكة بالقرب من الحدود، مما يحسن الدقة حيث تكون التدرجات عالية.

ب. الإجراء الفرعي الكمي (HHL)

الخوارزمية: تُستخدم خوارزمية HHL لحل $Ap = b$ ، حيث $b$ هو تباعد السرعة الوسيطة و $p$ هو الضغط.
التنفيذ:
- الترميز: يتم تفكيك المصفوفة $A$ إلى سلاسل باولي ( $A = \sum c_k P_k$ ).
- محاكاة هاميلتون (Hamiltonian Simulation): يتم تقريب التطور الزمني $e^{-iAt}$ باستخدام صيغة "لي-تروتر" للمنتج من الدرجة الأولى مع $R=150$ خطوة.
- عكس القيم الذاتية (Eigenvalue Inversion): يتم عبر دوران متحكم فيه على كيوبت مساعد باستخدام دالة ExactReciprocal من مكتبة Qiskit.
- تحضير الحالة: يفترض المؤلفون وجود "أوراكل" (Oracle) مثالي لتحميل المتجه المدخل $|b\rangle$ ، مع الإقرار بأن هذا يمثل عائقاً حالياً.

ج. تصوير الحالة الكمية التقريبي (QST)

للتغلب على عقبة القراءة دون إعادة بناء الحالة بالكامل، استخدم المؤلفون نهج QST القائم على كثيرات حدود تشيبيشيف:

الإسقاط الطيفي: بدلاً من قياس كل سعة، يتم إسقاط الحالة الكمية $|p\rangle$ على قاعدة مختزلة من $m$ من كثيرات حدود تشيبيشيف.
اختبار هادامارد (Hadamard Test): يتم تقدير معاملات كثيرات الحدود باستخدام اختبار هادامارد.
الكفاءة: يقلل هذا من تعقيد القراءة من التعقيد الأسي إلى $O(md) $، حيث$ m $هو عدد كثيرات الحدود و$ d$ هو البعد.
المعايرة: بما أن HHL تنتج حالة معيرة (normalized)، يقوم المؤلفون بإجراء معايرة طاقة عند كل خطوة زمنية عن طريق مطابقة معيار $L_2$ للضغط الكمي مع مرجع كلاسيكي لاستعادة المقدار الفيزيائي.

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

محلل هجين متكامل (نهاية إلى نهاية): أول عرض لمحلل متكامل (كمي-كلاسيكي) يحل معادلات نافييه-ستوكس، حيث يستبدل الإجراء الفرعي لـ HML محل حل بواسون الكلاسيكي.
دمج QST التقريبي: تنفيذ ناجح لـ QST القائم على تشيبيشيف لاستخراج بيانات ديناميكا الموائع من حالة كمية، متجاوزاً الحاجة إلى قياسات أسية.
التحقق من المعايير (Benchmark Validation): التحقق من صحة المحلل مقابل معايير CFD القياسية (تجويف الغطاء المتحرك وتدفق تايلور-غرين) باستخدام إطار عمل IBM Qiskit.
تحليل رسم الشبكة: استقصاء المقايضة بين تمدد الشبكة المنحنية (الذي يساعد دقة QST) وحالة المصفوفة (التي تؤثر على دقة HHL).

4. النتائج

تم اختبار المحلل على مشكلتين معياريتين بشبكة $16 \times 16$ و 8 كيوبتات ساعة ( $n_c=8$ ).

محلل معادلة بواسون:
- حقق متوسط خطأ نسبي (ARE) بنسبة ~2.3% (في 1D) و ~2.5% (في 2D) مقارنة بالحلول الكلاسيكية.
- كانت الدقة حساسة لعدد كيوبتات الساعة، حيث لوحظ تحسن حاد عند $n_c=8$ .
تدفق التجويف ذي الغطاء المتحرك ( $Re=100$ ):
- مجال السرعة: نجح المحلل الهجين في رصد ديناميكيات الدوامة العالمية بمتوسط خطأ نسبي (ARE) قدره 8.17%.
- تدرج الضغط: لوحظت أخطاء عالية (تصل إلى 416%) عند الزوايا السفلية للتجويف حيث يقترب الضغط من الصفر. يُعزى ذلك إلى صعوبة قدرة تقريب تشيبيشيف على التقاط القيم العالية والقيم القريبة من الصفر في آن واحد باستخدام قاعدة محدودة.
- التقارب: تقاربت ملفات السرعة على الخطوط المركزية نحو معيار غيا (Ghia benchmark).
تدفق تايلور-غرين (TGV):
- استخدم هذا الاختبار حلاً تحليلياً دقيقاً للتحقق من الصحة.
- النتائج: حقق المحلل دقة عالية بمتوسط خطأ نسبي (ARE) بنسبة ~1% للسرعة و ~4% للضغط.
- تؤكد النتائج أن النهج الهجين الكمي يمكنه نمذجة الاضمحلال الأسي للدوامات الناتج عن التبديد اللزوجي بدقة.

5. الأهمية والتوجهات المستقبلية

إثبات المفهوم: يثبت هذا العمل أن الإجراءات الكمية الفرعية يمكن دمجها بنجاح في سير عمل CFD التكراري دون التسبب في تباعد كارثي، حتى مع قيود المحاكاة الحالية.
معالجة عقبة القراءة: يوفر استخدام QST القائم على تشيبيشيف مساراً قابلاً للتطبيق لاستخراج بيانات فيزيائية مفيدة من الحالات الكمية، وهي خطوة حاسمة لتحقيق التفوق الكمي العملي في ديناميكا الموائع.
القيود والتحديات:
- تحضير الحالة: الاعتماد على "أوراكل" مثالي لتحميل المدخلات يظل عائقاً نظرياً.
- التوسع (Scaling): تقتصر المحاكاة الحالية على شبكات صغيرة ( $16 \times 16$ ) بسبب قيود الكيوبتات.
- الخطأ عند الحدود: يواجه منهج QST صعوبة في المناطق ذات القيم القريبة من الصفر عندما توجد قيم ذات مقدار عالٍ في نفس النطاق.
العمل المستقبلي:
- التوسع إلى معادلات نافييه-ستوكس ثلاثية الأبعاد (3D).
- استبدال تحضير الحالة المثالي بخوارزميات أكثر كفاءة (مثل شبكات التنسور أو HHL المعدلة).
- الانتقال من HHL القائم على التسامح مع الأخطاء إلى حلول الخطية الكمية التباينية (VQLS) للأجهزة الكمية في عصر NISQ.
- تطوير مخططات كمية كاملة تستبدل خطوات الإدفق-الانتشار الكلاسيكية بالتطور الهاميلتوني.

في الختام، تضع هذه الورقة إطاراً تأسيسياً لـ ديناميكا الموائع الكمية (Quantum CFD)، حيث تثبت أن البنى الهجينة يمكنها حل مشكلات معقدة في ديناميكا الموائع، وتوفر مساراً قوياً لدمج الحوسبة الكمية في عمليات محاكاة أرق دوامات (High-Reynolds number) مستقبلاً.