المعادلات الحاكمة الأساسية هي معادلة الاستمرارية غير القابلة للانضغاط ومعادلات الزخم لنافيير-ستوكس. بالنسبة لمائع نيوتوني غير قابل للانضغاط ومتساوي الحرارة بكثافة ثابتة (ρ) ولزوجة ديناميكية ثابتة (μ)، تكون المعادلات كالتالي:
[ \nabla \cdot \mathbf{v}=0 ]
و
[ \rho\left(\frac{\partial \mathbf{v}}{\partial t} +\mathbf{v}\cdot\nabla\mathbf{v}\right)
-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]
تحدد هذه المعادلات معاً مكونات السرعة الثلاث ومجال الضغط، وذلك بناءً على هندسة المعدات المناسبة، وظروف المدخل والمخرج والجدران والظروف الابتدائية.
الخلاصة الجوهرية: تتم نمذجة التدفق النيوتوني غير القابل للانضغاط ومتساوي الحرارة في معدات المحطات التجريبية بشكل أساسي من خلال حفظ الكتلة وحفظ الزخم. لا يتم إضافة معادلات الطاقة أو الأنواع الكيميائية أو التفاعلات أو الأطوار المتعددة إلا عندما تتضمن وحدة التشغيل تأثيرات لا يمكن تمثيلها بهاتين المعادلتين وحدهما.
مجموعة المعادلات الحاكمة
حفظ الكتلة: معادلة الاستمرارية
معادلة حفظ الكتلة العامة هي:
[ \frac{\partial \rho}{\partial t}+\nabla\cdot(\rho\mathbf{v})=0. ]
بالنسبة لمائع غير قابل للانضغاط ذي كثافة ثابتة،
[ \boxed{\nabla\cdot\mathbf{v}=0} ]
ينص هذا الشرط على أن عنصر المائع لا يتمدد ولا ينكمش أثناء حركته. في الإحداثيات الديكارتية،
[ \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z} =0, ]
حيث (u) و(v) و(w) هي مكونات السرعة في اتجاهات (x) و(y) و(z).
حفظ الزخم: معادلات نافيير-ستوكس
توازن الزخم الخطي هو:
[ \boxed{ \rho\frac{D\mathbf{v}}{Dt}
-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g} } ]
حيث المشتقة المادية هي:
[ \frac{D\mathbf{v}}{Dt}
\frac{\partial\mathbf{v}}{\partial t} + (\mathbf{v}\cdot\nabla)\mathbf{v}. ]
وبالتالي، تمثل معادلة الزخم توازناً بين:
- تسارع المائع: (\rho D\mathbf{v}/Dt)
- قوى الضغط: (-\nabla p)
- قوى اللزوجة: (\mu\nabla^2\mathbf{v})
- قوى الجسم: (\rho\mathbf{g})، وعادة ما تكون الجاذبية
يعد الحد الحملي ((\mathbf{v}\cdot\nabla)\mathbf{v}) مهماً في الانحناءات، والتقلصات، والتوسعات، والخلاطات، والمناطق الأخرى التي تتغير فيها السرعة مكانياً.
صيغة المكونات الديكارتية
بالنسبة للسرعة (\mathbf{v}=(u,v,w))، تكون معادلات الزخم:
[ \rho\left( \frac{\partial u}{\partial t} +u\frac{\partial u}{\partial x} +v\frac{\partial u}{\partial y} +w\frac{\partial u}{\partial z} \right)
-\frac{\partial p}{\partial x} +\mu\nabla^2u+\rho g_x, ]
[ \rho\left( \frac{\partial v}{\partial t} +u\frac{\partial v}{\partial x} +v\frac{\partial v}{\partial y} +w\frac{\partial v}{\partial z} \right)
-\frac{\partial p}{\partial y} +\mu\nabla^2v+\rho g_y, ]
[ \rho\left( \frac{\partial w}{\partial t} +u\frac{\partial w}{\partial x} +v\frac{\partial w}{\partial y} +w\frac{\partial w}{\partial z} \right)
-\frac{\partial p}{\partial z} +\mu\nabla^2w+\rho g_z. ]
هذه المعادلات الثلاث، مجتمعة مع الاستمرارية، تعطي أربع معادلات لأربعة مجالات أساسية: (u) و(v) و(w) و(p).
لماذا يهم السلوك النيوتوني؟
قانون نيوتن للزوجة
يحتوي المائع النيوتوني على إجهادات لزوجية تتغير خطياً مع معدل التشوه. في صيغة الموتر،
[ \boldsymbol{\tau}
2\mu\mathbf{D}, ]
للتدفق غير القابل للانضغاط، حيث
[ \mathbf{D}
\frac{1}{2} \left[ \nabla\mathbf{v} + (\nabla\mathbf{v})^{T} \right] ]
هو موتر معدل التشوه.
بالنسبة للزوجة الثابتة والتدفق غير القابل للانضغاط، يختزل تباعد الإجهاد اللزوجي إلى
[ \nabla\cdot\boldsymbol{\tau}
\mu\nabla^2\mathbf{v}. ]
هذا الاختزال هو ما ينتج الحد اللزوجي في معادلة نافيير-ستوكس القياسية.
التفسير الفيزيائي
يصف الحد اللزوجي انتشار الزخم. وهو يفسر انتقال الزخم الناتج عن تدرجات السرعة، مثل تلك الموجودة بين جدار الأنبوب الثابت والمائع المتحرك في مركز الأنبوب.
الشروط اللازمة لإغلاق النموذج
شرط عدم الانزلاق عند الجدار
عند جدار صلب ثابت،
[ \mathbf{v}=0. ]
هذا الشرط أساسي للتنبؤ بإجهاد قص الجدار، وفقدان الضغط، ومقاطع السرعة في الأنابيب والقنوات والمبادلات الحرارية والمعدات المعبأة.
بالنسبة للجدار المتحرك، تساوي سرعة المائع عند الجدار سرعة الجدار.
شروط المدخل
يمكن تحديد المدخل من خلال مقطع سرعة،
[ \mathbf{v}=\mathbf{v}_{\text{in}}, ]
أو من خلال معدل تدفق الكتلة،
[ \dot{m}
\int_A \rho,\mathbf{v}\cdot\mathbf{n},dA. ]
غالباً ما يُستخدم مقطع مدخل موحد للدراسات المبسطة، بينما تكون المقاطع المتطورة أو المقاسة أكثر ملاءمة لتحليل المحطات التجريبية التفصيلي.
شروط المخرج
تشمل مواصفات المخرج النموذجية ضغطاً معروفاً،
[ p=p_{\text{out}}, ]
أو شرط تدفق خارج مناسب عندما يكون التدفق في اتجاه المصب متطوراً بشكل كافٍ. يمكن أن يؤثر شرط المخرج المختار بشكل مادي على مجالات الضغط المتوقعة وإعادة التدوير.
الشرط الابتدائي
بالنسبة للتشغيل العابر (غير المستقر)، يجب تحديد مجال السرعة البدائي:
[ \mathbf{v}(\mathbf{x},0)=\mathbf{v}_0(\mathbf{x}). ]
مطلوب أيضاً مرجع للضغط لأن ضغط المائع غير القابل للانضغاط يتم تحديده عموماً فقط حتى ثابت اختياري.
التبسيطات الشائعة في وحدات المحطات التجريبية
التدفق المستقر
إذا كانت ظروف التشغيل لا تتغير بشكل كبير مع الزمن،
[ \frac{\partial\mathbf{v}}{\partial t}=0. ]
تصبح معادلة الزخم
[ \rho(\mathbf{v}\cdot\nabla)\mathbf{v}
-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]
يُستخدم هذا التقريب بشكل شائع لتحليلات الأنابيب والقنوات والمبادلات الحرارية المستقرة.
التدفق الداخلي المتطور بالكامل
بالنسبة للتدفق المتطور بالكامل في قناة أو أنبوب مستقيم، لا يتغير مقطع السرعة في الاتجاه المحوري. يمكن بعد ذلك اختزال المعادلات الحاكمة إلى توازن بين تدرج الضغط والمقاومة اللزوجية.
على سبيل المثال، بالنسبة للتدفق أحادي البعد بين لوحين متوازيين،
[ 0=-\frac{dp}{dx}+\mu\frac{d^2u}{dy^2}. ]
يكون مقطع السرعة الناتج مكافئاً (parabolic) عندما يكون تدرج الضغط ثابتاً.
القصور الذاتي المهمل
عند رقم رينولدز منخفض بما فيه الكفاية، تكون تأثيرات القصور الذاتي صغيرة مقارنة بالتأثيرات اللزوجية. تختزل معادلة الزخم إلى شكل تدفق ستوكس،
[ 0=-\nabla p+\mu\nabla^2\mathbf{v}+\rho\mathbf{g}. ]
هذا التقريب مفيد للموائع عالية اللزوجة، والقنوات الضيقة، والوحدات الدقيقة أو المختبرية ذات معدل التدفق المنخفض.
الجاذبية والضغط الهيدروستاتيكي
إذا كانت السرعة مهملة، تختزل معادلة الزخم إلى
[ \nabla p=\rho\mathbf{g}. ]
هذا يصف تغير الضغط الهيدروستاتيكي. في أنظمة التدفق، قد تظل الجاذبية مهمة في الأنابيب الرأسية، والأعمدة المعبأة، والخزانات، والتكوينات المتعلقة بالحمل الحراري الطبيعي.
كيفية تطبيق المعادلات على وحدات التشغيل
شبكات الأنابيب والمضخات
تصف المعادلات توزيعات السرعة وتدرجات الضغط في جميع أنحاء الأنابيب. تُستخدم بعد ذلك الصيغ المتكاملة أو المختزلة لتقدير فقدان الضغط وتحديد متطلبات المضخة.
بالنسبة لمحطة تجريبية، يمكن حل معادلات المجال الكاملة حاسوبياً، بينما غالباً ما يتم تمثيل أقسام الأنابيب الأبسط من خلال ارتباطات معامل الاحتكاك المشتقة من هذه التوازنات.
المبادلات الحرارية
بالنسبة لنمذجة تدفق الموائع متساوية الحرارة بدقة، تحدد الاستمرارية والزخم الديناميكا المائية، بما في ذلك توزيع السرعة، وقص الجدار، وانخفاض الضغط.
إذا كان الهدف يشمل تغير درجة الحرارة أو واجب نقل الحرارة، فيجب أيضاً إضافة معادلة الطاقة. توازن التيار المشترك هو
[ Q=\dot{m}C_p\Delta T, ]
بينما غالباً ما يرتبط أداء المبادل من خلال
[ Q=UA\Delta T_{\mathrm{lm}}F. ]
هذه علاقات حرارية إضافية، وليست بدائل لمعادلات التدفق.
المفاعلات ومعدات الخلط
في الخزانات المحركة ومفاعلات التدفق، تتنبأ معادلات نافيير-ستوكس بالدوران، وإعادة التدوير، ومجالات السرعة المتعلقة بالخلط.
تتطلب التفاعلات، ونقل الأنواع، وسلوك زمن البقاء معادلات حفظ إضافية. لا تتنبأ معادلات تدفق الموائع وحدها بالتحويل أو التركيز.
الأسرة والأعمدة المعبأة
في المعدات المعبأة، يمكن استخدام معادلات نافيير-ستوكس على نطاق المسام للنمذجة التفصيلية، ولكن نماذج النطاق الهندسي تستبدلها عادةً بعلاقات انخفاض ضغط الوسائط المسامية.
تظل المبادئ الأساسية هي حفظ الكتلة والزخم، مع تمثيل التعبئة من خلال حدود مقاومة إضافية.
وحدات الموائع الثنائية أو متعددة الأطوار
بالنسبة لمائعين غير قابلين للامتزاج، يتطلب كل طور معادلات الحفظ الخاصة به، إلى جانب شروط الواجهة. عند واجهة سائل-سائل، تشمل الشروط النموذجية
[ \mathbf{v}_1=\mathbf{v}_2 ]
لاستمرارية السرعة المماسية و
[ \boldsymbol{\tau}_1\cdot\mathbf{n}
\boldsymbol{\tau}_2\cdot\mathbf{n} ]
لاستمرارية الإجهاد اللزوجي، عندما يتم التعامل مع تأثيرات الواجهة بشكل مناسب.
لذلك، تكون معادلات المائع الواحد غير القابل للانضغاط غير كافية عندما تكون حدود الطور، أو انتقال الزخم بين الأطوار، أو تغير الطور مهمة.
فهم المقايضات
متساوي الحرارة لا يعني أن الطاقة غير ذات صلة دائماً
يفترض النموذج متساوي الحرارة أن درجة الحرارة ثابتة أو أن التأثيرات الحرارية لا تؤثر على حل التدفق. وهو لا يوفر توزيعات درجات الحرارة، أو معاملات نقل الحرارة، أو الواجب الحراري.
بالنسبة للمبادلات الحرارية، أو المفاعلات المسخنة، أو الأنظمة ذات التغير القوي في اللزوجة، قد تكون معادلة الطاقة ضرورية حتى لو استخدم تحليل التدفق الأولي تقريباً متساوي الحرارة.
عدم الانضغاط له نطاق محدد من الصلاحية
غالباً ما تتم نمذجة السوائل على أنها غير قابلة للانضغاط لأن تغيرات الكثافة صغيرة. يمكن أيضاً التعامل مع الغازات على أنها غير قابلة للانضغاط عندما تكون السرعة أقل بكثير من سرعة الصوت وتكون تغيرات الكثافة الناتجة عن الضغط مهملة.
يجب إعادة النظر في التقريب لتدفق الغاز عالي السرعة، أو تغيرات الضغط الكبيرة، أو التغير القوي في درجة الحرارة، أو الأنظمة التي تنطوي على تأثيرات مدفوعة بالانضغاط.
اللزوجة الثابتة قد تكون تقريباً ضعيفاً
تفترض المعادلة القياسية (μ) ثابتة. هذا معقول للعديد من السوائل عبر نطاقات معتدلة من درجة الحرارة ومعدل القص، ولكن ليس للموائع شديدة الاعتماد على درجة الحرارة أو المواد غير النيوتونية.
إذا تغيرت اللزوجة مكانياً، فيجب عموماً الاحتفاظ بالحد اللزوجي في شكل تباعد الإجهاد بدلاً من كتابته ببساطة كـ (μ∇²v).
افتراضات التدفق الصفائحي ليست مطلوبة من قبل نافيير-ستوكس
تنطبق معادلات نافيير-ستوكس على التدفق النيوتوني سواء كان صفائحياً أو مضطرباً. ومع ذلك، فإن الحل المباشر لمقاييس الاضطراب يتطلب جهداً حاسوبياً، لذا غالباً ما تكون نماذج الاضطراب أو الارتباطات التي تم التحقق منها تجريبياً مطلوبة للمعدات على نطاق تجريبي.
"تبدد اللزوجة" ينتمي إلى توازن الطاقة
تؤثر الإجهادات اللزوجية على معادلة الزخم مباشرة. تظهر الحرارة الناتجة عن تبدد اللزوجة في معادلة الطاقة ويتم حذفها من صياغة متساوية الحرارة بدقة.
يمنع هذا التمييز الخطأ الشائع المتمثل في التعامل مع تبدد اللزوجة كحد إضافي في معادلة الزخم.
اتخاذ الخيار الصحيح لهدفك
يعتمد النموذج الحاكم الصحيح على ما يجب أن يتنبأ به تحليل المحطة التجريبية.
- إذا كان تركيزك الأساسي هو مقاطع السرعة وفقدان الضغط: استخدم استمرارية غير قابلة للانضغاط بالإضافة إلى معادلات زخم نافيير-ستوكس، مع كثافة ولزوجة مناسبة لظروف التشغيل.
- إذا كان تركيزك الأساسي هو أداء نقل الحرارة: أضف معادلة طاقة واستخدم حل التدفق لتحديد الحمل الحراري، وانخفاض الضغط، وسلوك نقل الحرارة.
- إذا كان تركيزك الأساسي هو تحويل المفاعل أو توزيع الأنواع: أضف توازنات الأنواع أو التفاعلات إلى معادلات الديناميكا المائية.
- إذا كان تركيزك الأساسي هو تشغيل سائل-سائل أو غاز-سائل: استخدم معادلات حفظ خاصة بكل طور وقم بتضمين شروط السرعة البينية، والإجهاد، وانتقال الكتلة.
- إذا كان تركيزك الأساسي هو التصميم السريع على نطاق المعدات: استخدم نماذج مختزلة تم التحقق منها، مثل ارتباطات معامل الاحتكاك، أو الوسائط المسامية، أو زمن البقاء، مع إدراك حدود تشغيلها.
مع ذكر الافتراضات بوضوح وتحديد الشروط الحدودية بشكل صحيح، توفر معادلات الحفظ هذه أساساً موثوقاً لتحليل سلوك تدفق الموائع في المحطات التجريبية.
جدول الملخص:
| عنصر النموذج | العلاقة الحاكمة | الغرض الأساسي |
|---|---|---|
| حفظ الكتلة | ∇ · v = 0 | يضمن استمرارية التدفق غير القابل للانضغاط |
| حفظ الزخم | ρ(∂v/∂t + v · ∇v) = −∇p + μ∇²v + ρg | يتنبأ بتأثيرات السرعة، والضغط، واللزوجة، وقوى الجسم |
| اللزوجة النيوتونية | τ = 2μD | يربط الإجهاد اللزوجي بمعدل التشوه |
| الشروط الحدودية | المدخل، المخرج، الجدار، والشروط الابتدائية | يغلق النموذج الرياضي |
| فيزياء إضافية | معادلات الطاقة، الأنواع، التفاعل، أو الأطوار المتعددة | يوسع النموذج لنقل الحرارة، أو التفاعلات، أو الأطوار المتعددة |
هل تحتاج إلى معدات موثوقة للمحطات التجريبية لتعليم تدفق الموائع، أو تطوير العمليات، أو توسيع النطاق؟ اتصل بـ LABPARK لمناقشة متطلباتك. توفر LABPARK محطات تجريبية لوحدات التشغيل التعليمية والمهنية للهندسة الكيميائية، والعمليات الحيوية والتكنولوجيا الحيوية، وتطبيقات معالجة البيئة والمياه. تساعد أنظمتنا النمطية الجامعات ومعاهد البحوث والمؤسسات على دراسة انخفاض الضغط، ومقاطع السرعة، والخلط، ونقل الحرارة، والفصل، والتفاعل، والتحكم في العمليات باستخدام معدات عملية تركز على التطبيق. أخبرنا بعمليتك المستهدفة، وظروف التشغيل، وأهداف التدريس أو البحث، والمقياس المطلوب حتى يتمكن فريقنا من التوصية بتكوين مناسب ودعم تخطيط محطتك التجريبية.
المنتجات ذات الصلة
- وحدة عمليات التدريب النموذجية لإنتاج مستحضرات التجميل للأغراض العامة
- وحدة تدريب عمليات وحدات نقل السوائل بمضخات متعددة منشأة تجريبية
- وحدة تجارب تعليمية متعددة الوظائف للتجفيف
- وحدة تجريبية لعمليات التفاعل الكيميائي في طبقة ثابتة وإزالة الغبار والقطران من الغاز
- محطة تجريبية لعمليات الهدرجة التعليمية المستمرة بالدورة المغلقة بسعة 100 لتر
يسأل الناس أيضًا
- كيف يمكن لمحطات تجريبية لعمليات الوحدات إثبات الإزالة المتزامنة للشوائب الخفيفة والثقيلة؟ قم بتهيئة إعداد معياري.
- كيف تكمل المحاكاة المستقرة والديناميكية عمل وحدات التشغيل التجريبية؟ التحقق من التوسيع وقابلية التحكم
- لماذا يُرجع محلّل محطة تجريبية خطأ المصفوفة المفردة، وكيف يمكن حلّه؟
- كيف تؤثر الانحرافات في تقدير الحرارة الكامنة على الأنظمة الحرارية لمحطة التجربة؟ تجنب تحديد أحجام الأجهزة بشكل غير صحيح.
- لماذا يعتبر جدول بدء تشغيل المصنع الكيميائي أمرًا حاسمًا؟ تقليل مخاطر التوسع باستخدام مصانع النماذج الأولية.