A Positivity-Preserving Expectation Scheme for Hamilton--Jacobi--Bellman Equations with Oblique Robin Boundary Conditions

ملخص تقني: مخطط توقع يحافظ على الإيجابية لمعادلات هاملتون-جاكوبي-بلمان (HJB) مع شروط روبن المائلة (Oblique Robin Boundary Conditions)

إعداد المسألة
تتناول الورقة الحل العددي لمعادلات هاملتون-جاكوبي-بلمان (HJB) الخاضعة لشروط روبن المائلة. تمت صياغة المسألة في شكل عكسي على نطاق $O_T = [0, T) \times O$:
$$\begin{cases} -\partial_t u + H(t, x, Du, D^2u, u) = 0, & (t, x) \in O_T, \\ L(t, x, Du, u) = 0, & (t, x) \in [0, T) \times \partial O, \\ u(T, x) = \Psi(x), & x \in O, \end{cases}$$
حيث تتضمن الهاملتونيان $H$ انتشاراً متباين الخواص أو متدهوراً (يتحكم فيه $\sigma$) ودفعاً ($\mu$)، ويتضمن مؤثر الحدود $L$ حداً مشتقاً مائلاً مقترناً بحد من الدرجة صفر $k(t,x,b)u$ (حد روبن) ومصدراً $g$. غالباً ما تفشل طرق الفروق المحدودة أو العناصر المحدودة القياسية في الحفاظ على عدم السالبية (الرتابة) للانتشار متباين الخواص ما لم يتم استيفاء شروط هيمنة قطرية صارمة أو شروط شبكية محددة، وهو أمر إشكالي في الحالات المتدهورة.

المنهجية
يقترح المؤلفون مخططاً رقمياً كاملاً يعتمد على إطار التوقع الشرطي المستمد من تمثيل فاينمان-كاك (Feynman–Kac)، بدلاً من التفكيك الجبري للقوالب (stencils).

1. العقد الداخلية (هوية فاينمان-كاك المنعكسة):

   • يتم تمثيل عملية الانتشار باستخدام تقريب "ويك-أويلر" (weak-Euler) ذي خطوة واحدة مدفوع بمتجه رادماخر (Rademacher vector) ذي أبعاد $m$، مما يولد $P=2^m$ فرعاً ذا احتمالية متساوية.
   • الانعكاس المرآتي: إذا وقعت نقطة تجريبية $X^*$ خارج النطاق، يتم عكسها حول إسقاطها المائل $p_{\gamma_b}(X^*)$ على الحدود. النقطة المنعكسة هي $\tilde{X} = X^* - 2d_{\gamma_b}(X^*)\gamma_b$.
   • الزمن المحلي للحدود: تعمل مسافة الذهاب والإياب $e_d = 2d_{\gamma_b}(X^*)$ كزيادة منفصلة للزمن المحلي للحدود.
   • التوهين (Attenuation): يدخل معامل روبن $k$ في المخطط حصرياً من خلال عوامل التوهين $e^{-k e_d}$ (لقيمة الاستمرار) و $e_d e^{-k e_d/2} g$ (لمصدر الحدود). ومن الأهمية بمكان أن احتمالات الفروع تظل بالضبط $1/P$، مما يحافظ على إيجابية الأوزان دون الحاجة إلى شروط الهيمنة القطرية.
   • الاستيفاء (Interpolation): يتم استرداد القيم عند النقاط المنعكسة باستخدام استيفاء $P1$ (خطي قطعة قطعة) على الشبكة المكانية.

2. عقد الحدود (الإغلاق عند نفس المستوى):

   • على عكس العقد الداخلية، يتم تحديث عقد الحدود باستخدام معادلة منفصلة عند نفس المستوى عبر إزاحة مكانية $\ell_h \asymp h$.
   • يتم تمثيل شرط الحدود كالتالي:
     $$\frac{U_i - I_h[U](x_i - \ell_h \gamma(x_i))}{\ell_h} + k(t_n, x_i)U_i - g(t_n, x_i) = 0.$$
   • هذا النظام يكون ضمنياً بشكل عام. بالنسبة للبيانات الخطية الثابتة، فإنه ينتج مصفوفة M (M-matrix) متفرقة وغير منفردة. بالنسبة للمسائل المتحكم بها، يكون خريطة الحدود عبارة عن تطبيق تقليص (contraction mapping)، مما يضمن القابلية للحل.

3. التعامل مع التحكم:

   • يتعامل المخطط مع مجموعة التحكم $A \times B$ عن طريق أخذ الحد الأدنى (أو الحد الأعلى، اعتماداً على الصيغة) عبر الفروع بعد تكوين متوسط الفرع. يتم اختيار عناصر التحكم من جديد عند كل عقدة وخطوة زمنية.

المساهمات الرئيسية

1. الإيجابية دون الحاجة للهيمنة القطرية: يولد المخطط معاملات غير سالبة مباشرة من البناء الاحتمالي (احتمالات متساوية وعوامل توهين غير سالبة). يمكنه التعامل مع الانتشار متباين الخواص والمتدهور دون فرض شروط هيمنة قطرية لقوالب مدمجة أو علاقات من نوع CFL بين $\Delta t$ و $h$.
2. المعالجة الهندسية لروبن: يعالج المخطط شروط روبن المائلة عبر الانعكاس المرآتي والزمن المحلي المنفصل دون الحاجة إلى معامل إزاحة داخلي إضافي (على عكس مناهج "سي-لاغرانجيان" شبه المستقرة). ترتبط بيانات روبن فقط عبر التوهين، مما يترك احتمالات الفروع دون تغيير.
3. الرتابة والاستقرار: المخطط يحافظ على الرتبة (monotone) ومستقر في معيار $L^\infty$. يشكل إغلاق الحدود تطبيق تقليص، بينما يكون التحديث الداخلي صريحاً.
4. نظرية التقارب: تحت شرط التكرار $\Delta t, h \to 0$ مع $h^2/\Delta t \to 0$، يتقارب المخطط إلى الحل الفريد (viscosity solution) عبر حجة الحدود شبه المريحة (Barles–Souganidis framework).
5. تقديرات الخطأ: بالنسبة للحلول الكلاسيكية $u \in C^3_b$ ذات بقايا حدودية منتظمة من نوع $C^{1,1}$، يتم اشتقاق حد الخطأ كالتالي:
   $$\max |U^n_i - u(t_n, x_i)| \le C \left( \Delta t^{1/2} + \frac{h^2}{\Delta t} \right).$$
   موازنة الحدود باستخدام $\Delta t \asymp h^{4/3}$ تعطي رتبة قدرها $O(h^{2/3})$، بينما الاختيار القياسي $\Delta t \asymp h$ يعطي $O(h^{1/2})$.

النتائج العددية
تسجل الورقة أربع تجارب ثنائية الأبعاد على قرص الوحدة ومستطيل مختلط النطاق:

• المقارنة ($k=0$): تمت مقارنة المخطط مع تنفيذ مرجعي معاير (Calzola et al.) لمسألة $k=0$. يحقق المنهج المقترح دقة تنافسية (رتب $\approx 0.96$) دون الحاجة إلى معلمة إزاحة داخلية قابلة للضبط، مما يظهر سلوكاً أكثر توازناً في الشبكات الدقيقة.
• روبن متغير ($k \neq 0$): تؤكد اختبارات معاملات روبن المتغيرة استقرار المخطط. جدول الخطوات الزمنية "المتوازن" ($\Delta t \propto h^{4/3}$) يقلل الأخطاء بنسبة $\approx 55\%$ مقارنة بـ $\Delta t = h$ ويعطي رتباً مرصودة قدرها $1.242\text{--}1.261$ (في المرحلة قبل التقارب)، وهي أدق من الحد النظري $O(h^{2/3})$.
• الحساسية الزاوية: اضمحلال الخطأ بالنسبة لعدد عناصر التحكم الزاوية $K$ يتوافق مع $O(K^{-2})$.
• الحالات غير السلسة/المختلطة: تظهر تجارب الحل لزوجة المخروط (الانتشار المتدهور) ومسألة روبن-ديريكلت المختلطة اضمحلال خطأ رتيب، وإن كانت معدلات التقارب أقل ولا تغطيها نظريات النطاق السلس.

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

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