Identify dominant parameters of comprehensive pyrolysis model - Part I: Variance-based sensitivity analysis of wooden materials

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

वुहान यूनिवर्सिटी ऑफ टेक्नोलॉजी के शोधकर्ताओं की एक टीम ने लकड़ी के जलने के अध्ययन पर एक कठोर सांख्यिकीय पद्धति लागू करके इस प्रश्न का उत्तर देने का प्रयास किया। उन्होंने दो सामान्य सामग्रियों पर ध्यान केंद्रित किया: बिर्च की लकड़ी और मध्यम घनत्व वाला फाइबरबोर्ड, जो एक प्रकार का इंजीनियर लकड़ी उत्पाद है। यह अनुमान लगाने के बजाय कि कौन से कारक महत्वपूर्ण थे, उन्होंने 'वैरिएंस-बेस्ड सेंसिटिविटी एनालिसिस' (variance-based sensitivity analysis) नामक एक तकनीक का उपयोग किया। कल्पना कीजिए कि एक मॉडल एक 'ब्लैक बॉक्स' है जो कई अलग-अलग सेटिंग्स लेता है और एक आउटपुट देता है, जैसे कि लकड़ी के वजन घटने की दर या उसके केंद्र का तापमान। शोधकर्ताओं ने व्यवस्थित रूप से प्रत्येक इनपुट पैरामीटर को एक वास्तविक सीमा के भीतर बदला और देखा कि आउटपुट में कितना परिवर्तन हुआ। यह मापकर कि आउटपुट में अनिश्चितता को प्रत्येक विशिष्ट इनपुट के माध्यम से कितना ट्रैक किया जा सकता था, वे इस प्रक्रिया के प्रमुख चालकों की पहचान कर सके। उनका लक्ष्य आवश्यक भौतिकी को शोर (noise) से अलग करना था, यह निर्धारित करना था कि रासायनिक पथों के सटीक विवरण या ऊष्मा हस्तांतरण के व्यापक पहलू पायरोलिसिस प्रक्रिया के वास्तविक स्वामी हैं।

जांच ने लकड़ी के रासायनिक गुणों के बीच एक उल्लेखनीय पदानुक्रम का खुलासा किया। शोधकर्ताओं ने पाया कि सक्रियण ऊर्जा (activation energy), जो रासायनिक प्रतिक्रिया शुरू करने के लिए आवश्यक तापमान सीमा का एक माप है, अत्यधिक प्रभावशाली कारक थी। उनके सिमुलेशन में, इस तापमान सीमा की संवेदनशीलता प्रतिक्रिया क्रम (reaction order) की तुलना में, जो कि एक पैरामीटर है जो रासायनिक प्रतिक्रिया के विशिष्ट आकार या पथ का वर्णन करता है, सामान्यतः दस गुना से अधिक थी। यह निष्कर्ष बताता है कि लकड़ी के जलने की भविष्यवाणी करने के उद्देश्य से, वह सटीक तापमान सीमा जिसमें सामग्री टूटती है, प्रतिक्रिया पथ के सूक्ष्म विवरणों की तुलना में कहीं अधिक महत्वपूर्ण है। जबकि प्रतिक्रिया क्रम द्रव्यमान हानि (mass loss) दिखाने वाले वक्र के आकार को थोड़ा बदल सकता है, लेकिन यह समग्र परिणाम को महत्वपूर्ण रूप से नहीं बदलता है। इसका तात्पर्य यह है कि मॉडलर्स को प्रतिक्रिया क्रमों के लिए पूर्ण मानों के पीछे अत्यधिक समय व्यर्थ करने की आवश्यकता नहीं है; सही तापमान सीमा पर ध्यान केंद्रित करने से सटीकता के लिए बहुत बेहतर प्रतिफल (return on investment) प्राप्त होता है।

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

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

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

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