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

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

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

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