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

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

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

$$\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 $$ अब समीकरण के दाएँ पक्ष का मूल्यांकन करने के लिए हम एक वेटिंग पैरामीटर का उपयोग करते हैं $$ \theta $$ 0 और 1 के बीच, और हम का एकीकरण लिखते हैं $$ 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 $$ और समाधान को एक और समय चरण तक आगे बढ़ाने के लिए प्रक्रिया को दोहराया जाता है। यह देखा जा सकता है कि सभी गुणांक सकारात्मक हैं, जो समय के किसी भी आकार के लिए अंतर्निहित योजना को बिना शर्त स्थिर बनाता है। चूंकि योजना की सटीकता समय में केवल प्रथम-क्रम है, इसलिए परिणामों की सटीकता सुनिश्चित करने के लिए छोटे समय के कदमों की आवश्यकता होती है। इसकी मजबूती और बिना शर्त स्थिरता के कारण सामान्य प्रयोजन क्षणिक गणना के लिए अंतर्निहित विधि की सिफारिश की जाती है