2025-11-12T12:07:10.593714

A localized reduced basis approach for unfitted domain methods on parameterized geometries

Chasapi, Antolin, Buffa

This work introduces a reduced order modeling (ROM) framework for the solution of parameterized second-order linear elliptic partial differential equations formulated on unfitted geometries. The goal is to construct efficient projection-based ROMs, which rely on techniques such as the reduced basis method and discrete empirical interpolation. The presence of geometrical parameters in unfitted domain discretizations entails challenges for the application of standard ROMs. Therefore, in this work we propose a methodology based on i) extension of snapshots on the background mesh and ii) localization strategies to decrease the number of reduced basis functions. The method we obtain is computationally efficient and accurate, while it is agnostic with respect to the underlying discretization choice. We test the applicability of the proposed framework with numerical experiments on two model problems, namely the Poisson and linear elasticity problems. In particular, we study several benchmarks formulated on two-dimensional, trimmed domains discretized with splines and we observe a significant reduction of the online computational cost compared to standard ROMs for the same level of accuracy. Moreover, we show the applicability of our methodology to a three-dimensional geometry of a linear elastic problem.

academic

অপরিমাপিত ডোমেইন পদ্ধতিতে প্যারামিটারাইজড জ্যামিতির জন্য স্থানীয়করণ হ্রাসকৃত ভিত্তি পদ্ধতি

মৌলিক তথ্য

পত্র আইডি: 2212.11934
শিরোনাম: অপরিমাপিত ডোমেইন পদ্ধতিতে প্যারামিটারাইজড জ্যামিতির জন্য স্থানীয়করণ হ্রাসকৃত ভিত্তি পদ্ধতি
লেখক: মার্গারিটা চাসাপি, পাবলো অ্যান্টোলিন, অ্যানালিসা বুফা (ইকোল পলিটেকনিক ফেডারেল ডি লজান ও সিএনআর)
শ্রেণীবিভাগ: math.NA cs.NA
প্রকাশিত জার্নাল: কম্পিউটার মেথডস ইন অ্যাপ্লাইড মেকানিক্স এন্ড ইঞ্জিনিয়ারিং (গৃহীত)
পত্র লিঙ্ক: https://arxiv.org/abs/2212.11934

সারসংক্ষেপ

এই পত্রটি একটি হ্রাসকৃত ক্রম মডেলিং (ROM) কাঠামো প্রস্তাব করে যা অপরিমাপিত জ্যামিতিতে সংজ্ঞায়িত প্যারামিটারাইজড দ্বিতীয় ক্রমের রৈখিক উপবৃত্তাকার আংশিক ডিফারেনশিয়াল সমীকরণ সমাধানের জন্য। এই পদ্ধতিটি হ্রাসকৃত ভিত্তি পদ্ধতি এবং বিচ্ছিন্ন অভিজ্ঞতামূলক ইন্টারপোলেশনের মতো কৌশল ব্যবহার করে দক্ষ প্রজেকশন-ভিত্তিক ROM তৈরি করার লক্ষ্য রাখে। জ্যামিতিগত প্যারামিটারগুলি অপরিমাপিত ডোমেইন বিচ্ছিন্নকরণে উপস্থিত থাকার কারণে মান ROM-এর প্রয়োগে চ্যালেঞ্জ তৈরি হয়, এই পত্রটি (i) পটভূমি জালে স্ন্যাপশট সম্প্রসারণ এবং (ii) হ্রাসকৃত ভিত্তি ফাংশনের সংখ্যা হ্রাস করার জন্য স্থানীয়করণ কৌশলের উপর ভিত্তি করে পদ্ধতিবিদ্যা প্রস্তাব করে। ফলস্বরূপ পদ্ধতিটি গণনামূলকভাবে দক্ষ এবং নির্ভুল, যখন অন্তর্নিহিত বিচ্ছিন্নকরণ পছন্দের প্রতি অজ্ঞেয়বাদী।

গবেষণা পটভূমি এবং প্রেরণা

মূল সমস্যা

অপরিমাপিত ডোমেইন পদ্ধতির চ্যালেঞ্জ: কম্পিউটার-সহায়ক ডিজাইন (CAD) এবং আইসোজিওমেট্রিক বিশ্লেষণ (IGA)-তে, জটিল জ্যামিতি সাধারণত ট্রিমিং অপারেশনের মাধ্যমে পরিচালিত হয়, যার ফলে অপরিমাপিত জাল তৈরি হয় যেখানে জ্যামিতি সমাধানের বিচ্ছিন্নকরণ থেকে আলাদা।
প্যারামিটারাইজড জ্যামিতির হ্রাসকৃত ক্রম মডেলিং কঠিনতা: যখন জ্যামিতিগত প্যারামিটার পরিবর্তিত হয়, সক্রিয় ভিত্তি ফাংশনের সেট পরিবর্তিত হয়, যার ফলে স্ন্যাপশট সমাধান ভেক্টরগুলি বিভিন্ন দৈর্ঘ্যের হতে পারে, যা মান হ্রাসকৃত ভিত্তি পদ্ধতির সরাসরি প্রয়োগে বাধা দেয়।
অ-সম্বন্ধিত প্যারামিটার নির্ভরতা: জ্যামিতিগত প্যারামিটারগুলি প্যারামিটারের প্রতি ডিফারেনশিয়াল অপারেটরের অ-সম্বন্ধিত নির্ভরতা সৃষ্টি করে, যা হ্রাসকৃত ভিত্তি পদ্ধতির মৌলিক অনুমানগুলি ভেঙে দেয়।

গবেষণার গুরুত্ব

শিল্প প্রাসঙ্গিকতা: ট্রিমিং CAD-তে জটিল আকৃতি প্রতিনিধিত্বের একটি প্রধান কৌশল, একীভূত ডিজাইন-থেকে-বিশ্লেষণ কাঠামো বাস্তবায়নের জন্য গুরুত্বপূর্ণ
বহু-প্রশ্ন পরিস্থিতি: ডিজাইন অপ্টিমাইজেশন, অনিশ্চয়তা পরিমাণীকরণ, নিয়ন্ত্রণ এবং অন্যান্য প্রয়োগে প্যারামিটারাইজড PDE একাধিকবার সমাধান করার প্রয়োজন
গণনামূলক দক্ষতার প্রয়োজন: গণনামূলক ত্বরণ অর্জনের জন্য দক্ষ অফলাইন/অনলাইন পদ্ধতির প্রয়োজন

বিদ্যমান পদ্ধতির সীমাবদ্ধতা

মান ROM-এর জন্য একটি রেফারেন্স ডোমেইন ম্যাপিং প্রয়োজন, যা নির্দিষ্ট সমস্যার উপর অত্যন্ত নির্ভরশীল
বিদ্যমান অপরিমাপিত ডোমেইন ROM কাজ অ-সম্বন্ধিত নির্ভরতা সমস্যা সমাধান করে না, অতি-হ্রাসকরণ কৌশল অনুপস্থিত
বৈশ্বিক হ্রাসকৃত ভিত্তি সম্প্রসারিত সমাধানের অ-রৈখিক বহুগুণ অনুমানের জন্য দুর্বল, বড় সংখ্যক ভিত্তি ফাংশন প্রয়োজন

মূল অবদান

প্যারামিটারাইজড অপরিমাপিত জ্যামিতির জন্য সম্পূর্ণ হ্রাসকরণ কাঠামো প্রস্তাব করা, অতি-হ্রাসকরণ কৌশল একীভূত করে দক্ষ অফলাইন/অনলাইন বিভাজন অর্জন করা
প্যারামিটার-চালিত স্থানীয়করণ কৌশল বিকাশ করা, k-মানস ক্লাস্টারিং ব্যবহার করে একাধিক ছোট-মাত্রার স্থানীয় সাবস্পেস তৈরি করা, হ্রাসকৃত ভিত্তি ফাংশনের সংখ্যা উল্লেখযোগ্যভাবে হ্রাস করা
বিচ্ছিন্নকরণ পদ্ধতি-অজ্ঞেয়বাদী সার্বজনীন কাঠামো ডিজাইন করা, সীমিত উপাদান এবং আইসোজিওমেট্রিক পদ্ধতি উভয়ের জন্য প্রযোজ্য
রেডিয়াল ভিত্তি ফাংশন (RBF) ইন্টারপোলেশন-ভিত্তিক অ-আক্রমণাত্মক অনলাইন সহগ গণনা বাস্তবায়ন করা, উচ্চ-নির্ভুলতা সমাবেশ রুটিনে আক্রমণাত্মক পরিবর্তন এড়ানো
সংখ্যাগত পরীক্ষায় উল্লেখযোগ্য গণনামূলক ত্বরণ প্রদর্শন করা: মান ROM-এর তুলনায়, একই নির্ভুলতা স্তরে 17.6 গুণ থেকে 130 গুণ অনলাইন গণনা ত্বরণ অর্জন করা

পদ্ধতি বিস্তারিত

কাজের সংজ্ঞা

প্যারামিটারাইজড অপরিমাপিত ডোমেইন Ω(μ)-তে সংজ্ঞায়িত রৈখিক উপবৃত্তাকার PDE সমাধান করা:

Lu = f in Ω(μ)

যেখানে μ ∈ P ⊂ R^M হল জ্যামিতিগত প্যারামিটার ভেক্টর, ডোমেইন কাটিং অপারেশনের মাধ্যমে সংজ্ঞায়িত:

Ω(μ) = Ω̂₀(μ) \ ⋃ᵢ₌₁ᴷ Ω̂ᵢ(μ)

মূল পদ্ধতি স্থাপত্য

1. অপরিমাপিত ডোমেইন বিচ্ছিন্নকরণ

পটভূমি ডোমেইন সেটআপ: প্যারামিটারাইজড ডোমেইন Ω(μ) একটি প্যারামিটার-স্বাধীন পটভূমি ডোমেইন Ω₀-তে এম্বেড করা
বিচ্ছিন্ন স্পেস নির্মাণ: পটভূমি ডোমেইনে বিচ্ছিন্ন ফাংশন স্পেস V_{h,0} সংজ্ঞায়িত করা, সক্রিয় স্পেস V_h(μ) সমর্থন Ω(μ)-এর সাথে ছেদ করে এমন ভিত্তি ফাংশন দ্বারা বিস্তৃত
সম্প্রসারণ কৌশল: বিভিন্ন প্যারামিটারের অধীনে সমাধান ভেক্টরগুলি পটভূমি ডোমেইনে সম্প্রসারিত করা, অ-সক্রিয় অঞ্চলে শূন্য সেট করা

2. বিচ্ছিন্ন অভিজ্ঞতামূলক ইন্টারপোলেশন পদ্ধতি (DEIM)

সম্বন্ধিত বিয়োজন: কঠোরতা ম্যাট্রিক্স এবং ডান দিকের ভেক্টর সম্বন্ধিত ফর্মে অনুমান করা:
```
Â(μ) ≈ ∑_{q=1}^{Q_a} θ_q^a(μ)Âq
f̂(μ) ≈ ∑_{q=1}^{Q_f} θ_q^f(μ)f̂q
```
RBF ইন্টারপোলেশন: অনলাইনে দ্রুত প্যারামিটার-নির্ভর সহগ গণনা করতে রেডিয়াল ভিত্তি ফাংশন ব্যবহার করা:
```
θ_q^a(μ) ≈ ∑_{j=1}^{N_s^d} ω_{q,j}^a φ_{q,j}(‖μ-μⱼ‖₂)
```

3. স্থানীয়করণ কৌশল

প্যারামিটার স্পেস ক্লাস্টারিং: প্যারামিটার স্পেসকে N_c সাব-অঞ্চলে বিভক্ত করতে k-মানস অ্যালগরিদম ব্যবহার করা
স্থানীয় ভিত্তি নির্মাণ: প্রতিটি ক্লাস্টারের জন্য আলাদাভাবে DEIM অনুমান এবং হ্রাসকৃত ভিত্তি তৈরি করা
অনলাইন সুইচিং: প্যারামিটার থেকে ক্লাস্টার কেন্দ্রের দূরত্বের উপর ভিত্তি করে সংশ্লিষ্ট স্থানীয় ভিত্তি নির্বাচন করা

প্রযুক্তিগত উদ্ভাবন পয়েন্ট

জ্যামিতি-চালিত ক্লাস্টারিং কৌশল: সীমানা দূরত্বের উপর ভিত্তি করে প্যারামিটার ক্লাস্টারিং, সমান জ্যামিতি কনফিগারেশন গ্রুপিং নিশ্চিত করা
বহু-স্তরের স্থানীয়করণ: DEIM অনুমান এবং হ্রাসকৃত ভিত্তি উভয়ের জন্য একযোগে স্থানীয়করণ
অ-আক্রমণাত্মক বাস্তবায়ন: RBF ইন্টারপোলেশনের মাধ্যমে অন্তর্নিহিত সমাবেশ রুটিনে পরিবর্তন এড়ানো
সার্বজনীন কাঠামো ডিজাইন: নির্দিষ্ট বিচ্ছিন্নকরণ পদ্ধতি থেকে আলাদা, একাধিক অপরিমাপিত ডোমেইন কৌশল সমর্থন করা

পরীক্ষামূলক সেটআপ

পরীক্ষার সমস্যা

পয়সন সমীকরণ:
- 1D প্যারামিটারাইজেশন: বৃত্তাকার গর্ত সহ বর্গাকার ডোমেইন, গর্তের কেন্দ্র কর্ণ বরাবর চলমান
- 2D প্যারামিটারাইজেশন: গর্তের কেন্দ্র অবস্থান এবং ব্যাসার্ধ একসাথে পরিবর্তিত
রৈখিক স্থিতিস্থাপকতা সমস্যা:
- ছিদ্রযুক্ত ত্রৈমাসিক সিলিন্ডার: একাধিক বৃত্তাকার গর্ত ব্যাসার্ধ প্যারামিটারাইজেশন
- 3D ঘনক গোলাকার অন্তর্ভুক্তি: গোলক ব্যাসার্ধ প্যারামিটারাইজেশন

বিচ্ছিন্নকরণ সেটআপ

ভিত্তি ফাংশন: দ্বিতীয় এবং তৃতীয় ক্রমের C¹-C² ক্রমাগত B-স্প্লাইন
জাল: 32×32 কার্টেসিয়ান জাল (2D), 8×8×8 জাল (3D)
স্বাধীনতার ডিগ্রি: 1225-3000 DOF

মূল্যায়ন সূচক

নির্ভুলতা: L² এবং H¹ নর্ম আপেক্ষিক ত্রুটি
দক্ষতা: অনলাইন CPU সময় এবং ত্বরণ অনুপাত
হ্রাসকরণ প্রভাব: একবচন মূল্য ক্ষয় এবং ভিত্তি ফাংশন সংখ্যা

প্যারামিটার সেটিংস

POD সহনশীলতা: ε_ = 10⁻⁵, ε_^d = 10⁻⁷
প্রশিক্ষণ নমুনা: N_s = 100-1000, N_s^d = 250-1000
ক্লাস্টার সংখ্যা: N_c = 4-16, N_c^d = 8-16
RBF প্রকার: ঘন রেডিয়াল ভিত্তি ফাংশন

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

1. DEIM অনুমান কার্যকারিতা

বৈশ্বিক পদ্ধতি: 349টি ম্যাট্রিক্স পদ এবং 124টি ভেক্টর পদ প্রয়োজন (1D ক্ষেত্রে)
স্থানীয় পদ্ধতি: সর্বাধিক 31-38টি ম্যাট্রিক্স পদ এবং 11-15টি ভেক্টর পদ (16টি ক্লাস্টার)
নির্ভুলতা: 10⁻⁵ স্তরের ত্রুটি অর্জন করা

2. হ্রাসকৃত ভিত্তি মাত্রা তুলনা

সমস্যা	বৈশ্বিক ROM	স্থানীয় ROM	ক্লাস্টার সংখ্যা
পয়সন 1D	N=182	N=35	4
পয়সন 2D	N=201	N=17	16
স্থিতিস্থাপকতা ছিদ্রযুক্ত	N=198	N=16	10

3. গণনামূলক ত্বরণ অনুপাত

সমস্যা	বৈশ্বিক ROM ত্বরণ	স্থানীয় ROM ত্বরণ
পয়সন 2D	9.1×	17.6×
স্থিতিস্থাপকতা ছিদ্রযুক্ত	16.5×	130×

বিলোপন পরীক্ষা

ক্লাস্টার সংখ্যা প্রভাব: ক্লাস্টার সংখ্যা বৃদ্ধি ভিত্তি ফাংশন সংখ্যা উল্লেখযোগ্যভাবে হ্রাস করে, নির্ভুলতা উন্নত করে
প্রশিক্ষণ নমুনা আকার: DEIM নির্ভুলতা নিশ্চিত করতে যথেষ্ট সমৃদ্ধ প্রশিক্ষণ সেট প্রয়োজন (N_s^d ≥ 100)
RBF ইন্টারপোলেশন নির্ভুলতা: ঘন RBF ভাল মসৃণতা প্রদান করে কোনো মিথ্যা দোলন ছাড়াই

3D যাচাইকরণ

3D ঘনক গোলাকার অন্তর্ভুক্তি সমস্যায়:

DEIM ভিত্তি ফাংশন 180 থেকে 35-এ হ্রাস
হ্রাসকৃত ভিত্তি মাত্রা 35 থেকে 14-এ হ্রাস
উচ্চ-মাত্রার জ্যামিতিতে পদ্ধতির কার্যকারিতা যাচাই করা

উপসংহার এবং আলোচনা

প্রধান উপসংহার

কার্যকারিতা যাচাইকরণ: প্রস্তাবিত পদ্ধতি নির্ভুলতা বজায় রেখে অনলাইন গণনা খরচ উল্লেখযোগ্যভাবে হ্রাস করে
সার্বজনীনতা: কাঠামো অন্তর্নিহিত বিচ্ছিন্নকরণ পদ্ধতি থেকে স্বাধীন, একাধিক অপরিমাপিত ডোমেইন কৌশলের জন্য প্রযোজ্য
সম্প্রসারণযোগ্যতা: জটিল 3D জ্যামিতিতে সফলভাবে সম্প্রসারিত, ভাল স্কেলেবিলিটি প্রদর্শন করা

সীমাবদ্ধতা

ক্লাস্টার সংখ্যা নির্বাচন: ক্লাস্টার সংখ্যা পূর্ব-নির্ধারণের প্রয়োজন, যদিও বৈচিত্র্য মানদণ্ড প্রদান করা হয় তবুও অভিজ্ঞতামূলক বিচার প্রয়োজন
প্যারামিটার মাত্রা: উচ্চ-মাত্রার প্যারামিটার স্পেস RBF ইন্টারপোলেশন শর্ত সংখ্যা সমস্যা হতে পারে
জ্যামিতি জটিলতা: অত্যন্ত জটিল জ্যামিতি পরিবর্তনের জন্য, আরও ক্লাস্টার প্রয়োজন হতে পারে

ভবিষ্যত দিকনির্দেশনা

লোভী অ্যালগরিদম: স্থানীয় হ্রাসকৃত ভিত্তি এবং ত্রুটি সার্টিফিকেশন তৈরিতে লোভী অ্যালগরিদম প্রয়োগ করা
চতুর্থ ক্রম PDE: আরও জটিল PDE সমস্যায় সম্প্রসারণ করা
অপ্টিমাইজেশন প্রয়োগ: ডিজাইন, আকৃতি এবং টপোলজি অপ্টিমাইজেশনে প্রয়োগ করা

গভীর মূল্যায়ন

শক্তি

পদ্ধতি উদ্ভাবন শক্তিশালী: অপরিমাপিত ডোমেইন ROM-এ স্থানীয়করণ কৌশল প্রথমবারের মতো সিস্টেমেটিকভাবে প্রয়োগ করা, গুরুত্বপূর্ণ প্রযুক্তিগত সমস্যা সমাধান করা
তাত্ত্বিক কাঠামো সম্পূর্ণ: স্ন্যাপশট সম্প্রসারণ থেকে স্থানীয়করণ কৌশল পর্যন্ত, সম্পূর্ণ পদ্ধতিবিদ্যা সিস্টেম গঠন করা
পরীক্ষামূলক যাচাইকরণ পর্যাপ্ত: 2D/3D, বিভিন্ন PDE প্রকার, একাধিক প্যারামিটারাইজেশন পরিস্থিতি জুড়ে ব্যাপক যাচাইকরণ
ব্যবহারিক মূল্য উচ্চ: উল্লেখযোগ্য গণনা ত্বরণ এবং বিচ্ছিন্নকরণ পদ্ধতি-অজ্ঞেয়বাদী বৈশিষ্ট্য এটিকে বিস্তৃত প্রয়োগ সম্ভাবনা প্রদান করে

অপূর্ণতা

সীমিত তাত্ত্বিক বিশ্লেষণ: সংগ্রহ এবং ত্রুটি অনুমানের কঠোর তাত্ত্বিক বিশ্লেষণের অভাব
প্যারামিটার নির্ভরতা: পদ্ধতি কার্যকারিতা ক্লাস্টারিং কৌশল এবং প্যারামিটার নির্বাচনের উপর অত্যন্ত নির্ভরশীল
জটিলতা ভারসাম্য: যদিও অনলাইন দক্ষতা উচ্চ, অফলাইন পর্যায়ের গণনা এবং সংরক্ষণ খরচ উল্লেখযোগ্যভাবে বৃদ্ধি পায়

প্রভাব

একাডেমিক অবদান: অপরিমাপিত ডোমেইন ROM-এর জন্য নতুন গবেষণা দিক খোলা, পরবর্তী গবেষণা উদ্দীপিত করার প্রত্যাশা
প্রকৌশল মূল্য: CAD-বিশ্লেষণ একীকরণ এবং শিল্প জ্যামিতি অপ্টিমাইজেশনের জন্য গুরুত্বপূর্ণ
পুনরুৎপাদনযোগ্যতা: বিস্তারিত অ্যালগরিদম বর্ণনা এবং প্যারামিটার সেটিংস প্রদান করা, পুনরুৎপাদন এবং প্রয়োগ সহজতর করা

প্রযোজ্য পরিস্থিতি

জটিল জ্যামিতি পরিবর্তন জড়িত প্যারামিটারাইজড ডিজাইন অপ্টিমাইজেশন
ট্রিমড জ্যামিতির বহু-প্রশ্ন বিশ্লেষণ
রিয়েল-টাইম সিমুলেশন এবং নিয়ন্ত্রণ প্রয়োগ
অনিশ্চয়তা পরিমাণীকরণে জ্যামিতি প্যারামিটার গবেষণা

সংদর্ভ

পত্রটি সমৃদ্ধ সম্পর্কিত সাহিত্য উদ্ধৃত করে, যার মধ্যে রয়েছে:

অপরিমাপিত ডোমেইন পদ্ধতির ক্লাসিক কাজ (পেস্কিন, হিউজেস ইত্যাদি)
আইসোজিওমেট্রিক বিশ্লেষণের ভিত্তি সাহিত্য
হ্রাসকৃত ভিত্তি পদ্ধতি এবং অভিজ্ঞতামূলক ইন্টারপোলেশনের তাত্ত্বিক ভিত্তি
স্থানীয় ROM-এর অগ্রগামী গবেষণা

এই পত্রটি অপরিমাপিত ডোমেইনে প্যারামিটারাইজড PDE হ্রাসকৃত ক্রম মডেলিংয়ের জন্য গুরুত্বপূর্ণ পদ্ধতিবিদ্যা অবদান প্রদান করে, তাত্ত্বিক উদ্ভাবন এবং ব্যবহারিক প্রয়োগের মধ্যে ভাল ভারসাম্য অর্জন করে।