Tracer Diffusion in Granular Suspensions: Testing the Enskog Kinetic Theory with DSMC and Molecular Dynamics

यहाँ एंटोनियो एम. पुएर्टास और रूबेन गोमेज़ गोंजालेज के शोध पत्र "Tracer Diffusion in Granular Suspensions: Testing the Enskog Kinetic Theory with DSMC and Molecular Dynamics" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

यह शोध पत्र एक विस्कस कैरियर फ्लूइड (विलायक) में डूबे हुए ग्रेन्युलर गैस के भीतर एक एकल घुसपैठिया (ट्रेशर) कण के प्रसार (diffusion) की जांच करता है। यह प्रणाली, जिसे ग्रेन्युलर सस्पेंशन के रूप में जाना जाता है, एक जटिल भौतिक चुनौती प्रस्तुत करती है क्योंकि इसमें दो अलग-अलग तंत्रों का सह-अस्तित्व शामिल है:

1. विघटनकारी टकराव (Dissipative collisions): ठोस कणों के बीच अलोचदार (inelastic) टकराव, जिससे ऊर्जा की हानि होती है।
2. विलायक अंतःक्रिया (Solvent interaction): तरल एक विस्कस ड्रैग फोर्स (viscous drag force) और एक स्टोकेस्टिक (तापीय) बल लगाता है, जो एक थर्मोस्टेट के रूप में कार्य करता है।

केंद्रीय समस्या यह है कि इस गैर-संतुलन स्थिर अवस्था (non-equilibrium steady state) में ट्रेशर डिफ्यूजन गुणांक ($D$) की भविष्यवाणी करने के लिए एनस्कॉग काइनेटिक थ्योरी (विशेष रूप से दूसरे सोनाइन सन्निकटन तक चेपमैन-एनस्कॉग विस्तार) को मान्य करना है। जबकि पिछले सैद्धांतिक कार्यों ने विश्लेषणात्मक व्यंजक स्थापित किए थे, व्यापक रेंज में मापदंडों के विरुद्ध उनकी मजबूती—विशेष रूप से विलायक बलों की उपस्थिति में "आणविक अराजकता" (molecular chaos) की धारणा और मैक्सवेलियन सन्निकटन की वैधता के संबंध में—का व्यवस्थित रूप से परीक्षण नहीं किया गया था।

2. कार्यप्रणाली (Methodology)

लेखक एक बहु-आयामी दृष्टिकोण अपनाते हैं जो एक कठोर सत्यापन सुनिश्चित करने के लिए विश्लेषणात्मक सिद्धांत को दो अलग-अलग संख्यात्मक सिमुलेशन तकनीकों के साथ जोड़ता है:

• सैद्धांतिक ढांचा (Theoretical Framework):

  • सिस्टम को एक विलायक में चिकने हार्ड स्फेयर्स (ग्रेन्स और एक एकल घुसपैठिया) के द्विआधारी मिश्रण के रूप में मॉडल किया गया है।
  • विलायक अंतःक्रिया को एक लैंग्विन समीकरण (Langevin equation) के माध्यम से मॉडल किया गया है जिसमें एक रैखिक ड्रैग फोर्स ($-\gamma v$) और एक गाऊसी व्हाइट नॉइज़ टर्म (स्टोकेस्टिक बल) शामिल है जो फ्लक्चुएशन-डिसिपेशन प्रमेय का पालन करता है।
  • घुसपैठिया वितरण फलन (distribution function) के लिए एनस्कॉग-लोरेंत्ज़ काइनेटिक समीकरण व्युत्पन्न किया गया है।
  • वेग वितरण फलन को संभालने के लिए प्रथम और द्वितीय सोनाइन बहुपद सन्निकटन (Sonine polynomial approximations) का उपयोग करते हुए चेपमैन-एनस्कॉग (CE) विस्तार लागू किया गया है।

• संख्यात्मक सिमुलेशन (Numerical Simulations):

  • मॉलिक्यूलर डायनेमिक्स (MD): न्यूटन के गति के समीकरणों को हल करने वाला एक पूर्णतः नियतात्मक (deterministic) सिमुलेशन, जिसमें $N=1000$ (और 8000 तक) कण शामिल हैं। यह विधि सभी सूक्ष्म सहसंबंधों (microscopic correlations), परिमित-आकार प्रभावों (finite-size effects), और "आणविक अराजकता" परिकल्पना से विचलन को कैप्चर करती है। यह "ग्राउंड ट्रुथ" के रूप में कार्य करता है।
  • डायरेक्ट सिमुलेशन मोंटे कार्लो (DSMC): आणविक अराजकता की धारणा पर आधारित एक स्टोकेस्टिक एल्गोरिदम। यह एनस्कॉग समीकरण को संख्यात्मक रूप से हल करता है। DSMC, कण प्रक्षेपवक्रों (trajectories) की विशिष्ट गतिशीलता से सिद्धांत की मान्यताओं को अलग करने के लिए एक मध्यवर्ती संदर्भ के रूप में कार्य करता है।

• प्रमुख अवलोकन (Key Observables):

  • घुसपैठिया तापमान ($T_0$): ट्रेशर की स्थिर-अवस्था गतिज ऊर्जा।
  • वेलोसिटी ऑटोकोरिलेशन फंक्शन (VACF): $C_{vv}(t) = \langle \mathbf{v}(t) \cdot \mathbf{v}(0) \rangle$, जिसका उपयोग ग्रीन-कुबो संबंध के माध्यम से $D$ की गणना के लिए किया जाता है।
  • ट्रेशर डिफ्यूजन गुणांक ($D$): VACF के समय समाकल (time integral) और मीन स्क्वेयर्ड डिस्प्लेसमेंट (MSD) दोनों से गणना की गई है।

3. प्रमुख योगदान (Key Contributions)

• व्यवस्थित सत्यापन: यह शोध पत्र ग्रेन्युलर सस्पेंशन में ट्रेशर डिफ्यूजन के लिए एनस्कॉग-लैंग्विन काइनेटिक थ्योरी के पहले व्यापक संख्यात्मक सत्यापन को प्रदान करता है, जो विश्लेषणात्मक सिद्धांत और सूक्ष्म गतिशीलता के बीच के अंतर को पाटता है।
• दोहरी-सिमुलेशन रणनीति: MD (जिसमें सहसंबंध शामिल हैं) की तुलना DSMC (जो आणविक अराजकता मानता है) से करके, लेखक उन विशिष्ट क्षेत्रों की पहचान कर सकते हैं जहाँ सहसंबंधित कण गति के कारण काइनेटिक थ्योरी विफल हो जाती है।
• सोनाइन सन्निकटन का परिशोधन: अध्ययन प्रदर्शित करता है कि जबकि प्रथम सोनाइन सन्निकटन अक्सर अपर्याप्त होता है, द्वितीय सोनाइन सन्निकटन सैद्धांतिक भविष्यवाणियों की सटीकता में काफी सुधार करता है, जिससे वे DSMC और MD परिणामों के साथ निकट सहमति में आ जाते हैं।
• पैरामीटर स्पेस अन्वेषण: अध्ययन ने सिद्धांत की वैधता की सीमाओं को मैप करने के लिए प्रतिस्थापन गुणांक ($\alpha$), वॉल्यूम फ्रैक्शन ($\phi$), और घुसपैठिया द्रव्यमान अनुपात ($m_0/m$) को व्यवस्थित रूप से बदला है।

4. मुख्य परिणाम (Key Results)

A. घर्षण और ड्रैग का प्रभाव

• सिद्धांत कम-से-मध्यम ड्रैग गुणांकों के लिए डिफ्यूजन गुणांकों की सटीक भविष्यवाणी करता है।
• बहुत उच्च ड्रैग (ब्राउनियन मोशन द्वारा प्रभावी) की सीमा में, सिद्धांत टकरावों द्वारा प्रसार के दमन को पकड़ने में विफल रहता है, क्योंकि यहाँ सिंगल-पार्टिकल डिफ्यूजन लिमिट के करीब पहुँच जाता है।
• एक विशिष्ट क्षेत्र ($\gamma \sigma / \sqrt{m T_b} = 1$) की पहचान की गई जहाँ टकराव प्रभाव और विलायक प्रभाव संतुलित होते हैं, जो सिद्धांत के परीक्षण के लिए इष्टतम स्थितियाँ प्रदान करता है।

B. अलोचकता (Inelasticity - $\alpha$) का प्रभाव

• तापमान: बाथ तापमान ($T_b$) की तुलना में ग्रेन्युलर तापमान ($T$) कम होता है क्योंकि ऊर्जा का क्षय होता है। $T/T_b$ के लिए सैद्धांतिक भविष्यवाणियाँ, मजबूत अलोचकता ($\alpha \approx 0.3$) के लिए भी, MD और DSMC परिणामों के साथ उत्कृष्ट रूप से मेल खाती हैं।
• वेग वितरण: विलायक की उपस्थिति (स्टोकेस्टिक फोर्सिंग) उच्च-वेग वाली पूंछों (high-velocity tails) को दबा देती है जो अक्सर मुक्त रूप से ठंडे होने वाले ग्रेन्युलर गैसों में देखी जाती हैं। वेग वितरण मैक्सवेलियन के करीब रहता है, जो काइनेटिक थ्योरी की एक प्रमुख धारणा को मान्य करता है।
• डिफ्यूजन गुणांक:
  • सिद्धांत और सिमुलेशन उच्च घनत्व पर $\alpha$ के फलन के रूप में $D$ का गैर-मोनोटोनिक व्यवहार दिखाते हैं।
  • विसंगति: उच्च घनत्व और तीव्र अलोचकता पर, MD सिमुलेशन सैद्धांतिक भविष्यवाणी की तुलना में थोड़े उच्च डिफ्यूजन गुणांक दिखाते हैं। इसे सहसंबंधित कण गति (आणविक अराजकता का उल्लंघन) के कारण माना जाता है जिसे एनस्कॉग थ्योरी अनदेखा करती है।
  • समाधान: द्वितीय सोनाइन सन्निकटन प्रथम सन्निकटन की तुलना में इस विसंगति को काफी कम कर देता है।

C. घुसपैठिया द्रव्यमान ($m_0$) का प्रभाव

• प्रत्यास्थ मामला ($\alpha=1$): घुसपैठिया द्रव्यमान केवल क्षय के समय पैमाने को प्रभावित करता है; इक्विपार्टिशन बना रहता है ($T_0 = T_b$)।
• अलोचदार मामला ($\alpha < 1$): इक्विपार्टिशन टूट जाता है। भारी घुसपैठिए ($m_0 > m$) उच्च गतिज ऊर्जा बनाए रखते हैं और धीमी वेग डिकोरिलेशन (लंबी मेमोरी) प्रदर्शित करते हैं, जिससे हल्के घुसपैठियों की तुलना में उन्नत प्रसार होता है।
• सैद्धांतिक मॉडल पूरे रेंज ($0.01m$ से $100m$) में घुसपैठिया तापमान और डिफ्यूजन गुणांक दोनों की द्रव्यमान निर्भरता को सफलतापूर्वक पुनरुत्पादित करता है।

D. सहसंबंधित गति (Correlated Motion)

• कण विस्थापन के स्थानिक सहसंबंध ($\Delta(r)$) के विश्लेषण से पता चलता है कि जैसे-जैसे अलोचकता बढ़ती है, लंबी दूरी के सहसंबंध उभरते हैं।
• ये सहसंबंध ही MD परिणामों और एनस्कॉग थ्योरी के बीच विचलन का प्राथमिक स्रोत हैं। हालांकि, विलायक घर्षण ($\gamma$) को बढ़ाने से कणों की गति रैंडमाइज (randomize) हो जाती है, जिससे सहसंबंध दब जाते हैं और MD परिणामों को वापस थ्योरी के साथ सहमति में ले आता है।

5. महत्व और निष्कर्ष (Significance and Conclusion)

यह कार्य पुष्टि करता है कि एनस्कॉग काइनेटिक थ्योरी, जब विलायक बलों (लैंग्विन डायनेमिक्स) को शामिल करने के लिए विस्तारित की जाती है और इसे द्वितीय सोनाइन सन्निकटन का उपयोग करके हल किया जाता है, तो ग्रेन्युलर सस्पेंशन में ट्रेशर डिफ्यूजन का वर्णन करने के लिए एक मजबूत और सटीक ढांचा है।

• वैधता: यह सिद्धांत मध्यम घनत्व ($\phi \lesssim 0.2$) और अलोचकता की एक विस्तृत श्रृंखला के लिए अच्छी तरह से काम करता है, बशर्ते कि विलायक घर्षण सिस्टम को स्थिर अवस्था बनाए रखने और लंबी दूरी के सहसंबंधों को दबाने के लिए पर्याप्त हो।
• भौतिक अंतर्दृष्टि: अध्ययन स्पष्ट करता है कि विलायक न केवल एक थर्मोस्टेट के रूप में कार्य करता है बल्कि एक ऐसे तंत्र के रूप में भी कार्य करता है जो कणों के वेग को डिकोरिलेट करके "आणविक अराजकता" को बहाल करता है, जिससे शुष्क ग्रेन्युलर गैसों (dry granular gases) की तुलना में जहाँ क्लस्टरिंग और नॉन-गौसियन व्यवहार होता है, वहां भी काइनेटिक थ्योरी की प्रयोज्यता बढ़ जाती है।
• भविष्य की दिशाएं: लेखक सुझाव देते हैं कि यह प्रमाणित किया गया आइसोट्रोपिक बेसलाइन शीयर्ड ग्रेन्युलर सस्पेंशन (sheared granular suspensions) के भविष्य के अध्ययनों के लिए आवश्यक है, जहाँ अनिसोट्रोपिक ट्रांसपोर्ट टेंसर और गैर-संतुलन स्थिर अवस्थाओं की महत्वपूर्ण भूमिका होने की उम्मीद है।

संक्षेप में, यह शोध पत्र जटिल बहु-चरणीय प्रणालियों में सूक्ष्म कण गतिशीलता और स्थूल परिवहन सिद्धांत (macroscopic transport theory) के बीच के अंतर को सफलतापूर्वक पाटता है, जो औद्योगिक और प्राकृतिक तरल वातावरण में ग्रेन्युलर सामग्रियों के व्यवहार की भविष्यवाणी करने के लिए एक विश्वसनीय उपकरण प्रदान करता है।