وظيفة التدفق

بالنسبة لحقل متجه سرعة التدفق غير القابل للانضغاط (باللون الأحمر، في الأعلى)، يمكن حساب خطوط التدفق الخاصة به (المتقطعة) على أنها خطوط كفاف دالة التدفق (في الأسفل).

في ديناميكا الموائع ، يتم تعريف نوعين من دالة التدفق (أو دالة التدفق ):

إن خصائص دوال التدفق تجعلها مفيدة لتحليل التدفقات وتوضيحها بيانياً.

يصف الجزء المتبقي من هذه المقالة دالة التدفق ثنائية الأبعاد.

دالة التدفق ثنائية الأبعاد

الافتراضات

تعتمد دالة التدفق ثنائية الأبعاد على الافتراضات التالية:

  • يمكن وصف مجال التدفق بأنه تدفق مستوٍ ثنائي الأبعاد، مع متجه سرعة
u=[u(x،y،ت)v(x،y،ت)0].{\displaystyle \quad \mathbf {u} ={\begin{bmatrix}u(x,y,t)\\v(x,y,t)\\0\end{bmatrix}}.}
u=0.{\displaystyle \quad \nabla \cdot \mathbf {u} =0.}
  • لا يحتوي المجال على ثقوب، أو يحتوي فقط على ثقوب ليس لها تدفق صافٍ إلى الداخل أو إلى الخارج.

على الرغم من أن دالة التدفق لا تتطلب من حيث المبدأ استخدام نظام إحداثيات معين، إلا أن الوصف المقدم هنا يستخدم نظام إحداثيات ديكارتية يمينية مع إحداثيات(x،y،z){\displaystyle (x,y,z)}.

الاشتقاق

سطح الاختبار

ضع في اعتبارك نقطتينأ{\displaystyle A}وP{\displaystyle P}فيxy{\displaystyle xy}مستوى، ومنحنى متصلأP{\displaystyle AP}، وكذلك فيxy{\displaystyle xy}المستوى الذي يربط بينهما حيث كل نقطة على المنحنىأP{\displaystyle AP}لديهz{\displaystyle z}تنسيقz=0{\displaystyle z=0}لنفترض أن الطول الكلي للمنحنىأP{\displaystyle AP}يكونل{\displaystyle L}.

لنفترض أنه تم إنشاء سطح على شكل شريط عن طريق تمديد المنحنىأP{\displaystyle AP}لأعلى إلى المستوى الأفقيz=ب{\displaystyle z=b}(ب>0){\displaystyle (b>0)}، أينب{\displaystyle b}يمثل سمك التدفق. عندئذٍ يكون للسطح طولل{\displaystyle L}، عرضب{\displaystyle b}والمنطقةبل{\displaystyle b\,L}لنسمي هذا سطح الاختبار .

التدفق عبر سطح الاختبار

التدفق الحجمي عبر سطح الاختبار الذي يربط النقاطأ{\displaystyle A}وP.{\displaystyle P.}

إجمالي التدفق الحجمي عبر سطح الاختبار هو

سؤال(x،y،ت)=0ب0لuن^دsدz{\displaystyle Q(x,y,t)=\int _{0}^{b}\int _{0}^{L}\mathbf {u} \cdot {\hat {\mathbf {n} }}\,\mathrm {d} s\,\mathrm {d} z}

أينs{\displaystyle s}هو مُعامل طول القوس المُعرَّف على المنحنىأP{\displaystyle AP}، معs=0{\displaystyle s=0}عند النقطةأ{\displaystyle A}وs=ل{\displaystyle s=L}عند النقطةP{\displaystyle P}. هنان^{\displaystyle {\hat {\mathbf {n} }}}هو متجه الوحدة العمودي على سطح الاختبار، أي

ن^دs=-Rدر=[دy-دx0]{\displaystyle {\hat {\mathbf {n} }}\,\mathrm {d} s=-R\,\mathrm {d} \mathbf {r} ={\begin{bmatrix}\mathrm {d} y\\-\mathrm {d} x\\0\end{bmatrix}}}

أينR{\displaystyle R}هو3×3{\displaystyle 3\times 3}مصفوفة الدوران المقابلة لـ أ90{\displaystyle 90^{\circ }}دوران عكس اتجاه عقارب الساعة حول النقطة الموجبةz{\displaystyle z}محور:

R=Rz(90)=[0-10100001].{\displaystyle R=R_{z}(90^{\circ })={\begin{bmatrix}0&-1&0\\1&0&0\\0&0&1\end{bmatrix}}.}

الدالة التكاملية في التعبير عنسؤال{\displaystyle Q}مستقل عنz{\displaystyle z}لذا يمكن حساب التكامل الخارجي للحصول على

سؤال(x،y،ت)=بأP(uدy-vدx){\displaystyle Q(x,y,t)=b\,\int _{A}^{P}\left(u\,\mathrm {d} y-v\,\mathrm {d} x\right)}

التعريف الكلاسيكي

يُعرّف لامب وباتشيلور دالة التدفقψ{\displaystyle \psi }على النحو التالي. [ 3 ]

ψ(x،y،ت)=أP(uدy-vدx){\displaystyle \psi (x,y,t)=\int _{A}^{P}\left(u\,\mathrm {d} y-v\,\mathrm {d} x\right)}

باستخدام الصيغة المشتقة أعلاه للتدفق الحجمي الكلي،سؤال{\displaystyle Q}ويمكن كتابة ذلك على النحو التالي

ψ(x،y،ت)=سؤال(x،y،ت)ب{\displaystyle \psi (x,y,t)={\frac {Q(x,y,t)}{b}}}.

بعبارة أخرى، دالة التدفقψ{\displaystyle \psi }هو التدفق الحجمي عبر سطح الاختبار لكل وحدة سمك، حيث يتم قياس السمك بشكل عمودي على مستوى التدفق.

النقطةأ{\displaystyle A}هي نقطة مرجعية تحدد المكان الذي تكون فيه دالة التدفق صفراً تماماً. يتم اختيار موضعها بشكل عشوائي إلى حد ما، وبمجرد اختيارها، تظل ثابتة عادةً.

تحول متناهي الصغردP=(دx،دy){\displaystyle \mathrm {d} P=(\mathrm {d} x,\mathrm {d} y)}في موضع النقطةP{\displaystyle P}ينتج عن ذلك التغيير التالي في دالة التدفق:

دψ=uدy-vدx{\displaystyle \mathrm {d} \psi =u\,\mathrm {d} y-v\,\mathrm {d} x}.

من التفاضل التام

دψ=ψxدx+ψyدy،{\displaystyle \mathrm {d} \psi ={\frac {\partial \psi }{\partial x}}\,\mathrm {d} x+{\frac {\partial \psi }{\partial y}}\,\mathrm {d} y,}

لذا فإن مكونات سرعة التدفق بالنسبة لدالة التدفقψ{\displaystyle \psi }يجب أن يكون

u=ψy،v=-ψx.{\displaystyle u={\frac {\partial \psi }{\partial y}},\qquad v=-{\frac {\partial \psi }{\partial x}}.}

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

تأثير تغير موضع نقطة المرجع

لنفترض حدوث تغيير في موضع نقطة المرجع، على سبيل المثال منأ{\displaystyle A}لأ{\displaystyle A'}. يتركψ{\displaystyle \psi '}تشير إلى دالة التدفق بالنسبة إلى نقطة مرجعية مُزاحةأ{\displaystyle A'}:

ψ(x،y،ت)=أP(uدy-vدx).{\displaystyle \psi '(x,y,t)=\int _{A'}^{P}\left(u\,\mathrm {d} y-v\,\mathrm {d} x\right).}

ثم يتم إزاحة دالة التدفق بمقدار

Δψ(ت)=ψ(x،y،ت)-ψ(x،y،ت)=أأ(uدy-vدx)،{\displaystyle {\begin{aligned}\Delta \psi (t)&=\psi '(x,y,t)-\psi (x,y,t)\\&=\int _{A'}^{A}\left(u\,\mathrm {d} y-v\,\mathrm {d} x\right),\end{aligned}}}

وهذا يعني ما يلي:

  • يؤدي تغيير موضع نقطة المرجع فعليًا إلى إضافة ثابت (للتدفق المستقر) أو دالة تعتمد على الزمن فقط (للتدفق غير المستقر) إلى دالة التدفقψ{\displaystyle \psi }في كل نقطةP{\displaystyle P}.
  • التحول في دالة التدفق،Δψ{\displaystyle \Delta \psi }، يساوي التدفق الحجمي الكلي، لكل وحدة سمك، عبر السطح المتصل الذي يمتد من النقطةأ{\displaystyle A'}للإشارةأ{\displaystyle A}. بالتاليΔψ=0{\displaystyle \Delta \psi =0}إذا وفقط إذاأ{\displaystyle A}وأ{\displaystyle A'}تقع على نفس الخط الانسيابي.

من حيث دوران المتجه

السرعةu{\displaystyle \mathbf {u} }يمكن التعبير عنها بدلالة دالة التدفقψ{\displaystyle \psi }مثل

u=-Rψ{\displaystyle \mathbf {u} =-R\,\nabla \psi }

أينR{\displaystyle R}هو3×3{\displaystyle 3\times 3}مصفوفة الدوران المقابلة لـ أ90{\displaystyle 90^{\circ }}دوران عكس اتجاه عقارب الساعة حول النقطة الموجبةz{\displaystyle z}المحور. حل المعادلة أعلاه لـψ{\displaystyle \nabla \psi }ينتج الشكل المكافئ

ψ=Ru.{\displaystyle \nabla \psi =R\,\mathbf {u} .}

يتضح من هذه الأشكال على الفور أن المتجهاتu{\displaystyle \mathbf {u} }وψ{\displaystyle \nabla \psi }نكون

  • عمودي: uψ=0{\displaystyle \mathbf {u} \cdot \nabla \psi =0}
  • بنفس الطول: |u|=|ψ|{\displaystyle |\mathbf {u} |=|\nabla \psi |}.

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

من حيث الجهد المتجهي وأسطح التدفق

بشكل عام، مجال خالٍ من التباعد مثلu{\displaystyle \mathbf {u} }، والمعروف أيضًا باسم حقل متجه لولبي ، يمكن دائمًا تمثيله على أنه التفاف لبعض الكمون المتجهيأ{\displaystyle {\boldsymbol {A}}}:

u=×أ.{\displaystyle \mathbf {u} =\nabla \times {\boldsymbol {A}}.}

دالة التدفقψ{\displaystyle \psi }يمكن فهم ذلك على أنه يوفر قوة جهد متجه موجه عموديًا على المستوى: [ 4 ]

أ(x،y،ت)=[00ψ(x،y،ت)]،{\displaystyle {\boldsymbol {A}}(x,y,t)={\begin{bmatrix}0\\0\\\psi (x,y,t)\end{bmatrix}},}

بعبارة أخرىأ=ψz^{\displaystyle {\boldsymbol {A}}=\psi {\hat {\mathbf {z} }}}، أينz^{\displaystyle {\hat {\mathbf {z} }}}متجه الوحدة يشير إلى الاتجاه الموجبz{\displaystyle z}اتجاه.

ويمكن كتابة ذلك أيضًا على شكل الضرب الاتجاهي للمتجهات

u=ψ×z^{\displaystyle \mathbf {u} =\nabla \psi \times {\hat {\mathbf {z} }}}

حيث استخدمنا متطابقة حساب المتجهات

×(ψz^)=ψ×z^+ψ×z^.{\displaystyle \nabla \times \left(\psi {\hat {\mathbf {z} }}\right)=\psi \nabla \times {\hat {\mathbf {z} }}+\nabla \psi \times {\hat {\mathbf {z} }}.}

مع ملاحظة أنz^=z{\displaystyle {\hat {\mathbf {z} }}=\nabla z}، وتحديدϕ=z{\displaystyle \phi =z}يمكن التعبير عن حقل السرعة على النحو التالي:

u=ψ×ϕ.{\displaystyle \mathbf {u} =\nabla \psi \times \nabla \phi .}

يوضح هذا الشكل أن الأسطح المستوية لـψ{\displaystyle \psi }والأسطح المستوية لـz{\displaystyle z}(أي المستويات الأفقية) تشكل نظامًا من أسطح التدفق المتعامدة .

تعريف بديل (علامة معاكسة)

وهناك تعريف بديل، يُستخدم أحيانًا في علم الأرصاد الجوية وعلم المحيطات ، وهو

ψ=-ψ.{\displaystyle \psi '=-\psi .}

العلاقة بالدوامية

في التدفق المستوي ثنائي الأبعاد، يُعرف متجه الدوامة على النحو التالي:ω=×u{\displaystyle {\boldsymbol {\omega }}=\nabla \times \mathbf {u} }، يختزل إلىωz^{\displaystyle \omega \,{\hat {\mathbf {z} }}}، أين

ω=-2ψ{\displaystyle \omega =-\nabla ^{2}\psi }

أو

ω=+2ψ{\displaystyle \omega =+\nabla ^{2}\psi '}

هذه هي أشكال معادلة بواسون .

العلاقة بالخطوط الانسيابية

لنفترض جريانًا مستويًا ثنائي الأبعاد بنقطتين متقاربتين للغايةP=(x،y،z){\displaystyle P=(x,y,z)}وP=(x+دx،y+دy،z){\displaystyle P'=(x+dx,y+dy,z)}تقعان في نفس المستوى الأفقي. من حساب التفاضل والتكامل، يكون الفرق المتناهي الصغر المقابل بين قيم دالة التدفق عند النقطتين هو

دψ(x،y،ت)=ψ(x+دx،y+دy،ت)-ψ(x،y،ت)=ψxدx+ψyدy=ψدر{\displaystyle {\begin{aligned}\mathrm {d} \psi (x,y,t)&=\psi (x+\mathrm {d} x,y+\mathrm {d} y,t)-\psi (x,y,t)\\&={\partial \psi \over \partial x}\mathrm {d} x+{\partial \psi \over \partial y}\mathrm {d} y\\&=\nabla \psi \cdot \mathrm {d} \mathbf {r} \end{aligned}}}

يفترضψ{\displaystyle \psi }يأخذ نفس القيمة، على سبيل المثالج{\displaystyle C}، عند النقطتينP{\displaystyle P}وP{\displaystyle P'}ثم هذا يعطي

0=ψدر،{\displaystyle 0=\nabla \psi \cdot \mathrm {d} \mathbf {r} ,}

مما يعني أن المتجهψ{\displaystyle \nabla \psi }عمودي على السطحψ=ج{\displaystyle \psi =C}. لأنuψ=0{\displaystyle \mathbf {u} \cdot \nabla \psi =0}في كل مكان (انظر على سبيل المثال، من حيث دوران المتجهات )، يتوافق كل خط انسيابي مع تقاطع سطح تيار معين ومستوى أفقي معين. وبالتالي، في ثلاثة أبعاد، يتطلب التحديد الدقيق لأي خط انسيابي معين تحديد القيم المقابلة لكل من دالة التيار والارتفاع (z{\displaystyle z}تنسيق).

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

حالة الوجود

من السهل إثبات ذلك بالنسبة للتدفق المستوي ثنائي الأبعادu{\displaystyle \mathbf {u} }يحقق معادلة الالتفاف والتباعد

(u)z^=-×(Ru){\displaystyle (\nabla \cdot \mathbf {u} )\,{\hat {\mathbf {z} }}=-\nabla \times (R\,\mathbf {u} )}

أينR{\displaystyle R}هو3×3{\displaystyle 3\times 3}مصفوفة الدوران المقابلة لـ أ90{\displaystyle 90^{\circ }}دوران عكس اتجاه عقارب الساعة حول النقطة الموجبةz{\displaystyle z}المحور. هذه المعادلة صحيحة بغض النظر عما إذا كان التدفق غير قابل للانضغاط أم لا.

إذا كان التدفق غير قابل للانضغاط (أي،u=0{\displaystyle \nabla \cdot \mathbf {u} =0}ثم تعطي معادلة الالتفاف والتباعد

0=×(Ru){\displaystyle \mathbf {0} =\nabla \times (R\,\mathbf {u} )}.

ثم بحسب نظرية ستوكس، فإن التكامل الخطي لـRu{\displaystyle R\,\mathbf {u} }يختفي فوق كل حلقة مغلقة

Σ(Ru)دΓ=0.{\displaystyle \oint _{\partial \Sigma }(R\,\mathbf {u} )\cdot \mathrm {d} \mathbf {\Gamma } =0.}

وبالتالي، فإن التكامل الخطي لـRu{\displaystyle R\,\mathbf {u} }مستقلة عن المسار. وأخيرًا، بحسب عكس نظرية التدرج ، فإن الدالة العدديةψ(x،y،ت){\displaystyle \psi (x,y,t)}يوجد بحيث

Ru=ψ{\displaystyle R\,\mathbf {u} =\nabla \psi }.

هناψ{\displaystyle \psi }يمثل دالة التدفق.

وعلى العكس من ذلك، إذا كانت دالة التدفق موجودة، فإنRu=ψ{\displaystyle R\,\mathbf {u} =\nabla \psi }. بتعويض هذه النتيجة في معادلة الالتفاف والتباعد، نحصل علىu=0{\displaystyle \nabla \cdot \mathbf {u} =0}(أي أن التدفق غير قابل للانضغاط).

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

التدفق المحتمل

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

ملخص الخصائص

يمكن تلخيص الخصائص الأساسية لوظائف التدفق ثنائية الأبعاد على النحو التالي:

  1. يتم تحديد مركبتي السرعة x و y عند نقطة معينة من خلال المشتقات الجزئية لدالة التدفق عند تلك النقطة.
  2. تكون قيمة دالة التدفق ثابتة على طول كل خط تدفق (تمثل خطوط التدفق مسارات الجسيمات في التدفق المستقر). أي أن كل خط تدفق في بعدين هو منحنى مستوى لدالة التدفق.
  3. إن الفرق بين قيم دالة التدفق عند أي نقطتين يعطي التدفق الحجمي عبر السطح الرأسي الذي يربط النقطتين.

دالة التدفق ثنائية الأبعاد للتدفقات ذات الكثافة الثابتة مع الزمن

إذا كانت كثافة السائل ثابتة مع الزمن عند جميع النقاط داخل التدفق، أي

ρت=0{\displaystyle {\frac {\partial \rho }{\partial t}}=0}،

عندئذٍ تصبح معادلة الاستمرارية (انظر على سبيل المثال معادلة الاستمرارية#ديناميكا الموائع ) للجريان المستوي ثنائي الأبعاد كما يلي

(ρu)=0.{\displaystyle \nabla \cdot (\rho \,\mathbf {u} )=0.}

في هذه الحالة دالة التدفقψ{\displaystyle \psi }يتم تعريفها بحيث

ρu=ψy،ρv=-ψx{\displaystyle \rho \,u={\frac {\partial \psi }{\partial y}},\quad \rho \,v=-{\frac {\partial \psi }{\partial x}}}

ويمثل التدفق الكتلي (بدلاً من التدفق الحجمي) لكل وحدة سمك عبر سطح الاختبار.

انظر أيضاً

مراجع

الاقتباسات

  1. ^ لاغرانج، ج.-ل. (1868)، “Mémoire sur la théorie du mouvement des Fluides (in: Nouveaux Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Berlin، année 1781)”، Oevres de Lagrange ، المجلد.  المجلد الرابع، ص 695 – 748 
  2. ستوكس، جي جي (1842)، "حول الحركة المستقرة للسوائل غير القابلة للانضغاط"، معاملات الجمعية الفلسفية في كامبريدج ، 7 : 439-453 ، رمز Bibcode : 1848TCaPS...7..439Sأُعيد طبعه في: ستوكس، جي جي (1880)، أوراق رياضية وفيزيائية، المجلد الأول ، مطبعة جامعة كامبريدج ، الصفحات 1-16 
  3. لامب (1932 ، ص 62-63) وباتشيلور (1967 ، ص 75-79)  
  4. كاتوبوديس، نيكولاوس د. (2019). "تدفق الموائع اللزجة". تدفق السطح الحر . ص 324-426 . doi : 10.1016/B978-0-12-815489-2.00005-8 . ISBN  978-0-12-815489-2.

مصادر

  • باتشيلور، جي كي (1967)، مقدمة في ديناميكا الموائع ، مطبعة جامعة كامبريدج، رقم ISBN 0-521-09817-3
  • لامب، هـ. (1932)، الديناميكا المائية (  الطبعة السادسة)، مطبعة جامعة كامبريدج، أعيد نشرها بواسطة منشورات دوفر، رقم ISBN 0-486-60256-7{{citation}}عدم توافق رقم ISBN / التاريخ ( مساعدة )
  • ماسي، بي إس؛ وارد-سميث، جيه. (1998)، ميكانيكا الموائع (  الطبعة السابعة)، المملكة المتحدة: نيلسون ثورنز
  • وايت، إف إم (2003)، ميكانيكا الموائع (  الطبعة الخامسة)، نيويورك: ماكجرو هيل
  • جاملين، تي دبليو (2001)، التحليل المركب ، نيويورك: سبرينغر، رقم ISBN 0-387-95093-1
  • "وظيفة التدفق" ، معجم الأرصاد الجوية التابع للجمعية الأمريكية للأرصاد الجوية ، تم الاطلاع عليه بتاريخ 30 يناير 2014