2025-11-18T02:58:13.538030

Reduced order modelling of nonaffine problems on parameterized NURBS multipatch geometries

Chasapi, Antolin, Buffa

This contribution explores the combined capabilities of reduced basis methods and IsoGeometric Analysis (IGA) in the context of parameterized partial differential equations. The introduction of IGA enables a unified simulation framework based on a single geometry representation for both design and analysis. The coupling of reduced basis methods with IGA has been motivated in particular by their combined capabilities for geometric design and solution of parameterized geometries. In most IGA applications, the geometry is modelled by multiple patches with different physical or geometrical parameters. In particular, we are interested in nonaffine problems characterized by a high-dimensional parameter space. We consider the Empirical Interpolation Method (EIM) to recover an affine parametric dependence and combine domain decomposition to reduce the dimensionality. We couple spline patches in a parameterized setting, where multiple evaluations are performed for a given set of geometrical parameters, and employ the Static Condensation Reduced Basis Element (SCRBE) method. At the common interface between adjacent patches a static condensation procedure is employed, whereas in the interior a reduced basis approximation enables an efficient offline/online decomposition. The full order model over which we setup the RB formulation is based on NURBS approximation, whereas the reduced basis construction relies on techniques such as the Greedy algorithm or proper orthogonal decomposition (POD). We demonstrate the developed procedure using an illustrative model problem on a three-dimensional geometry featuring a multi-dimensional geometrical parameterization.

academic

पैरामीटराइज्ड NURBS मल्टीपैच ज्यामिति पर अनअफाइन समस्याओं का कम किया गया क्रम मॉडलिंग

मूल जानकारी

पेपर ID: 2211.07348
शीर्षक: पैरामीटराइज्ड NURBS मल्टीपैच ज्यामिति पर अनअफाइन समस्याओं का कम किया गया क्रम मॉडलिंग
लेखक: मार्गरिटा चसापी, पाब्लो एंटोलिन, अन्नालिसा बुफा
वर्गीकरण: math.NA cs.NA
प्रकाशन समय: 14 नवंबर 2022 (arXiv सबमिशन)
पेपर लिंक: https://arxiv.org/abs/2211.07348

सारांश

यह पेपर पैरामीटराइज्ड आंशिक अवकल समीकरणों के समाधान में कम किए गए आधार विधियों (Reduced Basis Methods) और समतापीय विश्लेषण (IsoGeometric Analysis, IGA) के संयुक्त अनुप्रयोग की खोज करता है। समतापीय विश्लेषण का परिचय एकल ज्यामितीय प्रतिनिधित्व के आधार पर डिजाइन और विश्लेषण के एकीकृत सिमुलेशन ढांचे को लागू करता है। विशेष रूप से कई पैच द्वारा मॉडल की गई जटिल ज्यामिति के लिए, यह पेपर उच्च-आयामी पैरामीटर स्पेस विशेषताओं वाली अनअफाइन समस्याओं पर ध्यान केंद्रित करता है। अनुभवजन्य प्रक्षेप विधि (EIM) का उपयोग करके अफाइन पैरामीटर निर्भरता को पुनः प्राप्त किया जाता है, जिसे डोमेन विघटन के साथ आयाम को कम करने के लिए जोड़ा जाता है। पैरामीटराइज्ड सेटिंग में स्पलाइन पैच को युग्मित किया जाता है, स्थिर संघनन कम किए गए आधार तत्व (SCRBE) विधि को नियोजित किया जाता है। आसन्न पैच के सामान्य इंटरफेस पर स्थिर संघनन प्रक्रिया लागू की जाती है, जबकि आंतरिक रूप से कम किए गए आधार सन्निकटन का उपयोग करके कुशल ऑफलाइन/ऑनलाइन विघटन को लागू किया जाता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या परिभाषा

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

अनुसंधान प्रेरणा

एकीकृत ढांचा: IGA का उपयोग करके ज्यामितीय डिजाइन और संख्यात्मक विश्लेषण का एकीकृत प्रतिनिधित्व
गणनात्मक दक्षता: कम किए गए आधार विधियों के माध्यम से पैरामीटराइज्ड समस्याओं का तीव्र समाधान
औद्योगिक प्रासंगिकता: जटिल मल्टीपैच ज्यामिति समस्याओं को संभालना जिनका औद्योगिक महत्व है

मूल योगदान

विधि संलयन: पहली बार IGA और कम किए गए आधार विधियों को अनअफाइन मल्टीपैच ज्यामिति समस्याओं के लिए व्यवस्थित रूप से संयोजित किया गया
EIM-IGA एकीकरण: स्पलाइन सन्निकटन ढांचे के तहत अनुभवजन्य प्रक्षेप विधि को पुनः तैयार किया गया, सटीक ज्यामितीय मानचित्रण का उपयोग करके अफाइन विघटन को पुनः प्राप्त किया गया
SCRBE विस्तार: स्थिर संघनन कम किए गए आधार तत्व विधि को NURBS मल्टीपैच ज्यामिति तक विस्तारित किया गया
डोमेन विघटन रणनीति: उच्च-आयामी पैरामीटर स्पेस को संभालने के लिए पैच-आधारित डोमेन विघटन विधि प्रस्तावित की गई
संख्यात्मक सत्यापन: त्रि-आयामी मल्टीपैच ज्यामिति पर विधि की प्रभावशीलता को सत्यापित किया गया

विधि विवरण

कार्य परिभाषा

पैरामीटराइज्ड अण्डाकार आंशिक अवकल समीकरण को हल करना: $-\Delta u = f \text{ in } \Omega(\mu)$ $u = 0 \text{ on } \Gamma_D(\mu)$ $\frac{\partial u}{\partial n} = 0 \text{ on } \Gamma_N(\mu)$

जहाँ $\Omega(\mu) \subset \mathbb{R}^d$ पैरामीटराइज्ड डोमेन है, $\mu \in \mathcal{P} \subset \mathbb{R}^P$ ज्यामितीय पैरामीटर हैं।

मॉडल आर्किटेक्चर

1. NURBS मल्टीपैच ज्यामिति प्रतिनिधित्व

पैरामीटर डोमेन: $\hat{\Omega} = [0,1]^{\hat{d}}$
ज्यामितीय मानचित्रण: $F(\xi;\mu) = \sum_i B_{i,p}(\xi)P_i(\mu)$
मल्टीपैच विघटन: $\Omega(\mu) = \bigcup_{k=1}^{n_p} \Omega_k(\mu^{(k)})$

2. अनुभवजन्य प्रक्षेप विधि (EIM)

अनअफाइन पैरामीटर निर्भरता को संदर्भ डोमेन पर प्रक्षेप समस्या में परिवर्तित करना:

रूपांतरित कमजोर रूप: $\alpha(u_h, v_h;\mu) = \int_{\hat{\Omega}} (\nabla \hat{u}_h DF^{-1}) \cdot (\nabla \hat{v}_h DF^{-1})|\det DF| d\hat{\Omega}$

EIM सन्निकटन: $g_M^\alpha(\xi;\mu) = I_{\xi}^M g^\alpha(\xi;\mu) = \sum_{m=1}^{M^\alpha} \theta_m^\alpha(\mu)\phi_m^\alpha(\xi)$

3. SCRBE विधि

स्थिर संघनन विघटन:

बुलबुला स्पेस: $V_{h,b}^k = \{v_h \in V_h^k : v_h|_{\gamma_{k,j}} = 0\}$
पोर्ट स्पेस: $V_{h,\gamma}^{k,j} = \{v_h \in V_h^k : v_h|_{\gamma_{k,j}} \neq 0\}$

समाधान का प्रतिनिधित्व: $u_h(\mu)|_{\Omega_k} = b_{h,f}^k(\mu^{(k)}) + \sum_{j=1}^{n_\gamma^k} \sum_{r=1}^{N_{\gamma}^{k,j}} \hat{u}_{G_k(j),r}(\mu) \phi_{h}^{k,j,r}(\mu^{(k)})$

तकनीकी नवाचार बिंदु

ज्यामितीय मानचित्रण उपयोग: NURBS के सटीक ज्यामितीय मानचित्रण का सीधे उपयोग करके चर परिवर्तन किया जाता है, जो परंपरागत परिमित तत्वों में जटिल मानचित्रण निर्माण से बचता है
स्तरीय आयाम में कमी:
- पैच के अंदर: बुलबुला कार्यों को सन्निकटित करने के लिए कम किए गए आधार का उपयोग
- पैच इंटरफेस: पोर्ट मोड का निर्माण करने के लिए POD का उपयोग
अफाइनिटी पुनः प्राप्ति: पैरामीटर डोमेन पर EIM प्रक्षेप के माध्यम से ज्यामितीय पैरामीटराइजेशन की अफाइन विघटन को पुनः प्राप्त किया जाता है

प्रयोगात्मक सेटअप

ज्यामितीय मॉडल

डोमेन: त्रि-आयामी वक्र मल्टीपैच ज्यामिति
पैरामीटर: 8 ज्यामितीय पैरामीटर $\mu = (\mu_1, \mu_2, ..., \mu_8)$
पैच वितरण: 4 पैच, प्रत्येक पैच में 2 पैरामीटर
जाल: प्रत्येक पैच में 64 तत्व, कुल 1225 स्वतंत्रता की डिग्री
बहुपद क्रम: $p = 3$

सीमा शर्तें

डिरिचलेट सीमा: $u = 0$ on $\Gamma_D$
न्यूमैन सीमा: $\frac{\partial u}{\partial n} = 0$ on $\Gamma_N$
स्रोत पद: $\tilde{f} = 2xyz$

प्रशिक्षण पैरामीटर

EIM प्रशिक्षण सेट: $|P_{train}^{EIM}| = 250$
RB प्रशिक्षण सेट: $|P_{train}| = 250$
POD स्नैपशॉट: $N_s = 25$ (लैटिन हाइपरक्यूब नमूनाकरण)
सहनशीलता: EIM $\varepsilon = 10^{-7}$ , RB $\varepsilon = 10^{-5}$

प्रयोगात्मक परिणाम

EIM सन्निकटन प्रदर्शन

पैच	$M^\alpha$	$M^f$
1	28	13
2	33	14
3	19	10
4	23	11

कठोरता मैट्रिक्स और दाहिने पक्ष के EIM आधार कार्यों की संख्या कम है ( $M^\alpha \leq 33, M^f \leq 14$ )
$10^{-7}$ सटीकता पर तीव्र अभिसरण

कम किए गए आधार प्रदर्शन

पैच	$N_b$	ऑफलाइन समय (मिनट)	ऑनलाइन समय (मिली सेकंड)
1	17	2.16	44.2
2	16	1.72	41.2
3	10	1.01	46.4
4	11	1.70	36.8

पोर्ट मोड विश्लेषण

विलक्षण मान तीव्र क्षय: सभी इंटरफेस के POD विलक्षण मान तीव्र क्षय दिखाते हैं
पोर्ट आयाम: प्रत्येक इंटरफेस को ≤25 पोर्ट मोड तक कम किया जाता है
इंटरफेस RB: प्रत्येक पोर्ट मोड के लिए औसतन 16 आधार कार्य

गणनात्मक दक्षता

औसत ऑनलाइन मूल्यांकन: प्रत्येक पैच के लिए 42 मिली सेकंड
पोर्ट मोड ऑनलाइन समय: 29.3 मिली सेकंड
समग्र त्वरण अनुपात: पूर्ण-क्रम मॉडल की तुलना में महत्वपूर्ण सुधार

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

विधि व्यवहार्यता: अनअफाइन मल्टीपैच ज्यामिति समस्याओं को संभालने के लिए IGA और RB विधियों को सफलतापूर्वक संयोजित किया गया
गणनात्मक दक्षता: महत्वपूर्ण आयाम में कमी और त्वरण प्रभाव प्राप्त किए गए
सटीकता संरक्षण: छोटे कम किए गए आधार आयाम पर अच्छी सन्निकटन सटीकता बनाई गई
विस्तारशीलता: विधि ढांचे में अधिक जटिल औद्योगिक समस्याओं तक विस्तारित होने की क्षमता है

सीमाएं

समस्या जटिलता: वर्तमान में केवल अपेक्षाकृत सरल मॉडल समस्याओं पर सत्यापित
पैरामीटर आयाम: हालांकि 8-आयामी पैरामीटर को संभाला गया है, उच्च-आयामी पैरामीटर स्पेस पर प्रदर्शन को आगे के अनुसंधान की आवश्यकता है
ज्यामितीय जटिलता: अधिक जटिल औद्योगिक ज्यामिति पर प्रयोज्यता को सत्यापित करने की आवश्यकता है
अरैखिक विस्तार: वर्तमान ढांचा मुख्य रूप से रैखिक समस्याओं के लिए है

भविष्य की दिशाएं

औद्योगिक अनुप्रयोग: जटिल ज्यामिति समस्याओं तक विस्तार जिनका औद्योगिक प्रासंगिकता है
पैरामीटर अनुकूलन: पैरामीटर अनुकूलन प्रक्रिया में एकीकरण
अरैखिक समस्याएं: अरैखिक आंशिक अवकल समीकरणों तक विस्तार
समानांतर कंप्यूटिंग: अधिक कुशल समानांतर कार्यान्वयन विकसित करना

गहन मूल्यांकन

लाभ

सैद्धांतिक पूर्णता: IGA-RB-EIM के सैद्धांतिक ढांचे को व्यवस्थित रूप से स्थापित किया गया
विधि नवाचार: NURBS ज्यामितीय मानचित्रण का कुशलतापूर्वक उपयोग करके EIM कार्यान्वयन को सरल बनाया गया
व्यावहारिक मूल्य: जटिल ज्यामिति के पैरामीटराइज्ड समाधान के लिए प्रभावी मार्ग प्रदान किया गया
सत्यापन पर्याप्त: संख्यात्मक प्रयोगों के माध्यम से विभिन्न घटकों की प्रभावशीलता को सत्यापित किया गया

कमियाँ

अनुप्रयोग सीमा: सत्यापन मामले अपेक्षाकृत सरल हैं, जटिल औद्योगिक समस्याओं का परीक्षण अभाव है
सैद्धांतिक विश्लेषण: अभिसरण और त्रुटि सीमा का कठोर सैद्धांतिक विश्लेषण अभाव है
तुलनात्मक अध्ययन: अन्य कम किए गए क्रम विधियों के साथ व्यवस्थित तुलना अपर्याप्त है
गणनात्मक जटिलता: विस्तृत गणनात्मक जटिलता विश्लेषण प्रदान नहीं किया गया है

प्रभाव

शैक्षणिक योगदान: IGA और RB विधियों के संयोजन के लिए महत्वपूर्ण सैद्धांतिक आधार प्रदान किया गया
व्यावहारिक मूल्य: CAD/CAE एकीकरण के लिए महत्वपूर्ण है
पुनरुत्पादनशीलता: खुले स्रोत सॉफ्टवेयर पैकेज का उपयोग करके, अच्छी पुनरुत्पादनशीलता है
विस्तार क्षमता: बाद के अनुसंधान के लिए ठोस पद्धति संबंधी आधार प्रदान किया गया

प्रयोज्य परिदृश्य

ज्यामितीय डिजाइन अनुकूलन: CAD पैरामीटर-संचालित डिजाइन अनुकूलन
बहु-प्रश्न सिमुलेशन: बड़ी संख्या में पैरामीटर मूल्यांकन की आवश्यकता वाली इंजीनियरिंग समस्याएं
वास्तविक समय सिमुलेशन: गणनात्मक समय पर कठोर आवश्यकताओं वाले अनुप्रयोग
आकार अनुकूलन: ज्यामितीय पैरामीटर पर आधारित संरचनात्मक अनुकूलन समस्याएं

संदर्भ

यह पेपर कम किए गए आधार विधियों, समतापीय विश्लेषण और डोमेन विघटन के क्षेत्रों से महत्वपूर्ण साहित्य का हवाला देता है, जिसमें शामिल हैं:

Hughes आदि 2005: IGA आधार सिद्धांत
Quarteroni आदि 2016: कम किए गए आधार विधियों का सर्वेक्षण
Barrault आदि 2004: अनुभवजन्य प्रक्षेप विधि
Huynh आदि 2013: SCRBE विधि आधार

समग्र मूल्यांकन: यह कम्प्यूटेशनल गणित के क्षेत्र में महत्वपूर्ण महत्व का एक पेपर है, जो जटिल ज्यामितीय पैरामीटराइजेशन समस्याओं के कुशल समाधान को व्यवस्थित रूप से हल करता है। विधि में ठोस सैद्धांतिक आधार और अच्छी व्यावहारिक संभावनाएं हैं, जो IGA और RB विधियों के संयोजन के लिए अनुसंधान की नई दिशा खोलता है।