फ्लक्स सीमक

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

ध्यान दें कि फ्लक्स अवरोधक को प्रवणता अवरोधक के रूप में भी जाना जाता है। क्योंकि उन दोनों का गणितीय रूप समान है, और दोनों शॉक या असंतोष के समीप समाधान प्रवणता को सीमित करने का प्रभाव होता है। अतः सामान्य रूप से, फ्लक्स अवरोधक शब्द का उपयोग तब किया जाता है। जब अवरोधक प्रणाली फ्लक्स पर कार्य करता है, और प्रवणता अवरोधक का उपयोग तब किया जाता है जब अवरोधक प्रणाली स्तिथि (जैसे दबाव, वेग आदि) पर कार्य करता है।

कार्य विधि

इस प्रकार से फ्लक्स अवरोधक योजनाओं के निर्माण के पीछे मुख्य विचार स्थानिक व्युत्पन्नों को यथार्थवादी मान - वैज्ञानिक और इंजीनियरिंग समस्याओं के लिए ही सीमित करना है। इसका अर्थ सामान्यतः भौतिक रूप से प्राप्य और सार्थक मान हैं। और इनका उपयोग पीडीई द्वारा वर्णित समस्याओं को हल करने के लिए उच्च संकल्प योजनाओं में किया जाता है। और यह केवल तभी परिचालन में आते हैं। जब तीव्र तरंग लहर उपस्तिथ होते हैं। जिसे सुचारू रूप से परिवर्तित तरंगों के लिए, फ्लक्स अवरोधक संचालित नहीं होते हैं। और स्थानिक व्युत्पन्नों को स्पूरियस दोलनों को प्रस्तुत किए बिना उच्च क्रम सन्निकटन द्वारा दर्शाया जा सकता है। इस प्रकार से नीचे दी गई 1D अर्ध-असतत योजना पर विचार करें,

${\frac {du_{i}}{dt}}+{\frac {1}{\Delta x_{i}}}\left[F\left(u_{i+{1}/{2}}\right)-F\left(u_{i-{1}/{2}}\right)\right]=0,$

जहाँ, $F\left(u_{i+{1}/{2}}\right)$ और $F\left(u_{i-1/2}\right)$ i-th सेल के लिए एज फ्लक्स का प्रतिनिधित्व करते है। यदि इन किनारे के फ्लक्स को निम्न और उच्च संकल्प योजनाओं द्वारा दर्शाया जा सकता है, तब फ्लक्स अवरोधक विशेष सेल के समीप प्रवणता के आधार पर इन योजनाओं के मध्य स्विच कर सकता है, निम्नानुसार:

$$F\left(u_{i+1/2}\right)=f_{i+1/2}^{\text{low}}-\phi \left(r_{i}\right)\left(f_{i+1/2}^{\text{low}}-f_{i+1/2}^{\text{high}}\right),$$

$$F\left(u_{i-1/2}\right)=f_{i-1/2}^{\text{low}}-\phi \left(r_{i-1}\right)\left(f_{i-1/2}^{\text{low}}-f_{i-1/2}^{\text{high}}\right),$$

जहाँ

• $f^{\text{low}}$ निम्न विभेदन प्रवाह है,
• $f^{\text{high}}$ उच्च विभेदन प्रवाह है,
• $\phi \ (r)$ फ्लक्स अवरोधक फलन है, और
• $r$ समाधान मेश पर क्रमिक प्रवणता के अनुपात का प्रतिनिधित्व करता है, अर्थात,
  $$r_{i}={\frac {u_{i}-u_{i-1}}{u_{i+1}-u_{i}}}.$$

इस प्रकार से अवरोधक फलन शून्य से अधिक या उसके समान होने के लिए बाध्य है, अर्थात, $\phi \ (r)\geq 0$. इसलिए, जब अवरोधक शून्य (तीव्र प्रवणता, विपरीत प्रवणता या शून्य प्रवणता) के समान होता है, तो फ्लक्स को कम संकल्प योजना द्वारा दर्शाया जाता है। इसी प्रकार, जब अवरोधक 1 (सुचारू समाधान) के समान होता है, तो इसे उच्च संकल्प योजना द्वारा दर्शाया जाता है। और विभिन्न सीमाओं में भिन्न-भिन्न स्विचिंग विशेषताएँ होती हैं और उन्हें विशेष समस्या और समाधान योजना के अनुसार चुना जाता है। सभी समस्याओं के लिए उचित कार्य करने वाला कोई विशेष अवरोधक नहीं पाया गया है, और विशेष विकल्प सामान्यतः परीक्षण और त्रुटि के आधार पर बनाया जाता है।

अवरोधक कार्य

जहाँ $\phi (r)$,फ़्लक्स/प्रवणता अवरोधक फलन के सामान्य रूप निम्नलिखित हैं:

• आकर्षण [दूसरे क्रम का टीवीडी नहीं]
  $$\phi _{cm}(r)={\begin{cases}{\frac {r\left(3r+1\right)}{\left(r+1\right)^{2}}},&r>0,&\lim _{r\to \infty }\phi _{cm}(r)=3\\0\quad \quad \,,&r\leq 0\end{cases}}$$
• एचसीयूएस [दूसरा क्रम टीवीडी नहीं]
  $$\phi _{hc}(r)={\frac {1.5\left(r+\left|r\right|\right)}{\left(r+2\right)}};\quad \lim _{r\to \infty }\phi _{hc}(r)=3.$$
• त्वरित [दूसरे क्रम के टीवीडी नोट्स]
  $$\phi _{hq}(r)={\frac {2\left(r+\left|r\right|\right)}{\left(r+3\right)}};\quad \lim _{r\to \infty }\phi _{hq}(r)=4.$$
• कोरेन - पर्याप्त रूप से सुचारू डेटा के लिए तृतीय क्रम का स्पष्ट
  $$\phi _{kn}(r)=\max \left[0,\min \left(2r,\min \left({\dfrac {(1+2r)}{3}},2\right)\right)\right];\quad \lim _{r\to \infty }\phi _{kn}(r)=2.$$
• मिनमोड - सममित
  $$\phi _{mm}(r)=\max \left[0,\min \left(1,r\right)\right];\quad \lim _{r\to \infty }\phi _{mm}(r)=1.$$
• मोनोटोनाइज्ड सेंट्रल (एमसी) - सममित
  $$\phi _{mc}(r)=\max \left[0,\min \left(2r,0.5(1+r),2\right)\right];\quad \lim _{r\to \infty }\phi _{mc}(r)=2.$$
• ओशर
  $$\phi _{os}(r)=\max \left[0,\min \left(r,\beta \right)\right],\quad \left(1\leq \beta \leq 2\right);\quad \lim _{r\to \infty }\phi _{os}(r)=\beta .$$
• ओस्प्रे - सममित
  $$\phi _{op}(r)={\frac {1.5\left(r^{2}+r\right)}{\left(r^{2}+r+1\right)}};\quad \lim _{r\to \infty }\phi _{op}(r)=1.5\,.$$
• स्मार्ट [दूसरे क्रम का टीवीडी नहीं]
  $$\phi _{sm}(r)=\max \left[0,\min \left(2r,\left(0.25+0.75r\right),4\right)\right];\quad \lim _{r\to \infty }\phi _{sm}(r)=4.$$
• सुपरबी - सममित
  $$\phi _{sb}(r)=\max \left[0,\min \left(2r,1\right),\min \left(r,2\right)\right];\quad \lim _{r\to \infty }\phi _{sb}(r)=2.$$
• स्वेबी - सममित
  $$\phi _{sw}(r)=\max \left[0,\min \left(\beta r,1\right),\min \left(r,\beta \right)\right],\quad \left(1\leq \beta \leq 2\right);\quad \lim _{r\to \infty }\phi _{sw}(r)=\beta .$$
• यूएमआईएसटी - सममित
  $$\phi _{um}(r)=\max \left[0,\min \left(2r,\left(0.25+0.75r\right),\left(0.75+0.25r\right),2\right)\right];\quad \lim _{r\to \infty }\phi _{um}(r)=2.$$
• वैन अल्बाडा 1 - सममित
  $$\phi _{va1}(r)={\frac {r^{2}+r}{r^{2}+1}};\quad \lim _{r\to \infty }\phi _{va1}(r)=1.$$
• वैन अल्बाडा 2 - वैकल्पिक रूप [2रे क्रम का टीवीडी नहीं] उच्च स्थानिक क्रम योजनाओं पर उपयोग किया जाता है
  $$\phi _{va2}(r)={\frac {2r}{r^{2}+1}};\quad \lim _{r\to \infty }\phi _{va2}(r)=0.$$
• वैन लीयर - सममित
  $$\phi _{vl}(r)={\frac {r+\left|r\right|}{1+\left|r\right|}};\quad \lim _{r\to \infty }\phi _{vl}(r)=2.$$
• सममिति के रूप में संकेतित उपरोक्त सभी सीमाएं निम्नलिखित समरूपता गुण प्रदर्शित करती हैंː
  $${\frac {\phi \left(r\right)}{r}}=\phi \left({\frac {1}{r}}\right).$$

अतः यह वांछनीय गुण है क्योंकि यह सुनिश्चित करता है, कि आगे और पीछे के प्रवणता के लिए सीमित क्रियाएं समान विधि से संचालित होती हैं।

दूसरे क्रम की टीवीडी योजनाओं के लिए स्वीकार्य अवरोधक क्षेत्र।

जब तक इसके विपरीत संकेत न दिया जाए, उपरोक्त अवरोधक कार्य दूसरे क्रम के टीवीडी हैं। इसका अर्थ यह है कि उन्हें इस तरह से डिज़ाइन किया गया है कि वे योजना की स्थिरता की प्रत्याभूत के लिए समाधान के निश्चित क्षेत्र से निकलते हैं, जिसे टीवीडी क्षेत्र के रूप में जाना जाता है। इस प्रकार से दूसरे क्रम के, टीवीडी अवरोधक कम से कम निम्नलिखित मानदंडों को पूर्ण करते हैं:

• $r\leq \phi (r)\leq 2r,\left(0\leq r\leq 1\right)\$,
• $1\leq \phi (r)\leq r,\left(1\leq r\leq 2\right)\$,
• $1\leq \phi (r)\leq 2,\left(r>2\right)\$,
• $\phi (1)=1\$,

दूसरे क्रम की टीवीडी योजनाओं के लिए स्वीकार्य अवरोधक क्षेत्र स्वेबी आरेख में विपरीत दिखाया गया है, और टीवीडी क्षेत्र पर अवरोधक फ़ंक्शंस को दिखाने वाले प्लॉट नीचे दिखाए गए हैं। इस छवि में, ओशर और स्वेबी अवरोधक के लिए प्लॉट $\beta =1.5$ का उपयोग करके तैयार किए गए हैं .

अवरोधक फलन दूसरे क्रम के टीवीडी क्षेत्र पर मढ़ा हुआ है।

सामान्यीकृत मिनमॉड अवरोधक

एक अतिरिक्त अवरोधक जिसका रोचक रूप है, वैन-लीयर का मिनमॉड अवरोधक का एक-पैरामीटर वर्ग है। इसे इस प्रकार परिभाषित किया गया है।

$$\phi _{mg}(u,\theta )=\max \left(0,\min \left(\theta r,{\frac {1+r}{2}},\theta \right)\right),\quad \theta \in \left[1,2\right].$$

टिप्पणी: $\phi _{mg}$, $\theta =1,$ के लिए सर्वाधिक विघटनकारी है जब यह घटकर $\phi _{mm},$ हो जाता है और $\theta =2$ के लिए अधिक कम विघटनकारी है।

यह भी देखें

• गोडुनोव का प्रमेय
• उच्च संकल्प योजना
• एमयूएससीएल योजना
• सर्गेई के. गोडुनोव
• कुल भिन्नता कम होना।

संदर्भ

• Hirsch, C. (1990), Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows, Wiley, ISBN 978-0-471-92452-4
• Leonard, B.P.; Leschziner, M.A.; McGuirk, J. (1978), "The QUICK algorithm: a uniformly 3rd-order finite-difference method for highly convective flows", Proc. 1st Conf. on Numerical Methods in Laminar & Turbulent Flow, Swansea, p. 807

अग्रिम पठन

• Laney, Culbert B. (1998), Computational Gasdynamics, Cambridge University Press, ISBN 978-0-521-57069-5
• LeVeque, Randall (1990), Numerical Methods for Conservation Laws, ETH Lectures in Mathematics Series, Birkhauser-Verlag, ISBN 3-7643-2464-3
• LeVeque, Randall (2002), Finite Volume Methods for Hyperbolic Problems, Cambridge University Press, ISBN 0-521-00924-3
• Toro, E.F. (1999), Riemann Solvers and Numerical Methods for Fluid Dynamics (2nd ed.), Springer-Verlag, ISBN 3-540-65966-8
• Tannehill, John C.; Anderson, Dale Arden; Pletcher, Richard H. (1997), Computational Fluid Mechanics and Heat Transfer (2nd ed.), Taylor and Francis, ISBN 1-56032-046-X
• Wesseling, Pieter (2001), Principles of Computational Fluid Dynamics, Springer-Verlag, ISBN 3-540-67853-0