2025-11-23T03:55:16.576493

A finite element method using a bounded auxiliary variable for solving the Richards equation

Benfanich, Bourgault, Beljadid
The Richards equation, a nonlinear elliptic parabolic equation, is widely used to model infiltration in porous media. We develop a finite element method for solving the Richards equation by introducing a new bounded auxiliary variable to eliminate unbounded terms in the weak formulation of the method. This formulation is discretized using a semi-implicit scheme and the resulting nonlinear system is solved using Newton's method. Our approach eliminates the need of regularization techniques and offers advantages in handling both dry and fully saturated zones. In the proposed techniques, a non-overlapping Schwarz domain decomposition method is used for modeling infiltration in layered soils. We apply the proposed method to solve the Richards equation using the Havercamp and van Genuchten models for the capillary pressure. Numerical experiments are performed to validate the proposed approach, including tests such as modeling flows in fibrous sheets where the initial medium is totally dry, two cases with fully saturated and dry regions, and an infiltration problem in layered soils. The numerical results demonstrate the stability and accuracy of the proposed numerical method. The numerical solutions remain positive in the presence of totally dry zones. The numerical investigations clearly demonstrated the capability of the proposed method to effectively predict the dynamics of flows in unsaturated soils.
academic

طريقة العناصر المحدودة باستخدام متغير مساعد محدود لحل معادلة ريتشاردز

المعلومات الأساسية

  • معرّف الورقة: 2510.13012
  • العنوان: طريقة العناصر المحدودة باستخدام متغير مساعد محدود لحل معادلة ريتشاردز
  • المؤلفون: Abderrahmane Benfanich (جامعة أوتاوا)، Yves Bourgault (جامعة أوتاوا)، Abdelaziz Beljadid (جامعة أوتاوا وجامعة محمد السادس متعددة التخصصات)
  • التصنيف: math.NA cs.NA physics.comp-ph
  • تاريخ النشر: 14 أكتوبر 2025 (نسخة أولية على arXiv)
  • رابط الورقة: https://arxiv.org/abs/2510.13012

الملخص

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

الخلفية البحثية والدافع

تعريف المشكلة

تصف معادلة ريتشاردز التدفق غير المشبع في الوسائط المسامية وهي نموذج أساسي في الهيدروجيولوجيا والهندسة البيئية وعلوم التربة. تتخذ المعادلة عدة أشكال:

  1. صيغة رأس الضغط (Ψ): مناسبة للتدفق المشبع بالكامل والجزئي غير المشبع، لكنها غير مستقرة عددياً عندما S→0
  2. صيغة التشبع (S): تؤدي أداءً أفضل في التربة الجافة، لكنها تواجه صعوبات عددية في المناطق المشبعة (S=1)
  3. الصيغة المختلطة (S,Ψ): تحسن خصائص حفظ الكتلة، لكنها تزيد التعقيد

قيود الطرق الموجودة

  1. طرق تبديل المتغيرات الرئيسية: قد تنتج حلولاً غير فيزيائية عند واجهات التشبع/عدم التشبع
  2. تقنيات التنظيم: تتطلب اختيار معاملات نمذجة إضافية
  3. طرق العناصر المحدودة التقليدية: تواجه صعوبات عددية في التعامل مع المناطق الجافة تماماً (S=0) والمشبعة بالكامل (S=1)

الدافع البحثي

تطوير طريقة عددية موحدة قادرة على:

  • تجنب تبديل المتغيرات والتنظيم
  • معالجة المناطق الجافة والمشبعة بشكل مستقر
  • الحفاظ على خصائص حفظ الكتلة الجيدة
  • التطبيق على الوسائط المسامية غير المتجانسة

المساهمات الأساسية

  1. اقتراح صيغة u جديدة: إدخال متغير مساعد محدود u لحذف الحدود غير المحدودة في الصيغة الضعيفة
  2. تطوير صيغة عددية مستقرة: دمج التقسيم الزمني شبه الضمني وتكرار نيوتن
  3. تنفيذ طريقة تحليل المجال: معالجة مشاكل الوسائط غير المتجانسة مثل التربة الطبقية
  4. التحقق من فعالية الطريقة: إثبات استقرار وتدقيق الطريقة من خلال تجارب عددية متعددة

شرح الطريقة

تعريف المهمة

حل معادلة ريتشاردز: ∂θ∂t+∇⋅q=0\frac{\partial\theta}{\partial t} + \nabla \cdot q = 0

حيث θ هي محتوى الماء، و q هي تدفق المياه وفقاً لقانون دارسي-باكينغهام: q=−Ks(x)Kr(x,S)∇(Ψ(x,S)+z)q = -K_s(x)K_r(x,S)\nabla(\Psi(x,S) + z)

بناء المتغير الجديد u

بالنسبة لجميع النماذج، يمكن كتابة مشتقة دالة ليفيريت على النحو التالي: J′(S)=C(x)S−a(1−Sc)−bJ'(S) = C(x)S^{-a}(1-S^c)^{-b}

تعريف المتغير المساعد: u=∫0S(1−sc)−bdsu = \int_0^S (1-s^c)^{-b} ds

التعبير باستخدام دالة بيتا غير الكاملة: u=1cB(Sc,1c,1−b)u = \frac{1}{c}B(S^c, \frac{1}{c}, 1-b)

صيغة u لمعادلة ريتشاردز

في الوسط المتجانس، صيغة u هي: ϕ∂S(u)∂t−∇⋅(KsKr(S(u))(hcapCS(u)−a∇u+ez))=0\phi \frac{\partial S(u)}{\partial t} - \nabla \cdot \left(K_sK_r(S(u))\left(h_{cap}CS(u)^{-a}\nabla u + e_z\right)\right) = 0

الصيغة المتغيرة

البحث عن u ∈ L²(I;H¹(Ω)) بحيث: ∫Ωϕ∂S(u)∂tvdx+∫ΩKsKr(S(u))(hcapS(u)−a∇u+ez)⋅∇vdx+حدود الحدود=0\int_\Omega \phi \frac{\partial S(u)}{\partial t} v dx + \int_\Omega K_sK_r(S(u))\left(h_{cap}S(u)^{-a}\nabla u + e_z\right) \cdot \nabla v dx + \text{حدود الحدود} = 0

التقسيم الزمني

استخدام صيغة أويلر شبه الضمنية: ∫ΩϕS(uhn+1)−S(uhn)Δtvhdx+∫ΩKsKr(S(uhn))(hcapS(uhn)−a∇uhn+1+ez)⋅∇vhdx=0\int_\Omega \phi \frac{S(u_h^{n+1}) - S(u_h^n)}{\Delta t} v_h dx + \int_\Omega K_sK_r(S(u_h^n))\left(h_{cap}S(u_h^n)^{-a}\nabla u_h^{n+1} + e_z\right) \cdot \nabla v_h dx = 0

تكرار نيوتن

بالنسبة للحدود غير الخطية S(u^{n+1,m+1})، استخدام تطور تايلور من الدرجة الأولى: S(un+1,m+1)=S(un+1,m)+∂S∂u(un+1,m)(un+1,m+1−un+1,m)+O((un+1,m+1−un+1,m)2)S(u^{n+1,m+1}) = S(u^{n+1,m}) + \frac{\partial S}{\partial u}(u^{n+1,m})(u^{n+1,m+1} - u^{n+1,m}) + O((u^{n+1,m+1} - u^{n+1,m})^2)

طريقة تحليل المجال

بالنسبة للتربة الطبقية، استخدام طريقة شوارتز غير المتداخلة:

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

إعداد التجارب

التنفيذ الرقمي

  • منصة البرمجيات: FreeFEM++ و FENICSx
  • حلالات خطية: UMFPACK و PETSc و MUMPS
  • الشبكات: شبكات مثلثة موحدة، عناصر محدودة خطية (k=1)

مؤشرات التقييم

خطأ L²: Eh,Δt(T)2=∥u(⋅,T)−uh,Δt(⋅,T)∥L2(Ω)2E_{h,\Delta t}(T)^2 = \|u(\cdot,T) - u_{h,\Delta t}(\cdot,T)\|_{L^2(\Omega)}^2

ترتيب التقارب:

  • ترتيب التقارب المكاني: p=log⁡2(Eh,Δt(T)/Eh/2,Δt(T))p = \log_2(E_{h,\Delta t}(T)/E_{h/2,\Delta t}(T))
  • ترتيب التقارب الزمني: q=log⁡2(Eh,Δt(T)/Eh,Δt/2(T))q = \log_2(E_{h,\Delta t}(T)/E_{h,\Delta t/2}(T))

حالات الاختبار

  1. شريط ألياف 1D غير مشبع بالكامل: التحقق من الاستقرار في ظروف الجفاف
  2. منطقة مختلطة مشبعة-غير مشبعة 1D: اختبار القدرة على التعامل مع القيم الابتدائية غير المتصلة
  3. التحقق من الحل المصنع: حساب ترتيب التقارب
  4. تدفق مشبع بالكامل 2D: المقارنة مع برنامج Hydrus
  5. مشكلة التربة الطبقية: التحقق من طريقة تحليل المجال
  6. وسط غير متجانس يحتوي على أجسام غريبة: اختبار معالجة الأشكال الهندسية المعقدة

نتائج التجارب

النتائج الرئيسية

اختبار شريط الألياف 1D

  • الشبكة: 500 نقطة شبكة موحدة
  • خطوة زمنية: Δt = 1s
  • النتيجة: تطابق كامل بين نتائج صيغة u وصيغة S، توافق جيد مع البيانات التجريبية

منطقة مختلطة مشبعة-غير مشبعة

  • الشبكة: 5000 نقطة شبكة موحدة
  • خطوة زمنية: Δt = 10⁻⁵ ساعة
  • تقارب نيوتن: متوسط 1-5 تكرارات
  • الاستقرار: يحافظ الحل على الرتابة، بدون تذبذب عددي

تحليل تقارب الحل المصنع

التقارب المكاني (Δt = 5×10⁻⁵ ثابت):

hخطأ L²ترتيب التقارب pوقت الحساب (s)
1.43×10⁻²6.04×10⁻³-25.1
7.14×10⁻³8.99×10⁻⁴2.7547.2
3.57×10⁻³1.39×10⁻⁴2.6982.8
1.79×10⁻³3.48×10⁻⁵2.00151.0

التقارب الزمني (h = 1×10⁻³ ثابت):

Δtخطأ L²ترتيب التقارب qوقت الحساب (s)
1×10⁻²8.23×10⁻³-6.9
5×10⁻³3.90×10⁻³1.0811.1
2.5×10⁻³1.92×10⁻³1.0219.3
1.25×10⁻³9.60×10⁻⁴1.0033.3

مقارنة التدفق المشبع 2D

  • الشبكة: 95,142 عنصر مثلثي، 47,972 رأس
  • خطوة زمنية: Δt = 0.001 ساعة
  • النتيجة: توافق عالي مع نتائج برنامج Hydrus

اختبار التربة الطبقية

  • الشبكة: 89,306 عنصر مثلثي، 45,064 رأس
  • خطوة زمنية: Δt = 0.05 ساعة
  • التقارب: تفاوت 10⁻³، بحد أقصى 10 تكرارات، λ = 25
  • النتيجة: متطابقة مع النتائج في المرجع 20

مقارنة الطرق

في اختبار الوسط غير المتجانس يحتوي على أجسام غريبة:

  • صيغة u: مستقرة ومتقاربة حتى عندما ε=0، متوسط 4 تكرارات
  • صيغة Ψ: تتطلب أكثر من 1000 تكرار عندما ε=10⁻⁵، تفاوت البواقي لا يزال 0.04

الأعمال ذات الصلة

الاتجاهات البحثية الرئيسية

  1. طرق تبديل المتغيرات: Diersch & Perrochet (1999)، Forsyth et al. (1995)
  2. تقنيات التنظيم: Schweizer (2007)، Pop & Schweizer (2011)
  3. الصيغ المختلطة: Celia et al. (1990)، Maina & Ackerer (2017)
  4. تحليل المجال: Lions (1990)، Manzini & Ferraris (2004)

مزايا هذه الورقة

  1. المعالجة الموحدة: بدون الحاجة إلى تبديل المتغيرات أو التنظيم
  2. الاستقرار العددي: مستقر في النطاق الكامل S∈0,1
  3. التوافق الفيزيائي: يحافظ على استمرارية رأس الضغط
  4. الكفاءة الحسابية: تقارب نيوتن سريع

الاستنتاجات والمناقشة

الاستنتاجات الرئيسية

  1. فعالية صيغة u: نجحت في حذف الحدود غير المحدودة في الصيغ التقليدية
  2. الاستقرار العددي: يحافظ على الاستقرار في المناطق الجافة والمشبعة بالكامل
  3. أداء التقارب: تقارب من الدرجة الثانية مكانياً، من الدرجة الأولى زمنياً
  4. الجدوى العملية: قابلة للتطبيق على مشاكل الوسائط غير المتجانسة المعقدة

القيود

  1. الاعتماد على النموذج: تتطلب تحقيق الشرط lim⁡S→0K(S)S−a<∞\lim_{S\to 0} K(S)S^{-a} < \infty
  2. حل غير خطي: تتطلب تكرار نيوتن في المناطق المشبعة
  3. معالجة الواجهات: تتطلب تقنيات تحليل المجال للوسائط غير المتجانسة
  4. حساسية المعاملات: معامل روبن λ في تحليل المجال يتطلب ضبطاً

الاتجاهات المستقبلية

  1. التحليل النظري: إثبات صارم للتقارب والاستقرار العددي
  2. الفيزياء المتعددة المقترنة: التوسع إلى نقل المذاب والحرارة
  3. التقنيات التكيفية: تطوير شبكات تكيفية وخطوات زمنية
  4. الحوسبة المتوازية: تنفيذ متوازي على نطاق واسع

التقييم المتعمق

المزايا

  1. الابتكار القوي: بناء متغير u ذكي، أساس نظري متين
  2. الاستقرار العددي: تجنب الصعوبات العددية للطرق التقليدية
  3. التجارب الشاملة: تغطي من 1D إلى 2D، من المتجانس إلى غير المتجانس
  4. الجدوى الهندسية: المقارنة مع برنامج Hydrus التجاري تعزز المصداقية

أوجه القصور

  1. نقص التحليل النظري: غياب إثبات صارم للتقارب والاستقرار
  2. تحليل التعقيد الحسابي: تحليل تكلفة حساب تكرار نيوتن غير كافٍ
  3. اختيار المعاملات: اختيار معامل روبن λ يفتقر إلى التوجيه النظري
  4. التحقق ثلاثي الأبعاد: غياب التحقق من مشاكل ثلاثية الأبعاد

التأثير

  1. القيمة الأكاديمية: توفر منظوراً جديداً لحل معادلة ريتشاردز عددياً
  2. القيمة العملية: آفاق تطبيق واسعة في الهيدروجيولوجيا وعلوم التربة
  3. قابلية التكرار: توفير تفاصيل تنفيذ مفصلة تسهل التكرار

السيناريوهات المعمول بها

  1. تدفق الوسائط المسامية: تدفق المياه الجوفية، حركة الماء في التربة
  2. التطبيقات الهندسية: تسرب السدود، تحليل استقرار المنحدرات
  3. المشاكل البيئية: هجرة الملوثات، استصلاح المياه الجوفية
  4. علوم المواد: تسرب السوائل في مواد الألياف

المراجع

المراجع الرئيسية

  1. Richards, L.A. (1931). Capillary conduction of liquids through porous mediums. Physics 1(5), 318-333.
  2. Celia, M.A., et al. (1990). A general mass-conservative numerical solution for the unsaturated flow equation. Water Resour. Res. 26(7), 1483-1496.
  3. van Genuchten, M.T. (1980). A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 44(5), 892-898.
  4. Manzini, G., Ferraris, S. (2004). Mass-conservative finite volume methods on 2-d unstructured grids for the Richards' equation. Adv. Water Resour. 27, 1199-1215.

التقييم الشامل: هذه ورقة ذات أهمية ابتكارية مهمة في مجال الحل العددي لمعادلة ريتشاردز. من خلال إدخال متغير مساعد محدود u، تحل الورقة بذكاء الصعوبات العددية للطرق التقليدية في التعامل مع المناطق الجافة والمشبعة. تتمتع الطريقة بأساس نظري جيد وقيمة عملية، وتتحقق التجارب الرقمية الشاملة من فعالية الطريقة. على الرغم من أن التحليل النظري يحتاج إلى تعزيز، إلا أن العمل بشكل عام ذو جودة عالية.