🔬 condensed matter

Rheology of dense suspensions of granular spherocylinders by particle-based simulation

यह शोध पत्र एक कण-आधारित सिमुलेशन मॉडल प्रस्तुत करता है जो शियर फ्लो (shear flow) के तहत ग्रेन्युलर स्फेरोसिलिंडर्स (granular spherocylinders) के सघन सस्पेंशन की श्यानता (viscosity) और सूक्ष्मसंरचना (microstructure) की सफलतापूर्वक भविष्यवाणी करता है, जो यह प्रकट करता है कि वॉल्यूम फ्रैक्शन और पहलू अनुपात (aspect ratio) के साथ ये गुण कैसे विकसित होते हैं और सीमित प्रयोगात्मक डेटा की पुष्टि करते हैं।

मूल लेखक: Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

प्रकाशित 2026-07-31

📖 5 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

एक ऐसी दुनिया की कल्पना करें जहाँ आप जो कुछ भी डालते, पंप करते या मिलाते हैं, वह केवल पानी या रेत नहीं है, बल्कि तरल में तैरती हुई छोटी, कठोर छड़ियों का एक अराजक मिश्रण (chaotic soup) है। यह रियोलॉजी (rheology) का क्षेत्र है, जो इस बात का विज्ञान है कि चीजें कैसे बहती हैं। आपने संभवतः इसे क्रिया में देखा होगा: जैसे गिरे हुए पेड़ों से भरा हुआ एक नदी का हिस्सा जो एक विशाल "लॉग जैम" (लकड़ियों का जाम) बना देता है, या कागज बनाने के लिए उपयोग किया जाने वाला गाढ़ा, चिपचिपा घोल। इन स्थितियों में, कण इतने घने होते हैं कि वे लगातार एक-दूसरे से टकराते हैं, और उनका आकार उतना ही महत्वपूर्ण होता है जितना कि उनकी संख्या। यदि आप कंचों (marbles) से भरी बाल्टी को हिलाने की कोशिश करते हैं, तो वे आसानी से एक-दूसरे के पास से लुढ़क जाते हैं। लेकिन यदि आप लंबी, पतली छड़ियों से भरी बाल्टी को हिलाने की कोशिश करते हैं, तो वे आपस में उलझ जाती हैं, एक साथ लॉक हो जाती हैं, और अचानक पूरा मिश्रण एक ठोस की तरह व्यवहार करने लगता है। वैज्ञानिक लंबे समय से समझते आए हैं कि गोल गेंदें इन भीड़ में कैसे व्यवहार करती हैं, लेकिन लंबी, छड़ के आकार के कणों के लिए नियम थोड़े रहस्यमय रहे हैं, खासकर जब वे बहुत सघन रूप से पैक होते हैं। इस बात को समझना अत्यंत महत्वपूर्ण है क्योंकि ये "चिपचिपी छड़ियों" वाले मिश्रण हर जगह दिखाई देते हैं, औद्योगिक विनिर्माण से लेकर ज्वालामुखियों के भीतर पिघली हुई चट्टानों के प्रवाह तक।

इस शोध पत्र में, शोधकर्ताओं की एक टीम ने इन घनी छड़ वाले सस्पेंशन के लिए एक वर्चुअल प्रयोगशाला के रूप में एक अत्यंत स्मार्ट कंप्यूटर सिमुलेशन बनाया। लैब में वास्तविक रसायनों को मिलाने के बजाय, उन्होंने हजारों "स्फेरोसिलिंडर्स" (spherocylinders)—ऐसे कण जो कैप्सूल या गोलियों की तरह दिखते हैं (एक सिलेंडर जिसके दोनों सिरों पर अर्धगोला है)—से भरी एक डिजिटल दुनिया बनाई। उन्होंने इन कैप्सूलों को एक आभासी तरल में तैराया और फिर उन्हें "शीयर" (shear) करना शुरू किया, जो कि एक फैंसी तरीका है यह कहने का कि उन्होंने तरल की ऊपरी परत को बगल की ओर खींचा जबकि निचली परत को स्थिर रखा, ठीक वैसे ही जैसे टोस्ट पर मक्खन फैलाया जाता है। लक्ष्य यह देखना था कि जैसे-जैसे वे अधिक छड़ियों को स्थान में पैक करते हैं और जैसे-जैसे वे छड़ियों को लंबा करते हैं, मिश्रण का गाढ़ापन (viscosity) कैसे बदलता है।

सिमुलेशन ने कुछ दिलचस्प व्यवहारों का खुलासा किया। सबसे पहले, जब उन्होंने 'शीयर' शुरू किया, तो मिश्रण एक क्षण के लिए अविश्वसनीय रूप से गाढ़ा हो गया—एक "विस्कोसिटी स्पाइक" (viscosity spike)—इससे पहले कि वह स्थिर हो जाए। ऐसा इसलिए हुआ क्योंकि छड़ें शुरू में यादृच्छिक दिशाओं में बिखरी हुई थीं, जो एक 'मोश पिट' (mosh pit) की तरह एक-दूसरे से टकरा रही थीं। एक बार प्रवाह शुरू होने के बाद, छड़ें एक कतार में मार्च करने वाले सैनिकों की तरह एक पंक्ति में व्यवस्थित होने लगीं, और मिश्रण को हिलाना आसान हो गया। शोधकर्ताओं ने पाया कि वे जितनी अधिक छड़ियों को पैक करते हैं (आयतन अंश, $\phi$ ), मिश्रण उतना ही गाढ़ा होता जाता है, और अंततः एक ऐसे बिंदु पर पहुँच जाता है जहाँ यह जाम हो जाता है और पूरी तरह से बहना बंद कर देता है। दिलचस्प बात यह है कि छड़ें जितनी लंबी थीं (पहलू अनुपात, $A$ , जो 1 के गोले से लेकर 20 की बहुत लंबी छड़ियों तक है), जाम होने की प्रक्रिया उतनी ही जल्दी हुई। उदाहरण के लिए, लंबाई 5 वाली छड़ियों के लिए, जैमिंग पॉइंट आयतन अंश 0.57 के आसपास था, लेकिन लंबाई 20 वाली छड़ियों के लिए, मिश्रण लगभग 0.41 के बहुत कम घनत्व पर जाम हो गया।

टीम ने यह भी खोजा कि छड़ें केवल पंक्ति में ही नहीं आतीं; वे विशेष रूप से प्रवाह के साथ संरेखित (align) होती हैं। उन्होंने संरेखण को मापने के लिए एक "ऑर्डर पैरामीटर" ( $S$ ) का उपयोग किया। कम घनत्व पर, छड़ें डगमगाती थीं और लूप में घूमती थीं (जिसे जेफ्री ऑर्बिट्स/Jeffery orbits कहा जाता है), लेकिन जैसे-जैसे भीड़ घनी होती गई, वे प्रवाह की दिशा में एक सीधी रेखा में लॉक हो गईं। हालांकि, यदि भीड़ बहुत अधिक घनी हो गई, तो छड़ें वास्तव में परेशान होने लगीं और पूरी तरह से संरेखित रहने में असमर्थ रहीं, जिससे संरेखण (order) थोड़ा गिर गया। शोधकर्ताओं ने काम करने वाले बलों का भी विश्लेषण किया, यह दिखाते हुए कि जब मिश्रण पतला था, तो तरल का अपना प्रतिरोध सबसे महत्वपूर्ण था, लेकिन जैसे-जैसे यह घना होता गया, छड़ियों के बीच भौतिक टकराव और घर्षण प्रमुख बल बन गए, जिसने मिश्रण को गाढ़ा बनाने का कार्य संभाल लिया।

संपर्क, घर्षण और तरल ड्रैग (fluid drag) के भौतिकी को संतुलित करने वाला एक मॉडल बनाकर, लेखकों ने यह अनुमान लगाने के लिए एक नया उपकरण प्रदान किया कि ये जटिल मिश्रण कैसे व्यवहार करते हैं। उन्होंने केवल अनुमान नहीं लगाया; उन्होंने अलग-अलग कणों की संख्या, अलग-अलग छड़ की लंबाई और अलग-अलग पैकिंग घनत्व के साथ हजारों सिमुलेशन चलाए। उनके परिणाम बताते हैं कि हालांकि "अधिक सामग्री = गाढ़ा तरल" का मूल विचार सही है, लेकिन कणों का आकार उस बिंदु को नाटकीय रूप से बदल देता है जिस पर मिश्रण लॉक हो जाता है। यह कार्य यह दावा नहीं करता कि उन्होंने द्रव यांत्रिकी (fluid mechanics) के हर रहस्य को सुलझा लिया है, लेकिन यह नदियों में लॉग जैम को रोकने से लेकर बेहतर औद्योगिक स्लरी (slurry) डिजाइन करने तक, वास्तविक दुनिया में इन कठिन, छड़-भरे तरल पदार्थों को प्रबंधित करने के लिए एक ठोस, सिम्युलेटेड आधार प्रदान करता है।

तकनीकी सारांश: कण-आधारित सिमुलेशन द्वारा सघन गोलाकार बेलनाकार कणों (Granular Spherocylinders) की श्यानता (Rheology)

समस्या विवरण
छड़ के आकार के कणों के सघन निलंबन (Dense suspensions) प्राकृतिक और औद्योगिक संदर्भों में व्यापक रूप से मौजूद हैं, जिनमें नदियों में लकड़ी के जाम से लेकर मैग्मा में क्रिस्टलीय निलंबन और कागज निर्माण में फाइबर नेटवर्क तक शामिल हैं। जबकि गोलाकार कणों की श्यानता (rheology) अच्छी तरह से परिभाषित है, विशेष रूप से बॉयलर (Boyer) आदि द्वारा गैर-ब्राउनियन निलंबन और शुष्क कण प्रवाह के एकीकरण के माध्यम से, लम्बे कणों के व्यवहार को पूरी तरह से समझा नहीं जा सका है। कण अभिविन्यास (orientation) से जुड़ी अतिरिक्त स्वतंत्रता (degree of freedom) जटिल सूक्ष्म संरचनात्मक घटनाओं को जन्म देती है, जैसे कि जेफ़री ऑर्बिट्स (Jeffery orbits) और शियर-प्रेरित संरेखण (shear-induced alignment), जो प्रभावी श्यानता और स्ट्रेस रिस्पॉन्स को महत्वपूर्ण रूप से बदल देते हैं। सघन शासन (dense regimes) (उच्च ठोस आयतन अंश, $\phi$ ) के लिए मौजूदा प्रयोगात्मक डेटा तैयारी और माप की कठिनाइयों के कारण सीमित है। इसके अलावा, पिछले सिमुलेशन दृष्टिकोणों में ट्रेड-ऑफ (trade-offs) देखने को मिलते हैं: स्टोक्सियन डायनेमिक्स (Stokesian Dynamics - SD) कई-शरीर हाइड्रोडायनामिक्स के लिए आवश्यक मैट्रिक्स व्युत्क्रमों (matrix inversions) के कारण सघन प्रणालियों के लिए अत्यधिक गणनात्मक रूप से महंगा है, जबकि डिस्क्रीट एलीमेंट मेथड (Discrete Element Method - DEM) मॉडल अक्सर कण ज्यामिति (जैसे, गोलों की श्रृंखला) का सरलीकरण करते हैं या आवश्यक हाइड्रोडायनामिक बलों की उपेक्षा करते हैं। सघन, गैर-ब्राउनियन गोलाकार बेलनाकार निलंबन के प्रभुत्व वाले सूक्ष्म भौतिकी को सटीक रूप से पकड़ने के लिए एक गणनात्मक रूप से सुलभ कण-आधारित मॉडल की आवश्यकता है।

कार्यप्रणाली (Methodology)
लेखकों ने सरल शियर फ्लो (simple shear flow) के तहत गैर-बेंडेबल, गैर-ब्राउनियन गोलाकार बेलनाकार कणों (अर्धगोलाकार कैप वाले सिलेंडर) के सघन निलंबन के लिए एक कण-आधारित सिमुलेशन मॉडल विकसित किया है। यह मॉडल उस शासन में कार्य करता है जहाँ कण जड़त्व (particle inertia) मौजूद है लेकिन नगण्य है।

गति के समीकरण (Equations of Motion): कणों के प्रक्षेपवक्र (trajectories) को वेलोसिटी-वेलेट (Velocity-Verlet) एल्गोरिदम का उपयोग करके अपडेट किया जाता है, जो स्थानांतरीय (translational) और घूर्णी (rotational) गति के समीकरणों को एकीकृत करता है। ये समीकरण युग्म-वार संपर्क बल (pairwise contact forces), युग्म-वार लुब्रिकेशन बल, और एकल-शरीर हाइड्रोडायनामिक बलों (ड्रैग, लिफ्ट और टॉर्क) को ध्यान में रखते हैं।
हाइड्रोडायनामिक्स: एकल-शरीर बलों को गोलाकार बेलनाकार कणों के अनुमान के रूप में प्रोलैट स्फेरॉइड्स (prolate spheroids) के विश्लेषणात्मक समाधानों का उपयोग करके मॉडल किया गया है। इसमें ओरिएंटेशन-निर्भर ड्रैग और शियर-प्रेरित लिफ्ट शामिल है। हाइड्रोडायनामिक टॉर्क जेफ़री ऑर्बिट्स को पुनरुत्पादित करता है।
संपर्क और लुब्रिकेशन (Contact and Lubrication):
- संपर्क: हुकियन (Hookean) के रूप में मॉडल किया गया है जिसमें कूलम्बिक घर्षण (Coulombic friction) शामिल है, जिसमें सामान्य (normal) और स्पर्शरेखीय (tangential) स्प्रिंग-डैम्पर घटक शामिल हैं।
- लुब्रिकेशन: गोलाकार बेलनाकार कणों के बीच लुब्रिकेशन बलों की गणना करने के लिए एक नवीन योजना विकसित की गई थी। एंड-शाफ्ट जंक्शन पर वक्रता (curvature) की विच्छिन्नता (discontinuity) को संबोधित करने के लिए, मॉडल तीन अंतःक्रिया घटकों (शाफ्ट-शाफ्ट, एंड-शाफ्ट, और एंड-एंड) को एक साथ गणना करता है। इन्हें कणों की सापेक्ष स्थिति के आधार पर भारित (weighted) किया जाता है ताकि अंतःक्रिया प्रकारों के बीच सुचारू संक्रमण सुनिश्चित हो सके, जिससे पिछले मॉडलों (जैसे, यामाने एट अल., बटलर और शाफकेह) में मौजूद विच्छिन्नता को हल किया जा सके।
गतिशील टाइमस्टेप (Dynamic Timestep): एक प्रमुख पद्धतिगत नवाचार एक गतिशील टाइमस्टेप गणना है। टाइमस्टेप को संपर्क अंतःक्रियाओं, लुब्रिकेशन बलों और लुब्रिकेशन टॉर्क के बीच सबसे छोटे विशिष्ट काल (timescale) द्वारा निर्धारित किया जाता है। यह संख्यात्मक स्थिरता सुनिश्चित करता है जबकि अनावश्यक रूप से छोटे टाइमस्टेप से बचता है, जिससे गणनात्मक दक्षता में सुधार होता है।
सिमुलेशन सेटअप: सिस्टम फ्लो और ग्रेडिएंट दिशाओं में आवधिक सीमाओं (periodic boundaries) और शियर के लिए लीज़-एडवर्ड्स (Lees-Edwards) सीमाओं का उपयोग करता है। क्रिस्टलीकरण को रोकने के लिए त्रिज्याओं (1:1.4) का एक बिडिस्पर्स मिश्रण उपयोग किया गया है। सिमुलेशन विभिन्न पहलू अनुपातों (aspect ratios, $A$ ) $1 $से$ 20 $तक और आयतन अंशों ($ \phi$) को विरल (dilute) से जैमिंग के करीब तक कवर करते हैं।

मुख्य परिणाम

श्यानता व्यवहार (Viscosity Behavior): मॉडल आयतन अंश ( $\phi$ $ϕ$ ) और पहलू अनुपात ( $A$ $A$ ) दोनों के साथ स्थिर-अवस्था श्यानता ( $\eta_s$ $η_{s}$ ) में व्यवस्थित वृद्धि की भविष्यवाणी करता है।
- क्षणिक प्रतिक्रिया (Transient Response): शियर स्टार्ट-अप पर, एक श्यानता स्पाइक (viscosity spike) देखा जाता है, जिसका कारण कणों के प्रारंभिक यादृच्छिक वितरण को माना जाता है। यह स्पाइक जैसे ही कण संरेखित होते हैं, कम हो जाता है, जिससे एक स्थिर अवस्था प्राप्त होती है। स्थिर अवस्था की तुलना में इस स्पाइक का परिमाण $\phi$ के साथ बढ़ता है।
- स्थिर अवस्था (Steady State): स्थिर-अवस्था श्यानता एक पावर-लॉ डाइवर्जेंस $\eta_s \sim (1 - \phi/\phi_m)^{-\varsigma}$ का अनुसरण करती है। जैमिंग आयतन अंश $\phi_m$ , $A$ बढ़ने के साथ घटता है ( $A=5$ के लिए $\approx 0.57$ से उच्च $A$ के लिए कम मान तक)। मॉडल स्थिर-अवस्था श्यानता की दर-स्वतंत्रता (rate-independence) की पुष्टि करता है, जो ग्रैनुलर सस्पेंशन के ढांचे के अनुरूप है।
- बल योगदान (Force Contributions): निम्न $\phi$ पर, श्यानता फ्लूइड स्ट्रेस और लुब्रिकेशन बलों द्वारा संचालित होती है। जैसे-जैसे $\phi$ बढ़ता है, संपर्क बल (सामान्य और स्पर्शरेखीय दोनों) प्रभावी हो जाते हैं, जो $\phi \approx 0.5$ के आसपास लुब्रिकेशन योगदान से अधिक हो जाते हैं।
सूक्ष्म संरचना और संरेखण (Microstructure and Alignment):
- ऑर्डर पैरामीटर: ओरिएंटेशनल ऑर्डर पैरामीटर ( $S$ ), $\phi$ के साथ एक महत्वपूर्ण बिंदु तक बढ़ता है, जिसके बाद पूर्ण संरेखण को रोकने वाली ज्यामितीय बाधाओं (geometric frustration) के कारण थोड़ा कम हो जाता है।
  - संरेखण दिशा (Alignment Direction):** कण अधिमानतः प्रवाह दिशा ( $x$ -अक्ष) के साथ संरेखित होते हैं, जिसमें प्रवाह-ग्रेडिएंट तल के सापेक्ष एक छोटा झुकाव कोण ( $\phi_v$ ) होता है। यह झुकाव कोण $\phi$ के साथ एक दहलीज तक बढ़ता है और फिर घट जाता है।

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

लेखक कहते हैं कि उनका कार्य (जो पहले से ही गोलों के लिए स्थापित है) को गोलाकार कणों के व्यापक वर्ग तक विस्तृत करता है। यह मॉडल जैमिंग पॉइंट की सूक्ष्म संरचनात्मक विवरणों पर निर्भरता को मैप करने, नॉर्मल स्ट्रेस की जांच करने और विषम स्थितियों (inhomogeneous conditions) के तहत व्यवहार का पता लगाने के लिए एक मंच प्रदान करता है। यह पेपर अपने योगदान को इन सामग्रियों के पूर्ण फ्लुइड मैकेनिकल विवरण की ओर एक कदम के रूप में विनम्रतापूर्वक प्रस्तुत करता है, जो संरेखण और घर्षण के संबंध में जैमिंग पॉइंट के व्यापक मात्रात्मक विवरण को भविष्य के कार्यों के लिए छोड़ देता है। यह सीधे तौर पर किसी विशिष्ट औद्योगिक प्रसंस्करण समस्याओं को हल करने का दावा नहीं करता है, बल्कि इन प्रणालियों के अंतर्निहित भौतिकी की जांच करने के लिए बुनियादी सिमुलेशन इंफ्रास्ट्रक्चर प्रदान करता है।

इसी तरह के और पेपर