High-Resolution Weighted Essentially Non-Oscillatory Compact Least-Squares Schemes with Implicit Time Integration for Compressible Navier-Stokes Equations on Curvilinear Grids

कल्पना कीजिए कि आप कंप्यूटर पर एक तूफानी समुद्र का एक आदर्श चित्र बनाने की कोशिश कर रहे हैं। आप लहरों की कोमल, लुढ़कती हुई गति को अविश्वसनीय विवरण के साथ पकड़ना चाहते हैं, लेकिन आपको उस तीखे, ऊबड़-खाबड़ किनारे को भी बनाना है जहाँ एक विशाल लहर चट्टान से टकराती है। यह उन वैज्ञानिकों की दैनिक चुनौती है जो यह सिम्युलेट करते हैं कि हवा और गैस उच्च गति पर कैसे चलती है, जैसे कि एक सुपरसोनिक जेट के ऊपर से गुजरने वाली हवा या एक रॉकेट इंजन से निकलने वाली दहकती हुई गैस। इस क्षेत्र को कम्प्यूटेशनल फ्लूइड डायनेमिक्स (computational fluid dynamics) कहा जाता है। पेचीदा बात यह है कि कंप्यूटर उन कलाकारों की तरह हैं जो चिकनी वक्र रेखाएं (smooth curves) बनाने में तो माहिर हैं लेकिन तीखे कोनों को बनाने में बहुत खराब हैं; जब वे एक चिकनी रेखा को अचानक रुकने के लिए मजबूर करते हैं, तो वे अक्सर हिलने-डुलने और कंपन करने लगते हैं, जिससे एक अस्त-व्यस्त, टेढ़ा-मेढ़ा दृश्य बन जाता है जो पूरे चित्र को खराब कर देता है। दशकों से, वैज्ञानिक एक ऐसा "जादुई ब्रश" आविष्कार करने की कोशिश कर रहे हैं जो चिकनी लहरों और तीखी चट्टानों दोनों को बना सके बिना कंप्यूटर को भ्रमित किए या उसे हिलाए बिना।

इस शोध पत्र में, लेखकों ने, जिनका नेतृत्व निंगबो विश्वविद्यालय और सिंघुआ विश्वविद्यालय के शोधकर्ताओं द्वारा किया गया है, एक नए प्रकार के "जादुई ब्रश" को विकसित किया है। वे इसे "वेटेड कॉम्पैक्ट लीस्ट-स्क्वायर्स" (Weighted Compact Least-Squares) स्कीम कहते हैं, लेकिन आइए इसे "स्मार्ट स्केचर" (Smart Sketcher) कहें। पुराना तरीका ऐसा था जैसे किसी मुड़े हुए कागज को एक साथ दबाकर उसे चिकना करने की कोशिश करना; इससे अक्सर चिकने हिस्से लहरदार हो जाते थे और तीखे हिस्से धुंधले पड़ जाते थे। नया तरीका अधिक समझदार है। पूरे चित्र को एक साथ चिकना करने के बजाय, स्मार्ट स्केचर कागज को देखता है और कहता है, "ठीक है, यह हिस्सा एक चिकनी लहर है, मैं इसे उच्च सटीकता के साथ बनाऊंगा। लेकिन यह हिस्सा एक तीखी चट्टान है जहाँ हवा टकराती है, इसलिए मैं यहाँ एक चिकनी रेखा को जबरदस्ती थोपने की कोशिश करना बंद कर दूंगा और बस एक साफ, तीखा किनारा बनाऊंगा।" ऐसा करके, कंप्यूटर उस बिखराव और कंपन से बच जाता है और हवा की सूक्ष्म लहरों और शॉकवेव्स (shockwaves) के तीखे किनारों दोनों को स्पष्ट रूप से देख पाता है।

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

तो, उन्होंने वास्तव में क्या किया? उन्होंने एक नया गणितीय नुस्खा बनाया कि कंप्यूटर हवा के आकार को कैसे "पुनर्गठित" (reconstruct) करता है। अतीत में, कंप्यूटर पूरे ग्रिड में एक एकल, चिकनी गणितीय वक्र रेखा फिट करने की कोशिश करते थे, भले ही उस वक्र को अचानक आने वाली शॉकवेव के ऊपर से कूदना पड़ता था। बिना त्रुटियों के ऐसा करना गणितीय रूप से असंभव है, ठीक वैसे ही जैसे एक सीधी रेखा खींचने की कोशिश करना जो अचानक एक समकोण में बदल जाए बिना पेन उठाए; पेन फिसल जाएगा और लकीरें बना देगा। लेखकों ने इसे बदलकर खेल के नियम बदल दिए। उन्होंने एक ऐसी प्रणाली पेश की जो स्वचालित रूप से पता लगा लेती है कि चिकनी हवा कहाँ समाप्त होती है और शॉक कहाँ शुरू होता है। जब कंप्यूटर एक शॉक देखता है, तो वह एक चिकना संबंध बनाने की कोशिश करना बंद कर देता है और इसके बजाय एक विशेष "पेनल्टी" (penalty) लागू करता है जो चित्र को स्थिर और गैर-दोलनशील (non-oscillatory) रखता है।

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

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

अपने सिमुलेशन में, लेखकों ने अपने नए स्मार्ट स्केचर की तुलना WENO (Weighted Essentially Non-Oscillatory) जैसे पुराने, प्रसिद्ध तरीकों से की। परिणाम स्पष्ट थे: नया तरीका "हेल्महोल्ट्ज़ अस्थिरता" (Helmholtz instabilities)—वे सूक्ष्म लहरें जो तब बनती हैं जब हवा की दो परतें एक-दूसरे के बगल से गुजरती हैं—को बहुत बेहतर तरीके से पकड़ता है। एक उच्च-गति गैस जेट से जुड़े परीक्षण में, पुराने तरीके जेट के सिर के बारीक विवरणों को धुंधला करने लगे, जबकि नए तरीके ने संरचना को स्पष्ट और साफ रखा। यहाँ तक कि एक परीक्षण में जहाँ गैस मैक 2000 (ध्वनी की गति से बीस गुना) की गति से चल रही थी, नया तरीका टूटा नहीं या बेकाबू त्रुटियां पैदा नहीं की, जबकि कुछ पुराने तरीके पूरी तरह विफल हो गए।

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

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