Predicting cement microstructure and mechanical properties in hydrating cement paste with a Phase-Field model

यहाँ शोध पत्र "Predicting cement microstructure and mechanical properties in hydrating cement paste with a Phase-Field model" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

सीमेंट-आधारित सामग्रियों के यांत्रिक गुणों (mechanical properties) का विकास हाइड्रेशन के दौरान उनकी सूक्ष्म संरचना (microstructure) के विकास से आंतरिक रूप से जुड़ा हुआ है। जबकि अनुभवजन्य (empirical) और अर्ध-अनुभवजन्य (semi-empirical) मॉडल मौजूद हैं, वे अक्सर आदर्श सूक्ष्म संरचनाओं या सरलीकृत प्रतिक्रिया गतिकी (reaction kinetics) पर निर्भर करते हैं जो सूक्ष्म स्तर पर विलेयता (dissolution) और अवक्षेपण (precipitation) की जटिल ऊष्मागतिकी (thermodynamics) को पकड़ने में विफल रहते हैं।

मौजूदा भौतिक-आधारित दृष्टिकोण, जैसे कि सेलुलर ऑटोमेटा (CA) मॉडल (जैसे, CEMHYD3D), व्यापक रूप से उपयोग किए जाते हैं लेकिन वे निम्नलिखित सीमाओं से ग्रस्त हैं:

• वोक्सेल रिज़ॉल्यूशन (Voxel Resolution): वे ग्रिड रिज़ॉल्यूशन द्वारा बाधित होते हैं, जिससे "स्टेयर-स्टेप" (सीढ़ीनुमा) इंटरफेस उत्पन्न होते हैं।
• भौतिक विसंगतियाँ (Physical Inconsistencies): सीमेंट हाइड्रेशन के पिछले फेज-फील्ड (PF) फॉर्मूलेशन में ऊष्मागतिक दोष थे, विशेष रूप से:
  1. साम्यावस्था (equilibrium) पर भी स्वतः जेल अवक्षेपण (gel precipitation) की अनुमति देना।
  2. यह मानना कि स्रोत विलेयता (source dissolution) और जेल अवक्षेपण (gel precipitation) के लिए समान साम्यावस्था सांद्रता होती है, जबकि स्थापित हाइड्रेशन सिद्धांत निर्देश देता है कि ये मान भिन्न होते हैं।

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

2. कार्यप्रणाली (Methodology)

A. अनुकूलित फेज-फील्ड (PF) फॉर्मूलेशन

लेखक पूर्व साहित्य में पाई गई ऊष्मागतिक विसंगतियों को हल करने के लिए एक संशोधित PF मॉडल प्रस्तावित करते हैं। यह मॉडल तीन प्राथमिक चरणों (phases) को ट्रैक करता है:

1. स्रोत ($C_3S$): अनहाइड्रस सीमेंट क्लिंकर।
2. पेस्ट ($C\text{-}S\text{-}H$): हाइड्रेशन उत्पाद (कैल्शियम सिलिकेट हाइड्रेट जेल)।
3. विलेय ($c$): पोर वॉटर (pore water) में घुले हुए प्रजाति।

प्रमुख नवाचार:

• संशोधित मुक्त-ऊर्जा क्षमता (Revised Free-Energy Potential): एक 'डबल-वेल पोटेंशियल' का उपयोग किया जाता है, लेकिन इसे स्थानीय सांद्रता विचलन के आधार पर प्रतिक्रिया को संचालित करने के लिए विशिष्ट स्रोत और सिंक टर्म्स ($E_d$ विलेयता के लिए, $E_p$ अवक्षेपण के लिए) का उपयोग करके "टिल्ट" (tilt) किया गया है।
• विशिष्ट साम्यावस्था स्थिरांक (Distinct Equilibrium Constants): मॉडल स्रोत ($c_{eq, C3S}$) और जेल ($c_{eq, CSH}$) के लिए अलग-अलग साम्यावस्था सांद्रता पेश करता है। यह साम्यावस्था पर स्वतः अवक्षेपण को रोकता है और न्यूक्लिएशन (nucleation) के लिए आवश्यक सुपरसैचुरेशन (supersaturation) को सही ढंग से मॉडल करता है।
• द्रव्यमान संरक्षण (Mass Conservation): युग्मित एलेन-कैन (Allen-Cahn) समीकरण चरण विकास को नियंत्रित करते हैं, जबकि एक विसरण (diffusion) समीकरण विलेय परिवहन को ट्रैक करता है। द्रव्यमान संरक्षण को संतृप्ति सीमाओं (saturation limits) से प्राप्त रूपांतरण गुणांकों ($\alpha$) के माध्यम से लागू किया जाता है।
• डिफ्यूसिविटी मॉडलिंग (Diffusivity Modeling): आयनों के भौतिक ट्रैपिंग को दर्शाने के लिए जेल चरण के भीतर विलेय प्रसारशीलता (solute diffusivity) को कम किया जाता है, जो हाइड्रेशन गतिकी को धीमा करने में एक महत्वपूर्ण कारक है।

B. संख्यात्मक कार्यान्वयन (Numerical Implementation)

• सॉल्वर (Solver): युग्मित आंशिक अंतर समीकरणों (coupled partial differential equations) के सिस्टम को MOOSE परिमित तत्व ढांचे (finite element framework) का उपयोग करके हल किया जाता है।
• डोमेन (Domain): सिमुलेशन 2D आवधिक डोमेन ($100 \mu m$) पर $0.2 \mu m$ मेश आकार के साथ किए जाते हैं।
• परिदृश्य (Scenarios): दो जल-से-सीमेंट ($w/c$) अनुपात परीक्षण किए गए हैं: $0.3$(अति-संतृप्त) और$0.5$ (असंतृप्त)।
• प्रारंभिक स्थितियाँ (Initial Conditions): एक संशोधित कण आकार वितरण (PF स्थिरता सीमाओं के लिए समायोजित) के आधार पर गोलाकार स्रोत कण उत्पन्न किए जाते हैं, जिसमें जेल वृद्धि शुरू करने के लिए स्रोत के पास यादृच्छिक न्यूक्लिएशन बीज (nucleation seeds) रखे जाते हैं।

C. कम्प्यूटेशनल होमोजेनाइजेशन (Computational Homogenization)

एक बार जब सूक्ष्म संरचना एक विशिष्ट हाइड्रेशन डिग्री तक विकसित हो जाती है, तो मॉडल मैकेनिकल विश्लेषण में बदल जाता है:

• लोडिंग (Loading): सूक्ष्म संरचना पर तन्यता (tensile) और कर्तन (shear) सीमा स्थितियाँ लागू की जाती हैं।
• सामग्री गुण (Material Properties): चरणों को आइसोट्रोपिक इलास्टिक गुणों (यंग का मापांक और पॉइसन अनुपात) को सौंपा जाता है। पोर चरण को तरल व्यवहार का अनुकरण करने के लिए लगभग शून्य मापांक दिया जाता है।
• आउटपुट: होमोजेनाइज्ड यंग का मापांक ($E_h$) और शियर मापांक ($G_h$) स्ट्रेस-स्ट्रेन औसत के माध्यम से गणना किए जाते हैं।

3. प्रमुख योगदान (Key Contributions)

1. ऊष्मागतिक रूप से सुसंगत PF मॉडल: शोध पत्र एक संशोधित मुक्त-ऊर्जा परिदृश्य और विशिष्ट साम्यावस्था स्थिरांक पेश करता है जो पिछले PF मॉडलों की भौतिक विसंगतियों को ठीक करता है, जिससे यथार्थवादी विलेयता/अवक्षेपण व्यवहार सुनिश्चित होता है।
2. स्पष्ट सूक्ष्म संरचना-मैकेनिक्स लिंक: ढांचा सफलतापूर्वक कम्प्यूटेशनल होमोजेनाइजेशन के माध्यम से सूक्ष्म संरचनात्मक विकास को सीधे यांत्रिक गुण भविष्यवाणी से जोड़ता है, जिससे मध्यवर्ती अनुभवजन्य सहसंबंधों की आवश्यकता समाप्त हो जाती है।
3. साहित्य के विरुद्ध सत्यापन: मॉडल का कठोरता से निम्नलिखित के विरुद्ध सत्यापन किया गया है:
   • प्रयोगात्मक डेटा (हाइड्रेशन कर्व्स, सरंध्रता/porosity)।
   • व्यापक रूप से उपयोग किया जाने वाला CEMHYD3D (CA) मॉडल।
   • स्कैनिंग इलेक्ट्रॉन माइक्रोस्कोप (SEM) अवलोकन।
4. मात्रात्मक रूपात्मक विश्लेषण (Quantitative Morphological Analysis): अध्ययन उन्नत रूपात्मक डिस्क्रिप्टर्स (two-point correlation functions, perimeter evolution, feature counting) का उपयोग यह प्रदर्शित करने के लिए करता है कि PF मॉडल CA मॉडलों की तुलना में पृथक समुच्चय (isolated aggregates) से एक निरंतर नेटवर्क (percolating network) में संक्रमण को अधिक सटीक रूप से कैप्चर करता है।

4. परिणाम (Results)

सूक्ष्म संरचनात्मक विकास (Microstructural Evolution)

• गतिकी (Kinetics): मॉडल विशिष्ट दो-चरणीय हाइड्रेशन वक्र (तेज प्रारंभिक विलेयता के बाद प्रसार-नियंत्रित धीमी अवस्था) को पुनरुत्पादित करता है।
• चरण वितरण (Phase Distribution): PF मॉडल CEMHYD3D की तुलना में कम सरंध्रता (porosity) और उच्च पेस्ट वॉल्यूम अंश की भविष्यवाणी करता है। पाया गया कि CA मॉडल छिद्र स्थान (pore space) को बढ़ा-चढ़ाकर दिखाता है, जिससे कृत्रिम रूप से कमजोर सूक्ष्म संरचना बनती है।
• इंटरफेस गुणवत्ता (Interface Quality): CA मॉडलों के "स्टेयर-स्टेप" इंटरफेस के विपरीत, PF मॉडल चिकने, निरंतर चरण सीमाओं को उत्पन्न करता है जो भौतिक वास्तविकता को बेहतर ढंग से दर्शाते हैं।
• परकोलेशन (Percolation): मॉडल उस टोपोलॉजिकल संक्रमण को सफलतापूर्वक कैप्चर करता है जहाँ अलग-थलग जेल द्वीप एक निरंतर, भार-वहन करने वाले नेटवर्क में विलीन हो जाते हैं।

यांत्रिक गुण (Mechanical Properties)

• कठोरता विकास (Stiffness Evolution): यंग और शियर दोनों मापांक हाइड्रेशन के साथ बढ़ते हैं। PF-अनुमानित मापांक प्रयोगात्मक डेटा और वास्तविक SEM-आधारित होमोजेनाइजेशन परिणामों के करीब संरेखित होते हैं।
• CA के साथ तुलना: CEMHYD3D-व्युत्पन्न सूक्ष्म संरचनाओं ने सरंध्रता और विच्छिन्न चरण सीमाओं को बढ़ा-चढ़ाकर दिखाने के कारण काफी कम यांत्रिक मापांक दिए। PF मॉडल कठोरता का अधिक भौतिक रूप से सुसंगत अनुमान प्रदान करता है।
• $w/c$ अनुपात का प्रभाव: $w/c = 0.3$ मिश्रण ने $w/c = 0.5$ की तुलना में उच्च कठोरता दिखाई क्योंकि इसमें अनहाइड्रस, कठोर स्रोत समावेशन (inclusions) मौजूद थे।

5. महत्व और निष्कर्ष (Significance and Conclusion)

यह कार्य सीमेंट हाइड्रेशन के लिए एक भौतिक-आधारित डिजिटल प्रयोगशाला स्थापित करता है। फेज-फील्ड मॉडलिंग में ऊष्मागतिक विसंगतियों को सुधारकर, लेखक एक ऐसा उपकरण प्रदान करते हैं जो:

• भविष्यवाणी सटीकता में सुधार करता है: यह अर्ध-अनुभवजन्य या CA-आधारित तरीकों की तुलना में रासायनिक संरचना, सूक्ष्म संरचना और यांत्रिक प्रदर्शन के बीच अधिक विश्वसनीय लिंक प्रदान करता है।
• मल्टीफिजिक्स विस्तार को सक्षम बनाता है: इस ढांचे को तनाव-निर्भर हाइड्रेशन, थर्मल प्रभाव और तरल परिवहन जैसे युग्मित घटनाओं को शामिल करने के लिए विस्तारित किया जा सकता है।
• अनुभववाद पर निर्भरता को कम करता है: यह प्रथम-सिद्धांत सिमुलेशन (first-principles simulation) की ओर क्षेत्र को ले जाता है, जिससे कम पर्यावरणीय प्रभाव के साथ सीमेंट फॉर्मूलेशन का अनुकूलन संभव होता है।

हालांकि वर्तमान कार्यान्वयन 2D है और गणनात्मक रूप से गहन है (~10 घंटे प्रति सिमुलेशन), लेखक तर्क देते हैं कि यह भविष्य के 3D विस्तार और अगली पीढ़ी की सीमेंट सामग्री के विकास के लिए एक आवश्यक आधार प्रदान करता है। सिम्युलेटेड सूक्ष्म संरचना विकास से सीधे यांत्रिक व्यवहार की भविष्यवाणी करने की क्षमता यह सुझाव देती है कि यह दृष्टिकोण पारंपरिक 28-दिवसीय परीक्षण प्रोटोकॉल के पूरक या अंततः उसका विकल्प बन सकता है।