طريقة التجزئة بالتدرج

الحل الدقيق u¯(x)=34(0.54/3-|x-0.5|4/3){\displaystyle {\overline {u}}(x)={\frac {3}{4}}\left({0.5}^{4/3}-|x-0.5|^{4/3}\right)} مسألة لابلاس من الرتبة p-(|u¯|2u¯)(x)=1{\displaystyle -(|{\overline {u}}'|^{2}{\overline {u}}')'(x)=1}على المجال [0,1] معu¯(0)=u¯(1)=0{\displaystyle {\overline {u}}(0)={\overline {u}}(1)=0} (الخط الأسود) والتقريبي (الخط الأزرق) المحسوب باستخدام طريقة غاليركين غير المتصلة من الدرجة الأولى الموصولة بـ GDM (شبكة موحدة مع 6 عناصر).

في الرياضيات العددية ، تُعدّ طريقة التقطيع بالتدرج ( GDM ) إطارًا يتضمن مخططات عددية كلاسيكية وحديثة لمسائل الانتشار بأنواعها المختلفة: الخطية وغير الخطية، والمستقرة أو المتغيرة مع الزمن. قد تكون هذه المخططات متوافقة أو غير متوافقة، وقد تعتمد على شبكات مضلعة أو متعددة السطوح عامة جدًا (أو قد تكون بدون شبكة).

تتطلب بعض الخصائص الأساسية إثبات تقارب طريقة التوزيع المعمم (GDM). تُمكّن هذه الخصائص الأساسية من تقديم براهين كاملة لتقارب طريقة التوزيع المعمم للمسائل الإهليلجية والقطع المكافئ، سواء كانت خطية أو غير خطية. بالنسبة للمسائل الخطية، سواء كانت مستقرة أو عابرة، يمكن تحديد تقديرات الخطأ بناءً على ثلاثة مؤشرات خاصة بطريقة التوزيع المعمم [ 1 ] (الكمياتجد{\displaystyle C_{D}}،Sد{\displaystyle S_{D}}ودبليود{\displaystyle W_{D}}( انظر أدناه ). بالنسبة للمسائل غير الخطية، تستند البراهين إلى تقنيات التراص ولا تتطلب أي افتراضات انتظام قوية غير فيزيائية على الحل أو بيانات النموذج. [ 2 ] تشمل النماذج غير الخطية التي تم فيها إجراء برهان تقارب طريقة GDM ما يلي: مسألة ستيفان التي تُنمذج مادة منصهرة، وتدفقات ثنائية الطور في الأوساط المسامية، ومعادلة ريتشاردز لتدفق المياه الجوفية، ومعادلات ليري-ليونز غير الخطية تمامًا. [ 3 ]

من المعروف أن أي مخطط يدخل في إطار عمل GDM يتقارب مع جميع هذه المشكلات. وينطبق هذا بشكل خاص على العناصر المحدودة المتوافقة ، والعناصر المحدودة المختلطة ، والعناصر المحدودة غير المتوافقة ، وفي حالة المخططات الأحدث، طريقة غاليركين غير المتصلة ، وطريقة المحاكاة المختلطة الهجينة، وطريقة الفروق المحدودة للمحاكاة العقدية ، وبعض مخططات الحجم المحدود للازدواجية المنفصلة، ​​وبعض مخططات تقريب التدفق متعدد النقاط.

مثال على مسألة الانتشار الخطي

لنفترض معادلة بواسون في مجال مفتوح محدودΩRد{\displaystyle \Omega \subset \mathbb {R} ^{d}}، مع شرط حدودي متجانس من نوع ديريشليه

أينول2(Ω){\displaystyle f\in L^{2}(\Omega )}المعنى المعتاد للحل الضعيف [ 4 ] لهذا النموذج هو:

باختصار، تتألف طريقة التقطيع التدرجي (GDM) لمثل هذا النموذج من اختيار فضاء محدود الأبعاد ومؤثرين لإعادة البناء (أحدهما للدوال والآخر للتدرجات) واستبدال هذه العناصر المنفصلة بالعناصر المستمرة في المعادلة (2). وبشكل أدق، تبدأ طريقة التقطيع التدرجي بتعريف التقطيع التدرجي (GD)، وهو عبارة عن ثلاثية.د=(Xد،0،Πد،د){\displaystyle D=(X_{D,0},\Pi _{D},\nabla _{D})}، أين:

  • مجموعة المجاهيل المنفصلةXد،0{\displaystyle X_{D,0}}هو فضاء متجهي حقيقي ذو أبعاد محدودة،
  • إعادة بناء الوظيفةΠد : Xد،0ل2(Ω){\displaystyle \Pi _{D}~:~X_{D,0}\to L^{2}(\Omega )}هي عملية تحويل خطية تعيد بناء عنصر منXد،0{\displaystyle X_{D,0}}، دالة علىΩأوميغا،
  • إعادة بناء التدرجد : Xد،0ل2(Ω)د{\displaystyle \nabla _{D}~:~X_{D,0}\to L^{2}(\Omega )^{d}}هي عملية تحويل خطية تعيد بناء عنصر منXد،0{\displaystyle X_{D,0}}، "التدرج" (دالة ذات قيم متجهة) علىΩأوميغايجب اختيار إعادة بناء التدرج هذا بحيثدل2(Ω)د{\displaystyle \Vert \nabla _{D}\cdot \Vert _{L^{2}(\Omega )^{d}}}هو أمر طبيعي فيXد،0{\displaystyle X_{D,0}}.

يُعطى مخطط التدرج ذي الصلة لتقريب (2) كما يلي: أوجدuXد،0{\displaystyle u\in X_{D,0}}بحيث

في هذه الحالة، تُعدّ طريقة GDM طريقةً غير متوافقة لتقريب المعادلة (2)، والتي تتضمن طريقة العناصر المحدودة غير المتوافقة. تجدر الإشارة إلى أن العكس غير صحيح، بمعنى أن إطار عمل GDM يتضمن طرقًا بحيث تكون الدالةدu{\displaystyle \nabla _{D}u}لا يمكن حسابها من الدالةΠدu{\displaystyle \Pi _{D}u}.

ينطبق تقدير الخطأ التالي، المستوحى من اللمة الثانية لـ ج. سترانج، [ 5 ]

و

تعريف:

والذي يقيس الإكراه (ثابت بوانكاريه المنفصل)،

والذي يقيس خطأ الاستيفاء،

والذي يقيس عيب المطابقة.

لاحظ أنه يمكن اشتقاق الحدود العليا والسفلى التالية لخطأ التقريب:

ثم إن الخصائص الأساسية اللازمة والكافية لتقارب الطريقة، بالنسبة لمجموعة من طرق التدرج المعمم، هي خصائص الإكراه، والاتساق مع التدرج المعمم، والتوافق الحدي، كما هو موضح في القسم التالي. وبشكل أعم، تكفي هذه الخصائص الأساسية الثلاث لإثبات تقارب طريقة التدرج المعمم للمسائل الخطية ولبعض المسائل غير الخطية مثل...ص{\displaystyle p}مسألة لابلاس. بالنسبة للمسائل غير الخطية مثل الانتشار غير الخطي، والمسائل القطعية المتدهورة...، نضيف في القسم التالي خاصيتين أساسيتين أخريين قد تكونان مطلوبتين.

الخصائص الأساسية التي تسمح بتقارب نموذج GDM

يترك(دم)مشمال{\displaystyle (D_{m})_{m\in \mathbb {N} }}أن تكون عائلة من GDs، معرفة كما سبق (مرتبطة بشكل عام بتسلسل من الشبكات المنتظمة التي يميل حجمها إلى 0).

الإكراه

التسلسل(جدم)مشمال{\displaystyle (C_{D_{m}})_{m\in \mathbb {N} }}(المحدد بواسطة ( 6 )) يبقى محدودًا.

اتساق GD

للجميعφح01(Ω){\displaystyle \varphi \in H_{0}^{1}(\Omega )}،ليممSدم(φ)=0{\displaystyle \lim _{m\to \infty }S_{D_{m}}(\varphi )=0}(محدد بواسطة ( 7 )).

المطابقة الحدية

للجميعφحdiv(Ω){\displaystyle \varphi \in H_{\operatorname {div} }(\Omega )}،ليممدبليودم(φ)=0{\displaystyle \lim _{m\to \infty }W_{D_{m}}(\varphi )=0}(المحددة بواسطة ( 8 )). هذه الخاصية تستلزم خاصية الإكراه.

التراص (مطلوب لبعض المسائل غير الخطية)

لكل التسلسل(uم)مشمال{\displaystyle (u_{m})_{m\in \mathbb {N} }}بحيثuمXدم،0{\displaystyle u_{m}\in X_{D_{m},0}}للجميع مشمال{\displaystyle m\in \mathbb {N} }و (uمدم)مشمال{\displaystyle (\Vert u_{m}\Vert _{D_{m}})_{m\in \mathbb {N} }}إذا كانت محدودة، فإن المتتالية(Πدمuم)مشمال{\displaystyle (\Pi _{D_{m}}u_{m})_{m\in \mathbb {N} }}يتميز بصغر حجمه نسبياًل2(Ω){\displaystyle L^{2}(\Omega )}(هذه الخاصية تستلزم خاصية الإكراه).

إعادة بناء ثابتة جزئياً (مطلوبة لبعض المسائل غير الخطية)

يترك د=(Xد،0،Πد،د){\displaystyle D=(X_{D,0},\Pi _{D},\nabla _{D})}يكون تجزئة تدرجية كما هو مُعرَّف أعلاه. المؤثرΠد{\displaystyle \Pi _{D}}تكون إعادة البناء ثابتة جزئياً إذا وُجد أساس(هـأنا)أناب{\displaystyle (e_{i})_{i\in B}}لXد،0{\displaystyle X_{D,0}}وعائلة من المجموعات الفرعية المنفصلة (Ωأنا)أناب{\displaystyle (\Omega _{i})_{i\in B}}لΩ{\displaystyle \Omega }بحيث Πدu=أنابuأناχΩأنا{\textstyle \Pi _{D}u=\sum _{i\in B}u_{i}\chi _{\Omega _{i}}}للجميعu=أنابuأناهـأناXد،0{\textstyle u=\sum _{i\in B}u_{i}e_{i}\in X_{D,0}}، أينχΩأنا{\displaystyle \chi _{\Omega _{i}}}هي الوظيفة المميزة لـΩأنا{\displaystyle \Omega _{i}}.

بعض المسائل غير الخطية مع براهين تقارب كاملة لنموذج GDM

نستعرض بعض المشاكل التي يمكن إثبات تقارب طريقة GDM فيها عندما يتم استيفاء الخصائص الأساسية المذكورة أعلاه.

مشاكل الانتشار الثابت غير الخطي

-div(Λ(u¯)u¯)=و{\displaystyle -\operatorname {div} (\Lambda ({\overline {u}})\nabla {\overline {u}})=f}

في هذه الحالة، يتقارب نموذج GDM في ظل خصائص الإكراه، والاتساق GD، والتوافق الحدي، والتراص.

مسألة لابلاس من الرتبة p عندما تكون p > 1

-div(|u¯|ص-2u¯)=و{\displaystyle -\operatorname {div} \left(|\nabla {\overline {u}}|^{p-2}\nabla {\overline {u}}\right)=f}

في هذه الحالة، يجب كتابة الخصائص الأساسية، واستبدالهال2(Ω){\displaystyle L^{2}(\Omega )} بواسطةلص(Ω){\displaystyle L^{p}(\Omega )}،ح01(Ω){\displaystyle H_{0}^{1}(\Omega )}بواسطةدبليو01،ص(Ω){\displaystyle W_{0}^{1,p}(\Omega )} وحdiv(Ω){\displaystyle H_{\operatorname {div} }(\Omega )}بواسطةدبليوdivص(Ω){\displaystyle W_{\operatorname {div} }^{p'}(\Omega )}مع1ص+1ص=1{\textstyle {\frac {1}{p}}+{\frac {1}{p'}}=1}، ولا يتقارب نموذج GDM إلا في ظل خصائص الإكراه، والاتساق مع GD، والتوافق الحدي.

معادلة الحرارة الخطية وغير الخطية

تu¯-div(Λ(u¯)u¯)=و{\displaystyle \partial _{t}{\overline {u}}-\operatorname {div} (\Lambda ({\overline {u}})\nabla {\overline {u}})=f}

في هذه الحالة، يتقارب نموذج GDM في ظل خصائص الإكراه، والاتساق GD (المكيف مع مشاكل المكان والزمان)، والتوافق الحدي، والتراص (للحالة غير الخطية).

مسائل القطع المكافئ المتدهورة

افترض أنβ{\displaystyle \beta }وζ{\displaystyle \zeta }هي دوال متصلة غير متناقصة وفقًا لشرط ليبشيتز:

تβ(u¯)-Δζ(u¯)=و{\displaystyle \partial _{t}\beta ({\overline {u}})-\Delta \zeta ({\overline {u}})=f}

لاحظ أنه بالنسبة لهذه المشكلة، فإن خاصية إعادة البناء الثابتة جزئيًا مطلوبة، بالإضافة إلى خصائص الإكراه، والاتساق GD (المكيف مع مشاكل المكان والزمان)، والتوافق الحدي، والتراص.

مراجعة لبعض الطرق العددية التي تُعدّ من طرق إدارة البيانات العامة

جميع الطرق المذكورة أدناه تفي بالخصائص الأساسية الأربع الأولى لـ GDM (الإكراه، والاتساق GD، والتوافق الحدي، والتراص)، وفي بعض الحالات بالخاصية الخامسة (إعادة البناء الثابتة جزئيًا).

طرق غاليركين وطرق العناصر المحدودة المطابقة

يتركVحح01(Ω){\displaystyle V_{h}\subset H_{0}^{1}(\Omega )}يتم توليدها بواسطة الأساس المحدود(ψأنا)أناأنا{\displaystyle (\psi _{i})_{i\in I}}طريقة غاليركين فيVح{\displaystyle V_{h}}وهو مطابق لنموذج إدارة البيانات العامة (GDM) حيث يتم تعريف

  • Xد،0={u=(uأنا)أناأنا}=Rأنا،{\displaystyle X_{D,0}=\{u=(u_{i})_{i\in I}\}=\mathbb {R} ^{I},}
  • Πدu=أناأناuأناψأنا{\displaystyle \Pi _{D}u=\sum _{i\in I}u_{i}\psi _{i}}
  • دu=أناأناuأناψأنا.{\displaystyle \nabla _{D}u=\sum _{i\in I}u_{i}\nabla \psi _{i}.}

في هذه الحالة،جد{\displaystyle C_{D}}هو الثابت المُتضمن في متباينة بوانكاريه المستمرة، و، بالنسبة لجميعφحdiv(Ω){\displaystyle \varphi \in H_{\operatorname {div} }(\Omega )}،دبليود(φ)=0{\displaystyle W_{D}(\varphi )=0}(المحددة بواسطة ( 8 )). ثم يتم استنتاج ( 4 ) و( 5 ) من خلال ليمّة سيا .

"المُجمّع"P1{\displaystyle P^{1}}تدخل حالة العناصر المحدودة في إطار عمل GDM، لتحل محلΠدu{\displaystyle \Pi _{D}u}بواسطة Π~دu=أناأناuأناχΩأنا{\textstyle {\widetilde {\Pi }}_{D}u=\sum _{i\in I}u_{i}\chi _{\Omega _{i}}}، أينΩأنا{\displaystyle \Omega _{i}}هي خلية مزدوجة متمركزة على الرأس المفهرس بواسطة أناأنا{\displaystyle i\in I}. استخدام التجميع الكتلي يسمح بالحصول على خاصية إعادة البناء الثابتة جزئيًا.

عنصر محدود غير متوافق

على شبكةتي{\displaystyle T}وهي مجموعة متوافقة من الأشكال البسيطة لـRد{\displaystyle \mathbb {R} ^{d}}غير المتوافقP1{\displaystyle P^{1}}تُعرَّف العناصر المحدودة بواسطة الأساس (ψأنا)أناأنا{\displaystyle (\psi _{i})_{i\in I}}من الدوال التي تكون خطية في أي كتي{\displaystyle K\in T}، وقيمتها عند مركز ثقل وجه معين من الشبكة تساوي 1 و0 عند جميع الوجوه الأخرى (تُستخدم هذه العناصر المحدودة في [Crouzeix et al ] [ 6 ] لتقريب معادلات ستوكس ونافير-ستوكس ). عندئذٍ، تدخل الطريقة في إطار GDM بنفس التعريف كما في حالة طريقة غاليركين، باستثناء حقيقة أن ψأنا{\displaystyle \nabla \psi _{i}}يجب فهمها على أنها "التدرج المكسور" لـψأنا{\displaystyle \psi _{i}}، بمعنى أنها دالة ثابتة جزئياً تساوي في كل سيمبلكس تدرج الدالة الأفينية في السيمبلكس.

عنصر محدود مختلط

تعتمد طريقة العناصر المحدودة المختلطة على تعريف فضاءين منفصلين، أحدهما لتقريبu¯{\displaystyle \nabla {\overline {u}}}وواحد آخر لـu¯{\displaystyle {\overline {u}}}[ 7 ] يكفي استخدام العلاقات المنفصلة بين هذه التقريبات لتعريف نموذج GDM. يسمح استخدام دوال أساس رافيارت-توماس منخفضة الدرجة بالحصول على خاصية إعادة البناء الثابتة جزئيًا.

طريقة غاليركين غير المتصلة

تعتمد طريقة غاليركين غير المتصلة على تقريب المسائل بدالة متعددة الحدود مجزأة، دون اشتراطات على الانتقالات من عنصر إلى آخر. [ 8 ] يتم دمجها في إطار عمل GDM عن طريق تضمين حد قفز في التدرج المنفصل، والذي يعمل كتنظيم للتدرج بمعنى التوزيع.

طريقة الفروق المحدودة المحاكية وطريقة الفروق المحدودة المحاكية العقدية

تم تقديم هذه المجموعة من الطرق بواسطة [Brezzi et al ] [ 9 ] واكتملت في [Lipnikov et al ] [ 10 ] . وهي تتيح تقريب المسائل الإهليلجية باستخدام فئة واسعة من الشبكات متعددة السطوح. وقد تم إثبات اندماجها في إطار GDM في [Droniou et al ] [ 2 ] .

انظر أيضاً

مراجع

  1. ر. إيمارد، س. غيشارد، و ر. هيربين . مخططات ثلاثية الأبعاد باستخدام قوالب صغيرة لتدفقات الانتشار في الأوساط المسامية. M2AN، 46:265–290، 2012.
  2. 1 2 ج. درونيو، ر. إيمارد، ت. غالويت، و ر. هيربين . مخططات التدرج: إطار عام لتقسيم المعادلات الإهليلجية والقطع المكافئ الخطية وغير الخطية وغير المحلية. نماذج رياضية وطرق العلوم التطبيقية (M3AS)، 23(13):2395–2432، 2013.
  3. ^ ج. ليري و ج. ليونز. بعض نتائج Višik على المشكلات البيضاوية غير الخطية باستخدام طرق Minty-Browder. ثور. شركة نفط الجنوب. الرياضيات. فرنسا، 93: 97-107، 1965.
  4. هـ. بريزيس . التحليل الوظيفي، فضاءات سوبوليف والمعادلات التفاضلية الجزئية. سلسلة Universitext. سبرينغر، نيويورك، 2011.
  5. ج. سترانج. الجرائم التباينية في طريقة العناصر المحدودة. في الأسس الرياضية لطريقة العناصر المحدودة مع تطبيقات على المعادلات التفاضلية الجزئية (وقائع الندوة، جامعة ماريلاند، بالتيمور، ماريلاند، 1972) ، الصفحات 689-710. دار النشر الأكاديمية، نيويورك، 1972.
  6. م. كروزيكس وب. أ. رافيارت. طرق العناصر المحدودة المطابقة وغير المطابقة لحل معادلات ستوكس الثابتة. الجزء الأول. المجلة الفرنسية للأتمتة والمعلوماتية وبحوث العمليات، السلسلة الحمراء، 7(R-3):33–75، 1973.
  7. بي. أ. رافيارت وجيه إم توماس. طريقة العناصر المحدودة المختلطة لمسائل القطع الناقص من الرتبة الثانية. في الجوانب الرياضية لطرق العناصر المحدودة (وقائع المؤتمر، المجلس الوطني للبحوث (CNR)، روما، 1975) ، الصفحات 292-315. سلسلة محاضرات في الرياضيات، المجلد 606. سبرينغر، برلين، 1977.
  8. دي بيترو وإرن. الجوانب الرياضية لطرق غاليركين غير المتصلة، المجلد 69 من سلسلة الرياضيات والتطبيقات (برلين). سبرينغر، هايدلبرغ، 2012.
  9. ف. بريتزي، ك. ليبنيكوف، و م. شاشكوف. تقارب طريقة الفروق المحدودة المحاكية لمسائل الانتشار على الشبكات متعددة السطوح. مجلة SIAM للتحليل العددي، 43(5):1872–1896، 2005.
  10. ك. ليبنيكوف، ج. مانزيني، و م. شاشكوف. طريقة الفروق المحدودة المحاكية. مجلة الفيزياء الحاسوبية، 257-الجزء ب: 1163-1227، 2014.