Rheology of dense vibrated granular flows: non-monotonic response controlled by granular temperature

यहाँ शोध पत्र "Rheology of dense vibrated granular flows: non-monotonic response controlled by granular temperature" का विस्तृत तकनीकी सारांश दिया गया है।

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

ग्रैनुलर सामग्रियाँ जटिल रियोलॉजिकल व्यवहार प्रदर्शित करती हैं, जो बाहरी बल के आधार पर ठोस जैसी और तरल जैसी अवस्थाओं के बीच परिवर्तित होती हैं। जबकि गैर-कंपित (non-vibrated) ग्रैनुलर प्रवाह का रियोलॉजी अच्छी तरह से परिभाषित है (जैसे कि $\mu(I)$ रियोलॉजी द्वारा), ऊर्ध्वाधर कंपन (vertical vibration) के तहत घने ग्रैनुलर पदार्थों का व्यवहार कम समझ में आया है, विशेष रूप से कंपन मापदंडों (आयाम $A$, आवृत्ति $f$) और विन्यास दबाव (confining pressure $p$) के बीच परस्पर क्रिया के संबंध में।

पिछले अध्ययनों ने विरोधाभासी या अपूर्ण रुझान पहचाने हैं:

• कुछ कार्य बताते हैं कि घर्षण क्षीणन (friction weakening) $A^2/p$ के अनुपात द्वारा नियंत्रित होता है, जो बल श्रृंखलाओं (force chains) को बाधित करने के लिए एक सीमा का संकेत देता है।
• अन्य अध्ययन कंपन आवृत्ति $f$ पर श्यानता (viscosity) और घर्षण की गैर-एकदिष्ट निर्भरता (non-monotonic dependence) देखते हैं, जहाँ प्रवाह क्षमता एक अधिकतम तक बढ़ती है और फिर उच्च आवृत्तियों पर घट जाती है।
• एक एकीकृत ढांचे का अभाव है जो यह समझा सके कि ऊर्जा इंजेक्शन, अपव्यय (dissipation) और विन्यास (confinement) मिलकर पूर्ण आवृत्ति रेंज में प्रभावी श्यानता ($\eta$) को कैसे निर्धारित करते हैं, विशेष रूप से उस शासन (regime) में जहाँ उच्च आवृत्तियों पर श्यानता फिर से बढ़ जाती है।

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

लेखकों ने एक वास्तविक प्रयोगात्मक सेटअप, जिसे स्ट्रेस-इम्पोज़्ड वेन रियोमीटर (stress-imposed vane rheometer) के रूप में जाना जाता है, को दोहराने के लिए डिस्क्रीट एलिमेंट मेथड (DEM) संख्यात्मक सिमुलेशन का उपयोग किया।

• सिमुलेशन सेटअप:
  • प्रणाली: एक बेलनाकार पात्र में सीमित एक मोनोडिस्पर्सिड ग्रैनुलर सिस्टम ($N=2600$ कण, त्रिज्या $r=2$ mm), जिसमें एक शंक्वाकार आधार और द्रव्यमान $M$ वाला एक गोलाकार ढक्कन है।
  • बल (Forcing): पात्र ऊर्ध्वाधर साइनुसोइडल विस्थापन $\Delta z(t) = A \sin(2\pi f t)$ का अनुभव करता है।
  • शियर माप (Shear Measurement): एक आयताकार वेन (vane) जो कणों में डूबा हुआ है, उसे एक निरंतर टॉर्क $T$ दिया जाता है। परिणामी कोणीय वेग $\Omega$ को प्रभावी श्यानता $\eta = T / |\langle \Omega \rangle|$ की गणना के लिए मापा जाता है।
  • भौतिक मॉडल: कण-कण, कण-वेन और कण-दीवार इंटरैक्शन का वर्णन करने के लिए हर्ट्ज़-मिंडलिन (Hertz-Mindlin) संपर्क मॉडल का उपयोग किया गया था, जो लोचदार और विसरणात्मक बलों (elastic and dissipative forces) को ध्यान में रखता है।
• परिवर्तित किए गए पैरामीटर:
  • कंपन आयाम ($A$): $6$से$58$ $\mu$m।
  • कंपन आवृत्ति ($f$): पुनर्गठित त्वरण $\Gamma = (2\pi f)^2 A / g$ को बदलने के लिए $f$ को परिवर्तित किया गया।
  • विन्यास दबाव ($p$): ढक्कन के द्रव्यमान को बदलकर $p$ को परिवर्तित किया गया।
• प्रमुख मेट्रिक्स:
  • प्रभावी श्यानता ($\eta$)।
  • ग्रैनुलर तापमान ($K$): कणों की औसत गतिज ऊर्जा के रूप में परिभाषित ($K \propto \langle v^2 \rangle$)।
  • ऊर्जा अनुपात ($\bar{K}$): कण उत्तेजना ऊर्जा और विन्यास द्वारा निर्धारित ऊर्जा पैमाने के बीच का अनुपात।

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

यह शोध पत्र एक एकीकृत भौतिक ढांचा प्रदान करता है जो पिछले अवलोकनों—श्यानता क्षीणन और सुदृढ़ीकरण—को सुसंगत बनाता है। मुख्य योगदान हैं:

1. शासी चर (Governing Variable) की पहचान: लेखक प्रदर्शित करते हैं कि प्रभावी श्यानता सीधे $A$, $f$, या $p$ द्वारा नियंत्रित नहीं होती है, बल्कि ग्रैनुलर तापमान ($K$) और विन्यास ऊर्जा पैमाने ($K_p$) के अनुपात द्वारा नियंत्रित होती है।
   • इसे $\bar{K} = K / K_p$ के रूप में परिभाषित किया गया है, जहाँ $K_p \propto p$।
2. गैर-एकदिष्टता (Non-Monotonicity) की व्याख्या: शोध पत्र स्पष्ट करता है कि आवृत्ति के संबंध में श्यानता की गैर-एकदिष्ट व्यवहार ऊर्जा हस्तांतरण की दक्षता द्वारा संचालित है।
   • निम्न आवृत्तियों पर, $f$ बढ़ाने से ऊर्जा इंजेक्शन बढ़ता है, जिससे $K$ बढ़ता है और $\eta$ घटता है (द्रवीकरण/fluidization)।
   • उच्च आवृत्तियों पर, ऊर्जा अवशोषण की तुलना में अपव्यय दर (टकराव दर) तेजी से बढ़ती है, जिससे $K$ गिर जाता है, जो बदले में $\eta$ को बढ़ाता है (सघनता/thickening)।
3. स्केलिंग कानूनों का सामंजस्य: अध्ययन पिछले घर्षण अध्ययनों में देखे गए $A^2/p$ स्केलिंग को मान्य करता है, लेकिन इसकी वैधता को केवल न्यूनतम श्यानता शासन (द्रवित अवस्था) तक सीमित करता है। यह दिखाता है कि यह स्केलिंग उच्च आवृत्तियों पर व्यवहार की भविष्यवाणी करने में विफल रहती है जहाँ गैर-एकदिष्ट रुझान हावी होता है।
4. न्यूनतम मॉडल (Minimal Model): एक सैद्धांतिक दो-ब्लॉक मॉडल (परिशिष्ट C) विकसित किया गया है जो $A^2/p$ स्केलिंग को पुनरुत्पादित करता है, यह दर्शाता है कि यह केवल बल श्रृंखला टूटने के बजाय गैर-रेखीय लोचदार इंटरैक्शन (Hertzian contacts) से उत्पन्न होता है।

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

• श्यानता बनाम आवृत्ति ($\eta$ बनाम $f$):

  • प्रभावी श्यानता एक गैर-एकदिष्ट रुझान प्रदर्शित करती है: आवृत्ति बढ़ने पर यह एक न्यूनतम ($\eta_{min}$) तक घटती है, और फिर उच्च आवृत्तियों पर बढ़ती है।
  • यह रुझान ग्रैनुलर तापमान ($K$) के साथ पूरी तरह से सहसंबंधित है, जो एक मध्यवर्ती आवृत्ति $f^*$ पर चरम (peak) पर पहुँचता है।
  • इस संबंध को एक पावर लॉ द्वारा वर्णित किया गया है: $\eta \propto K^\alpha$ जहाँ $\alpha \in [-2, -1.5]$।

• श्यानता बनाम विन्यास दबाव ($\eta$ बनाम $p$):

  • विन्यास दबाव $p$ बढ़ाने से प्रभावी श्यानता बढ़ती है (प्रणाली को स्थिर करती है)।
  • हालाँकि, $p$ का ग्रैनुलर तापमान $K$ पर नगण्य प्रभाव पड़ता है। $p$ का प्रभाव शुद्ध रूप से कणों को हिलाने के लिए आवश्यक ऊर्जा अवरोध (energy barrier) के माध्यम से है।

• मास्टर कर्व ($\eta$ बनाम $\bar{K}$):

  • जब $\bar{K} = K / (p \pi r^3)$ के विरुद्ध $\eta$ को प्लॉट किया जाता है, तो विभिन्न आयामों, आवृत्तियों और दबावों के लिए सभी डेटा बिंदु एक एकल मास्टर कर्व पर आ जाते हैं।
  • यह पुष्टि करता है कि रियोलॉजिकल अवस्था उत्तेजना ऊर्जा (कणों के पिंजरों को तोड़ने की कोशिश) और विन्यास ऊर्जा (पिंजरों को बरकरार रखने की कोशिश) के बीच प्रतिस्पर्धा द्वारा निर्धारित होती है।

• $A^2/p$ की सीमाएँ:

  • जबकि $\eta_{min}$, $A^2/p$ के साथ स्केल करता है, यह चर उच्च आवृत्तियों (बढ़ते हुए भाग) पर श्यानता की भविष्यवाणी नहीं कर सकता है।
  • $A^2/p$ स्केलिंग एक विशिष्ट द्रवित अवस्था तक पहुँचने का परिणाम है, लेकिन यह आवृत्ति-निर्भर ऊर्जा अपव्यय तंत्र को कैप्चर नहीं करती है।

5. महत्व और निहितार्थ (Significance and Implications)

• एकीकृत ढांचा: यह शोध पत्र $A^2/p$ (घर्षण क्षीणन) पर केंद्रित अध्ययनों और गैर-एकदिष्ट आवृत्ति निर्भरता को देखने वाले अध्ययनों के बीच के स्पष्ट विरोधाभास को हल करता है। यह स्थापित करता है कि ऊर्जा इंजेक्शन बनाम अपव्यय ही मौलिक चालक है।
• नियंत्रण पैरामीटर के रूप में ग्रैनुलर तापमान: यह कार्य "ग्रैनुलर तापमान" को वाइब्रेटेड सिस्टम के रियोलॉजी के लिए एक वर्णनात्मक मीट्रिक से ऊपर उठाकर एक भविष्य कहनेवाला नियंत्रण चर (predictive control variable) के रूप में स्थापित करता है, जो आणविक प्रणालियों में तापीय तापमान के समान है, हालांकि इसके गैर-संतुलन गतिकी (non-equilibrium dynamics) भिन्न हैं।
• व्यावहारिक अनुप्रयोग: गैर-एकदिष्ट प्रतिक्रिया को समझना औद्योगिक अनुप्रयोगों (जैसे हॉपर फ्लो, फार्मास्युटिकल मिक्सिंग) और भूभौतिकीय घटनाओं (जैसे भूकंप-प्रेरित लैंडस्लाइड) के लिए महत्वपूर्ण है। यह सुझाव देता है कि केवल कंपन आवृत्ति बढ़ाने से हमेशा प्रवाह क्षमता में सुधार नहीं होता है; एक इष्टतम आवृत्ति रेंज होती है, जिसके बाद सामग्री प्रवाह के प्रति अधिक प्रतिरोधी हो जाती है।
• सैद्धांतिक अंतर्दृष्टि: निष्कर्षों से पता चलता है कि गैर-संतुलन ग्रैनुलर प्रणालियों में, ड्राइविंग मापदंडों और सिस्टम प्रतिक्रिया के बीच का संबंध जटिल ऊर्जा अपव्यय तंत्रों द्वारा मध्यस्थ होता है, जो सरल स्केलिंग कानूनों को सभी शासनों में लागू होने से रोकता है।

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