Unification of The Navier-Stokes Equation for Fluid Mechanics and Aeroacoustics
تقترح هذه الورقة صياغة موحدة لمعادلة نافييه-ستوكس لميكانيكا الموائع والصوتيات الهوائية من خلال تقديم لزوجة كينماتيكية معممة ومنهجية جديدة لتقدير معامل اللزوجة الثاني، مما يلغي الحاجة إلى فرضية ستوكس.
الصوت هو همس من الضغط يتحرك عبر الهواء، تموج من الانضغاط والتمدد يمكن لآذاننا رصده حتى عندما يكون الاهتزاز أصغر بمليار مرة من وزن الغلاف الجوي الذي يضغط علينا. لفهم كيفية انتقال هذه التموجات، اعتمد العلماء لفترة طويلة على مجموعة من القواعد تُعرف بمعادلات نافييه-ستوكس، والتي تصف كيفية تحرك السوائل مثل الهواء والماء. تأخذ هذه القواعد في الاعتبار كيفية مقاومة السائل للتدفق، وهي خاصية تسمى اللزوجة، وكيفية تخزينه للطاقة. ومع ذلك، وعلى مدار قرنين تقريبًا، كان هناك افتراض محدد مدمج في هذه القواعد للسوائل القابلة للانضغاط مثل الهواء: الفكرة القائلة بأن مقاومة السائل للانضغاط والتمدد متوازنة تمامًا مع مقاومته للانزلاق فوق نفسه. هذا الافتراض، المعروف باسم فرضية ستوكس، يخبر المعادلات فعليًا بتجاهل الطاقة المفقودة عند انضغاط السائل وتمدده، ومعاملتها كما لو كانت صفرًا. وبينما نجح هذا التبسيط في العديد من المشكلات الهندسية، إلا أنه يخلق تناقضًا رياضيًا عند تطبيقه على الموجات الصوتية، والتي تتعلق جوهريًا بعملية الانضغاط والتمدد تلك.
لقد أعاد باحث مراجعة هذه القواعد التأسيسية، جازمًا بأن الافتراض القديم ليس مجرد تبسيط غير ضروري، بل هو استحالة فيزيائية تؤدي إلى تنبؤات غير مستقرة وخاطئة. ومن خلال تجريد المعادلات من الافتراض القائل بأن مقاومة السائل للانضغاط هي صفر، استنتج نسخة جديدة من المعادلات الحاكمة التي تعامل فقدان الطاقة أثناء الانضغاط كقوة حقيقية قابلة للقياس. وفي عمله، قدم مفهومًا جديدًا يسميه "اللزوجة الكينماتيكية المعممة". فكر في هذا كقياس واحد مدمج لمدى فقدان السائل للطاقة سواء عند انزلاقه أو عند انضغاطه. هذا القياس الجديد يحل محل النسخة القديمة غير المكتملة التي تجاهلت تأثير الانضغاط. ووجد الباحث أنه بدون هذا التصحيح، تتنبأ المعادلات بأن الموجات الصوتية يجب أن تسلك سلوكيات تخالف الفيزياء، مما يشير أساسًا إلى أن السائل سيكتسب طاقة من العدم بدلاً من فقدانها بسبب الاحتكاك والحرارة.
ولاختبار هذا الإطار الجديد، توجه الباحث إلى مشكلة فيزيائية كلاسيكية تُعرف باسم "عدم استقرار رايلي-تايلور". يحدث هذا عندما يستقر سائل ثقيل فوق سائل أخف، وهو وضع غير مستقر بطبيعته. وفي دراسته، قام بمحاكاة مكعب من الهواء حيث وُضع هواء بارد وكثيف فوق هواء دافئ وأخف. وعندما تمت إزالة الحاجز الذي يفصل بينهما فجأة، بدأ الهواء الثقيل في الهبوط وبدأ الهواء الخفيف في الصعود. وما لاحظه الباحث في عمليات المحاكاة الحاسوبية عالية الدقة هو أنه قبل أن تبدأ السوائل في الاختلاط في دوامة فوضوية، أدى التحرر المفاجئ إلى خلق نبضة ضغط أحادية البعد انطلقت من الواجهة. انتقلت هذه النبضة عبر الهواء مثل موجة، حاملة طيفًا من الترددات من المنخفضة جدًا إلى العالية جدًا. وكان سلوك هذه النبضة هو المفتاح لفك اللغز؛ فمن خلال مراقبة كيفية انتشار النبضة وفقدانها للطاقة أثناء انتقالها، استطاع الباحث حساب القيمة الدقيقة للزوجة المعممة الجديدة.
أكدت نتائج هذه المحاكاة أن الافتراض القديم كان بالفعل معيبًا. فعندما استخدم الباحث الطريقة التقليدية التي تتجاهل الطاقة المفقودة بسبب الانضغاط، فشل النموذج في رصد الحركات المتموجة الصغيرة التي تظهر في بداية عملية الاختلاط. كما قلل من تقدير سرعة نمو طبقات الاختلاط. ومع ذلك، عندما استخدم المعادلات الجديدة مع اللزوجة المعممة، طابقت عمليات المحاكاة الواقع الفيزيائي لفقدان الطاقة بدقة أكبر بكثير. وأشار الباحث إلى أن فقدان الطاقة في هذه الموجات الصوتية لا يعود فقط إلى انزلاق السائل فوق نفسه، بل يعود بشكل كبير إلى انضغاط السائل وتمدده. وهذا التمييز أمر بالغ الأهمية لأن هذا يعني أنه من أجل تنبؤات دقيقة بانتشار الصوت واختلاط السوائل، يجب أن تأخذ المعادلات في الاعتبار معامل اللزوجة الثاني، وهي القيمة التي تم اعتبارها صفرًا فعليًا لأكثر من قرن.
يخلص البحث إلى أن المسار المستقبلي يتطلب نهجًا موحدًا يعامل ميكانيكا السوائل والصوتيات الهوائية كنظام واحد متسق. ويقترح الباحث أنه يمكن تقدير اللزوجة الكينماتيكية المعممة الجديدة من خلال تجارب دقيقة أو طرق عددية متقدمة، مثل تتبع كيفية انتقال نبضة صوتية محددة عبر وسط ساكن. وبينما لا تزال القيمة الدقيقة لهذه اللزوجة للسوائل الشائعة مثل الهواء والماء تحت ظروف مختلفة بحاجة إلى تحديد من خلال مزيد من الاختبارات، فإن الإطار النظري بات جاهزًا الآن. وتثبت الدراسة أنه من خلال إزالة افتراض عفا عليه الزمن، يمكن للعلماء أخيرًا وصف كيفية تفاعل الصوت وحركة السوائل بمستوى من الدقة كان مستحيلاً في السابق، مما يفتح الباب أمام نماذج أكثر دقة لكل شيء، بدءًا من الصوتيات الجوية وصولاً إلى الاختلاط المعقد للسوائل في البيئات القاسية.
بيان المشكلة تتناول الورقة عدم اتساق جوهري في الصيغة الحالية لمعادلات نافير-ستوكس (NSE) عند تطبيقها على كل من ميكانيكا الموائع والصوتيات الهوائية. وتنتقد الورقة الاعتماد على فرضية ستوكس، التي تفترض أن اللزوجة الحجمية (μb) تتلاشى في التدفقات الانضغاطية (μb=λ+32μ=0). ويجادل المؤلف بأن هذه الفرضية تؤدي إلى عواقب غير فيزيائية:
إنها تجبر معامل اللزوجة الثاني (λ) أن يكون سالباً (λ=−2μ/3)، مما يجعله مضاداً للانتشار (anti-diffusive). وهذا يتناقض مع الواقع الفيزيائي الذي يقضي بأن الموجات الصوتية، والتي تتضمن انضغاطاً وتمدداً متتاليين، يجب أن تشهد فقداناً في الطاقة (تبديد)، مما يعني أن λ يجب أن يكون معامل انتشار موجباً.
من خلال مساواة الضغوط الميكانيكية والثيرموديناميكية لإلغاء λ، تفشل معادلة نافير-ستوكس القياسية في تفسير التوهن الصوتي الناتج عن اللزوجة الحجمية، حيث تعزو الخسائر اللزوجية حصرياً إلى لزوجة القص (μ).
تخلق هذه الفرضية انفصالاً بين المعادلات الحاكمة المستخدمة في ميكانيكا الموائع العامة وتلك المطلوبة للدقة في انتشار الإشارات الصوتية الهوائية، لا سيما فيما يتعلق بالتشتت والتوهن.
المنهجية يمضي الاستقصاء من خلال إعادة اشتقاق نظري للمعادلات الحاكمة واستراتيجية تحقق عددية:
الاشتقاق النظري: بدءاً من حفظ الكتلة والزخم للمائع النيوتوني الانضغاطي دون استدعاء فرضية ستوكس، يشتق المؤلف معادلة ضغط الاضطراب. يستخدم هذا الاشتقاق لزوجة كينماتيكية معممة جديدة (νl)، مُعرفة بـ νl=(λ+2μ)/ρ0.
تصنيف المعادلات: يتم تحليل معادلة اضطراب الضغط الناتجة باستخدام التحليل الطيفي العالمي (GSA) لتحديد طبيعتها الرياضية (زائدية أم مكافئة) وخصائصها التشتتية.
التحقق العددي (DNS): تشير الدراسة إلى نتائج المحاكاة العددية المباشرة (DNS) لعدم استقرار رايلي-تايلور (RTI) في نطاق مكعب مغلق. في هذه المحاكاة، يستقر مائع ثقيل فوق مائع أخف في حالة سكون، وعند الإزالة الاندفاعية للفاصل، يولد النظام نبضات صوتية.
استراتيجية التقدير: تقترح الورقة منهجية لتقدير νl (وبالتالي λ) من خلال تحليل نبضات الضغط أحادية البعد، وتوازن بين طيف متعدد الترددات معقد لبداية عدم استقرار RTI وبين إعداد نظري مبسط يتضمن حزمة موجية غاوسية اندفاعية زمنية، وهو ما يكون أكثر ملاءمة للتقدير التحليلي والعددي لمعاملات اللزوجة.
المساهمات الرئيسية
صيغة موحدة: تقدم الورقة إطار عمل معدلاً لمعادلات نافير-ستوكس يدمج كلاً من لزوجة القص (ν) واللزوجة الكينماتيكية المعممة (νl). هذه الصيغة قابلة للتطبيق على كل من ميكانيكا الموائع والصوتيات الهوائية دون الحاجة لفرضية ستوكس.
معادلة الزخم المعدلة: تم اشتقاق معادلة حفظ زخم جديدة باستخدام التدوين المؤشر (Equation 3.2) تحتفظ صراحة بمعامل اللزوجة الثاني (λ) عبر الحد νl. وتعمل هذه المعادلة بشكل صحيح على نمذجة انتشار الزخم في التدفقات الانضغاطية كدالة لكل من تأثيرات القص والكتلة.
معادلة اضطراب الضغط: يشتق المؤلف معادلة حاكمة لاضطراب الضغط (Equation 2.5) تتضمن حداً للانتشار اللزوجي يعتمد على νl. وتوضح هذه المعادلة أن انتشار الإشارات الصوتية هو بطبيعته تشتتي وتبديدي، حيث تصبح سرعة الصوت معتمدة على الطول والمقياس الزمني، بدلاً من كونها عملية ثابتة الكظومة (adiabatic).
التبرير الفيزيائي: يقدم العمل حجة فيزيائية ضد فرضية ستوكس، مشيراً إلى أن λ السالب يعني عدم استقرار فيزيائي، بينما يتطلب الانتشار الصوتي تبديداً موجباً.
النتائج
التشتت والتبديد: تُظهر معادلة الضغط المشتقة (Equation 2.5) أن النظام تشتتي. ويشير علاقة التشتت (Equation 4.2) إلى أنه بالنسبة للأعداد الموجية العالية (k>kc=2c/νl)، تنتقل المعادلة الحاكمة من الطبيعة الزائدية (الموجية) إلى الطبيعة المكافئة البحتة، وهو سلوك يعتمد على νl.
ديناميكيات بداية RTI: يكشف تحليل بيانات المحاكة العددية المباشرة (DNS) عالية الدقة لـ RTI (بالإشارة إلى Sengupta et al., 2022) أن بداية عدم الاستقرار من حالة السكون مدفوعة بنبضات ضغط أحادية البعد تنتشر عمودياً على الفاصل، وتغطي هذه النبضات نطاقاً ترددياً من الترددات تحت الصوتية إلى فوق الصوتية.
عدم كفاية البيانات الحالية: بينما تم وضع الإطار النظري، تشير الورقة إلى أن البيانات التجريبية لمعامل اللزوجة الثاني (λ) نادرة، مع الاعتماد الأكبر على تحليلات الانحدار لبيانات التوهن الصوتي من Ash et al. (1991) و Zuckerwar & Ash (2009).
التحقق من الانتشار أحادي البعد: تؤكد الدراسة أنه رغم كون RTI ظاهرة ثلاثية الأبعاد، إلا أن الإثارة الصوتية الأولية تسلك سلوك نبضة ضغط أحادية البعد، مما يجعلها مرشحاً مناسباً لتقدير νl إذا كان الطيف أقل تعقيداً.
الأهمية والادعاءات تدعي الورقة تقديم "نهج موحد" لميكانيكا الموائع والصوتيات الهوائية من خلال إزالة القيد الاصطناعي لفرضية ستوكس. تكمن الأهمية الرئيسية في:
تصحيح التناقضات الفيزيائية: إزالة الطبيعة المضادة للانتشار غير الفيزيائية لمعامل اللزوجة الثاني في صيغ نافير-ستوكس الحالية.
تمكين الصوتيات الهوائية الدقيقة: توفير إطار عمل حيث يتم نمذجة توهن وتشتت الصوت بشكل طبيعي عبر νl، بدلاً من تجاهلها أو عزوها بشكل غير صحيح.
مسار تقدير جديد: اقتراح إمكانية تقدير νl تجريبياً أو عددياً عن طريق تتبع انتشار نبضات الضغط المستوية، مما يسمح بالتحديد الكامل لمعامل اللزوجة الثاني دون الاعتماد على فرضية تلاشي اللزوجة الحجمية.
تخلص الورقة إلى أنه بينما الإطار النظري مكتمل، فإن المهمة المتبقية الحاسمة هي التقدير التجريبي أو العددي الدقيق لـ νl للموائع الشائعة مثل الهواء والماء تحت ظروف محيطة متغيرة لتحقيق الإمكانات الكاملة لهذه الصيغة الموحدة.