एचपी-एफईएम: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
 
(9 intermediate revisions by 3 users not shown)
Line 1: Line 1:
'''एचपी-एफईएम''' परिमित अवयव विधि (एफईएम) का सामान्य संस्करण होता है | जो खंड अनुसार-बहुपद सन्निकटन के आधार पर आंशिक अंतर समीकरणों का समाधान करने के लिए [[संख्यात्मक विश्लेषण]] विधि होती है | और जो वैरिएबल आकार ''(h)'' और [[बहुपद की डिग्री]] ''(p)'' के अवयवों को नियोजित करता है। एचपी-एफईएम की उत्पत्ति बार्ना A सज़ाबो और इवो बाबूस्का के अग्रणी कार्य से हुई है | <ref>[[Barna Szabó|B. A. Szabó]], A. K. Mehta: p-Convergent Finite Element Approximations in Fracture Mechanics, Int. J. Num. Meth. Engng, Volume 12, pp. 551-560, 1978.</ref> <ref>[[Ivo Babuška|I. Babuška]], [[Barna Szabó|B. A. Szabó]] and I. N. Katz: The p-Version of the Finite Element Method, SIAM J. Numer. Anl., Volume 18, pp. 515-544, 1981.</ref> <ref>[[Ivo Babuška|I. Babuška]], [[Barna Szabó|B. A. Szabó]], On the Rates of Convergence of the Finite Element Method, Int. J. Numer. Meth.Engng., Volume 18, pp. 323-341, 1982.</ref> <ref>[[Ivo Babuška|I. Babuška]]: The p- and hp-Versions of the Finite Element Method: the State of the Art, Finite Elements: Theory and Applications, edited by D. L. Dwoyer, M. Y. Hussaini and R. G. Voigt, New York, Springer-Verlag, 1988.</ref> <ref>[[Barna Szabó|B. A. Szabó]], [[Ivo Babuška|I. Babuška]]: Finite Element Analysis, John Wiley & Sons, {{ISBN|978-0-471-50273-9}}, 1991.</ref> <ref>[[Ivo Babuška|I. Babuška]], B.Q. Guo: The h, p and h-p version of the finite element method: basis theory and applications, Advances in Engineering Software, Volume 15, Issue 3-4, 1992.</ref> जिन्होंने इसमें पाया कि परिमित अवयव विधि शीघ्रता से परिवर्तित होती है। और जालक को h-शोधन (अवयवों को लघु भागों में विभाजित करना) होता हैं | इस प्रकार p-शोधन (उनकी बहुपद डिग्री को बढ़ाने) के उपयुक्त संयोजन का उपयोग करके परिष्कृत किया जाता है। यह घातीय अभिसरण अधिकांश अन्य परिमित अवयव विधियों की तुलना में विधि को बहुत आकर्षक बनाता है | जिससे यह सिर्फ बीजगणितीय दर के साथ अभिसरण करता है। एचपी-एफईएम के घातीय अभिसरण का पूर्वानुमान न सिर्फ सैद्धांतिक रूप से किया गया था, किंतु यह अनेक स्वतंत्र शोधकर्ताओं द्वारा भी देखा गया था। <ref>J.M. Melenk: hp-Finite Element Methods for Singular Perturbations, Springer, 2002</ref> <ref>C. Schwab: p- and hp- Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics, Oxford University Press, 1998</ref> <ref>P. Solin: Partial Differential Equations and the Finite Element Method, J. Wiley & Sons, 2005</ref>
==मानक एफईएम से अंतर==
एचपी-एफईएम अनेक पहलुओं में मानक (निम्नतम-क्रम) एफईएम से भिन्न होते है।<ref>P. Solin, K. Segeth, I. Dolezel: Higher-Order Finite Element Methods, Chapman & Hall/CRC Press, 2003</ref>
*उच्च-क्रम आकार कार्यों का चयन उदाहरण आवश्यक: अवयवों में उच्च-डिग्री बहुपद को आकार कार्यों के विभिन्न समुच्चयो का उपयोग करके उत्पन्न किया जा सकता है। ऐसे समुच्चय का चयन कठोरता आव्युह की कंडीशनिंग और उसी स्थान में संपूर्ण समाधान प्रक्रिया को नाटकीय रूप से प्रभावित कर सकता है। इस समस्या को सबसे पहले बाबुस्का एट अल द्वारा प्रलेखित किया गया था।<ref>I. Babuska, M. Griebel and J. Pitkaranta, The problem of selecting the shape functions for a p-type finite element, Internat. J. Numer. Methods Engrg. (1989), pp. 1891–1908</ref>
*स्वचालित एचपी-अनुकूलन: एचपी-एफईएम में, अवयव को अनेक भिन्न-भिन्न विधियों से एचपी-परिष्कृत किया जा सकता है, जैसे: इसे स्पेस में उप-विभाजित किए बिना इसकी बहुपद डिग्री बढ़ाना, या अवयव को ज्यामितीय रूप से उप-विभाजित करना होता हैं | जहां विभिन्न बहुपद डिग्री को उप-अवयवों पर प्रयुक्त किया जा सकता है। और अवयव शोधन प्रत्याशी की संख्या सरलता से दो आयामों में 100 और तीन आयामों में 1000 तक पहुंच जाती है। किसी अवयव में त्रुटि के आकार को सांकेतिक करने वाली संख्या स्वचालित hp-अनुकूलता को निर्देशित करने के लिए पर्याप्त नहीं होती है | यह (मानक एफईएम में अनुकूलता के विपरीत) होती हैं । इस प्रकार यह प्रत्येक अवयव में त्रुटि के आकार के बारे में अधिक सूचना प्राप्त करने के लिए संदर्भ समाधान या विश्लेषणात्मक विचार जैसी अन्य तकनीकों को नियोजित किया जाना चाहिए। <ref>L. Demkowicz, W. Rachowicz, and Ph. Devloo:  A Fully Automatic hp-Adaptivity, Journal of Scientific Computing, 17, Nos 1–3 (2002), 127–155</ref>
* संयोजन और समाधान सीपीयू समय का अनुपात: मानक एफईएम में, कठोरता आव्युह होता हैं | यह सामान्यतः शीघ्रता से एकत्रित किया जाता है किन्तु यह अधिक विस्तृत होता है। यह सामान्यतः, असतत समस्या के समाधान में कुल कंप्यूटिंग समय का सबसे विस्तृत भाग व्यय होता है। और इसके विपरीत, hp- में एफईएम में कठोरता आव्युह सामान्यतः बहुत लघु होते हैं | किन्तु यह (समान आव्युह आकार के लिए) होता हैं और उनमें एकत्रित मानक एफईएम की तुलना में अधिक समय लगता है। यह मुख्य रूप से संख्यात्मक चतुर्भुज की कम्प्यूटेशनल निवेश के कारण होता है | जिसमें शीघ्र अभिसरण दरों का लाभ उठाने के लिए मानक एफईएम की तुलना में उच्च परिशुद्धता होनी चाहिए | और इसलिए यह उच्च क्रम का होना चाहिए।
*विश्लेषणात्मक चुनौतियाँ: hp-एफईएम को सामान्यतः मानक एफईएम की तुलना में विश्लेषणात्मक दृष्टिकोण से समझना अधिक कठिन माना जाता है। जिसके अनुसार यह अनेक तकनीकों से संबंधित होता है, जैसे वृत्ताकार समस्याओं के लिए असतत अधिकतम सिद्धांत (डीएमपी) होता हैं। यह परिणाम बताते हैं कि, सामान्यतः जालक पर कुछ सीमित धारणाओं के साथ, खंड अनुसार-बहुपद एफईएम सन्निकटन अंतर्निहित वृत्ताकार पीडीई के समान अधिकतम सिद्धांत का पालन करता है। ऐसे परिणाम बहुत महत्वपूर्ण होते हैं क्योंकि वे गारंटी देते हैं कि सन्निकटन भौतिक रूप से स्वीकार्य रहता है | जिससे नकारात्मक घनत्व, नकारात्मक एकाग्रता, या नकारात्मक निरपेक्ष तापमान की गणना करने की कोई संभावना नहीं रहती है। डीएमपी निम्नतम-क्रम एफईएम के लिए अधिक अच्छी तरह से समझा जाता है किन्तु दो या दो से अधिक आयामों में hp-एफईएम के लिए पूर्ण तरह से अज्ञात होता है। इस प्रकार यह स्थानिक आयाम में प्रथम डीएमपी वर्तमान में तैयार किया गया था। <ref>P. Solin, T. Vejchodsky: A Weak Discrete Maximum Principle for hp-FEM, J. Comput. Appl. Math. 209 (2007) 54–65</ref> <ref>T. Vejchodsky, P. Solin: Discrete Maximum Principle for Higher-Order Finite Elements in 1D, Math. Comput. 76 (2007), 1833–1846</ref>
* प्रोग्रामिंग चुनौतियाँ: मानक एफईएम कोड की तुलना में hp-एफईएम सॉल्वर को प्रयुक्त करना बहुत कठिन होता है। जिनमें अनेक विवादों को दूर करने की आवश्यकता होती है | यह उनमें सम्मिलित होता हैं (किन्तु यह सिर्फ यहीं तक सीमित नहीं होता हैं) | उच्च-क्रम चतुर्भुज सूत्र, उच्च-क्रम आकार फ़ंक्शन, भौतिक डोमेन में आधार कार्यों के साथ संदर्भ डोमेन पर आकार कार्यों से संबंधित कनेक्टिविटी और अभिविन्यास सूचना आदि होते हैं। <ref>L. Demkowicz, J. Kurtz, D. Pardo, W. Rachowicz, M. Paszynski, A. Zdunek: Computing with hp-Adaptive Finite Elements, Chapman & Hall/CRC Press, 2007</ref>




==फ़िचेरा समस्या                                                                                                                                ==


 
फिचेरा समस्या (जिसे फिचेरा कॉर्नर समस्या भी कहा जाता है) | यह अनुकूल एफईएम कोड के लिए मानक बेंचमार्क समस्या होती है। कोई इसका उपयोग मानक एफईएम और hp-एफईएम के प्रदर्शन में नाटकीय अंतर दिखाने के लिए कर सकता है। यह समस्या ज्यामिति घन होता है जिसका कॉर्नर लुप्त होता है।इसमें स्पष्ट समाधान के केंद्र में विलक्षण स्लोप (अनंत तनाव का सादृश्य) होता है। इसमें स्पष्ट समाधान का ज्ञान सन्निकटन त्रुटि की स्पष्ट गणना करना हैं और इस प्रकार यह विभिन्न संख्यात्मक विधियों की तुलना को संभव बनाता है। उदाहरण के लिए, समस्या को अनुकूली एफईएम के तीन भिन्न-भिन्न संस्करणों का उपयोग करके समाधान किया गया था | जिसमे यह रैखिक अवयवों, द्विघात अवयवों और hp-एफईएम के साथ होता हैं।                                                           
एचपी-एफईएम परिमित तत्व विधि (एफईएम) का सामान्य संस्करण है, जो टुकड़े-टुकड़े-बहुपद सन्निकटन के आधार पर आंशिक अंतर समीकरणों को हल करने के लिए [[संख्यात्मक विश्लेषण]] विधि है जो चर आकार के तत्वों को नियोजित करता है। ''(एच)'' और [[बहुपद की डिग्री]] ''(पी)''। एचपी-एफईएम की उत्पत्ति बार्ना ए. सज़ाबो और इवो बाबुस्का के अग्रणी कार्य से हुई है।<ref>[[Barna Szabó|B. A. Szabó]], A. K. Mehta: p-Convergent Finite Element Approximations in Fracture Mechanics, Int. J. Num. Meth. Engng, Volume 12, pp. 551-560, 1978.</ref><ref>[[Ivo Babuška|I. Babuška]], [[Barna Szabó|B. A. Szabó]] and I. N. Katz: The p-Version of the Finite Element Method, SIAM J. Numer. Anl., Volume 18, pp. 515-544, 1981.</ref><ref>[[Ivo Babuška|I. Babuška]], [[Barna Szabó|B. A. Szabó]], On the Rates of Convergence of the Finite Element Method, Int. J. Numer. Meth.Engng., Volume 18, pp. 323-341, 1982.</ref><ref>[[Ivo Babuška|I. Babuška]]: The p- and hp-Versions of the Finite Element Method: the State of the Art, Finite Elements: Theory and Applications, edited by D. L. Dwoyer, M. Y. Hussaini and R. G. Voigt, New York, Springer-Verlag, 1988.</ref><ref>[[Barna Szabó|B. A. Szabó]], [[Ivo Babuška|I. Babuška]]: Finite Element Analysis, John Wiley & Sons, {{ISBN|978-0-471-50273-9}}, 1991.</ref><ref>[[Ivo Babuška|I. Babuška]], B.Q. Guo: The h, p and h-p version of the finite element method: basis theory and applications, Advances in Engineering Software, Volume 15, Issue 3-4, 1992.</ref> किसने पता लगाया कि परिमित तत्व विधि तेजी से तेजी से परिवर्तित होती है
जाल को एच-शोधन के उपयुक्त संयोजन का उपयोग करके परिष्कृत किया जाता है
(तत्वों को छोटे-छोटे भागों में विभाजित करना) और पी-एफईएम|पी-शोधन (उन्हें बढ़ाना)।
बहुपद डिग्री)। घातीय अभिसरण अधिकांश अन्य परिमित तत्व विधियों की तुलना में विधि को बहुत आकर्षक बनाता है, जो केवल बीजगणितीय दर के साथ अभिसरण करता है। घातीय अभिसरण
एचपी-एफईएम की न केवल सैद्धांतिक रूप से भविष्यवाणी की गई, बल्कि इसका अवलोकन भी किया गया कई स्वतंत्र शोधकर्ताओं द्वारा।<ref>J.M. Melenk: hp-Finite Element Methods for Singular Perturbations, Springer, 2002</ref><ref>C. Schwab: p- and hp- Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics, Oxford University Press, 1998</ref><ref>P. Solin: Partial Differential Equations and the Finite Element Method, J. Wiley & Sons, 2005</ref>
 
 
==मानक FEM से अंतर==
एचपी-एफईएम कई पहलुओं में मानक (निम्नतम-क्रम) एफईएम से भिन्न है।<ref>P. Solin, K. Segeth, I. Dolezel: Higher-Order Finite Element Methods, Chapman & Hall/CRC Press, 2003</ref>
* उच्च-क्रम आकार के कार्यों का विकल्: तत्वों में उच्च-डिग्री बहुपद को आकार कार्यों के विभिन्न सेटों का उपयोग करके उत्पन्न किया जा सकता है। ऐसे सेट का चुनाव कठोरता मैट्रिक्स की कंडीशनिंग और बदले में संपूर्ण समाधान प्रक्रिया को नाटकीय रूप से प्रभावित कर सकता है। इस समस्या को सबसे पहले बाबुस्का एट अल द्वारा प्रलेखित किया गया था।<ref>I. Babuska, M. Griebel and J. Pitkaranta, The problem of selecting the shape functions for a p-type finite element, Internat. J. Numer. Methods Engrg. (1989), pp. 1891–1908</ref>
* स्वचालित एचपी-अनुकूलन: एचपी-एफईएम में, तत्व को कई अलग-अलग तरीकों से एचपी-परिष्कृत किया जा सकता है, जैसे: इसे अंतरिक्ष में उप-विभाजित किए बिना इसकी बहुपद डिग्री बढ़ाना, या तत्व को ज्यामितीय रूप से उप-विभाजित करना, जहां विभिन्न बहुपद डिग्री को उप-तत्वों पर लागू किया जा सकता है। तत्व शोधन उम्मीदवारों की संख्या आसानी से दो आयामों में 100 और तीन आयामों में 1000 तक पहुंच जाती है। किसी तत्व में त्रुटि के आकार को इंगित करने वाली संख्या स्वचालित एचपी-अनुकूलता को निर्देशित करने के लिए पर्याप्त नहीं है (मानक एफईएम में अनुकूलता के विपरीत)। प्रत्येक तत्व में ''त्रुटि के आकार'' के बारे में अधिक जानकारी प्राप्त करने के लिए ''संदर्भ समाधान'' या ''विश्लेषणात्मकता विचार'' जैसी अन्य तकनीकों को नियोजित किया जाना चाहिए।<ref>L. Demkowicz, W. Rachowicz, and Ph. Devloo:  A Fully Automatic hp-Adaptivity, Journal of Scientific Computing, 17, Nos 1–3 (2002), 127–155</ref>
* असेंबलिंग और समाधान सीपीयू समय का अनुपात: मानक एफईएम में, कठोरता मैट्रिक्स आमतौर पर जल्दी से असेंबल किया जाता है लेकिन यह काफी बड़ा होता है। आमतौर पर, असतत समस्या के समाधान में कुल कंप्यूटिंग समय का सबसे बड़ा हिस्सा खर्च होता है। इसके विपरीत, एचपी-एफईएम में कठोरता मैट्रिक्स आमतौर पर बहुत छोटे होते हैं, लेकिन (समान मैट्रिक्स आकार के लिए) उनकी असेंबली में मानक एफईएम की तुलना में अधिक समय लगता है। यह मुख्य रूप से संख्यात्मक चतुर्भुज की कम्प्यूटेशनल लागत के कारण है, जिसमें तेज अभिसरण दरों का लाभ उठाने के लिए मानक एफईएम की तुलना में उच्च परिशुद्धता होनी चाहिए, और इसलिए उच्च क्रम का होना चाहिए।
* विश्लेषणात्मक चुनौतियाँ: एचपी-एफईएम को आम तौर पर मानक एफईएम की तुलना में विश्लेषणात्मक दृष्टिकोण से समझना अधिक कठिन माना जाता है। यह कई तकनीकों से संबंधित है, जैसे अण्डाकार समस्याओं के लिए असतत अधिकतम सिद्धांत (डीएमपी)। ये परिणाम बताते हैं कि, आमतौर पर जाल पर कुछ सीमित धारणाओं के साथ, टुकड़े-टुकड़े-बहुपद एफईएम सन्निकटन अंतर्निहित अण्डाकार पीडीई के समान अधिकतम सिद्धांतों का पालन करता है। ऐसे परिणाम बहुत महत्वपूर्ण हैं क्योंकि वे गारंटी देते हैं कि सन्निकटन भौतिक रूप से स्वीकार्य रहता है, जिससे नकारात्मक घनत्व, नकारात्मक एकाग्रता, या नकारात्मक निरपेक्ष तापमान की गणना करने की कोई संभावना नहीं बचती है। डीएमपी निम्नतम-क्रम एफईएम के लिए काफी अच्छी तरह से समझा जाता है लेकिन दो या दो से अधिक आयामों में एचपी-एफईएम के लिए पूरी तरह से अज्ञात है। स्थानिक आयाम में पहला डीएमपी हाल ही में तैयार किया गया था।<ref>P. Solin, T. Vejchodsky: A Weak Discrete Maximum Principle for hp-FEM, J. Comput. Appl. Math. 209 (2007) 54–65</ref><ref>T. Vejchodsky, P. Solin: Discrete Maximum Principle for Higher-Order Finite Elements in 1D, Math. Comput. 76 (2007), 1833–1846</ref>
* प्रोग्रामिंग चुनौतियाँ: मानक FEM कोड की तुलना में hp-FEM सॉल्वर को लागू करना बहुत कठिन है। जिन कई मुद्दों को दूर करने की आवश्यकता है उनमें शामिल हैं (लेकिन केवल यहीं तक सीमित नहीं हैं): उच्च-क्रम चतुर्भुज सूत्र, उच्च-क्रम आकार फ़ंक्शन, भौतिक डोमेन में आधार कार्यों के साथ संदर्भ डोमेन पर आकार कार्यों से संबंधित कनेक्टिविटी और अभिविन्यास जानकारी, आदि।<ref>L. Demkowicz, J. Kurtz, D. Pardo, W. Rachowicz, M. Paszynski, A. Zdunek: Computing with hp-Adaptive Finite Elements, Chapman & Hall/CRC Press, 2007</ref>
 
 
==फ़िचेरा समस्या==
 
फिचेरा समस्या (जिसे फिचेरा कॉर्नर समस्या भी कहा जाता है) अनुकूली FEM कोड के लिए मानक बेंचमार्क समस्या है। कोई इसका उपयोग मानक FEM और hp-FEM के प्रदर्शन में नाटकीय अंतर दिखाने के लिए कर सकता है। समस्या ज्यामिति घन है जिसका कोना लुप्त है। सटीक समाधान के केंद्र में विलक्षण ढाल (अनंत तनाव का सादृश्य) है। सटीक समाधान का ज्ञान सन्निकटन त्रुटि की सटीक गणना करना और इस प्रकार विभिन्न संख्यात्मक तरीकों की तुलना करना संभव बनाता है। उदाहरण के लिए, समस्या को अनुकूली FEM के तीन अलग-अलग संस्करणों का उपयोग करके हल किया गया था: रैखिक तत्वों, द्विघात तत्वों और hp-FEM के साथ।


<gallery>
<gallery>
Line 27: Line 17:
Image:conv fichera.png|समस्या चार्ट: अभिसरण तुलना।
Image:conv fichera.png|समस्या चार्ट: अभिसरण तुलना।
</gallery>
</gallery>
अभिसरण ग्राफ स्वतंत्रता की डिग्री (डीओएफ) की संख्या के फ़ंक्शन के रूप में सन्निकटन त्रुटि दिखाते हैं। डीओएफ अज्ञात मापदंडों को संदर्भित करता है जो सन्निकटन को परिभाषित करने के लिए आवश्यक हैं, और डीओएफ की संख्या कठोरता मैट्रिक्स के आकार के बराबर है। पाठक ग्राफ़ में देख सकते हैं कि एचपी-एफईएम का अभिसरण अन्य दोनों तरीकों के अभिसरण की तुलना में बहुत तेज़ है। प्रदर्शन अंतर इतना बड़ा है कि रैखिक एफईएम बिल्कुल भी (उचित समय में) अभिसरण नहीं कर सकता है और द्विघात एफईएम को उस सटीकता तक पहुंचने के लिए सैकड़ों हजारों या शायद लाखों डीओएफ की आवश्यकता होगी जो एचपी-एफईएम ने लगभग 17,000 डीओएफ के साथ प्राप्त की थी। स्वतंत्रता की अपेक्षाकृत कुछ डिग्री का उपयोग करके बहुत सटीक परिणाम प्राप्त करना एचपी-एफईएम की मुख्य ताकत है।
अभिसरण ग्राफ स्वतंत्रता की डिग्री (डीओएफ) की संख्या के फ़ंक्शन के रूप में सन्निकटन त्रुटि दिखाते हैं। और डीओएफ अज्ञात मापदंडों को संदर्भित करता है | जो सन्निकटन को परिभाषित करने के लिए आवश्यक होते हैं | और इसमें डीओएफ की संख्या कठोरता आव्युह के आकार के सामान्य होता है। इसको रीडर ग्राफ़ में देख सकते हैं कि hp-एफईएम का अभिसरण अन्य दोनों विधियों के अभिसरण की तुलना में बहुत शीघ्र होता है। इसमें प्रदर्शन अंतर इतना विस्तृत होता है कि रैखिक एफईएम पूर्णतया सभी (उचित समय में) अभिसरण नहीं कर सकते है और द्विघात एफईएम को उस स्पष्टता तक पहुंचने के लिए सैकड़ों हजारों या संभवतः लाखों डीओएफ की आवश्यकता होती हैं जो एचपी-एफईएम ने लगभग 17,000 डीओएफ के साथ प्राप्त की थी। यह स्वतंत्रता की अपेक्षाकृत कुछ डिग्री का उपयोग करके बहुत स्पष्ट परिणाम प्राप्त करना एचपी-एफईएम की मुख्य शक्ति होती है।
 
==एचपी-एफईएम की दक्षता                                                                                                          ==
लघु -रेखीय अवयवों की तुलना में विस्तृत उच्च-क्रम वाले अवयवों का उपयोग करके सुचारू कार्यों का अधिक कुशलता से अनुमान लगाया जा सकता है। इसे नीचे दिए गए चित्र में दर्शाया गया है | जहां दो भिन्न-भिन्न जालकों पर शून्य डिरिचलेट सीमा स्थितियों के साथ आयामी पॉइसन समीकरण समाधान किया गया है। यह स्पष्ट समाधान साइन फ़ंक्शन होता है।


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


* बाएँ: दो रैखिक तत्वों से युक्त जाल।
[[Image:grad sin h.png|टुकड़े-टुकड़े-रैखिक सन्निकटन।]][[Image:conv sin hp.png|द्विघात सन्निकटन.]]
* दाएँ: द्विघात तत्व से युक्त जाल।


[[Image:grad sin h.png|टुकड़े-टुकड़े-रैखिक सन्निकटन।]]
जबकि दोनों स्तिथियों में अज्ञात की संख्या समान है यह (1 डीओएफ), संबंधित मानदंड में त्रुटियां क्रमशः 0.68 और 0.20 हैं। इसका कारण यह है कि द्विघात सन्निकटन खंड-रेखीय सन्निकटन की तुलना में लगभग 3.5 गुना अधिक कुशल था। जब हम कदम आगे बढ़ाते हैं और (a) चार रैखिक अवयवों की तुलना (b) चतुर्थक अवयव (p=4) से करते हैं, तब दोनों भिन्न-भिन्न समस्याओं में तीन डीओएफ होंते हैं | किन्तु चतुर्थक सन्निकटन लगभग 40 गुना अधिक कुशल होता हैं।
[[Image:conv sin hp.png|द्विघात सन्निकटन.]]जबकि दोनों मामलों में अज्ञात की संख्या समान है (1 डीओएफ), संबंधित मानदंड में त्रुटियां क्रमशः 0.68 और 0.20 हैं। इसका मतलब यह है कि द्विघात सन्निकटन टुकड़ा-रेखीय सन्निकटन की तुलना में लगभग 3.5 गुना अधिक कुशल था। जब हम कदम आगे बढ़ते हैं और (ए) चार रैखिक तत्वों की तुलना (बी) चतुर्थक तत्व (पी=4) से करते हैं, तो दोनों अलग-अलग समस्याओं में तीन डीओएफ होंगे, लेकिन चतुर्थक सन्निकटन लगभग 40 गुना अधिक कुशल होगा।


इसके विपरीत, छोटे निम्न-क्रम वाले तत्व बड़े उच्च-क्रम वाले तत्वों की तुलना में छोटे पैमाने की विशेषताओं जैसे विलक्षणताओं को बेहतर तरीके से पकड़ सकते हैं। एचपी-एफईएम इन दो दृष्टिकोणों के इष्टतम संयोजन पर आधारित है जो घातांकीय अभिसरण की ओर ले जाता है। ध्यान दें कि यह घातीय अभिसरण त्रुटि की धुरी बनाम स्वतंत्रता की डिग्री में व्यक्त किया गया है। वास्तविक जीवन के अनुप्रयोगों के लिए, हम आम तौर पर सटीकता के समान स्तर तक पहुंचने के लिए आवश्यक कम्प्यूटेशनल समय पर विचार करते हैं। इस प्रदर्शन संकेतक के लिए एच- और एचपी-शोधन समान परिणाम प्रदान कर सकते हैं, उदाहरण के लिए। अंतिम आंकड़ा यहां देखें <ref>{{Cite web|url=http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html|title = Microwave Oven — Hermes Examples Guide}}</ref> (वेब आर्काइव लिंक <ref>{{cite web |url=http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html |title=Microwave Oven — Hermes Examples Guide |website=hpfem.org |access-date=12 January 2022 |archive-url=https://web.archive.org/web/20180807173436/http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html |archive-date=7 August 2018 |url-status=dead}}</ref>). जैसे ही एच-एफईएम की तुलना में एचपी-एफईएम को प्रोग्राम करना और [[समानांतर कंप्यूटिंग]] करना कठिन हो जाता है, एचपी-शोधन की अभिसरण उत्कृष्टता अव्यावहारिक हो सकती है।
इसके विपरीत, लघु निम्न-क्रम वाले अवयव विस्तृत उच्च-क्रम वाले अवयवों की तुलना में लघु मापदंड की विशेषताओं जैसे विलक्षणताओं को उत्तम विधियों से पकड़ सकते हैं। hp-एफईएम इन दो दृष्टिकोणों के इष्टतम संयोजन पर आधारित होता है जो घातांकीय अभिसरण की ओर ले जाता है। ध्यान दें कि यह घातीय अभिसरण त्रुटि की धुरी और स्वतंत्रता की डिग्री में व्यक्त किया गया है। वास्तविक जीवन के अनुप्रयोगों के लिए, हम सामान्यतः स्पष्टता के समान स्तर तक पहुंचने के लिए आवश्यक कम्प्यूटेशनल समय पर विचार करते हैं। इस प्रदर्शन संकेतक के लिए h- और hp-शोधन समान परिणाम प्रदान कर सकते हैं | उदाहरण के लिए <ref>{{Cite web|url=http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html|title = Microwave Oven — Hermes Examples Guide}}</ref> (वेब आर्काइव लिंक <ref>{{cite web |url=http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html |title=Microwave Oven — Hermes Examples Guide |website=hpfem.org |access-date=12 January 2022 |archive-url=https://web.archive.org/web/20180807173436/http://hpfem.org/wp-content/uploads/doc-web/doc-examples/src/hermes2d/examples/maxwell/microwave-oven.html |archive-date=7 August 2018 |url-status=dead}}</ref>) पर अंतिम आंकड़ा देखते हैं | जैसे ही h-एफईएम की तुलना में hp-एफईएम को प्रोग्राम करना और [[समानांतर कंप्यूटिंग]] करना कठिन हो जाता है | और hp-शोधन की अभिसरण उत्कृष्टता अव्यावहारिक हो सकती है।


==एचपी-अनुकूलन==
==एचपी-अनुकूलन                                                                                                               ==


कुछ एफईएम साइटें एचपी-अनुकूलता को एच-अनुकूलता (उनकी बहुपद डिग्री को स्थिर रखते हुए अंतरिक्ष में तत्वों को विभाजित करना) और पी-अनुकूलता (केवल उनकी बहुपद डिग्री को बढ़ाना) के संयोजन के रूप में वर्णित करती हैं।. यह पूरी तरह से सटीक नहीं है, क्योंकि एचपी-अनुकूलता एच- और पी-अनुकूलता दोनों से काफी भिन्न है क्योंकि किसी तत्व का एचपी-शोधन कई अलग-अलग तरीकों से किया जा सकता है। पी-शोधन के अलावा, तत्व को अंतरिक्ष में उप-विभाजित किया जा सकता है (जैसा कि एच-अनुकूलता में), लेकिन उप-तत्वों पर बहुपद डिग्री के लिए कई संयोजन हैं। यह दाहिनी ओर के चित्र में दर्शाया गया है। उदाहरण के लिए, यदि त्रिकोणीय या चतुर्भुज तत्व को चार उप-तत्वों में विभाजित किया जाता है, जहां बहुपद डिग्री को अधिकतम दो तक भिन्न होने की अनुमति होती है, तो इससे 3^4 = 81 शोधन उम्मीदवार मिलते हैं (बहुपद अनिसोट्रोपिक उम्मीदवारों पर विचार नहीं किया जाता है)। अनुरूप रूप से, हेक्साहेड्रोन को आठ उप-तत्वों में विभाजित करना और अधिकतम दो द्वारा उनकी बहुपद डिग्री को बदलना 3^8 = 6,561 शोधन उम्मीदवार प्राप्त करता है। प्रति तत्व स्थिर संख्या प्रदान करने वाला मानक FEM त्रुटि अनुमान स्वचालित एचपी-अनुकूलन का मार्गदर्शन करने के लिए पर्याप्त नहीं है।
कुछ एफईएम साइटें एचपी-अनुकूलता को h-अनुकूलता (उनकी बहुपद डिग्री को स्थिर रखते हुए स्पेस में अवयवों को विभाजित करता) हैं | और p-अनुकूलता (सिर्फ उनकी बहुपद डिग्री को बढ़ाना) के संयोजन के रूप में वर्णित करती हैं। यह पूर्णता से स्पष्ट नहीं होती है | क्योंकि hp-अनुकूलता h- और p-अनुकूलता दोनों से अधिक भिन्न होती है क्योंकि किसी अवयव का hp-शोधन अनेक भिन्न-भिन्न विधियों से किया जा सकता है। p-शोधन के अतिरिक्त, अवयव को स्पेस में उप-विभाजित किया जा सकता है (जैसा कि h-अनुकूलता में) हैं, किन्तु उप-अवयवों पर बहुपद डिग्री के लिए अनेक संयोजन होते हैं। यह दाहिनी ओर के चित्र में दर्शाया गया है। उदाहरण के लिए, यदि त्रिकोणीय या चतुर्भुज अवयव को चार उप-अवयवों में विभाजित किया जाता है | जहां बहुपद डिग्री को अधिकतम दो तक भिन्न होने की अनुमति होती है | तब इससे 3^4 = 81 शोधन प्रत्याशी मिलते हैं | इसमें (बहुपद अनिसोट्रोपिक प्रत्याशी पर विचार नहीं किया जाता है)। इस रूप से, हेक्साहेड्रोन को आठ उप-अवयवों में विभाजित करना होता हैं और अधिकतम दो द्वारा उनकी बहुपद डिग्री को परिवर्तित 3^8 = 6,561 शोधन प्रत्याशी प्राप्त करता है। यह प्रति अवयव स्थिर संख्या प्रदान करने वाला मानक एफईएम त्रुटि अनुमान स्वचालित hp-अनुकूलन का मार्गदर्शन करने के लिए पर्याप्त नहीं होता है।


==उच्च-क्रम आकार के कार्य==
==उच्च-क्रम आकार के कार्य                                                                                       ==
मानक FEM में केवल ग्रिड से जुड़े आकार कार्यों के साथ काम करता है
मानक एफईएम में सिर्फ ग्रिड शीर्षों (तथाकथित शीर्ष कार्यों) से जुड़े आकार कार्यों के साथ कार्य करता है। इसके विपरीत, hp-एफईएम का उपयोग करते समय, व्यक्ति एज के कार्यों (अवयव एज से जुड़े), फेस के कार्यों (अवयव फेसों के अनुरूप - केवल 3 डी), और बबल कार्यों (उच्च-क्रम बहुपद जो अवयव सीमाओं पर लुप्त हो जाते हैं) यह इसका भी ध्यान रखता है। इस प्रकार निम्नलिखित छवियां इन कार्यों को दिखाती हैं | और यह (एकल अवयव तक सीमित) होते हैं |
शीर्ष (तथाकथित शीर्ष फलन)। इसके विपरीत, एचपी-एफईएम का उपयोग करते समय, व्यक्ति एज फ़ंक्शंस (संबंधित) पर भी ध्यान देता है
तत्व किनारों), चेहरे के कार्य (तत्व चेहरों के अनुरूप - केवल 3डी),
और बबल फ़ंक्शंस (उच्च-क्रम बहुपद जो गायब हो जाते हैं
तत्व सीमाएँ)। निम्नलिखित छवियां इन कार्यों को दिखाती हैं
(एकल तत्व तक सीमित):


<gallery>
<gallery>
Line 58: Line 44:
Image:bubble new.jpg|बुलबुला समारोह.
Image:bubble new.jpg|बुलबुला समारोह.
</gallery>
</gallery>
ध्यान दें: ये सभी फ़ंक्शन संपूर्ण तत्व इंटीरियर में परिभाषित हैं।
ध्यान दें: यह सभी फ़ंक्शन संपूर्ण अवयव इंटीरियर में परिभाषित होते हैं।


==ओपन सोर्स एचपी-एफईएम कोड==
==ओपन सोर्स एचपी-एफईएम कोड                                         ==
* डील.II: डील.II परिमित तत्व विधि का उपयोग करके आंशिक अंतर समीकरणों को हल करने के लिए निःशुल्क, ओपन-सोर्स लाइब्रेरी है।
* डील.II: डील.II परिमित अवयव विधि का उपयोग करके आंशिक अंतर समीकरणों का समाधान करने के लिए निःशुल्क, ओपन-सोर्स लाइब्रेरी है।
* [http://www.concepts.math.ethz.ch/ अवधारणाएं]: अण्डाकार समीकरणों के लिए C/C++ hp-FEM/DGFEM/BEM लाइब्रेरी SAM, ETH ज्यूरिख (स्विट्जरलैंड) और बर्लिन के तकनीकी विश्वविद्यालय (जर्मनी) में K. श्मिट के समूह में विकसित की गई।
*[http://www.concepts.math.ethz.ch/ अवधारणाएं]: एसएएम, ईटीएच ज्यूरिख (स्विट्जरलैंड) और टीयू बर्लिन (जर्मनी) में के. श्मिट के समूह में वृत्ताकार समीकरणों के लिए C/C++ hp-एफईएम/डीजीएफईएम/बीईएम लाइब्रेरी विकसित की गई हैं।
* 2dhp90, 3dhp90: अण्डाकार समस्याओं और मैक्सवेल के समीकरणों के लिए फोरट्रान कोड, ICES, UT ऑस्टिन में एल. डेमकोविज़ द्वारा विकसित।
*2dhp90, 3dhp90: वृत्ताकार समस्याओं और मैक्सवेल के समीकरणों के लिए फोरट्रान कोड आईसीईएस, यूटी ऑस्टिन में एल. डेमकोविज़ द्वारा विकसित होता हैं।
* पीएचएएमएल: समानांतर पदानुक्रमित अनुकूली बहु-स्तरीय परियोजना। अनुकूली जाल शोधन और मल्टी-ग्रिड समाधान तकनीकों का उपयोग करके वितरित मेमोरी समानांतर कंप्यूटर और मल्टी-कोर कंप्यूटर पर 2 डी अण्डाकार आंशिक अंतर समीकरणों के संख्यात्मक समाधान के लिए, संयुक्त राज्य अमेरिका के राष्ट्रीय मानक और प्रौद्योगिकी संस्थान में परिमित तत्व सॉफ्टवेयर विकसित किया गया।
* पीएचएएमएल: समानांतर पदानुक्रमित अनुकूल बहु-स्तरीय परियोजना हैं। जिसमे अनुकूल जालक शोधन और मल्टी-ग्रिड समाधान तकनीकों का उपयोग करके वितरित मेमोरी समानांतर कंप्यूटर और मल्टी-कोर कंप्यूटर पर 2 डी वृत्ताकार आंशिक अंतर समीकरणों के संख्यात्मक समाधान के लिए,राष्ट्रीय मानक और प्रौद्योगिकी संस्थान, संयुक्त राज्य अमेरिका में परिमित तत्व सॉफ्टवेयर विकसित किया गया है।
* [[हर्मीस परियोजना]]: पीडीई और मल्टीफिजिक्स पीडीई सिस्टम की विशाल विविधता के लिए अंतरिक्ष और अंतरिक्ष-समय अनुकूली एचपी-एफईएम सॉल्वरों के तेजी से प्रोटोटाइप के लिए सी/सी++/पायथन लाइब्रेरी, नेवादा विश्वविद्यालय, रेनो (यूएसए), थर्मो-मैकेनिक्स संस्थान, प्राग (चेक गणराज्य) और पिल्सेन (चेक गणराज्य) में वेस्ट बोहेमिया विश्वविद्यालय में एचपी-एफईएम समूह द्वारा विकसित - शीर्ष पर निर्मित [[एग्रोस2डी]] इंजीनियरिंग सॉफ्टवेयर के साथ हर्मीस पुस्तकालय का.
* [[हर्मीस परियोजना]]: पीडीई और मल्टीफिजिक्स पीडीई सिस्टम की विशाल विविधता के लिए स्पेस और स्पेस-समय अनुकूल hp-एफईएम सॉल्वरों के शीघ्रता से प्रोटोटाइप के लिए C/C++/पायथन लाइब्रेरी, नेवादा विश्वविद्यालय, रेनो (यूएसए), थर्मो-मैकेनिक्स संस्थान, प्राग (चेक गणराज्य) और पिल्सेन (चेक गणराज्य) में वेस्ट बोहेमिया विश्वविद्यालय में hp-एफईएम समूह द्वारा विकसित - [[एग्रोस2डी]] इंजीनियरिंग सॉफ्टवेयर के शीर्ष पर निर्मित हर्मीस पुस्तकालय हैं |
* [http://lsec.cc.ac.cn/phg/index_en.htm PHG]: PHG समानांतर अनुकूली परिमित तत्व प्रोग्राम विकसित करने के लिए टूलबॉक्स है। यह h-, p- और hp-fem के लिए उपयुक्त है। पीएचजी वर्तमान में वैज्ञानिक और इंजीनियरिंग कंप्यूटिंग की राज्य प्रमुख प्रयोगशाला, कम्प्यूटेशनल गणित संस्थान और चीनी विज्ञान अकादमी (एलएसईसी, सीएएस, चीन) के वैज्ञानिक/इंजीनियरिंग कंप्यूटिंग संस्थान में सक्रिय विकास के अधीन है। पीएचजी अनुरूप टेट्राहेड्रल जाल से संबंधित है और संदेश भेजने के लिए अनुकूली स्थानीय जाल शोधन और एमपीआई के लिए द्विभाजन का उपयोग करता है। पीएचजी में ऑब्जेक्ट-ओरिएंटेड डिज़ाइन है जो समानांतर विवरण छुपाता है और अमूर्त तरीके से मेष और परिमित तत्व कार्यों पर सामान्य संचालन प्रदान करता है, जिससे उपयोगकर्ताओं को अपने संख्यात्मक एल्गोरिदम पर ध्यान केंद्रित करने की अनुमति मिलती है।
* [http://lsec.cc.ac.cn/phg/index_en.htm पीएचजी]: पीएचजी समानांतर अनुकूल परिमित अवयव प्रोग्राम विकसित करने के लिए टूलबॉक्स होता है। यह h-, p- और hp-फेम के लिए उपयुक्त होता है। और पीएचजी वर्तमान में वैज्ञानिक और इंजीनियरिंग कंप्यूटिंग की स्थान प्रमुख प्रयोगशाला हैं | कम्प्यूटेशनल गणित संस्थान और चीनी विज्ञान अकादमी (एलएसईसी, सीएएस, चीन) के वैज्ञानिक/इंजीनियरिंग कंप्यूटिंग संस्थान में सक्रिय विकास पर निर्भर होती है। पीएचजी अनुरूप टेट्राहेड्रल जालक से संबंधित होता है | और संदेश भेजने के लिए अनुकूल स्थानीय जालक शोधन और एमपीआई के लिए द्विभाजन का उपयोग करता है। पीएचजी में ऑब्जेक्ट-ओरिएंटेड डिज़ाइन होता है जो समानांतर विवरण गुप्त रखता है | और अमूर्त विधियों से मेष और परिमित अवयव कार्यों पर सामान्य संचालन प्रदान करता है, जिससे उपयोगकर्ताओं को अपने संख्यात्मक एल्गोरिदम पर ध्यान केंद्रित करने की अनुमति मिलती है।
* [http://mofem.eng.gla.ac.uk MoFEM] परिमित तत्व विश्लेषण कोड है जो बहु-भौतिकी समस्याओं के समाधान के लिए मनमाने ढंग से अनुमान के स्तर, जाल शोधन के विभिन्न स्तरों और उच्च-प्रदर्शन कंप्यूटिंग के लिए अनुकूलित है। इसे L2,H1, H-div और H-कर्ल स्थानों के लिए सन्निकटन के विषम क्रम से संबंधित जटिलताओं का प्रबंधन करने में सक्षम होने के लिए डिज़ाइन किया गया है।
* [http://mofem.eng.gla.ac.uk एमओएफईएम] परिमित अवयव विश्लेषण कोड है जो बहु-भौतिकी समस्याओं के समाधान के लिए अनेैतिक रूप से अनुमान के स्तर, जालक शोधन के विभिन्न स्तरों और उच्च-प्रदर्शन कंप्यूटिंग के लिए अनुकूलित किया गया है। और इसे L2,H1, H-डीआईवी और H-कर्ल स्थानों के लिए सन्निकटन के विषम क्रम से संबंधित सम्मिश्रों का प्रबंधन करने में सक्षम होने के लिए डिज़ाइन किया गया है।
* [http://www.sparselizard.org स्पार्सेलिज़ार्ड] बहु-भौतिकी, एचपी-अनुकूली, उपयोगकर्ता के अनुकूल, ओपन-सोर्स सी++ परिमित तत्व पुस्तकालय है जिसे वर्तमान में टाम्परे विश्वविद्यालय, फिनलैंड में विकसित किया गया है। यह सामान्य स्थैतिक और क्षणिक एचपी-एफईएम के लिए मनमाना क्रम पदानुक्रमित एच 1 और एच-कर्ल फ़ंक्शन रिक्त स्थान के साथ 3 डी टेट्राहेड्रल और 2 डी त्रिकोण / चतुर्भुज अनुरूप अनुकूली जाल शोधन को जोड़ता है।
* [http://www.sparselizard.org स्पार्सेलिज़ार्ड] बहु-भौतिकी, hp-अनुकूली, उपयोगकर्ता के अनुकूल, ओपन-सोर्स C++ परिमित अवयव पुस्तकालय है जिसे वर्तमान में टाम्परे विश्वविद्यालय, फिनलैंड में विकसित किया गया है। यह सामान्य स्थैतिक और क्षणिक hp-एफईएम के लिए इच्छानुसार क्रम पदानुक्रमित H1 और H-कर्ल फ़ंक्शन रिक्त स्थान के साथ 3 डी टेट्राहेड्रल और 2 डी त्रिकोण / चतुर्भुज अनुरूप अनुकूली जालक शोधन को जोड़ता है।


==वाणिज्यिक एचपी-एफईएम सॉफ्टवेयर==
==वाणिज्यिक एचपी-एफईएम सॉफ्टवेयर                                                                                           ==
* [https://www.esrd.com/ StressCheck] विस्तृत संरचनात्मक विश्लेषण की ओर उन्मुख एचपी-क्षमताओं वाला सीमित तत्व विश्लेषण उपकरण है।
* [https://www.esrd.com/ स्ट्रेसचेक] विस्तृत संरचनात्मक विश्लेषण की ओर उन्मुख hp-क्षमताओं वाला सीमित अवयव विश्लेषण उपकरण है।


==संदर्भ==
==संदर्भ==
Line 77: Line 63:


{{Numerical PDE}}
{{Numerical PDE}}
[[Category: सीमित तत्व विधि]]


[[Category: Machine Translated Page]]
[[Category:Collapse templates]]
[[Category:Created On 23/07/2023]]
[[Category:Created On 23/07/2023]]
[[Category:Machine Translated Page]]
[[Category:Navigational boxes| ]]
[[Category:Navigational boxes without horizontal lists]]
[[Category:Pages with script errors]]
[[Category:Sidebars with styles needing conversion]]
[[Category:Template documentation pages|Documentation/doc]]
[[Category:Templates Vigyan Ready]]
[[Category:Templates generating microformats]]
[[Category:Templates that are not mobile friendly]]
[[Category:Templates using TemplateData]]
[[Category:Wikipedia metatemplates]]
[[Category:सीमित तत्व विधि]]

Latest revision as of 17:08, 1 August 2023

एचपी-एफईएम परिमित अवयव विधि (एफईएम) का सामान्य संस्करण होता है | जो खंड अनुसार-बहुपद सन्निकटन के आधार पर आंशिक अंतर समीकरणों का समाधान करने के लिए संख्यात्मक विश्लेषण विधि होती है | और जो वैरिएबल आकार (h) और बहुपद की डिग्री (p) के अवयवों को नियोजित करता है। एचपी-एफईएम की उत्पत्ति बार्ना A सज़ाबो और इवो बाबूस्का के अग्रणी कार्य से हुई है | [1] [2] [3] [4] [5] [6] जिन्होंने इसमें पाया कि परिमित अवयव विधि शीघ्रता से परिवर्तित होती है। और जालक को h-शोधन (अवयवों को लघु भागों में विभाजित करना) होता हैं | इस प्रकार p-शोधन (उनकी बहुपद डिग्री को बढ़ाने) के उपयुक्त संयोजन का उपयोग करके परिष्कृत किया जाता है। यह घातीय अभिसरण अधिकांश अन्य परिमित अवयव विधियों की तुलना में विधि को बहुत आकर्षक बनाता है | जिससे यह सिर्फ बीजगणितीय दर के साथ अभिसरण करता है। एचपी-एफईएम के घातीय अभिसरण का पूर्वानुमान न सिर्फ सैद्धांतिक रूप से किया गया था, किंतु यह अनेक स्वतंत्र शोधकर्ताओं द्वारा भी देखा गया था। [7] [8] [9]

मानक एफईएम से अंतर

एचपी-एफईएम अनेक पहलुओं में मानक (निम्नतम-क्रम) एफईएम से भिन्न होते है।[10]

  • उच्च-क्रम आकार कार्यों का चयन उदाहरण आवश्यक: अवयवों में उच्च-डिग्री बहुपद को आकार कार्यों के विभिन्न समुच्चयो का उपयोग करके उत्पन्न किया जा सकता है। ऐसे समुच्चय का चयन कठोरता आव्युह की कंडीशनिंग और उसी स्थान में संपूर्ण समाधान प्रक्रिया को नाटकीय रूप से प्रभावित कर सकता है। इस समस्या को सबसे पहले बाबुस्का एट अल द्वारा प्रलेखित किया गया था।[11]
  • स्वचालित एचपी-अनुकूलन: एचपी-एफईएम में, अवयव को अनेक भिन्न-भिन्न विधियों से एचपी-परिष्कृत किया जा सकता है, जैसे: इसे स्पेस में उप-विभाजित किए बिना इसकी बहुपद डिग्री बढ़ाना, या अवयव को ज्यामितीय रूप से उप-विभाजित करना होता हैं | जहां विभिन्न बहुपद डिग्री को उप-अवयवों पर प्रयुक्त किया जा सकता है। और अवयव शोधन प्रत्याशी की संख्या सरलता से दो आयामों में 100 और तीन आयामों में 1000 तक पहुंच जाती है। किसी अवयव में त्रुटि के आकार को सांकेतिक करने वाली संख्या स्वचालित hp-अनुकूलता को निर्देशित करने के लिए पर्याप्त नहीं होती है | यह (मानक एफईएम में अनुकूलता के विपरीत) होती हैं । इस प्रकार यह प्रत्येक अवयव में त्रुटि के आकार के बारे में अधिक सूचना प्राप्त करने के लिए संदर्भ समाधान या विश्लेषणात्मक विचार जैसी अन्य तकनीकों को नियोजित किया जाना चाहिए। [12]
  • संयोजन और समाधान सीपीयू समय का अनुपात: मानक एफईएम में, कठोरता आव्युह होता हैं | यह सामान्यतः शीघ्रता से एकत्रित किया जाता है किन्तु यह अधिक विस्तृत होता है। यह सामान्यतः, असतत समस्या के समाधान में कुल कंप्यूटिंग समय का सबसे विस्तृत भाग व्यय होता है। और इसके विपरीत, hp- में एफईएम में कठोरता आव्युह सामान्यतः बहुत लघु होते हैं | किन्तु यह (समान आव्युह आकार के लिए) होता हैं और उनमें एकत्रित मानक एफईएम की तुलना में अधिक समय लगता है। यह मुख्य रूप से संख्यात्मक चतुर्भुज की कम्प्यूटेशनल निवेश के कारण होता है | जिसमें शीघ्र अभिसरण दरों का लाभ उठाने के लिए मानक एफईएम की तुलना में उच्च परिशुद्धता होनी चाहिए | और इसलिए यह उच्च क्रम का होना चाहिए।
  • विश्लेषणात्मक चुनौतियाँ: hp-एफईएम को सामान्यतः मानक एफईएम की तुलना में विश्लेषणात्मक दृष्टिकोण से समझना अधिक कठिन माना जाता है। जिसके अनुसार यह अनेक तकनीकों से संबंधित होता है, जैसे वृत्ताकार समस्याओं के लिए असतत अधिकतम सिद्धांत (डीएमपी) होता हैं। यह परिणाम बताते हैं कि, सामान्यतः जालक पर कुछ सीमित धारणाओं के साथ, खंड अनुसार-बहुपद एफईएम सन्निकटन अंतर्निहित वृत्ताकार पीडीई के समान अधिकतम सिद्धांत का पालन करता है। ऐसे परिणाम बहुत महत्वपूर्ण होते हैं क्योंकि वे गारंटी देते हैं कि सन्निकटन भौतिक रूप से स्वीकार्य रहता है | जिससे नकारात्मक घनत्व, नकारात्मक एकाग्रता, या नकारात्मक निरपेक्ष तापमान की गणना करने की कोई संभावना नहीं रहती है। डीएमपी निम्नतम-क्रम एफईएम के लिए अधिक अच्छी तरह से समझा जाता है किन्तु दो या दो से अधिक आयामों में hp-एफईएम के लिए पूर्ण तरह से अज्ञात होता है। इस प्रकार यह स्थानिक आयाम में प्रथम डीएमपी वर्तमान में तैयार किया गया था। [13] [14]
  • प्रोग्रामिंग चुनौतियाँ: मानक एफईएम कोड की तुलना में hp-एफईएम सॉल्वर को प्रयुक्त करना बहुत कठिन होता है। जिनमें अनेक विवादों को दूर करने की आवश्यकता होती है | यह उनमें सम्मिलित होता हैं (किन्तु यह सिर्फ यहीं तक सीमित नहीं होता हैं) | उच्च-क्रम चतुर्भुज सूत्र, उच्च-क्रम आकार फ़ंक्शन, भौतिक डोमेन में आधार कार्यों के साथ संदर्भ डोमेन पर आकार कार्यों से संबंधित कनेक्टिविटी और अभिविन्यास सूचना आदि होते हैं। [15]


फ़िचेरा समस्या

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

  • समस्या फ़ाइल: एकवचन ग्रेडिएंट.

    समस्या फ़ाइल: एकवचन ग्रेडिएंट.

  • समस्या चार्ट: अभिसरण तुलना।

    समस्या चार्ट: अभिसरण तुलना।

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

एचपी-एफईएम की दक्षता

लघु -रेखीय अवयवों की तुलना में विस्तृत उच्च-क्रम वाले अवयवों का उपयोग करके सुचारू कार्यों का अधिक कुशलता से अनुमान लगाया जा सकता है। इसे नीचे दिए गए चित्र में दर्शाया गया है | जहां दो भिन्न-भिन्न जालकों पर शून्य डिरिचलेट सीमा स्थितियों के साथ आयामी पॉइसन समीकरण समाधान किया गया है। यह स्पष्ट समाधान साइन फ़ंक्शन होता है।

  • बाएँ: दो रैखिक अवयवों से युक्त जालक हैं।
  • दाएँ: द्विघात अवयव से युक्त जालक हैं।

टुकड़े-टुकड़े-रैखिक सन्निकटन।द्विघात सन्निकटन.

जबकि दोनों स्तिथियों में अज्ञात की संख्या समान है यह (1 डीओएफ), संबंधित मानदंड में त्रुटियां क्रमशः 0.68 और 0.20 हैं। इसका कारण यह है कि द्विघात सन्निकटन खंड-रेखीय सन्निकटन की तुलना में लगभग 3.5 गुना अधिक कुशल था। जब हम कदम आगे बढ़ाते हैं और (a) चार रैखिक अवयवों की तुलना (b) चतुर्थक अवयव (p=4) से करते हैं, तब दोनों भिन्न-भिन्न समस्याओं में तीन डीओएफ होंते हैं | किन्तु चतुर्थक सन्निकटन लगभग 40 गुना अधिक कुशल होता हैं।

इसके विपरीत, लघु निम्न-क्रम वाले अवयव विस्तृत उच्च-क्रम वाले अवयवों की तुलना में लघु मापदंड की विशेषताओं जैसे विलक्षणताओं को उत्तम विधियों से पकड़ सकते हैं। hp-एफईएम इन दो दृष्टिकोणों के इष्टतम संयोजन पर आधारित होता है जो घातांकीय अभिसरण की ओर ले जाता है। ध्यान दें कि यह घातीय अभिसरण त्रुटि की धुरी और स्वतंत्रता की डिग्री में व्यक्त किया गया है। वास्तविक जीवन के अनुप्रयोगों के लिए, हम सामान्यतः स्पष्टता के समान स्तर तक पहुंचने के लिए आवश्यक कम्प्यूटेशनल समय पर विचार करते हैं। इस प्रदर्शन संकेतक के लिए h- और hp-शोधन समान परिणाम प्रदान कर सकते हैं | उदाहरण के लिए [16] (वेब आर्काइव लिंक [17]) पर अंतिम आंकड़ा देखते हैं | जैसे ही h-एफईएम की तुलना में hp-एफईएम को प्रोग्राम करना और समानांतर कंप्यूटिंग करना कठिन हो जाता है | और hp-शोधन की अभिसरण उत्कृष्टता अव्यावहारिक हो सकती है।

एचपी-अनुकूलन

कुछ एफईएम साइटें एचपी-अनुकूलता को h-अनुकूलता (उनकी बहुपद डिग्री को स्थिर रखते हुए स्पेस में अवयवों को विभाजित करता) हैं | और p-अनुकूलता (सिर्फ उनकी बहुपद डिग्री को बढ़ाना) के संयोजन के रूप में वर्णित करती हैं। यह पूर्णता से स्पष्ट नहीं होती है | क्योंकि hp-अनुकूलता h- और p-अनुकूलता दोनों से अधिक भिन्न होती है क्योंकि किसी अवयव का hp-शोधन अनेक भिन्न-भिन्न विधियों से किया जा सकता है। p-शोधन के अतिरिक्त, अवयव को स्पेस में उप-विभाजित किया जा सकता है (जैसा कि h-अनुकूलता में) हैं, किन्तु उप-अवयवों पर बहुपद डिग्री के लिए अनेक संयोजन होते हैं। यह दाहिनी ओर के चित्र में दर्शाया गया है। उदाहरण के लिए, यदि त्रिकोणीय या चतुर्भुज अवयव को चार उप-अवयवों में विभाजित किया जाता है | जहां बहुपद डिग्री को अधिकतम दो तक भिन्न होने की अनुमति होती है | तब इससे 3^4 = 81 शोधन प्रत्याशी मिलते हैं | इसमें (बहुपद अनिसोट्रोपिक प्रत्याशी पर विचार नहीं किया जाता है)। इस रूप से, हेक्साहेड्रोन को आठ उप-अवयवों में विभाजित करना होता हैं और अधिकतम दो द्वारा उनकी बहुपद डिग्री को परिवर्तित 3^8 = 6,561 शोधन प्रत्याशी प्राप्त करता है। यह प्रति अवयव स्थिर संख्या प्रदान करने वाला मानक एफईएम त्रुटि अनुमान स्वचालित hp-अनुकूलन का मार्गदर्शन करने के लिए पर्याप्त नहीं होता है।

उच्च-क्रम आकार के कार्य

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

  • वर्टेक्स फ़ंक्शन.

    वर्टेक्स फ़ंक्शन.

  • किनारे का कार्य।

    किनारे का कार्य।

  • चेहरे का कार्य.

    चेहरे का कार्य.

  • बुलबुला समारोह.

    बुलबुला समारोह.

ध्यान दें: यह सभी फ़ंक्शन संपूर्ण अवयव इंटीरियर में परिभाषित होते हैं।

ओपन सोर्स एचपी-एफईएम कोड

  • डील.II: डील.II परिमित अवयव विधि का उपयोग करके आंशिक अंतर समीकरणों का समाधान करने के लिए निःशुल्क, ओपन-सोर्स लाइब्रेरी है।
  • अवधारणाएं: एसएएम, ईटीएच ज्यूरिख (स्विट्जरलैंड) और टीयू बर्लिन (जर्मनी) में के. श्मिट के समूह में वृत्ताकार समीकरणों के लिए C/C++ hp-एफईएम/डीजीएफईएम/बीईएम लाइब्रेरी विकसित की गई हैं।
  • 2dhp90, 3dhp90: वृत्ताकार समस्याओं और मैक्सवेल के समीकरणों के लिए फोरट्रान कोड आईसीईएस, यूटी ऑस्टिन में एल. डेमकोविज़ द्वारा विकसित होता हैं।
  • पीएचएएमएल: समानांतर पदानुक्रमित अनुकूल बहु-स्तरीय परियोजना हैं। जिसमे अनुकूल जालक शोधन और मल्टी-ग्रिड समाधान तकनीकों का उपयोग करके वितरित मेमोरी समानांतर कंप्यूटर और मल्टी-कोर कंप्यूटर पर 2 डी वृत्ताकार आंशिक अंतर समीकरणों के संख्यात्मक समाधान के लिए,राष्ट्रीय मानक और प्रौद्योगिकी संस्थान, संयुक्त राज्य अमेरिका में परिमित तत्व सॉफ्टवेयर विकसित किया गया है।
  • हर्मीस परियोजना: पीडीई और मल्टीफिजिक्स पीडीई सिस्टम की विशाल विविधता के लिए स्पेस और स्पेस-समय अनुकूल hp-एफईएम सॉल्वरों के शीघ्रता से प्रोटोटाइप के लिए C/C++/पायथन लाइब्रेरी, नेवादा विश्वविद्यालय, रेनो (यूएसए), थर्मो-मैकेनिक्स संस्थान, प्राग (चेक गणराज्य) और पिल्सेन (चेक गणराज्य) में वेस्ट बोहेमिया विश्वविद्यालय में hp-एफईएम समूह द्वारा विकसित - एग्रोस2डी इंजीनियरिंग सॉफ्टवेयर के शीर्ष पर निर्मित हर्मीस पुस्तकालय हैं |
  • पीएचजी: पीएचजी समानांतर अनुकूल परिमित अवयव प्रोग्राम विकसित करने के लिए टूलबॉक्स होता है। यह h-, p- और hp-फेम के लिए उपयुक्त होता है। और पीएचजी वर्तमान में वैज्ञानिक और इंजीनियरिंग कंप्यूटिंग की स्थान प्रमुख प्रयोगशाला हैं | कम्प्यूटेशनल गणित संस्थान और चीनी विज्ञान अकादमी (एलएसईसी, सीएएस, चीन) के वैज्ञानिक/इंजीनियरिंग कंप्यूटिंग संस्थान में सक्रिय विकास पर निर्भर होती है। पीएचजी अनुरूप टेट्राहेड्रल जालक से संबंधित होता है | और संदेश भेजने के लिए अनुकूल स्थानीय जालक शोधन और एमपीआई के लिए द्विभाजन का उपयोग करता है। पीएचजी में ऑब्जेक्ट-ओरिएंटेड डिज़ाइन होता है जो समानांतर विवरण गुप्त रखता है | और अमूर्त विधियों से मेष और परिमित अवयव कार्यों पर सामान्य संचालन प्रदान करता है, जिससे उपयोगकर्ताओं को अपने संख्यात्मक एल्गोरिदम पर ध्यान केंद्रित करने की अनुमति मिलती है।
  • एमओएफईएम परिमित अवयव विश्लेषण कोड है जो बहु-भौतिकी समस्याओं के समाधान के लिए अनेैतिक रूप से अनुमान के स्तर, जालक शोधन के विभिन्न स्तरों और उच्च-प्रदर्शन कंप्यूटिंग के लिए अनुकूलित किया गया है। और इसे L2,H1, H-डीआईवी और H-कर्ल स्थानों के लिए सन्निकटन के विषम क्रम से संबंधित सम्मिश्रों का प्रबंधन करने में सक्षम होने के लिए डिज़ाइन किया गया है।
  • स्पार्सेलिज़ार्ड बहु-भौतिकी, hp-अनुकूली, उपयोगकर्ता के अनुकूल, ओपन-सोर्स C++ परिमित अवयव पुस्तकालय है जिसे वर्तमान में टाम्परे विश्वविद्यालय, फिनलैंड में विकसित किया गया है। यह सामान्य स्थैतिक और क्षणिक hp-एफईएम के लिए इच्छानुसार क्रम पदानुक्रमित H1 और H-कर्ल फ़ंक्शन रिक्त स्थान के साथ 3 डी टेट्राहेड्रल और 2 डी त्रिकोण / चतुर्भुज अनुरूप अनुकूली जालक शोधन को जोड़ता है।

वाणिज्यिक एचपी-एफईएम सॉफ्टवेयर

  • स्ट्रेसचेक विस्तृत संरचनात्मक विश्लेषण की ओर उन्मुख hp-क्षमताओं वाला सीमित अवयव विश्लेषण उपकरण है।

संदर्भ

  1. ↑ B. A. Szabó, A. K. Mehta: p-Convergent Finite Element Approximations in Fracture Mechanics, Int. J. Num. Meth. Engng, Volume 12, pp. 551-560, 1978.
  2. ↑ I. Babuška, B. A. Szabó and I. N. Katz: The p-Version of the Finite Element Method, SIAM J. Numer. Anl., Volume 18, pp. 515-544, 1981.
  3. ↑ I. Babuška, B. A. Szabó, On the Rates of Convergence of the Finite Element Method, Int. J. Numer. Meth.Engng., Volume 18, pp. 323-341, 1982.
  4. ↑ I. Babuška: The p- and hp-Versions of the Finite Element Method: the State of the Art, Finite Elements: Theory and Applications, edited by D. L. Dwoyer, M. Y. Hussaini and R. G. Voigt, New York, Springer-Verlag, 1988.
  5. ↑ B. A. Szabó, I. Babuška: Finite Element Analysis, John Wiley & Sons, ISBN 978-0-471-50273-9, 1991.
  6. ↑ I. Babuška, B.Q. Guo: The h, p and h-p version of the finite element method: basis theory and applications, Advances in Engineering Software, Volume 15, Issue 3-4, 1992.
  7. ↑ J.M. Melenk: hp-Finite Element Methods for Singular Perturbations, Springer, 2002
  8. ↑ C. Schwab: p- and hp- Finite Element Methods: Theory and Applications in Solid and Fluid Mechanics, Oxford University Press, 1998
  9. ↑ P. Solin: Partial Differential Equations and the Finite Element Method, J. Wiley & Sons, 2005
  10. ↑ P. Solin, K. Segeth, I. Dolezel: Higher-Order Finite Element Methods, Chapman & Hall/CRC Press, 2003
  11. ↑ I. Babuska, M. Griebel and J. Pitkaranta, The problem of selecting the shape functions for a p-type finite element, Internat. J. Numer. Methods Engrg. (1989), pp. 1891–1908
  12. ↑ L. Demkowicz, W. Rachowicz, and Ph. Devloo: A Fully Automatic hp-Adaptivity, Journal of Scientific Computing, 17, Nos 1–3 (2002), 127–155
  13. ↑ P. Solin, T. Vejchodsky: A Weak Discrete Maximum Principle for hp-FEM, J. Comput. Appl. Math. 209 (2007) 54–65
  14. ↑ T. Vejchodsky, P. Solin: Discrete Maximum Principle for Higher-Order Finite Elements in 1D, Math. Comput. 76 (2007), 1833–1846
  15. ↑ L. Demkowicz, J. Kurtz, D. Pardo, W. Rachowicz, M. Paszynski, A. Zdunek: Computing with hp-Adaptive Finite Elements, Chapman & Hall/CRC Press, 2007
  16. ↑ "Microwave Oven — Hermes Examples Guide".
  17. ↑ "Microwave Oven — Hermes Examples Guide". hpfem.org. Archived from the original on 7 August 2018. Retrieved 12 January 2022.
Retrieved from ‘https://www.vigyanwiki.in/index.php?title=एचपी-एफईएम&oldid=232879’