💻 computer science

Adaptation of the hybrid fictitious domain-immersed boundary method for Reynolds-averaged turbulence modeling

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

मूल लेखक: Lucie Kubíčková, Martin Isoz

प्रकाशित 2026-06-05

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Lucie Kubíčková, Martin Isoz

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

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

यह शोध पत्र इस तरह के सिमुलेशन का एक नया, स्मार्ट तरीका पेश करता है, जिसे openHFDIBRANS कहा जाता है। यह कैसे काम करता है, इसके लिए सरल उपमाओं (analogies) का उपयोग किया गया है:

1. "घोस्ट वॉल" ट्रिक (इमर्स्ड बाउंड्री)

हर नए आकार के लिए नया कागज का टुकड़ा (मेश) काटने के बजाय, शोधकर्ता पूरे क्षेत्र को कवर करने वाले एक विशाल, सरल ग्रिड (जैसे चेकरबोर्ड) का उपयोग करते हैं। वे ग्रिड को बदलते नहीं हैं। इसके बजाय, वे एक "जादुई मार्कर" का उपयोग करके कंप्यूटर को बताते हैं: "हे, ये विशिष्ट वर्ग ठोस दीवारें हैं, और ये खाली हवा हैं।"

यदि आप वस्तु का आकार बदलना चाहते हैं, तो आप बस ग्रिड पर "जादुई मार्कर" को स्थानांतरित करते हैं। आपको पूरा ग्रिड फिर से बनाने की आवश्यकता नहीं है। यह एक डिजिटल सैंडबॉक्स के साथ खेलने जैसा है जहाँ आप पूरे रेत के टेबल को फिर से बनाए बिना तुरंत एक महल का आकार बदल सकते हैं।

2. "ट्रैफिक पुलिस" की समस्या (टर्बुलेंस)

जटिल हिस्सा यह है कि तरल पदार्थ (जैसे हवा या पानी) अव्यवस्थित और अशांत (turbulent) हो जाते हैं, खासकर दीवारों के पास। वास्तविक दुनिया में, दीवार के ठीक बगल में हवा बहुत धीमी गति से चलती है, जबकि थोड़ी दूर की हवा तेज चलती है। इसे "बाउंड्री लेयर" कहा जाता है।

आमतौर पर, कंप्यूटरों को दीवार के ठीक बगल में बहुत विस्तृत, उच्च-रिज़ॉल्यूशन वाले मानचित्र की आवश्यकता होती है ताकि वे इन सूक्ष्म विवरणों को देख सकें। लेकिन चूंकि हमारा "जादुई मार्कर" ग्रिड सरल है और दीवार से पूरी तरह नहीं चिपकता है, इसलिए कंप्यूटर उन सूक्ष्म विवरणों को नहीं देख पाता है।

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

3. "रस्सी पर संतुलन" की चुनौती (स्थिरता)

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

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

4. क्या यह काम आया? (परिणाम)

टीम ने अपने नए टूल का कई क्लासिक चुनौतियों पर परीक्षण किया:

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

एक कमी:
कुछ मामलों में, नए टूल ने "वेक" (वस्तु के पीछे छोड़े गए अव्यवस्थित, घूमते हुए हवा के निशान) के आकार को थोड़ा "अधिक आंक लिया" (overestimated)। यह ऐसा है जैसे टूल ने सोचा कि वस्तु वास्तव में जितनी बड़ी है उससे थोड़ी बड़ी थी, जिससे अशांति का एक बड़ा निशान पीछे छूट गया। हालांकि, लेखकों का कहना है कि यह एक अनुमानित व्यवहार है और यह परिणामों को खराब नहीं करता है।

यह क्यों महत्वपूर्ण है?

इस कार्य का मुख्य लक्ष्य ऑप्टिमाइज़ेशन (अनुकूलन) है। क्योंकि इस नए टूल को हर बार आकार बदलने पर ग्रिड को फिर से बनाने की आवश्यकता नहीं होती है, इसलिए यह टोपोलॉजिकल ऑप्टिमाइज़ेशन के लिए एकदम सही है। यह एक ऐसी प्रक्रिया है जहाँ कंप्यूटर स्वचालित रूप से हजारों अलग-अलग आकारों को आज़माता है ताकि सबसे अच्छा आकार (जैसे, सबसे हल्का पुल या सबसे वायुगतिकीय कार का हिस्सा) पाया जा सके।

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

संक्षेप में: उन्होंने तरल प्रवाह के सिमुलेशन का एक तेज़, अधिक लचीला तरीका बनाया है जो आकारों को डिजाइन करने के लिए बहुत अच्छा काम करता है, बशर्ते आप प्रवाह के "पूंछ" (tail) में होने वाली एक छोटी, अनुमानित त्रुटि को स्वीकार करें।

तकनीकी सारांश: RANS मॉडलिंग के लिए हाइब्रिड फिक्टिशियस डोमेन-इमर्स्ड बाउंड्री विधि का अनुकूलन

समस्या विवरण
अभियांत्रिकी अनुप्रयोगों में अक्सर तरल घटकों के आकार या टोपोलॉजी अनुकूलन (TO) की आवश्यकता होती है, जहाँ कम्प्यूटेशनल फ्लूइड डायनेमिक्स (CFD) लागत फलन (cost function) मूल्यांकन के आधार के रूप में कार्य करता है। CFD-आधारित TO फ्रेमवर्क में एक महत्वपूर्ण बाधा बार-बार रीमेशिंग (remeshing) की आवश्यकता है क्योंकि ज्यामिति बदलती रहती है। जबकि इमर्स्ड बाउंड्री मेथड्स (IBMs) सरल, नॉन-कॉन्फॉर्मिंग मेश का उपयोग करके और स्केलर इंडिकेटर फील्ड के माध्यम से ज्यामिति को प्रदर्शित करके एक समाधान प्रदान करते हैं, उनका रेनॉल्ड्स-एवरेज्ड नेवियर-स्टोक्स (RANS) टर्बुलेंस मॉडल के साथ एकीकरण चुनौतीपूर्ण बना हुआ है। प्राथमिक कठिनाइयों में अनिसोट्रोपिक मेश रिफाइनमेंट के बिना बाउंड्री लेयर को मॉडल करना, उच्च रेनॉल्ड्स नंबर पर दबाव और वॉल शियर स्ट्रेस में स्प्यूरियस ऑसिलेशन (spurious oscillations) को रोकना, और अस्थिरता से बचने के लिए मोमेंटम और टर्बुलेंस वेरिएबल्स के बीच एक कड़ा जुड़ाव सुनिश्चित करना शामिल है। इसके अलावा, मौजूदा IBM-RANS सॉल्वर अक्सर ओपन-सोर्स नहीं होते हैं या उनमें मानक CFD प्रथाओं का अभाव होता है।

कार्यप्रणाली
लेखक openHFDIBRANS का प्रस्ताव करते हैं, जो OpenFOAM लाइब्रेरी पर निर्मित एक ओपन-सोर्स सॉल्वर है जो स्थिर-अवस्था (steady-state), इनकम्प्रेसिबल RANS सिमुलेशन के लिए हाइब्रिड फिक्टिशियस डोमेन-इमर्स्ड बाउंड्री (HFDIB) विधि को अनुकूलित करता है। यह कार्यप्रणाली निम्नलिखित प्रमुख घटकों को एकीकृत करती है:

शासी समीकरण (Governing Equations): सॉल्वर चार टू-इक्वेशन टर्बुलेंस मॉडल्स: $k$ - $\epsilon$ , $k$ - $\omega$ , $k$ - $\omega$ SST, और रियलाइजेबल $k$ - $\epsilon$ के साथ युग्मित स्थिर-अवस्था RANS समीकरणों को संबोधित करता है। समीकरणों को एक फेज इंडिकेटर फील्ड $\lambda$ द्वारा संशोधित किया गया है ताकि समाधान को फ्लूइड डोमेन तक सीमित रखा जा सके और बाउंड्री कंडीशंस को लागू करने के लिए सोर्स टर्म्स ( $f_{ib}$ और $S_{ib}$ ) शामिल किए गए हैं।
मास्किंग और सेल वर्गीकरण: कम्प्यूटेशनल डोमेन को $\lambda$ फील्ड के आधार पर इन-सॉलिड सेल्स ( $\Omega^h_{ic}$ ), बाउंड्री सेल्स ( $\Omega^h_{bc}$ ), और फ्री-स्ट्रीम सेल्स ( $\Omega^h_{fc}$ ) में विभाजित किया गया है। मास्किंग फील्ड्स ( $\alpha_u, \alpha_k, \alpha_\omega$ ) प्रासंगिक क्षेत्रों में सोर्स टर्म्स के अनुप्रयोग को सीमित करते हैं।
वॉल मॉडलिंग और बाउंड्री ट्रीटमेंट: बाउंड्री लेयर को हल किए बिना संभालने के लिए, विधि वॉल फंक्शन्स का उपयोग करती है।
- फ्रिक्शन वेलोसिटी: मोमेंटम और टर्बुलेंस के बीच स्थिरता और कड़े जुड़ाव को सुनिश्चित करने के लिए, फ्रिक्शन वेलोसिटी ( $u_\tau$ ) की गणना स्थानीय वेग प्रोफाइल के बजाय टर्बुलेंस वेरिएबल्स का उपयोग करके की जाती है।
- इम्पोज़्ड फील्ड्स: उन बाउंड्री सेल्स में जहाँ बाउंड्री लेयर अनरिज़ॉल्व्ड है ( $y^+_{ib} > y^+_{lam}$ ), वेग को इंटरपोलेशन के माध्यम से पुनर्गठित नहीं किया जाता है, बल्कि केवल टर्बुलेंट विस्कोसिटी द्वारा सुधारा जाता है। रिज़ॉल्व्ड क्षेत्रों में ( $y^+_{ib} \leq y^+_{lam}$ ), वॉल फंक्शन्स के आधार पर वेग और टर्बुलेंस काइनेटिक एनर्जी ( $k$ ) को लागू करने के लिए पॉलिनोमियल इंटरपोलेशन का उपयोग किया जाता है।
- टर्बुलेंस वेरिएबल्स: विशिष्ट डिसिपेशन रेट ( $\omega$ ) के लिए विशिष्ट सोर्स टर्म्स नहीं जोड़े जाते हैं क्योंकि यह दीवारों के पास डाइवर्ज (diverge) हो जाता है। इसके बजाय, डिस्क्रीटाइजेशन मैट्रिक्स में हेरफेर करके प्रभावित सेल्स में $\omega$ और $k$ के प्रोडक्शन ( $P_k$ ) को सीधे असाइन किया जाता है, जो प्रभावी रूप से डिDirichlet-जैसे बाउंड्री कंडीशंस को लागू करता है।
समाधान रणनीति (Solution Strategy): सॉल्वर एक विस्तारित SIMPLE एल्गोरिदम का उपयोग करता है। महत्वपूर्ण रूप से, मोमेंटम प्रेडिक्टर से पहले, इम्पोज़्ड फील्ड्स की गणना की जाती है। मोमेंटम समीकरण को हल करने के बाद, बाउंड्री कंडीशंस को सीधे लागू करने के लिए एक करेक्शन स्टेप ( $u = u + \alpha_u(u_{ib} - u)$ ) लगाया जाता है। यह स्टेप स्थिर-अवस्था सिमुलेशन में अत्यंत महत्वपूर्ण है जहाँ अंडर-रिलैक्सेशन की आवश्यकता होती है। इनकम्प्रेसिबिलिटी कंडीशन को केवल फ्लूइड क्षेत्र में लागू किया जाता है ताकि IBM फोर्सिंग के साथ संघर्ष से बचा जा सके।

प्रमुख योगदान

ओपन-सोर्स कार्यान्वयन: यह कार्य व्यापक रूप से उपयोग किए जाने वाले OpenFOAM फ्रेमवर्क के भीतर एक IBM-RANS सॉल्वर का पूर्ण ओपन-सोर्स कार्यान्वयन प्रदान करता है, जिसमें मानक CFD प्रथाओं को शामिल किया गया है।
स्थिरता-उन्मुख डिज़ाइन: लेखक टर्बुलेंस वेरिएबल्स के माध्यम से फ्रिक्शन वेलोसिटी का मूल्यांकन करके और मोमेंटम एवं टर्बुलेंस ट्रांसपोर्ट समीकरणों के बीच एक कड़ा संबंध स्थापित करके सामान्य IBM-RANS अस्थिरताओं को संबोधित करते हैं।
अनुकूलन के लिए मजबूती: विधि को विशेष रूप से ज्यामिति अनुकूलन अनुप्रयोगों के लिए डिज़ाइन किया गया है, जो बिना रीमेशिंग के विभिन्न रेनॉल्ड्स नंबरों ( $10^1$ से $10^6$ ) और जटिल ज्यामिति में मजबूती प्रदर्शित करती है।
व्यापक बेंचमार्किंग: सॉल्वर को छह परिदृश्यों में मानक बॉडी-फिटेड सिमुलेशन (simpleFoam) और प्रयोगात्मक डेटा के विरुद्ध मान्य किया गया है: 2D/3D पाइप फ्लो, बैकवर्ड-फेसिंग स्टेप, सिलेंडर के ऊपर प्रवाह, NACA-0009 एयरफॉइल, और अहमद बॉडी।

परिणाम
openHFDIBRANS के प्रदर्शन का मूल्यांकन मानक बॉडी-फिटेड CFD (simpleFoam) और प्रयोगात्मक डेटा के विरुद्ध किया गया:

सटीकता: सॉल्वर ने विभिन्न रेनॉल्ड्स नंबरों और विभिन्न टर्बुलेंस मॉडल्स के तहत simpleFoam के अनुरूप परिणाम दिए। वेग, टर्बुलेंट विस्कोसिटी और रेनॉल्ड्स स्ट्रेस के लिए $L_2$ नॉर्म एरर आम तौर पर कम थे, हालांकि टर्बुलेंट विस्कोसिटी में अचानक बदलाव वाले क्षेत्रों (जैसे $k$ - $\omega$ SST और रियलाइजेबल $k$ - $\epsilon$ मॉडल्स) में थोड़े अधिक थे।
संवेदनशीलता: विधि ने मेश के सापेक्ष सॉलिड-फ्लुइड इंटरफ़ेस की स्थिति के संबंध में संतोषजनक मजबूती दिखाई। हालांकि मेश फेसेस के साथ संरेखण (alignment) ने सटीकता में सुधार किया, लेकिन सॉल्वर मिस-अलाइन्ड इंटरफेस के साथ भी कार्यात्मक रहा।
ड्रैग और लिफ्ट: सिलेंडर और NACA-0009 एयरफॉइल के लिए, सॉल्वर ने प्रवाह पैटर्न को कैप्चर किया और संदर्भ डेटा के साथ अच्छे सामंजस्य के साथ लिफ्ट और ड्रैग गुणांकों की भविष्यवाणी की।
सीमाएँ: विशिष्ट मामलों में (जैसे बैकवर्ड-फेसिंग स्टेप और NACA-0009), सॉल्वर ठोस वस्तुओं के पीछे के वेक (wake) क्षेत्र को थोड़ा अधिक आंकने की प्रवृत्ति रखता है। इसके अतिरिक्त, अन्य RAS मॉडल्स की तरह, यह सिलेंडर के लिए ड्रैग क्राइसिस को हल नहीं कर सका, जिसके लिए अधिक उन्नत टू-लेयर मॉडल्स की आवश्यकता होती है।

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

1. "घोस्ट वॉल" ट्रिक (इमर्स्ड बाउंड्री)

2. "ट्रैफिक पुलिस" की समस्या (टर्बुलेंस)

3. "रस्सी पर संतुलन" की चुनौती (स्थिरता)

4. क्या यह काम आया? (परिणाम)

यह क्यों महत्वपूर्ण है?

इसी तरह के और पेपर