अस्थिर प्रवाह के लिए परिमित आयतन विधि

अस्थिर प्रवाह को ऐसे प्रवाह के रूप में जाना जाता है जिसमें तरल पदार्थ के गुण समय पर निर्भर होते हैं। यह समीकरण संचालन में प्रतिबिंबित होता है क्योंकि गुणों का अवकलज समय अनुपस्थित है। अस्थिर प्रवाह के लिए परिमित-मात्रा विधि का अध्ययन करने के लिए कुछ नियामक समीकरण हैं >

समीकरण संचालन
अस्थिर प्रवाह में अदिश के परिवहन के लिए संरक्षण समीकरण का सामान्य रूप इस प्रकार है

$$\frac{\partial \rho \phi }{\partial t} + \operatorname{div}\left(\rho \phi \upsilon\right) = \operatorname{div}\left(\Gamma \operatorname{grad} \phi\right) + S_\phi$$

$$\rho$$ घनत्व है और $$ \phi $$ सभी द्रव प्रवाह का अपरिवर्तनवादी रूप है,

$$\Gamma$$ प्रसार गुणांक है और $$S$$ स्रोत पद है। $$\operatorname{div}\left(\rho \phi \upsilon\right)$$तरल पदार्थ तत्व (संवहन) से $$ \phi $$ के प्रवाह की परिष्कृत दर है, $$\operatorname{div}\left(\Gamma \operatorname{grad} \phi\right) $$ की वृद्धि दर है $$ \phi $$ प्रसार के कारण,

$$ S_\phi$$ स्रोतों के कारण $$\phi$$ की वृद्धि की दर है। $$\frac{\partial \rho \phi }{\partial t} $$ द्रव तत्व के $$ \phi $$ की वृद्धि की दर (क्षणिक) है,

समीकरण का पहला पद प्रवाह की अस्थिरता को दर्शाता है और स्थिर प्रवाह के मामले में अनुपस्थित है। समीकरण संचालन का परिमित आयतन एकीकरण एक नियंत्रण आयतन और एक सीमित समय चरण ∆t पर भी किया जाता है।

$$\int\limits_{cv} \!\!\!\int_t^ {t+\Delta t} \left(\frac{\partial \rho \phi }{\partial t} \,\mathrm{d}t\right)\,\mathrm{d}V + \int_t^ {t+\Delta t}\!\!\!\int\limits_A \left(n.{\rho \phi u} \,\mathrm{d}A\right)\,\mathrm{d}t = \int_t^ {t+\Delta t}\!\!\!\int\limits_A \left(n \cdot \left(\Gamma \operatorname{grad} \phi\right)\,\mathrm{d}A\right)\,\mathrm{d}t +\int_t^ {t+\Delta t} \!\!\!\int\limits_{cv} S_\phi\,\mathrm{d}V\,\mathrm{d}t $$

समीकरण के स्थिर भाग का नियंत्रण आयतन एकीकरण स्थिर अवस्था शासी समीकरण के एकीकरण के समान है। हमें समीकरण के अस्थिर घटक के एकीकरण पर ध्यान केंद्रित करने की आवश्यकता है। एकीकरण तकनीक का एहसास पाने के लिए, हम एक-आयामी अस्थिर ताप चालन समीकरण का संदर्भ लेते हैं।

$$ \rho c \frac{\partial T} {\partial t} = \frac{\partial \frac{ k \partial T} {\partial x}} {\partial x} + S $$

$$\int_t^ {t+\Delta t} \!\!\!\int\limits_{cv} \rho c \frac{\partial T} {\partial t}\,\mathrm{d}V\,\mathrm{d}t = \int_t^ {t+\Delta t} \!\!\!\int\limits_{cv} \frac{\partial \frac{ k \partial T} {\partial x}} {\partial x}\,\mathrm{d}V\,\mathrm{d}t + \int_t^ {t+\Delta t} \!\!\!\int\limits_{cv} S\,\mathrm{d}V\,\mathrm{d}t$$

$$\int_e^w \!\!\!\int_t^ {t+\Delta t} \left(\rho c \frac{\partial T} {\partial t}\,\mathrm{d}t\right)\,\mathrm{d}V = \int_t^ {t+\Delta t} \left[ \left(k A \frac{\partial T} {\partial x}\right)_e - \left(k A \frac{\partial T} {\partial x}\right)_w\right]\,\mathrm{d}t + \int_t^ {t+\Delta t} \bar S\Delta V \,\mathrm{d}t $$

अब, संपूर्ण नियंत्रण आयतन में प्रचलित नोड पर तापमान की धारणा को ध्यान में रखते हुए, समीकरण के बाईं ओर को के रूप में लिखा जा सकता है।

$$\int\limits_{cv} \!\!\!\int_t^ {t+\Delta t} \left(\rho c \frac{\partial T} {\partial t}\,\mathrm{d}t\right)\,\mathrm{d}V = \rho c\left(T_P - {T_P}^O\right) \Delta V $$

पहले क्रम की पश्चगामी अवकलन योजना का उपयोग करके, हम समीकरण के दाहिने हाथ को इस प्रकार लिख सकते हैं

$$ \rho c \left(T_P - {T_P}^0\right) \Delta V = \int_t^{t+\Delta t} \left[\left( K_e A \frac {T_E - T_P} {\delta x_{PE}}\right) - \left( K_w A \frac {T_P - T_W} { \delta x_{WP}}\right)\right] \,\mathrm{d}t + \int_t^{t+\Delta t} \bar S\Delta V \,\mathrm{d}t $$

अब समीकरण के दाहिने पक्ष का मूल्यांकन करने के लिए हम 0 और 1 के बीच एक वेटिंग पैरामीटर $$ \theta $$ का उपयोग करते हैं, और हम $$ T_P $$ का एकीकरण लिखते हैं।

$$ I_T = \int_t^{t+\Delta t} T_P \,\mathrm{d}t = \left[ \theta T_P - \left(1 - \theta \right) {T_P}^0 \right] \Delta t $$

अब, अंतिम पृथक समीकरण का सटीक रूप $$ \Theta $$ के मूल्य पर निर्भर करता है। चूंकि $$ \Theta $$ का विचरण 0< $$ \Theta $$<1 है, $$ T_P $$ की गणना करने के लिए उपयोग की जाने वाली योजना $$ \Theta $$ के मान पर निर्भर करती है।

विभिन्न योजनाएँ
1. स्पष्ट योजना स्पष्ट योजना में स्रोत शब्द को इस प्रकार रैखिक किया गया है $$ b = S_u + {S_P}{T_P}^0 $$. हम स्थानापन्न करते हैं $$ \theta = 0 $$ स्पष्ट विवेक प्राप्त करने के लिए अर्थात:

$$ a_P T_P = a_w {T_w}^0 + a_e {T_e}^0 + \left[ {a_P}^0 - \left( a_w + a_e - S_P \right)\right] {T_P}^0 + S_u $$ कहाँ $$ a_P = {a_P}^0 $$. ध्यान देने योग्य एक बात यह है कि दाईं ओर पुराने समय के चरण के मान शामिल हैं और इसलिए बाईं ओर की गणना समय में आगे मिलान करके की जा सकती है। यह योजना बैकवर्ड डिफरेंसिंग पर आधारित है और इसकी टेलर श्रृंखला ट्रंकेशन त्रुटि समय के संबंध में पहले क्रम की है। सभी गुणांक सकारात्मक होने चाहिए. निरंतर k और समान ग्रिड रिक्ति के लिए, $$ \delta x_{PE} = \delta x_{WP} = \Delta x $$ इस शर्त को इस प्रकार लिखा जा सकता है

$$ \rho c \frac { \Delta x } { \Delta t } > \frac {2K} { \Delta x } $$ यह असमानता अधिकतम समय कदम पर एक कठोर शर्त निर्धारित करती है जिसका उपयोग किया जा सकता है और योजना पर एक गंभीर सीमा का प्रतिनिधित्व करता है। स्थानिक सटीकता में सुधार करना बहुत महंगा हो जाता है क्योंकि अधिकतम संभव समय चरण को वर्ग के रूप में कम करने की आवश्यकता होती है $$ \Delta x $$ 2. क्रैंक-निकोलसन योजना: क्रैंक-निकोलसन विधि सेटिंग से उत्पन्न होती है $$ \theta = \frac {1}{2}$$. विवेचित अस्थिर ऊष्मा चालन समीकरण बन जाता है

$$ a_P T_P = a_E \left[ \frac {T_E + {T_E}^0} {2}\right] + a_W \left[ \frac {T_W + {T_W}^0} {2}\right] + \left[ {a_P}^0 - \frac {a_E} {2} - \frac {a_W} {2}\right] {T_P}^0 + b $$ कहाँ $$ a_P = \frac {a_W + a_E} {2} + {a_P}^0 - \frac {S_P} {2} $$ चूंकि नए समय स्तर पर टी के एक से अधिक अज्ञात मान समीकरण में मौजूद हैं, इसलिए विधि अंतर्निहित है और प्रत्येक समय चरण पर सभी नोड बिंदुओं के लिए एक साथ समीकरणों को हल करने की आवश्यकता है। हालाँकि योजनाओं के साथ $$ \frac {1}{2} < \theta < 1 $$ क्रैंक-निकोलसन योजना सहित, समय चरण के सभी मूल्यों के लिए बिना शर्त स्थिर हैं, यह सुनिश्चित करना अधिक महत्वपूर्ण है कि सभी गुणांक शारीरिक रूप से यथार्थवादी और सीमित परिणामों के लिए सकारात्मक हैं। यह मामला है यदि का गुणांक $$ {T_P}^0$$ निम्नलिखित शर्त को पूरा करता है

$$ {a_P}^0 = \left[ \frac {a_E + a_W} {2} \right]$$ जिससे होता है

$$ \Delta t < \rho c \frac { \Delta x^2} {K} $$ क्रैंक-निकोलसन केंद्रीय भिन्नता पर आधारित है और इसलिए समय में दूसरा क्रम सटीक है। गणना की समग्र सटीकता स्थानिक भिन्नता अभ्यास पर भी निर्भर करती है, इसलिए क्रैंक-निकोलसन योजना का उपयोग आम तौर पर स्थानिक केंद्रीय भिन्नता के साथ संयोजन में किया जाता है

3. पूरी तरह से अंतर्निहित योजना जब Ѳ का मान 1 पर सेट किया जाता है तो हमें पूरी तरह से अंतर्निहित योजना मिलती है। विच्छेदित समीकरण है:

$$ a_P T_P = a_W T_W + a_E T_E + {a_P}^0 {T_P}^0 + S_u $$

$$ a_P = {a_P}^0 + a_W + a_E - S_P $$ समीकरण के दोनों पक्षों में नए समय चरण पर तापमान होता है, और प्रत्येक समय स्तर पर बीजगणितीय समीकरणों की एक प्रणाली को हल किया जाना चाहिए। टाइम मार्चिंग प्रक्रिया तापमान के दिए गए प्रारंभिक क्षेत्र से शुरू होती है $$ T^0 $$. समय चरण का चयन करने के बाद समीकरणों की प्रणाली को हल किया जाता है $$ \Delta t $$. अगला समाधान $$ T $$ को सौंपा गया है $$ T^0 $$ और समाधान को एक और समय चरण तक आगे बढ़ाने के लिए प्रक्रिया को दोहराया जाता है। यह देखा जा सकता है कि सभी गुणांक सकारात्मक हैं, जो समय के किसी भी आकार के लिए अंतर्निहित योजना को बिना शर्त स्थिर बनाता है। चूंकि योजना की सटीकता समय में केवल प्रथम-क्रम है, इसलिए परिणामों की सटीकता सुनिश्चित करने के लिए छोटे समय के कदमों की आवश्यकता होती है। इसकी मजबूती और बिना शर्त स्थिरता के कारण सामान्य प्रयोजन क्षणिक गणना के लिए अंतर्निहित विधि की सिफारिश की जाती है