Enucleated incompressible red blood cells in shear flow: theoretical analysis of shape instabilities

तकनीकी सारांश: शियर फ्लो (Shear Flow) में एन्यूक्लिएटेड (Enucleated) असंपीड्य लाल रक्त कोशिकाएं: आकार अस्थिरता का सैद्धांतिक विश्लेषण

समस्या विवरण
लाल रक्त कोशिकाएं (RBCs) ऑक्सीजन परिवहन के लिए आवश्यक हैं, जो प्रवाह में महत्वपूर्ण विरूपण (deformation) करने की अपनी क्षमता पर निर्भर करती हैं। जबकि रिजिड-बॉडी मॉडल (जैसे जेफ़री ऑर्बिट्स) निम्न शियर दरों पर RBC गतिकी (tumbling, tank-treading, swinging) का सफलतापूर्वक वर्णन करते हैं, वे मध्यवर्ती शियर दरों ($100\text{--}1000 \text{ s}^{-1}$) पर जटिल, बहु-लोब वाले आकारों और सिमेट्री-ब्रोकन मॉर्फोलॉजी (जैसे स्टोमेटोसाइट्स) के उद्भव को समझाने में विफल रहते हैं। इन आकार अस्थिरताओं को चलाने वाले भौतिक तंत्रों—विशेष रूप से झिल्ली लोच (membrane elasticity), बेंडिंग प्रतिरोध, झिल्ली तनाव, अतिरिक्त क्षेत्रफल और श्यानता अनुपात (viscosity ratios) के बीच परस्पर क्रिया—पर बहस जारी है। इसके अलावा, मौजूदा सैद्धांतिक ढांचे अक्सर RBC को एक स्थिर दीर्घवृत्त (ellipsoid) या कैप्सूल के रूप में मानते हैं, जिससे एन्यूक्लिएशन की परिपक्वता प्रक्रिया के परिणामस्वरूप होने वाले विशिष्ट ज्यामितीय परिवर्तनों की उपेक्षा होती है, जहाँ कोशिका केंद्रक (nucleus) को बाहर निकाल दिया जाता है, जिससे स्ट्रेस-फ्री झिल्ली क्षेत्रफल और अंतर्निहित आयतन के बीच संबंध बदल जाता है।

कार्यप्रणाली
लेखक एक न्यूटनियन शियर फ्लो में निलंबित एक बंद, तरल-भरे, असंपीड्य झिल्ली के रूप में RBC को मॉडल करने के लिए एक परिवर्तक सैद्धांतिक ढांचे (perturbative theoretical framework) का उपयोग करते हैं। यह दृष्टिकोण पूर्व में विलाहोस्का एट अल. [26] द्वारा स्थापित रैखिकीकृत गैर-रैखिक ढांचे का विस्तार करता है, जिसमें स्पष्ट रूप से दो अलग-अलग अतिरिक्त क्षेत्रों को शामिल किया गया है:

1. $\Delta_0$ (तनाव-मुक्त अतिरिक्त क्षेत्रफल): एन्यूक्लिएशन से पहले RBC के दीर्घवृत्ताकार आकार के अनुरूप।
2. $\Delta$ (वर्तमान अतिरिक्त क्षेत्रफल): एन्यूक्लिएशन के बाद डिस्कोसाइट आकार के अनुरूप, जहाँ आयतन घट जाता है ($V < V_0$) जबकि सतह का क्षेत्रफल संरक्षित रहता है, जिससे आंतरिक तनाव लॉक हो जाता है।

RBC के आकार को गोलाकार हार्मोनिक्स ($Y_{l,m}$) में विस्तारित किया जाता है, जो वर्तमान विन्यास $\rho$ और तनाव-मुक्त विन्यास $\rho_0$ को अलग करता है। गति के शासन समीकरणों को फ्लुइड स्ट्रेस (स्टोक्स फ्लो) और झिल्ली बलों के बीच संतुलन बनाकर व्युत्पन्न किया जाता है, जिन्हें इलास्टिक (स्कालाक मॉडल), बेंडिंग (हेल्फ्रिक ऊर्जा) और झिल्ली तनाव योगदानों में विभाजित किया गया है। झिल्ली तनाव ($\Sigma$) स्थानीय क्षेत्रफल को स्थिर रखने के लिए एक लैग्रेंज मल्टीप्लायर के रूप में कार्य करता है और विभिन्न स्थानिक मोड (spatial modes) को जोड़ता है।

विश्लेषण दो चरणों में आगे बढ़ता है:

1. न्यूनीकरण प्रणाली ($l=2$): सिस्टम को दीर्घवृत्ताकार मोड तक सीमित किया जाता है ताकि टंबलिंग, टैंक-ट्रेडिंग और स्विंगिंग की गतिकी को फिर से देखा जा सके, जिसमें $\Delta$ और $\Delta_0$ के बीच के अंतर को स्पष्ट रूप से शामिल किया गया है।
2. पूर्ण परिवर्तक प्रणाली ($l \geq 2$): एक रैखिक अस्थिरता मानदंड प्राप्त करने के लिए अतिरिक्त स्थानिक मोड पेश किए जाते हैं। सिस्टम की स्थिरता का विश्लेषण आयामहीन मापदंडों: बेंडिंग नंबर ($B$), कैपिलरी नंबर ($C$), श्यानता अनुपात ($\lambda$), और अतिरिक्त क्षेत्रों के फलन के रूप में मोड के विकास दर की जांच करके किया जाता है।

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

• सैद्धांतिक ढांचे का विस्तार: यह अध्ययन एन्यूक्लिएशन द्वारा प्रेरित ज्यामितीय और स्ट्रेस अवस्था परिवर्तनों को स्पष्ट रूप से ध्यान में रखते हुए मौजूदा परिवर्तक दृष्टिकोण का सामान्यीकरण करता है। तनाव-मुक्त आकार ($\Delta_0$) और वर्तमान आकार ($\Delta$) के बीच अंतर करके, मॉडल कोशिका की परिपक्वता प्रक्रिया की "स्मृति" (memory) को पकड़ता है, जो कि फ्लूइड वेसिकल मॉडल में अनुपस्थित है।
• अतिरिक्त क्षेत्रफल की न्यूनीकरण गतिकी में भूमिका: $l=2$ सबस्पेस में, लेखक प्रदर्शित करते हैं कि $\Delta$ और $\Delta_0$ के बीच का अंतर प्रभावी इलास्टिक स्थिरांक ($S$) को पुनर्रचना (renormalize) करता है। यह संदर्भ आकार के अतिरिक्त क्षेत्रफल में परिवर्तन के कारण होने वाले गतिकी संक्रमणों (जैसे, स्विंगिंग से टंबलिंग तक) को संचालित करने की अनुमति देता है, भले ही शियर दर स्थिर हो। प्रारंभिक ओरिएंटेशन (इन-प्लेन बनाम आउट-ऑफ-प्लेन) परिणामी प्रक्षेपवक्रों को महत्वपूर्ण रूप से प्रभावित करता है।
• अस्थिरता मानदंड और झिल्ली तनाव: आकार संक्रमणों के लिए एक रैखिक अस्थिरता मानदंड का व्युत्पन्न करना प्राथमिक योगदान है। विश्लेषण से पता चलता है कि आकार की अस्थिरता (बहु-लोब वाले आकारों का उद्भव) झिल्ली तनाव $\sigma(t)$ पर महत्वपूर्ण रूप से निर्भर है। एक अस्थिरता तब होती है जब झिल्ली तनाव एक महत्वपूर्ण सीमा ($\sigma < \sigma^*$) से नीचे गिर जाता है, जो शियर इलास्टिसिटी और बेंडिंग स्टिफनेस के अनुपात ($\alpha = B/C$) द्वारा निर्धारित होता है।
  • अस्थिरता तब शुरू होती है जब बेंडिंग और इलास्टिक प्रतिक्रिया बाहरी प्रवाह और आंतरिक तनावों का मुकाबला करने में अपर्याप्त होती है।
  • सबसे अस्थंत मोड ($l^*$) और महत्वपूर्ण तनाव $\alpha$, श्यानता अनुपात $\lambda$, और अतिरिक्त क्षेत्रों पर निर्भर करता है।
• संतृप्ति और पुनर्वितरण: अध्ययन अस्थिरता के बाद दीर्घकालिक व्यवहार का विश्लेषण करता है। अस्थिर मोड का विकास झिल्ली तनाव में वृद्धि करता है, जो अंततः सिस्टम को स्थिर करता है। यह प्रक्रिया अतिरिक्त क्षेत्रफल के पुनर्वितरण का कारण बनती है, जो प्रमुख $l=2$ मोड से उच्च-क्रम मोड (जैसे, $l=3, 4$) में स्थानांतरित हो जाता है, जिसके परिणामस्वरूप एक नया, स्थिर, गैर-दीर्घवृत्ताकार आकार विन्यास प्राप्त होता है।
• फेज डायग्राम (Phase Diagrams): युग्मित समीकरणों के संख्यात्मक एकीकरण से विश्लेषणात्मक भविष्यवाणियों की पुष्टि होती है। लेखक $(B, C)$ तल में फेज डायग्राम का निर्माण करते हैं, जो उन क्षेत्रों की पहचान करते हैं जहाँ विशिष्ट मोड ($l=3, 4, \dots$) अस्थिर हो जाते हैं। ये डायग्राम दिखाते हैं कि श्यानता अनुपात $\lambda$ और अतिरिक्त क्षेत्रफल $\Delta$ को बढ़ाने से उच्च कैपिलरी नंबर (कमजोर इलास्टिसिटी) पर RBC अस्थिर हो जाते हैं, जिससे जटिल आकारों का उद्भव होता है।

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