رقم غراشوف

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

تعريف

انتقال الحرارة

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

رقم غراشوف هو:

جيرل=زβ(تيs-تي)ل3ν2{\displaystyle \mathrm {Gr} _{L}={\frac {g\beta (T_{s}-T_{\infty })L^{3}}{\nu ^{2}}}\,}للألواح المسطحة الرأسية
جيرد=زβ(تيs-تي)د3ν2{\displaystyle \mathrm {Gr} _{D}={\frac {g\beta (T_{s}-T_{\infty })D^{3}}{\nu ^{2}}}\,}للأنابيب والأجسام الصلبة

أين:

تشير الرموز السفلية L و D إلى أساس مقياس الطول لرقم غراشوف.

يحدث الانتقال إلى التدفق المضطرب في النطاق 10⁸ < Gr L < 10⁹ للحمل الحراري الطبيعي من الألواح المسطحة الرأسية. عند أرقام غراشوف الأعلى، تكون الطبقة الحدية مضطربة؛ وعند أرقام غراشوف الأدنى ، تكون الطبقة الحدية انسيابية، أي في النطاق 10³ < Gr L < 10⁶ .

انتقال الكتلة

يوجد شكل مماثل لرقم غراشوف يُستخدم في مسائل انتقال الكتلة بالحمل الحراري الطبيعي . في حالة انتقال الكتلة، ينشأ الحمل الحراري الطبيعي عن تدرجات التركيز وليس عن تدرجات درجة الحرارة. [ 2 ]

جيرج=زβ*(جأ،s-جأ،أ)ل3ν2{\displaystyle \mathrm {Gr} _{c}={\frac {g\beta ^{*}(C_{a,s}-C_{a,a})L^{3}}{\nu ^{2}}}}

أين

β*=-1ρ(ρجأ)تي،ص{\displaystyle \beta ^{*}=-{\frac {1}{\rho }}\left({\frac {\partial \rho }{\partial C_{a}}}\right)_{T,p}}

و:

  • g هو التسارع الناتج عن الجاذبية الأرضية
  • يمثل C a,s تركيز النوع a على السطح
  • C a,a هو تركيز النوع a في الوسط المحيط
  • L هو الطول المميز
  • ν هي اللزوجة الحركية
  • ρ هي كثافة السائل
  • C a هو تركيز النوع a
  • T هي درجة الحرارة (ثابتة)
  • p هو الضغط (ثابت).

العلاقة بالأعداد الأخرى عديمة الأبعاد

يُعدّ عدد رايلي ، الموضح أدناه، عددًا لا بُعديًا يُميّز مسائل الحمل الحراري في انتقال الحرارة. توجد قيمة حرجة لعدد رايلي ، يتجاوزها يحدث تحرك المائع. [ 3 ]

Rأx=جيرxPر{\displaystyle \mathrm {Ra} _{x}=\mathrm {Gr} _{x}\mathrm {Pr} }

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

Rأنا=جيرRهـ21تجاهل الحمل الحراري القسري{\displaystyle \mathrm {Ri} ={\frac {\mathrm {Gr} }{\mathrm {Re} ^{2}}}\gg 1\يتضمن {\text{تجاهل الحمل القسري}}}
Rأنا=جيرRهـ21الحمل الحراري القسري والحمل الحراري الطبيعي معًا{\displaystyle \mathrm {Ri} ={\frac {\mathrm {Gr} }{\mathrm {Re} ^{2}}}\approx 1\implies {\text{الحمل الحراري القسري والطبيعي المدمج}}}
Rأنا=جيرRهـ21تجاهل الحمل الحراري الطبيعي{\displaystyle \mathrm {Ri} ={\frac {\mathrm {Gr} }{\mathrm {Re} ^{2}}}\ll 1\يتضمن {\text{تجاهل الحمل الحراري}}}

الاشتقاق

تتمثل الخطوة الأولى لاستنتاج رقم غراشوف في معالجة معامل تمدد الحجم،β{\displaystyle \mathrm {\beta } }على النحو التالي.

β=1v(vتي)ص=-1ρ(ρتي)ص\beta = \frac{1}{v} \left( \frac{\partial v}{\partial T} \right)_p = \frac{-1}{\rho} \left( \frac{\partial \rho}{\partial T} \right)_p

الv{\displaystyle v}إن المعادلة أعلاه، التي تمثل الحجم النوعي ، ليست هي نفسهاv{\displaystyle v}في الأقسام اللاحقة من هذا الاشتقاق، والتي ستمثل سرعة. هذه العلاقة الجزئية لمعامل تمدد الحجم،β{\displaystyle \mathrm {\beta } }فيما يتعلق بكثافة السائل،ρ{\displaystyle \mathrm {\rho } }يمكن إعادة كتابة المعادلة، عند ثبات الضغط، على النحو التالي:

ρ=ρ0(1-βΔتي){\displaystyle \rho =\rho _{0}(1-\beta \Delta T)}

أين:

  • ρ0{\displaystyle \rho _{0}}هي كثافة السائل الكلية
  • ρ{\displaystyle \rho }هي كثافة الطبقة الحدية
  • Δتي=(تي-تي0){\displaystyle \Delta T=(T-T_{0})}، الفرق في درجة الحرارة بين الطبقة الحدية والسائل الرئيسي.

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

معادلة الطاقة

تتناول هذه المناقشة معادلة الطاقة في سياق التدفق المتناظر دورانيًا. وسيأخذ هذا التحليل في الاعتبار تأثير تسارع الجاذبية على التدفق وانتقال الحرارة. تنطبق المعادلات الرياضية التالية على كلٍ من التدفق المتناظر دورانيًا والتدفق المستوي ثنائي الأبعاد.

s(ρuر0ن)+y(ρvر0ن)=0{\displaystyle {\frac {\partial }{\partial s}}(\rho ur_{0}^{n})+{\frac {\partial }{\partial y}}(\rho vr_{0}^{n})=0}

أين:

  • s{\displaystyle s}هو اتجاه الدوران، أي الاتجاه الموازي للسطح
  • u{\displaystyle u}هي السرعة المماسية، أي السرعة الموازية للسطح
  • y{\displaystyle y}هو الاتجاه المستوي، أي الاتجاه العمودي على السطح
  • v{\displaystyle v}هي السرعة العمودية، أي السرعة العمودية على السطح
  • ر0{\displaystyle r_{0}}هو نصف القطر.

في هذه المعادلة، يُستخدم الرمز العلوي n للتمييز بين التدفق المتناظر دورانيًا والتدفق المستوي. وتتحقق الخصائص التالية لهذه المعادلة.

  • ن{\displaystyle n}= 1: تدفق متناظر دورانيًا
  • ن{\displaystyle n}= 0: تدفق ثنائي الأبعاد مستوٍ
  • ز{\displaystyle g}التسارع الجاذبي

تتوسع هذه المعادلة إلى ما يلي بإضافة خصائص السوائل الفيزيائية:

ρ(uus+vuy)=y(μuy)-دصدs+ρز.{\displaystyle \rho \left(u{\frac {\partial u}{\partial s}}+v{\frac {\partial u}{\partial y}}\right)={\frac {\partial }{\partial y}}\left(\mu {\frac {\partial u}{\partial y}}\right)-{\frac {dp}{ds}}+\rho g.}

ومن هنا يمكننا تبسيط معادلة الزخم بشكل أكبر عن طريق ضبط سرعة السائل الكلية على 0 (u=0{\displaystyle u=0}).

دصدs=ρ0ز{\displaystyle {\frac {dp}{ds}}=\rho _{0}g}

تُظهر هذه العلاقة أن تدرج الضغط هو ببساطة ناتج ضرب كثافة المائع الكلية في تسارع الجاذبية الأرضية. والخطوة التالية هي إدخال تدرج الضغط في معادلة الزخم.

(uus+vuy)=ν(2uy2)+زρ-ρ0ρ=ν(2uy2)-ρ0ρزβ(تي-تي0){\displaystyle \left(u{\frac {\partial u}{\partial s}}+v{\frac {\partial u}{\partial y}}\right)=\nu \left({\frac {\partial ^{2}u}{\partial y^{2}}}\right)+g{\frac {\rho -\rho _{0}}{\rho }}=\nu \left({\frac {\partial ^{2}u}{\partial y^{2}}}\right)-{\frac {\rho _{0}}{\rho }}g\beta (T-T_{0})}

حيث العلاقة بين معامل التمدد الحجمي والكثافةρ-ρ0=-ρ0β(تي-تي0){\displaystyle \rho -\rho _{0}=-\rho _{0}\beta (T-T_{0})}تم العثور على ما سبق وعلاقة اللزوجة الحركيةν=μρ{\displaystyle \nu ={\frac {\mu }{\rho }}}تم استبدالها في معادلة الزخم.u(us)+v(vy)=ν(2uy2)-ρ0ρزβ(تي-تي0){\displaystyle u\left({\frac {\partial u}{\partial s}}\right)+v\left({\frac {\partial v}{\partial y}}\right)=\nu \left({\frac {\partial ^{2}u}{\partial y^{2}}}\right)-{\frac {\rho _{0}}{\rho }}g\beta (T-T_{0})}

لإيجاد عدد غراشوف من هذه النقطة، يجب تحويل المعادلة السابقة إلى صيغة لا بُعدية. وهذا يعني أن كل متغير في المعادلة يجب أن يكون بلا بُعد، وأن يكون بدلاً من ذلك نسبة مميزة لهندسة المسألة وإعدادها. ويتم ذلك بقسمة كل متغير على الكميات الثابتة المناظرة له. تُقسم الأطوال على طول مميز.لج{\displaystyle L_{c}}يتم تقسيم السرعات على سرعات مرجعية مناسبة.V{\displaystyle V}، وهو ما يعطي، مع الأخذ في الاعتبار رقم رينولدزV=Rهـلνلج{\displaystyle V={\frac {\mathrm {Re} _{L}\nu }{L_{c}}}}يتم تقسيم درجات الحرارة على فرق درجة الحرارة المناسب.(تيs-تي0){\displaystyle (T_{s}-T_{0})}تبدو هذه المعاملات عديمة الأبعاد كما يلي:

  • s*=sلج{\displaystyle s^{*}={\frac {s}{L_{c}}}}،
  • y*=yلج{\displaystyle y^{*}={\frac {y}{L_{c}}}}،
  • u*=uV{\displaystyle u^{*}={\frac {u}{V}}}،
  • v*=vV{\displaystyle v^{*}={\frac {v}{V}}}، و
  • تي*=(تي-تي0)(تيs-تي0){\displaystyle T^{*}={\frac {(T-T_{0})}{(T_{s}-T_{0})}}}.

تمثل العلامات النجمية معلمات لا بُعدية. وبدمج هذه المعادلات اللا بُعدية مع معادلات الزخم، نحصل على المعادلة المبسطة التالية.

u*u*s*+v*u*y*=-[ρ0زβ(تيs-تي0)لج3ρν2Rهـل2]تي*+1Rهـل2u*y*2{\displaystyle u^{*}{\frac {\partial u^{*}}{\partial s^{*}}}+v^{*}{\frac {\partial u^{*}}{\partial y^{*}}}=-\left[{\frac {\rho _{0}g\beta (T_{s}-T_{0})L_{c}^{3}}{\rho \nu ^{2}\mathrm {Re} _{L}^{2}}}\right]T^{*}+{\frac {1}{\mathrm {Re} _{L}}}{\frac {\partial ^{2}u^{*}}{\partial {y^{*}}^{2}}}}
=-(ρ0ρ)[زβ(تيs-تي0)لج3ν2]تي*Rهـل2+1Rهـل2u*y*2{\displaystyle =-\left({\frac {\rho _{0}}{\rho }}\right)\left[{\frac {g\beta (T_{s}-T_{0})L_{c}^{3}}{\nu ^{2}}}\right]{\frac {T^{*}}{\mathrm {Re} _{L}^{2}}}+{\frac {1}{\mathrm {Re} _{L}}}{\frac {\partial ^{2}u^{*}}{\partial {y^{*}}^{2}}}}

أين:

تيs{\displaystyle T_{s}}هي درجة حرارة السطح
تي0{\displaystyle T_{0}}درجة حرارة السائل الكلية
لج{\displaystyle L_{c}}هو الطول المميز.

يُعرف المعامل عديم الأبعاد الموجود بين قوسين في المعادلة السابقة باسم رقم غراشوف:

جير=زβ(تيs-تي0)لج3ν2.{\displaystyle \mathrm {Gr} ={\frac {g\beta (T_{s}-T_{0})L_{c}^{3}}{\nu ^{2}}}.}

نظرية باكنغهام باي

هناك شكل آخر من أشكال التحليل البُعدي الذي ينتج عنه عدد غراشوف يُعرف باسم نظرية باكنغهام باي . تأخذ هذه الطريقة في الاعتبار قوة الطفو لكل وحدة حجم.Fب{\displaystyle F_{b}}بسبب اختلاف الكثافة في الطبقة الحدية والسائل الرئيسي.

Fب=(ρ-ρ0)ز{\displaystyle F_{b}=(\rho -\rho _{0})g}

يمكن التلاعب بهذه المعادلة للحصول على،

Fب=-βزρ0Δتي.{\displaystyle F_{b}=-\beta g\rho _{0}\Delta T.}

فيما يلي قائمة المتغيرات المستخدمة في طريقة باكنغهام π، إلى جانب رموزها وأبعادها.

عاملرمزأبعاد
طول كبيرل{\displaystyle L}ل{\displaystyle \mathrm {L} }
لزوجة السائلμ{\displaystyle \mu }ملت{\displaystyle \mathrm {\frac {M}{Lt}} }
السعة الحرارية للسوائلجص{\displaystyle c_{p}}سؤالمتي{\displaystyle \mathrm {\frac {Q}{MT}} }
الموصلية الحرارية للسوائلك{\displaystyle k}سؤاللتتي{\displaystyle \mathrm {\frac {Q}{LtT}} }
معامل تمدد الحجمβ{\displaystyle \beta }1تي{\displaystyle \mathrm {\frac {1}{T}} }
التسارع الناتج عن الجاذبيةز{\displaystyle g}لت2{\displaystyle \mathrm {\frac {L}{t^{2}}} }
فرق درجة الحرارةΔتي{\displaystyle \Delta T}تي{\displaystyle \mathrm {T} }
معامل انتقال الحرارةح{\displaystyle h}سؤالل2تتي{\displaystyle \mathrm {\frac {Q}{L^{2}tT}} }

بالرجوع إلى نظرية باكنغهام باي، يوجد 9 - 5 = 4 مجموعات لا بُعدية. اختر L ،μ،{\displaystyle \mu ,}ك ، ج وβ{\displaystyle \beta }باعتبارها المتغيرات المرجعية. وبالتالي فإنπ{\displaystyle \pi }المجموعات هي كالتالي:

π1=لأμبكجβدزهـجص{\displaystyle \pi _{1}=L^{a}\mu ^{b}k^{c}\beta ^{d}g^{e}c_{p}}،
π2=لوμزكحβأنازجρ{\displaystyle \pi _{2}=L^{f}\mu ^{g}k^{h}\beta ^{i}g^{j}\rho }،
π3=لكμلكمβنزoΔتي{\displaystyle \pi _{3}=L^{k}\mu ^{l}k^{m}\beta ^{n}g^{o}\Delta T}،
π4=لqμركsβتزuح{\displaystyle \pi _{4}=L^{q}\mu ^{r}k^{s}\beta ^{t}g^{u}h}.

حل هذه المشكلاتπ{\displaystyle \pi }المجموعات تعطي:

π1=μ(جص)ك=Pر{\displaystyle \pi _{1}={\frac {\mu (c_{p})}{k}}=\mathrm {Pr} }،
π2=ل3زρ2μ2{\displaystyle \pi _{2}={\frac {l^{3}g\rho ^{2}}{\mu ^{2}}}}،
π3=βΔتي{\displaystyle \pi _{3}=\beta \Delta T}،
π4=حلك=شمالu{\displaystyle \pi _{4}={\frac {hL}{k}}=\mathrm {Nu} }

من المجموعتينπ2{\displaystyle \pi _{2}}وπ3،{\displaystyle \pi _{3},}يشكل الناتج رقم غراشوف:

π2π3=βزρ2Δتيل3μ2=جير.{\displaystyle \pi _{2}\pi _{3}={\frac {\beta g\rho ^{2}\Delta TL^{3}}{\mu ^{2}}}=\mathrm {Gr} .}

أخذν=μρ{\displaystyle \nu ={\frac {\mu }{\rho }}}وΔتي=(تيs-تي0){\displaystyle \Delta T=(T_{s}-T_{0})}يمكن صياغة المعادلة السابقة على أنها نفس النتيجة من اشتقاق رقم غراشوف من معادلة الطاقة.

جير=βزΔتيل3ν2{\displaystyle \mathrm {Gr} ={\frac {\beta g\Delta TL^{3}}{\nu ^{2}}}}

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

التفكير الفيزيائي

من الممكن أيضاً استنتاج رقم غراشوف من خلال التعريف الفيزيائي للرقم على النحو التالي:

جير=بuoyأنجy FoرجهـFرأناجتأناoن Foرجهـ=مزτأ=ل3ρβ(Δتي)زμ(V/ل)ل2=ل2β(Δتي)زνV{\displaystyle \mathrm {Gr} ={\frac {\mathrm {Buoyancy~Force} }{\mathrm {Friction~Force} }}={\frac {mg}{\tau A}}={\frac {L^{3}\rho \beta (\Delta T)g}{\mu (V/L)L^{2}}}={\frac {L^{2}\beta (\Delta T)g}{\nu V}}}

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

μ=ρVل{\displaystyle \mathrm {\mu } =\rho VL}الكتابة فوق المقياس بالجرام تعطي؛

جير=ل3β(Δتي)زν2{\displaystyle \mathrm {Gr} ={\frac {L^{3}\beta (\Delta T)g}{\nu ^{2}}}}يُساعد التفكير الفيزيائي على فهم معنى الرقم. من جهة أخرى، يمكن استخدام تعريف السرعة التالي كقيمة مميزة للسرعة لجعل بعض السرعات بلا أبعاد.

V=ل2β(Δتي)زνجير{\displaystyle \mathrm {V} ={\frac {L^{2}\beta (\Delta T)g}{\nu Gr}}}

تأثير رقم غراشوف على تدفق السوائل المختلفة

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

ملحوظات

  1. على الرغم من أن مصطلح "رقم غراشوف" كان مستخدمًا بالفعل، إلا أنه لم يُطلق عليه هذا الاسم إلا في حوالي عام 1921، أي بعد 28 عامًا من وفاة فرانز غراشوف. ولا يزال سبب تسمية هذه المجموعة باسمه غير واضح. [ 1 ]

مراجع

  1. ساندر، سي جيه؛ هولمان، جيه بي (1972). "فرانز غراشوف ورقم غراشوف". المجلة الدولية لانتقال الحرارة والكتلة . 15 (3): 562-563 . Bibcode : 1972IJHMT..15..562S . doi : 10.1016/0017-9310(72)90220-7 .
  2. 1 2 3 إنكروبيرا، فرانك (2007). أساسيات انتقال الحرارة والكتلة (الطبعة السادسة ). هوبوكين، نيوجيرسي: وايلي. الصفحات 408 ، 599، 629. ISBN   9780471457282. OCLC 288958608 . 
  3. 1 2 بيرد، ر. بايرون؛ ستيوارت، وارن إي.؛ لايتفوت، إدوين ن. (2002). ظواهر النقل (الطبعة الثانية ). نيويورك: جون وايلي. ص 318 ، 359. ISBN   9780471410775. OCLC 471520548 . 
  4. شاه، نهاد علي؛ أنيماسون، إلينوي؛ ابراهيم، ر.ع. باباتوندي، HA؛ سانديب، ن.؛ بوب، آي. (2018). "دراسة تأثير رقم جراشوف على جريان الموائع المختلفة المدفوعة بالحمل الحراري على الأسطح المختلفة". مجلة السوائل الجزيئية . 249 : 980– 990. دوى : 10.1016/j.molliq.2017.11.042 . ردمك 0167-7322 . 

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

  • سينجل، يونس أ. (2003). انتقال الحرارة والكتلة: منهج عملي (  الطبعة الثالثة). بوسطن: ماكجرو هيل.
  • إيكرت، إرنست آر جي ؛ دريك، روبرت إم. (1972). تحليل انتقال الحرارة والكتلة . نيويورك: ماكجرو هيل.
  • جالوريا، يوجيش (1980). انتقال الحرارة والكتلة بالحمل الحراري الطبيعي . نيويورك: دار بيرغامون للنشر.{{cite book}}: CS1 maint: موقع الناشر مفقود ( رابط )
  • ويلتي، جيمس ر. (1976). أساسيات نقل الزخم والحرارة والكتلة . نيويورك: جون وايلي وأولاده.