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

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

في ديناميكا الموائع ، يشير التدفق الكامن أو التدفق غير الدوراني إلى التدفق المثالي عديم الاحتكاك للمائع. ويتم تصوير نوعين من التدفقات على النحو التالي:

  1. تدفق سائل غير لزج
  2. تدفق سائل ذي لزوجة منخفضة ، في مناطق لا تحتوي على طبقة حدية . انظر فرضية براندتل .

يصف التدفق الكامن حقل السرعة على أنه تدرج دالة قياسية: جهد السرعة . ونتيجة لذلك، يتميز التدفق الكامن بحقل سرعة غير دوراني ، وهو تقريب صالح للعديد من التطبيقات. ويعود عدم دوران التدفق الكامن إلى أن التفاف تدرج الدالة القياسية يساوي صفرًا دائمًا.

في حالة التدفق غير القابل للانضغاط، يحقق جهد السرعة معادلة لابلاس ، ويمكن تطبيق نظرية الجهد . مع ذلك، استُخدمت تدفقات الجهد أيضًا لوصف التدفقات القابلة للانضغاط وتدفقات هيلي-شو . يُستخدم نهج تدفق الجهد في نمذجة كلٍ من التدفقات الثابتة وغير الثابتة. تشمل تطبيقات تدفق الجهد: مجال التدفق الخارجي للأجنحة الهوائية ، وأمواج الماء ، والتدفق الكهروأسموزي ، وتدفق المياه الجوفية .

يمكن وصف المنطقة التي تحتوي على قوى احتكاك (تُعرف بقوى القص وقوى اللزوجة) بأنها منطقة تحتوي على دوامات . بالنسبة للتدفقات (أو أجزاء منها) ذات التأثيرات الدوامية القوية، لا يُمكن تطبيق تقريب التدفق الكامن. في مناطق التدفق التي تُعرف فيها أهمية الدوامات، مثل مناطق التخلف والطبقات الحدية ، لا تستطيع نظرية التدفق الكامن تقديم تنبؤات معقولة للتدفق. [ 1 ] ومع ذلك، غالبًا ما توجد مناطق واسعة من التدفق يكون فيها افتراض انعدام الدوامات صحيحًا، مما يسمح باستخدام التدفق الكامن في تطبيقات متنوعة؛ تشمل هذه التطبيقات التدفق حول الطائرات ، وتدفق المياه الجوفية ، وعلم الصوتيات ، وأمواج الماء ، والتدفق الكهروأسموزي . [ 2 ]

الوصف والخصائص

يتم إنشاء التدفق المحتمل عن طريق إضافة تدفقات أولية بسيطة وملاحظة النتيجة.
خطوط انسيابية لتدفق الكامن غير القابل للانضغاط حول أسطوانة دائرية في تدفق منتظم.

في التدفق الكامن أو غير الدوراني، يكون حقل متجه الدوامة صفراً، أي

ω×v=0،{\displaystyle {\boldsymbol {\omega }}\equiv \nabla \times \mathbf {v} =0,}

أينv(x،ت){\displaystyle \mathbf {v} (\mathbf {x} ,t)}هو حقل السرعة وω(x،ت){\displaystyle {\boldsymbol {\omega }}(\mathbf {x} ,t)}يمثل حقل الدوامة . ومثل أي حقل متجه ذي دوران صفري، يمكن التعبير عن حقل السرعة على أنه تدرج كمية قياسية معينة، على سبيل المثالφ(x،ت){\displaystyle \varphi (\mathbf {x} ,t)}وهذا ما يُسمى بجهد السرعة ، لأن التفاف التدرج يساوي صفرًا دائمًا. لذلك لدينا [ 3 ]

v=φ.{\displaystyle \mathbf {v} =\nabla \varphi .}

إن جهد السرعة ليس محددًا بشكل فريد، إذ يمكن إضافة دالة زمنية اختيارية إليه، على سبيل المثالو(ت){\displaystyle f(t)}، دون التأثير على الكمية الفيزيائية ذات الصلة وهيv{\displaystyle \mathbf {v} }عادةً ما يتم التخلص من عدم التفرد عن طريق اختيار الشروط الأولية أو الحدودية المناسبة التي تحققهاφ{\displaystyle \varphi }وبالتالي، قد يختلف الإجراء من مشكلة إلى أخرى.

في التدفق الكامن، الدورانΓ{\displaystyle \Gamma }حول أي محيط متصل ببساطةج{\displaystyle C}يساوي صفرًا. ويمكن إثبات ذلك باستخدام نظرية ستوكس .

Γجvدل=ωدو=0{\displaystyle \Gamma \equiv \oint _{C}\mathbf {v} \cdot d\mathbf {l} =\int {\boldsymbol {\omega }}\cdot d\mathbf {f} =0}

أيندل{\displaystyle d\mathbf {l} }هو عنصر الخط على المحيط ودو{\displaystyle d\mathbf {f} }هي عنصر المساحة لأي سطح محدد بالمحيط. في الفضاء متعدد الاتصالات (على سبيل المثال، حول محيط يحيط بجسم صلب في بعدين أو حول محيط يحيط بحلقة في ثلاثة أبعاد) أو في وجود دوامات مركزة (على سبيل المثال، في ما يسمى بالدوامات غير الدورانية أو الدوامات النقطية، أو في حلقات الدخان)، فإن الدورانΓ{\displaystyle \Gamma }ليس بالضرورة أن تكون القيمة صفرًا. في الحالة الأولى، لا يمكن تطبيق نظرية ستوكس، وفي الحالة الثانية،ω{\displaystyle {\boldsymbol {\أوميغا }}}لا يساوي الصفر داخل المنطقة المحصورة بالمحيط. حول محيط يحيط بأسطوانة صلبة لا نهائية الطول يلتف حولها المحيطشمال{\displaystyle N}في بعض الأحيان، لدينا

Γ=شمالκ{\displaystyle \Gamma =N\kappa }

أينκ{\displaystyle \kappa }هو ثابت دوري. ينتمي هذا المثال إلى فضاء متصل ثنائيًا. فين{\displaystyle n}- مساحة متصلة بشكل متماثل، هناكن-1{\displaystyle n-1}هذه الثوابت الدورية، وهي:κ1،κ2،...،κن-1.{\displaystyle \kappa _{1},\kappa _{2},\dots ,\kappa _{n-1}.}

التدفق غير القابل للانضغاط

في حالة التدفق غير القابل للانضغاط - على سبيل المثال السائل ، أو الغاز عند أرقام ماخ المنخفضة ؛ ولكن ليس بالنسبة للموجات الصوتية - فإن السرعة v لها تباعد صفري : [ 3 ]

v=0،{\displaystyle \nabla \cdot \mathbf {v} =0\,,}

استبدال هناv=φ{\displaystyle \mathbf {v} =\nabla \varphi }يُظهر ذلك أنφ{\displaystyle \varphi }يحقق معادلة لابلاس [ 3 ]

2φ=0،{\displaystyle \nabla ^{2}\varphi =0\,,}

حيث ∇² = ⋅ ∇ هو مؤثر لابلاس (ويُكتب أحيانًا Δ ). بما أن حلول معادلة لابلاس هي دوال توافقية ، فإن كل دالة توافقية تمثل حلاً لتدفق كامن. وكما هو واضح، في حالة التدفق غير القابل للانضغاط، يُحدد حقل السرعة كليًا من حركيته : بافتراض عدم الدوران وانعدام تباعد التدفق. أما الديناميكا المرتبطة بمعادلات الزخم، فلا يلزم تطبيقها إلا لاحقًا، إذا كان الهدف هو حساب حقل الضغط: على سبيل المثال، للتدفق حول الأجنحة باستخدام مبدأ برنولي .

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

μ2v=μ(v)-μ×ω=0{\displaystyle \mu \nabla ^{2}\mathbf {v} =\mu \nabla (\nabla \cdot \mathbf {v} )-\mu \nabla \times {\boldsymbol {\omega }}=0}

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

في بعدين، بمساعدة الدالة التوافقيةφ{\displaystyle \varphi }ودالتها التوافقية المترافقةψ{\displaystyle \psi }(دالة التدفق)، يختزل التدفق الكامن غير القابل للانضغاط إلى نظام بسيط للغاية يتم تحليله باستخدام التحليل المركب (انظر أدناه).

التدفق القابل للانضغاط

تدفق ثابت

يمكن أيضًا استخدام نظرية التدفق الكامن لنمذجة التدفق الانضغاطي غير الدوراني. اشتقاق المعادلة الحاكمة لـφ{\displaystyle \varphi }إن استنتاج معادلات أويلر أمرٌ بسيطٌ للغاية. تُعطى معادلات الاستمرارية ومعادلات الزخم (التدفق الكامن) للتدفقات المستقرة كما يلي:

ρv+vρ=0،(v)v=-1ρص=-ج2ρρ// }

حيث تُستنتج المعادلة الأخيرة من حقيقة أن الإنتروبيا ثابتة بالنسبة لجزيء سائل وأن مربع سرعة الصوت هوج2=(ص/ρ)s{\displaystyle c^{2}=(\partial p/\partial \rho )_{s}}. القضاء علىρ{\displaystyle \nabla \rho }ينتج عن المعادلتين الحاكمتين ما يلي:

ج2v-v(v)v=0.{\displaystyle c^{2}\nabla \cdot \mathbf {v} -\mathbf {v} \cdot (\mathbf {v} \cdot \nabla )\mathbf {v} =0.}

تظهر النسخة غير القابلة للانضغاط في النهايةج{\displaystyle c\to \infty }. استبدال هناv=φ{\displaystyle \mathbf {v} =\nabla \varphi }النتائج في [ 4 ] [ 5 ]

(ج2-φx2)φxx+(ج2-φy2)φyy+(ج2-φz2)φzz-2(φxφyφxy+φyφzφyz+φzφxϕzx)=0// _ {xy}+\varphi _ {y}\varphi _ {z}\varphi _ {yz}+\varphi _ {z}\varphi _ {x}\phi _ {zx}) = 0}

أينج=ج(v){\displaystyle c=c(v)}يتم التعبير عنها كدالة لمقدار السرعةv2=(ϕ)2{\displaystyle v^{2}=(\nabla \phi )^{2}}بالنسبة للغاز متعدد الخواص ،ج2=(γ-1)(ح0-v2/2){\displaystyle c^{2}=(\gamma -1)(h_{0}-v^{2}/2)}، أينγ{\displaystyle \gamma }هي نسبة الحرارة النوعية وح0{\displaystyle h_{0}}هي إنثالبي الركود . في بعدين، تُبسط المعادلة إلى

(ج2-φx2)φxx+(ج2-φy2)φyy-2φxφyφxy=0.{\displaystyle (c^{2}-\varphi _{x}^{2})\varphi _{xx}+(c^{2}-\varphi _{y}^{2})\varphi _{yy}-2\varphi _{x}\varphi _{y}\varphi _{xy}=0.}

الصلاحية: في وضعها الحالي، تُعدّ المعادلة صالحة لأي تدفقات كامنة غير لزجة، بغض النظر عما إذا كان التدفق دون صوتي أو فوق صوتي (مثل تدفق براندتل-ماير ). مع ذلك، في التدفقات فوق الصوتية والتدفقات الانتقالية، قد تحدث موجات صدمية تُدخل إنتروبيا ودوامية في التدفق، مما يجعله دورانيًا. ومع ذلك، توجد حالتان يسود فيهما التدفق الكامن حتى في وجود الموجات الصدمية، وهما مُفسَّرتان من معادلة الزخم (التي ليست بالضرورة كامنة) المكتوبة بالشكل التالي.

(ح+v2/2)-v×ω=تيs{\displaystyle \nabla (h+v^{2}/2)-\mathbf {v} \times {\boldsymbol {\omega }}=T\nabla s}

أينح{\displaystyle h}هي المحتوى الحراري النوعي ،ω{\displaystyle {\boldsymbol {\omega }}}هو حقل الدوامة ،تي{\displaystyle T}درجة الحرارة وs{\displaystyle s}هي الإنتروبيا النوعية. وبما أن لدينا تدفقًا كامنًا أمام الموجة الصدمية الأمامية، فإن معادلة برنولي تُظهر أنح+v2/2{\displaystyle h+v^{2}/2}ثابت، وهو ثابت أيضًا عبر موجة الصدمة ( شروط رانكين-هوغونيو ) وبالتالي يمكننا كتابة [ 4 ]

v×ω=-تيs{\displaystyle \mathbf {v} \times {\boldsymbol {\omega }}=-T\nabla s}

1) عندما تكون شدة الموجة الصدمية ثابتة، يكون عدم استمرارية الإنتروبيا عبر الموجة الصدمية ثابتًا أيضًا، أيs=0{\displaystyle \nabla s=0}وبالتالي، يكون إنتاج الدوامات معدومًا. تتميز موجات الصدمة عند الحافة الأمامية المدببة للوتد ثنائي الأبعاد أو المخروط ثلاثي الأبعاد ( تدفق تايلور-ماكول ) بشدة ثابتة. 2) بالنسبة لموجات الصدمة الضعيفة، تكون قفزة الإنتروبيا عبر موجة الصدمة كمية من الدرجة الثالثة بالنسبة لقوة موجة الصدمة، وبالتاليs{\displaystyle \nabla s}يمكن إهمالها. تقع موجات الصدمة في الأجسام النحيلة بشكل شبه موازٍ للجسم وهي ضعيفة.

التدفقات شبه المتوازية: عندما يكون التدفق أحادي الاتجاه في الغالب مع انحرافات صغيرة كما هو الحال في التدفق حول الأجسام النحيلة، يمكن تبسيط المعادلة الكاملة بشكل أكبر. لنفترضيوهـx{\displaystyle U\mathbf {e} _{x}}لنفترض أن السرعة هي التيار السائد، ولنأخذ في الاعتبار الانحرافات الطفيفة عن مجال السرعة هذا. ويمكن كتابة جهد السرعة المقابل على النحو التالي:φ=xيو+ϕ{\displaystyle \varphi =xU+\phi }أينϕ{\displaystyle \phi }يُحدد هذا الانحراف الطفيف عن التدفق المنتظم، ويُحقق الصيغة الخطية للمعادلة الكاملة. ويُعطى هذا بواسطة

(1-م2)2ϕx2+2ϕy2+2ϕz2=0{\displaystyle (1-M^{2}){\frac {\partial ^{2}\phi }{\partial x^{2}}}+{\frac {\partial ^{2}\phi }{\partial y^{2}}}+{\frac {\partial ^{2}\phi }{\partial z^{2}}}=0}

أينم=يو/ج{\displaystyle M=U/c_{\infty }}يمثل عدد ماخ الثابت المقابل للتدفق المنتظم. وتكون هذه المعادلة صحيحة بشرطم{\displaystyle M}لا يقترب من الوحدة. عندما|م-1|{\displaystyle |M-1|}إذا كانت صغيرة (تدفق فوق صوتي)، فلدينا المعادلة غير الخطية التالية [ 4 ]

2α*ϕx2ϕx2=2ϕy2+2ϕz2{\displaystyle 2\alpha _{*}{\frac {\partial \phi }{\partial x}}{\frac {\partial ^{2}\phi }{\partial x^{2}}}={\frac {\partial ^{2}\phi }{\partial y^{2}}}+{\frac {\partial ^{2}\phi }{\partial z^{2}}}}

أينα*{\displaystyle \alpha _{*}}القيمة الحرجة لمشتقة لاندوα=(ج4/2υ3)(2υ/ص2)s{\displaystyle \alpha =(c^{4}/2\upsilon ^{3})(\partial ^{2}\upsilon /\partial p^{2})_{s}}[ 6 ] [ 7 ] وυ=1/ρ{\displaystyle \upsilon =1/\rho }يمثل الحجم النوعي . يتميز التدفق فوق الصوتي تمامًا بمعامل واحد.α*{\displaystyle \alpha _{*}}، والتي تأخذ القيمة بالنسبة للغاز متعدد الخواصα*=α=(γ+1)/2{\displaystyle \alpha _{*}=\alpha =(\gamma +1)/2}. تحت تحويل الهودوغراف ، تصبح المعادلة فوق الصوتية في بعدين معادلة أويلر-تريكومي .

تدفق غير مستقر

تُعطى معادلات الاستمرارية ومعادلات الزخم (التدفق الكامن) للتدفقات غير المستقرة كما يلي:

ρت+ρv+vρ=0،vت+(v)v=-1ρص=-ج2ρρ=-ح.{\displaystyle {\frac {\partial \rho }{\partial t}}+\rho \nabla \cdot \mathbf {v} +\mathbf {v} \cdot \nabla \rho =0,\quad {\frac {\partial \mathbf {v} }{\partial t}}+(\mathbf {v} \cdot \nabla )\mathbf {v} =-{\frac {1}{\rho }}\nabla p=-{\frac {c^{2}}{\rho }}\nabla \rho =-\nabla h.}

يُعطى التكامل الأول لمعادلة الزخم (التدفق الكامن) بالصيغة التالية:

φت+v22+ح=و(ت)،حت=-2φت2-12v2ت+دودت{\displaystyle {\frac {\partial \varphi }{\partial t}}+{\frac {v^{2}}{2}}+h=f(t),\quad \Rightarrow \quad {\frac {\partial h}{\partial t}}=-{\frac {\partial ^{2}\varphi }{\partial t^{2}}}-{\frac {1}{2}}{\frac {\partial v^{2}}{\partial t}}+{\frac {df}{dt}}}

أينو(ت){\displaystyle f(t)}هي دالة اختيارية. وبدون فقدان للعمومية، يمكننا أن نضعو(ت)=0{\displaystyle f(t)=0}منذφ{\displaystyle \varphi }لا يتم تعريفها بشكل فريد. بدمج هذه المعادلات، نحصل على

2φت2+v2ت=ج2v-v(v)v.{\displaystyle {\frac {\partial ^{2}\varphi }{\partial t^{2}}}+{\frac {\partial v^{2}}{\partial t}}=c^{2}\nabla \cdot \mathbf {v} -\mathbf {v} \cdot (\mathbf {v} \cdot \nabla )\mathbf {v} .}

استبدال هناv=φ{\displaystyle \mathbf {v} =\nabla \varphi }النتائج في

φتت+(φx2+φy2+φz2)ت=(ج2-φx2)φxx+(ج2-φy2)φyy+(ج2-φz2)φzz-2(φxφyφxy+φyφzφyz+φzφxϕzx).{\displaystyle \varphi _{tt}+(\varphi _{x}^{2}+\varphi _{y}^{2}+\varphi _{z}^{2})_{t}=(c^{2}-\varphi _{x}^{2})\varphi _{xx}+(c^{2}-\varphi _{y}^{2})\varphi _{yy}+(c^{2}-\varphi _{z}^{2})\varphi _{zz}-2(\varphi _{x}\varphi _{y}\varphi _{xy}+\varphi _{y}\varphi _{z}\varphi _{yz}+\varphi _{z}\varphi _{x}\phi _{zx}).}

التدفقات شبه المتوازية: كما في السابق، بالنسبة للتدفقات شبه المتوازية، يمكننا كتابة (بعد إدخال وقت مُعاد حسابه)τ=جت{\displaystyle \tau =c_{\infty }t})

2ϕτ2+2م2ϕxτ=(1-م2)2ϕx2+2ϕy2+2ϕz2{\displaystyle {\frac {\partial ^{2}\phi }{\partial \tau ^{2}}}+2M{\frac {\partial ^{2}\phi }{\partial x\partial \tau }}=(1-M^{2}){\frac {\partial ^{2}\phi }{\partial x^{2}}}+{\frac {\partial ^{2}\phi }{\partial y^{2}}}+{\frac {\partial ^{2}\phi }{\partial z^{2}}}}

بشرط ثبات عدد ماخم{\displaystyle M}لا يقترب من الوحدة. عندما|م-1|{\displaystyle |M-1|}إذا كانت صغيرة (تدفق فوق صوتي)، فلدينا المعادلة غير الخطية التالية [ 4 ]

2ϕτ2+22ϕxτ=-2α*ϕx2ϕx2+2ϕy2+2ϕz2.{\displaystyle {\frac {\partial ^{2}\phi }{\partial \tau ^{2}}}+2{\frac {\partial ^{2}\phi }{\partial x\partial \tau }}=-2\alpha _{*}{\frac {\partial \phi }{\partial x}}{\frac {\partial ^{2}\phi }{\partial x^{2}}}+{\frac {\partial ^{2}\phi }{\partial y^{2}}}+{\frac {\partial ^{2}\phi }{\partial z^{2}}}.}

الموجات الصوتية: في الموجات الصوتية، مقدار السرعةv{\displaystyle v}(أو عدد ماخ) صغير جدًا، على الرغم من أن الحد غير المستقر أصبح الآن قابلاً للمقارنة مع الحدود الرئيسية الأخرى في المعادلة. وبالتالي، بإهمال جميع الحدود التربيعية والحدود ذات الرتب الأعلى، وملاحظة أنه في نفس التقريب،ج{\displaystyle c}ثابت (على سبيل المثال، في الغاز متعدد الخواص)ج2=(γ-1)ح0{\displaystyle c^{2}=(\gamma -1)h_{0}}لدينا [ 8 ] [ 4 ]

2φت2=ج22φ،{\displaystyle {\frac {\partial ^{2}\varphi }{\partial t^{2}}}=c^{2}\nabla ^{2}\varphi ,}

وهي معادلة موجية خطية لكمون السرعة φ . وكما في السابق، يرتبط الجزء التذبذبي من متجه السرعة v بكمون السرعة بالعلاقة v = ∇ φ ، بينما Δ هو مؤثر لابلاس ، و c هي متوسط ​​سرعة الصوت في الوسط المتجانس . تجدر الإشارة إلى أن الأجزاء التذبذبية للضغط p والكثافة ρ ، كلٌ على حدة، تحقق المعادلة الموجية في هذا التقريب.

نطاق التطبيق والقيود

لا يشمل التدفق الكامن جميع خصائص التدفقات التي نصادفها في العالم الحقيقي. لا يمكن تطبيق نظرية التدفق الكامن على التدفقات الداخلية اللزجة ، [ 1 ] باستثناء التدفقات بين الصفائح المتقاربة . اعتبر ريتشارد فاينمان التدفق الكامن غير فيزيائي لدرجة أن السائل الوحيد الذي يفي بافتراضاته هو "الماء الجاف" (نقلاً عن جون فون نيومان ). [ 9 ] كما يقدم التدفق الكامن غير القابل للانضغاط عددًا من التنبؤات غير الصحيحة، مثل مفارقة دالمبير ، التي تنص على أن مقاومة أي جسم يتحرك عبر سائل لا نهائي ساكن تساوي صفرًا. [ 10 ] وبشكل أدق، لا يستطيع التدفق الكامن تفسير سلوك التدفقات التي تتضمن طبقة حدية . [ 1 ] ومع ذلك، يُعد فهم التدفق الكامن مهمًا في العديد من فروع ميكانيكا الموائع. على وجه الخصوص، تمتلك التدفقات الكامنة البسيطة (التي تُسمى التدفقات الأولية ) مثل الدوامة الحرة والمصدر النقطي حلولًا تحليلية جاهزة. يمكن دمج هذه الحلول لإنشاء تدفقات أكثر تعقيدًا تُحقق مجموعة متنوعة من الشروط الحدية. تُحاكي هذه التدفقات بدقة تدفقات الحياة الواقعية في جميع مجالات ميكانيكا الموائع؛ بالإضافة إلى ذلك، تبرز العديد من الأفكار القيّمة عند دراسة الانحراف (الذي غالبًا ما يكون طفيفًا) بين التدفق المرصود والتدفق الكامن المقابل. يُستخدم التدفق الكامن في العديد من التطبيقات في مجالات مثل تصميم الطائرات. على سبيل المثال، في ديناميكا الموائع الحسابية ، تتمثل إحدى التقنيات في ربط حل التدفق الكامن خارج الطبقة الحدية بحل معادلات الطبقة الحدية داخلها. يعني غياب تأثيرات الطبقة الحدية إمكانية استبدال أي خط انسيابي بحدود صلبة دون أي تغيير في مجال التدفق، وهي تقنية تُستخدم في العديد من مناهج التصميم الديناميكي الهوائي. ومن التقنيات الأخرى استخدام أجسام ريابوشينسكي الصلبة .

تحليل التدفق غير القابل للانضغاط ثنائي الأبعاد

يُعدّ تحليل التدفق الكامن في بُعدين بسيطًا باستخدام التحويلات المطابقة ، وذلك من خلال استخدام تحويلات المستوى المركب . ومع ذلك، لا يُشترط استخدام الأعداد المركبة، كما هو الحال في التحليل الكلاسيكي لتدفق الموائع حول أسطوانة. ولا يُمكن حلّ مسألة التدفق الكامن باستخدام الأعداد المركبة في ثلاثة أبعاد. [ 11 ]

تعتمد الفكرة الأساسية على استخدام دالة تحليلية (أو دالة هولومورفية ) أو دالة ميرومورفية f ، والتي تُحوّل المجال الفيزيائي ( x , y ) إلى المجال المُحوَّل ( φ , ψ ) . وبينما تكون x و y و φ و ψ جميعها ذات قيم حقيقية ، فمن الملائم تعريف الكميات المركبة.

z=x+أناy، و w=φ+أناψ.{\displaystyle {\begin{aligned}z&=x+iy\,,{\text{ and }}&w&=\varphi +i\psi \,.\end{aligned}}}

الآن، إذا كتبنا الدالة f على النحو التالي [ 11 ]

و(x+أناy)=φ+أناψ، أو و(z)=w.{\displaystyle {\begin{aligned}f(x+iy)&=\varphi +i\psi \,,{\text{ or }}&f(z)&=w\,.\end{aligned}}}

ثم، ولأن f دالة هولومورفية أو ميرومورفية، فيجب أن تحقق معادلات كوشي-ريمان [ 11 ].

φx=ψy،φy=-ψx.{\displaystyle {\begin{aligned}{\frac {\partial \varphi }{\partial x}}&={\frac {\partial \psi }{\partial y}}\,,&{\frac {\partial \varphi }{\partial y}}&=-{\frac {\partial \psi }{\partial x}}\,.\end{aligned}}}

يمكن الحصول على مركبات السرعة ( u , v ) في الاتجاهين ( x , y ) على التوالي، مباشرة من f عن طريق التفاضل بالنسبة إلى z . أي [ 11 ]

دودz=u-أناv{\displaystyle {\frac {df}{dz}}=u-iv}

لذا فإن حقل السرعة v = ( u , v ) يتم تحديده بواسطة [ 11 ]

u=φx=ψy،v=φy=-ψx.{\displaystyle {\begin{aligned}u&={\frac {\partial \varphi }{\partial x}}={\frac {\partial \psi }{\partial y}},&v&={\frac {\partial \varphi }{\partial y}}=-{\frac {\partial \psi }{\partial x}}\,.\end{aligned}}}

وبالتالي فإن كل من φ و ψ تحققان معادلة لابلاس : [ 11 ]

Δφ=2φx2+2φy2=0، و Δψ=2ψx2+2ψy2=0.{\displaystyle {\begin{aligned}\Delta \varphi &={\frac {\partial ^{2}\varphi }{\partial x^{2}}}+{\frac {\partial ^{2}\varphi }{\partial y^{2}}}=0\,,{\text{ and }}&\Delta \psi &={\frac {\partial ^{2}\psi }{\partial x^{2}}}+{\frac {\partial ^{2}\psi }{\partial y^{2}}}=0\,.\end{aligned}}}

لذا يمكن تعريف φ على أنها جهد السرعة، وتسمى ψ دالة التدفق . [ 11 ] تُعرف خطوط ψ الثابتة بخطوط التدفق ، وتُعرف خطوط φ الثابتة بخطوط تساوي الجهد (انظر سطح تساوي الجهد ).

خطوط الانسياب وخطوط تساوي الجهد متعامدة مع بعضها البعض، لأن [ 11 ]

φψ=φxψx+φyψy=ψyψx-ψxψy=0.{\displaystyle \nabla \varphi \cdot \nabla \psi ={\frac {\partial \varphi }{\partial x}}{\frac {\partial \psi }{\partial x}}+{\frac {\partial \varphi }{\partial y}}{\frac {\partial \psi }{\partial y}}={\frac {\partial \psi }{\partial y}}{\frac {\partial \psi }{\partial x}}-{\frac {\partial \psi }{\partial x}}{\frac {\partial \psi }{\partial y}}=0\,.}

وبالتالي يحدث التدفق على طول خطوط ψ الثابتة وبزوايا قائمة على خطوط φ الثابتة . [ 11 ]

يتحقق الشرط Δψ = 0 أيضًا، وهذه العلاقة مكافئة للعلاقة × v = 0. لذا فإن التدفق غير دوراني . ويعطي الشرط التلقائي ∂²ψ / ∂x / ∂y = ∂²ψ / ∂y / ∂x قيد عدم الانضغاط· v = 0 .

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

يمكن استخدام أي دالة قابلة للتفاضل لـ f . تستخدم الأمثلة التالية مجموعة متنوعة من الدوال الأولية ؛ كما يمكن استخدام دوال خاصة . تجدر الإشارة إلى أنه يمكن استخدام دوال متعددة القيم مثل اللوغاريتم الطبيعي ، ولكن يجب التركيز على سطح ريمان واحد .

قوانين القوة

أمثلة على التحويلات المطابقة لقانون القوة w = Az n
أمثلة على التحويلات المطابقة لقانون القوة w = Az n ، لقيم مختلفة للقوة n . يظهر المستوى z ، موضحًا خطوط الجهد الثابت φ ودالة التدفق ψ ، بينما w = φ + .

في حالة تطبيق التحويل المطابق لقانون القوة التالي ، من z = x + iy إلى w = φ + : [ 12 ]

w=أzن،{\displaystyle w=Az^{n}\,,}

ثم، بكتابة z في الإحداثيات القطبية على النحو التالي: z = x + iy = re ، لدينا [ 12 ]

φ=أرنكوسنθوψ=أرنالخطيئةنθ.{\displaystyle \varphi =Ar^{n}\cos n\theta \qquad {\text{and}}\qquad \psi =Ar^{n}\sin n\theta \,.}

تُظهر الأشكال على اليمين أمثلةً لعدة قيم لـ n . الخط الأسود هو حدود التدفق، بينما الخطوط الزرقاء الداكنة هي خطوط الانسياب، والخطوط الزرقاء الفاتحة هي خطوط تساوي الجهد. بعض قوى n المهمة هي: [ 12 ]

  • n = 1 / 2 : هذا يتوافق مع التدفق حول صفيحة شبه لانهائية ،
  • n = 2 / 3 : التدفق حول زاوية يمين ،
  • n = 1 : حالة بسيطة من التدفق المنتظم،
  • n = 2 : التدفق عبر زاوية، أو بالقرب من نقطة ركود، و
  • n = −1 : التدفق الناتج عن مصدر مزدوج

الثابت A هو معلمة قياس: قيمته المطلقة | A | تحدد المقياس، بينما وسيطه arg( A ) يُدخل دورانًا (إذا كان غير صفري).

قوانين القوة مع n = 1 : تدفق منتظم

إذا كانت w = Az 1 ، أي قانون قوة حيث n = 1 ، فإن خطوط الانسياب (أي خطوط ψ الثابتة ) هي نظام من الخطوط المستقيمة الموازية للمحور x . ويسهل فهم ذلك بالكتابة بدلالة الأجزاء الحقيقية والخيالية.

و(x+أناy)=أ(x+أناy)=أx+أناأy{\displaystyle f(x+iy)=A\,(x+iy)=Ax+iAy}

وبذلك نحصل على φ = Ax و ψ = Ay . ويمكن تفسير هذا التدفق على أنه تدفق منتظم موازٍ للمحور x .

قوانين القوة مع n = 2

إذا كانت n = 2 ، فإن w = Az وخطوط التدفق المقابلة لقيمة معينة لـ ψ هي تلك النقاط التي تحقق

ψ=أر2الخطيئة2θ،{\displaystyle \psi =Ar^{2}\sin 2\theta \,,}

وهو نظام من القطوع الزائدة المستطيلة . ويمكن ملاحظة ذلك بإعادة كتابة المعادلة بدلالة الأجزاء الحقيقية والخيالية. بملاحظة أن sin²θ = 2 sinθ cosθ ، وإعادة كتابة sinθ = y / r و cosθ = x / r ، يتضح (بعد التبسيط) أن خطوط الانسياب تُعطى بالعلاقة التالية :

ψ=2أxy.{\displaystyle \psi =2Axy\,.}

يُعطى حقل السرعة بواسطة φ ، أو

(uv)=(φxφy)=(+ψy-ψx)=(+2أx-2أy).{\displaystyle {\begin{pmatrix}u\\v\end{pmatrix}}={\begin{pmatrix}{\frac {\partial \varphi }{\partial x}}\\[2px]{\frac {\partial \varphi }{\partial y}}\end{pmatrix}}={\begin{pmatrix}+{\partial \psi \over \partial y}\\[2px]-{\partial \psi \over \partial x}\end{pmatrix}}={\begin{pmatrix}+2Ax\\[2px]-2Ay\end{pmatrix}}\,.}

في ديناميكا الموائع، يُمثل مجال التدفق بالقرب من نقطة الأصل نقطة ركود . لاحظ أن المائع عند نقطة الأصل يكون ساكنًا (وهذا ناتج عن اشتقاق الدالة f (z) = عند z = 0 ). يُعد خط الانسياب ψ = 0 مثيرًا للاهتمام بشكل خاص: إذ يتفرع إلى فرعين (أو أربعة فروع) باتجاه محوري الإحداثيات، أي x = 0 و y = 0. وبما أنه لا يوجد تدفق للمائع عبر المحور x ، فيمكن اعتباره حدًا صلبًا. وبالتالي، يُمكن تجاهل التدفق في النصف السفلي من المستوى حيث y < 0 والتركيز على التدفق في النصف العلوي. وفقًا لهذا التفسير، يكون التدفق عبارة عن نفاثة موجهة رأسيًا تصطدم بلوحة أفقية مسطحة. كما يُمكن تفسير التدفق على أنه تدفق باتجاه زاوية قائمة (90 درجة) إذا تم تجاهل المناطق المحددة بـ (على سبيل المثال) x ، y < 0 .

قوانين القوة مع n = 3

إذا كانت n = 3 ، فإن التدفق الناتج يكون أشبه بنسخة سداسية من حالة n = 2 المذكورة أعلاه. تُعطى خطوط الانسياب بالمعادلة ψ = 3 x 2 yy 3 ، ويمكن تفسير التدفق في هذه الحالة على أنه تدفق نحو زاوية 60 درجة.

قوانين القوة مع n = −1 : ثنائية

إذا كانت n = −1 ، فإن خطوط الانسياب تُعطى بواسطة

ψ=-أرالخطيئةθ.{\displaystyle \psi =-{\frac {A}{r}}\sin \theta .}

يمكن تفسير ذلك بسهولة أكبر من حيث المكونات الحقيقية والخيالية: ψ=-أyر2=-أyx2+y2،x2+y2+أyψ=0،x2+(y+أ2ψ)2=(أ2ψ)2.{\displaystyle {\begin{aligned}\psi ={\frac {-Ay}{r^{2}}}&={\frac {-Ay}{x^{2}+y^{2}}}\,,\\x^{2}+y^{2}+{\frac {Ay}{\psi }}&=0\,,\\x^{2}+\left(y+{\frac {A}{2\psi }}\right)^{2}&=\left({\frac {A}{2\psi }}\right)^{2}\,.\end{aligned}}}

وبالتالي، فإن خطوط الانسياب عبارة عن دوائر مماسية للمحور السيني عند نقطة الأصل. لذا، تتدفق الدوائر في النصف العلوي من المستوى في اتجاه عقارب الساعة، بينما تتدفق تلك الموجودة في النصف السفلي عكس اتجاه عقارب الساعة. لاحظ أن مركبات السرعة تتناسب مع r⁻² ، وقيمتها عند نقطة الأصل لانهائية. يُشار عادةً إلى نمط التدفق هذا باسم " ثنائي القطب" ، ويمكن تفسيره على أنه مزيج من مصدر ومستقبل ذي قوة لانهائية، يفصل بينهما مسافة متناهية الصغر. يُعطى حقل السرعة بالعلاقة التالية:

(u،v)=(ψy،-ψx)=(أy2-x2(x2+y2)2،-أ2xy(x2+y2)2).{\displaystyle (u,v)=\left({\frac {\partial \psi }{\partial y}},-{\frac {\partial \psi }{\partial x}}\right)=\left(A{\frac {y^{2}-x^{2}}{\left(x^{2}+y^{2}\right)^{2}}},-A{\frac {2xy}{\left(x^{2}+y^{2}\right)^{2}}}\right)\,.}

أو بالإحداثيات القطبية:

(uر،uθ)=(1رψθ،-ψر)=(-أر2كوسθ،-أر2الخطيئةθ).{\displaystyle (u_{r},u_{\theta })=\left({\frac {1}{r}}{\frac {\partial \psi }{\partial \theta }},-{\frac {\partial \psi }{\partial r}}\right)=\left(-{\frac {A}{r^{2}}}\cos \theta ,-{\frac {A}{r^{2}}}\sin \theta \right)\,.}

قوانين القوى مع n = −2 : رباعي الأقطاب

إذا كانت n = −2 ، فإن خطوط الانسياب تُعطى بواسطة

ψ=-أر2الخطيئة2θ.{\displaystyle \psi =-{\frac {A}{r^{2}}}\sin 2\theta \,.}

هذا هو مجال التدفق المرتبط بـ رباعي الأقطاب . [ 13 ]

مصدر ومصب الخط

مصدر أو مصب للقوةسؤال{\displaystyle Q}(سؤال>0{\displaystyle Q>0}للحصول على المصدر وسؤال<0{\displaystyle Q<0}يُعطى (للحوض) بواسطة الجهد

w=سؤال2πlnz{\displaystyle w={\frac {Q}{2\pi }}\ln z}

أينسؤال{\displaystyle Q}في الواقع، يمثل التدفق الحجمي لكل وحدة طول عبر سطح يحيط بالمصدر أو المصرف. حقل السرعة في الإحداثيات القطبية هو

uر=سؤال2πر،uθ=0{\displaystyle u_{r}={\frac {Q}{2\pi r}},\quad u_{\theta }=0}

أي، تدفق شعاعي بحت.

دوامة خطية

دوامة خطية من القوةΓ{\displaystyle \Gamma }يُعطى بواسطة

w=Γ2πأناlnz{\displaystyle w={\frac {\Gamma }{2\pi i}}\ln z}

أينΓ{\displaystyle \Gamma }هي الدوران حول أي محيط مغلق بسيط يحيط بالدوامة. حقل السرعة في الإحداثيات القطبية هو

uر=0،uθ=Γ2πر{\displaystyle u_{r}=0,\quad u_{\theta }={\frac {\Gamma }{2\pi r}}}

أي، تدفق دائري بحت.

تحليل التدفقات غير القابلة للانضغاط ثلاثية الأبعاد

بالنسبة للتدفقات ثلاثية الأبعاد، لا يمكن الحصول على الجهد المركب.

مصادر ومصارف نقطية

جهد السرعة لمصدر نقطي أو مصب ذي قوةسؤال{\displaystyle Q}(سؤال>0{\displaystyle Q>0}للحصول على المصدر وسؤال<0{\displaystyle Q<0}بالنسبة للحوض، يُعطى في الإحداثيات القطبية الكروية بالصيغة التالية:

ϕ=-سؤال4πر{\displaystyle \phi =-{\frac {Q}{4\pi r}}}

أينسؤال{\displaystyle Q}في الواقع، يمثل التدفق الحجمي عبر سطح مغلق يحيط بالمصدر أو المصرف. حقل السرعة في الإحداثيات القطبية الكروية هو

uر=سؤال4πر2،uθ=0،uϕ=0.{\displaystyle u_{r}={\frac {Q}{4\pi r^{2}}},\quad u_{\theta }=0,\quad u_{\phi }=0.}

انظر أيضاً

ملحوظات

  1. 1 2 3 باتشيلور (1973) ص. 378–380.
  2. كيربي، بي جيه (2010)، ميكانيكا الموائع على المستويين الميكروي والنانوي: النقل في الأجهزة الميكروفلويدية. ، مطبعة جامعة كامبريدج، رقم ISBN 978-0-521-11903-0
  3. 1 2 3 باتشيلور (1973) ص. 99–101.
  4. 1 2 3 4 5 لاندو، إل دي؛ ليفشيتز، إي إم (2013). "114". ميكانيكا الموائع . لاندو وليفشيتز: دورة في الفيزياء النظرية. المجلد 6. إلسيفير. ص 436.  
  5. أندرسون، جيه دي (2002). التدفق الانضغاطي الحديث . ماكجرو هيل. ص 358-359 . ISBN  0-07-242443-5.
  6. 1942، لاندو، إل دي "حول الموجات الصدمية" مجلة الفيزياء في الاتحاد السوفيتي 6 229-230
  7. تومسون، ب.أ. (1971). مشتق أساسي في ديناميكا الغازات. فيزياء السوائل، 14(9)، 1843-1849.
  8. لامب (1994) §287، ص 492-495.
  9. فاينمان، آر بي ؛ لايتون، آر بي ؛ ساندز، إم. (1964)، محاضرات فاينمان في الفيزياء ، المجلد 2، أديسون-ويسلي ، ص 40-3. الفصل 40 يحمل العنوان: تدفق المياه الجافة .
  10. باتشيلور (1973) ص. 404–405.
  11. 1 2 3 4 5 6 7 8 9 باتشيلور (1973) ص. 106–108.
  12. 1 2 3 باتشيلور (1973) ص. 409–413.
  13. كيرالا، أ. (1972). الدوال التطبيقية لمتغير مركب . وايلي-إنترساينس. ص 116-117 . ISBN  9780471511298.

مراجع

  • باتشيلور، جي كي (1973)، مقدمة في ديناميكا الموائع ، مطبعة جامعة كامبريدج، رقم ISBN 0-521-09817-3
  • شانسون، هـ. (2009)، الديناميكا المائية التطبيقية: مقدمة في تدفقات الموائع المثالية والحقيقية ، مطبعة سي آر سي، مجموعة تايلور وفرانسيس، ليدن، هولندا، 478 صفحة، رقم ISBN 978-0-415-49271-3
  • لامب، هـ. (1994) [1932]، الديناميكا المائية (  الطبعة السادسة)، مطبعة جامعة كامبريدج، رقم ISBN 978-0-521-45868-9
  • ميلن-تومسون، إل إم (1996) [1968]، الديناميكا المائية النظرية (  الطبعة الخامسة)، دوفر، رقم ISBN 0-486-68970-0

للمزيد من القراءة

  • Chanson، H. (2007)، “Le Pottiel de vitesse pour les écoulements de Fluides réels: la مساهمة جوزيف لويس لاغرانج [إمكانية السرعة في تدفقات السوائل الحقيقية: مساهمة جوزيف لويس لاغرانج]” ، لا هويل بلانش (بالفرنسية)، 93 (5): 127– 131، بيب كود : 2007LHBl...93..127C ، دوى : 10.1051/lhb:2007072
  • ويهاوزن، جيه في ؛ لايتون، إي في (1960)، "موجات السطح"، في فلوغه، إس ؛ تروسديل، سي (محرران)، موسوعة الفيزياء ، المجلد  التاسع، سبرينغر فيرلاغ، الصفحات 446-778 ، مؤرشف من الأصل في 5 يناير 2009 ، تم استرجاعه في 29 مارس 2009