2025-11-22T10:49:15.885853

Transient rod-climbing in an Oldroyd-B fluid

Ruangkriengsin, Brandão, Wu et al.
The Weissenberg effect, or rod-climbing phenomenon, occurs in non-Newtonian fluids where the fluid interface ascends along a rotating rod. Despite its prominence, theoretical insights into this phenomenon remain limited. In earlier work, Joseph \& Fosdick (\emph{Arch. Rat. Mech. Anal.}, vol. 49, 1973, pp. 321--380) employed domain perturbation methods for second-order fluids to determine the equilibrium interface height by expanding solutions based on the rotation speed. In this work, we investigate the time-dependent interface height through asymptotic analysis with dimensionless variables and equations using the Oldroyd-B model. We begin by neglecting surface tension and inertia to focus on the interaction between gravity and viscoelasticity. In the small-deformation scenario, the governing equations indicate the presence of a boundary layer in time, where the interface rises rapidly over a short time scale before gradually approaching a steady state. By employing a stretched time variable, we derive the transient velocity field and corresponding interface profile on this short time scale and recover the steady-state profile on a longer time scale. Subsequently, we reintroduce small but finite inertial effects to investigate their interplay with viscoelasticity and propose a criterion for determining the conditions under which rod-climbing occurs.
academic

Oldroyd-B द्रव में क्षणिक रॉड-क्लाइम्बिंग

मूल जानकारी

  • पेपर ID: 2501.01016
  • शीर्षक: Oldroyd-B द्रव में क्षणिक रॉड-क्लाइम्बिंग
  • लेखक: Tachin Ruangkriengsin, Rodolfo Brandão, Katie Wu, Jonghyun Hwang, Evgeniy Boyko, Howard A. Stone
  • वर्गीकरण: physics.flu-dyn (द्रव गतिविज्ञान)
  • प्रकाशन तिथि: 2 जनवरी 2025 (arXiv को प्रस्तुत)
  • पत्रिका: जर्नल ऑफ फ्लूइड मैकेनिक्स (समीक्षाधीन)
  • पेपर लिंक: https://arxiv.org/abs/2501.01016

सारांश

यह पेपर गैर-न्यूटोनियन द्रव में Weissenberg प्रभाव (रॉड-क्लाइम्बिंग घटना) का अध्ययन करता है, जहाँ द्रव इंटरफेस घूर्णनशील रॉड के साथ ऊपर उठता है। लेखकों ने Oldroyd-B मॉडल का उपयोग करके, स्पर्शोन्मुख विश्लेषण और विमाहीन चर के माध्यम से, समय-निर्भर इंटरफेस ऊंचाई का अध्ययन किया। अनुसंधान से पता चलता है कि छोटी विकृति के मामले में, नियंत्रण समीकरण समय सीमा परत के अस्तित्व को दर्शाते हैं, जहाँ इंटरफेस छोटे समय पैमाने पर तेजी से ऊपर उठता है, फिर धीरे-धीरे स्थिर अवस्था के करीब पहुंचता है। खिंचाव समय चर को प्रस्तुत करके, लेखकों ने क्षणिक वेग क्षेत्र और संबंधित इंटरफेस समोच्च प्राप्त किए, और छोटे लेकिन परिमित जड़त्व प्रभावों को पुनः प्रस्तुत किया ताकि उनकी viscoelasticity के साथ परस्पर क्रिया का अध्ययन किया जा सके।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

रॉड-क्लाइम्बिंग घटना (Weissenberg प्रभाव) गैर-न्यूटोनियन द्रव के सबसे विशेषता वाले व्यवहारों में से एक है। जब रॉड viscoelastic द्रव में घूमता है, तो द्रव की सतह रॉड के साथ ऊपर की ओर चढ़ती है। यह घटना सैद्धांतिक रूप से अभी भी गहन समझ की कमी है।

महत्व

  1. सैद्धांतिक महत्व: यह घटना जटिल द्रव व्यवहार को समझने की एक मौलिक अवधारणा है, जिसका शिक्षण और अनुसंधान में महत्वपूर्ण स्थान है
  2. व्यावहारिक मूल्य: रियोलॉजी माप में महत्वपूर्ण अनुप्रयोग है, जिसका उपयोग सामग्री पैरामीटर जैसे सामान्य तनाव गुणांक और विश्राम समय निर्धारित करने के लिए किया जा सकता है
  3. इंजीनियरिंग अनुप्रयोग: बहुलक समाधान और अन्य जटिल द्रव के औद्योगिक अनुप्रयोग के लिए मार्गदर्शन प्रदान करता है

मौजूदा विधियों की सीमाएं

  1. सैद्धांतिक सीमा: पूर्ववर्ती अनुसंधान मुख्य रूप से स्थिर अवस्था के व्यवहार पर केंद्रित है, द्वितीय-क्रम द्रव मॉडल का उपयोग करते हुए
  2. मॉडल समस्या: द्वितीय-क्रम द्रव मॉडल केवल छोटे Deborah संख्या या Weissenberg संख्या के तहत वैध है, और अत्यधिक अस्थिर है
  3. समय निर्भरता: क्षणिक गतिविज्ञान के सैद्धांतिक विश्लेषण की कमी है

अनुसंधान प्रेरणा

लेखकों ने द्वितीय-क्रम द्रव मॉडल की गैर-स्थिर प्रवाह समस्याओं से बचने के लिए Oldroyd-B मॉडल चुना, और पहली बार रॉड-क्लाइम्बिंग घटना की समय-निर्भरता का व्यवस्थित रूप से अध्ययन किया।

मुख्य योगदान

  1. पहला सैद्धांतिक विश्लेषण: Oldroyd-B द्रव में रॉड-क्लाइम्बिंग घटना का पहला संपूर्ण क्षणिक सैद्धांतिक विश्लेषण प्रदान करता है
  2. समय सीमा परत सिद्धांत: समय सीमा परत संरचना की खोज और विश्लेषण, तेजी से ऊपर उठने और धीरे-धीरे स्थिर अवस्था की ओर जाने की दो-चरणीय प्रक्रिया को प्रकट करता है
  3. विमाहीन विश्लेषण ढांचा: Weissenberg संख्या के आधार पर एक व्यवस्थित विमाहीन विश्लेषण विधि स्थापित करता है
  4. जड़त्व प्रभाव अनुसंधान: पहली बार छोटे लेकिन परिमित जड़त्व प्रभाव और viscoelasticity की परस्पर क्रिया का अध्ययन करता है
  5. क्लाइम्बिंग मानदंड: रॉड-क्लाइम्बिंग की घटना की शर्तें निर्धारित करने के लिए एक सैद्धांतिक मानदंड प्रस्तावित करता है: Re<4βpWiRe < 4\beta_p Wi

विधि विवरण

कार्य परिभाषा

अनंत लंबे वृत्ताकार रॉड (त्रिज्या aa) के कोणीय वेग Ω\Omega पर viscoelastic बहुलक समाधान में घूमने पर, द्रव मुक्त सतह की ऊंचाई h(r,t)h(r,t) के समय विकास नियम का अध्ययन करना।

गणितीय मॉडल

नियंत्रण समीकरण

Oldroyd-B संवैधानिक मॉडल का उपयोग:

  • विन्यास टेंसर विकास समीकरण: At+uAAu(u)TA=1λ(AI)\frac{\partial\mathbf{A}}{\partial t} + \mathbf{u} \cdot \nabla\mathbf{A} - \mathbf{A} \cdot \nabla\mathbf{u} - (\nabla\mathbf{u})^T \cdot \mathbf{A} = -\frac{1}{\lambda}(\mathbf{A} - \mathbf{I})
  • तनाव टेंसर: σ=pI+2μsE+G(AI)\boldsymbol{\sigma} = -p\mathbf{I} + 2\mu_s\mathbf{E} + G(\mathbf{A} - \mathbf{I})
  • गति समीकरण: σ=ρgez\nabla \cdot \boldsymbol{\sigma} = \rho g\mathbf{e}_z

विमाहीनकरण

मुख्य विमाहीन पैरामीटर प्रस्तुत करना:

  • Weissenberg संख्या: Wi=λΩWi = \lambda\Omega
  • गुरुत्व पैरामीटर: G=ρga(μs+μp)ΩG = \frac{\rho ga}{(\mu_s + \mu_p)\Omega}
  • Reynolds संख्या: Re=ρa2Ωμs+μpRe = \frac{\rho a^2\Omega}{\mu_s + \mu_p}

विशेषता सीमा

लेखकों ने मुख्य विशेषता सीमा की पहचान की: G=WiGG = Wi \cdot \mathcal{G}, जहाँ G=O(1)\mathcal{G} = O(1), यह द्वि-पैरामीटर स्पर्शोन्मुख विस्तार को एकल-पैरामीटर विस्तार में सरल करता है।

तकनीकी नवाचार बिंदु

1. समय पैमाने पृथक्करण

दो समय पैमानों की खोज:

  • बाहरी परत समय: T=O(1)T = O(1), स्थिर अवस्था के व्यवहार के अनुरूप
  • आंतरिक परत समय: τ=T/Wi\tau = T/Wi, क्षणिक सीमा परत के अनुरूप

2. डोमेन विक्षोभ विधि

छोटी इंटरफेस विकृति के लिए, सीमा शर्तों को Z=H(R,T)Z = H(R,T) से Z=0Z = 0 तक Taylor विस्तारित करना, मुक्त सतह सीमा शर्तों को सरल करता है।

3. Weber रूपांतरण

जड़त्व प्रभाव को संभालते समय, अर्ध-अनंत डोमेन पर Bessel अवकल संचालक को संभालने के लिए Weber रूपांतरण का उपयोग करना, आंशिक अवकल समीकरण को सामान्य अवकल समीकरण में परिवर्तित करता है।

प्रायोगिक सेटअप

सैद्धांतिक विश्लेषण ढांचा

यह पेपर शुद्ध सैद्धांतिक अनुसंधान है, स्पर्शोन्मुख विश्लेषण विधि का उपयोग करता है:

  1. छोटा पैरामीटर: Wi1Wi \ll 1 (कमजोर viscoelasticity)
  2. गुरुत्व श्रेणी: GWi1G \gg Wi^{-1} (मजबूत गुरुत्व प्रभाव)
  3. छोटी विकृति धारणा: इंटरफेस विकृति अपेक्षाकृत छोटी है

विश्लेषण रणनीति

  1. अनुभाग 4: बाहरी परत विश्लेषण - स्थिर अवस्था इंटरफेस समोच्च
  2. अनुभाग 5: आंतरिक परत विश्लेषण - क्षणिक इंटरफेस समोच्च
  3. अनुभाग 6: छोटे जड़त्व प्रभाव सहित विश्लेषण

प्रायोगिक परिणाम

मुख्य परिणाम

1. स्थिर अवस्था इंटरफेस समोच्च

Hsteady(R)=2βpWi2GR4H_{steady}(R) = \frac{2\beta_p Wi^2}{\mathcal{G}R^4}

यह Joseph & Fosdick (1973) द्वारा द्वितीय-क्रम द्रव मॉडल का उपयोग करके प्राप्त परिणाम के अनुरूप है।

2. क्षणिक इंटरफेस समोच्च

Htransient(R,τ)=2βp(1(τ+1)eτ)Wi2GR4H_{transient}(R,\tau) = \frac{2\beta_p(1-(\tau+1)e^{-\tau})Wi^2}{\mathcal{G}R^4}

जहाँ τ=T/Wi\tau = T/Wi खिंचाव समय चर है।

3. क्लाइम्बिंग मानदंड

जब जड़त्व प्रभाव को शामिल किया जाता है, तो रॉड-क्लाइम्बिंग की शर्त है: 4Elβp>14El\beta_p > 1 जहाँ El=Wi/ReEl = Wi/Re लोच संख्या है।

समय विकास विशेषताएं

  1. तेजी से ऊपर उठने का चरण: इंटरफेस τ6\tau \approx 6 समय में तेजी से ऊपर उठता है
  2. स्थिर अवस्था के करीब: क्षणिक समाधान लंबे समय की सीमा में स्थिर अवस्था समाधान में परिवर्तित होता है
  3. गैर-एकरस व्यवहार: जड़त्व को शामिल करते समय, इंटरफेस ऊंचाई पहले ऊपर जा सकती है फिर नीचे आ सकती है

मुख्य निष्कर्ष

  1. परिधीय तनाव: बहुलक खिंचाव द्वारा उत्पन्न परिधीय तनाव क्लाइम्बिंग को चलाने का मुख्य तंत्र है
  2. जड़त्व प्रतिस्पर्धा: जड़त्व द्वारा उत्पन्न अपकेंद्र बल viscoelasticity के परिधीय तनाव के साथ प्रतिस्पर्धा करता है
  3. सीमा परत संरचना: समय सीमा परत का अस्तित्व तेजी से क्षणिक और धीमी स्थिर अवस्था के दो-चरणीय व्यवहार को समझाता है

संबंधित कार्य

ऐतिहासिक विकास

  1. Serrin (1959): पहला सैद्धांतिक प्रयास, द्वितीय-क्रम द्रव मॉडल का उपयोग करते हुए
  2. Giesekus (1961): जड़त्व प्रभाव को अनदेखा करते हुए व्यापक संवैधानिक समीकरण
  3. Joseph & Fosdick (1973): व्यवस्थित डोमेन विक्षोभ विधि ढांचा
  4. Oldroyd (1950, 1951): प्रारंभिक Oldroyd-B मॉडल अनुप्रयोग

इस पेपर के लाभ

  1. समय निर्भरता: पहला संपूर्ण क्षणिक विश्लेषण
  2. मॉडल स्थिरता: द्वितीय-क्रम द्रव मॉडल की अस्थिरता समस्या से बचाता है
  3. विमाहीन ढांचा: व्यवस्थित विमाहीन पैरामीटर विश्लेषण
  4. जड़त्व प्रभाव: पहली बार जड़त्व और viscoelasticity की परस्पर क्रिया का अध्ययन

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. सीमा परत संरचना: रॉड-क्लाइम्बिंग घटना में स्पष्ट समय सीमा परत होती है, इंटरफेस तेजी से ऊपर उठता है फिर धीरे-धीरे स्थिर अवस्था की ओर जाता है
  2. पैमाने पृथक्करण: क्षणिक समय पैमाना O(Wi)O(Wi) है, स्थिर अवस्था समय पैमाना O(1)O(1) है
  3. क्लाइम्बिंग तंत्र: बहुलक खिंचाव द्वारा उत्पन्न परिधीय तनाव क्लाइम्बिंग को चलाने का मुख्य तंत्र है
  4. जड़त्व प्रभाव: छोटा जड़त्व प्रभाव क्लाइम्बिंग व्यवहार को महत्वपूर्ण रूप से प्रभावित कर सकता है, एक महत्वपूर्ण लोच संख्या मौजूद है

सीमाएं

  1. छोटी विकृति धारणा: केवल छोटी इंटरफेस विकृति के मामले में लागू
  2. कमजोर viscoelasticity: Wi1Wi \ll 1 श्रेणी में सीमित
  3. मजबूत गुरुत्व: GWi1G \gg Wi^{-1} शर्त की आवश्यकता है
  4. सरलीकृत मॉडल: सतह तनाव और संपर्क कोण प्रभाव को अनदेखा करता है

भविष्य की दिशाएं

  1. उच्च Weissenberg संख्या: Wi1Wi \gg 1 स्थिति में व्यवहार का अध्ययन
  2. जटिल संवैधानिक मॉडल: FENE-P, Giesekus आदि मॉडल तक विस्तार
  3. परिमित विकृति: बड़ी इंटरफेस विकृति के मामले को संभालना
  4. प्रायोगिक सत्यापन: प्रायोगिक डेटा के साथ तुलना सत्यापन

गहन मूल्यांकन

लाभ

  1. सैद्धांतिक कठोरता: गणितीय व्युत्पत्ति सुदृढ़ है, स्पर्शोन्मुख विश्लेषण विधि का उचित अनुप्रयोग
  2. भौतिक अंतर्दृष्टि: समय सीमा परत संरचना और क्लाइम्बिंग तंत्र को प्रकट करता है
  3. विधि नवाचार: Weber रूपांतरण और डोमेन विक्षोभ विधि का चतुर संयोजन
  4. व्यावहारिक मूल्य: रियोलॉजी माप के लिए सैद्धांतिक आधार प्रदान करता है

कमियां

  1. पैरामीटर श्रेणी: लागू पैरामीटर श्रेणी अपेक्षाकृत सीमित है
  2. प्रायोगिक तुलना: प्रायोगिक डेटा के साथ सीधी तुलना की कमी है
  3. संख्यात्मक सत्यापन: कुछ परिणाम केवल संख्यात्मक गणना दिए गए हैं, विश्लेषणात्मक अभिव्यक्ति की कमी है

प्रभाव

  1. शैक्षणिक योगदान: क्षणिक रॉड-क्लाइम्बिंग सिद्धांत में अंतराल को भरता है
  2. विधि मूल्य: प्रदान की गई विश्लेषण विधि अन्य जटिल द्रव समस्याओं तक विस्तारित की जा सकती है
  3. इंजीनियरिंग अनुप्रयोग: रियोलॉजी यंत्र डिजाइन और सामग्री लक्षण वर्णन के लिए सैद्धांतिक समर्थन प्रदान करता है

लागू परिदृश्य

  1. पतला बहुलक समाधान: Boger द्रव जैसे स्थिर कतरनी चिपचिपापन viscoelastic द्रव
  2. कम घूर्णन गति: छोटी Weissenberg संख्या के प्रायोगिक शर्तें
  3. रियोलॉजी माप: सामान्य तनाव गुणांक और विश्राम समय का निर्धारण

संदर्भ

  1. Joseph, D.D. & Fosdick, R.L. (1973). विभिन्न गति से घूमने वाले सिलेंडरों के बीच द्रव की मुक्त सतह। Arch. Rat. Mech. Anal., 49, 321-380.
  2. Oldroyd, J.G. (1950). रियोलॉजिकल अवस्था समीकरणों के निर्माण पर। Proc. R. Soc. A, 200, 523-541.
  3. Bird, R.B., Armstrong, R.C. & Hassager, O. (1987). बहुलक द्रव की गतिविज्ञान, खंड 1, 2nd संस्करण। John Wiley and Sons.

समग्र मूल्यांकन: यह एक उच्च गुणवत्ता का सैद्धांतिक द्रव गतिविज्ञान पेपर है, जो पहली बार Oldroyd-B द्रव में रॉड-क्लाइम्बिंग घटना के क्षणिक व्यवहार का व्यवस्थित विश्लेषण करता है। पेपर गणितीय व्युत्पत्ति में कठोर है, भौतिक अंतर्दृष्टि गहन है, विधि नवीन है, और जटिल द्रव सिद्धांत और रियोलॉजी अनुप्रयोग के लिए महत्वपूर्ण योगदान देता है। यद्यपि पैरामीटर श्रेणी सीमाएं हैं, इसका सैद्धांतिक ढांचा और विश्लेषण विधि बहुत अच्छी विस्तार क्षमता रखते हैं।