🔢 mathematics

Modelling Capillary Rise with a Slip Boundary Condition: Well-posedness and Long-time Dynamics of Solutions to Washburn's Equation

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

मूल लेखक: Isidora Rapajić, Srboljub Simić, Endre Süli

प्रकाशित 2026-04-10

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

CC BY 4.0

मूल लेखक: Isidora Rapajić, Srboljub Simić, Endre Süli

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

कल्पना कीजिए कि आपके पास एक बहुत ही पतली स्ट्रॉ (नली) है जो पानी के गिलास में डूबी हुई है। आप जानते हैं कि पानी वहां सिर्फ बैठा नहीं रहता; वह अपने आप स्ट्रॉ में ऊपर चढ़ता है, जो कुछ समय के लिए गुरुत्वाकर्षण को चुनौती देता है। इसे केशिका उत्थान (capillary rise) कहा जाता है, और यही वह सिद्धांत है जो पेड़ों को अपनी जड़ों से पत्तियों तक पानी खींचने में मदद करता है या यह कि कैसे एक पेपर टॉवल किसी फैले हुए तरल को सोख लेता है।

एक सदी से भी अधिक समय से, वैज्ञानिक वॉशबर्न के समीकरण (Washburn's Equation) नामक एक प्रसिद्ध सूत्र का उपयोग यह अनुमान लगाने के लिए करते आए हैं कि तरल कितनी तेजी से और कितनी ऊंचाई तक चढ़ेगा। हालाँकि, इस पुराने सूत्र में कुछ समस्याएँ थीं:

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

रपाज़िक (Rapajić), सिमिक (Simić) और सुली (Süli) का यह शोध पत्र इन समस्याओं को ठीक करता है। यहाँ इसका एक सरल विवरण दिया गया है, जिसमें रोजमर्रा के उदाहरणों का उपयोग किया गया है।

1. "फिसलन भरी" दीवार (नया भौतिक विज्ञान)

पुराना मॉडल मानता था कि पाइप के किनारे पर मौजूद पानी के अणु अपनी जगह पर चिपके रहते हैं, जैसे ब्रेक लगे हुए कार की स्थिति हो। लेकिन वास्तव में, विशेष रूप से छोटी नलियों में, तरल दीवार के साथ थोड़ा "फिसल" सकता है, जैसे बर्फ के टुकड़े पर चलती कार।

लेखकों ने समीकरण में एक "स्लिप पैरामीटर" (slip parameter) जोड़ा है। इसे एक "फिसलन वाले डायल" की तरह समझें।

डायल 0 पर: दीवार चिपचिपी है (कोई फिसलन नहीं)।
डायल ऊपर की ओर: दीवार फिसलन भरी है।
परिणाम: जब दीवार फिसलन भरी होती है, तो तरल तेजी से ऊपर चढ़ता है क्योंकि उसे पीछे रोकने के लिए घर्षण कम होता है। लेखकों ने सिद्ध किया कि इस नए "फिसलन" वाले कारक के साथ भी, भौतिक विज्ञान अभी भी तर्कसंगत बना रहता है।

2. "जीरो सेकंड" की समस्या को ठीक करना (गणित)

पुराने समीकरण के साथ सबसे बड़ी सिरदर्दी यह थी कि समय = 0 पर क्या होता है। यदि आप तरल को बिल्कुल नीचे (ऊंचाई 0) और वेग 0 से शुरू करते हैं, तो गणित त्वरण (acceleration) की गणना करने की कोशिश करता है और चिल्ला उठता है, "मैं यह नहीं कर सकता!" यह बिल्कुल वैसा ही है जैसे बिना इंजन की शक्ति के, कार के बिल्कुल उसी क्षण की गति की गणना करना जब वह स्थिर अवस्था से चलना शुरू करती है।

लेखकों ने महसूस किया कि पुराना प्रमाण कि यह समीकरण काम करता है, "लीकी" (कमजोर) था। उन्होंने इसे निम्नलिखित तरीकों से ठीक किया:

शुरुआत को सुचारू बनाना: उन्होंने तरल को एक बहुत ही मामूली, गैर-शून्य ऊंचाई (जैसे कि नीचे पहले से ही बैठा एक छोटा सा قط्रा) से शुरू होने की अनुमति दी, न कि एक पूर्ण गणितीय बिंदु से।
"रेगुलराइजेशन" (Regularization) का तरीका: उन्होंने समीकरण में अस्थायी रूप से थोड़ा सा "गणितीय गोंद" जोड़ा ताकि यह व्यवहार कर सके, फिर इसे हल किया, और फिर धीरे-धीरे उस गोंद को हटा दिया। इससे यह सिद्ध हुआ कि एक समाधान मौजूद है और वह अद्वितीय है—अर्थात तरल अनुमानित तरीके से व्यवहार करता है, अराजक रूप से नहीं।

3. फिनिश लाइन की दौड़ (स्थिरता)

एक बार जब तरल ऊपर चढ़ना शुरू कर देता है, तो क्या वह रुक जाता है? या क्या वह लक्ष्य से आगे निकल जाता है, वापस गिरता है, और एक बाउंस बॉल की तरह ऊपर-नीचे होता रहता है?

लेखकों ने पूछा: क्या तरल अंततः एक विशिष्ट ऊंचाई पर आकर रुक जाएगा?

गंतव्य: एक "कंफर्ट ज़ोन" ऊंचाई (जिसे साम्यावस्था ऊंचाई कहा जाता है) होती है जहाँ सतह तनाव (surface tension) का खिंचाव पानी के स्तंभ के वजन को संतुलित करता है।
यात्रा: तरल के "भारीपन" (जड़त्व/inertia) और उसकी "गाढ़ेपन" (श्यानता/viscosity) के आधार पर, तरल:
- एक पत्ते तक पहुँचने वाले घोंघे की तरह सुचारू रूप से चढ़ सकता है।
- एक पेंडुलम की तरह लक्ष्य से आगे निकलकर डगमगा सकता है और फिर स्थिर हो सकता है।
प्रमाण: लेखकों ने एक गणितीय "ऊर्जा जाल" (जिसे लयापुनोव फंक्शन कहा जाता है) बनाया। एक कटोरे के अंदर रखी मार्बल (कंचे) की कल्पना करें। आप मार्बल को कहीं भी रखें (जब तक कि वह कटोरे के अंदर है), वह अंततः बिल्कुल नीचे की ओर ही लुढ़केगी। उन्होंने सिद्ध किया कि किसी भी शुरुआती ऊंचाई के लिए (जब तक कि वह अत्यधिक ऊँची न हो), तरल उस मार्बल की तरह है, और वह अंततः सही ऊंचाई पर आकर ही स्थिर होगा। वह कटोरे से बाहर नहीं निकलेगा।

4. यात्रा पर "फिसलन" का प्रभाव

एक दिलचस्प खोज यह है कि हालांकि "फिसलन भरी दीवार" तरल को तेजी से ऊपर चढ़ाती है, लेकिन यह उसके गंतव्य को नहीं बदलती। तरल अभी भी उसी अंतिम ऊंचाई पर रुकता है, चाहे पाइप कितना भी फिसलन भरा क्यों न हो। यह बस एक अलग शैली के साथ वहां पहुँचता है (शायद तेज़ डगमगाहट या एक सुचारू फिसलन के साथ)।

सारांश: यह क्यों महत्वपूर्ण है

इस शोध पत्र को छोटी नलियों में तरल पदार्थों को समझने के लिए "यूजर मैनुअल अपडेट" की तरह समझें।

पहले: मैनुअल में एक चेतावनी थी कि "यदि आप शून्य से शुरू करते हैं, तो गणित विफल हो जाता है," और "हमें 100% यकीन नहीं है कि तरल रुकेगा।"
अब: नया मैनुअल कहता है, "हमने गणित को ठीक कर दिया है ताकि यह पहले सेकंड से ही काम करे। हमने सिद्ध किया है कि तरल हमेशा सही ऊंचाई पर आकर रुकेगा, चाहे पाइप चिपचिपा हो या फिसलन भरा।"

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

यह शोध पत्र "Modelling Capillary Rise with a Slip Boundary Condition: Well-posedness and Long-time Dynamics of Solutions to Washburn's Equation" (रपाजीक, सिमिक और सुली द्वारा) का एक विस्तृत तकनीकी सारांश है।

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

यह शोध पत्र एक ऊर्ध्वाधर बेलनाकार पाइप में तरल के स्वतः स्फूर्त केशिका उत्थान (spontaneous capillary rise) के गणितीय मॉडलिंग को संबोधित करता है। शास्त्रीय वॉशबर्न समीकरण (1921) इस घटना का वर्णन करता है लेकिन पारंपरिक रूप से पाइप की दीवार पर नो-स्लिप बाउंडरी कंडीशन (no-slip boundary condition) मानता है। यह धारणा हु-स्क्रीवन विरोधाभास (Huh–Scriven paradox) की ओर ले जाती है, जहाँ दीवार पर तरल का वेग शून्य होता है, फिर भी संपर्क रेखा (liquid, solid, और gas के बीच का इंटरफेस) गतिमान होती है।

इसे हल करने के लिए, लेखक एक स्लिप बाउंडरी कंडीशन (जिसे स्लिप लंबाई $L$ या आयामहीन पैरामीटर $\beta$ द्वारा दर्शाया गया है) को शामिल करते हैं। इसके अलावा, यह शोध पत्र वॉशबर्न समीकरण के समाधानों के वैश्विक अस्तित्व (global existence), विशिष्टता (uniqueness) और वेल-पोज़डनेस (well-posedness) के संबंध में पिछले गणितीय प्रमाणों की कमियों को कठोरता से संबोधित करने का लक्ष्य रखता है, विशेष रूप से गैर-शून्य प्रारंभिक डेटा और अनंत समय अंतराल $t \in [0, \infty)$ के लिए।

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

A. भौतिक व्युत्पत्ति (Physical Derivation)

लेखक वॉशबर्न समीकरण को केवल एक कंट्रोल वॉल्यूम पर न्यूटन के द्वितीय नियम को लागू करने के बजाय, प्रथम सिद्धांतों (continuum mechanics) से व्युत्पन्न करते हैं।

मान्यताएँ: असंपीड्य न्यूटोनियन तरल (incompressible Newtonian fluid), स्थिर घनत्व $\rho$ और श्यानता (viscosity) $\mu$ , स्थिर पृष्ठ तनाव $\gamma$ , और संपर्क कोण $\theta$ ।
प्रवाह प्रोफ़ाइल (Flow Profile): वे एक अक्षीय सममित (axially symmetric), समानांतर प्रवाह मानते हैं जो एक स्लिप पैरामीटर द्वारा संशोधित प्वाइज़ुली वेग प्रोफ़ाइल (Poiseuille velocity profile) है।
संवेग संतुलन (Momentum Balance): वे मेनिस्कस ऊंचाई $h(t)$ के लिए गवर्निंग साधारण अवकल समीकरण (ODE) व्युत्पन्न करने के लिए संवेग संतुलन के वैश्विक (integral) रूप का उपयोग करते हैं।
परिणामी समीकरण:
$\frac{d}{dt}\left[\rho h(t) \frac{dh}{dt}\right] + \frac{8\mu\beta}{R^2} h(t) \frac{dh}{dt} + \rho g h(t) = \frac{2\gamma \cos \theta}{R}$
जहाँ $\beta = (1 + 4L/R)^{-1}$ स्लिप पैरामीटर है ( $\beta=1$ के लिए no-slip)।

B. आयामी विश्लेषण और मॉडल न्यूनीकरण (Dimensional Analysis and Model Reduction)

समीकरण को साम्यावस्था ऊंचाई $h_e$ और एक विशिष्ट समय पैमाने $\tau$ का उपयोग करके गैर-आयामी (non-dimensionalized) बनाया जाता है। इससे निम्नलिखित आयामहीन समीकरण प्राप्त होता है:
$\omega(HH')' + \beta HH' + H = 1$
जहाँ $H$ आयामहीन ऊंचाई है, $T$ आयामहीन समय है, और $\omega$ ओहनेसोरज (Ohnesorge) और बॉन्ड नंबरों से संबंधित एक पैरामीटर है।

लेखक समय और ऊंचाई के लिए पावर-लॉ स्केलिंग (power-law scalings) पेश करके एक व्यवस्थित मॉडल न्यूनीकरण विश्लेषण करते हैं।
वे $\omega$ $ω$ (जड़त्व बनाम श्यानता बनाम गुरुत्वाकर्षण) के सापेक्ष परिमाण के आधार पर चार अलग-अलग प्रवाह शासन (flow regimes) की पहचान करते हैं:
1. नगण्य गुरुत्वाकर्षण (Inertia + Viscosity)।
2. नगण्य जड़त्व (Viscosity + Gravity)।
3. नगण्य गुरुत्वाकर्षण और जड़त्व (Viscous dominated, जो शास्त्रीय $\sqrt{t}$ वॉशबर्न नियम को पुनः प्राप्त करता है)।
4. नगण्य श्यानता (Inertia + Gravity, जो दोलन व्यवहार (oscillatory behavior) की ओर ले जाता है)।

C. गणितीय विश्लेषण (वेल-पोज़डनेस)

मुख्य गणितीय योगदान रूपांतरित समीकरण के लिए समाधानों के अस्तित्व और विशिष्टता को सिद्ध करने में शामिल है:
$u'' + \frac{\beta}{\sqrt{\omega}} u' + \sqrt{2u} = 1$
जहाँ $u = \frac{1}{2}H^2$ ।

रूपांतरण (Transformation): प्रतिस्थापन $u = H^2/2$ के माध्यम से अग्रणी पद (leading term) की नॉन-लिनियरिटी को हटा दिया जाता है।
अस्तित्व (Existence): लेखक एक रेगुलराइजेशन दृष्टिकोण का उपयोग करते हैं। वे नॉन-लीनियर टर्म में एक छोटा $\epsilon$ जोड़कर समीकरण के रेगुलराइज्ड संस्करण को हल करते हैं ताकि लिप्सचिट्ज़ (Lipschitz) स्थिति बनी रहे। वे फिर स्थानीय अस्तित्व के लिए पिकार्ड-लिंडेलोफ प्रमेय (Picard–Lindelöf theorem) लागू करते हैं और $\epsilon \to 0$ की सीमा तक जाने के लिए अर्ज़ेला-एस्कोली कॉम्पैक्टनेस तर्क (Arzelà–Ascoli compactness arguments) का उपयोग करते हैं, जिससे वैश्विक अस्तित्व स्थापित होता है।
विशिष्टता (Uniqueness): वे ऊर्जा फलन (energy functional) का विश्लेषण करके और धनात्मक प्रारंभिक डेटा के लिए ग्रोनवाल के लेम्मा (Grönwall's lemma) का उपयोग करके विशिष्टता सिद्ध करते हैं, तथा शून्य प्रारंभिक डेटा के लिए एक फिक्स्ड-पॉइंट प्रमेय (Precup's Theorem 11.3) का उपयोग करते हुए, उस सिंगुलैरिटी को सावधानीपूर्वक संबोधित करते हैं जहाँ $t=0$ पर लिप्सचिट्ज़ स्थिति विफल हो जाती है।
वेल-पोज़डनेस (Well-posedness): वे अधिकतम नॉर्म (maximum norm) में प्रारंभिक डेटा पर निरंतर निर्भरता प्रदर्शित करते हैं, जो हैडेमैन (Hadamard) के मानदंडों को पूरा करता है।

D. स्थिरता विश्लेषण (Stability Analysis)

रैखिक स्थिरता (Linear Stability): सिस्टम को साम्यावस्था बिंदु $(u_e, v_e) = (1/2, 0)$ के आसपास लीनियर किया जाता है। लेखक निर्धारित करते हैं कि साम्यावस्था स्पर्शोन्मुख रूप से स्थिर (asymptotically stable) है। अभिसरण (convergence) का प्रकार (एकदिष्ट बनाम दोलनकारी) महत्वपूर्ण मान $\omega^* = \beta^2/4$ पर निर्भर करता है।
वैश्विक स्थिरता (Global Stability): विशिष्ट प्रारंभिक डेटा ( $H(0)=\alpha, H'(0)=0$ ) के लिए वैश्विक अभिसरण को सिद्ध करने के लिए, वे एक लियापुनोव फलन (Lyapunov function) का निर्माण करते हैं:
$V(u, v) = \frac{1}{2}v^2 - u + \frac{2\sqrt{2}}{3}u^{3/2} + \frac{1}{6}$
लासेल के अपरिवर्तनीयता सिद्धांत (LaSalle's invariance principle) का उपयोग करते हुए, वे आकर्षण क्षेत्र (basin of attraction) $\Omega_0$ की पहचान करते हैं जिसमें स्पष्ट रूप से भौतिक रूप से प्रासंगिक प्रारंभिक अवस्थाएं ( $0 \le \alpha \le 3/2$ ) शामिल हैं।

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

कठोर व्युत्पत्ति (Rigorous Derivation): संरक्षण नियमों (conservation laws) से स्लिप के साथ वॉशबर्न समीकरण की व्युत्पत्ति, जो सामान्यीकरण के लिए एक ठोस आधार प्रदान करती है।
पिछले प्रमाणों में सुधार: यह शोध पत्र वैश्विक अस्तित्व के पिछले प्रमाण (Ref [9]) में त्रुटियों की पहचान करता है और उन्हें सुधारता है, विशेष रूप से फिक्स्ड-पॉइंट ऑपरेटर के सेल्फ-मैपिंग गुण को सत्यापित करने में विफलता और $t=0$ पर मोनोटोनिसिटी थ्योरम के गलत अनुप्रयोग के संबंध में।
सामान्य प्रारंभिक डेटा: प्रमाण शून्य प्रारंभिक ऊंचाई के बजाय, प्रारंभिक डेटा की एक व्यापक श्रेणी ( $0 \le \alpha \le 3/2$ ) के लिए मान्य हैं।
स्लिप पैरामीटर विश्लेषण: यह प्रदर्शन करना कि जबकि स्लिप पैरामीटर $\beta$ अभिसरण की दर और दोलन व्यवहार के महत्वपूर्ण दहलीज को प्रभावित करता है, यह गुणात्मक दीर्घकालिक व्यवहार (साम्यावस्था की ओर अभिसरण) या स्वयं साम्यावस्था ऊंचाई को नहीं बदलता है।
वैश्विक स्थिरता प्रमाण: एक लियापुनोव फलन का निर्माण जो कठोरता से सिद्ध करता है कि सिस्टम केशिका उत्थान प्रयोगों में उपयोग किए जाने वाले विशिष्ट प्रारंभिक डेटा से साम्यावस्था की ओर अभिसरण करता है, जो कि पहले गणितीय प्रमाण के अभाव में था।

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

अस्तित्व और विशिष्टता: किसी भी $\beta > 0, \omega > 0$ और प्रारंभिक ऊंचाई अनुपात $\alpha \in [0, 3/2]$ के लिए, सभी $T \in [0, \infty)$ के लिए परिभाषित एक अद्वितीय, बाध्य (bounded), धनात्मक समाधान $H(T)$ मौजूद है।
वेल-पोज़डनेस: प्रारंभिक मान समस्या (initial value problem) हैडेमैन के अर्थ में वेल-पोज़ड है (अस्तित्व, विशिष्टता और डेटा पर निरंतर निर्भरता)।
साम्यावस्था (Equilibrium): अद्वितीय साम्यावस्था ऊंचाई $H_e = 1$ (आयामहीन) है, जो स्लिप पैरामीटर $\beta$ से स्वतंत्र है।
गतिकी (Dynamics):
- यदि $\omega < \beta^2/4$ : समाधान एकदिष्ट रूप से (monotonically) साम्यावस्था की ओर बढ़ता है (स्थिर नोड)।
- यदि $\omega > \beta^2/4$ : समाधान दोलनकारी रूप से (oscillatorily) साम्यावस्था की ओर बढ़ता है (स्थिर स्पाइरल)।
- स्लिप पैरामीटर $\beta$ महत्वपूर्ण दहलीज $\omega^*$ को स्थानांतरित करता है, जो प्रभावी रूप से बढ़ते स्लिप (कम $\beta$ ) के साथ दोलनों को अधिक मजबूती से कम (damping) करता है।
अभिसरण (Convergence): विश्राम की स्थिति से शुरू होने वाले सभी समाधान ( $0 \le H(0) \le 1.5 H_e$ ) $t \to \infty$ होने पर साम्यावस्था ऊंचाई $H_e$ की ओर अभिसत होते हैं।

5. महत्व (Significance)

यह कार्य भौतिक मॉडलिंग और कठोर गणितीय विश्लेषण के बीच के अंतर को पाटता है।

भौतिक प्रासंगिकता: स्लिप स्थिति को शामिल करके, मॉडल हु-स्क्रीवन विरोधाभास को हल करता है, जिससे यह प्रवाह के प्रारंभिक चरण के लिए अधिक भौतिक रूप से सटीक बनता है जहाँ संपर्क रेखा गतिमान होती है।
गणितीय कठोरता: यह वॉशबर्न समीकरण के लिए पूर्ण और कठोर वैश्विक वेल-पोज़डनेस का पहला पूर्ण प्रमाण प्रदान करता है, जो अपूर्ण तर्कों को सुधारता है।
भविet भविष्य कहनेवाला शक्ति (Predictive Power): आकर्षण क्षेत्र की पहचान यह सुनिश्चित करती है कि मॉडल उन प्रयोगात्मक स्थितियों के लिए मजबूत है जहाँ तरल एक छोटी ऊंचाई से विश्राम की स्थिति से शुरू होता है।
भावी दिशाएँ: यह शोध पत्र परिवर्तनशील स्लिप मापदंडों, नकारात्मक स्लिप (संपर्क कोण हिस्टेरेसिस), और दीर्घकालिक सिमुलेशन के लिए स्थिर संख्यात्मक विधियों के विकास के विश्लेषण के लिए मार्ग प्रशस्त करता है।