Volumetric Growth in Linear Elasticity Driven by an Optimality Criterion

तकनीकी सारांश: एक रैखिकीकृत लोचशीलता ढांचे (Linearized Elasticity Framework) में अनुकूलन-संचालित आयतनिक वृद्धि

समस्या का निरूपण
यह शोध पत्र लोचदार पिंडों (elastic bodies) में आयतनिक वृद्धि के मॉडलिंग को संबोधित करता है, जो एक ऐसी प्रक्रिया है जहाँ द्रव्यमान पूरे आयतन में उत्पन्न होता है, जिससे ज्यामिति और यांत्रिक प्रतिक्रिया बदल जाती है। पारंपरिक निरंतरता यांत्रिकी (continuum mechanics) दृष्टिकोण आमतौर पर प्रतिबल (stress), विकृति (strain), या जैव-रासायनिक उद्दीपनों पर निर्भर घटनात्मक गतिक नियमों (phenomenological kinetic laws) का उपयोग करके ग्रोथ टेंसर ($E_g$) के विकास को निर्धारित करते हैं। इसके विपरीत, यह कार्य एक ऐसा ढांचा प्रस्तावित करता है जहाँ ग्रोथ टेंसर को पूर्व-निर्धारित नहीं किया जाता है, बल्कि प्रत्येक विविक्त समय चरण (discrete time step) पर एक प्रतिबन्धित अनुकूलन समस्या (constrained optimization problem) के समाधान के रूप में अप्रत्यक्ष रूप से निर्धारित किया जाता है।

सेटिंग दो आयामी रैखिकीकृत लोचशीलता (linearized elasticity) की है (जिसे तीन आयामों में आसानी से विस्तारित किया जा सकता है)। पिंड एक डोमेन $\Omega$ में स्थित है जिसकी सीमा को डिरिचलेट ($\partial_D \Omega$) और न्यूमैन ($\partial_N \Omega$) भागों में विभाजित किया गया है। वृद्धि को स्ट्रेन टेंसर के योगात्मक अपघटन (additive decomposition) के माध्यम से मॉडल किया गया है, जहाँ कुल स्ट्रेन, लोचदार स्ट्रेन (elastic strain) और एक अलोचदार वृद्धि स्ट्रेन ($E_g$) का योग है। संतुलन समीकरणों को इस प्रकार संशोधित किया गया है कि $E_g$ एक आंतरिक स्रोत पद (internal source term) के रूप में कार्य करता है।

मुख्य समस्या प्रत्येक समय चरण $i$ पर विस्थापन क्षेत्र $u^{(i)}$ और ग्रोथ टेंसर क्षेत्र $E_g^{(i)}$ पर एक वृद्धिशील न्यूनीकरण (incremental minimization) है। अनुकूलन तीन प्राथमिक प्रतिबंधों के अधीन है:

1. संतुलन: रैखिकीकृत लोचशीलता के संतुलन समीकरणों को संतुष्ट किया जाना चाहिए।
2. द्रव्यमान संतुलन: वृद्धि के कारण कुल (या स्थानीय) आयतन परिवर्तन एक निर्धारित वृद्धि $\Gamma^{(i)}$ के बराबर होना चाहिए।
3. अपरिवर्तनीयता (संचय/Accretion): वृद्धि विशुद्ध रूप से संचयी होनी चाहिए, जिसका अर्थ है कि ग्रोथ टेंसर का अंतर धनात्मक अर्ध-निश्चित (positive semi-definite) होना चाहिए ($\lambda_{\min}(E_g^{(i)} - E_g^{(i-1)}) \geq 0$)।

ऑब्जेक्टिव फंक्शन $\Phi$ प्रेरक तंत्र (driving mechanism) को दर्शाता है। यह शोध पत्र दो विशिष्ट मानदंडों की जांच करता है:

• बाह्य कार्य का न्यूनीकरण: यह एक यांत्रिक रूप से संचालित अनुकूलन का प्रतिनिधित्व करता है जहाँ पिंड कठोरता (stiffness) बढ़ाने का प्रयास करता है।
• परिधि का न्यूनीकरण: यह एक ज्यामितीय रूप से संचालित विकास (आइसोपिरिमेट्रिक सिद्धांत) का प्रतिनिधित्व करता है।

एक द्विघातीय नियमितीकरण पद (quadratic regularization term), जो $|E_g^{(i)} - E_g^{(i-1)}|^2$ के समानुपाती है, शामिल किया गया है ताकि पिछले विन्यास से बड़े विचलन को दंडित किया जा सके। डी जॉर्गी (De Giorgi) के 'मिनिमाइजिंग मूवमेंट्स' के सिद्धांत से प्रेरित यह पद, लौकिक निरंतरता (temporal continuity) सुनिश्चित करता है और निरंतर सीमा में एक ग्रेडिएंट-फ्लो संरचना की ओर ले जाता है।

कार्यप्रणाली
लेखक निरंतर समस्या को विविक्त करने के लिए परिमित तत्व विधि (Finite Element Method - FEM) का उपयोग करते हैं।

• विवक्तीकरण (Discretization): विस्थापन क्षेत्र को निरंतर, खंडिक रैखिक (piecewise affine) आधार फलनों (basis functions) का उपयोग करके अनुमानित किया जाता है, जबकि ग्रोथ टेंसर को त्रिकोणीय तत्वों पर खंडिक स्थिर (piecewise constant) आधार फलनों का उपयोग करके अनुमानित किया जाता है।
• न्यूनीकरण (Reduction): विस्थापन क्षेत्र को ग्रोथ टेंसर के संदर्भ में रैखिक संतुलन समीकरणों को स्पष्ट रूप से हल करके, लेखक विस्थापन चरों को हटा देते हैं। यह समस्या को केवल ग्रोथ चरों वाले एक परिमित-आयामी प्रतिबन्धित न्यूनीकरण समस्या में बदल देता है।
• अनुकूलन: परिणामी विविक्त समस्या को MATLAB के fmincon रूटीन के साथ 'सीक्वेंशियल क्वाड्रेटिक प्रोग्रामिंग' (SQP) एल्गोरिदम का उपयोग करके संख्यात्मक रूप से हल किया जाता है। ऑब्जेक्टिव फंक्शन और प्रतिबंधों के लिए विश्लेषणात्मक ग्रेडिएंट प्रदान किए गए हैं ताकि अभिसरण (convergence) सुनिश्चित हो सके।
• विश्लेषणात्मक सन्निकटन: अंतर्निहित संरचना को प्रकट करने के लिए, लेखक गैर-रैखिक आइजनवैल्यू प्रतिबंध (अपरिवर्तनीयता) को अस्थायी रूप से अनदेखा करके एक अनुमानित विश्लेषणात्मक समाधान व्युत्पन्न करते हैं। यह ग्रोथ वृद्धि के लिए एक बंद-रूप अभिव्यक्ति (closed-form expression) के व्युत्पन्न की अनुमति देता है, जो यह प्रदर्शित करता है कि विकास अलोचदार स्ट्रेन के स्थान में एक 'प्रोजेक्टेड ग्रेडिएंट फ्लो' का अनुसरण करता है।

प्रमुख परिणाम
संख्यात्मक प्रयोग तीन प्रतिनिधि मामलों पर किए गए: एक डबल-क्लैम्प्ड बीम, एक कैंटिलीवर बीम (दोनों बाहरी भार के तहत), और एक मुक्त पिंड (परिधि न्यूनीकरण हेतु)।

• बीम अनुकूलन: बाहरी भार के तहत, बीम संरचनात्मक कठोरता बढ़ाने के लिए विकसित होते हैं। ग्रोथ वितरण गैर-समान है, जो उन क्षेत्रों में केंद्रित होता है जहाँ भौतिक फाइबर लंबे होते हैं। डबल-क्लैम्प्ड बीम मुख्य रूप से संपीड़न अवशिष्ट प्रतिबल (compressive residual stresses) विकसित करता है, जबकि कैंटिलीवर वृद्धि को विस्तार (elongation) के माध्यम से समायोजित करता है, जिससे प्रतिबल के परिमाण में कमी आती है।
• परिधि न्यूनीकरण: मुक्त पिंड एक गोलाकार आकृति की ओर विकसित होता है, जो शास्त्रीय आइसोपिरिमेट्रिक सिद्धांत के अनुरूप है।
• प्रतिबल पैटर्न: परिधि-संचालित मामले में, अवशिष्ट प्रतिबल क्षेत्र एक संपीड़ित कोर (compressive core) और एक तन्यता बाहरी खोल (tensile outer shell) प्रदर्शित करता है। लेखक नोट करते हैं कि यह पैटर्न चूहों और मनुष्यों के ट्यूमर में देखे गए प्रयोगात्मक और संख्यात्मक अवलोकनों को दर्शाता है, जहाँ ट्यूमर का कोर संपीड़ित होता है और बाहरी खोल तनाव में होता है।
• विश्लेषणात्मक बनाम संख्यात्मक: अनुमानित विश्लेषणात्मक समाधान (जो अपरिवर्तनीयता प्रतिबंध के बिना प्राप्त किए गए थे) बीम मामलों के लिए पूर्ण संख्यात्मक परिणामों के साथ निकटता से मेल खाते हैं, जिनमें ऑब्जेक्टिव फंक्शन मानों और मिनिमाइज़र में नगण्य अंतर है। परिधि मामले के लिए, यदि नियमितीकरण पैरामीटर पर्याप्त रूप से बड़ा है, तो विश्लेषणात्मक समाधान आइजनवैल्यू प्रतिबंध के प्रारंभिक उल्लंघन के बावजूद संख्यात्मक परिणाम की ओर अभिसरित होता है।

महत्व और दावे
यह शोध पत्र यह दावा करने का प्रयास करता है कि आयतनिक वृद्धि को एक अनुकूलन-संचालित प्रक्रिया के रूप में तैयार किया जा सकता है जहाँ ग्रोथ टेंसर को किसी घटनात्मक नियम के बजाय एक यांत्रिक चयन सिद्धांत द्वारा चुना जाता है।

• सैद्धांतिक बदलाव: यह कार्य इस बात से अलग है कि यह वृद्धि को केवल एक वर्णनात्मक गतिक इनपुट के रूप में नहीं, बल्कि यांत्रिक बाधाओं के अधीन एक विशिष्ट फंक्शनल के न्यूनीकरण के परिणाम के रूप में मानता है।
• संरचनात्मक अंतर्दृष्टि: मॉडल यह प्रकट करता है कि विकास को अलोचदार स्ट्रेन के स्थान में एक 'प्रोजेक्टेड ग्रेडिएंट फ्लो' के रूप में व्याख्यायित किया जा सकता है। लेखक इस बात पर जोर देते हैं कि कंप्लायंस फंक्शनल (compliance functional) की उत्तलता (convexity) अत्यंत महत्वपूर्ण है; जब उत्तलता समाप्त हो जाती है, तो गैर-विशिष्टता (non-uniqueness) और स्थानीयकरण (localization) उत्पन्न हो सकते हैं।
• प्रूफ ऑफ कॉन्सेप्ट: रैखिकीकृत सिद्धांत एक 'प्रूफ ऑफ कॉन्सेप्ट' के रूप में कार्य करता है, जो यह दर्शाता है कि विशिष्ट जैविक विकास नियमों को माने बिना, केवल यांत्रिक बाधाओं वाले अनुकूलन ही जटिल, वृद्धि-प्रेरित रूपात्मकता (जैसे आकार अनुकूलन और अवशिष्ट प्रतिबल पैटर्न) उत्पन्न करने के लिए पर्याप्त है।
• भविष्य की संभावना: लेखक कहते हैं कि यह ढांचा पूर्णतः गैर-रैखिक लोचशीलता और प्रतिस्पर्धी ऊर्जा उद्देश्यों वाले अधिक यथार्थवादी तंत्रों में इस दृष्टिकोण के विस्तार का मार्ग प्रशस्त करता है, हालांकि वे वर्तमान कार्य में इन जटिल जैविक प्रणालियों को हल करने का दावा नहीं करते हैं।

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