निमज्जित सीमा विधि

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



\rho \left(\frac{\partial{u}({x},t)}{\partial{t}} + {u}\cdot\nabla{u}\right) = -\nabla p + \mu\, \Delta u(x,t) + f(x,t) $$ और यदि प्रवाह असम्पीडित होता है, तो हमारे पास आगे के नियम इस प्रकार है



\nabla \cdot u = 0. \, $$ निमज्जित संरचनाओं को सामान्यतः एक-आयामी फाइबर के संग्रह के रूप में प्रदर्शित किया जाता है, जिसे $$ \Gamma $$ द्वारा प्रदर्शित किया जाता है। प्रत्येक फाइबर को पैरामीट्रिक वक्र $$ X(s,t) $$ के रूप में देखा जा सकता है जहाँ $$ s $$ फाइबर के साथ लैग्रेंजियन समन्वय होता है और $$t $$ समय होता है। फाइबर की भौतिकी को फाइबर बल वितरण फ़ंक्शन $$ F(s,t) $$ के माध्यम से प्रदर्शित किया जाता है। स्प्रिंग बल, बेन्डिंग प्रतिरोध या किसी अन्य प्रकार का व्यवहार इस शब्द में निर्मित किया जा सकता है। द्रव पर संरचना द्वारा लगाए गए बल को फिर संवेग समीकरण में स्रोत शब्द के रूप में प्रक्षेपित किया जाता है



f(x,t) = \int_\Gamma F(s,t) \, \delta\big(x - X(s,t)\big) \, ds, $$ जहाँ $$ \delta $$ डिराक डेल्टा फ़ंक्शन है। इलास्टिक सतहों या त्रि-आयामी ठोस पदार्थों को प्रतिरूप करने के लिए बल को कई आयामों तक बढ़ाया जा सकता है। एक द्रव्यमान रहित संरचना मानते हुए, इलास्टिक फाइबर स्थानीय द्रव वेग के साथ चलता है और डेल्टा फ़ंक्शन के माध्यम से प्रक्षेपित किया जा सकता है



\frac{\partial X(s,t)}{\partial t} = u(X,t) = \int_\Omega u(x,t) \, \delta\big(x - X(s,t)\big) \, dx, $$ जहाँ $$ \Omega $$ संपूर्ण द्रव कार्यक्षेत्र को प्रदर्शित करता है। इन समीकरणों का डिसक्रेटीजेशन द्रव पर एक यूलेरियन ग्रिड और फाइबर पर एक अलग लैग्रेंजियन ग्रिड मानकर किया जा सकता है। सुचारू कार्यों द्वारा डेल्टा वितरण का प्राक्कलन हमें दो ग्रिडों के मध्य अंतरण करने की अनुमति देता है। निमज्जित सीमा समीकरणों को हल करने के लिए किसी भी उपस्थित द्रव सॉल्वर को फाइबर समीकरणों के सॉल्वर के साथ जोड़ा जा सकता है। इस मूलभूत दृष्टिकोण के भिन्न रूप को इलास्टिक संरचनाओं से युक्त विभिन्न प्रकार की यांत्रिक प्रणालियों को अनुकरण करने के लिए प्रयुक्त किया जाता है जो द्रव प्रवाह के साथ अन्तःक्रिया करते हैं।

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

यह भी देखें

 * स्टोकेस्टिक यूलेरियन लैग्रेंजियन विधि
 * स्टोकेशियन गतिकी
 * तरल पदार्थ की मात्रा विधि
 * लेवल-सेट विधि
 * मार्कर-और-सेल विधि

सॉफ्टवेयर: न्यूमेरिकल कोड

 * FloEFD: वाणिज्यिक सीएफडी आईबीएम कोड
 * उन्नत सिमुलेशन लाइब्रेरी
 * मैंगो-सेल्म: निमज्जित सीमा विधियाँ और एसईएलएम सिमुलेशन, 3डी पैकेज, (पायथन इंटरफ़ेस, एलएएमएमपीएस एमडी इंटीग्रेशन), पी. एट्ज़बर्गर, यूसीएसबी
 * 3डी में स्टोकेस्टिक निमज्जित सीमा विधि, पी. एट्ज़ बर्गर, यूसीएसबी
 * 2डी में यूनिफार्म लैटिस के लिए निमज्जित सीमा विधि, ए. फोगेलसन, यूटा
 * IBAMR: 3डी में अडाप्टिव मेशेस के लिए निमज्जित सीमा विधि, बी. ग्रिफ़िथ, एनवाईयू।
 * IB2d: 60+ उदाहरणों के साथ 2डी में मैटलैब और पायथन के लिए निमज्जित सीमा विधि, एन.ए. बतिस्ता, टीसीएनजे
 * ESPResSo: सॉफ्ट इलास्टिक वस्तुओं के लिए निमज्जित सीमा विधि
 * OpenFoam पर आधारित CFD IBM कोड
 * sdfibm: ओपनफोम पर आधारित एक और CFD IBM कोड
 * सिमस्केल: क्लाउड में द्रव यांत्रिकी और संयुग्मी हीट हस्तांतरण ट्रान्सफर के लिए निमज्जित सीमा विधि