🔢 mathematics

Mass-Lumped Virtual Element Method with Strong Stability-Preserving Runge-Kutta Time Stepping for Two-Dimensional Parabolic Problems

यह शोध पत्र सामान्य बहुभुज मेषों (polygonal meshes) पर द्वि-आयामी परवलयिक समस्याओं (parabolic problems) के लिए एक मास-लम्पड वर्चुअल एलीमेंट मेथड (mass-lumped Virtual Element Method) को स्पष्ट स्ट्रॉन्ग स्टेबिलिटी-प्रिजर्विंग रनगे-कुट्टा (explicit Strong Stability-Preserving Runge-Kutta) टाइम स्टेपिंग के साथ प्रस्तुत करता है, जो एक शास्त्रीय $\mathcal{O}(h^2)$ CFL स्थिति के तहत सैद्धांतिक स्थिरता स्थापित करता है और संख्यात्मक प्रयोगों के माध्यम से इष्टतम अभिसरण दरों (optimal convergence rates) और मेष विरूपण (mesh distortion) के विरुद्ध सुदृढ़ता का प्रदर्शन करता है।

मूल लेखक: Paulo Akira F. Enabe, Rodrigo Provasi

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

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

CC BY 4.0

मूल लेखक: Paulo Akira F. Enabe, Rodrigo Provasi

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

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

आमतौर पर, ये पहेली के टुकड़े सटीक वर्ग या त्रिकोण होते हैं। लेकिन वास्तविक दुनिया में, चीजें अव्यवस्थित होती हैं। चट्टानें, कोशिकाएं, या फटी हुई सतहें सटीक वर्गों में फिट नहीं बैठतीं। यहीं पर वर्चुअल एलीमेंट मेथड (VEM) काम आता है। यह एक सुपर-लचीले पहेली समाधानकर्ता की तरह है जो किसी भी आकार के टुकड़े—त्रिकोण, पंचकोण, अजीबोगरीब आकार, या विकृत वर्ग—को बिना भ्रमित हुए संभाल सकता है।

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

मुख्य विचार: "मास लंप्ड" (Mass Lumping) और "SSP-RK"

यह शोध पत्र इस प्रक्रिया को तेज़, स्थिर और सटीक बनाने के लिए एक चतुर शॉर्टकट पेश करता है, भले ही आकार अनियमित हों। वे दो मुख्य तरकीबों को मिलाते हैं:

1. "लम्प्ड मास" की तरकीब (पहेली को सरल बनाना)

कल्पना कीजिए कि "मास मैट्रिक्स" एक बहीखाता (ledger) है जो ट्रैक करता है कि प्रत्येक पहेली के टुकड़े में कितना "सामग्री" (स्याही या ऊष्मा) है और वह अपने पड़ोसियों से कैसे जुड़ता है।

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

2. "SSP-RK" टाइम स्टेपिंग (स्मार्ट रनर)

अब जबकि गणित सरल हो गया है, तो आपको बिना क्रैश हुए समय में आगे बढ़ने की रणनीति की आवश्यकता है।

समस्या: यदि आप बहुत बड़े कदम लेते हैं, तो सिमुलेशन फट जाएगा (स्याही गायब हो जाएगी या नकारात्मक संख्याओं में बदल जाएगी)। यदि आप बहुत छोटे कदम लेते हैं, तो इसमें बहुत समय लगेगा।
समाधान: वे स्ट्रॉन्ग स्टेबिलिटी-प्रिजर्विंग रंज (Strong Stability-Preserving Runge-Kutta - SSP-RK) विधियों का उपयोग करते हैं।
- उपमा: कल्पना कीजिए कि एक धावक (सिमुलेशन) एक खड़ी ढलान (डिफ्यूजन प्रक्रिया) से नीचे दौड़ने की कोशिश कर रहा है।
  - फॉरवर्ड यूलर (शुरुआती): शुरुआती बहुत छोटे, सुरक्षित कदम लेता है। वे गिरेंगे नहीं, लेकिन वे धीमे हैं।
  - SSP-RK (प्रो): प्रो धावक जानता है कि बड़े, तेज़ कदम (उच्च-क्रम सटीकता) कैसे लिए जाते हैं, लेकिन उसके पास एक विशेष नियम है: "कभी भी ऐसा कदम न लें जिससे आपका संतुलन बिगड़ जाए।"
- यह शोध पत्र सिद्ध करता है कि यदि शुरुआती (फॉरवर्ड यूलर) एक निश्चित स्टेप साइज के तहत सुरक्षित है, तो प्रो (SSP-RK) बड़े कदम ले सकता है और साथ ही यह गारंटी दे सकता है कि वह उतना ही सुरक्षित रहेगा। यह सिमुलेशन को जोखिम के बिना तेज़ बनाता है।

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

लेखकों ने इसे तीन प्रकार के "अव्यवस्थित" मानचित्रों पर परखा:

विकृत वर्ग (Distorted Squares): ऐसे वर्ग जिन्हें दबाया या खींचा गया है।
सेरेंडिपिटी शेप्स (Serendipity Shapes): विशेष आकार जो कुशल होने के लिए डिज़ाइन किए गए हैं।
वोरोनोई सेल्स (Voronoi Cells): यादृच्छिक, जैविक आकार (जैसे झाग में बुलबुले या पत्ती में कोशिकाएं)।

परिणाम:

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

निष्कर्ष

यह शोध पत्र एक बंपी, ऑफ-रोड इलाके के लिए जीपीएस (GPS) बनाने जैसा है।

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

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

यहाँ शोध पत्र "Mass-Lumped Virtual Element Method with Strong Stability-Preserving Runge-Kutta Time Stepping for Two-Dimensional Parabolic Problems" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

यह शोध पत्र सामान्य बहुभुज मेष (general polygonal meshes) पर द्वि-आयामी समय-निर्भर परवलयिक समस्याओं (विशेष रूप से स्केलर डिफ्यूजन समीकरणों) के संख्यात्मक समाधान को संबोधित करता है। नियंत्रित समीकरण है:
$\frac{\partial u}{\partial t} - \Delta u = f \quad \text{in } \Omega \times (0, T]$
जिसमें होमोजेनियस डिरिचलेट बाउंड्री कंडीशंस और एक दी गई इनिशियल कंडीशन है।

मुख्य चुनौतियाँ:

मेष जनरैलिटी (Mesh Generality): मानक फाईनाइट एलीमेंट मेथड्स (FEM) गैर-उत्तल (non-convex) या अत्यधिक विकृत बहुभुज तत्वों के साथ संघर्ष करते हैं। वर्चुअल एलीमेंट मेथड (VEM) इन्हें स्वाभाविक रूप से संभालता है लेकिन आमतौर पर घने मास मैट्रिसेस (dense mass matrices) का परिणाम देता है।
एक्सप्लिसिट टाइम स्टेपिंग (Explicit Time Stepping): एक्सप्लिसिट स्कीम्स समानांतरकरण (parallelization) और ग्लोबल लीनियर सोल्व्स से बचने के लिए आकर्षक हैं, लेकिन उन्हें प्रत्येक स्टेप में मास मैट्रिक्स को इनवर्ट करने की आवश्यकता होती है। मानक VEM में, मास मैट्रिक्स विकर्ण (diagonal) नहीं होता है, जिससे इनवर्जन महंगा हो जाता है।
स्थिरता (Stability): डिफ्यूजन समस्याओं के लिए एक्सप्लिसिट स्कीम्स सख्त कॉोरेंट-फ्रेड्रिक्स-लेवी (CFL) स्थितियों ( $\Delta t \propto h^2$ ) के अधीन होती हैं। मास लंपिंग के साथ मनमाने बहुभुज मेषों पर स्थिरता और धनात्मकता (strong stability) सुनिश्चित करना गैर-तुच्छ कार्य है।

2. कार्यप्रणाली (Methodology)

प्रस्तावित ढांचा तीन मुख्य घटकों को जोड़ता है: स्थानिक विविक्तीकरण (spatial discretization) के लिए Virtual Element Method (VEM), विकर्ण बनाने के लिए Mass Lpping, और समय एकीकरण (time integration) के लिए Strong Stability-Preserving Runge-Kutta (SSP-RK)।

A. स्थानिक विविक्तीकरण (Spatial Discretization - VEM)

लेखक क्रम $k$ के कन्फॉर्मिंग VEM (विशेष रूप से प्रयोगों में $k=1$ ) का उपयोग करते हैं।

स्पेस (Spaces): स्थानीय स्पेस $V_{h,k}^E$ में किनारों पर बहुपद ट्रेस (polynomial traces) वाले फलन और $P_{k-2}(E)$ में लैप्लासियन शामिल हैं।
प्रोजेक्टर्स (Projectors): कंप्यूट करने योग्य द्विपद रूपों (bilinear forms) का निर्माण किया गया है:
- एनर्जी प्रोजेक्टर ( $\Pi^\nabla_k$ ): स्टिफनेस मैट्रिक्स (डिफ्यूजन टर्म) के लिए।
- $L^2$ प्रोजेक्टर ( $\Pi^0_k$ ): मास मैट्रिक्स और लोड टर्म्स के लिए।
स्थिरीकरण (Stabilization): प्रोजेक्टर के कर्नेल में बहुपद मोड को नियंत्रित करने के लिए स्टिफनेस फॉर्म में एक स्थिरीकरण पद जोड़ा जाता है।

B. मास लंपिंग रणनीति (Mass Lumping Strategy)

एक्सप्लिसिट टाइम स्टेपिंग को सक्षम करने के लिए, कंसिस्टेंट मास मैट्रिक्स $M_h$ को एक विकर्ण मैट्रिक्स $\hat{M}_h$ द्वारा प्रतिस्थापित किया जाता है।

रो-सम लंपिंग (Row-Sum Lumping): विकर्ण प्रविष्टियों की गणना कंसिस्टेंट मास मैट्रिक्स की पंक्तियों के योग के रूप में की जाती है।
लुप्त स्थिरीकरण (Vanishing Stabilization): एक प्रमुख सैद्धांतिक निष्कर्ष यह है कि रो-समिंग के तहत मास मैट्रिक्स में स्थिरीकरण पद स्वतः ही शून्य हो जाते हैं क्योंकि स्थिरीकरण प्रोजेक्टर के कर्नेल पर कार्य करता है, जबकि रो सम प्रोजेक्टर के इमेज (image) में मौजूद स्थिरांकों के एकीकरण के अनुरूप होते हैं।
फ्लोरिंग वेट्स (Floored Weights): $k \ge 2$ के लिए, सरल रो सम के परिणामस्वरूप शून्य या ऋणात्मक विकर्ण प्रविष्टियाँ प्राप्त हो सकती हैं। लेखक एक फ्लोरिंग मैकेनिज्म (भार को एक न्यूनतम धनात्मक मान $\delta |E|/N_{dof}$ पर सेट करना) पेश करते हैं ताकि तत्वों के किनारों की संख्या से स्वतंत्र, एकरूप धनात्मकता (uniform positivity) और $L^2$ -तुल्यता सुनिश्चित की जा सके।
दक्षता (Efficiency): लुम्प्ड वेट्स को एक छोटे स्थानीय बहुपद सिस्टम $G_E w_E = c_E$ को हल करके कंप्यूट किया जा सकता है जिसका खर्च $O(N_k^3)$ है, जहाँ $N_k \ll N_{dof}$ ।

C. समय एकीकरण (Time Integration - SSP-RK)

सेमी-डिस्क्रीट सिस्टम $\hat{M}_h \dot{u}_h + K_h u_h = f_h$ को एक्सप्लिसिट Strong Stability-Preserving Runge-Kutta (SSP-RK) स्कीम्स का उपयोग करके हल किया जाता है।

मैकेनिज्म: SSP-RK विधियाँ (विशेष रूप से 3rd-order SSP-RK3 और 4th-order SSP-RK(5,4)) फॉरवर्ड यूलर स्टेप्स के उत्तल संयोजन (convex combinations) हैं।
स्थिरता विरासत (Stability Inheritance): यदि फॉरवर्ड यूलर स्टेप एक विशिष्ट टाइम-स्टेप प्रतिबंध के तहत स्थिर (contractive) है, तो SSP-RK स्कीम इस स्थिरता गुण को विरासत में लेती है, जिससे मोनोटोनिसिटी और ऊर्जा क्षय (energy decay) बना रहता है।
CFL स्थिति: फॉरवर्ड यूलर स्टेप की स्थिरता $(\hat{M}_h)^{-1}K_h$ के स्पेक्ट्रल रेडियस से जुड़ी है।

3. प्रमुख योगदान (Key Contributions)

मास-लम्प्ड VEM फॉर्मूलेशन:
- एक विकर्ण मास मैट्रिक्स वाला Vented VEM का विकास जहाँ स्थिरीकरण पद रो-समिंग के दौरान स्वाभाविक रूप से लुप्त हो जाते हैं।
- एक "फ्लोरिंग वेट" तकनीक का परिचय जो सामान्य बहुभुजों पर मास मैट्रिक्स की एकरूप धनात्मकता सुनिश्चित करती है, जिससे एक सममित धनात्मक निश्चित (symmetric positive definite) डिस्क्रीट इनर प्रोडक्ट सुनिश्चित होता है।
मेष-रोबस्ट स्पेक्ट्रल बाउंड्स (Mesh-Robust Spectral Bounds):
- यह प्रमाण कि डिस्क्रीट ऑपरेटर का सबसे बड़ा आइजनवैल्यू $\lambda_{\max} \le C h^{-2}$ के रूप में स्केल करता है।
- महत्वपूर्ण रूप से, स्थिरांक $C$ केवल स्पेस डायमेंशन, बहुपद डिग्री और मेष रेगुलैरिटी (chunkiness) पर निर्भर करता है, और यह तत्व के किनारों/चेहरों की संख्या से स्वतंत्र है। यह जटिल मेषों पर मानक डिफ्यूजन-टाइप CFL स्थितियों के उपयोग को मान्य करता है।
एक्सप्लिसिट स्कीम्स के लिए स्थिरता विश्लेषण:
- $\hat{M}_h$ द्वारा प्रेरित डिस्क्रीट एनर्जी नॉर्म में फॉरवर्ड यूलर और SSP-RK स्कीम्स के लिए कठोर स्थिरता अनुमानों का व्युत्पन्न।
- टाइम-स्टेप प्रतिबंध $\Delta t \le C_{SSP} \Delta t_{FE}, \text{ जहाँ } \Delta t_{FE} \propto h^2$ की स्थापना।
कंप्यूटेशनल दक्षता:
- यह प्रदर्शन कि लुम्प्ड वेट्स को स्थानीय रूप से और कुशलतापूर्वक ( $O(N_k^3)$ प्रति तत्व) कंप्यूट किया जा सकता है, जो ग्लोबल सोल्व्स से बचता है।

4. संख्यात्मक परिणाम (Numerical Results)

लेखकों ने विकृत चतुर्भुज (distorted quadrilaterals), सेरेन्डिपिटी एलिमेंट्स (serendipity elements) और वोरोनोई मेष (Voronoi meshes) पर व्यापक संख्यात्मक प्रयोगों के माध्यम से सिद्धांत को सत्यापित किया।

अभिसरण दरें (Convergence Rates):
- $k=1$ के लिए, विधि इष्टतम अभिसरण प्राप्त करती है: $H^1$ -सेमी-नॉर्म में $O(h)$ और $L^2$ -नॉर्म में $O(h^2)$ ।
- मेष विरूपण और डायगोनल मास सन्निकटन के बावजूद सटीकता बनी रहती है।
- उच्च-क्रम वाली SSP-RK (RK3, RK4) स्कीम्स इन स्थानिक दरों को बनाए रखते हुए एरर कांस्टेंट को कम करती हैं।
स्थिरता और CFL:
- संख्यात्मक प्रयोग अनुमानित $\Delta t \propto h^2$ स्केलिंग की पुष्टि करते हैं।
- विधि अत्यधिक अनियमित वोरोनोई मेष पर भी स्थिर रहती है।
- हेट्रोजेनियस एनिसोट्रोपिक डिफ्यूजन परीक्षणों में (विच्छिन्न और रोटेटेड टेंसर के साथ), अधिकतम स्थिर टाइम स्टेप डिफ्यूजन टेंसर के सबसे बड़े आइजनवैल्यू के व्युत्क्रमानुपाती स्केल करता है ( $\Delta t_{max} \propto 1/\lambda_{\max}(K)$ ), जो मजबूती की पुष्टि करता है।
दक्षता तुलना (Explicit vs. Implicit):
- एक मानक कंसिस्टेंट-मास इमप्लिसिट यूलर स्कीम की तुलना में, एक्सप्लिसिट लम्प्ड VEM को काफी अधिक टाइम स्टेप्स की आवश्यकता होती है ( $\Delta t \propto h^2$ बनाम $\Delta t \propto h$ )।
- हालांकि, क्योंकि VEM में डोमिनेंट कॉस्ट मैट्रिक्स असेंबली (प्रोजेक्शन कंप्यूट करना) है न कि लीनियर सोल्व्स, दोनों परीक्षण किए गए परिदृश्यों में कुल वॉल-क्लॉक टाइम तुलनीय है।
- एक्सप्लिसिट विधि एक सरल कार्यान्वयन प्रदान करती है (कोई ग्लोबल लीनियर सॉल्वर नहीं) और स्वाभाविक रूप से समानांतर करने योग्य (parallelizable) है।

5. महत्व (Significance)

यह कार्य सामान्य बहुभुज मेषों पर परवलयिक समस्याओं के लिए एक प्रमाणित रूप से स्थिर, एक्सप्लिसिट टाइम-स्टेपिंग फ्रेमवर्क प्रदान करके वर्चुअल एलीमेंट मेथड के साहित्य में एक महत्वपूर्ण अंतर को पाटता है।

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

संक्षेप में, यह शोध पत्र स्थापित करता है कि Mass-Lumped VEM + SSP-RK, मनमाने बहुभुज ग्रिडों पर क्षणिक डिफ्यूजन समस्याओं के लिए इमप्लिसिट स्कीम्स के एक कठोर, सटीक और कंप्यूटेशनल रूप से प्रतिस्पर्धी विकल्प के रूप में कार्य करता है।