لهب مخلوط مسبقًا

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

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

انتشار اللهب المخلوط مسبقًا

رقائقي

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

سرعة الاحتراق الصفائحي

بالنسبة للتفاعلات الكيميائية غير القابلة للانعكاس التي تتم في خطوة واحدة، أيνFF+νيايا2Pرoدuجتs{\displaystyle \nu _{F}{\rm {{F}+\nu _{O}{\rm {{O}_{2}\rightarrow {\rm {Products}}}}}}}للهب المستوي الأديباتي تعبير صريح لسرعة الاحتراق مشتق من سلوك طاقة التنشيط التقاربي عندما يكون عدد زيلدوفيتشβ1.{\displaystyle \beta \gg 1.}معدل التفاعلω{\displaystyle \omega }يُفترض أن (عدد مولات الوقود المستهلكة لكل وحدة حجم لكل وحدة زمنية) يأخذ شكل أرهينيوس ،

ω=ب(ρYFدبليوF)م(ρYيا2دبليويا2)نهـ-هـأ/Rتي،\displaystyle \omega =B\left({\frac {\rho Y_{F}}{W_{F}}}\right)^{m}\left({\frac {\rho Y_{O_{2}}}{W_{O_{2}}}}\right)^{n}e^{-E_{a}/RT},}

أينب{\displaystyle B}هو العامل الأسي المسبق ،ρ{\displaystyle \rho }الكثافة ،YF{\displaystyle Y_{F}}يمثل الكسر الكتلي للوقود ،Yيا2{\displaystyle Y_{O_{2}}}يمثل الكسر الكتلي للمؤكسد ،هـأ{\displaystyle E_{a}}هي طاقة التنشيط ،R{\displaystyle R}هو ثابت الغازات العام ،تي{\displaystyle T}هي درجة الحرارة ،دبليوF و دبليويا2{\displaystyle W_{F}\ \&\ W_{O_{2}}}تمثل هذه الأوزان الجزيئية للوقود والمؤكسد على التوالي، وم و ن{\displaystyle m\ \&\ n}تمثل هذه رتب التفاعل. ولنرمز إلى الظروف غير المحترقة البعيدة عن اللهب بالرمز السفلي.u{\displaystyle u}وبالمثل، فإن ظروف الغاز المحترق بواسطةب{\displaystyle b}ثم يمكننا تعريف نسبة التكافؤϕ{\displaystyle \phi }بالنسبة للخليط غير المحترق كما

ϕ=νيا2دبليويا2νFدبليوFYF،uYيا2،u{\displaystyle \phi ={\frac {\nu _{O_{2}}W_{O_{2}}}{\nu _{F}W_{F}}}{\frac {Y_{F,u}}{Y_{O_{2},u}}}}.

ثم سرعة الاحتراق الصفائحي المستوي للخليط الغني بالوقود (ϕ>1{\displaystyle \phi >1}) معطى بواسطة [ 2 ] [ 3 ]

Sل={2بλبρبم+نνFمYيا2،uم+ن-1جي(ن،م،أ)جص،بρu2νيا2دبليويا2م+ن-1βم+ن+1لهـيا2-نلهـF-م}1/2هـ-هـأ/2Rتيب+يا(β-1)،{\displaystyle S_{L}=\left\{{\frac {2B\lambda _{b}\rho _{b}^{m+n}\nu _{F}^{m}Y_{O_{2},u}^{m+n-1}G(n,m,a)}{c_{p,b}\rho _{u}^{2}\nu _{O_{2}}W_{O_{2}}^{m+n-1}\beta ^{m+n+1}\mathrm {Le} _{O_{2}}^{-n}\mathrm {Le} _{F}^{-m}}}\right\}^{1/2}e^{-E_{a}/2RT_{b}}+O(\beta ^{-1}),}

أين

جي(ن،م،أ)=0yن(y+أ)م دy{\displaystyle G(n,m,a)=\int _{0}^{\infty }y^{n}(y+a)^{m}\ dy}

وأ=β(ϕ-1)/لهـF{\displaystyle a=\beta (\phi -1)/\mathrm {Le} _{F}}. هناλ{\displaystyle \lambda }هي الموصلية الحرارية ،جص{\displaystyle c_{p}}هي الحرارة النوعية عند ضغط ثابت ولهـ{\displaystyle \mathrm {Le} }هو رقم لويس . وبالمثل، يمكن كتابة صيغة الهزال.ϕ<1{\displaystyle \phi <1}المخاليط. وقد حصل ت. ميتاني على هذه النتيجة لأول مرة عام 1980. [ 4 ] وقد اشتق فورمان أ. ويليامز وزملاؤه في ثمانينيات القرن الماضي تصحيحًا من الدرجة الثانية لهذه الصيغة مع خصائص نقل أكثر تعقيدًا. [ 5 ] [ 6 ] [ 7 ]

تنشأ اختلافات سرعة الانتشار الموضعية للهب الصفائحي نتيجة لما يُعرف بتمدد اللهب. قد يحدث تمدد اللهب بسبب الإجهاد الناتج عن مجال سرعة التدفق الخارجي أو انحناء اللهب؛ ويُعد الفرق في سرعة الانتشار عن سرعة اللهب الصفائحي المقابلة دالةً لهذه التأثيرات، ويمكن التعبير عنه كما يلي: [ 8 ] [ 9 ]

Sتي=Sل+مجدلتال(Sل-vن)ن-متدلتالتvت{\displaystyle S_{T}=S_{L}+{\mathcal {M}}_{c}\delta _{L}(S_{L}-\mathbf {v} \cdot \mathbf {n} )\nabla \cdot \mathbf {n} -{\mathcal {M}}_{t}\delta _{L}\nabla _{t}\cdot \mathbf {v} _{t}}

أينن{\displaystyle \mathbf {n} }الوحدة عمودية على سطح اللهب (متجهة نحو جانب الغاز المحترق)،v{\displaystyle \mathbf {v} }هل يتم تقييم مجال سرعة التدفق عند سطح اللهب وتvت{\displaystyle \nabla _{t}\cdot \mathbf {v} _{t}}هو تباعد السطح للسرعة المماسيةvت=(أنا-نن)v{\displaystyle \mathbf {v} _{t}=(\mathbf {I} -\mathbf {n} \otimes \mathbf {n} )\mathbf {v} }مع كون التدفق غير قابل للانضغاط خارج اللهب،تvت=-نن:v-(vن)ن{\displaystyle \nabla _{t}\cdot \mathbf {v} _{t}=-\mathbf {n} \otimes \mathbf {n} :\nabla \mathbf {v} -(\mathbf {v} \cdot \mathbf {n} )\nabla \cdot \mathbf {n} } . علاوة على ذلك،مج{\displaystyle {\mathcal {M}}_{c}}ومت{\displaystyle {\mathcal {M}}_{t}}يمثلان رقمي ماركشتاين ، المرتبطين بالانحناء والإجهاد المماسي. [ 10 ]

مضطرب

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

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

يمكن تحليل انتشار هذا اللهب المختلط مسبقًا باستخدام معادلة المجال المسماة معادلة G [ 12 ] [ 13 ] لكمية قياسيةجي{\displaystyle G}مثل:

جيت+vجي=Sتي|جي|{\displaystyle {\frac {\partial G}{\partial t}}+\mathbf {v} \cdot \nabla G=S_{T}|\nabla G|}،

والتي تُعرَّف بحيث تمثل مجموعات المستويات لـ G الواجهات المختلفة داخل اللهب الممزوج مسبقًا والذي ينتشر بسرعة محليةSتي{\displaystyle S_{T}}. ومع ذلك، فإن هذا ليس هو الحال عادةً حيث أن سرعة انتشار الواجهة (بالنسبة للخليط غير المحترق) تختلف من نقطة إلى أخرى بسبب التمدد الديناميكي الهوائي الناتج عن التدرجات في مجال السرعة.

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

تكوين اللهب المخلوط مسبقًا

يؤثر شكل تدفق الغازات المخلوطة مسبقًا على استقرار وخصائص احتراق

لهب بنسن

في لهب بنسن، يُوفَّر معدل تدفق ثابت يُطابق سرعة اللهب لتحقيق استقراره. إذا كان معدل التدفق أقل من سرعة اللهب، سيتحرك اللهب عكس اتجاه التدفق حتى ينفد الوقود أو حتى يصطدم بحاجز . أما إذا كان معدل التدفق مساويًا لسرعة اللهب، فنتوقع جبهة لهب مسطحة ثابتة عمودية على اتجاه التدفق. وإذا كان معدل التدفق أعلى من سرعة اللهب، ستصبح جبهة اللهب مخروطية الشكل بحيث يكون مُركِّب متجه السرعة العمودي على جبهة اللهب مساويًا لسرعة اللهب.

لهيب الركود

هنا، تتدفق الغازات المختلطة مسبقًا بطريقة تشكل منطقة ركود (سرعة صفرية) حيث يمكن تثبيت اللهب.

لهب كروي

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

التطبيقات

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

انظر أيضاً

مراجع

  1. لويس، برنارد؛ إلبه، غونتر فون (2012). احتراق الغازات ولهيبها وانفجاراتها . إلسيفير. ISBN 9780323138024.
  2. ويليامز، ف. أ. (2018). نظرية الاحتراق. مطبعة سي آر سي.
  3. لينان، أ.، وويليامز، ف.أ. (1993). الجوانب الأساسية للاحتراق.
  4. ميتاني، ت. (1980). سرعات انتشار اللهب ثنائي المتفاعل. علوم وتكنولوجيا الاحتراق، 21(3-4)، 175-177.
  5. روغ، ب.، وويليامز، ف.أ. (1985). التحليل التقاربي لانتشار اللهب الصفائحي بمعاملات نقل متغيرة. علم وتكنولوجيا الاحتراق، 42(5-6)، 301-316.
  6. تشيليا، إتش كيه، وويليامز، إف إيه (1987). التحليل التقاربي للهب ثنائي المتفاعل ذي الخصائص المتغيرة ونقل ستيفان-ماكسويل. علم وتكنولوجيا الاحتراق، 51(4-6)، 129-144.
  7. روغ، ب. (1986). حول دقة التنبؤات التقاربية لسرعة اللهب في اللهب ثنائي المتفاعلات. علم وتكنولوجيا الاحتراق، 45(5-6)، 317-329.
  8. كلافين، ب.، وغرانيا-أوتيرو، ج.س. (2011). اللهب المنحني والممتد: رقما ماركشتاين. مجلة ميكانيكا الموائع، 686، 187-217.
  9. كلافين، بول، وجيوف سيربي. موجات الاحتراق والجبهات في التدفقات: اللهب، والصدمات، والانفجارات، وجبهات التآكل، وانفجار النجوم. مطبعة جامعة كامبريدج، 2016.
  10. كلافين، بول؛ غرانيا-أوتيرو، خوسيه سي. (10-11-2011). "اللهب المنحني والممتد: رقما ماركشتاين". مجلة ميكانيكا الموائع . 686 : 187-217 . Bibcode : 2011JFM...686..187C . doi : 10.1017/jfm.2011.318 . ISSN 0022-1120 . 
  11. بيترز، نوربرت (2000). الاحتراق المضطرب . مطبعة جامعة كامبريدج. ISBN 9780511612701. OCLC 56066895 . 
  12. ويليامز، ف. أ. (1985). الاحتراق المضطرب. في رياضيات الاحتراق (ص 97-131). جمعية الرياضيات الصناعية والتطبيقية.
  13. كيرستين، آلان ر. (1988-01-01). "معادلة المجال لانتشار السطح البيني في مجال تدفق متجانس غير مستقر". مجلة Physical Review A. 37 ( 7): 2728–2731 . Bibcode : 1988PhRvA..37.2728K . doi : 10.1103/PhysRevA.37.2728 . PMID 9899999 .