फ्लक्स सीमक

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

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

वे कैसे काम करते हैं

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

${\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)$ आई-वें सेल के लिए एज फ्लक्स का प्रतिनिधित्व करें। यदि इन किनारे के फ्लक्स को निम्न और उच्च रिज़ॉल्यूशन योजनाओं द्वारा दर्शाया जा सकता है, तो एक फ्लक्स लिमिटर विशेष सेल के करीब ग्रेडिएंट के आधार पर इन योजनाओं के बीच स्विच कर सकता है, निम्नानुसार:

$$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 .$$
• UMIST - सममित
  $$\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$.

Limiter functions overlaid onto second-order TVD region.

सामान्यीकृत मिनमॉड सीमक

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

$$\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$.

यह भी देखें

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

संदर्भ

• 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{{citation}}: CS1 maint: location missing publisher (link)

अग्रिम पठन

• 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