🔬 condensed matter

Richards' equation as a hydrodynamic limit: Chapman--Enskog reduction of the continuum kinetic equation for unsaturated soil water

यह शोध पत्र एक निरंतर गतिज समीकरण (continuum kinetic equation) के प्रथम-क्रम चैपमैन-एनस्कोग (Chapman-Enskog) न्यूनीकरण के रूप में असंतृप्त मृदा जल प्रवाह के लिए रिचर्ड्स के समीकरण को व्युत्पन्न करता है, जो पोर-स्केल गतिकी से इसके चालकता जैसे स्थिरांकों को स्थापित करता है और ड्यूल-परमीबिलिटी मॉडलों के लिए एक कठोर प्रथम-सिद्धांत आधार प्रदान करता है।

मूल लेखक: Riccardo Rigon

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

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

CC BY 4.0

मूल लेखक: Riccardo Rigon

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

मिट्टी में पानी का अदृश्य नृत्य

कल्पना कीजिए कि आपके पैरों के नीचे की ज़मीन एक ठोस मिट्टी के ब्लॉक के रूप में नहीं, बल्कि अरबों नन्हे, अदृश्य सुरंगों से बने एक विशाल, त्रि-आयामी स्पंज के रूप में है। यह मृदा भौतिकी (soil physics) की दुनिया है, जो इस बात को समझने के लिए समर्पित है कि पानी इन सूक्ष्म भूलभुलैयाओं के माध्यम से कैसे चलता है। लगभग एक सदी से, वैज्ञानिक इस गति की भविष्यवाणी करने के लिए रिचर्ड्स समीकरण (Richards' equation) नामक एक प्रसिद्ध नियम का उपयोग कर रहे हैं। इस समीकरण को मिट्टी के लिए मौसम के पूर्वानुमान की तरह समझें: यह हमें बताता है कि मिट्टी कितनी गीली होगी और पानी कितनी तेज़ी से निकल जाएगा। यह मिट्टी को एक चिकने, निरंतर पदार्थ के रूप में मानकर काम करता है, जिसमें दो मुख्य सामग्रियां शामिल हैं: मिट्टी में वर्तमान में कितना पानी है (जल सामग्री/water content) और मिट्टी उस पानी को कितनी ज़ोर से खींच रही है (मैट्रिक क्षमता/matric potential)।

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

शोध पत्र की बड़ी खोज: ज़मीन से ही नियमों को व्युत्पन्न करना

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

इस शोध पत्र का मूल भाग चैपमैन-एनस्कोग (Chapman–Enskog) रिडक्शन नामक एक गणितीय यात्रा है। इसे एक अराजक, तेज़ी से चलती भीड़ को एक चिकनी, बहती हुई नदी में बदलने के तरीके के रूप में समझें। लेखक समस्या को दो भागों में विभाजित करता है: मिट्टी के एक अकेले समूह के भीतर पानी की नन्ही, तेज़ गतिविधियाँ ("माइक्रो" स्केल) और पूरे खेत में पानी की धीमी, व्यापक-चित्र वाली गति ("मैक्रो" स्केल)। शोध पत्र यह सिद्ध करता है कि प्रसिद्ध रिचर्ड्स समीकरण केवल तभी मान्य होता है जब मिट्टी के अंदर का पानी मौसम के बदलने की तुलना में बहुत तेज़ी से खुद को व्यवस्थित कर सकता है। लेखक एक विशिष्ट संख्या पेश करता है, डैमकोहलर नंबर (Damköhler number - Da), जो इस दौड़ के लिए एक स्पीडोमीटर की तरह कार्य करता है। यदि पानी मुकाबला करने के लिए पर्याप्त तेज़ी से चलता है (छोटा Da), तो पुराना नियम पूरी तरह से काम करता है। लेकिन यदि मौसम बहुत तेज़ी से बदलता है या मिट्टी बहुत अधिक अवरुद्ध है, तो पानी तालमेल नहीं बिठा पाता, "चिकनी नदी" की धारणा टूट जाती है, और पुराना समीकरण विफल हो जाता है।

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

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

अंततः, यह शोध पत्र रिचर्ड्स समीकरण को एक "नियम" से बदलकर एक मौलिक भौतिक नियम में बदल देता है, ठीक वैसे ही जैसे तरल गतिकी (fluid dynamics) के नियमों को व्यक्तिगत गैस अणुओं की गति से व्युत्पन्न किया जाता है। यह पुष्टि करता है कि पुराना समीकरण सही परिस्थितियों में सही है, लेकिन यह सटीक गणितीय उपकरण भी प्रदान करता है कि उन परिस्थितियों के टूटने पर क्या होगा और वास्तविक दुनिया की जटिलताओं के लिए बेहतर मॉडल कैसे बनाए जाएं। मिट्टी के पानी को एक काइनेटिक सिस्टम के रूप में मानकर, लेखक ने एक सदी पुराने अनुभवजन्य सूत्र को एक व्युत्पन्न सत्य में बदल दिया है, जिसके साथ एक स्पष्ट मानचित्र भी है कि यह कहाँ तक जाता है और कहाँ इसकी सीमा समाप्त होती है।

तकनीकी सारांश: रिचर्ड्स के समीकरण का एक हाइड्रोडायनामिक सीमा (Hydrodynamic Limit) के रूप में निरूपण

समस्या विवरण
रिचर्ड्स का समीकरण लगभग एक सदी से असंतृप्त जल प्रवाह (unsaturated water flow) के लिए मानक मैक्रोस्कोपिक मॉडल के रूप में कार्य कर रहा है। हालाँकि, इसे पारंपरिक रूप से मैक्रोस्कोपिक द्रव्यमान संरक्षण और डार्सी-बकमिंगहैम फ्लक्स नियम से अनुमानित किया जाता है, न कि किसी सूक्ष्म सिद्धांत (microscopic theory) से व्युत्पन्न किया जाता है। फलस्वरूप, इसकी वैधता का क्षेत्र अनिर्दिष्ट बना हुआ है, और यह प्रिफरेंशियल फ्लो (preferential flow), हिस्टेरेसिस (hysteresis), और डायनेमिक नॉन-इक्विलिब्रियम प्रभावों जैसी घटनाओं को पकड़ने में विफल रहता है। ये विफलताएं इस धारणा से उत्पन्न होती हैं कि आंतरिक पोर-स्केल (pore-scale) विन्यास स्थायी रूप से स्थानीय संतुलन में है, जिससे जल सामग्री ( $\theta$ ) या मैट्रिक्स क्षमता ( $\psi$ ) पर्याप्त अवस्था चर (state variables) बन जाते हैं। यह शोध पत्र मृदा जल के लिए एक निरंतरता गतिज समीकरण (continuum kinetic equation - CKE) के हाइड्रोडायनामिक सीमा के रूप में रिचर्ड्स के समीकरण को कठोरता से व्युत्पन्न करके सूक्ष्म पोर-स्केल भौतिकी और मैक्रोस्कोपिक निरंतरता विवरण के बीच के अंतर को संबोधित करता है।

कार्यप्रणाली
व्युत्पत्ति एक दो-चरणीय सीमा प्रक्रिया (two-step limiting process) का उपयोग करती है, जो स्थानिक समरूपीकरण (spatial homogenization) को टेम्पोरल स्केल सेपरेशन (temporal scale separation) से अलग करती है, जो बोल्ट्ज़मैन समीकरण से नेवियर-स्टोक्स समीकरणों के व्युत्पत्ति के समान है।

स्थानिक सीमा ( $\varepsilon \to 0$ ): प्रतिनिधि प्राथमिक आयतन (Representative Elementary Volume - REV) को एक बिंदु में संकुचित किया जाता है। यह विशुद्ध रूप से काइनेमैटिक चरण मेसोस्केल मास्टर समीकरण को एक बंद निरंतरता गतिज समीकरण (CKE) में परिवर्तित करता है:
$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$
यहाँ, $g(r, x, t)$ त्रिज्या $r$ वाले छिद्रों का अंश है जो जल से भरे हुए हैं। फ्लक्स $\mathbf{F}$ एक "एंटैंगल्ड" (entangled) प्री-क्लोजर फ्लक्स है जहाँ परिवहन गुणांक और ड्राइविंग ग्रेडिएंट अभी अलग नहीं किए गए हैं, और $C[g]$ इंट्रा-REV पुनर्वितरण ऑपरेटर है।
टेम्पोरल सीमा (चैपमैन-एनस्कॉग रिडक्शन): यह शोध पत्र डैमोक्लर संख्या (Damköhler number), $Da = \tau_{redis}/\tau_{forcing}$ , यानी पोर-स्केल पुनर्वितरण समय और मैक्रोस्कोपिक फोर्सिंग समय के अनुपात के नियंत्रण में एक चैपमैन-एनस्कॉग (Chapman–Enskog - CE) विस्तार करता है। व्युत्पत्ति मानती है कि $Da \ll 1$ , जिसका अर्थ है कि मैक्रोस्कोपिक फोर्सिंग की तुलना में इंट्रा-REV पुनर्वितरण तीव्र है।
- $g$ का विस्तार $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$ के रूप में किया जाता है।
- शून्य क्रम ( $g^{(0)}$ ) स्थानीय संतुलन स्टेप फंक्शन (सभी छिद्र एक महत्वपूर्ण त्रिज्या $r^*$ से छोटे हैं तो पूर्ण हैं, अन्य खाली हैं) के अनुरूप है।
- प्रथम क्रम ( $g^{(1)}$ ) एक रैखिककृत पुनर्वितरण ऑपरेटर $J = C'|_{g_{eq}}$ को उलटने (inverting) द्वारा निर्धारित होता है।

मुख्य योगदान और परिणाम

रिचर्ड्स के समीकरण का व्युत्पत्ति: यह शोध पत्र सिद्ध करता है कि जब CKE पर CE विस्तार लागू किया जाता है, तो रिचर्ड्स का समीकरण ही एकमात्र मैक्रोस्कोपिक समीकरण प्राप्त होता है जब $Da \ll 1$ की स्थिति होती है। प्रथम-क्रम सुधार की समाधान योग्यता स्थिति (solvability condition) (फ्रेडहोम अल्टरनेटिव) मैक्रोस्कोपिक द्रव्यमान संरक्षण समीकरण प्रदान करती है, जो सिस्टम को क्लोज करती है।
हाइड्रोलिक कंडक्टिविटी एक परिवहन गुणांक के रूप में: हाइड्रोलिक कंडक्टिविटी $K(\theta)$ को केवल एक अनुभवजन्य (empirical) रचना संबंध के रूप में नहीं, बल्कि एक प्रथम-क्रम परिवहन गुणांक के रूप में व्युत्पन्न किया गया है, जो गैस काइनेटिक थ्योरी में श्यानता (viscosity) के संरचनात्मक रूप से समान है। इसे इस प्रकार व्यक्त किया गया है:
$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$
जहाँ $\kappa$ पोर-स्केल कंडक्टेंस है, $J^{-1}$ रैखिककृत पुनर्वितरण ऑपरेटर का व्युत्क्रम है, और $S_0 C_w$ स्रोत पद (source term) है। यह सूत्र $K$ को माध्यम के रिलैक्सेशन स्पेक्ट्रम के गुण के रूप में प्रकट करता है।
शास्त्रीय मॉडलों की रिकवरी:
- रिटेंशन कर्व्स: इक्विलिब्रियम स्टेप फंक्शन यंग-लाप्लास नियम के माध्यम से रिटेंशन कर्व $\psi(\theta)$ को स्वाभाविक रूप से प्राप्त करता है। ब्रूक्स-कोरी (Brooks–Corey) और वैन गेन्युचेन (van Genuchten) जैसे विशिष्ट रूपों को अंतर्निहित पोर-साइज डिस्ट्रीब्यूशन (PSD) के कोऑर्डिनेट ट्रांसफॉर्मेशन के रूप में पहचाना गया है।
- मुआलेम का विषमता दंड (Mualem's Heterogeneity Penalty): यह शोध पत्र लॉग-नॉर्मल PSDs के लिए मुआलेम फैक्टर $\exp(-4\sigma^2)$ को CE परिवहन गुणांक के सीरियल-पाथ करेक्शन (न्यूमैन सीरीज़) के एक सटीक रिसम्प्शन (resummation) के रूप में व्युत्पन्न करता है, न कि एक अनुभवजन्य फिट के रूप में।
ड्यूल- और मल्टीपल-परमीएबिलिटी मॉडल: यह ढांचा उन मामलों के लिए विस्तारित है जहाँ रिलैक्सेशन स्पेक्ट्रम में कई अंतराल (जैसे द्वि-मोडल PSDs) होते हैं। CKE को अलग-अलग स्पेक्ट्रल बैंड्स (तेज "बॉडी" पोर और धीमे "थ्रोट" पोर) पर प्रोजेक्ट करके और प्रत्येक बैंड के भीतर CE लागू करके, जबकि $O(1)$ पर क्रॉस-बैंड रिलैक्सेशन को बनाए रखा जाता है, यह शोध पत्र युग्मित रिचर्ड्स समीकरणों को व्युत्पन्न करता है। यह जेर्के-वैन गेन्युचेन ड्यूल-परमीएबिलिटी मॉडल, वीलर के IN3M, और मोबाइल-इमोबाइल मॉडलों को प्रथम सिद्धांतों से प्राप्त करता है। महत्वपूर्ण रूप से, इंटर-बैंड एक्सचेंज गुणांक को फिट करने के बजाय पोर कनेक्टिविटी से गणना किया जाता है।
ब्रेकडाउन के तंत्र: यह शोध पत्र रिडक्शन की सीमाओं को स्पष्ट रूप से वर्णित करता है। विस्तार तब विफल हो जाता है जब:
1. $Da \sim 1$ (काइनेटिक फोर्सिंग), जिससे तीव्र फ्रंट और नॉन-इक्विलिब्रियम अवस्थाएं उत्पन्न होती हैं।
2. स्पेक्ट्रल गैप पर्कोलेशन थ्रेशोल्ड ( $\theta \to \theta_c$ ) पर बंद हो जाता है, जिससे रिलैक्सेशन टाइम अनंत हो जाता है और कंडक्टिविटी शून्य हो जाती है।
3. स्केल सेपरेशन समाप्त हो जाता है (कनेक्टिविटी लंबाई REV आकार के करीब पहुँच जाती है)।

महत्व
यह शोध पत्र दावा करता है कि रिचर्ड्स का समीकरण केवल एक मॉडल नहीं है, बल्कि मृदा जल के काइनेटिक थ्योरी का विशिष्ट हाइड्रोडायनामिक लिमिट है, जो सख्ती से तेज़ स्थानीय पुनर्वितरण ( $Da \ll 1$ ) की स्थिति के तहत मान्य है। स्थानिक और टेम्पोरल सीमाओं को अलग करके, यह व्युत्पत्ति उन संरचनात्मक धारणाओं (स्थानीय संतुलन, सांख्यिकीय समरूपता) को स्पष्ट करती है जो समीकरण के मान्य होने के लिए आवश्यक हैं।

यह कार्य उन घटनाओं के लिए एक प्रथम-सिद्धांत आधार प्रदान करता है जिन्हें अक्सर घटनात्मक (phenomenological) रूप से माना जाता है। यह प्रदर्शित करता है कि ड्यूल-परमीएबिलिटी मॉडल पोर नेटवर्क में मल्टी-बैंड स्पेक्ट्रल संरचनाओं से स्वाभाविक रूप से उत्पन्न होते हैं और हिस्टेरेसिस का मूल कारण काइनेटिक थ्योरी में वेटिंग (wetting) और ड्राइंग (drying) ऑपरेटर्स की नॉन-कम्यूटेटिविटी (non-commutativity) में निहित है। इसके अतिरिक्त, यह हाइड्रोलिक कंडक्टिविटी को रिलैक्सेशन स्पेक्ट्रम के एक संरचनात्मक गुण के रूप में पहचानता है, जो पर्कोलेशन थ्रेशोल्ड के पास इसके व्यवहार और पोर ज्योमेट्री पर इसकी निर्भरता के लिए एक कठोर स्पष्टीकरण प्रदान करता है। यह व्युत्पत्ति शास्त्रीय मृदा भौतिकी परिणामों (रिटेंशन कर्व्स, कंडक्टिविटी इंटीग्रल्स) को अंतर्निहित काइनेटिक ऑपरेटर आर्किटेक्चर के परिणाम के रूप में एकीकृत करती है, ठीक वैसे ही जैसे बोल्ट्ज़मैन समीकरण से नेवियर-स्टोक्स समीकरण उभरते हैं।

मिट्टी में पानी का अदृश्य नृत्य

शोध पत्र की बड़ी खोज: ज़मीन से ही नियमों को व्युत्पन्न करना

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