تصميم المبلور: العائد وضبط التشبع الفائق
تصميم المبلور هو قبل كل شيء تمرين في توازن الكتلة: في المثال المحلول أدناه، تستعيد مرور تبريد واحدة من 60 إلى 25 °م نسبة 48.8% من كبريتات النحاس المذابة في صورة بلورات CuSO4·5H2O. العائد يحدده منحنى الذوبانية والتوازن، وحجم البلورة يحدده ضبط التشبع الفائق وزمن الإقامة. يعرض هذا المقال طريقة الحساب ونوافذ التشغيل وأنماط العطل التي تفصل مصنعًا يحقق هدف الاستعادة عن آخر لا يحققه أبدًا.
أهم النقاط
- عائد البلورات يحدده توازن الكتلة ومنحنى الذوبانية، لا اختيار المعدات.
- التشبع الفائق، معبَّرًا عنه بـΔC أو بالتشبع النسبي σ، هو القوة الدافعة للتنوّي والنمو معًا.
- σ = 0.2–2% هو النطاق العملي للحالة المستقرة في معظم أنظمة الأملاح غير العضوية.
- في مبلورات MSMPR يتغير الحجم المتوسط تغيرًا ضعيفًا مع زمن الإقامة، وإدارة النواعم هي الرافعة الأقوى للحجم.
- معظم نقص العائد يعود إلى تحولات الهيدرات وإلى حركية غير مقيسة، وكلاهما يُغلق في الاختبار التجريبي.
الأساسيات: الذوبانية والتشبع الفائق والعائد
منحنى الذوبانية هو الحد الذي يفصل المحلول غير المشبع عن المشبع في التوازن: يعطي لكل درجة حرارة تركيز C* للمذاب (محسوبًا على وحدة الماء الحر، أي الماء الباقي مذيبًا دون الماء المحتجز في الشبكة البلورية) الذي يكون عنده المحلول مشبعًا تمامًا. في كبريتات النحاس خماسية الماء يرتفع C* من نحو 0.207 كغ CuSO4 لامائي لكل كغ ماء حر عند 20 °م إلى نحو 0.400 عند 60 °م. وفي كبريتات الصوديوم ينقطع المنحنى عند 32.4 °م حيث يتحول الصلب المستقر من الميرابيليت (Na2SO4·10H2O) إلى الثينارديت (Na2SO4 اللامائي)، وهو تحول أربك أكثر من مخطط تدفق واحد.
التشبع الفائق حالة ترموديناميكية تنشأ عندما يتجاوز التركيز الفعلي C قيمة التوازن C* عند الحرارة نفسها. وهو القوة الدافعة الوحيدة المتاحة للتنوّي وللنمو. وهو عابر دائمًا بمعنى أن المنظومة تستهلكه: كل غرام من الملح يترسب على سطح بلوري يخفض C عائدًا نحو C*. لذا فالتصميم العملي هو فن توليد التشبع الفائق بمعدل مضبوط وإتلافه عبر مساحة سطح البلورات.
عرض المنطقة الميتاستابلة (MSZW) فترة تركيز أو حرارة تمتد من منحنى الذوبانية إلى النقطة التي يبدأ عندها التنوّي التلقائي (الأولي). وتحت MSZW يمكن للمحلول النظيف أن يبقى مشبعًا فائقًا بلا نهاية دون تنوّي؛ وداخلها يصبح التنوّي محتملًا ثم عنيفًا. ويمكن تشغيل المحاليل المبذورة عند قيم ΔC أقل بكثير من حد MSZW، وهذه بالضبط طريقة تشغيل المعدات الصناعية: المنطقة الميتاستابلة هي غرفة العمل، وحد التنوّي هو جدارها.
وأخيرًا توزيع حجوم البلورات (CSD) إحصاء تجمعات يصف كيف تتوزع الكتلة البلورية على فئات الحجم. العائد يخبرك بكم طنًا تبيع، وCSD يخبرك هل المنتج قابل للطرد المركزي وسهل الانسياب أم طين ناعم يعمي الأقمشة ويحمل المحلول الأم إلى المنتج. ويجري الحساب كله على أساس جاف لامائي — كيلوغرامات مذاب لامائي لكل كيلوغرام ماء حر — لأن البلورات المائية تأخذ ماءً إضافة إلى الملح، والمحاسبة المتراخية لماء التبلور هي أشهر خطأ حسابي في هذا المجال.
حساب العائد: طريقة توازن الكتلة
يُستنتج عائد مرحلة التبلور من توازن للمذاب وتوازن للماء. ولأن الهيدرات تسحب الاثنين، فالتوازنان مترابطان ويُحلان معًا.
الخطوة 1 — تثبيت الأساس. اختر كمية مريحة من محلول التغذية (هنا 1,000 كغ) وجدولها على أساس جاف: كتلة الماء الحر وكتلة المذاب اللامائي الذائب.
الخطوة 2 — أخذ الذوبانية عند درجتي الحرارة. اقرأ C1 (حرارة التغذية) وC2 (الحرارة النهائية) بوحدة كغ مذاب لامائي لكل كغ ماء حر. في المثال: CuSO4 عند 60 °م، C1 ≈ 0.400؛ وعند 25 °م، C2 ≈ 0.230.
الخطوة 3 — حساب نسبة الهيدرات. R = M(هيدرات)/M(لامائي). في CuSO4·5H2O: 249.7 مقابل 159.6 غ/مول، أي R = 1.5646. كل كيلوغرام منتج يحتوي 1/R = 0.6392 كغ CuSO4 لامائي و0.3608 كغ ماء تبلور.
الخطوة 4 — حل التوازن المترابط. لمرور تبريد بلا تبخير (V = 0) يكون عائد بلورات الهيدرات Y = S1·(C1 − C2) / [1/R − C2·(1 − 1/R)]، حيث S1 كتلة الماء الحر الابتدائية. ويصحح المقام للماء الذي يسحبه الهيدرات من الدوران. وإذا كان المرور يبخر أيضًا V كغ من الماء الحر، فاستبدل البسط بـC1·S1 − C2·(S1 − V).
الخطوة 5 — التحويل إلى استعادة. استعادة الملح المذاب = (Y/R)/(C1·S1) كنسبة مئوية.
الخطوة 6 — التحقق من التوازنين. اجمع المذاب اللامائي والماء الحر في المنتج والمحلول الأم؛ ويجب أن يطابق كل منهما التغذية. إذا لم يُطابق توازن الماء فحدّ ماء التبلور خطأ.
المثال المحلول. الأساس: 1,000 كغ محلول منقّى عند 60 °م يحتوي 28.57% وزنية CuSO4 لامائي (C1 = 0.400؛ S1 = 714.3 كغ). التبريد إلى 25 °م عند C2 = 0.230، بلا تبخير ولا حمل، والمنتج CuSO4·5H2O نقي. البسط: 714.3 × (0.400 − 0.230) = 121.4. المقام: 0.6392 − 0.230 × 0.3608 = 0.5562. Y = 121.4 / 0.5562 = 218.3 كغ بلورات CuSO4·5H2O. الاستعادة = (218.3 × 0.6392) / 285.7 = 139.5 / 285.7 = 48.8% من كبريتات النحاس المذابة.
| التيار | ماء حر (كغ) | CuSO4 لامائي (كغ) | الإجمالي (كغ) |
|---|---|---|---|
| محلول التغذية، 60 °م | 714.3 | 285.7 | 1,000.0 |
| بلورات المنتج | 78.8 (ماء هيدرات) | 139.5 | 218.3 |
| المحلول الأم، 25 °م | 635.5 | 146.2 | 781.7 |
يتطابق التوازنان (714.3 = 78.8 + 635.5؛ 285.7 = 139.5 + 146.2؛ 218.3 + 781.7 = 1,000.0)، ويعيد تركيز المحلول الأم إنتاج C2: 146.2 / 635.5 = 0.230. والحساسية: دفع المحلول نفسه إلى 20 °م (C2 ≈ 0.207) يعطي Y = 714.3 × 0.193 / 0.5645 = 244.3 كغ، أي استعادة 54.7% — نحو ست نقاط مئوية إضافية مقابل 5 K تبريد إضافية. والأرقام التي يستند إليها التصميم جيدة بقدر جودة منحنى الذوبانية خلفها، وفي محاليل المصنع — مع الحديد والحمض والمغنيسيوم — ينزاح المنحنى؛ ولهذا نتعامل مع ذوبانية الأدبيات كتقدير ابتدائي لا كمدخل تصميمي.
كيف يُولَّد التشبع الفائق ويُقاس
بما أن العائد يأتي من التوازن، فإن الشيء الوحيد المتبقي لضبطه هو معدل ومسار توليد التشبع الفائق واستهلاكه. التشبع المطلق ΔC = C − C*، والنسبي σ = ΔC/C*. مثال مصغّر: عند 50 °م تكون C* للنحاس ≈ 0.336 كغ/كغ؛ ومحلول عند 0.340 له ΔC = 0.004 كغ/كغ وσ = 0.004/0.336 = 1.2%. والممارسة الصناعية للأملاح غير العضوية تبقي σ في نطاق 0.2–2% في الحالة المستقرة: عالٍ بما يكفي لمعدلات نمو مقبولة ومنخفض بما يكفي لتفادي جدار التنوّي.
| الطريقة | كيف يُولَّد ΔC | σ النمطي في الحالة المستقرة | الخدمة الهيدروميتالورجية النمطية |
|---|---|---|---|
| تبريد سطحي | تنخفض الذوبانية مع الحرارة عبر مبادل حراري | 0.5–2% | CuSO4·5H2O وMgSO4·7H2O والشبّة |
| تبخير | يُزال المذيب عند حرارة شبه ثابتة | 0.2–1% | Na2SO4 وMgCl2 وأملاح نيكل ذات منحنى مسطح |
| تبريد تحت تفريغ (أدياباتي) | يزيل التبخر الوميضي الحرارة المحسوسة والمذيب | 0.5–2% | الميرابيليت والبوراكس |
| التمليح بالإزاحة | مذيب مضاد قابل للامتزاج يخفض C* | 1–5%، قصير الأمد | أنظمة كحول/محلول ملحي في تكرير Ni–Co |
| تفاعل | يتولد المذاب في الموقع أسرع مما يمكن أن يترسب | 1–5%، عابر | الجبس وMg(OH)2 وترسيب النيكل |
وللمسار التبخيري يكون المعدات المختارة عادةً المبخر-المبلور بالدوران القسري، حيث يحفظ محيط الدوران سطح نقل الحرارة فقيرًا بالتشبع الفائق ويترك ΔC ينمو في منطقة الغليان. وللمنتجات الحساسة للحجم في خدمات التبريد يحفظ مبلور Oslo طبقة مميعة من البلورات الخشنة مغمورة في محلول قليل التشبع الفائق، مقايضًا السعة بجودة التوزيع.
وتستحق MSZW فقرة إضافية لأنها تترجم مباشرة إلى سياسة ضبط. القيم المقيسة للأنظمة الكبريتية الشائعة تمتد نحو 5–25 K من التبريد عند معدلات 5–20 K/ساعة بلا بذور (أرقام من رتبة الأدبيات وتعتمد على النظام والشوائب). والتبذير يمكن أن يقلّص المدى القابل للاستخدام إلى النصف. والحكم الذي لا يأتي إلا من الممارسة: نبذّر عند نحو ثلث MSZW المقيسة لا عند نقطة التشبع الترموديناميكية، ولا نسمح أبدًا للعابرين الابتدائيين بأن تحدد نغمة تنوّي التشغيلة. والوعاء الذي ينوّي بشكل سيئ في الثلاثين دقيقة الأولى يقضي بقية الحملة في إنتاج النواعم.
توازن التجمعات في MSMPR: التنوّي مقابل النمو
يوصف معظم المبلورات الصناعية المستمرة وصفًا جيدًا بمثالية MSMPR. وفي افتراضاتها يعطي توازن التجمعات الكثافة العددية للبلورات n(L) = n⁰·exp(−L/Gτ)، حيث n⁰ كثافة التنوّي وG سرعة النمو الخطية وτ زمن إقامة المنتج.
يهيمن قانونان حركيان على السلوك. النمو يتبع G = kg·σ^g حيث g عادةً 1–2؛ والتنوّي (ثانوي في الغالب، تلامسي على بلورات قائمة) يتبع B⁰ = kb·σ^b حيث b عادةً 2–4، أي أن التنوّي دالة أكثر حدة في التشبع الفائق. وهذه اللامساواة في الأسّ هي الحقيقة المركزية في ضبط التشبع الفائق: عندما يرتفع σ تظهر بلورات جديدة أسرع مما تنمو القائمة، فيتوزع إنتاج الملح نفسه على بلورات أكثر وأصغر.
والنتيجة الكلاسيكية للحجم: الحجم المتوسط المرجح بالكتلة L43 = 3.67·G·τ. مثال مصغّر: عند G = 1.0 × 10⁻⁸ م/ث وτ = 6 ساعات، L43 = 3.67 × 1.0 × 10⁻⁸ × 21,600 = 0.79 مم. وهل مضاعفة زمن الإقامة إلى 12 ساعة تضاعفه؟ لا، وتوازن التجمعات يفسر السبب: عند كثافة معلّق ثابتة يشترط توازن الكتلة B⁰·G³·τ⁴ = ثابت؛ وبالتعويض بالحركية نحصل على σ ∝ τ^(−4/(b+3g))، ومع b = 2g يكون L43 ∝ τ^0.2. فمضاعفة τ ترفع L43 نحو 2^0.2 = 1.15 مرة فقط، إلى ≈ 0.91 مم. وهذه القاعدة الثانية في التصميم من أجل الحجم: الأوعية الأكبر تشتري قليلًا بشكل مفاجئ، والروافع القوية هي خفض التنوّي وإزالة النواعم.
| المعامل | الرمز | النطاق النمطي (أملاح غير عضوية) | الأثر عند الزيادة |
|---|---|---|---|
| زمن إقامة المنتج | τ | 2–12 ساعة | L43 أكبر لكن بضعف (∝ τ^0.2 عند b = 2g)؛ رأس مال أكبر للطن |
| سرعة النمو الخطية | G | 10⁻⁹–10⁻⁷ م/ث | إنتاجية أعلى عند σ نفسه |
| سرعة التنوّي | B⁰ | 10⁵–10⁸ عدد/(م³·ث) | نواعم أكثر وL43 أصغر |
| جزيء النواعم المُزال | Rf | 5–20% من الماغما الدائرة | النواعم المذابة تعود نموًا؛ توزيع أخشن وأضيق |
| كثافة المعلّق | MT | 50–400 كغ/م³ | مساحة أكبر لكل وحدة ΔC؛ تلطّف σ وقد ترفع التنوّي التلامسي |
واستغلال هذه الروافع في مخطط مستمر أو دفعي يعتمد على المقياس وتغير التغذية ومواصفة الحجم؛ وقيمة إطار MSMPR العملية تناسبية: يقول أي مقبض يحرك التوزيع وبأي قوة، قبل مطابقة أي بيانات تجريبية.
أكثر أخطاء التصميم شيوعًا
تتكرر في مصانع الأملاح، وكل بند مذكور كعرض — سبب — علاج.
- عائد محسوب من التوازن الحراري وحده. العرض: تخفق استعادة المصنع في بلوغ رقم التصميم بـ5–15 نقطة مئوية. السبب: استُخدمت ذوبانية الماء النقي أو أُهمل ماء التبلور وتحولات الهيدرات. العلاج: إعادة توازن منحنى الذوبانية على تركيب المحلول الفعلي، على أساس جاف، مع محاسبة صريحة لماء التبلور.
- إهمال تحولات الهيدرات. العرض: مصنع كبريتات الصوديوم ينتج نواعم طينية تحت 32.4 °م أو يتبدد العائد في موسم دافئ. السبب: تغيّر الميرابيليت/الثينارديت يجعل الذوبانية وR والطلب على الماء تتغير تغيرًا متقطعًا. العلاج: اختيار نطاق الحرارة عمدًا — تبلور الميرابيليت باردًا لعائد أعلى (R أعلى) أو الثينارديت ساخنًا لمنتج سهل الانسياب — وتصميم المسار الحراري حول الطور المختار.
- ضبط منخفض جدًا لتصريف النواعم. العرض: حجم المنتج المتوسط ينزلق أسبوعًا بعد أسبوع وتنسدّ أقمشة الطاردة. السبب: حلقة إزالة النواعم صغيرة أو لم تُضبط قيمتها أبدًا مقابل معدل التنوّي الفعلي. العلاج: قياس التوزيع ثم ضبط التصريف بحيث تزيل الحلقة مساحة كافية وتحوّل النمو إلى البلورات الباقية.
- قشرة على أسطح نقل الحرارة بسبب فرض ΔC. العرض: تضمحل الحمولة الحرارية خلال أيام وتقصر دورات التنظيف. السبب: يُولد التشبع الفائق عند الجدار لا في الكتلة. العلاج: خفض القوة الدافعة، ورفع سرعة الدوران، ونقل التبخير بعيدًا عن المبادل إلى منطقة الغليان.
- انحراف الحجم بعد تغيّر التغذية. العرض: توزيع ناعم بعد ترقية في المرحلة الأعلى ترفع المحتوى. السبب: إنتاجية أعلى عند هندسة ثابتة تتطلب ΔC أكبر، ولأن B⁰ يستجيب أحدّ من G فإن التنوّي يفوز. العلاج: إعادة فحص توازن التجمعات عند المحتوى الجديد وإضافة سعة تبذير أو تدمير نواعم.
- غياب قياس MSZW قبل التصميم التفصيلي. العرض: الوعاء ينوّي بعنف عند بدء التشغيل وأسابيع من منتج خارج المواصفة. السبب: افتُرضت بيانات أدبيات، لكن شوائب محلول المصنع ضيّقت MSZW. العلاج: قياس MSZW على المحلول الحقيقي، مبذورًا وغير مبذور، عند معدلات تبريد ممثلة.
وحكم من الميدان مرة أخرى: معظم هذه الإخفاقات ليست أخطاء حساب بل أخطاء محاسبة: ماء الهيدرات، وتيارات التدوير، والذوبانية المنزاحة بالشوائب. ولذا يجب أن تبدأ كل مراجعة تصميم وتنتهي عند توازن الكتلة، الذي لا يغفر شيئًا ولا ينسى شيئًا.
سياق التكلفة والطاقة: التبخير مقابل التبريد
تفصل الطاقة بين مساري التوليد أكثر مما يفصل رأس المال. يجب أن يوفر التبخير الحرارة الكامنة للماء: نحو 2,260 كيلوجول/كغ عند 100 °م ونحو 2,407 عند 40 °م تحت تفريغ. ويستهلك المبخر أحادي التأثير نحو 1.1–1.3 كغ بخار لكل كغ ماء مبخر؛ وخط ثلاثي التأثير 0.35–0.45 كغ/كغ؛ ووحدة إعادة الضغط الميكانيكي (MVR) 20–45 كيلوواط·ساعة لكل طن ماء مبخر، بافتراض ΔT فعال 4–8 K وارتفاع نقطة الغليان (BPE) 3–8 K لمحاليل كبريتية مشبعة وكفاءة ضاغط 75–80%. وللمقارنة، فإن 30 كيلوواط·ساعة/طن تمثل 4–5% فقط من المحتوى الحراري لتلك المهمة بالبخار، ولهذا تهيمن MVR على تصاميم التبخير الجديدة حيث تتوفر الكهرباء.
والتبلور بالتبريد أرخص ترموديناميكيًا بكثير. الحمولة حرارة محسوسة زائد حرارة تبلور؛ وفي CuSO4·5H2O تبلغ الأخيرة ≈ 78 كيلوجول/مول، أي ≈ 0.31 ميجاجول لكل كغ هيدرات. وفي المثال: الحمولة المحسوسة ≈ 1,000 كغ × 3.5 كيلوجول/(كغ·K) × 35 K = 122.5 ميجاجول، وحرارة التبلور ≈ 218.3 × 0.31 = 67.7 ميجاجول، فيرى المبرد ≈ 190 ميجاجول لكل طن تغذية — نحو 53 كيلوواط·ساعة — وهي مقاربة للطن لحمولة MVR لكن مقدمة كبرودة منخفضة الحرارة لا كشغل محور. واقتصاد الاختيار يتوقف على ميل الذوبانية: المنحنيات الحادة (كبريتات النحاس) تكافئ التبريد، والمسطحة تفرض التبخير.
وتحذير بشأن BPE: لأنها تقلص ΔT الفعال المتاح للمبادل وللضاغط معًا، فإن التقليل من تقديرها في محلول مركّز سبب كلاسيكي لسعة تبخير أصغر من اللازم. ولا شيء من هذه الأرقام عرض سعر؛ إنها نطاقات سياق تصميمي.
ما العوائد المتوقعة للنحاس والنيكل وكبريتات الصوديوم
ثلاث حالات استرشادية بالطريقة نفسها، وكلها موسومة بوضوح بأنها استرشادية (ذوبانيات من رتبة الأدبيات، ومحاليل نقية، وبلا تدوير).
كبريتات النحاس (تبريد): كما حُلّ أعلاه، 1,000 كغ محلول 28.57% وزنية مبردة من 60 إلى 25 °م تعطي 218.3 كغ CuSO4·5H2O واستعادة 48.8%؛ وتعميق التبريد إلى 20 °م يرفعها إلى 54.7%. المنحنى الحاد يجعل كبريتات النحاس من الأملاح الهيدروميتالورجية القليلة التي يعيد فيها التبريد وحده نصف قيمة المعدن المذاب في مرور واحد.
كبريتات النيكل (تبخير): أساس 1,000 كغ محلول منقّى عند 30% وزنية NiSO4 (أساس جاف)، تبخير 250 كغ ماء حر، وانتهاء مشبع عند 30 °م بـC2 استرشادي 0.50 كغ/كغ. وبـR = 262.85/154.75 = 1.699 للهيدرات السداسي: Y = (0.4286 × 700 − 0.50 × 450) / (0.5887 − 0.50 × 0.4112) = 75 / 0.3832 = 195.7 كغ NiSO4·6H2O، أي 38.4% من النيكل المذاب. والمنحنى المسطح نسبيًا هو سبب كون مصانع كبريتات النيكل تبخيرية.
كبريتات الصوديوم (تبريد تحت تفريغ): أساس 1,000 كغ محلول 25% وزنية Na2SO4، مبرد من 40 إلى 20 °م بلا تبخير صافٍ. وبـR = 322.20/142.04 = 2.269 للميرابيليت: Y = 750 × (0.3333 − 0.195) / (0.4409 − 0.195 × 0.5591) = 103.7 / 0.3318 = 312.7 كغ Na2SO4·10H2O، باستعادة 55.1% — فنسبة الهيدرات العالية (أكثر من نصف البلورة ماء) تدفع العائد لأعلى. وأي تبخر وميضي في المبرد التفريغي يضيف استعادة فوق هذا الحد الأدنى.
ولاحظ التباين: 38–55% استعادة في مرور واحد بالطريقة نفسها، ولهذا يكون تصميم المبلور خاصًا بالملح دائمًا.
الأسئلة الشائعة
ما الفرق بين التشبع الفائق المطلق والنسبي؟
ΔC هو القوة الدافعة المطلقة: الفرق بين التركيز الفعلي والذوبانية المتوازنة C* عند الحرارة نفسها، بوحدة كغ/كغ ماء حر. والتشبع النسبي σ = ΔC/C* يعبّر عن هذا الفرق كجزء من الذوبانية، وهو ما ترتبط به فعليًا سرعتا التنوّي والنمو. استخدم ΔC للمحاسبة وتحديد أحجام المعدات، وσ لمقارنة الأنظمة ووضع حدود الضبط.
هل يزيد التبذير عائد البلورات؟
لا. العائد يحدده توازن الكتلة ومنحنى الذوبانية، فالتبذير لا يغيّر كمية الملح التي تخرج من المحلول في التوازن. ما يغيّره هو كيف يظهر ذلك العائد: يستهلك التشبع الفائق المتاح بالنمو على أسطح قائمة بدل التنوّي التلقائي، فينتج بلورات أقل عددًا وأكبر حجمًا وبتوزيع أضيق. وتُعدّ كمية البذور وحجمها وتوقيتها من أكثر المتغيرات التشغيلية حساسية.