The Frenet apparatus is a new framework for constructing high order geometry-conforming immersed finite element functions for interface problems. In this report, we present a procedure for constructing the local IFE bases in some detail as well as a new approach for constructing orthonormal bases using the singular value decomposition of the local generalized Vandermonde matrix. A sample implementation in MATLAB is provided to showcase the simplicity and extensionability of the framework.
গবেষণাপত্র ID : 2510.12018শিরোনাম : Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Methodলেখক : Slimane Adjerid (Virginia Tech), Tao Lin (Virginia Tech), Haroun Meghaichi (Ohio State University)শ্রেণীবিভাগ : math.NA cs.NAপ্রকাশনার সময় : অক্টোবর ১৫, ২০২৫গবেষণাপত্র লিংক : https://arxiv.org/abs/2510.12018 এই গবেষণাপত্রে Frenet যন্ত্রের উপর ভিত্তি করে একটি নতুন কাঠামো প্রস্তাব করা হয়েছে, যা সীমানা সমস্যা সমাধানের জন্য উচ্চ-ক্রমের জ্যামিতি সামঞ্জস্যপূর্ণ নিমজ্জিত সীমান্ত উপাদান ফাংশন নির্মাণের জন্য। গবেষণাপত্রে স্থানীয় IFE ভিত্তি ফাংশন নির্মাণের প্রক্রিয়া এবং স্থানীয় সাধারণীকৃত Vandermonde ম্যাট্রিক্সের একবচন মান বিয়োজন ব্যবহার করে অর্থোনরমাল ভিত্তি নির্মাণের একটি নতুন পদ্ধতি বিস্তারিতভাবে বর্ণনা করা হয়েছে। MATLAB উদাহরণ বাস্তবায়ন প্রদান করা হয়েছে, যা এই কাঠামোর সরলতা এবং সম্প্রসারণযোগ্যতা প্রদর্শন করে।
এই গবেষণা উপবৃত্তাকার সীমানা সমস্যা সমাধান করে:
− ∇ ⋅ ( β ∇ u ) = f , on Ω − ∪ Ω + -\nabla \cdot (\beta\nabla u) = f, \text{ on } \Omega^- \cup \Omega^+ − ∇ ⋅ ( β ∇ u ) = f , on Ω − ∪ Ω +
যেখানে গণনা ক্ষেত্র Ω ⊂ R 2 \Omega \subset \mathbb{R}^2 Ω ⊂ R 2 সীমানা Γ \Gamma Γ দ্বারা দুটি উপক্ষেত্র Ω − \Omega^- Ω − এবং Ω + \Omega^+ Ω + তে বিভক্ত, এবং বিস্তার সহগ β \beta β একটি খণ্ডশ: ধ্রুবক ফাংশন।
বিদ্যমান IFE পদ্ধতির সীমাবদ্ধতা : ঐতিহ্যবাহী নিমজ্জিত সীমান্ত উপাদান পদ্ধতি শুধুমাত্র সীমানা লাফ শর্ত আনুমানিক সন্তুষ্ট করতে পারেউচ্চ-ক্রম নির্ভুলতার প্রয়োজন : যেকোনো ক্রমের বহুপদীর জন্য IFE ফাংশন নির্মাণের প্রয়োজনজ্যামিতি সামঞ্জস্যতার প্রয়োজনীয়তা : সীমানা লাফ শর্ত সঠিকভাবে সন্তুষ্ট করে এমন IFE ফাংশনের প্রয়োজনগণনামূলক স্থিতিশীলতা : ভালো শর্ত সংখ্যা সহ ভিত্তি ফাংশন সিস্টেমের প্রয়োজনসীমানা সমস্যা বহু-পর্যায়ের প্রবাহ, উপকরণ বিজ্ঞান, জৈব চিকিৎসা এবং অন্যান্য ক্ষেত্রে ব্যাপকভাবে দেখা যায়। উচ্চ-নির্ভুলতা সংখ্যাসূচক সমাধান পদ্ধতি প্রকৌশল প্রয়োগের জন্য গুরুত্বপূর্ণ।
Frenet যন্ত্রের উপর ভিত্তি করে GC-IFE নির্মাণ কাঠামো প্রস্তাব : অন্তরক জ্যামিতিতে Frenet ফ্রেম ব্যবহার করে সীমানা লাফ শর্ত সরল করাদুটি ভিত্তি ফাংশন নির্মাণ পদ্ধতি উন্নয়ন :
বিশেষ ফর্ম্যাটের উপর ভিত্তি করে প্রাথমিক নির্মাণ পদ্ধতি সম্প্রসারণ ধারণার উপর ভিত্তি করে সাধারণীকৃত নির্মাণ পদ্ধতি দুটি অর্থোগোনালাইজেশন পুনর্নির্মাণ অ্যালগরিদম প্রস্তাব :
ভর ম্যাট্রিক্স SVD এর উপর ভিত্তি করে পদ্ধতি (পদ্ধতি ১) সাধারণীকৃত Vandermonde ম্যাট্রিক্স SVD এর উপর ভিত্তি করে পদ্ধতি (পদ্ধতি ২) সম্পূর্ণ MATLAB বাস্তবায়ন প্রদান : পদ্ধতির ব্যবহারিকতা এবং সম্প্রসারণযোগ্যতা প্রদর্শনমূল ধারণা হল Frenet রূপান্তর ব্যবহার করে ভৌত স্থানাঙ্ক ( x , y ) (x,y) ( x , y ) কে স্থানীয় স্থানাঙ্ক ( η , ξ ) (\eta,\xi) ( η , ξ ) তে রূপান্তরিত করা:
[ η ξ ] → P Γ [ x ( η , ξ ) y ( η , ξ ) ] = g ( ξ ) + η n ( ξ ) \begin{bmatrix} \eta \\ \xi \end{bmatrix} \xrightarrow{P_\Gamma} \begin{bmatrix} x(\eta,\xi) \\ y(\eta,\xi) \end{bmatrix} = g(\xi) + \eta n(\xi) [ η ξ ] P Γ [ x ( η , ξ ) y ( η , ξ ) ] = g ( ξ ) + η n ( ξ )
যেখানে g ( ξ ) g(\xi) g ( ξ ) সীমানা বক্রের প্যারামিটারকরণ এবং n ( ξ ) n(\xi) n ( ξ ) স্বাভাবিক ভেক্টর।
স্থানীয় স্থানাঙ্ক সিস্টেমে, সীমানা লাফ শর্ত হয়ে ওঠে:
⟦ u ^ ⟧ Γ ^ K F = 0 \llbracket \hat{u} \rrbracket_{\hat{\Gamma}_{K_F}} = 0 [ [ u ^ ] ] Γ ^ K F = 0 ⟦ β ^ u ^ η ⟧ Γ ^ K F = 0 \llbracket \hat{\beta}\hat{u}_\eta \rrbracket_{\hat{\Gamma}_{K_F}} = 0 [ [ β ^ u ^ η ] ] Γ ^ K F = 0 ⟦ β ^ ∂ j ∂ η j L ( u ^ ) ⟧ Γ ^ K F = 0 , j = 0 , 1 , … , m − 2 \llbracket \hat{\beta}\frac{\partial^j}{\partial\eta^j}\mathcal{L}(\hat{u}) \rrbracket_{\hat{\Gamma}_{K_F}} = 0, j = 0,1,\ldots,m-2 [ [ β ^ ∂ η j ∂ j L ( u ^ ) ] ] Γ ^ K F = 0 , j = 0 , 1 , … , m − 2 নিম্নলিখিত ফর্মের ভিত্তি ফাংশন নির্মাণ করুন:
ϕ ^ i , j ( η , ξ ) = 1 β ^ ( η , ξ ) q j ( η η h ) p i ( ξ − ξ m i d ξ h ) \hat{\phi}_{i,j}(\eta,\xi) = \frac{1}{\hat{\beta}(\eta,\xi)}q_j\left(\frac{\eta}{\eta_h}\right)p_i\left(\frac{\xi-\xi_{mid}}{\xi_h}\right) ϕ ^ i , j ( η , ξ ) = β ^ ( η , ξ ) 1 q j ( η h η ) p i ( ξ h ξ − ξ mi d )
1 ≤ j ≤ m , 0 ≤ i ≤ m 1 \leq j \leq m, 0 \leq i \leq m 1 ≤ j ≤ m , 0 ≤ i ≤ m এর জন্য, এই ফাংশনগুলি স্বয়ংক্রিয়ভাবে লাফ শর্ত সন্তুষ্ট করে।
আরও সাধারণ ফর্ম গ্রহণ করুন:
λ ^ j ( η , ξ ) = { λ ^ j − ( η , ξ ) = ∑ i = 1 ( m + 1 ) 2 C i , j − R i ( η , ξ ) , η < 0 λ ^ j + ( η , ξ ) = ∑ i = 1 ( m + 1 ) 2 C i , j + R i ( η , ξ ) , η > 0 \hat{\lambda}_j(\eta,\xi) = \begin{cases}
\hat{\lambda}_j^-(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^- R_i(\eta,\xi), & \eta < 0 \\
\hat{\lambda}_j^+(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^+ R_i(\eta,\xi), & \eta > 0
\end{cases} λ ^ j ( η , ξ ) = { λ ^ j − ( η , ξ ) = ∑ i = 1 ( m + 1 ) 2 C i , j − R i ( η , ξ ) , λ ^ j + ( η , ξ ) = ∑ i = 1 ( m + 1 ) 2 C i , j + R i ( η , ξ ) , η < 0 η > 0
রৈখিক সমীকরণ সিস্টেম সমাধানের মাধ্যমে সহগ ম্যাট্রিক্স C ± C^{\pm} C ± নির্ধারণ করুন।
M q = V 1 Λ V 1 T M_q = V_1\Lambda V_1^T M q = V 1 Λ V 1 T সেট করুন, পুনর্নির্মাণের পরে ভর ম্যাট্রিক্স একক ম্যাট্রিক্স হওয়ার জন্য Q 1 = V 1 Λ − 1 / 2 Q_1 = V_1\Lambda^{-1/2} Q 1 = V 1 Λ − 1/2 নির্বাচন করুন।
V ~ = W V \tilde{V} = \sqrt{W}V V ~ = W V এর SVD বিয়োজন করুন: V ~ = U 2 Σ V 2 T \tilde{V} = U_2\Sigma V_2^T V ~ = U 2 Σ V 2 T , Q 2 = V 2 Σ − 1 Q_2 = V_2\Sigma^{-1} Q 2 = V 2 Σ − 1 নির্বাচন করুন।
তাত্ত্বিক বিশ্লেষণ দেখায় যে পদ্ধতি ২ সংখ্যাসূচক বাস্তবায়নে আরও স্থিতিশীল, কারণ এটি শর্ত সংখ্যার বর্গ সহ ম্যাট্রিক্সের SVD বিয়োজন এড়ায়।
ক্ষেত্র Ω = [ − 1 , 1 ] 2 \Omega = [-1,1]^2 Ω = [ − 1 , 1 ] 2 বিবেচনা করুন, সীমানা x 2 + y 2 = r 0 2 x^2 + y^2 = r_0^2 x 2 + y 2 = r 0 2 (r 0 = 1 / 3 r_0 = 1/\sqrt{3} r 0 = 1/ 3 ), পরীক্ষা ফাংশন:
u ( x , y ) = { 1 β + cos ( 2 π r 2 ) , r > r 0 1 β − cos ( 2 π r 2 ) + cos ( 2 π r 0 2 ) ( 1 β + − 1 β − ) , r < r 0 u(x,y) = \begin{cases}
\frac{1}{\beta^+}\cos(2\pi r^2), & r > r_0 \\
\frac{1}{\beta^-}\cos(2\pi r^2) + \cos(2\pi r_0^2)\left(\frac{1}{\beta^+} - \frac{1}{\beta^-}\right), & r < r_0
\end{cases} u ( x , y ) = { β + 1 cos ( 2 π r 2 ) , β − 1 cos ( 2 π r 2 ) + cos ( 2 π r 0 2 ) ( β + 1 − β − 1 ) , r > r 0 r < r 0
L 2 L^2 L 2 প্রজেকশন ত্রুটির সংগ্রহ ক্রমভর ম্যাট্রিক্স শর্ত সংখ্যা ভিত্তি ফাংশন নির্মাণের সংখ্যাসূচক স্থিতিশীলতা বিভিন্ন বহুপদী ক্রম m = 1 , 2 , 3 , 4 m = 1,2,3,4 m = 1 , 2 , 3 , 4 এর জন্য, L 2 L^2 L 2 প্রজেকশন ত্রুটি সর্বোত্তম সংগ্রহ ক্রম প্রদর্শন করে:
m m m N = 16 N = 16 N = 16 N = 32 N = 32 N = 32 N = 64 N = 64 N = 64 N = 128 N = 128 N = 128 N = 256 N = 256 N = 256 সংগ্রহ ক্রম 1 8.14E-02 2.08E-02 5.23E-03 1.31E-03 3.28E-04 ≈2.0 2 9.29E-03 1.19E-03 1.50E-04 1.88E-05 2.36E-06 ≈3.0 3 8.80E-04 5.74E-05 3.62E-06 2.27E-07 1.42E-08 ≈4.0 4 7.55E-05 2.39E-06 7.49E-08 2.34E-09 7.33E-11 ≈5.0
পুনর্নির্মাণের আগে এবং পরে ভর ম্যাট্রিক্স শর্ত সংখ্যার তুলনা দেখায়:
প্রাথমিক নির্মাণের শর্ত সংখ্যা বহুপদী ক্রমের সাথে সূচকীয়ভাবে বৃদ্ধি পায় পদ্ধতি ২ শর্ত সংখ্যা স্থিতিশীলভাবে ১.০ এ হ্রাস করতে পারে পদ্ধতি ১ উচ্চ ক্রমে সংখ্যাসূচক অস্থিতিশীলতা প্রদর্শন করে সহগ ম্যাট্রিক্স A A A এবং A ~ \tilde{A} A ~ এর জন্য, Jacobi পূর্ব-প্রক্রিয়াকারী এবং সারি স্বাভাবিকীকরণ পূর্ব-প্রক্রিয়াকারী উভয়ই শর্ত সংখ্যা উল্লেখযোগ্যভাবে উন্নত করতে পারে, বিশেষত উচ্চ ক্রম এবং ছোট জাল আকারের ক্ষেত্রে।
ঐতিহ্যবাহী IFE পদ্ধতি: শুধুমাত্র সীমানা শর্ত আনুমানিক সন্তুষ্ট করতে পারে সম্প্রসারিত IFE পদ্ধতি: Cauchy সম্প্রসারণের উপর ভিত্তি করে উচ্চ-ক্রম পদ্ধতি এই গবেষণাপত্রের GC-IFE: সীমানা শর্ত সঠিকভাবে সন্তুষ্ট করে এমন জ্যামিতি সামঞ্জস্যপূর্ণ পদ্ধতি ভার্চুয়াল উপাদান পদ্ধতিতে অর্থোগোনাল বহুপদী নির্মাণ DG পদ্ধতিতে অর্থোনরমাল ভিত্তি ফাংশন এই গবেষণাপত্র অর্থোগোনালাইজেশন ধারণা IFE কাঠামোতে প্রবর্তন করে তাত্ত্বিক অবদান : Frenet যন্ত্রের উপর ভিত্তি করে GC-IFE তাত্ত্বিক কাঠামো প্রতিষ্ঠাঅ্যালগরিদম উদ্ভাবন : দুটি নির্মাণ পদ্ধতি এবং দুটি অর্থোগোনালাইজেশন কৌশল প্রস্তাবসংখ্যাসূচক যাচাইকরণ : পদ্ধতির সর্বোত্তম সংগ্রহ এবং সংখ্যাসূচক স্থিতিশীলতা প্রমাণব্যবহারিক মূল্য : সম্পূর্ণ MATLAB বাস্তবায়ন প্রদানজ্যামিতি সীমাবদ্ধতা : সীমানা বক্র যথেষ্ট মসৃণ হওয়ার প্রয়োজন (C 3 C^3 C 3 )গণনামূলক জটিলতা : উচ্চ-ক্রম ক্ষেত্রে ম্যাট্রিক্স সমাবেশ খরচ বেশিসম্প্রসারণযোগ্যতা : বর্তমানে শুধুমাত্র দ্বি-মাত্রিক উপবৃত্তাকার সীমানা সমস্যায় প্রযোজ্যত্রি-মাত্রিক সমস্যা এবং জটিল জ্যামিতিতে সম্প্রসারণ সময়-নির্ভর সীমানা সমস্যায় প্রয়োগ স্থিতিস্থাপক বলবিদ্যা, Stokes সিস্টেম ইত্যাদিতে সম্প্রসারণ আরও দক্ষ বাস্তবায়ন অ্যালগরিদম উন্নয়ন তাত্ত্বিক কঠোরতা : অন্তরক জ্যামিতির Frenet যন্ত্র দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করেপদ্ধতি উদ্ভাবন : প্রথমবারের মতো সীমানা লাফ শর্ত সঠিকভাবে সন্তুষ্ট করে এমন উচ্চ-ক্রম IFE পদ্ধতি বাস্তবায়নসংখ্যাসূচক স্থিতিশীলতা : অর্থোগোনালাইজেশন পুনর্নির্মাণ গণনামূলক স্থিতিশীলতা উল্লেখযোগ্যভাবে উন্নত করেসম্পূর্ণ বাস্তবায়ন : বিস্তারিত অ্যালগরিদম বর্ণনা এবং MATLAB কোড প্রদানপর্যাপ্ত যাচাইকরণ : তাত্ত্বিক বিশ্লেষণ এবং সংখ্যাসূচক পরীক্ষা পরস্পর নিশ্চিত করেসীমিত প্রয়োগযোগ্যতা : বর্তমানে শুধুমাত্র দ্বি-মাত্রিক উপবৃত্তাকার সীমানা সমস্যায় সীমাবদ্ধগণনামূলক খরচ : উচ্চ-ক্রম ক্ষেত্রে গণনামূলক জটিলতা বিশ্লেষণ যথেষ্ট গভীর নয়কঠোর জ্যামিতি প্রয়োজনীয়তা : সীমানা বক্রের মসৃণতার জন্য উচ্চ প্রয়োজনীয়তাব্যবহারিক প্রয়োগ কেস : জটিল প্রকৌশল সমস্যার প্রয়োগ যাচাইকরণের অভাবএকাডেমিক মূল্য : IFE পদ্ধতির জন্য নতুন তাত্ত্বিক কাঠামো এবং প্রযুক্তিগত পথ প্রদান করেব্যবহারিক মূল্য : উচ্চ-নির্ভুলতা সীমানা সমস্যা সমাধানের জন্য কার্যকর সরঞ্জাম প্রদান করেপুনরুৎপাদনযোগ্যতা : সম্পূর্ণ কোড বাস্তবায়ন অন্যান্য গবেষকদের ব্যবহার এবং সম্প্রসারণ সহজতর করেউন্নয়ন সম্ভাবনা : পরবর্তী গবেষণার জন্য দৃঢ় ভিত্তি স্থাপন করেবহু-পর্যায়ের প্রবাহ : বিভিন্ন তরল সীমানার সংখ্যাসূচক অনুকরণউপকরণ বিজ্ঞান : যৌগিক উপকরণ সীমানা সমস্যার উচ্চ-নির্ভুলতা গণনাজৈব চিকিৎসা : টিস্যু সীমানার তাপ এবং ভর স্থানান্তর সমস্যাভূ-পদার্থবিজ্ঞান : ভূস্তর সীমানার তরঙ্গ প্রচার সমস্যাগবেষণাপত্র ২২টি সম্পর্কিত সংদর্ভ উদ্ধৃত করে, যা অন্তরক জ্যামিতি, সংখ্যাসূচক বিশ্লেষণ, সীমান্ত উপাদান পদ্ধতি এবং অন্যান্য ক্ষেত্রের গুরুত্বপূর্ণ কাজ অন্তর্ভুক্ত করে, এই গবেষণার জন্য দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।