DOI: https://doi.org/10.1038/s41598-026-42447-6
PMID: https://pubmed.ncbi.nlm.nih.gov/41813788
تاريخ النشر: 2026-03-11
المؤلف: Muhammad Zainul Abidin
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
في هذه الدراسة، نحقق في نظام الديناميكا المائية المغناطيسية (MHD) غير القابل للضغط والذي يتضمن مشتق كابوتو الزمني من الدرجة \(0 < \alpha < 1\). يتم إعادة صياغة معادلات الزخم-التحريض المترابطة باستخدام إسقاط مزدوج الالتواء، مما يزيل بشكل فعال كل من الضغط الديناميكي المائي والضغط الزائف المغناطيسي، مما ينتج عنه زوج خالٍ من التباين من مجالات السرعة والمجالات المغناطيسية. ثم نقوم بتفكيك هذه المشكلة المعاد صياغتها من خلال استخدام تقريب طيفي فورييه خالٍ من التباين في الفضاء جنبًا إلى جنب مع تفكيك الالتفاف L1 بخطوات زمنية متغيرة. يضمن هذا النهج أن الطريقة المنفصلة تلتزم بعدم المساواة في الطاقة الحركية-المغناطيسية الكسرية، وتحافظ على قيود التباين بدقة عالية، وتتقارب إلى قانون الطاقة الكلاسيكي MHD عندما \(\alpha \to 1^-\). تم تصميم طريقتنا العددية لتكون كاملة التفكيك وتحافظ على الهيكل، مما يدير بشكل فعال التفرد الأولي الضعيف من خلال شبكات زمنية متدرجة. نثبت الاستقرار غير المشروط ونستخرج حدود خطأ مثالية من الدرجة \(O(\tau^{2-\alpha} + N^{-m})\) على شبكات زمنية عشوائية. تضمن الطريقة أن الطاقة المنفصلة تتلاشى بشكل أحادي، بغض النظر عن عدم انتظام خطوات الزمن. من خلال تجارب عددية واسعة، بما في ذلك اختبارات على دوامات تايلور-جرين المغناطيسية الكسرية وأورزاغ-تانغ، نتحقق من معدلات التقارب النظرية ونظهر كفاءة استراتيجية الخطوات الزمنية التكيفية، مما يؤكد قوة نهجنا في الحفاظ على الخصائص الفيزيائية لنظام MHD.
مقدمة
في هذا القسم، يقدم المؤلفون مخططًا عدديًا يحافظ على الهيكل لنظام الديناميكا المائية المغناطيسية غير القابلة للضغط (TFMHD) ذو المشتقات الزمنية الكسرية، والذي يتم وصفه بمجموعة من المعادلات التي تحكم ديناميات مجالات السرعة \( u \) والمجالات المغناطيسية \( b \) تحت مشتقات زمنية كسرية. يتم طرح المعادلات على مستطيل دوري مزدوج \( \Omega \subset \mathbb{R}^2 \) مع شروط ابتدائية لـ \( u \) و \( b \)، وتدمج شروط حدود دورية جنبًا إلى جنب مع تطبيع للضغط \( p \). يتم استخدام مشتق كابوتو من الدرجة \( 0 < \alpha < 1 \) لنمذجة السلوك الزمني الكسر للنظام، والذي يعد ذا صلة لفهم ظواهر الاضطراب المغناطيسي MHD المتقطع. يسلط المؤلفون الضوء على أهمية نتائجهم من خلال الإشارة إلى دراسات حديثة تظهر صلة معادلات TFMHD في نمذجة اضطراب الرياح الشمسية. كما يلخصون النتائج التحليلية الحالية المتعلقة بوجود وحيدة الحلول، بالإضافة إلى خصائص التقارب لطرق عددية لنماذج ذات صلة. يحدد البحث خاصيتين هيكليتين حرجتين لنظام TFMHD: الحفاظ على قيود السولينويد على مر الزمن واستقرار الطاقة. تثبت النظريات المقدمة في القسم قيود التباين، وتبدد الطاقة المنفصلة، واستقرار الحلول العددية، جنبًا إلى جنب مع تقديرات الخطأ التي تضمن دقة المخطط الكامل المنفصل تحت ظروف معينة على البيانات الأولية وحجم خطوة الزمن.
نقاش
في هذا القسم، يناقش المؤلفون تطوير والتحقق من خوارزمية عددية مصممة للحفاظ على الهيكل الداخلي لمعادلات الديناميكا المائية المغناطيسية (MHD) ذات المشتقات الزمنية الكسرية. يظهرون أنه مع اقتراب المعامل $\alpha$ من 1، يتقارب الدالة الطاقية $E_\alpha$ إلى الطاقة الكلاسيكية لـ MHD، مما يضمن التوافق مع النماذج التقليدية. يتم التأكيد على أهمية الحفاظ على الثوابت على المستوى المنفصل كأمر حاسم لمحاكاة طويلة الأجل، ويتم الإشارة إلى خوارزميات مختلفة تحافظ على الهيكل لمشاكل التفاضل الكسرية، مع تسليط الضوء على مساهماتها في تلاشي الطاقة واستقرارها.
تستخدم الخوارزمية المقترحة إسقاط مزدوج الالتواء لإعادة صياغة المعادلات الحاكمة، مما يقضي بشكل فعال على الضغوط الديناميكية المائية والمغناطيسية لإنتاج زوج من السرعة والمغناطيسية خالٍ من التباين. يتم تحقيق التفكيك الزمني من خلال قاعدة الالتفاف L1 بخطوات متغيرة، بينما يستخدم التفكيك المكاني تقريب طيفي فورييه سولينويدي. يؤكد المؤلفون أن هذه هي أول خوارزمية تحافظ على الهيكل مع نظرية تقارب شاملة لنظام MHD ذو المشتقات الزمنية الكسرية، محققة دقة زمنية مثالية قدرها $2^{-\alpha}$ وتوفير حدود خطأ صارمة للسرعة، المجال المغناطيسي، والضغط. تؤكد التجارب العددية النتائج النظرية، مما يؤكد معدلات التقارب المتوقعة والتلاشي الأحادي للطاقة المنفصلة.
DOI: https://doi.org/10.1038/s41598-026-42447-6
PMID: https://pubmed.ncbi.nlm.nih.gov/41813788
Publication Date: 2026-03-11
Author(s): Muhammad Zainul Abidin
Primary Topic: Fractional Differential Equations Solutions
Overview
In this study, we investigate the incompressible magnetohydrodynamic (MHD) system incorporating a Caputo time-fractional derivative of order \(0 < \alpha < 1\). The coupled momentum-induction equations are reformulated using a double-curl projection, which effectively removes both hydrodynamic pressure and magnetic pseudopressure, resulting in a divergence-free pair of velocity and magnetic fields. We then discretize this reformulated problem by employing a divergence-free Fourier spectral approximation in space alongside a variable-step L1 convolution discretization in time. This approach ensures that the discrete method adheres to a fractional kinetic-magnetic energy inequality, maintains divergence constraints with high precision, and converges to the classical MHD energy law as \(\alpha \to 1^-\). Our numerical method is designed to be fully discrete and structure-preserving, effectively managing the weak initial singularity through graded time meshes. We establish unconditional stability and derive optimal error bounds of order \(O(\tau^{2-\alpha} + N^{-m})\) on arbitrary time meshes. The method guarantees that the discrete energy decays monotonically, independent of the nonuniformity of time steps. Through extensive numerical experiments, including tests on fractional magnetic Taylor-Green and Orszag-Tang vortices, we validate the theoretical convergence rates and demonstrate the efficiency of the adaptive time-stepping strategy, confirming the robustness of our approach in preserving the physical properties of the MHD system.
Introduction
In this section, the authors present a structure-preserving numerical scheme for the time-fractional incompressible magnetohydrodynamic (TFMHD) system, described by a set of equations governing the dynamics of velocity fields \( u \) and magnetic fields \( b \) under fractional time derivatives. The equations are posed on a doubly periodic rectangle \( \Omega \subset \mathbb{R}^2 \) with initial conditions for \( u \) and \( b \), and they incorporate periodic boundary conditions along with a normalization for pressure \( p \). The Caputo fractional derivative of order \( 0 < \alpha < 1 \) is utilized to model the time-fractional behavior of the system, which is relevant for understanding intermittent MHD turbulence phenomena. The authors highlight the significance of their findings by referencing recent studies that demonstrate the relevance of TFMHD equations in modeling solar-wind turbulence. They also summarize existing analytical results regarding the existence and uniqueness of solutions, as well as convergence properties of numerical methods for related models. The paper outlines two critical structural properties of the TFMHD system: the preservation of solenoidal constraints over time and energy stability. Theorems presented in the section establish the divergence constraints, discrete energy dissipation, and stability of the numerical solutions, alongside error estimates that ensure the accuracy of the fully discrete scheme under specific conditions on the initial data and time-step size.
Discussion
In this section, the authors discuss the development and validation of a numerical algorithm designed to preserve the intrinsic structure of time-fractional incompressible magnetohydrodynamic (MHD) equations. They demonstrate that as the parameter $\alpha$ approaches 1, the energy functional $E_\alpha$ converges to the classical MHD energy, ensuring compatibility with traditional models. The preservation of invariants at the discrete level is emphasized as crucial for long-term simulations, and various structure-preserving algorithms for fractional differential problems are referenced, highlighting their contributions to energy decay and stability.
The proposed algorithm employs a double-curl projection to reformulate the governing equations, effectively eliminating hydrodynamic and magnetic pressures to yield a divergence-free velocity-magnetic pair. Temporal discretization is achieved through a variable-step L1 convolution rule, while spatial discretization utilizes a solenoidal Fourier spectral approximation. The authors assert that this is the first structure-preserving algorithm with a comprehensive convergence theory for the time-fractional MHD system, achieving optimal temporal accuracy of $2^{-\alpha}$ and providing rigorous error bounds for the velocity, magnetic field, and pressure. Numerical experiments validate the theoretical results, confirming the expected rates of convergence and the monotonic decay of discrete energy.