🔬 materials science

A complete phase-field fracture model for brittle materials subjected to thermal shocks

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

मूल लेखक: Bo Zeng, John E. Dolbow

प्रकाशित 2026-06-19

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

CC BY 4.0

मूल लेखक: Bo Zeng, John E. Dolbow

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

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

लेकिन उस पुराने जासूस ने एक बहुत बड़ा सुराग मिस कर दिया: एक खरोंच होने से पहले ही वह सामग्री वास्तव में कितनी मजबूत है।

यह शोध पत्र इन टूट-फूटों को मॉडल करने का एक नया, "पूर्ण" तरीका पेश करता है। इसे एक ऐसे जासूस में अपग्रेड करने के रूप में सोचें जो केवल दरारों को देखने के बजाय सामग्री के पूरे व्यक्तित्व को समझता है।

नए मॉडल के तीन घटक

लेखकों का तर्क है कि यह भविष्यवाणी करने के लिए कि कोई चीज़ कब टूटती है, आपको तीन चीजों को अलग-अलग, स्वतंत्र घटकों के रूप में मानना होगा, जैसे कि एक केक बनाना जहाँ आप मैदा, चीनी और अंडे को स्वतंत्र रूप से समायोजित कर सकते हैं:

कठोरता (Elasticity/Stiffness): सामग्री कितनी खिंचती या दबती है?
सहनशीलता (Toughness/Fracture Resistance): एक बार दरार शुरू होने के बाद, उसे लंबा करने में कितनी कठिनाई होती है?
शक्ति (Strength/The Breaking Point): दरार शुरू होने से पहले सामग्री कितना तनाव झेल सकती है?

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

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

तीन परीक्षण (बेंचमार्क)

यह साबित करने के लिए कि उनका नया मॉडल काम करता है, लेखकों ने इसे तीन वास्तविक दुनिया के परिदृश्यों के विरुद्ध परखा, जैसे कि एक नई कार के इंजन को तीन अलग-अलग प्रकार के रास्तों पर चलाना:

1. ग्लास प्लेट रेस ("क्वेंचिंग" टेस्ट)

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

2. सिरेमिक डिस्क ("हीट लैंप" टेस्ट)

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

3. न्यूक्लियर फ्यूल पेलेट ("पावर सर्ज" टेस्ट)

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

मुख्य निष्कर्ष

शोध पत्र निष्कर्ष निकालता है कि अत्यधिक गर्मी के तहत भंगुर सामग्रियों (जैसे कांच, सिरेमिक या परमाणु ईंधन) के टूटने की भविष्यवाणी करने के लिए, आप केवल यह नहीं देख सकते कि वे कितने कठोर या कितने लचीले हैं। आपको कंप्यूटर को स्पष्ट रूप से यह भी बताना होगा कि वे कितने मजबूत हैं।

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

तकनीकी सारांश: थर्मल शॉक के अधीन भंगुर सामग्रियों के लिए एक पूर्ण फेज़-फील्ड फ्रैक्चर मॉडल

समस्या विवरण
थर्मल शॉक (तापीय आघात) के तहत भंगुर सामग्रियों के व्यवहार की भविष्यवाणी करना एक महत्वपूर्ण वैज्ञानिक चुनौती बनी हुई है। लीनियर इलास्टिक फ्रैक्चर मैकेनिक्स (LEFM) पर आधारित पारंपरिक दृष्टिकोण पूर्व-विद्यमान दरारों (pre-existing cracks) पर निर्भर करते हैं और दरार के उद्भव (nucleation) या ब्रांचिंग जैसे जटिल विकसित होते पथों का वर्णन करने में संघर्ष करते हैं। जबकि फेज़-फील्ड मॉडल ने तीक्ष्ण दरारों को विसरित क्षेत्रों (diffuse fields) से बदलकर एक सामान्य विकल्प प्रदान किया है, शास्त्रीय वैरिएशनल फेज़-फील्ड फॉर्मूलेशन में एक मौलिक भौतिक सीमा है: वे सामग्री की मजबूती (strength) को एक स्वतंत्र पदार्थ गुण (constitutive property) के रूप में नहीं देखते हैं। इन शास्त्रीय मॉडलों में, एकअक्षीय तनन शक्ति (uniaxial tensile strength) एक स्वतंत्र इनपुट नहीं है, बल्कि यह फ्रैक्चर टफनेस ( $G_c$ ), यंग्स मॉड्यूलस ( $E$ ), और रेगुलराइजेशन लंबाई ( $\ell$ ) के साथ $\sigma_{ts} \sim \sqrt{G_c E / \ell}$ के रूप में स्केल होती है। फलस्वरूप, ये मॉडल स्थानिक रूप से समान प्रतिबल अवस्थाओं (spatially uniform stress states) के तहत बल्क में दरार के उद्भव को सही ढंग से नहीं बता सकते, न ही शक्ति और ऊर्जा (energetics) के बीच परस्पर क्रिया द्वारा नियंत्रित विशिष्ट फ्रैक्चर शासनों (fracture regimes) के बीच संक्रमण को पकड़ सकते हैं।

कार्यप्रणाली
यह कार्य "पूर्ण" फेज़-फील्ड फ्रैक्चर मॉडल के विस्तार के रूप में है—जिसे मूल रूप से कुमार, फ्रांकोर्ट और लोपेज़-पामिएस द्वारा विशुद्ध रूप से यांत्रिक समस्याओं के लिए विकसित किया गया था—इसे युग्मित थर्मो-मैकेनिकल परिवेशों के लिए विस्तारित किया गया है। मुख्य नवाचार तीन स्थूल भौतिक गुणों का स्वतंत्र विनिर्देश है: प्रत्यास्थता (elasticity), फ्रैक्चर टफनेस, और सामग्री की मजबूती।

शासी समीकरण (Governing Equations): मॉडल एक 'वन-वे' थर्मल कपलिंग का उपयोग करता है जहाँ तापमान क्षेत्र तापीय प्रसार और यांत्रिक प्रतिबल को संचालित करता है, लेकिन यांत्रिक विरूपण ताप समीकरण में फीडबैक नहीं देता है। तापीय समस्या मानक ऊष्मा समीकरण द्वारा शासित होती है। यांत्रिक समस्या एक डिग्रेडेड कॉची स्ट्रेस टेंसर (degraded Cauchy stress tensor) का उपयोग करती है जो एक द्विघातीय डिग्रेडेशन फंक्शन ( $g(v)=v^2$ ) पर आधारित है।
पूर्ण मॉडल फॉर्मूलेशन: शास्त्रीय वैरिएशनल मॉडलों के विपरीत, जो एक वैश्विक ऊर्जा कार्यात्मक (global energy functional) को न्यूनतम करने से प्राप्त होते हैं, पूर्ण मॉडल फेज़-फील्ड विकास समीकरण को एक बाहरी माइक्रोफोर्स, $c_e$ $c_{e}$ के साथ संवर्धित करता है। यह माइक्रोफोर्स स्पष्ट रूप से एक सामग्री शक्ति सतह $F(\sigma) = 0$ $F (σ) = 0$ (जिसे यहाँ ड्रकर-प्रागर सतह के रूप में अपनाया गया है) द्वारा निर्धारित होता है।
- उद्भव (Nucleation): दरारें केवल वहीं उत्पन्न होती हैं जहाँ स्ट्रेंथ सरफेस (शक्ति सतह) का उल्लंघन होता है (एक आवश्यक शर्त)।
- प्रसार (Propagation): एक बार उत्पन्न होने के बाद, दरारें विरूपण और फ्रैक्चर सतह ऊर्जा के योग को न्यूनतम करने के लिए प्रसारित होती हैं (पर्याप्तता की शर्त)।
विविक्तकरण (Discretization): प्रणाली को स्थान में परिमित तत्वों (finite elements) और समय में परिमित अंतरों (finite differences) का उपयोग करके हल किया जाता है। युग्मित यांत्रिक और फेज़-फील्ड उप-समस्याओं के लिए एक स्टेगर्ड स्कीम का उपयोग किया जाता है, जिसमें एक इरेवर्सिबिलिटी कंस्ट्रेंट (अपरिवर्तनीयता बाधा) को प्राइमल-डुअल एक्टिव सेट रणनीति के माध्यम से लागू किया जाता है। विभिन्न मेश आकारों और रेगुलराइजेशन लंबाई के बीच प्रभावी फ्रैक्चर टफनेस को सुसंगत बनाए रखने के लिए पैरामीटर $\delta$ पर मेश-डिपेंडेंट सुधार लागू किए जाते हैं।

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

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

परिणाम
मॉडल की क्षमताओं का मूल्यांकन तीन प्रतिनिधि बेंचमार्क समस्याओं के माध्यम से किया गया:

ग्लास प्लेट्स का प्रोग्रेसिव क्वेंचिंग:
- मॉडल ने सीधी प्रसार से लेकर नियमित दोलन (oscillation), ब्रांचिंग और अराजक पैटर्न (chaotic patterns) तक दरार की आकृति विज्ञान (morphology) में संक्रमण को सफलतापूर्वक दोहराया, जैसे-जैसे थर्मल लोड (तापमान जंप) बढ़ता गया।
- परिणाम पेकलेट संख्या और थर्मल शॉक अनुपात द्वारा परिभाषित एक स्थापित फेज डायग्राम के अनुरूप थे।
- महत्वपूर्ण निष्कर्ष: पैरामीट्रिक अध्ययनों से पता चला कि बड़े पूर्व-विद्यमान दरारों वाले मामलों के लिए भी (जो पारंपरिक रूप से ऊर्जा-प्रधान माने जाते हैं), सामग्री की मजबूती की सतह अंतिम दरार पैटर्न को महत्वपूर्ण रूप से प्रभावित करती है। सिमुलेशन में सामग्री की मजबूती को कम करने से एक एकल दोलन करती दरार 'ट्यूनिंग-फोर्क' आकार में विभाजित हो गई, जो एक ऐसा व्यवहार है जिसे स्वतंत्र शक्ति के बिना मॉडलों द्वारा नहीं पकड़ा जा सकता।
सिरेमिक डिस्क (एलुमिना) का थर्मल शॉक:
- अध्ययन में नॉच (notch) के साथ और बिना नॉच वाली इन्फ्रारेड-हीटेड डिस्क का परीक्षण किया गया।
- नॉच्ड डिस्क: मॉडल ने डिस्क के पार नॉच टिप से सीधे दरार के प्रसार की सटीक भविष्यवाणी की।
- अछूती (Intact) डिस्क: मॉडल ने किनारे से दूर (लगभग 2-3 मिमी) दरार के उद्भव, उसके बाद सीधे प्रसार, ब्रांचिंग और मुड़ने के प्रयोगात्मक अवलोकन को सफलतापूर्वक दोहराया। यह स्थानिक रूप से परिवर्तनशील स्ट्रेंथ फील्ड और सूक्ष्म ज्यामितीय विक्षोभ (geometric perturbations) को पेश करके प्राप्त किया गया, जो दर्शाता है कि शक्ति और ऊर्जा के बीच की प्रतिस्पर्धा नॉच्ड बनाम अछूती नमूनों के बीच विशिष्ट फ्रैक्चर पैटर्न को नियंत्रित करती है।
- मॉडल ने संपीड़न क्षेत्रों (compressive regions) में कम वेग (लगभग $0.2c_R$ ) पर गतिशील दरार ब्रांचिंग को कैप्चर किया, जो तनन क्षेत्रों (tensile regions) में $0.6c_R$ पर ब्रांचिंग के पारंपरिक ज्ञान से भिन्न है।
तीव्र पावर पल्स के तहत सिरेमिक फ्यूल पेलेट्स:
- मॉडल ने क्षणिक न्यूट्रॉन और फोटॉन विकिरण के अधीन परमाणु ईंधन पेलेट्स (BeO–UO2) का अनुकरण किया।
- इसने कम ऊर्जा पल्स (200 MJ) पर अछूते पेलेट्स से उच्च पल्स (247 MJ) पर खंडित पेलेट्स में संक्रमण को सफलतापूर्वक दोहराया।
- एक सांख्यिकीय रूप से विक्षेपित स्ट्रेंथ फील्ड का उपयोग करके, मॉडल ने फ्रैक्चर परिणामों की परिवर्तनशीलता को कैप्चर किया और प्रयोगात्मक डेटा के अनुरूप दरार-कोण वितरण (crack-angle distribution) उत्पन्न किया।

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

लेखक का तर्क है कि यह फॉर्मूलेशन:

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

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

नए मॉडल के तीन घटक

तीन परीक्षण (बेंचमार्क)

मुख्य निष्कर्ष

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