🔬 materials science

Bridging perturbation and variational approaches in brittle fracture

यह शोध पत्र एक वेरिएशनल रिड्यूस्ड-ऑर्डर मॉडल प्रस्तुत करता है जो विषम माध्यमों में 3D भंगुर दरार प्रसार (brittle crack propagation) को कुशलतापूर्वक सिम्युलेट करने के लिए पर्टर्बेशन और वेरिएशनल फ्रैक्चर सिद्धांतों के बीच सेतु बनाता है, जिससे यह प्रकट होता है कि अव्यवस्था की तीव्रता (disorder intensity) और मोड मिक्सिटी (mode mixity), सुचारू से रुक-रुक कर होने वाली वृद्धि (intermittent growth) के संक्रमण और अव्यवस्था-प्रेरित दुर्बलता एवं सुदृढ़ीकरण (weakening and toughening) के बीच के क्रॉसओवर को कैसे नियंत्रित करते हैं।

मूल लेखक: Serafim Egorov, Antoine Sanner, Jean Sulem, Lars Pastewka, Mathias Lebihain

प्रकाशित 2026-05-14

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

CC BY 4.0

मूल लेखक: Serafim Egorov, Antoine Sanner, Jean Sulem, Lars Pastewka, Mathias Lebihain

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

कल्पना कीजिए कि आप जेली-ओ (Jell-O) के एक ब्लॉक के माध्यम से कांच के एक नुकीले, ऊबड़-खाबड़ टुकड़े को धकेलने की कोशिश कर रहे हैं, जिसके अंदर हार्ड कैंडी और सॉफ्ट मार्शमैलो के टुकड़े बेतरतीब ढंग से बिखरे हुए हैं। जैसे ही आप धकेलते हैं, कांच की दरार मक्खन के माध्यम से कटने वाले चाकू की तरह सुचारू रूप से आगे नहीं बढ़ती। इसके बजाय, यह हार्ड कैंडी पर अटक जाती है, दबाव बनाती है, और फिर अचानक अगले स्थान पर "कूद" (snap) जाती है, और फिर से अटक जाती है। वास्तविक, जटिल सामग्रियों जैसे चट्टान, कंक्रीट या हड्डी में दरारें इसी तरह चलती हैं।

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

यहाँ उनके काम का रोजमर्रा के उदाहरणों का उपयोग करके विवरण दिया गया है:

1. समस्या: "बहुत धीमा" बनाम "बहुत सरल" की दुविधा

वैज्ञानिकों के पास मॉडल बनाने के दो मुख्य तरीके हैं:

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

शोध पत्र का समाधान: लेखकों ने एक "हाइब्रिड" मॉडल बनाया। उन्होंने "एज-ओनली" (केवल किनारे वाले) दृष्टिकोण की गति को "मेगा-मेश" दृष्टिकोण के कठोर ऊर्जा नियमों के साथ जोड़ा। उन्होंने एक वेरिएशनल रिड्यूस्ड-ऑर्डर मॉडल बनाया। इसे एक ऐसे GPS के रूप में सोचें जो भीड़ के केवल अग्रणी किनारे को ट्रैक करता है लेकिन यह भविष्यवाणी करने के लिए जटिल ट्रैफिक नियमों का उपयोग करता है कि भीड़ कहाँ जाम होगी या बहेगी, बिना हर एक व्यक्ति का अनुकरण किए।

2. यह कैसे काम करता है: "ऊर्जा न्यूनीकरण" (Energy Minimization) का खेल

कंप्यूटर "न्यूनतम ऊर्जा" का खेल खेलता है।

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

इसके बाद उन्होंने एक स्मार्ट सर्च एल्गोरिदम ("न्यूटन कंजुगेट ग्रेडिएंट" के साथ "ट्रस्ट रीजन") का उपयोग किया ताकि एकदम सही आकार खोजा जा सके।

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

3. उन्होंने क्या खोजा: "116,000 सिमुलेशन"

टीम ने यह देखने के लिए कि अव्यवस्थपूर्ण, यादृच्छिक सामग्रियों में दरारें कैसे व्यवहार करती हैं, एक सिंगल कंप्यूटर कोर पर 116,000 सिमुलेशन चलाए। यहाँ उनके मुख्य निष्कर्ष दिए गए हैं:

सुचारू से झटकेदार (Smooth to Jerky): जब दरार छोटी होती है, तो यह सुचारू रूप से चलती है। लेकिन जैसे-जैसे यह बड़ी होती जाती है, यह अनियमित व्यवहार करने लगती है—कुछ समय के लिए अटकना, फिर अचानक आगे कूद जाना। इसे "इंटरमिटेंसी" (intermittency) कहा जाता है।
"शीयर" प्रभाव: अधिकांश पिछले अध्ययन केवल सामग्री को खींचकर अलग करने (मोड I) को देखते थे। इस शोध पत्र ने घुमाने और फिसलने (मोड II और III) को देखा। उन्होंने पाया कि जब आप सामग्री को घुमाते हैं, तो दरार का अगला हिस्सा गोल नहीं रहता; बल्कि यह एक क्वासी-एलिप्टिक (अंडे के आकार का) आकार ले लेता है।
आकार का महत्व ("क्रॉसओवर"):
- छोटी दरारें: एक अव्यवस्थपूर्ण सामग्री में, छोटी दरारें वास्तव में बढ़ने के लिए अधिक आसान पाती हैं (कमजोरी)। वे कठोर स्थानों के चारों ओर आसानी से लहरा सकती हैं।
- बड़ी दरारें: एक बार जब दरार पर्याप्त बड़ी हो जाती है, तो वह कठोर स्थानों द्वारा "पिन" (रोक दी जाती) हो जाती है। उसे तोड़ने के लिए भारी दबाव बनाना पड़ता है। यह सामग्री को उसके वास्तविक स्वरूप से अधिक मजबूत (tougher) दिखाता है।
- बदलाव (The Switch): एक विशिष्ट आकार होता है जहाँ सामग्री "अव्यवस्था" (messiness) के कारण कमजोर होने से बदलकर उससे "मजबूत" होने की ओर स्विच करती है।

4. यह क्यों महत्वपूर्ण है (शोध पत्र के अनुसार)

यह विधि वैज्ञानिकों को घंटों के भीतर एक सिंगल कंप्यूटर पर लाखों सूक्ष्म अशुद्धियों के साथ दरारों का अनुकरण करने की अनुमति देती है, जिसमें पहले दिनों या हफ्तों का समय लगता था।

उन्होंने अपने गणित को नए, हाथ से निकाले गए सूत्रों के विरुद्ध सत्यापित किया ताकि यह सिद्ध हो सके कि यह काम करता है। उन्होंने दिखाया कि उनका मॉडल सही ढंग से भविष्यवाणी करता है:

दरारें कैसे कूदती और रुकती हैं (इंटरमिटेंसी)।
ऊर्जा कैसे संग्रहीत और मुक्त की जाती है (जैसे स्प्रिंग का टूटना)।
कैसे सामग्री की "अव्यवस्था" दरार के आकार के आधार पर उसकी समग्र मजबूती को बदल देती है।

संक्षेप में: उन्होंने एक तेज़, सटीक "क्रैक सिम्युलेटर" बनाया है जो दरार के सामने को बाधाओं के क्षेत्र में चलते हुए एक लचीले रबर बैंड की तरह मानता है, और उन्नत गणित का उपयोग करता है ताकि यह सुनिश्चित हो सके कि यह कभी भी भौतिकी के नियमों को न तोड़े। यह हमें समझने में मदद करता है कि कुछ सामग्रियां अचानक विफल क्यों हो जाती हैं और अन्य तनाव के तहत क्यों टिकी रहती हैं।

तकनीकी सारांश: भंगुर फ्रैक्चर में परटर्बेशन और वेरिएशनल दृष्टिकोणों का समन्वय

समस्या विवरण
पूर्णतः भंगुर, विषम (heterogeneous) ठोसों में त्रि-आयामी (3D) दरार प्रसार (crack propagation) का मॉडलिंग करना गणनात्मक रूप से अत्यंत कठिन है, क्योंकि दरार के अग्रभाग (crack front) पर सिंगुलर स्ट्रेस कंसन्ट्रेशन और अव्यवस्थित माध्यमों में फ्रैक्चर की बहु-स्तरीय प्रकृति (multiscale nature) मौजूद होती है। मौजूदा विधियाँ महत्वपूर्ण समझौते (trade-offs) का सामना करती हैं:

फेज़-फील्ड विधियाँ (फ्रांकोर्ट और मैरिगो, 1998 पर आधारित) विस्कस रेगुलाइजेशन (viscous regularization) के बिना विस्कस मीडिया में 3D प्रसार के लिए सुदृढ़ हैं, लेकिन इन्हें महंगे वॉल्यूमेट्रिक मेशिंग की आवश्यकता होती है और ये एक विसरित क्षति लंबाई पैमाने ( $\ell$ ) को पेश करती हैं, जो विषमता के आकार के भौतिकी के साथ परस्पर क्रिया कर सकती है और उसे बदल सकती है।
परटर्बेशन दृष्टिकोण (राइस, 1985 पर आधारित) तीक्ष्ण दरारों को बनाए रखते हैं और केवल दरार के अग्रभाग को मेश करके उच्च गणनात्मक दक्षता प्रदान करते हैं। हालाँकि, वे पारंपरिक रूप से ग्रिफिथ के मानदंड के लिए 'एड हॉक' विस्कस रेगुलाइजेशन पर निर्भर करते हैं, जो फ्रैक्चर डायनेमिक्स को प्रभावित करने वाले संख्यात्मक आर्टिफैक्ट्स (numerical artifacts) उत्पन्न कर सकते हैं। इसके अतिरिक्त, अधिकांश परटर्बेशन अध्ययन मोड I में अर्ध-अनंत दरारों पर केंद्रित हैं, जिससे परिमित-आकार प्रभावों (finite-size effects), विषमता की तीव्रता और 3D दरारों पर मिश्रित-मोड (I+II+III) लोडिंग के प्रभाव को पर्याप्त रूप से समझा नहीं जा सका है।

कार्यप्रणाली (Methodology)
लेखक एक वेरिएशनल रिड्यूस्ड-ऑर्डर मॉडल प्रस्तावित करते हैं जो फ्रांकोर्ट और मैरिगो (1998) के वेरिएशनल फॉर्मूलेशन को राइस (1985) के परटर्बेशन थ्योरी के साथ जोड़ता है। मुख्य उद्देश्य विस्कस रेगुलाइजेशन का सहारा लिए बिना, सिस्टम की कुल ऊर्जा को न्यूनतम करके साम्यावस्था दरार-अग्रभाग (equilibrium crack-front) कॉन्फ़िगरेशन की गणना करना है, जिसे इलास्टिक पोटेंशियल एनर्जी और विसरित ऊर्जा (dissipated energy) के योग के रूप में परिभाषित किया गया है।

वेरिएशनल परटर्बेटिव फॉर्मूलेशन:
- कुल ऊर्जा $\Pi_{tot}$ को अपरिवर्तनीयता बाधा (irreversibility constraint - दरार ठीक नहीं हो सकती) के अधीन न्यूनतम किया जाता है।
- पोटेंशियल एनर्जी ( $\Pi_{pot}$ ): पारंपरिक परटर्बेशन विधियों की तरह ऊर्जा रिलीज रेट ( $G$ ) का विस्तार करने के बजाय, लेखक एक संदर्भ त्रिज्या $a_0$ से विचलित एक अर्ध-वृत्ताकार दरार अग्रभाग के लिए स्वयं पोटेंशियल एनर्जी का एक एसिम्प्टोटिक एक्सपेंशन (asymptotic expansion) प्राप्त करते हैं। यह व्युत्पत्ति ब्यूकनर-राइस वेट फंक्शन थ्योरी और मिश्रित-मोड स्ट्रेस इंटेंसिटी फैक्टर्स (SIFs) के लिए गाओ (1988) के लीनियर थ्योरी का उपयोग करती है।
- डिसिपेटेड एनर्जी ( $\Pi_{dis}$ ): इसे दरार की सतह पर विषम फ्रैक्चर ऊर्जा क्षेत्र $G_c$ को एकीकृत करके गणना किया जाता है।
- कुशल मूल्यांकन: पोटेंशियल एनर्जी और उसके डेरिवेटिव का मूल्यांकन फास्ट फूरियर ट्रांसफॉर्म (FFT) का उपयोग करके किया जाता है, जो इलास्टिक इंटरैक्शन की नॉन-लोकल प्रकृति (फ्रैक्शनल लैपलेसियन और हिलबर्ट ट्रांसफॉर्म ऑपरेटरों द्वारा प्रतिनिधित्व) का लाभ उठाता है।
संख्यात्मक कार्यान्वयन (Numerical Implementation):
- समस्या को एक नॉन-कॉन्वेक्स, बॉक्स-कंस्ट्रेंड मिनिमाइजेशन समस्या के रूप में तैयार किया गया है।
- सॉल्वर: एक मैट्रिक्स-फ्री न्यूटन-कंजुगेट ग्रेडिएंट (CG) एल्गोरिदम का उपयोग ट्रस्ट रीजन के साथ किया जाता है।
- बाधाएं (Constraints):
  - बॉक्स कंस्ट्रेंट्स दरार की अपरिवर्तनीयता ( $a_i \ge a_i^{prev}$ ) को लागू करते हैं।
  - ट्रस्ट रीजन कंस्ट्रेंट्स (त्रिज्या को विषमता सहसंबंध लंबाई $d$ द्वारा निर्धारित किया जाता है) सॉल्वर को ऊर्जा बाधाओं (energy barriers) के ऊपर कूदने से रोकते हैं, जिससे यह सुनिश्चित होता है कि सॉल्वर मनमाने वैश्विक न्यूनतम (global minima) के बजाय भौतिक रूप से सार्थक मेटास्टेबल इक्विलिब्रियम का चयन करे।
- प्रीकंडीशनिंग: अभिसरण (convergence) को तेज करने के लिए एक भौतिक-आधारित प्रीकंडीशनर लगाया जाता है, जो एक समतुल्य होमोजेनियस समस्या के इनवर्स हेसियन (inverse Hessian) से प्राप्त होता है, विशेष रूप से कम विषमता कंट्रास्ट के लिए।
- सॉफ्टवेयर: इसे petsc4py और ऑटोमैटिक डिफरेंशिएशन एवं जस्ट-इन-टाइम कंपाइलेशन के लिए JAX का उपयोग करके पायथन में कार्यान्वित किया गया है।

मुख्य योगदान

सैद्धांतिक व्युत्पत्ति: मिश्रित-मोड I+II+III लोडिंग के तहत एक अर्ध-वृत्ताकार दरार के लिए पोटेंशियल एनर्जी का एक कठोर एसिम्प्टोटिक एक्सपेंशन, जो एक क्लोज्ड-फॉर्म फॉर्मूला प्रदान करता है जो वेरिएशनल सिद्धांतों को परटर्बेशन थ्योरी से जोड़ता है।
एल्गोरिद्मिक फ्रेमवर्क: एक सुदृढ़, मैट्रिक्स-फ्री सॉल्वर जो भौतिक बाधाओं (अपरिवर्तनीयता, ऊर्जा बाधाएं, लॉन्ग-रेंज इंटरैक्शन) को लागू करते हुए भी गणनात्मक दक्षता ( $O(N \log N)$ जटिलता प्रति इटरेशन) बनाए रखता है।
वैलिडेशन: कार्यान्वयन को होमोजेनियस और विषम माध्यमों में शियर क्रैक्स के लिए नए प्राप्त विश्लेषणात्मक समाधानों के विरुद्ध मान्य किया गया है, जो फ्रंट विरूपण और फूरियर मोड विकास की भविष्यवाणी करने में सटीकता प्रदर्शित करता है।

परिणाम
लेखकों ने विषमता तीव्रता, परिमित-आकार प्रभावों और मोड मिक्सिटी के प्रभाव को मापने के लिए विस्कस मीडिया में टेंसाइल (मोड I) और शियर (मोड II+III) दरार प्रसार के 116,000 बड़े पैमाने के सिमुलेशन किए।

इंटरमिटेंसी की ओर संक्रमण: सिमुलेशन निरंतर दरार वृद्धि से इंटरमिटेंट ग्रोथ (पिंनिंग चरणों और एवलांच के बीच बारी-बारी से होने वाली वृद्धि) में संक्रमण को सफलतापूर्वक दर्शाते हैं। यह संक्रमण दरार के परिधि और लार्किन लेंथ (Larkin length) के अनुपात द्वारा नियंत्रित होता है, जो विषमता तीव्रता और विषमता स्केल पर निर्भर करता है।
मोड मिक्सिटी प्रभाव:
- मोड मिक्सिटी का इंटरमिटेंसी की शुरुआत पर सीमित प्रभाव पड़ता है।
- हालाँकि, गैर-शून्य पॉइसन रेशियो ( $\nu \neq 0$ ) के साथ मिश्रित मोड II+III लोडिंग में, दरार का अग्रभाग एक क्वासी-इलिप्सॉइडल (quasi-elliptic) आकार विकसित करता है, जहाँ मोड III सेक्टर, मोड II सेक्टर की तुलना में धीमी गति से आगे बढ़ते हैं।
ऊर्जा अपव्यय और रेडिएटेड एनर्जी: सिमुलेशन दर्शाते हैं कि इंटरमिटेंट डायनेमिक्स ऊर्जा अपव्यय और पोटेंशियल एनर्जी में मल्टीस्केल उतार-चढ़ाव की ओर ले जाते हैं। पोटेंशियल एनर्जी में गिरावट और डिसिपेटेड एनर्जी के स्पाइक्स के बीच का अंतर रेडिएटेड एनर्जी के प्रॉक्सी के रूप में कार्य करता है, जो भूकंपीय मोमेंट रिलीज (seismic moment release) के समान बर्स्ट-लाइक घटनाओं को दिखाता है।
आकार-निर्भर टफनिंग/वीकनिंग:
- छोटी दरारें: विषमताएं स्पष्ट वीकनिंग ( $\Delta S > 0$ ) प्रेरित करती हैं क्योंकि दरार का अग्रभाग लो-टफनेस क्षेत्रों में झुक जाता है।
- बड़ी दरारें: एक स्ट्रॉन्ग पिंनिंग रिजीम में प्रवेश करने पर टफनिंग ( $\Delta S < 0$ ) की ओर एक क्रॉसओवर होता है, जहाँ अस्थिरताएं दरार को प्राथमिकता से हाई-टफनेस क्षेत्रों को सैंपल करने के लिए मजबूर करती हैं।
- इस प्रभाव का परिमाण विषमता तीव्रता ( $\sigma/G_c^0$ ) के वर्ग के साथ स्केल करता है।

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

लेखक इस बात पर जोर देते हैं कि उनका दृष्टिकोण दरार की स्थिरता का सीधा आकलन और परिमित-आकार प्रभावों तथा मोड मिक्सिटी का मात्रात्मक मूल्यांकन करने की अनुमति देता है, जो पहले 3D विषम फ्रैक्चर में पूरी तरह से समझ में नहीं आए थे। परिणाम सूक्ष्म स्तर की पैटर्न संरचना (microscale patterning) के माध्यम से डैमेज-टोलरेंट सामग्रियों के डिजाइन के लिए एक तंत्र को उजागर करते हैं, जहाँ वीकनिंग से टफनिंग के संक्रमण को दरार के सापेक्ष आकार और विषमता वितरण के माध्यम से नियंत्रित किया जा सकता है। इस फ्रेमवर्क को अन्य ज्यामिति (सेमी-इनफिनिट, टनल क्रैक्स) के अनुकूल और संभावित रूप से डायनेमिक रप्चर और कोहेसिव फ्रैक्चर के लिए विस्तार योग्य प्रस्तुत किया गया है।

1. समस्या: "बहुत धीमा" बनाम "बहुत सरल" की दुविधा

2. यह कैसे काम करता है: "ऊर्जा न्यूनीकरण" (Energy Minimization) का खेल

3. उन्होंने क्या खोजा: "116,000 सिमुलेशन"

4. यह क्यों महत्वपूर्ण है (शोध पत्र के अनुसार)

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