🔬 physics

Similarity Solution of the Fractional Time Independent Advection-Diffusion Equation in Two Dimensions

تقدم هذه الورقة حلاً تشابهياً لمعادلة الانتشار والتحريك ثنائية الأبعاد ذات الزمن المستقل والكسري لنمذجة تركيزات ثاني أكسيد الكبريت في الغلاف الجوي تحت ظروف رياح وعدم استقرار متغيرة، مما يثبت أن هذه الطريقة تعطي تنبؤات دقيقة تقع إلى حد كبير ضمن عامل قدره اثنين من ملاحظات تجربة "برايري-جراس" (Prairie-Grass).

المؤلفون الأصليون: Khaled S. M. Essa, Ahmed M. Mosallem, Fawzia Mubarak

نُشر 2026-08-19

📖 3 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Khaled S. M. Essa, Ahmed M. Mosallem, Fawzia Mubarak

البحث الأصلي مرخَّص بموجب CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

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

ولاختبار نهجهم الجديد، لجأ الفريق إلى مجموعة شهيرة من البيانات الواقعية تُعرف بتجربة "برايري غراس" (Prairie-Grass)، التي أُجريت في أو نيل، نبراسكا، عام 1956. في تلك الاختبارات التاريخية، أطلق العلماء غاز ثاني أكسيد الكبريت من مدخنة قصيرة وقاسوا كمية الغاز التي وصلت إلى الأرض على مسافات مختلفة باتجاه الريح. استخدم الباحثون نموذجهم الكسري لمحاكاة هذه الظروف بدقة، وحسبوا التركيزات المتوقعة للغاز عند مسافات 50 و100 و200 و400 و800 متر. وقارنوا نتائجهم بالقياسات الفعلية التي أُخذت خلال التجارب الأصلية، مع تشغيل الأرقام بإعدادات مختلفة لنموذجهم الكسري لمعرفة أي نسخة تتوافق مع الواقع بشكل أفضل.

أظهرت النتياج أن النهج الكسري كان فعالاً للغاية، لا سيما عند ضبطه على إعداد محدد. ووجد الباحثون أن نموذجهم يمكنه التنبؤ بتركيز ثاني أكسيد الكبريت بدرجة عالية من الموثوقية عبر معظم المسافات التي تم اختبارها. وفي الغالبية العظمى من الحالات، وقعت القيم المتوقعة ضمن عامل اثنين من القياسات المرصودة، وهو معيار مرجعي للدقة في علوم الغلاف الجوي. وهذا يعني أن النموذج كان قريباً بما يكفي ليكون مفيداً لاتخاذ القرارات العملية. وكان الأداء الأفضل عندما استخدم النموذج قيمة كسرية محددة قدرها 0.85؛ فعند هذا الإعداد، طابق النموذج البيانات المرصودة بشكل استثنائي عند مسافات 50 و100 و200 و800 متر، مما سمح برصد سلوك السحابة بدقة أكبر من النماذج التقليدية غير الكسرية.

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

ملخص تقني: حل التشابه لمعادلة الانتشار والتحريك الكسرية المستقلة عن الزمن في بعدين

بيان المشكلة
تتناول الورقة البحثية تحدي نمذجة تشتت الملوثات، وتحديداً ثاني أكسيد الكبريت (SO₂)، في بيئة جوية غير مستقرة. وبينما تُستخدم معادلة الانتشار والتحريك (ADE) الكلاسيكية على نطاق واسع، يشير المؤلفون إلى أنها غالباً ما تفشل في رصد سلوكيات الانتشار الشاذة والمعقدة الموجودة في البيئات المختلطة. ولمعالجة ذلك، تركز الدراسة على معادلة الانتشار والتحريك الكسرية (FADE) في بعدين (المسافة في اتجاه الريح $x$ والارتفاع الرأسي $z$ ). وتتضمن المشكلة المحددة حل معادلة (FADE) الثابتة زمنياً تحت ظروف عدم استقرار جوي، حيث تتغير سرعة الرياح ( $u$ ) والانتشار الدوامي الرأسي ( $k_z$ ) مع الارتفاع. والهدف هو اشتقاق حل تشابه يأخذ في الاعتبار هذه التغيرات والتحقق من صحة النموذج مقابل البيانات التجريبية من تجربة "برايري غراس" (Prairie-Grass).

المنهجية
يستخدم المؤلفون طريقة حل التشابه لتقليل تعقيد المعادلة التفاضلية الجزئية (PDE) إلى معادلة تفاضلية عادية (ODE). ويتم بناء الإطار الرياضي بناءً على الخطوات التالية:

المعادلة الحاكمة: تستخدم الدراسة معادلة (ADE) الكسرية في حالة الاستقرار:
$u \frac{\partial^\alpha C}{\partial x^\alpha} = \frac{\partial}{\partial z} \left( k_z \frac{\partial C}{\partial z} \right)$
حيث يمثل $0 < \alpha < 1$ الرتبة الكسرية، و $C$ تركيز الملوث، و $u$ و $k_z$ هما سرعة الرياح والانتشار الدوامي المعتمدان على الارتفاع، على التوالي.
المعلمات الجوية: تحت الظروف غير المستقرة، يتم تحديد ملف سرعة الرياح $u(z)$ والانتشار الدوامي $k_z(z)$ باستخدام نظرية تشابه "مونين-أوبوكوف"، والتي تتضمن سرعة الاحتكاك ( $u_*$ )، وارتفاع الخشونة ( $z_0$ )، وطول "مونين-أوبوكوف" ( $L$ ).
تحويل التشابه: يقدم المؤلفون متغير تشابه $\mu = z/\delta(x)$ ودالة قياس $\delta(x)$ (سمك التشابه) لموازنة حد التحريك الكسري مع حد الانتشار الرأسي. يؤدي هذا إلى اختزال المعادلة الجزئية (PDE) إلى معادلة تفاضلية عادية (ODE) من الدرجة الثانية ذات معاملات ثابتة.
- يتم تقييم سرعة الرياح الفعالة ( $u_{eff}$ ) والانتشار الفعال ( $k_{eff}$ ) عند ارتفاع المدخنة ( $h_s$ ).
- يتم اشتقاق سمك التشابه كـ $\delta(x) \propto x^{\alpha/2}$ .
الشروط الحدية: يفرض الحل شروط التدفق الصفري (نيومان) عند الأرض ( $z=0$ ) وعند ارتفاع الخلط ( $z=h$ ). ويتم بناء الحل التحليلي النهائي باستخدام طريقة السلسلة الصورية (image series method) لتلبية هذه الحدود، مما ينتج عنه مجموع من النوى الغاوسية (Gaussian kernels).
التحقق من الصحة: تتم مقارنة التركيزات الكسرية المشتقة مع تركيزات SO₂ المرصودة من تجربة برايري غراس (Barad, 1958)، التي أجريت في أو نيل، نبراسكا. وتغطي المقارنة مسافات في اتجاه الريح تبلغ 50، 100، 200، 400، و800 متر تحت ظروف جوية غير مستقرة.

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

الاشتقاق التحليلي: توفر الورقة حلاً تحليلياً مغلقاً لحل التشابه لمعادلة (ADE) الكسرية في بعدين وثابتة زمنياً، مع مراعاة صريحة للتغير الرأسي في سرعة الرياح والانتشار الدوامي في الأجواء غير المستقرة.
تحليل الرتبة الكسرية: تقيم الدراسة أداء النموذج عبر رتب كسرية مختلفة ( $\alpha = 0.85, 0.95,$ والحالة الكلاسيكية $\alpha = 1$ ).
التقييم الإحصائي: يتم قياس دقة النموذج باستخدام مؤشرات إحصائية قياسية: الانحياز الجزئي (FB)، متوسط مربع الخطأ الموزون (NMSE)، معامل الارتباط (COR)، ومعيار عامل الاثنين (FAC2).

النتائج
تؤدي المقارنة بين التنبؤات الكسرية وبيانات تجربة "برايري غراس" إلى النتائج التالية:

الاتفاق العام: بالنسبة لمعظم المسافات، تقع التركيزات المتوقعة كسريةً ضمن عامل اثنين من القيم المرصودة، رغم وجود انحرافات في مسافات محددة.
شذوذ المسافة: يعمل النموذج بشكل جيد عند مسافات 50، 100، 200، و800 متر. ومع ذلك، عند مسافة 400 متر، تقع التركيزات المتوقعة خارج نطاق عامل الاثنين لجميع قيم $\alpha$ المختبرة.
الرتبة الكسرية المثلى: تشير النتائج إلى أن النموذج الكسري مع $\alpha = 0.85$ يوفر أفضل اتفاق مع البيانات المرصودة مقارنة بـ $\alpha = 0.95$ $α = 0.95$ والحالة الكلاسيكية ( $\alpha = 1$ $α = 1$ ).
- عند $x=50$ م، بلغت قيمة FAC2 لـ $\alpha=0.85$ هي 0.63.
- عند $x=100$ م، بلغت قيمة FAC2 لـ $\alpha=0.85$ هي 0.99.
- عند $x=200$ م، بلغت قيمة FAC2 لـ $\alpha=0.85$ هي 1.5 (كما ورد في الجدول 4، وهي تتجاوز النطاق القياسي [0,1]).
- عند $x=800$ م، بلغت قيمة FAC2 لـ $\alpha=0.85$ هي 0.84، مما يعني أن 16% من نقاط البيانات تقع خارج نطاق عامل الاثنين.
الأداء الإحصائي: تظهر المؤشرات الإحصائية (NMSE, FB, COR) عموماً أداءً أفضل للحالات الكسرية ( $\alpha < 1$ ) مقارنة بالحالة الكلاسيكية ( $\alpha = 1$ )، لا سيما عند المسافات القصيرة.

الأهمية والادعاءات
يخلص المؤلفون إلى أن حل التشابه العددي يعطي إجابات موثوقة بمستوى مقبول من الدقة لنمذجة تشتت الملوثات في الظروف الجوية غير المستقرة. وتدعي الدراسة أن النهج الكسري، وتحديداً مع $\alpha = 0.85$ ، يقدم تمثيلاً أكثر دقة لعمليات الانتشار الفيزيائية مقارنة بالنموذج ذي الرتبة الصحيحة الكلاسيكية ( $\alpha = 1$ ). وتؤكد الورقة أن طريقة التشابه تنجح في تقليل التعقيد الحسابي مع الحفاظ على دقة عالية، مما يجعلها أداة قابلة للتطبيق في تقييمات الأثر البيئي المتعلقة بتشتت SO₂. ويصرح المؤلفون صراحة بأن النتائج موثوقة، لكنهم يشيرون بتواضع إلى التباين عند مسافة 400 م، مما يوحي بأنه بينما الطريقة قوية، فإن الظروف المحلية المحددة عند تلك المسافة قد تتطلب مزيداً من الاستقصاء.

ملخص تقني: حل التشابه لمعادلة الانتشار والتحريك الكسرية المستقلة عن الزمن في بعدين

أبحاث مشابهة