نظرية تشابمان-إنسكوج

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

سميت النظرية نسبة إلى سيدني تشابمان وديفيد إنسكوج ، اللذين قدماها بشكل مستقل في عامي 1916 و1917. [ 1 ]

وصف

تُعد معادلة بولتزمان لدالة توزيع الجسيم الواحد نقطة البداية لنظرية تشابمان-إنسكوغو(ر،v،ت){\displaystyle f(\mathbf {r} ,\mathbf {v} ,t)}:

وت+vور+Fموv=ج^و،{\displaystyle {\frac {\partial f}{\partial t}}+\mathbf {v} \cdot {\frac {\partial f}{\partial \mathbf {r} }}+{\frac {\mathbf {F} }{m}}\cdot {\frac {\partial f}{\partial \mathbf {v} }}={\hat {C}}f,}

أينج^{\displaystyle {\hat {C}}}هو عامل تكاملي غير خطي يُمثل تطورو{\displaystyle f}تحت تأثير تصادمات الجسيمات. هذه اللاخطية تجعل حل معادلة بولتزمان الكاملة أمراً صعباً، وتحفز تطوير تقنيات تقريبية مثل تلك التي توفرها نظرية تشابمان-إنسكوغ.

انطلاقاً من هذه النقطة، فإن الافتراضات المختلفة التي تقوم عليها معادلة بولتزمان تنطبق أيضاً على نظرية تشابمان-إنسكوغ. وأبسط هذه الافتراضات يتطلب فصلاً في المقياس بين مدة التصادمτج{\displaystyle \tau _{\mathrm {c} }}ومتوسط ​​الوقت الحر بين التصادماتτو{\displaystyle \tau _{\mathrm {f} }}:τجτو{\displaystyle \tau _{\mathrm {c} }\ll \tau _{\mathrm {f} }}يضمن هذا الشرط أن تكون التصادمات أحداثًا محددة جيدًا في المكان والزمان، ويتحقق إذا كان المعامل عديم الأبعادγرج3ن{\displaystyle \gamma \equiv r_{\mathrm {c} }^{3}n}صغير، حيثرج{\displaystyle r_{\mathrm {c} }}هو نطاق التفاعلات بين الجسيمات ون{\displaystyle n}هي كثافة العدد . [ 2 ] بالإضافة إلى هذا الافتراض، تتطلب نظرية تشابمان-إنسكوج أيضًا أنτو{\displaystyle \tau _{\mathrm {f} }}أصغر بكثير من أي نطاقات زمنية خارجيةτext{\displaystyle \tau _{\text{ext}}}هذه هي الأطر الزمنية المرتبطة بالحدود الموجودة على الجانب الأيسر من معادلة بولتزمان، والتي تصف تغيرات حالة الغاز على أطوال ماكروسكوبية. عادةً ما تُحدد قيمها بواسطة الشروط الابتدائية/الحدودية و/أو المجالات الخارجية. هذا الفصل بين المقاييس يعني أن حد التصادم على الجانب الأيمن من معادلة بولتزمان أكبر بكثير من حدود التدفق على الجانب الأيسر. وبالتالي، يمكن إيجاد حل تقريبي من

ج^و=0.{\displaystyle {\hat {C}}f=0.}

يمكن إثبات أن حل هذه المعادلة هو توزيع غاوسي :

و=ن(ر،ت)(م2πكبتي(ر،ت))3/2خبرة[-م|v-v0(ر،ت)|22كبتي(ر،ت)]،{\displaystyle f=n(\mathbf {r} ,t)\left({\frac {m}{2\pi k_{\text{B}}T(\mathbf {r} ,t)}}\right)^{3/2}\exp \left[-{\frac {m{\left|\mathbf {v} -\mathbf {v} _{0}(\mathbf {r} ,t)\right|}^{2}}{2k_{\text{B}}T(\mathbf {r} ,t)}}\right],}

أينم{\displaystyle m}كتلة الجزيء وكب{\displaystyle k_{\text{B}}}هو ثابت بولتزمان . [ 3 ] يُقال إن الغاز في حالة توازن موضعي إذا حقق هذه المعادلة. [ 4 ] يؤدي افتراض التوازن الموضعي مباشرةً إلى معادلات أويلر ، التي تصف الموائع بدون تبديد، أي ذات موصلية حرارية ولزوجة تساوي0{\displaystyle 0}يتمثل الهدف الرئيسي لنظرية تشابمان-إنسكوغ في الحصول بشكل منهجي على تعميمات لمعادلات أويلر تتضمن التبديد. ويتحقق ذلك من خلال التعبير عن الانحرافات عن التوازن المحلي كسلسلة اضطرابية في عدد كنودسن.كن{\displaystyle {\text{Kn}}}وهو أمر صغير إذاτوτext{\displaystyle \tau _{\mathrm {f} }\ll \tau _{\text{ext}}}من الناحية النظرية، تصف المعادلات الهيدروديناميكية الناتجة التفاعل الديناميكي بين التدفق الحر والتصادمات بين الجسيمات. تميل الأخيرة إلى دفع الغاز نحو حالة التوازن الموضعي، بينما تعمل الأولى عبر عدم التجانس المكاني لدفع الغاز بعيدًا عن حالة التوازن الموضعي. [ 5 ] عندما يكون عدد كنودسن من رتبة 1 أو أكبر، لا يمكن وصف الغاز في النظام قيد الدراسة بأنه سائل.

للطلب الأول فيكن{\displaystyle {\text{Kn}}}نحصل على معادلات نافيير-ستوكس . وتؤدي المعادلات من الرتبة الثانية والثالثة، على التوالي، إلى معادلات بورنيت ومعادلات بورنيت الفائقة.

الصياغة الرياضية

بما أن عدد كنودسن لا يظهر صراحةً في معادلة بولتزمان، بل ضمنيًا من حيث دالة التوزيع والشروط الحدية، فإن متغيرًا وهميًاε{\displaystyle \varepsilon }تم إدخالها لتتبع الطلبات المناسبة في توسعة تشابمان-إنسكوج:

وت+vور+Fموv=1εج^و.{\displaystyle {\frac {\partial f}{\partial t}}+\mathbf {v\cdot } {\frac {\partial f}{\partial \mathbf {r} }}+{\frac {\mathbf {F} }{m}}\cdot {\frac {\partial f}{\partial \mathbf {v} }}={\frac {1}{\varepsilon }}{\hat {C}}f.}

صغيرε{\displaystyle \varepsilon }يشير إلى مصطلح التصادمج^و{\displaystyle {\hat {C}}f}يهيمن على مصطلح البثvور+Fموv{\displaystyle \mathbf {v\cdot } {\frac {\partial f}{\partial \mathbf {r} }}+{\frac {\mathbf {F} }{m}}\cdot {\frac {\partial f}{\partial \mathbf {v} }}}وهذا يعني أن عدد كنودسن صغير. وبالتالي، فإن الشكل المناسب لتوسيع تشابمان-إنسكوغ هو

و=و(0)+εو(1)+ε2و(2)+ .{\displaystyle f=f^{(0)}+\varepsilon f^{(1)}+\varepsilon ^{2}f^{(2)}+\cdots \ .}

تُعرف الحلول التي يمكن توسيعها رسميًا بهذه الطريقة بالحلول العادية لمعادلة بولتزمان. [ 6 ] يستثني هذا النوع من الحلول المساهمات غير الاضطرابية (مثلهـ-1/ε{\displaystyle e^{-1/\varepsilon }}والتي تظهر في الطبقات الحدية أو بالقرب من طبقات الصدمات الداخلية . وبالتالي، فإن نظرية تشابمان-إنسكوغ تقتصر على الحالات التي تكون فيها هذه الحلول ضئيلة.

باستبدال هذا التوسع ومساواة رتبε{\displaystyle \varepsilon }يؤدي إلى التسلسل الهرمي

ج(و(0)،و(0))=02ج(و(0)،و(ن))=(ت+vر+Fمv)و(ن-1)-م=1ن-1ج(و(م)،و(ن-م))،ن>0،{\displaystyle {\begin{aligned}J(f^{(0)},f^{(0)})&=0\\2J(f^{(0)},f^{(n)})&=\left({\frac {\partial }{\partial t}}+\mathbf {v\cdot } {\frac {\partial }{\partial \mathbf {r} }}+{\frac {\mathbf {F} }{m}}\cdot {\frac {\partial }{\partial \mathbf {v} }}\right)f^{(n-1)}-\sum _{m=1}^{n-1}J(f^{(m)},f^{(nm)}),\qquad n>0,\end{aligned}}}

أينج{\displaystyle J}هو مؤثر تكاملي، خطي في كلا وسيطيه، والذي يحققج(و،ز)=ج(ز،و){\displaystyle J(f,g)=J(g,f)}وج(و،و)=ج^و{\displaystyle J(f,f)={\hat {C}}f}حل المعادلة الأولى هو توزيع غاوسي:

و(0)=ن(ر،ت)(م2πكبتي(ر،ت))3/2خبرة[-م|v-v0(ر،ت)|22كبتي(ر،ت)].{\displaystyle f^{(0)}=n'(\mathbf {r} ,t)\left({\frac {m}{2\pi k_{\text{B}}T'(\mathbf {r} ,t)}}\right)^{3/2}\exp \left[-{\frac {m\left|\mathbf {v} -\mathbf {v} '_{0}(\mathbf {r} ,t)\right|^{2}}{2k_{\text{B}}T'(\mathbf {r} ,t)}}\right].}

بالنسبة لبعض الوظائفن(ر،ت){\displaystyle n'(\mathbf {r} ,t)}،v0(ر،ت){\displaystyle \mathbf {v} '_{0}(\mathbf {r} ,t)}، وتي(ر،ت){\displaystyle T'(\mathbf {r} ,t)}. التعبير عنو(0){\displaystyle f^{(0)}}يشير ذلك إلى وجود صلة بين هذه الوظائف والحقول الهيدروديناميكية الفيزيائية المحددة على أنها لحظاتو(ر،v،ت){\displaystyle f(\mathbf {r} ,\mathbf {v} ,t)}:

ن(ر،ت)=و(ر،v،ت)دvن(ر،ت)v0(ر،ت)=vو(ر،v،ت)دvن(ر،ت)تي(ر،ت)=م3كبv2و(ر،v،ت)دv.{\displaystyle {\begin{aligned}n(\mathbf {r} ,t)&=\int f(\mathbf {r} ,\mathbf {v} ,t)\,d\mathbf {v} \\n(\mathbf {r} ,t)\mathbf {v} _{0}(\mathbf {r} ,t)&=\int \mathbf {v} f(\mathbf {r} ,\mathbf {v} ,t)\,d\mathbf {v} \\n(\mathbf {r} ,t)T(\mathbf {r} ,t)&=\int {\frac {m}{3k_{\text{B}}}}v^{2}f(\mathbf {r} ,\mathbf {v} ,t)\,d\mathbf {v} .\end{aligned}}}

لكن من وجهة نظر رياضية بحتة، فإن مجموعتي الدوال ليستا بالضرورة متطابقتين بالنسبة لـε>0{\displaystyle \varepsilon >0}ε=0{\displaystyle \varepsilon =0}(هما متساويان بحكم التعريف). في الواقع، عند اتباع التسلسل الهرمي بشكل منهجي، يجد المرء أنه على غرارو(0){\displaystyle f^{(0)}}، كلو(ن){\displaystyle f^{(n)}}يحتوي أيضًا على وظائف اختيارية منر{\displaystyle \mathbf {r} }وت{\displaystyle t}والتي لا تُعرف علاقتها بالحقول الهيدروديناميكية الفيزيائية مسبقًا . أحد الافتراضات التبسيطية الرئيسية لنظرية تشابمان-إنسكوغ هو افتراض إمكانية كتابة هذه الدوال، التي تبدو اعتباطية، بدلالة الحقول الهيدروديناميكية الدقيقة وتدرجاتها المكانية. بعبارة أخرى، فإن اعتمادها على المكان والزمانو{\displaystyle f}يدخل ضمنيًا فقط من خلال الحقول الهيدروديناميكية. هذا البيان معقول فيزيائيًا لأن أرقام كنودسن الصغيرة تتوافق مع النظام الهيدروديناميكي، حيث تُحدد حالة الغاز فقط بواسطة الحقول الهيدروديناميكية. في حالةو(0){\displaystyle f^{(0)}}، الوظائفن(ر،ت){\displaystyle n'(\mathbf {r} ,t)}،v0(ر،ت){\displaystyle \mathbf {v} '_{0}(\mathbf {r} ,t)}، وتي(ر،ت){\displaystyle T'(\mathbf {r} ,t)}يُفترض أنها تساوي تمامًا الحقول الهيدروديناميكية الفيزيائية.

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

ن=1εنو(ن)دv=0=ن=1εنو(ن)v2دvن=1εنو(ن)vأنادv=0،أنا{x،y،z}.{\displaystyle {\begin{aligned}\int \sum _{n=1}^{\infty }\varepsilon ^{n}f^{(n)}\,d\mathbf {v} =0=\int \sum _{n=1}^{\infty }\varepsilon ^{n}f^{(n)}\mathbf {v} ^{2}\,d\mathbf {v} \\[1ex]\int \sum _{n=1}^{\infty }\varepsilon ^{n}f^{(n)}v_{i}\,d\mathbf {v} =0,\qquad i\in \{x,y,z\}.\end{aligned}}}

علاوة على ذلك، حتى لو وُجدت مثل هذه الحلول، يبقى السؤال الإضافي حول ما إذا كانت تشمل المجموعة الكاملة من الحلول الطبيعية لمعادلة بولتزمان، أي ما إذا كانت لا تمثل تقييدًا مصطنعًا للتوسع الأصلي فيε{\displaystyle \varepsilon }من أهم الإنجازات التقنية لنظرية تشابمان-إنسكوغ الإجابة على هذين السؤالين بالإيجاب. [ 6 ] وبالتالي، على الأقل على المستوى الرسمي، لا يوجد فقدان للعمومية في منهج تشابمان-إنسكوغ.

بعد تحديد هذه الاعتبارات الرسمية، يمكن للمرء أن يشرع في الحسابو(1){\displaystyle f^{(1)}}والنتيجة هي [ 1 ]

و(1)=[-1ن(2كبتيم)1/2أ(v)lnتي-2نب(v):v0]و(0)،{\displaystyle f^{(1)}=\left[-{\frac {1}{n}}\left({\frac {2k_{\text{B}}T}{m}}\right)^{1/2}\mathbf {A} (\mathbf {v} )\cdot \nabla \ln T-{\frac {2}{n}}\mathbb {B(\mathbf {v} )\colon \nabla } \mathbf {v} _{0}\right]f^{(0)},}

أينأ(v){\displaystyle \mathbf {A} (\mathbf {v} )}هو متجه وب(v){\displaystyle \mathbb {B} (\mathbf {v} )}موتر ، كل منها حل لمعادلة تكاملية خطية غير متجانسة يمكن حلها صراحةً بواسطة توسيع متعدد الحدود. هنا، تشير النقطتان الرأسيتان إلى الضرب النقطي المزدوج .تي:تي=أنا،جتيأناجتيجأنا{\textstyle \mathbb {T} :\mathbb {T'} =\sum _{i,j}T_{ij}T'_{ji}} للموتراتتي{\displaystyle \mathbb {T} }،تي{\displaystyle \mathbb {T'} }.

التوقعات

للحصول على الترتيب الأول في رقم كنودسن، فإن التدفق الحراريq=م2و(ر،v،ت)v2vدv{\textstyle \mathbf {q} ={\frac {m}{2}}\int f(\mathbf {r} ,\mathbf {v} ,t)\,v^{2}\mathbf {v} \,d\mathbf {v} }وُجد أنه يخضع لقانون فورييه لتوصيل الحرارة ، [ 7 ]

q=-λتي،{\displaystyle \mathbf {q} =-\lambda \nabla T,}

وموتر تدفق الزخمσ=م(v-v0)(v-v0)تيو(ر،v،ت)دv{\textstyle \mathbf {\sigma } =m\int (\mathbf {v} -\mathbf {v} _{0})(\mathbf {v} -\mathbf {v} _{0})^{\mathsf {T}}f(\mathbf {r} ,\mathbf {v} ,t)\,d\mathbf {v} }وهو ما يشبه السائل النيوتوني ، [ 7 ]

σ=صأنا-μ(v0+v0تي)+23μ(v0)أنا،{\displaystyle \mathbf {\sigma } =p\mathbb {I} -\mu \left(\nabla \mathbf {v_{0}} +\nabla \mathbf {v_{0}} ^{T}\right)+{\frac {2}{3}}\mu (\nabla \cdot \mathbf {v_{0}} )\mathbb {I} ,}

معأنا{\displaystyle \mathbb {I} }موتر الهوية. هنا،λ{\displaystyle \lambda }وμ{\displaystyle \mu }تمثل α الموصلية الحرارية وβ اللزوجة. ويمكن حسابهما بشكل صريح بدلالة المعاملات الجزيئية عن طريق حل معادلة تكاملية خطية؛ ويلخص الجدول أدناه النتائج لبعض النماذج الجزيئية المهمة (م{\displaystyle m}كتلة الجزيء وكب{\displaystyle k_{\text{B}}}(ثابت بولتزمان). [ 8 ]

الجدول 1: التعبيرات المتوقعة للتوصيل الحراري واللزوجة.
نموذجμ{\displaystyle \mu }λ{\displaystyle \lambda }ملحوظات
كرات مرنة صلبة قطرهاσ{\displaystyle \sigma }1.016516σ2(كبمتيπ)1/2{\displaystyle 1.016\cdot {\frac {5}{16\sigma ^{2}}}\left({\frac {k_{\text{B}}mT}{\pi }}\right)^{1/2}}2.52232كبمμ{\displaystyle 2.522\cdot {\frac {3}{2}}{\frac {k_{\text{B}}}{m}}\cdot \mu }صحح إلى 3 منازل عشرية.
جزيئات ذات قوة تنافريةκ/رν{\displaystyle \kappa /r^{\nu }}581أ2(ν)Γ(4-2ν-1)(كبمتيπ)1/2(2كبتيκ)2/(ν-1){\displaystyle {\frac {5}{8}}{\frac {1}{A_{2}(\nu )\Gamma {\left(4-{\frac {2}{\nu -1}}\right)}}}\left({\frac {k_{\text{B}}mT}{\pi }}\right)^{1/2}\left({\frac {2k_{\text{B}}T}{\kappa }}\right)^{2/(\nu -1)}}154كبمμ{\displaystyle {\frac {15}{4}}{\frac {k_{\text{B}}}{m}}\cdot \mu }Γ{\displaystyle \Gamma }يرمز إلى دالة غاما ، وأ2(ν){\displaystyle A_{2}(\nu )}هو عامل عددي. وقد ذكر تشابمان وكاولينغ عدة قيم لهذا العامل، على سبيل المثالأ2(5)=0.436{\displaystyle A_{2}(5)=0.436}وأ2(11)=0.319{\displaystyle A_{2}(11)=0.319}[ 9 ]
إمكانات لينارد-جونز :V(ر)=4ε[(σر)12-(σر)6]{\displaystyle V(r)=4\varepsilon \left[{\left({\frac {\sigma }{r}}\right)}^{12}-{\left({\frac {\sigma }{r}}\right)}^{6}\right]}516σ2(كبمتيπ)1/21دبليو1(2)(2){\displaystyle {\frac {5}{16\sigma ^{2}}}\left({\frac {k_{\text{B}}mT}{\pi }}\right)^{1/2}\cdot {\frac {1}{{\mathcal {W}}_{1}^{(2)}(2)}}}154كبمμ{\displaystyle {\frac {15}{4}}{\frac {k_{\text{B}}}{m}}\cdot \mu }دبليو1(2)(2){\displaystyle {\mathcal {W}}_{1}^{(2)}(2)}هي وظيفة منكبتي/ε{\displaystyle k_{\text{B}}T/\varepsilon }والتي يمكن حسابها عددياً. وهي تختلف من5.682{\displaystyle 5.682}لكبتي/ε=0.3{\displaystyle k_{\text{B}}T/\varepsilon =0.3}ل1.1738{\displaystyle 1.1738}لكبتي/ε=100{\displaystyle k_{\text{B}}T/\varepsilon =100}[ 10 ]

باستخدام هذه النتائج، يصبح من السهل الحصول على معادلات نافيير-ستوكس. وبأخذ عزوم السرعة لمعادلة بولتزمان، نحصل على معادلات التوازن الدقيقة للحقول الهيدروديناميكية.ن(ر،ت){\displaystyle n(\mathbf {r} ,t)}،v0(ر،ت){\displaystyle \mathbf {v} _{0}(\mathbf {r} ,t)}، وتي(ر،ت){\displaystyle T(\mathbf {r} ,t)}:

نت+(نv0)=0v0ت+v0v0-Fم+1نσ=0تيت+v0تي+23كبن(σ:v0+q)=0.{\displaystyle {\begin{aligned}{\frac {\partial n}{\partial t}}+\nabla \cdot \left(n\mathbf {v} _{0}\right)&=0\\{\frac {\partial \mathbf {v} _{0}}{\partial t}}+\mathbf {v} _{0}\cdot \nabla \mathbf {v} _{0}-{\frac {\mathbf {F} }{m}}+{\frac {1}{n}}\nabla \cdot \mathbf {\sigma } &=0\\{\frac {\partial T}{\partial t}}+\mathbf {v} _{0}\cdot \nabla T+{\frac {2}{3k_{\text{B}}n}}\left(\mathbf {\sigma :} \nabla \mathbf {v} _{0}+\nabla \cdot \mathbf {q} \right)&=0.\end{aligned}}}

كما في القسم السابق، تشير النقطتان الرأسيتان إلى حاصل الضرب النقطي المزدوج .تي:تي=أنا،جتيأناجتيجأنا{\textstyle \mathbb {T} :\mathbb {T'} =\sum _{i,j}T_{ij}T'_{ji}} . باستبدال تعبيرات تشابمان-إنسكوغ لـq{\displaystyle \mathbf {q} }وσ{\displaystyle \sigma }، ومن ثم نصل إلى معادلات نافيير-ستوكس.

مقارنة بالتجربة

من أهم تنبؤات نظرية تشابمان-إنسكوج أن اللزوجة،μ{\displaystyle \mu }وهي مستقلة عن الكثافة (يمكن ملاحظة ذلك لكل نموذج جزيئي في الجدول 1، ولكنها في الواقع مستقلة عن النموذج). تعود هذه النتيجة غير البديهية إلى جيمس كلارك ماكسويل ، الذي استنتجها عام 1860 بناءً على حجج حركية أبسط. [ 11 ] وقد تم التحقق منها تجريبياً بشكل جيد للغازات ذات الكثافات العادية.

الجدول 2: القيم المقاسة تجريبياً لـو=λ/μجv{\displaystyle f=\lambda /\mu c_{v}}بالنسبة للغازات النبيلة الخمسة الأولى. [ 12 ]
الهيليوم2.45
نيون2.52
الأرجون2.48
كريبتون2.535
زينون2.58

من ناحية أخرى، تتنبأ النظرية بأنμ{\displaystyle \mu }يعتمد ذلك على درجة الحرارة. بالنسبة للكرات المرنة الصلبة، يكون التناسب المتوقع هوμتي1/2{\displaystyle \mu \propto T^{1/2}}بينما تُظهر النماذج الأخرى عادةً تباينًا أكبر مع درجة الحرارة. على سبيل المثال، بالنسبة للجزيئات التي تتنافر مع بعضها البعض بقوةر-ν{\displaystyle \propto r^{-\nu }}القياس المتوقع هوμتيs{\displaystyle \mu \propto T^{s}}، أينs=1/2+2/(ν-1){\displaystyle s=1/2+2/(\nu -1)}أخذs=0.668{\displaystyle s=0.668}، بما يتوافق معν12.9{\displaystyle \nu \approx 12.9}يُظهر توافقًا معقولًا مع القياسات التجريبية المُلاحظة للهيليوم. أما بالنسبة للغازات الأكثر تعقيدًا، فإن التوافق ليس جيدًا، ويرجع ذلك على الأرجح إلى إهمال قوى التجاذب. [ 13 ] في الواقع، يمكن جعل نموذج لينارد-جونز ، الذي يتضمن قوى التجاذب، أكثر توافقًا مع التجربة (وإن كان ذلك على حساب زيادة الغموض).تي{\displaystyle T}(انظر مدخل لينارد-جونز في الجدول 1). [ 14 ] وللحصول على توافق أفضل مع البيانات التجريبية مقارنةً بما تم الحصول عليه باستخدام نموذج لينارد-جونز ، تم استخدام جهد مي الأكثر مرونة ، [ 15 ] تسمح المرونة الإضافية لهذا الجهد بالتنبؤ الدقيق بخصائص النقل لخلائط مجموعة متنوعة من الجزيئات المتناظرة كرويًا.

تتنبأ نظرية تشابمان-إنسكوج أيضاً بعلاقة بسيطة بين الموصلية الحرارية،λ{\displaystyle \lambda }واللزوجة،μ{\displaystyle \mu }، بالشكلλ=وμجv{\displaystyle \lambda =f\mu c_{v}}، أينجv{\displaystyle c_{v}}هي الحرارة النوعية عند حجم ثابت وو{\displaystyle f}هو عامل عددي بحت. بالنسبة للجزيئات المتناظرة كرويًا، من المتوقع أن تكون قيمته قريبة جدًا من2.5{\displaystyle 2.5}بطريقة تعتمد قليلاً على النموذج. على سبيل المثال، الكرات المرنة الصلبة لهاو2.522{\displaystyle f\approx 2.522}والجزيئات ذات القوة التنافريةر-13{\displaystyle \propto r^{-13}}يملكو2.511{\displaystyle f\approx 2.511}(تم تجاهل الانحراف الأخير في الجدول 1). الحالة الخاصة لجزيئات ماكسويل (قوة التنافر)ر-5{\displaystyle \propto r^{-5}}) لديهو=2.5{\displaystyle f=2.5} بالضبط. [ 16 ] منذλ{\displaystyle \lambda }،μ{\displaystyle \mu }، وجv{\displaystyle c_{v}}يمكن قياسها مباشرة في التجارب، ويتمثل أحد الاختبارات التجريبية البسيطة لنظرية تشابمان-إنسكوج في قياسو{\displaystyle f}بالنسبة للغازات النبيلة ذات التناظر الكروي . يوضح الجدول 2 وجود توافق معقول بين النظرية والتجربة. [ 12 ]

الإضافات

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

يمكن أيضًا تطبيق النظرية على رتبة أعلى في عدد كنودسن. وعلى وجه الخصوص، المساهمة من الرتبة الثانيةو(2){\displaystyle f^{(2)}}تم حسابها بواسطة بورنيت. [ 17 ] مع ذلك، في الظروف العامة، قد لا تُحسّن هذه التصحيحات ذات الرتبة الأعلى نظرية الرتبة الأولى تحسينًا موثوقًا، نظرًا لأن متسلسلة تشابمان-إنسكوغ لا تتقارب دائمًا. [ 18 ] (من ناحية أخرى، يُعتقد أن المتسلسلة على الأقل مقاربة لحلول معادلة بولتزمان، وفي هذه الحالة، يُعطي الاقتطاع عند الرتبة المنخفضة نتائج دقيقة.) [ 19 ] حتى لو حسّنت التصحيحات ذات الرتبة الأعلى نظامًا معينًا، فإن تفسير المعادلات الهيدروديناميكية المقابلة لا يزال محل نقاش. [ 20 ]

نظرية إنسكوج المنقحة

تم توسيع نظرية تشابمان-إنسكوغ للخلائط متعددة المكونات لتشمل الكثافات المرتفعة، ولا سيما الكثافات التي يكون فيها الحجم المشترك للخليط غير مهمل، في سلسلة من الأعمال التي قام بها إي جي دي كوهين وآخرون [ 21 ] [ 22 ] [ 23 ] [ 24 ] [ 25 ] ، وأُطلق عليها اسم نظرية إنسكوغ المنقحة (RET). وقد جاء الاشتقاق الناجح لنظرية إنسكوغ المنقحة بعد عدة محاولات سابقة في نفس المجال، والتي أسفرت عن نتائج تبين أنها غير متسقة مع الديناميكا الحرارية غير العكوسية . وتُعد الصيغة المعدلة لمعادلة بولتزمان نقطة الانطلاق لتطوير نظرية إنسكوغ المنقحة.s{\displaystyle s}دالة توزيع سرعة الجسيمات،

(ت+vأنار+Fأنامأناvأنا)وأنا=جSأناج(وأنا،وج){\displaystyle \left({\frac {\partial }{\partial t}}+\mathbf {v} _{i}\cdot {\frac {\partial }{\partial \mathbf {r} }}+{\frac {\mathbf {F} _{i}}{m_{i}}}\cdot {\frac {\partial }{\partial \mathbf {v} _{i}}}\right)f_{i}=\sum _{j}S_{ij}(f_{i},f_{j})}

أينvأنا(ر،ت){\displaystyle \mathbf {v} _{i}(\mathbf {r} ,t)}هي سرعة جسيمات النوعأنا{\displaystyle i}، في الموقعر{\displaystyle \mathbf {r} }والوقتت{\displaystyle t}،مأنا{\displaystyle m_{i}}كتلة الجسيم،Fأنا{\displaystyle \mathbf {F} _{i}}هي القوة الخارجية، و

Sأناج(وأنا،وج)=[زأناج(σأناجك)وأنا(ر)وج(ر+σأناجك)-زأناج(-σأناجك)وأنا(ر)وج(ر-σأناجك)]دτ{\displaystyle S_{ij}(f_{i},f_{j})=\iiint \left[g_{ij}(\sigma _{ij}\mathbf {k} )\,f_{i}'(\mathbf {r} )\,f_{j}'(\mathbf {r} +\sigma _{ij}\mathbf {k} )-g_{ij}(-\sigma _{ij}\mathbf {k} )\,f_{i}(\mathbf {r} )\,f_{j}(\mathbf {r} -\sigma _{ij}\mathbf {k} )\right]d\tau }

يكمن الاختلاف في هذه المعادلة عن نظرية تشابمان-إنسكوغ الكلاسيكية في عامل التدفق.Sأناج{\displaystyle S_{ij}}، والتي يتم من خلالها تقييم توزيع سرعة الجسيمين عند نقاط مختلفة في الفضاء، يفصل بينهماσأناجك{\displaystyle \sigma _{ij}\mathbf {k} }، أينك{\displaystyle \mathbf {k} }يمثل متجه الوحدة على طول الخط الواصل بين مركزي كتلة الجسيمين. وينشأ اختلاف هام آخر من إدخال العوامل.زأناج{\displaystyle g_{ij}}والتي تمثل احتمالية التصادمات المتزايدة بسبب الحجم المستبعد. يتم استعادة معادلات تشابمان-إنسكوغ الكلاسيكية عن طريق وضعσأناج=0{\displaystyle \sigma _{ij}=0}وزأناج(σأناجك)=1{\displaystyle g_{ij}(\sigma _{ij}\mathbf {k} )=1}.

يُعد اختيار العوامل نقطة مهمة لنجاح اختبار RETزأناج{\displaystyle g_{ij}}، والتي تُفسر على أنها دالة توزيع الأزواج المحسوبة عند مسافة التلامسσأناج{\displaystyle \sigma _{ij}}من العوامل المهمة التي يجب ملاحظتها هنا أنه من أجل الحصول على نتائج تتوافق مع الديناميكا الحرارية غير العكوسة ، فإنزأناج{\displaystyle g_{ij}}يجب التعامل معها كدوال لحقول الكثافة، بدلاً من كونها دوال للكثافة المحلية.

نتائج من نظرية إنسكوج المنقحة

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

صنكتي=1+2πن3أناجxأناxجσأناج3زأناج،{\displaystyle {\frac {p}{nkT}}=1+{\frac {2\pi n}{3}}\sum _{i}\sum _{j}x_{i}x_{j}\sigma _{ij}^{3}g_{ij},}

وهو ما يتوافق مع معادلة حالة كارناهان-ستارلينغ ، ويتحول إلى قانون الغاز المثالي في حالة التخفيف اللانهائي (أي عندمانأنا،جxأناxجσأناج31{\textstyle n\sum _{i,j}x_{i}x_{j}\sigma _{ij}^{3}\ll 1})

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

λ=(1+نαλ)λ0+ن2تي1/2λσ{\displaystyle \lambda =(1+n\alpha _{\lambda })\lambda _{0}+n^{2}T^{1/2}\lambda _{\sigma }}

أينαλ{\displaystyle \alpha _{\lambda }}وλσ{\displaystyle \lambda _{\sigma }}وهي دوال ضعيفة نسبياً للتركيب ودرجة الحرارة والكثافة، وλ0{\displaystyle \lambda _{0}}هي الموصلية الحرارية التي تم الحصول عليها من نظرية تشابمان-إنسكوج الكلاسيكية.

وبالمثل، يمكن كتابة التعبير الذي تم الحصول عليه للزوجة على النحو التالي:

μ=(1+نتيαμ)μ0+ن2تي1/2μσ{\displaystyle \mu =(1+nT\alpha _{\mu })\mu _{0}+n^{2}T^{1/2}\mu _{\sigma }}

معαμ{\displaystyle \alpha _{\mu }}وμσ{\displaystyle \mu _{\sigma }}الدوال الضعيفة للتركيب ودرجة الحرارة والكثافة، وμ0{\displaystyle \mu _{0}}القيمة المستمدة من نظرية تشابمان-إنسكوج الكلاسيكية.

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

د1ن،دتي1ن{\displaystyle D\sim {\frac {1}{n}},\quad D_{T}\sim {\frac {1}{n}}}

لجميع الكثافات، بل يتنبأ بأن المعاملات ستنخفض ببطء أكبر مع الكثافة عند الكثافات العالية، وهو ما يتوافق جيدًا مع التجارب. كما أن هذه التغيرات في اعتماد الكثافة تؤدي إلى تنبؤ RET باعتماد معامل سوريه على الكثافة .

Sتي=دتيد،(Sتين)تي0،{\displaystyle S_{T}={\frac {D_{T}}{D}},\quad \left({\frac {\partial S_{T}}{\partial n}}\right)_{T}\neq 0,}

بينما تتنبأ نظرية تشابمان-إنسكوج الكلاسيكية بأن معامل سوريه، مثل اللزوجة والتوصيل الحراري، مستقل عن الكثافة.

التطبيقات

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

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

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

انظر أيضاً

ملحوظات

  1. 1 2 تشابمان، سيدني؛ كاولينج، تي جي (1970)، النظرية الرياضية للغازات غير المتجانسة (  الطبعة الثالثة)، مطبعة جامعة كامبريدج
  2. باليسكو، رادو (1975)، الميكانيكا الإحصائية في حالات التوازن وعدم التوازن ، جون وايلي وأولاده، ISBN 978-0-471-04600-4
  3. سيرسينياني، كارلو (1975)، نظرية وتطبيق معادلة بولتزمان ، إلسيفير، ص 78-79 ، ISBN  978-0-444-19450-3
  4. باليسكو، ص 450
  5. باليسكو، ص 451
  6. 1 2 غراد، هارولد ( 1958)، "مبادئ النظرية الحركية للغازات"، في فلوغه، س. (محرر)، موسوعة الفيزياء ، المجلد الثاني عشر، سبرينغر-فيرلاغ، الصفحات 205-294  
  7. 1 2 بيرد، ر. برايون؛ أرمسترونج، روبرت سي؛ هاساجر، أولي (1987)، ديناميكا السوائل البوليمرية، المجلد 1: ميكانيكا الموائع ( الطبعة الثانية)، جون وايلي وأولاده، ص 10-11  
  8. تشابمان وكاولينغ، الفصل 10
  9. تشابمان وكاولينغ، ص 172
  10. تشابمان وكاولينغ، ص 185
  11. ماكسويل، جيمس (1860)، "5. توضيحات لنظرية الديناميكا للغازات. - الجزء الأول. حول حركات واصطدامات الكرات المرنة تمامًا"، المجلة الفلسفية ، 19 (124): 19-32 ، doi : 10.1080/14786446008642818
  12. 1 2 تشابمان وكاولينغ ص 249
  13. تشابمان وكاولينغ، الصفحات 230-232
  14. تشابمان وكاولينغ، الصفحات 235-237
  15. 1 2 3 جيرفيل، فيغارد ج.؛ ويلهلمسن، أويفيند (2023-06-08). "نظرية إنسكوغ المنقحة لسوائل مي: التنبؤ بمعاملات الانتشار، ومعاملات الانتشار الحراري، واللزوجة، والتوصيل الحراري" . مجلة الفيزياء الكيميائية . 158 (22). doi : 10.1063/5.0149865 . ISSN 0021-9606 . 
  16. تشابمان وكاولينغ، ص 247
  17. بورنيت، د. (1936)، "توزيع السرعات الجزيئية والحركة المتوسطة في غاز غير منتظم"، وقائع الجمعية الرياضية في لندن ، 40 : 382، doi : 10.1112/plms/s2-40.1.382
  18. سانتوس، أندريس؛ بري، ج. خافيير؛ دفتي، جيمس و. (1986)، "تباعد توسيع تشابمان-إنسكوغ"، رسائل المراجعة الفيزيائية ، 56 (15): 1571-1574 ، Bibcode : 1986PhRvL..56.1571S ، doi : 10.1103/PhysRevLett.56.1571 ، PMID 10032711 
  19. غراد، هارولد (1963)، "النظرية التقاربية لمعادلة بولتزمان"، فيزياء الموائع ، 6 (2): 147، Bibcode : 1963PhFl....6..147G ، doi : 10.1063/1.1706716
  20. غارسيا-كولين، إل إس؛ فيلاسكو، آر إم؛ أوريبي، إف جيه (2008)، "ما وراء معادلات نافيير-ستوكس: ديناميكا بورنيت المائية"، تقارير الفيزياء ، 465 (4): 149-189 ، رمز Bibcode : 2008PhR...465..149G ، doi : 10.1016/j.physrep.2008.04.010
  21. لوبيز دي هارو، م.؛ كوهين، إي. جي. دي.؛ كينكيد، ج. م. (1983-03-01). "نظرية إنسكوغ للخلائط متعددة المكونات. الجزء الأول: نظرية النقل الخطي" . مجلة الفيزياء الكيميائية . 78 (5): 2746-2759 . doi : 10.1063/1.444985 . ISSN 0021-9606 . 
  22. كينكيد، جيه إم؛ لوبيز دي هارو، إم؛ كوهين، إي جي دي (1983-11-01). "نظرية إنسكوغ للخلائط متعددة المكونات. الجزء الثاني: الانتشار المتبادل" . مجلة الفيزياء الكيميائية . 79 (9): 4509-4521 . doi : 10.1063/1.446388 . ISSN 0021-9606 . 
  23. لوبيز دي هارو، م.؛ كوهين، إي جي دي (1984-01-01). "نظرية إنسكوغ للخلائط متعددة المكونات. الجزء الثالث: خصائص النقل للخلائط الثنائية الكثيفة ذات المكون المتتبع الواحد" . مجلة الفيزياء الكيميائية . 80 (1): 408-415 . doi : 10.1063/1.446463 . ISSN 0021-9606 . 
  24. كينكيد، جيه إم؛ كوهين، إي جي دي؛ لوبيز دي هارو، إم. (15 يناير 1987). "نظرية إنسكوغ للخلائط متعددة المكونات. الجزء الرابع: الانتشار الحراري" . مجلة الفيزياء الكيميائية . 86 (2): 963-975 . doi : 10.1063/1.452243 . ISSN 0021-9606 . 
  25. فان بيجيرين، هـ.؛ إرنست، م.هـ. (مارس 1973). "معادلة إنسكوغ-بولتزمان غير الخطية" . رسائل الفيزياء أ . 43 (4): 367-368 . doi : 10.1016/0375-9601(73)90346-0 .

مراجع

الدراسة الكلاسيكية حول هذا الموضوع:

  • تشابمان، سيدني؛ كاولينج، تي جي (1970)، النظرية الرياضية للغازات غير المتجانسة (  الطبعة الثالثة)، مطبعة جامعة كامبريدج

يحتوي على مقدمة تقنية للحلول العادية لمعادلة بولتزمان:

  • غراد، هارولد ( 1958)، "مبادئ النظرية الحركية للغازات"، في فلوغه، س. (محرر)، موسوعة الفيزياء ، المجلد  الثاني عشر، سبرينغر-فيرلاغ، الصفحات 205-294