حاسبة تدفق الأنابيب

تحليل تدفق الأنابيب أمر أساسي لتصميم أنظمة السباكة والتكييف والتبريد والأنظمة الصناعية الكفؤة.

تشمل المتغيرات الرئيسية ما يلي:

• معدل التدفق
• السرعة
• انخفاض الضغط
• نظام التدفق

التدفق الصفائحي (Re < 2300) يمكن التنبؤ به ولكنه نادر في التطبيق العملي.

التدفق المضطرب (Re > 4000) شائع ويوفر خلطًا أفضل.

ينبغي تجنب المنطقة الانتقالية لضمان سلوك متوقع للنظام.

معدل التدفق الحجمي يساوي مساحة المقطع العرضي الداخلي للأنبوب مضروبة في متوسط سرعة التدفق: Q = A x v، حيث A = pi x d^2 / 4. بالنسبة لأنبوب دائري بقطر داخلي d وسرعة v، تصبح المعادلة Q = (pi/4) x d^2 x v.

مثال: أنبوب بقطر 50 مم (0.05 م) ينقل الماء بسرعة 2 م/ث تكون مساحته A = pi x 0.05^2 / 4 = 0.001963 m^2، بالتالي Q = 0.001963 x 2 = 0.003927 m^3/s، وهو ما يعادل 3.93 لتر/ثانية أو حوالي 62 جالونًا أمريكيًا في الدقيقة.

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

في وحدات النظام الدولي (SI):

• يُقاس معدل التدفق الحجمي بالمتر المكعب في الثانية (m^3/s)
• تُقاس السرعة بالمتر في الثانية (m/s)
• يُقاس الضغط وانخفاض الضغط بالباسكال (Pa)
• يُقاس قطر الأنبوب بالمتر

يُعرَّف الباسكال بواسطة BIPM وNIST بأنه نيوتن واحد لكل متر مربع (1 Pa = 1 N/m^2). تُفضّل الهندسة العملية اللتر في الثانية (L/s) أو الجالون الأمريكي في الدقيقة (GPM) للتدفق، والكيلو باسكال (kPa) أو البار للضغط، والمليمتر للقطر.

تحويلات مفيدة:

• 1 m^3/s = 1000 L/s
• 1 L/s = 15.85 US GPM
• 1 kPa = 1000 Pa

رقم رينولدز ومعامل احتكاك دارسي هما كميتان لا بعديتان، لذا تظل قيمتهما كما هي في أي نظام وحدات متسق.

رقم رينولدز هو نسبة لا بعدية بين قوى القصور الذاتي وقوى اللزوجة، ويُحسب بالمعادلة Re = (rho x v x d) / mu، حيث rho هي كثافة المائع (kg/m^3)، وv هي السرعة (m/s)، وd هو القطر الداخلي (m)، وmu هي اللزوجة الديناميكية (Pa.s).

بالنسبة للماء عند درجة حرارة 20 C، تكون rho حوالي 998 kg/m^3 وتكون mu حوالي 0.001002 Pa.s.

مثال: الماء بسرعة 2 م/ث في أنبوب بقطر 50 مم يعطي Re = (998 x 2 x 0.05) / 0.001002 = 99,600.

كما يوضح مشروع HyperPhysics بجامعة ولاية جورجيا، فإن قيمة Re أقل من 2,300 تشير إلى تدفق صفائحي، وأكبر من 4,000 تشير إلى تدفق مضطرب، بينما النطاق بينهما انتقالي وغير متوقع.

يُحسب انخفاض الضغط الناتج عن احتكاك الأنابيب بـ معادلة دارسي-وايسباخ: delta_P = f x (L/d) x (rho x v^2 / 2)، حيث f هو معامل احتكاك دارسي، وL هو طول الأنبوب (m)، وd هو القطر (m)، وrho هي الكثافة (kg/m^3)، وv هي السرعة (m/s). تكون النتيجة بالباسكال.

فقدان الضاغط يُحسب مباشرة كـ h_f = f x (L/d) x (v^2 / 2g)، مع g = 9.80665 m/s^2 (قيمة التسارع الأرضي القياسية المحددة بواسطة NIST).

يعتمد معامل الاحتكاك على رقم رينولدز والخشونة النسبية؛ بالنسبة للتدفق الصفائحي f = 64/Re، بينما يستخدم التدفق المضطرب معادلة كولبروك أو مخطط مودي.

تشير الموسوعة البريطانية (Encyclopaedia Britannica) إلى أن هذه العلاقة تشكل الأساس لجميع تصاميم الأنابيب الهيدروليكية تقريبًا.

يحدد معامل احتكاك دارسي (f) مقادير المقاومة للتدفق الناتجة عن خشونة جدار الأنبوب والقص اللزج.

بالنسبة للتدفق الصفيحي (Re < 2,300)، يكون المعامل دقيقًا: f = 64/Re. أما بالنسبة للتدفق المضطرب، فهو يعتمد على كل من رقم رينولدز والخشونة النسبية (epsilon/d) ويُحل ضمنيًا باستخدام معادلة كولبروك-وايت، أو يُقرأ من مخطط مودي.

تختلف الخشونة المطلقة بحسب المادة:

• المواد الملساء مثل النحاس المسحوب والبوليمر (PVC) تتميز بخشونة مطلقة منخفضة جدًا (نحو 0.0015 مم و0.0015-0.007 مم)
• الصلب التجاري نحو 0.045 مم
• حديد الزهر نحو 0.26 مم
• الخرسانة 0.3-3 مم

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

لتحديد حجم الأنبوب، أعِد ترتيب المعادلة Q = (pi/4) x d^2 x v لإيجاد القطر: d = sqrt(4Q / (pi x v)). حدد السرعة المستهدفة أولاً، وتكون عادة بين 1-2 م/ث (3-6 قدم/ث) للماء لتحقيق التوازن بين خسائر الاحتكاك والضوضاء والتآكل.

مثال: لنقل 5 لتر/ث (0.005 م^3/ث) بسرعة 1.5 م/ث، يكون d = sqrt((4 x 0.005) / (pi x 1.5)) = sqrt(0.004244) = 0.0651 م، أي نحو 65 مم، لذا تختار الحجم القياسي الأكبر التالي مباشرة.

نظرًا لأن هبوط الضغط يتناسب تقريبًا مع 1/d^5 عند تدفق ثابت، فإن تقليل القطر إلى النصف يضاعف خسائر الاحتكاك بنحو 32 مرة، مما يجعل تحديد الحجم الصحيح أهم قرار تصميمي على الإطلاق.

يُستخدم تحليل تدفق الأنابيب في العديد من التخصصات:

• يقوم السباكون ومهندسو المباني بتحديد أحجام خطوط التغذية والصرف لتلبية طلب التجهيزات دون ضوضاء مفرطة أو طرق مائي، وفقًا للقواعد المستمدة من ASHRAE.
• يحسب مصممو التدفئة والتهوية وتكييف الهواء (HVAC) تدفق المياه المبردة والدورات المائية لتحقيق نقل الحرارة المحدد.
• ينشئ مهندسو البلدية نماذج لشبكات توزيع المياه والمجاري لضمان ضغط كافٍ وسرعات تنظيف ذاتي.
• تستخدمها المصانع الكيميائية والتشغيلية لاختيار المضخات، وتحديد أحجام الصمامات، وأنظمة أمان تصريف الضغط، وغالبًا ما تعتمد على منهجية الورقة الفنية Crane 410.
• يحدد مهندسو الحماية من الحرائق أحجام خطوط الرشاشات الرئيسية للتدفق والضغط المتبقي المطلوبين بموجب معايير NFPA.

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

تغير خصائص المائع رقم رينولدز وهبوط الضغط حتى مع ثبات السرعة. تضطلع الكثافة (rho) بالحد العطالي rho x v^2 / 2 في معادلة دارسي-وايسباخ، بينما تحكم اللزوجة الديناميكية (mu) المقاومة اللزجة ومعامل الاحتكاك.

تكون كثافة الماء عند درجة حرارة 20 مئوية نحو 998 كجم/م^3 ولزوجته نحو 0.001002 باسكال.ثانية؛ وتسخينه إلى 60 مئوية يقلل اللزوجة إلى نحو 0.000466 باسكال.ثانية، مما يرفع رقم رينولدز ويدفع التدفق أكثر نحو النظام المضطرب.

تتميز الموائع اللزجة مثل زيت SAE 30 بلزوجة أعلى بكثير، لذا قد تظل صفيحية عند سرعات يكون فيها الماء مضطربًا.

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

تتكرر عدة أخطاء شائعة عند حساب تدفق الأنابيب:

• من الأخطاء المتكررة استخدام الحجم الاسمي للأنبوب بدلاً من القطر الداخلي الحقيقي؛ إذ يمكن لسمك الجدار أن يقلص مساحة التدفق بنسبة 10 بالمائة أو أكثر، مما يخل بمعادلة Q = A x v.
• خطأ آخر هو خلط أنظمة الوحدات، مثل تعويض المليمترات في صيغة تتطلب الأمتار، مما يشوه المساحة ورقم رينولدز.
• يطبق العديد من المستخدمين صيغة التدفق الصفيحي f = 64/Re على التدفق المضطرب، أو يتجاهلون التوصيلات والأكواع التي تضيف خسائر ثانوية فوق احتكاك الأنابيب المستقيمة.
• كما أن افتراض ثبات خصائص المائع مع تغير درجة الحرارة يؤدي إلى تشويه رقم رينولدز وهبوط الضغط.
• وأخيرًا، ينسى البعض أن هبوط الضغط يتناسب مع مربع السرعة وعكسيًا مع القطر مرفوعًا للقوة الخامسة، لذا فإن التغييرات الصغيرة في الأبعاد لها تأثيرات هائلة.

تأكد دائمًا من المدخلات والوحدات قبل الاعتماد على النتائج.