An adaptive perfectly matched layer finite element method for acoustic-elastic interaction in periodic structures

तकनीकी सारांश: आवधिक संरचनाओं में ध्वनिक-प्रत्यास्थ अंतःक्रिया के लिए एक अनुकूली पूर्णतः सुमेलित परत (PML) परिमित तत्व विधि (FEM)

1. समस्या विवरण

यह शोध पत्र एक अनबाउंड (unbounded) आवधिक सतह (एक विवर्तन ग्रेटिंग) द्वारा अभिलक्षित एक प्रत्यास्थ पिंड (elastic body) द्वारा एक समय-हार्मोनिक ध्वनिक समतलीय तरंग (acoustic plane wave) के प्रकीर्णन को संबोधित करता है। भौतिक प्रणाली में दो डोमेन शामिल हैं:

• द्रव डोमेन ($\Omega_u$): आवधिक सतह के ऊपर एक सजातीय, संपीड्य, अक्रिय तरल माध्यम।
• ठोस डोमेन ($\Omega_d$): सतह के नीचे एक समदैशिक (isotropic), रैखिक प्रत्यास्थ ठोस पदार्थ।

अंतःक्रिया को द्रव में हेल्महोल्ट्ज़ समीकरण (Helmholtz equation) और ठोस में नेवियर समीकरण (Navier equation) द्वारा नियंत्रित किया जाता है। इंटरफ़ेस ($\Gamma$) पर युग्मन (coupling) के लिए सामान्य वेग की निरंतरता (kinematic condition) और ट्रैक्शन की निरंतरता (dynamic condition) की आवश्यकता होती है। चूंकि संरचना आवधिक है, इसलिए समस्या का विश्लेषण एक एकल आवधिक सेल के भीतर किया जा सकता है, लेकिन ऊर्ध्वाधर दिशा में डोमेन अनबाउंड रहते हैं, जिसके लिए अनंत पर विकिरण स्थितियों (radiation conditions) को संभालने के लिए परिष्कृत ट्रंकेशन विधियों की आवश्यकता होती है।

2. कार्यप्रणाली

लेखक एक अनुकूली पूर्णतः सुमेलित परत (PML) परिमित तत्व विधि (FEM) का प्रस्ताव करते हैं। कार्यप्रणाली को तीन मुख्य घटकों में संरचित किया गया है:

• PML के माध्यम से डोमेन ट्रंकेशन: अनबाउंड समस्या को बाउंडेड (bounded) समस्या में बदलने के लिए, लेखक जटिल निर्देशांक विस्तार (complex coordinate stretching) का उपयोग करके एक PML तकनीक लागू करते हैं। एक प्रमुख तकनीकी उपलब्धि ध्वनिक और प्रत्यास्थ दोनों तरंगों के लिए समतुल्य पारदर्शी सीमा स्थितियों (transparent boundary conditions - TBCs) का एक साथ निर्माण करना है। यह बाहरी PML सीमाओं पर सरल होमोजेनियस डिरिचलेट सीमा स्थितियों का उपयोग करके ट्रंकेटेड समस्या को हल करने की अनुमति देता है।
• गणितीय विश्लेषण:
  • सुव्यवस्थितता (Well-posedness): लेखक ट्रंकेटेड PML वेरिएशनल फॉर्मूलेशन की अद्वितीय समाधान क्षमता (unique solvability) स्थापित करते हैं।
  • अभिसरण (Convergence): वे कठोरता से सिद्ध करते हैं कि PML सन्निकटन त्रुटि (approximation error), PML मापदंडों (परत की मोटाई और अवशोषण शक्ति) के संबंध में घातांकीय (exponentially) रूप से घटती है।
  • ए पोस्टीरियो त्रुटि अनुमान (A Posteriori Error Estimation): एक रेसिडुअल-प्रकार का त्रुटि अनुमानक व्युत्पन्न किया गया है। महत्वपूर्ण रूप से, यह अनुमानक दो अलग-अलग त्रुटि स्रोतों को स्पष्ट रूप से समाहित करता है: परिमित तत्व विविक्तीकरण त्रुटि (finite element discretization error) और PML ट्रंकेशन त्रुटि।
• अनुकूली एल्गोरिदम: त्रुटि अनुमानक का उपयोग करते हुए, लेखक एक अनुकूली परिशोधन रणनीति (adaptive refinement strategy) विकसित करते हैं। एल्गोरिदम एक "मार्किंग" रणनीति का उपयोग करता है ताकि उन क्षेत्रों में स्थानीय रूप से मेश (mesh) को परिष्कृत किया जा सके जहाँ समाधान विलक्षणताएं (singularities) प्रदर्शित करता है (जैसे कि रीएंट्रेंट कॉर्नर्स या गैर-चिकनी सतह प्रोफाइल), जिससे कंप्यूटेशनल संसाधनों को अनुकूलित किया जा सके।

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

1. बहुभौतिकी (Multiphysics) PML ढांचा: पिछले कार्यों के विपरीत जो एकल-भौतिकी प्रकीर्णन पर ध्यान केंद्रित करते हैं, यह शोध पत्र आवधिक संरचनाओं में युग्मित ध्वनिक-प्रत्यास्थ अंतःक्रिया के लिए एक कठोर सैद्धांतिक ढांचा प्रदान करता है।
2. कठोर त्रुटि विश्लेषण: एक रेसिडुअल-आधारित ए पोस्टीरियो त्रुटि अनुमानक का व्युत्पन्न, जो PML ट्रंकेशन त्रुटि को भी शामिल करता है, एक महत्वपूर्ण सैद्धांतिक प्रगति है।
3. विलक्षणताओं के लिए अनुकूली रणनीति: ध्वनिक-प्रत्यास्थ युग्मन की जटिलताओं के लिए विशेष रूप से तैयार किया गया एक अनुकूली FEM एल्गोरिदम, जो ज्यामितीय विलक्षणताओं को कुशलतापूर्वक संभालने की अनुमति देता है।

4. परिणाम और संख्यात्मक सत्यापन

प्रस्तावित विधि की प्रभावशीलता को चार अलग-अलग संख्यात्मक परिदृश्यों के माध्यम से सत्यापित किया गया:

• सटीकता सत्यापन: एक विश्लेषणात्मक समाधान वाले चिकने-इंटरफेस मामले में, PML-FEM ने सटीक समाधान के साथ उत्कृष्ट सहमति दिखाई और इष्टतम अभिसरण दर प्रदर्शित की।
• विलक्षणता हैंडलिंग: तीखे कोनों (ट्रैपेज़ॉइडल और एंगुलर ग्रेटिंग्स) वाले मामलों में, अनुकूली एल्गोरिदम ने विलक्षण बिंदुओं के आसपास मेश परिशोधन को सफलतापूर्वक केंद्रित किया। इसने समान परिशोधन (uniform refinement) की तुलना में बहुत कम डिग्री ऑफ फ्रीडम (DOF) के साथ काफी उच्च सटीकता प्राप्त की।
• उच्च-आवृत्ति मजबूती: विधि उच्च-वेव नंबर ($\kappa = 20$) समस्याओं के लिए स्थिर और प्रभावी रही, जो व्यावहारिक इंजीनियरिंग अनुप्रयोगों के लिए इसकी उपयोगिता को प्रदर्शित करती है।
• दक्षता: डिरिचले-टू-नेमैन (DtN) विधि के साथ तुलना से पता चला कि जबकि प्रति DOF में DtN थोड़ा अधिक सटीक हो सकता है, PML विधि को लागू करना बहुत सरल है क्योंकि यह गैर-स्थानीय सीमा स्थितियों (non-local boundary conditions) से बचता है।

5. महत्व

यह शोध अंडरवॉटर ध्वनिकी (underwater acoustics) और अल्ट्रासोनिक गैर-विनाशकारी मूल्यांकन (ultrasonic non-destructive evaluation) जैसे जटिल, आवधिक माध्यमों में तरंग प्रसार से जुड़े क्षेत्रों के लिए अत्यधिक महत्वपूर्ण है। ध्वनिक-प्रत्यास्थ अंतःक्रियाओं को सिम्युलेट करने के लिए एक गणितीय रूप से कठोर और गणनात्मक रूप से कुशल तरीका प्रदान करके, यह शोध पत्र ऑप्टिकल घटकों, सेंसरों और ध्वनिक मेटामटेरियल्स के डिजाइन के लिए एक शक्तिशाली उपकरण प्रदान करता है। ज्यामितीय विलक्षणताओं को अनुकूली परिशोधन के माध्यम से संभालने की क्षमता इसे वास्तविक दुनिया की इंजीनियरिंग संरचनाओं के लिए विशेष रूप से उपयोगी बनाती है जो शायद ही कभी पूरी तरह से चिकनी होती हैं।