🔬 condensed matter

Continuum granular flow model with restitution-derived viscoelastic damping

यह शोध पत्र दानेदार प्रवाह (granular flows) के लिए एक एकीकृत विस्कोइलास्टिक-विस्कोप्लास्टिक निरंतरता ढांचे (unified viscoelastic-visoplastic continuum framework) को प्रस्तुत करता है जो दर-निर्भर व्यवहार को मॉडल करने के लिए प्रत्यास्थता गुणांक (coefficient of restitution) और निरंतरता श्यानता (continuum viscosity) के बीच एक सीधा संबंध स्थापित करता है, जिसे शास्त्रीय $\mu(I)$ प्लास्टिक प्रवाह नियम को संरक्षित करते हुए, जटिल क्षणिक प्रक्रियाओं के अनुकरण के लिए मैटेरियल पॉइंट मेथड (material point method) के माध्यम से कार्यान्वित किया गया है।

मूल लेखक: Bodhinanda Chandra, Sachith Dunatunga, Ken Kamrin

प्रकाशित 2026-04-24

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

CC BY 4.0

मूल लेखक: Bodhinanda Chandra, Sachith Dunatunga, Ken Kamrin

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

रेत का एक ढेर, बॉल बेयरिंग की एक बाल्टी, या बर्फ का एक टीला कल्पना कीजिए। एक भौतिक विज्ञानी (physicist) के लिए, ये केवल चीजों के ढेर नहीं हैं; ये जटिल सामग्रियां हैं जो इस बात पर निर्भर करती हैं कि आप उन्हें कितनी तेज़ी से हिलाते हैं या उन्हें कितनी ज़ोर से मारते हैं—वे एक ठोस चट्टान की तरह व्यवहार कर सकती हैं, पानी की तरह बह सकती हैं, या गैस की तरह तैर भी सकती हैं।

यह शोध पत्र कंप्यूटर सिमुलेशन के लिए एक सार्वभौमिक "नियम पुस्तिका" (rulebook) बनाने के बारे में है जो सटीक रूप से भविष्यवाणी कर सके कि ये दानेदार सामग्रियां (granular materials) इन सभी विभिन्न अवस्थाओं में कैसे व्यवहार करती हैं। लेखकों, बोधिनांद चंद्र, सचित दुनातुंगा और केन काम्रिन ने एक नया गणितीय मॉडल बनाया है जो व्यक्तिगत दानों की सूक्ष्म दुनिया और हिमस्खलन (avalanches) व भूस्खलन (landslides) की बड़ी दुनिया के बीच के अंतर को पाटता है।

यहाँ उनके इस क्रांतिकारी आविष्कार का विवरण दिया गया है, जिसे रोजमर्रा के उदाहरणों के माध्यम से समझाया गया है।

1. समस्या: "बाउंस" बनाम "स्क्विशी" की दुविधा

अतीत में, कंप्यूटर मॉडलों को एक साथ दो विशिष्ट चीजों को संभालने में कठिनाई होती थी:

"बाउंस" (प्रतिघात/Restitution): जब दो दाने आपस में टकराते हैं, तो वे सुपरबॉल की तरह पूरी तरह से नहीं उछलते। वे कुछ ऊर्जा खो देते हैं। इसे "प्रतिघात गुणांक" (coefficient of restitution) कहा जाता है। यदि वे बहुत अधिक ऊर्जा खो देते हैं, तो वे चिपक जाते हैं; यदि वे कम ऊर्जा खोते हैं, तो वे उछलते हैं।
"प्रवाह" (प्लास्टिसिटी/Plasticity): जब आप रेत के ढेर को पर्याप्त ज़ोर से धकेलते हैं, तो यह शहद की तरह बहने लगता है। यह प्रवाह इस बात पर निर्भर करता है कि आप इसे कितनी तेज़ी से धकेल रहे हैं।

पुराने मॉडल एक ही बर्तन में दो अलग-अलग भोजन पकाने की कोशिश करने वाले शेफ की तरह थे। वे या तो "बाउंस" को अच्छी तरह से संभालते थे (गैस जैसी अवस्थाओं के लिए अच्छा) या "प्रवाह" को अच्छी तरह से संभालते थे (ठोस जैसी अवस्थाओं के लिए अच्छा), लेकिन दोनों को एक साथ बिना सिमुलेशन को अव्यवस्थित या अवास्तविक बनाए संभालना दुर्लभ था।

2. समाधान: एक "स्मार्ट शॉक एब्जॉर्बर"

लेखकों ने एक नया मॉडल बनाया है जो कार के स्मार्ट शॉक एब्जॉर्बर की तरह काम करता है।

उपमा: कल्पना कीजिए कि आप एक ऊबड़-खाबड़ सड़क पर गाड़ी चला रहे हैं।
- आपकी कार के स्प्रिंग्स दानों के वापस धकेलने की लोच (elasticity) का प्रतिनिधित्व करते हैं।
- शॉक एब्जॉर्बर (डैम्पर्स) श्यानता (viscosity) (घर्षण और ऊर्जा हानि) का प्रतिनिधित्व करते हैं।
- इंजन प्लास्टिक प्रवाह (plastic flow) (सामग्री का स्थायी रूप से विकृत होना) का प्रतिनिधित्व करता है।

यहाँ मुख्य नवाचार यह है कि उन्होंने यह पता लगाया कि शॉक एब्जॉर्बर को व्यक्तिगत दानों के उछाल (bounciness) के आधार पर कैसे ट्यून किया जाए। उन्होंने एक सीधा संबंध निकाला है: यदि आप जानते हैं कि एक एकल दाना कितना उछलता है (प्रतिघात गुणांक), तो आप गणितीय रूप से गणना कर सकते हैं कि पूरे ढेर को कितनी "विस्कस डैम्पिंग" (viscous damping) की आवश्यकता है।

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

3. "नो-टेंशन" नियम: अदृश्य धागा

दानेदार सामग्रियों की एक अजीब विशेषता होती है: उन्हें दबाया (squish) जा सकता है, लेकिन खींचा नहीं जा सकता।

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

4. तनाव के लिए "ट्रैफिक लाइट" प्रणाली

इन सिमुलेशन में सबसे बड़ी चुनौतियों में से एक यह सुनिश्चित करना है कि "बाउंस" (विस्कोइलास्टिसिटी) "प्रवाह" (प्लास्टिसिटी) में बाधा न डाले।

रूपक: एक व्यस्त चौराहे की कल्पना करें।
- इलास्टिक वेव्स (Elastic waves) चौराहे से गुजरने वाली कारों की तरह हैं।
- प्लास्टिक फ्लो (Plastic flow) एक निर्माण क्षेत्र (construction zone) की तरह है जहाँ सड़क स्थायी रूप से बदल रही है।
- लेखकों का मॉडल एक ट्रैफिक लाइट लगाता है। यह कहता है: "शॉक एब्जॉर्बर (डैम्पिंग) केवल कारों (तरंगों) को धीमा करते हैं। वे निर्माण दल (प्लास्टिक फ्लो) को धीमा नहीं करते हैं।"
यह क्यों मायने रखता है? यदि आप प्रवाह को बहुत अधिक डैम्प (धीमा) कर देते हैं, तो रेत बहने के बजाय रुक जाएगी। यदि आप तरंगों को डैम्प नहीं करते हैं, तो रेत सिमुलेशन में बेतहाशा कंपन करेगी और फट जाएगी। यह मॉडल दोनों प्रक्रियाओं को अलग रखते हुए भी उन्हें मिलकर काम करने में सक्षम बनाता है।

5. उन्होंने क्या सिद्ध किया? (प्रयोग)

उन्होंने अपने "स्मार्ट शॉक एब्जॉर्बर" मॉडल का परीक्षण पांच अलग-अलग परिदृश्यों के साथ किया, जैसे कि एक ड्राइवर अलग-अलग इलाकों में नई कार का परीक्षण करता है:

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

व्यापक परिदृश्य

यह शोध पत्र एक बड़ा कदम है क्योंकि यह टकराव (दानों का आपस में टकराना) के भौतिकी को प्रवाह (रेत का ढलान से बहना) के भौतिकी के साथ एकीकृत करता है।

इसे एक वीडियो गेम इंजन को अपग्रेड करने जैसा समझें। पहले, आपके पास एक ऐसा गेम हो सकता था जहाँ रेत वास्तविक रूप से बहती थी, या एक ऐसा गेम जहाँ रेत वास्तविक रूप से उछलती थी, लेकिन दोनों नहीं। अब, इस "रेस्टीट्यूशन-व्युत्पन्न विस्कोइलास्टिक डैम्पिंग" (restitution-derived viscoelastic damping) की मदद से, कंप्यूटर एक रेत के महल पर लहर के टकराने, पहाड़ से नीचे आते हिमस्खलन, या वाइब्रेटिंग बॉक्स में नाचते दानों को, एक सुसंगत सेट नियमों के साथ सिम्युलेट कर सकता है जो प्रत्येक व्यक्तिगत दाने के भौतिकी का सम्मान करता है।

यह एक एकल दाने के उछाल की सूक्ष्म दुनिया और लैंडस्लाइड की स्थूल दुनिया के बीच एक सेतु है, जो एक सरल, सुंदर गणितीय सूत्र द्वारा जुड़ी हुई है।

यहाँ चंद्रा, दुनाटुंगा और काम्रिन के शोध पत्र "Continuum granular flow model with restitution–derived viscoelastic damping" का विस्तृत तकनीकी सारांश दिया गया है।

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

शुष्क दानेदार प्रवाह (Dry granular flows) जटिल दर-निर्भर व्यवहार प्रदर्शित करते हैं जो दो अलग-अलग सूक्ष्म तंत्रों से उत्पन्न होते हैं:

माइक्रो-इन्ertia (Micro-inertia): यह सघन, प्लास्टिक शियर प्रवाह में प्रभावी होता है, जिससे $\mu(I)$ रियोलॉजी द्वारा वर्णित दर-निर्भर घर्षण उत्पन्न होता है।
पुनः स्थापन (Restitution - Collisional Dissipation): यह लोचदार कण टकरावों के दौरान ऊर्जा हानि को नियंत्रित करता है। यह बल्क इलास्टिक प्रतिक्रियाओं (जैसे कि स्ट्रेस वेव क्षीणन और सॉलिड-लाइक अवस्थाओं में दोलन डैम्पिंग) और गैस-लाइक अवस्थाओं में ऊर्जा हानि को प्रभावित करता है।

चुनौती: जबकि माइक्रो-इन्ertia को निरंतर मॉडल (जैसे कि $\mu(I)$ मॉडल के माध्यम से) में अच्छी तरह से शामिल किया गया है, इलास्टिक रेस्टिट्यूशन (Elastic restitution) को निरंतर यांत्रिकी (Continuum mechanics) में शामिल करने के लिए कोई व्यवस्थित ढांचा नहीं है, बिना प्लास्टिक प्रवाह व्यवहार को कृत्रिम रूप से बदले। मौजूदा दृष्टिकोण अक्सर कोलिशनल डिसिपेशन को सीधे प्लास्टिक विरूपण से जोड़ते हैं या एड-हॉक डैम्पिंग मापदंडों (जैसे कि रेली डैम्पिंग) पर निर्भर करते हैं जिनमें सूक्ष्म गुणांक ऑफ रेस्टिट्यूशन ( $e$ ) के साथ सीधा भौतिक संबंध नहीं होता है। लक्ष्य एक एकीकृत विस्कोइलास्टिक-विस्कोप्लास्टिक ढांचे को विकसित करना है जो सूक्ष्म गुणांक ऑफ रेस्टिट्यूशन को मैक्रोस्कोपिक कंटिनम विस्कोसिटी से जोड़ता है, जिससे यह सुनिश्चित हो सके कि विस्कस डैम्पिंग वेव प्रोपेगेशन और टकरावों को प्रभावित करती है, बिना प्लास्टिक फ्लो रूल को संशोधित किए।

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

A. सैद्धांतिक सूत्रीकरण (Theoretical Formulation)

लेखक निम्नलिखित सिद्धांतों पर आधारित एक एकीकृत विस्कोइलास्टिक-विस्कोप्लास्टिक कंटिनम फ्रेमवर्क प्रस्तावित करते हैं:

काइनेमैटिक डीकंपोजिशन (Kinematic Decomposition): वेग प्रवणता ( $L$ ) को विस्कोइलास्टिक ( $L_{ve}$ ) और प्लास्टिक ( $L_p$ ) घटकों में योगात्मक रूप से विभाजित किया जाता है।
स्ट्रेस डीकंपोजिशन (Stress Decomposition): केल्विन-वोइग्ट व्यवस्था (Kelvin-Voigt arrangement) का अनुसरण करते हुए, कुल स्ट्रेस को इलास्टिक ( $\sigma_e$ $σ_{e}$ ) और विस्कस ( $\sigma_v$ $σ_{v}$ ) भागों में विभाजित किया जाता है, जिन्हें फिर सीरीज व्यवस्था में प्लास्टिक स्ट्रेस ( $\sigma_p$ $σ_{p}$ ) के बराबर रखा जाता है।
- महत्वपूर्ण अंतर: विस्कस डैम्पिंग केवल इलास्टिक घटक ( $\sigma_v$ का $D_{ve}$ पर निर्भरता) पर लागू होती है। यह सुनिश्चित करता है कि विस्कस डिसिपेशन तरंगों और टकरावों को डैम्प करता है, लेकिन प्लास्टिक प्रवाह में कृत्रिम रूप से डैम्पिंग नहीं जोड़ता है, जो पूरी तरह से $\mu(I)$ रियोलॉजी द्वारा शासित रहता है।
नो-टेंशन कंडीशन (No-Tension Condition): मॉडल एक पृथक्करण नियम (Separation rule) को शामिल करता है जहाँ यदि बल्क डेंसिटी क्रिटिकल वैल्यू ( $\rho_c$ ) से नीचे गिर जाती है, तो सामग्री स्ट्रेस-फ्री ( $\sigma=0$ ) हो जाती है, जो ग्रेन सेपरेशन (गैस-लाइक अवस्था) का अनुकरण करती है।

B. रेस्टिट्यूशन-विस्कोसिटी लिंक का व्युत्पन्न (Derivation of Restitution-Viscosity Link)

एक केंद्रीय योगदान गुणांक ऑफ रेस्टिट्यूशन ( $e$ ) और कंटिनम बल्क विस्कोसिटी ( $\vartheta$ ) के बीच एक स्पष्ट संबंध का व्युत्पन्न है।

दृष्टिकोण: लेखक एक गोलाकार ग्रेनुलर असेंबली के वॉल्यूमेट्रिक कम्प्रेशन इम्पल्स के प्रति प्रतिक्रिया का विश्लेषण करते हैं। वे एक कंटिनम-स्तर के रेस्टिट्यूशन माप ( $E$ ) को वॉल्यूमेट्रिक स्ट्रेन रेट के अनुपात के रूप में परिभाषित करते हैं (टकराव से पहले और बाद में)।
कॉन्टैक्ट टाइम की परिभाषा: श्वागर और पोशेल (Schwager and Pöschel) का अनुसरण करते हुए, कॉन्टैक्ट टाइम ( $t_c$ ) को उस क्षण के रूप में परिभाषित किया गया है जब संपर्क बल (स्ट्रेस) शून्य हो जाता है, न कि जब विस्थापन शून्य होता है। यह ग्रेनुलर कॉन्टैक्ट की एकतरफा प्रकृति (कोई टेंशन नहीं) को ध्यान में रखता है।
परिणाम: गोलाकार असेंबली के लिए डैम्प्ड वेव इक्वेशन को हल करके, वे आयामहीन बल्क विस्कोसिटी ( $\tilde{\vartheta}$ $\tilde{ϑ}$ ) और कंटिनम रेस्टिट्यूशन ( $E$ $E$ ) के बीच एक क्लोज्ड-फॉर्म संबंध प्राप्त करते हैं।
- अनुमानित संबंध: $\tilde{\vartheta} \approx 0.237 \frac{|\ln E|}{\pi^2}$ ।
- यह प्रयोगिक गुणांक ऑफ रेस्टिट्यूशन ( $e$ ) को सीधे कंटिनम बल्क विस्कोसिटी ( $\vartheta$ ) में बदलने की अनुमति देता है।
शियर विस्कोसिटी (Shear Viscosity): शियर विस्कोसिटी ( $\eta$ ) को रेली-टाइप स्टिफनेस-प्रोपोर्शनल डैम्पिंग के माध्यम से बल्क विस्कोसिटी ( $\vartheta$ ) से जोड़ा गया है: $\eta = (G/K)\vartheta$ ।

C. संख्यात्मक कार्यान्वयन (Numerical Implementation)

सॉल्वर (Solver): मॉडल को मटेरियल पॉइंट मेथड (MPM) के भीतर लागू किया गया है, जो एक हाइब्रिड यूलेरियन-लैग्रेंजियन फ्रेमवर्क है और बड़े विरूपण, सामग्री पृथक्करण और पुनर्गठन के लिए आदर्श है।
स्ट्रेस अपडेट एल्गोरिदम: एक प्रेडिक्टर-करेक्टर स्कीम विकसित की गई है।
1. विस्कस टर्म्स सहित ट्रायल इलास्टिक स्ट्रेस की गणना करें।
2. यील्ड कंडीशन ( $\tau \leq \mu(I)p$ ) की जाँच करें।
3. यदि यील्डिंग होती है, तो प्लास्टिक स्ट्रेन रेट और स्ट्रेस को अपडेट करने के लिए एक रिटर्न-मैपिंग एल्गोरिदम का उपयोग करें, यह सुनिश्चित करते हुए कि प्लास्टिक फ्लो केवल $\mu(I)$ रियोलॉजी द्वारा संचालित है जबकि विस्कस टर्म्स इलास्टिक रिस्पॉन्स को डैम्प करते हैं।

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

एकीकृत ढांचा (Unified Framework): एक एकल कंटिनम विवरण जो सॉलिड-लाइक (इलास्टिक/प्लास्टिक), लिक्विड-लाइक (प्लास्टिक फ्लो), और गैस-लाइक (सेपरेटेड) अवस्थाओं को सहजता से संभालता है।
भौतिक लिंक (Physical Link): सूक्ष्म कोलिशनल रेस्टिट्यूशन और मैक्रोस्कोपिक कंटिनम विस्कोसिटी के बीच एक कठोर संबंध का व्युत्पन्न, जिससे एड-हॉक डैम्पिंग कैलिब्रेशन की आवश्यकता समाप्त हो जाती है।
डीकपल्ड डिसिपेशन (Decoupled Dissipation): एक कंटिन्यू स्ट्रैटेजी जहाँ विस्कस डैम्पिंग इलास्टिक वेव्स और टकरावों को कम करती है, बिना माइक्रो-इन्ertia द्वारा नियंत्रित दर-निर्भर प्लास्टिक फ्लो में हस्तक्षेप किए।
संख्यात्मक स्थिरता (Numerical Stability): भौतिक विस्कोइलास्टिक डैम्पिंग का समावेश उन संख्यात्मक अस्थिरताओं (स्प्यूरियस ऑसिलेशन) को कम करता है जो मानक एक्सप्लिसिट स्कीम्स का उपयोग करने वाले कोलिशनल ग्रेनुलर फ्लो में देखी जाती हैं।

4. परिणाम और सत्यापन (Results and Validation)

मॉडल को पाँच संख्यात्मक उदाहरणों के माध्यम से सत्यापित किया गया:

गोलाकार संपीड़न (Spherical Compaction): बल्क विस्कोसिटी और कंटिनम रेस्टिट्यूशन के बीच व्युत्पन्न विश्लेषणात्मक संबंध को सत्यापित किया गया। संख्यात्मक परिणाम वेग अनुपातों के एक्सपोनेंशियल डिके (exponential decay) के सैद्धांतिक अनुमान से मेल खाते हैं।
इनक्लाइंड प्लेन फ्लो (Inclined Plane Flow): इसने प्रदर्शित किया कि गुणांक ऑफ रेस्टिट्यूशन (और इस प्रकार विस्कोसिटी) को बदलने से स्थिर बेनॉल्ड वेलोसिटी प्रोफाइल (Bagnold velocity profile) में कोई बदलाव नहीं आता है। यह पुष्टि करता है कि मॉडल सही ढंग से प्लास्टिक फ्लो नियम से विस्कस डैम्पिंग को अलग करता है, जो $\mu(I)$ सिद्धांत के अनुरूप है।
फ्लैट-बॉटम साइलो फ्लो (Flat-Bottom Silo Flow): इसने दिखाया कि जबकि डिस्चार्ज रेट (जो प्लास्टिक फ्लो द्वारा नियंत्रित है) $e$ के प्रति असंवेदनशील है, रिकंसोलिडेशन व्यवहार अत्यधिक संवेदनशील है। कम रेस्टिट्यूशन (उच्च डैम्पिंग) ने अवास्तविक उछाल (unphysical bouncing) को कम किया और कम रनआउट दूरी और तीव्र डायनेमिक रिपोज एंगल का परिणाम दिया, जो DEM अवलोकनों से मेल खाता है।
ग्रेनुलर बेड इम्पैक्ट (Granular Bed Impact): एक इम्पैक्टर द्वारा ग्रेनुलर बेड से टकराने का अनुकरण किया गया। मॉडल ने उच्च-आवृत्ति वाले स्ट्रेस वेव्स और अनडैम्प्ड इलास्टिक वाइब्रेशन से जुड़े बीट फ्रीक्वेंसी को सफलतापूर्वक कम किया, जबकि सही प्लास्टिक विरूपण और इजेक्टा पैटर्न को बनाए रखा।
वाइब्रेटेड ग्रेनुलर पैटर्न्स (Vibrated Granular Patterns): वर्टिकल वाइब्रेटेड ग्रेनुलर लेयर्स में स्क्वायर और डायमंड पैटर्न के निर्माण को पुनरुत्पादित किया गया। महत्वपूर्ण रूप से, ये पैटर्न केवल तभी उभरे जब विस्कोइलास्टिक डैम्पिंग (परिमित $e < 1$ ) को शामिल किया गया, जो प्रायोगिक अवलोकनों से मेल खाता है। शुद्ध इलास्टिक सिमुलेशन ( $e=1$ ) स्थिर पैटर्न बनाने में विफल रहे।

5. महत्व (Significance)

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

सघन माध्यमों में वेव प्रोपेगेशन और क्षीणन (attenuation)।
क्षणिक प्रभाव (Transient impacts) और ऊर्जा का क्षय।
सॉलिड, लिक्विड और गैस अवस्थाओं के बीच फेज ट्रांजिशन।
वाइब्रेटेड सिस्टम में पैटर्न फॉर्मेशन।

यह ढांचा विशेष रूप से उन अनुप्रयोगों के लिए महत्वपूर्ण है जिनमें प्लास्टिक फ्लो भविष्यवाणियों की सटीकता से समझौता किए बिना सटीक ऊर्जा क्षय मॉडलिंग की आवश्यकता होती है, जैसे कि जियोटेक्निकल इंजीनियरिंग, प्लेनेटरी साइंस, और इंडस्ट्रियल पाउडर प्रोसेसिंग। MPM में इसका कार्यान्वयन इसे अत्यधिक विरूपण और सामग्री विखंडन वाली समस्याओं के लिए और अधिक उपयोगी बनाता है।