🔢 mathematics

A generalized vertical coordinate transformation based on SPH(2) for efficient free surface flow simulations

यह शोध पत्र जटिल निचली सीमाओं वाले मुक्त सतह प्रवाह (free surface flows) के अनुकरण के लिए तीन कुशल कण विधियों—BF-SPH, BFE-SPH, और $\sigma$ -SPH—का प्रस्ताव करता है, जो सीमा स्थितियों को सटीक रूप से लागू करने और गणनात्मक प्रदर्शन को अनुकूलित करने के लिए द्वितीय-क्रम सटीक SPH(2) पर आधारित एक सामान्य ऊर्ध्वाधर समन्वय रूपांतरण (Vertical Coordinate Transformation) का उपयोग करते हैं।

मूल लेखक: Shujiro Fujioka, Kumpei Tsuji, Naoto Mitsume, Mitsuteru Asai

प्रकाशित 2026-08-03

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

CC BY 4.0

मूल लेखक: Shujiro Fujioka, Kumpei Tsuji, Naoto Mitsume, Mitsuteru Asai

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

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

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

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

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

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

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

यह सुनिश्चित करने के लिए कि ये आकार बदलने वाली ट्रिक्स गणित को खराब न करें, लेखकों को SPH(2) नामक एक बहुत ही सटीक गणना उपकरण का उपयोग करना पड़ा। इन कण सिमुलेशन के मानक तरीके थोड़े धुंधले लेंस की तरह हैं; वे सरल चीजों के लिए ठीक काम करते हैं लेकिन जब आप स्थान को खींचने की कोशिश करते हैं तो वे गड़बड़ा जाते हैं। SPH(2) विधि एक हाई-डेफिनिशन लेंस है जो सटीकता खोए बिना खींचने और सिकोड़ने के जटिल गणित को संभाल सकती है। शोधकर्ताओं ने इन नई विधियों का परीक्षण कई सिमुलेशन में किया, जिसमें पानी का एक शांत टैंक, एक त्रिकोणीय उभार के ऊपर टूटता हुआ बांध, और यहाँ तक कि एक दीवार से टकराते पानी का 3D सिमुलेशन शामिल था।

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

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

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

तकनीकी सारांश: कुशल मुक्त सतह प्रवाह सिमुलेशन के लिए सामान्यीकृत ऊर्ध्वाधर समन्वय रूपांतरण (Generalized Vertical Coordinate Transformation)

समस्या विवरण
कण विधियाँ (Particle methods), जैसे कि स्मूद्थ्ड पार्टिकल हाइड्रोडायनामिक्स (SPH) और मूविंग पार्टिकल सेमी-इम्प्लिसिट (MPS), बड़े विरूपण (deformations) और गतिशील विच्छेदन (moving discontinuities) वाले मुक्त सतह प्रवाह के अनुकरण के लिए प्रभावी हैं। हालाँकि, उन्हें महत्वपूर्ण कम्प्यूटेशनल लागत चुनौतियों का सामना करना पड़ता है। पारंपरिक कण विधियों को पूरे डोमेन में समान कण अंतराल (uniform particle spacing) की आवश्यकता होती है ताकि सटीकता बनी रहे, जिससे उन क्षेत्रों में भी बड़ी संख्या में कणों की आवश्यकता होती है जहाँ उच्च रिज़ॉल्यूशन महत्वपूर्ण नहीं है (जैसे सुनामी सिमुलेशन में गहरे समुद्री जल में)। जबकि मेश-आधारित विधियाँ लागत कम करने के लिए ग्रिड आकारों को समायोजित कर सकती हैं, पारंपरिक कण विधियों में यह लचीलापन नहीं होता है। इसके अलावा, जटिल बॉटम सीमाओं (जैसे समुद्री तल, नदी के तल) को संभालना अक्सर जटिल कण व्यवस्थाओं या सीमा उपचारों की मांग करता है जो सटीकता या स्थिरता को कम कर सकते हैं। लेखक एक ऐसी विधि की पहचान करते हैं जो कण विधियों के लैग्रेंजियन लाभों से समझौता किए बिना गैर-समान कण रिज़ॉल्यूशन (विशेष रूप से ऊर्ध्वाधर दिशा में) की अनुमति देकर कम्प्यूटेशनल दक्षता में सुधार करती है और सीमा हैंडलिंग को सरल बनाती है।

कार्यप्रणाली (Methodology)
यह शोध पत्र एक ढांचे का प्रस्ताव करता है जो SPH विधि पर लागू सामान्यीकृत ऊर्ध्वाधर समन्वय रूपांतरणों (VCTs) पर आधारित है। मुख्य अवधारणा भौतिक स्थान (जटिल सीमाओं और परिवर्तनशील गहराई के साथ) को एक सामान्यीकृत "प्रक्षेपित स्थान" (projected space) पर प्रक्षेपित करना है जहाँ ऊर्ध्वाधर समन्वय को रूपांतरित किया जाता है। यह प्रक्षेपित स्थान में दीर्घवृत्ताकार कणों (ellipsoidal particles) या सरलीकृत सीमा निरूपणों के उपयोग की अनुमति देता है।

इसका कार्यान्वयन तीन प्रमुख तकनीकी घटकों पर निर्भर करता है:

समन्वय रूपांतरण (Coordinate Transformation): एक रूपांतरण परिभाषित किया गया है जहाँ भौतिक ऊर्ध्धर समन्वय $z$ को स्केलिंग कारकों $\alpha(x,y)$ और ऊँचाई सुधार $\beta(x,y)$ का उपयोग करके एक प्रक्षेपित समन्वय $b_z$ में मैप किया जाता है। इस रूपांतरित स्थान में गवर्निंग समीकरणों को सटीक रूप से हल करने के लिए क्रॉस-डेरिवेटिव्स सहित द्वितीय-क्रम के डेरिवेटिव्स की गणना आवश्यक है।
उच्च-सटीक विविक्तीकरण (High-Accuracy Discretization - SPH(2)): समन्वय रूपांतरण द्वारा पेश किए गए जटिल डेरिवेटिव पदों को संभालने के लिए, लेखक SPH(2) का उपयोग करते हैं, जो एक द्वितीय-क्रम सटीक SPH मॉडल है। यह मॉडल उच्च सटीकता के साथ द्वितीय-क्रम के डेरिवेटिव्स (क्रॉस-डेरिवेटिव्स सहित) की गणना करता है, जो VCTs की स्थिरता और सटीकता के लिए आवश्यक है।
आयतन संरक्षण (Volume Conservation): चूंकि समन्वय रूपांतरण जैकोबियन डिटरमिनेंट (और इस प्रकार कणों के प्रतिनिधि आयतन) को बदल देता है जब $\alpha$ परिवर्तित होता है, इसलिए एक नया आयतन संरक्षण सुधार तकनीक पेश की गई है। इसमें जैकोबियन को स्मूथ करना और संख्यात्मक अस्थिरता और आयतन हानि को रोकने के लिए प्रेशर पॉइसन समीकरण में एक स्थिरीकरण पद (stabilization term) के माध्यम से कणों के आयतन को समायोजित करना शामिल है।

तीन प्रस्तावित विधियाँ
लेखक VCT अवधारणा को तीन विशिष्ट विधियों में सामान्यीकृत करते हैं:

BF-SPH (बॉटम बाउंड्री-फिटेड SPH): यह फाइनाइट डिफरेंस विधियों में बॉडी-फिटेड समन्वय प्रणालियों से प्रेरित है। यह जटिल बॉटम सीमाओं को प्रक्षेपित स्थान में सपाट सतहों में बदल देता है ( $\beta = -h(x,y)$ , $\alpha=1$ )। यह बॉटम बाउंड्री कंडीशंस लगाने को सरल बनाता है।
BFE-SPH (बॉटम बाउंड्री-फिटेड दीर्घवृत्ताकार SPH): यह BF-SPH को एक दीर्घवृत्ताकार कण मॉडल के साथ जोड़ता है। यह भौतिक स्थान में कणों को ऊर्ध्वाधर रूप से लंबा करने के लिए एक स्थिर ऊर्ध्वाधर स्केलिंग कारक ( $\alpha > 1$ ) का उपयोग करता है, जिससे बॉटम-फिटेड लाभ बनाए रखते हुए कुल कणों की संख्या कम हो जाती है।
$\sigma$ -SPH: यह एक $\sigma$ -समन्वय प्रणाली (जो समुद्र विज्ञान में सामान्य है, जैसे प्रिंसटन ओशन मॉडल) का उपयोग करता है जहाँ ऊर्ध्वाधर स्केलिंग कारक $\alpha$ जल की गहराई के आधार पर बदलता है ( $\alpha = \bar{H} / (H - h(x,y))$ )। यह रिज़ॉल्यूशन में एक स्टेपवाइज ट्रांजिशन की अनुमति देता है, जो उथले तटीय क्षेत्रों में उच्च रिज़ॉल्यूशन और गहरे समुद्री क्षेत्रों में मोटा (coarser) रिज़ॉल्यूशन प्रदान करता है, जिससे कणों की कुल संख्या में काफी कमी आती है।

मुख्य परिणाम
इन विधियों को हाइड्रोस्टेटिक दबाव, डैम ब्रेक समस्याओं और ढलानों पर गतिशील प्रवाह के 2-D और 3-D सिमुलेशन के माध्यम से मान्य किया गया:

सटीकता: VCTs के साथ SPH(2) के उपयोग ने मानक ISPH की तुलना में सटीकता में काफी सुधार किया। हाइड्रोस्टेटिक परीक्षणों में, SPH(2) के लिए दबाव में $L_2$ त्रुटि मानक ISPH की तुलना में दो क्रम (two orders) छोटी थी। इन विधियों ने जटिल बॉटम ज्योमेट्री और परिवर्तनशील आस्पेक्ट रेश्यो के साथ भी हाइड्रोस्टेटिक दबाव वितरण को सफलतापूर्वक पुनरुत्पादित किया।
दक्षता: $\sigma$ -SPH विधि ने उच्चतम दक्षता प्रदर्शित की। एक हाइड्रोस्टेटिक दबाव समस्या में, $\sigma$ -SPH ने मानक BF-SPH विधि की तुलना में कणों की संख्या को लगभग 1/4 और कम्प्यूटेशनल समय को 1/6 तक कम कर दिया, जबकि तुलनीय सटीकता बनाए रखी।
स्थिरता और आयतन संरक्षण: प्रस्तावित आयतन संरक्षण तकनीक (डायनेमिक समस्या में कंडीशन C) ने लंबे सिमुलेशन समय (100s तक) के दौरान कुल आयतन त्रुटियों को $\pm 1\%$ के भीतर सफलतापूर्वक बनाए रखा, जबकि बिना वॉल्यूम स्मूथिंग या गलत स्मूथिंग पैरामीटर्स वाली विधियों में आयतन की हानि या विचलन देखा गया।
3-D क्षमता: एक 3-D डैम ब्रेक सिमुलेशन ने पुष्टि की कि $\sigma$ -SPH विधि वेव रन-अप, रिफ्लेक्शन और वोर्टेक्स फॉर्मेशन सहित जटिल 3-D प्रवाह घटनाओं को घनत्व त्रुटियों और वेलोसिटी डायवर्जेंस को कम रखते हुए मजबूती से संभाल सकती है।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह कार्य कण विधि (particle method) में $\sigma$ -समन्वय प्रणाली को लागू करने का पहला प्रयास है। लेखक इस बात पर जोर देते हैं कि इन विधियों की सफलता समन्वय रूपांतरण द्वारा आवश्यक द्वितीय-क्रम के डेरिवेटिव्स को सटीक रूप से गणना करने के लिए SPH(2) के उपयोग पर अत्यधिक निर्भर है।

प्राथमिक महत्व निम्नलिखित क्षमताओं में निहित है:

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

लेखक निष्कर्ष निकालते हैं कि हालांकि यह अध्ययन मुख्य रूप से 3-D में इन विधियों की व्यवहार्यता की पुष्टि करता है, $\sigma$ -SPH विधि गैर-बायस्ड विधियों की तुलना में पांच गुना बेहतर कम्प्यूटेशनल दक्षता प्रदान करती है, जो इसे तटीय और समुद्री इंजीनियरिंग के व्यावहारिक अनुप्रयोगों के लिए एक आशाजनक उपकरण बनाती है। वे यह भी नोट करते हैं कि ऊर्ध्वाधर रिज़ॉल्यूशन को नियंत्रित करने की क्षमता भविष्य में निकट-दीवार टर्बुलेंस (near-wall turbulence) को हल करने के लाभ प्रदान कर सकती है, हालांकि इसके लिए आगे के अनुसंधान की आवश्यकता है।

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