Lithium-ion cell modelling: an example of fractional order behavior modeling with distributed time delay models
تُثبت هذه الورقة أن نماذج تأخير الوقت الموزعة الكلاسيكية يمكنها محاكاة ديناميكيات قانون القوة وتأثيرات الذاكرة المرتبطة عادةً بالنماذج ذات الرتب الكسرية بفعالية، مما يوفر بديلاً أبسط وأكثر قابلية للتفسير الفيزيائي لتطبيقات مثل نمذجة خلايا بطاريات الليثيوم أيون.
في العالم الخفي داخل البطارية، لا تتدفق الكهرباء دائمًا مثل الماء عبر أنبوب. فأحيانًا، يشبه تحرك الجسيمات المشحونة بقعة تنتشر ببطء على قطعة من الورق، حيث تتغير سرعة الانتشار بمرور الوقت بطريقة تتحدى التنبؤ البسيط. يطلق العلماء على هذا السلوك اسم "الذاكرة الطويلة"، مما يعني أن النظام يتذكر حالاته الماضية ويسمح لهذا التاريخ بالتأثير على ما سيحدث بعد ذلك. لعقود من الزمن، استخدم الباحثون فرعًا متخصصًا من الرياضيات، يُعرف باسم الحساب الكسري (fractional calculus)، لوصف هذه التأثيرات المعقدة والمستمرة. وقد أصبح هذا الأداة الرياضية هي اللغة القياسية لتفسير سبب استجابة البطاريات، والأنسجة البيولوجية، وبعض المواد بالطريقة التي تقوم بها على فترات زمنية طويلة. ومع ذلك، فإن هذه اللغة الرياضية صعبة التفسير فيزيائيًا بشكل ملحوظ، مما يترك المهندسين غالبًا مع تنبؤات دقيقة ولكن دون صورة واضحة للآليات الفعلية في العمل.
تتحدى دراسة جديدة أجرتها جوسلين ساباتيه في جامعة بوردو فكرة أن هذه الرياضيات المعقدة هي الطريقة الوحيدة لوصف هذه الظواواهر. فقد سعى الباحثة لمعرفة ما إذا كان هناك فئة أقدم وأبسط بكثير من النماذج — وهي النماذج التي تأخذ في الاعتبار التأخيرات المنتشرة عبر الزمن بدلاً من لحظة واحدة — يمكن أن تنتج النتائج نفسها تمامًا. انصب التركيز على خلية بطارية ليثيوم-أيون، من النوع الموجود في المركبات الكهربائية والهواتف الذكية. فعندما يتم شحن البطارية أو تفريغها، ينخفض جهدها بنمط محدد يبدو كقانون القوة (power law)، وهو منحنى يعزوه العلماء تقليديًا إلى القواعد الغامضة للحساب الكسري. اقترحت ساباتيه أن هذا المنحنى نفسه قد يكون ببساطة نتيجة لعمليات انتشار متعددة تحدث في آن واحد، لكل منها سرعتها الخاصة، بدلاً من قاعدة رياضية واحدة غريبة.
ولاختبار ذلك، قام الفريق بتحليل استجابة الجهد لخلية بطارية عند تطبيق تيار ثابت. وراقبوا انخفاض الجهد بمرور الوقت، والذي ينتج عن مزيج من المقاومة الفورية، والتفاعلات السطحية، والانتشار البطيء لأيونات الليثيوم في عمق الجسيمات الصلبة للقطب الكهربائي. ثم بنى الباحثون نموذجين مختلفين للتنبؤ بهذا السلوك. الأول كان نموذجًا قياسيًا من الرتب الكسرية، وهو المفضل حاليًا في الصناعة، والذي يستخدم قوى غير صحيحة لالتقاط تأثير الذاكرة الطويلة. والثاني كان نموذج تأخير زمني موزع، والذي يعامل النظام كمجموعة من التأخيرات البسيطة، كل منها موزون بشكل مختلف، لمحاكاة انتشار الوقت.
كانت النتائج مذهلة. فعندما طابق الفريق كلا النموذجين مع البيانات التجريبية الفعلية، وجدوا أن نموذج التأخير الموزع الأكثر بساطة قد قدم أداءً يضاهي النموذج الكسري المعقد. فقد تنبأ كلا النموذجين بجهد البطارية بدقة متطوّرة ومتطابقة تقريبًا، محققين دقة مطابقة تزيد عن 94 بالمائة. كما تطلب نموذج التأخير الموزع نفس عدد المعلمات القابلة للضبط للوصول إلى هذا المستوى من الدقة، ومع ذلك فعل ذلك دون استدعاء المشتقات الكسرية. في الواقع، أظهرت الدراسة أن المعاملات الرياضية لنموذج التأخير تطورت طبيعيًا لتبدو مثل منحنيات قانون القوة التي يشتهر بها الحساب الكسري، مما يثبت أن الرياضيات المعقدة لم تكن ضرورية تمامًا لتوليد السلوك الملحوظ.
يشير هذا الاكتشاف إلى أن استجابات قانون القوة التي تُلاحظ في البطاريات والأنظمة الأخرى قد لا تكون دليلًا على قانون فيزيائي أساسي وغريب، بل هي النتيجة الطبيعية لتداخل العديد من العمليات الكلاسيكية. ويقدم نموذج التأخير الموزع ميزة واضحة في الوضوح: فبينما يصف الحساب الكسري الذاكرة من خلال مؤثرات مجردة يصعب تصورها، يصف نموذج التأخيرها من خلال توزيع لفترات التأخير التي يمكن ربطها مباشرة بأوقات النقل الفيزيائي، أو معدلات الانتشار، أو أوقات الإقامة. لا تدعي الدراسة أن النماذج الكسرية خاطئة أو عديمة الفائدة، لكنها تجادل بأنها قد تكون تعقيدًا غير ضروري عندما يمكن لنهج كلاسيكي أن يشرح نفس البيانات بدقة متساوية وبديهة فيزيائية أكبر.
من خلال إثبات أن إطارًا كلاسيكيًا يمكنه تكرار سلوك نموذج كسري، تدعو هذه الأبحاث إلى إعادة تقييم كيفية نمذجة العلماء للأنظمة المدفوعة بالذاكرة. وهي تشير إلى أن التعقيد المتزايد للأطر الكسرية الحديثة، مع العديد من تنويعاتها وتعاريفها، قد يكون يحل مشكلات غير موجودة. فإذا كان بإمكان توزيع بسيط للتأخيرات التقاط جوهر ظواهر الذاكرة الطويلة، فإن الطريق لفهم هذه الأنظمة قد لا يكمن في رياضيات أكثر تجريدًا، بل في أوصاف أفضل لكيفية تطور العمليات الفيزيائية بمرور الوقت. لا يعلن هذا العمل انتصارًا نهائيًا لطريقة على أخرى، ولكنه يثبت بقوة أن النهج الأكثر بساطة وقابلية للتفسير هو بديل حيوي وقوي لفهم الديناميكيات البطيئة والمستمرة للعالم الحقيقي.
ملخص تقني: نمذجة خلايا الليثيوم-أيون باستخدام نماذج التأخير الزمني الموزع
بيان المشكلة تعد الديناميكيات ذات الرتب الكسرية واستجابات قانون القوة (power-law) ظواهر منتشرة في الأنظمة المعقدة، بما في ذلك الأجهزة الكهروكيميائية مثل بطاريات الليثيوم-أيون، والمواد اللزجة المرنة، والأنسجة البيولوجية. تقليديًا، يتم نمذجة هذه السلوكيات باستخدام الحساب الكسري، الذي يولد بشكل طبيعي تأثيرات الذاكرة الطويلة واستجابات قانون القوة عبر المعادلات التفاضلية الكسرية. ومع ذلك، فإن التفسير الفيزيائي للحالات ذات الرتب الكسرية، وتحديد الشروط الأولية، ووجود مشغلين كسريين متنافسين (مثل ريمان-ليوفيل و كابوتو)، يفرض تحديات قائمة.
بالتزامن مع ذلك، تشكل نماذج التأخير الزمني الموزع فئة راسخة من الأنظمة ذات الأبعاد اللانهائية القادرة على وصف تأثيرات الذاركة والوراثة من خلال مستمر من قيم التأخير. وبينما تُستخدم هذه النماذج على نطاق واسع في ديناميكيات السكان والتحكم في العمليات، إلا أن قدرتها المحددة على إعادة إنتاج استجابات قانون القوة المميزة للرتب الكسرية لم تحظَ باهتمام كافٍ. وتتمثل المشكلة المركزية التي يتناولها هذا البحث في ما إذا كانت نماذج التأخير الزمني الموزع الكلاسيكية يمكنها إعادة إنتاج ديناميكيات قانون القوة التي تُنسب عادةً إلى المشغلات الكسرية، مما يوفر بديلًا قابلًا للتفسير الفيزيائي للنمذجة الكسرية.
المنهجية يقترح البحث اشتقاقًا نظريًا يربط بين استجابات قانون القوة ونماذج التأخير الزمني الموزع المنفصلة والمتصلة، يليه تحقق تجريبي باستخدام بيانات خلايا الليثيوم-أيون.
الاشتقاق النظري:
النموذج المنفصل: انطلاقًا من استجابة زمنية لقانون القوة تميز نظامًا يُظهر سلوكًا كسريًا (تحديدًا Q(t)∝tα)، يشتق المؤلف نموذج تأخير زمني موزع منفصل. ومن خلال جعل الصيغة التكاملية لقانون القوة منفصلة وتطبيق صيغة أويلر الصريحة، يوضح البحث أن المعاملات الناتجة تعيد إنتاج اضمحلال قانون القوة (m−ν) تقاربيًا. يؤدي هذا إلى معادلة تكرارية حيث تكون الحالة الحالية عبارة عن مجموع مرجح للحالات الماضية، مما يشكل فعليًا نموذج تأخير موزع.
النموذج المتصل: يشير البحث إلى طريقة تقريب مستمر حيث يتم تقريب التكامل الكسري بواسطة دالة نقل تتضمن نواة ηa(t). تم بناء هذه النواة باستخدام حاصل ضرب الأقطاب والأصفار لتقريب t−ν عبر فاصل زمني محدد. يسمح هذا بتمثيل السلوك الكسري كنظام تأخير زمني موزع ضمن النظرية الكلاسيكية للأنظمة الديناميكية.
التطبيق على خلايا الليثيوم-أيون:
يستخدم البحث بيانات تجريبية لخلية ليثيوم-أيون (مدخل التيار، مخرج الجهد) مقدمة من شركة PSA Peugeot-Citroën.
يتم نمذجة ديناميكيات الخلية باستخدام إطار عمل القطب الواحد حيث يتم تفكيك هبوط الجهد إلى مقاومة لحظية، وتأثيرات سعة الطبقة المزدوجة، والاستقطاب المدفوع بالانتشار.
تتم مقارنة ثلاثة نماذج:
النموذج ذو الرتبة الكسرية: يستخدم دوال نقل كسرية (H1(s) و H2(s)) مشتقة من معادلات الانتشار في الهندسة الكروية، مقربة بواسطة أقطاب وأصفار متداخلة.
نموذج التأخير الزمني الموزع المنفصل: ينفذ المعادلة التكرارية المنفصلة المشتقة.
نموذج التأخير الزمني الموزع المتصل: ينفذ صيغة التأخير المستمرة باستخدام تقريب النواة.
يتم إجراء تحديد المعلمات لجميع النماذج باستخدام خوارزمية تحسين "ليفنبرج-ماركواتد" غير الخطية لتقليل جذر متوسط مربع الخطأ (RMSE) بين الجهود المحاكية والمقاسة.
المساهمات الرئيسية
الرابط النظري: يثبت البحث وجود علاقة رياضية مباشرة تظهر أن نموذج التأخير الزمني الموزع ذو الرتبة الصحيحة يمكنه توليد استجابات تظهر نفس خصائص قانون القوة الخاصة بالحساب الكسري.
إطار عمل بديل: يوضح البحث أن نماذج التأخير الزمني الموزع يمكن أن تعمل كبديل قابل للتفسير الفيزيائي للنماذج ذات الرتب الكسرية للظواهر المدفوعة بالذاكرة.
التحقق المقارن: تقدم الدراسة مقارنة جنبًا إلى جنب بين النماذج الكسرية ونماذج التأخير الموزع المطبقة على بطارية ليثيوم-أيون حقيقية، مما يظهر أن كلا النهجين يحققان دقة مماثلة مع عدد مشابه من معاملات الضبط.
النتائج تشير نتائج النمذجة لخلايا الليثيوم-أيون إلى أن:
الدقة: أعادت كل من النماذج ذات الرتب الكسرية ونماذج التأخير الزمني الموزع (المنفصلة والمتصلة على حد سواء) إنتاج ديناميكيات الجهد المقاسة بدقة عالية.
حقق النموذج الكسري (H2) نسبة مطابقة (FIT) بلغت 95.15% وجذر متوسط مربع خطأ (RMSE) قدره 0.00500.
حقق نموذج التأخير الزمني الموزع المنفصل نسبة مطابقة بلغت 95.08% وRMSE قدره 0.00508.
حقق نموذج التأخير الزمني الموزع المتصل نسبة مطابقة بلغت 94.93% وRMSE قدره 0.00523.
كفاءة المعلمات: تطلبت جميع النماذج عددًا محدودًا ومتساويًا من معاملات الضبط (من 4 إلى 6 معاملات اعتمادًا على تعقيد النموذج).
تحليل الخطأ: كانت أخطاء النمذجة (ϵ) لجميع المنهجيات متقاربة بصريًا وإحصائيًا، مما يؤكد أن نماذج التأخير الموزع لا تضحي بالدقة في سبيل التفسير الفيزيائي.
الأهمية والادعاءات يزعم البحث أن سلوكيات قانون القوة الملحوظة في الأنظمة الكهروكيميائية لا ينبغي تفسيرها تلقائيًا كدليل على آلية كامنة ذات رتبة كسرية. بدلاً من ذلك، يمكن أن تنشأ هذه السلوكيات من تجمع عمليات موزعة كلاسيكية تعمل عبر مقاييس زمنية متعددة.
تكمن أهمية هذا العمل في النقاط التالية:
التفسير الفيزيائي: على عكس النماذج الكسرية، التي تعتمد على مشغلات تفاضلية مجردة وتواجه تحديات تتعلق بالشروط الأولية، تمثل نماذج التأخير الزمني الموزع صراحةً تأثير الحالات الماضية من خلال نوى التأخير. ويمكن ربط هذه النوى مباشرة بكميات فيزيائية مثل أوقات النقل، وأوقات الانتشار، وتوزيعات زمن الإقامة.
البساطة الرياضية: تعمل نماذج التأخير الموزع ضمن النظرية الكلاسيكية للأنظمة الديناميكية، متجنبة التعقيدات والتعريفات المتنافسة المتأصلة في الحساب الككسري.
تحول النموذج: تشير النتائج إلى أن نماذج التأخير الزمني الموزع ليست مجرد تقريبات عددية للنماذج الكسرية، بل تشكل نموذجًا متميزًا. فهي توفر وسيلة لالتقاط سلوكيات الذاكرة الطويلة وقانون القوة مع الحفاظ على اتصال واضح بالآليات الفيزيائية الأساسية.
إعادة تقييم التعقيد: يدعو البحث إلى فحص التزايد في تعقيد أطر النمذجة الكسرية (مثل الرتب المتغيرة أو الرتب الموزعة). فإذا كانت هياكل التأخير الموزع الكلاسيكية قادرة على التقاط الديناميكيات الجوهرية، فإن ضرورة هذه الصياغات الأكثر تطورًا لبعض التطبيقات قد تحتاج إلى إعادة تقييم.
في الختام، يفترض البحث أن نماذج التأخير الزمني الموزع تشكل بديلًا ذا صلة وجذابًا لوصف الظواهر المدفوعة بالانتشار والاعتماد على الذاكرة، مما يوفر توازنًا بين أداء النمذجة والتفسير الفيزيائي.