pith. sign in

arxiv: 2604.24030 · v2 · submitted 2026-04-27 · 🧮 math.NA · cs.NA

Numerical Analysis of a Variable-Order Time-Fractional Incompressible Magnetohydrodynamics System

Abdumauvlen Berdyshev , Dossan Baigereyev , Aibek Bakishev , Nurlana Alimbekova , Talgat Farkhadov This is my paper

Pith reviewed 2026-05-08 02:18 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords variable-order Caputo derivativestime-fractional MHDL1 approximationfinite element methoddiscrete fractional Grönwall theoremstability analysisconvergenceenergy estimates
0
0 comments X

The pith

A convergent finite-element scheme for variable-order time-fractional MHD shows that fractional orders alter energy and enstrophy evolution even at fixed Reynolds number.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a numerical method for incompressible magnetohydrodynamics equations in which classical time derivatives are replaced by variable-order Caputo fractional derivatives to capture time-varying memory. It introduces a fully discrete scheme that combines finite-element discretization in space with an L1-type approximation of the fractional operators in time. Stability estimates are derived and convergence is established by verifying that the approximation kernels meet the hypotheses of an abstract discrete fractional Grönwall theorem, which is then applied to the coupled system. Numerical tests examine temporal accuracy, consistency with the integer-order case, and the effects of different order profiles on energy and enstrophy histories in a periodic vortex benchmark. The results indicate that variable orders produce measurable differences in these quantities without any change in the Reynolds number.

Core claim

The kernels generated by the variable-order L1 approximation satisfy the assumptions of an abstract discrete fractional Grönwall theorem, which is applied to prove stability and convergence for the fully discrete finite-element solution of the time-fractional incompressible MHD system. Numerical experiments confirm consistency with the classical MHD equations as orders approach one and demonstrate that variable-order profiles noticeably affect the evolution of kinetic energy, enstrophy, and current enstrophy even when the Reynolds number remains fixed.

What carries the argument

Variable-order L1 approximation of Caputo derivatives, whose generated kernels are shown to fulfill the conditions of an abstract discrete fractional Grönwall theorem used to obtain stability and error bounds for the coupled finite-element MHD discretization.

If this is right

  • The fully discrete scheme satisfies a discrete stability estimate for admissible variable-order functions.
  • Solutions of the fractional scheme converge to those of the classical MHD system in appropriate norms as the orders approach one.
  • Different variable-order profiles produce distinct energy, enstrophy, and current-enstrophy histories in the divergence-free vortex benchmark.
  • Global indicators such as time-integrated diagnostics vary systematically with the parameters that define the variable orders.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The kernel-verification technique may extend to stability proofs for other coupled systems that incorporate variable-order fractional operators.
  • Adaptive choice of the order function during a simulation could be explored to match local solution features.
  • Physical experiments on fluids or plasmas with time-dependent memory could be compared against the numerical energy histories to test model relevance.

Load-bearing premise

The kernels produced by the variable-order L1 approximation satisfy the positivity, monotonicity, and summability conditions required by the abstract discrete fractional Grönwall theorem.

What would settle it

A specific variable-order function for which the computed discrete solution diverges or the observed temporal convergence rate falls below the predicted order in the representative-order temporal tests.

Figures

Figures reproduced from arXiv: 2604.24030 by Abdumauvlen Berdyshev, Aibek Bakishev, Dossan Baigereyev, Nurlana Alimbekova, Talgat Farkhadov.

Figure 1
Figure 1. Figure 1: shows the time evolution of ∥∇ · uh∥ and ∥∇ · Bh∥ for Cases 1–3. In all cases, the velocity divergence remains small, and for the magnetic field the cleaning step reduces the divergence substantially over the whole time interval. Hence, the divergence errors remain small in the convergence tests. This behavior is consistent with the use of H1 -conforming finite element spaces, for which the constraints ∇ ·… view at source ↗
Figure 2
Figure 2. Figure 2: shows that both E (ε) u (t) and E (ε) B (t) decrease monotonically as ε decreases. For larger ε, the deviations from the integer-order solution grow rapidly at early times and remain noticeable over the interval shown. As ε becomes smaller, the curves are shifted downward over the whole time interval, and for ε = 10−6 the discrepancies are reduced to the level of about 10−7–10−5 . This behavior is consiste… view at source ↗
Figure 3
Figure 3. Figure 3: shows the time evolution of |Kε(t) − K1 (t)| and |Mε(t) − M1 (t)|. In both panels, the discrepancies decrease as ε becomes smaller. For the magnetic energy, the curves remain or￾dered over the whole interval, with smaller ε giving uniformly smaller values of |Mε (t) − M1 (t)|. For the kinetic energy, the same overall trend is observed, although the curves pass through val￾ues close to zero near the middle … view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of the divergence norms in the linear ramp case (Case 2) used for selecting the stabilization parameters. (a) ∥∇ · uh∥ for several values of ζ with χ = 1; (b) ∥∇ · Bh∥ for several values of χ with ζ = 2000. Dashed lines correspond to the values before cleaning and solid lines correspond to the values after cleaning. The corresponding results are collected in view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of the diagnostics for the classical model (α (t) = β (t) ≡ 1) and the constant-order fractional cases, showing the effect of α on the evolution of the main energy and enstrophy measures: (a) Kinetic energy K (t); (b) Magnetic energy M (t); (c) Enstrophy Z (t); (d) Current enstrophy J (t). The kinetic energy K (t) behaves differently. The cases α = 0.6 and α = 0.75 decrease faster at early times… view at source ↗
Figure 6
Figure 6. Figure 6: compares the variable-order profiles (Cases 2–5) with the classical model using the same diagnostics. In all four cases, the magnetic energy M (t) and the current enstrophy J (t) decay much faster than in the classical solution. At the same time, the variable-order curves are not identical, which shows that the results depend not only on the values of α (t), but also on how α (t) changes in time view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of the diagnostics across classical model (α (t) = β (t) ≡ 1) and unsym￾metric variable-order fractional cases, showing the effect of α (t) = β (t) on energy levels and small-scale activity: (a) Kinetic energy K (t); (b) Magnetic energy M (t); (c) Enstrophy Z (t); (d) Current enstrophy J (t). For the constant-order cases, all four deviations decrease as α increases from 0.6 to 0.9. Thus, in the … view at source ↗
Figure 8
Figure 8. Figure 8: Comparison of the diagnostics across different values of the Reynolds numbers Re = Rm: (a) Kinetic energy K (t); (b) Magnetic energy M (t); (c) Enstrophy Z (t); (d) Current enstrophy J (t) view at source ↗
Figure 9
Figure 9. Figure 9: Time evolution of the divergence norms for several representative variable-order profiles. (a) Velocity divergence ∥∇ · uh∥. (b) Magnetic divergence ∥∇ · Bh∥, with dashed lines corresponding to the values before cleaning and solid lines to the values after cleaning. 27 view at source ↗
Figure 10
Figure 10. Figure 10: Heatmaps of the relative deviations of the time-integrated diagnostics for the linear-ramp profile (50). The horizontal and vertical axes correspond to the values of α0 and αT , respectively. (a) δIK, relative deviation of the time-integrated kinetic energy; (b) δIM, relative deviation of the time-integrated magnetic energy; (c) δIZ, relative deviation of the time-integrated enstrophy; (d) δIJ , relative … view at source ↗
read the original abstract

We consider an incompressible magnetohydrodynamics (MHD) model in which the classical first-order time derivatives in the momentum and magnetic induction equations are replaced by variable-order Caputo time-fractional derivatives. This formulation allows the memory effect to vary during the evolution and represents a time-fractional generalization of the incompressible MHD system with nonstationary memory. To approximate the problem, we use a fully discrete scheme combining a finite element discretization in space with an L1-type approximation of the variable-order Caputo operators in time. For this discretization, we establish a discrete stability estimate and also derive an auxiliary corrected discrete energy estimate for the fully discrete solution. Convergence is proved by showing that the kernels generated by the variable-order L1 approximation satisfy the assumptions of an abstract discrete fractional Gr\"{o}nwall theorem, which is then applied to the coupled MHD system. The numerical study consists of four parts. First, representative order profiles are used to examine temporal convergence. Second, consistency with the classical incompressible MHD equations is studied as the fractional orders approach one, using norms of solution differences and deviations in kinetic and magnetic energies. Third, the influence of the variable-order fractional terms on nonlinear evolution is investigated through the periodic divergence-free vortex benchmark, with comparisons based on energy and enstrophy histories, divergence errors, Reynolds-number dependence, and time-integrated diagnostics. Fourth, parameter-space maps show how the parameters defining the variable orders affect global indicators. The results show that the variable orders can noticeably affect the evolution of the energy, enstrophy, and current enstrophy even when the Reynolds number is fixed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 3 minor

Summary. The manuscript analyzes a variable-order time-fractional incompressible MHD system. It introduces a fully discrete scheme that combines finite-element discretization in space with an L1-type approximation of the variable-order Caputo operators in time. Discrete stability and a corrected energy estimate are derived for the scheme. Convergence is established by verifying that the discrete kernels generated by the variable-order L1 approximation satisfy the hypotheses of an abstract discrete fractional Grönwall theorem, which is then applied to obtain error bounds for the coupled velocity and magnetic-field variables. Numerical experiments examine temporal convergence rates, consistency with the classical MHD system as the orders approach unity, the effect of variable orders on energy/enstrophy evolution in a periodic vortex benchmark, and sensitivity to the parameters defining the order functions.

Significance. If the kernel verification and application of the abstract theorem hold, the work supplies a rigorous, reusable framework for the numerical analysis of variable-order fractional MHD systems. The demonstration that variable orders can visibly alter kinetic and magnetic energy histories even at fixed Reynolds number is a concrete, falsifiable observation that may inform modeling choices in applications with time-varying memory. The combination of an abstract convergence argument with targeted numerical diagnostics is efficient and extends existing fractional-calculus techniques to a coupled nonlinear system.

major comments (2)
  1. [Convergence analysis section] The central convergence argument rests on showing that the kernels arising from the variable-order L1 scheme satisfy the hypotheses of the abstract discrete fractional Grönwall theorem. Please supply the explicit verification (positivity, monotonicity, or the precise summation inequalities required by the theorem) for general continuous α(t), β(t) in the relevant section; a brief calculation or reference to the precise conditions used would strengthen the claim.
  2. [Stability and energy estimate section] In the corrected discrete energy estimate, the auxiliary terms introduced to handle the variable-order operators must be controlled uniformly with respect to the mesh parameters. Clarify how the constants in this estimate depend on the bounds of α(t) and β(t) and whether they remain independent of the time-step size.
minor comments (3)
  1. [Introduction] The four-part numerical study is clearly described in the abstract; a short roadmap paragraph at the end of the introduction would improve navigation for readers.
  2. [Model formulation] Notation for the variable-order functions α(t) and β(t) and the associated Caputo operators should be introduced once with all required regularity assumptions and then used consistently.
  3. [Numerical experiments] In the parameter-space maps, state the precise norms or integrated quantities plotted and indicate the range of Reynolds numbers considered.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation and the constructive comments, which have helped clarify the presentation of the convergence and stability results. We address each major comment below and have incorporated the suggested clarifications into the revised manuscript.

read point-by-point responses

  1. Referee: [Convergence analysis section] The central convergence argument rests on showing that the kernels arising from the variable-order L1 scheme satisfy the hypotheses of the abstract discrete fractional Grönwall theorem. Please supply the explicit verification (positivity, monotonicity, or the precise summation inequalities required by the theorem) for general continuous α(t), β(t) in the relevant section; a brief calculation or reference to the precise conditions used would strengthen the claim.

    Authors: We appreciate the referee's suggestion to make the kernel verification more explicit. The original manuscript indicates that the kernels satisfy the required hypotheses but provides only a sketch rather than the full calculation for arbitrary continuous α(t), β(t). In the revised version we have inserted a new lemma (Lemma 3.3) in the convergence section that supplies the explicit verification: we prove positivity of the weights, monotonicity of the discrete kernels, and the precise summation inequalities needed by the abstract theorem, under the standing assumption that α(t) and β(t) are continuous and bounded in a closed subinterval of (0,1). The argument follows directly from the definition of the variable-order L1 weights and the uniform continuity of the order functions on the time interval. revision: yes

  2. Referee: [Stability and energy estimate section] In the corrected discrete energy estimate, the auxiliary terms introduced to handle the variable-order operators must be controlled uniformly with respect to the mesh parameters. Clarify how the constants in this estimate depend on the bounds of α(t) and β(t) and whether they remain independent of the time-step size.

    Authors: We thank the referee for this observation. The constants appearing in the corrected discrete energy estimate depend on the essential infima and suprema of α(t) and β(t) (denoted α_min, α_max, β_min, β_max) but are independent of both the time-step size Δt and the spatial mesh size h. This uniformity follows from the standard bounds on the variable-order L1 kernels and the coercivity properties of the bilinear forms. We have added a short remark immediately after the statement of the energy estimate (Remark 2.4) that records this dependence explicitly and confirms independence from the discretization parameters. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The central proof proceeds by verifying that kernels generated by the variable-order L1 scheme satisfy the hypotheses of a pre-existing abstract discrete fractional Grönwall theorem, then invoking the theorem on the coupled MHD energy estimates to obtain stability and convergence. This is a standard, non-self-referential application of an external abstract result rather than a derivation that reduces to its own inputs by construction. No self-definitional steps, fitted parameters renamed as predictions, load-bearing self-citations, or ansatz smuggling appear in the described chain; the numerical experiments provide independent validation against the classical MHD limit and parameter variations.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim relies on the suitability of the variable-order functions and the applicability of the abstract theorem to the MHD system. No new entities are invented. Since only the abstract is available, details on any fitted parameters are not verifiable.

free parameters (1)
  • variable order functions α(t), β(t)
    Chosen as representative profiles in the numerical study; may function as model parameters.
axioms (1)
  • domain assumption The kernels of the variable-order L1 approximation satisfy the assumptions of the abstract discrete fractional Grönwall theorem.
    This is key to the convergence proof as stated in the abstract.

pith-pipeline@v0.9.0 · 5625 in / 1308 out tokens · 55106 ms · 2026-05-08T02:18:37.413837+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 39 canonical work pages · 1 internal anchor

  1. [1]

    Frac- tional analysis of magnetohydrodynamics Maxwell flow over an inclined plate with the effect of thermal radiation

    Abbas, S., Nisa, Z.U., Gilani, S.F.F., Nazar, M., Metwally, A.S.M., Jan, A.Z., 2024. Frac- tional analysis of magnetohydrodynamics Maxwell flow over an inclined plate with the effect of thermal radiation. International Journal of Theoretical Physics 63. doi:10.1007/ s10773-024-05654-3

    2024

  2. [2]

    An implicit difference scheme for a mixed problem of hyperbolic type with memory

    Abdiramanov, Z.A., Baishemirov, Z.D., Berdyshev, A.S., Shiyapov, K.M., 2024. An implicit difference scheme for a mixed problem of hyperbolic type with memory. Lobachevskii Journal of Mathematics 45, 569–577. doi:10.1134/s1995080224600249

    work page doi:10.1134/s1995080224600249 2024

  3. [3]

    Inverse problem for non-stationary system of mag- netohydrodynamics

    Abylkairov, U.U., Aitzhanov, S.E., 2015. Inverse problem for non-stationary system of mag- netohydrodynamics. Boundary Value Problems 2015. doi:10.1186/s13661-015-0438-x

    work page doi:10.1186/s13661-015-0438-x 2015

  4. [4]

    Analysis of variable-order fractional enzyme kinetics model with time delay

    Agilan, K., Naveen, S., Suganya, S., Parthiban, V., 2025. Analysis of variable-order fractional enzyme kinetics model with time delay. Scientific Reports 15. doi:10.1038/ s41598-025-16382-x

    2025

  5. [5]

    Investigation of MHD fractionalized viscous fluid and thermal memory with slip and Newtonian heating effect: a fractional model based on Mittag-Leffler kernel

    Ali, Q., Amir, M., Metwally, A.S.M., Younas, U., Jan, A.Z., Amjad, A., 2024. Investigation of MHD fractionalized viscous fluid and thermal memory with slip and Newtonian heating effect: a fractional model based on Mittag-Leffler kernel. Journal of Thermal Analysis and Calorimetry 149, 8257–8270. doi:10.1007/s10973-024-13205-5

    work page doi:10.1007/s10973-024-13205-5 2024

  6. [6]

    Numerical method for the variable- order fractional filtration equation in heterogeneous media

    Alimbekova, N., Bakishev, A., Berdyshev, A., 2024. Numerical method for the variable- order fractional filtration equation in heterogeneous media. Fractal and Fractional 8, 640. doi:10.3390/fractalfract8110640

    work page doi:10.3390/fractalfract8110640 2024

  7. [7]

    Numerical identification of a time-dependent coefficient in a time- fractional diffusion equation with integral constraints

    Altybay, A., 2026. Numerical identification of a time-dependent coefficient in a time- fractional diffusion equation with integral constraints. Zeitschrift für angewandte Math- ematik und Physik 77. doi:10.1007/s00033-025-02653-0

    work page doi:10.1007/s00033-025-02653-0 2026

  8. [8]

    A sample implementation for parallelizing divide-and-conquer algorithms on the GPU,

    Arif, M., Kumam, P., Seangwattana, T., Suttiarporn, P., 2023. A fractional model of mag- netohydrodynamics Oldroyd-B fluid with couple stresses, heat and mass transfer: A com- parison among non-Newtonian fluid models. Heliyon 9, e17642. doi:10.1016/j.heliyon. 2023.e17642

    work page doi:10.1016/j.heliyon 2023

  9. [9]

    Magnetohydrodynamics in Materials Processing

    Asai, S., 2011. Magnetohydrodynamics in Materials Processing. Springer Netherlands. pp. 49–86. doi:10.1007/978-94-007-2645-1_3

    work page doi:10.1007/978-94-007-2645-1_3 2011

  10. [10]

    Numerical investiga- tion of fractional Maxwell nano-fluids between two coaxial cylinders via the finite difference approach

    Asjad, M.I., Usman, M., Assiri, T.A., Ali, A., Tag-ElDin, E.M., 2023. Numerical investiga- tion of fractional Maxwell nano-fluids between two coaxial cylinders via the finite difference approach. Frontiers in Materials 9. doi:10.3389/fmats.2022.1050767

    work page doi:10.3389/fmats.2022.1050767 2023

  11. [11]

    AFiD-MHD: A finite difference method for magnetohydrody- namic flows

    Bader, S.H., Zhu, X., 2025. AFiD-MHD: A finite difference method for magnetohydrody- namic flows. Journal of Computational Physics 523, 113658. doi:10.1016/j.jcp.2024. 113658

    work page doi:10.1016/j.jcp.2024 2025

  12. [12]

    Efficient numerical im- plementation of the time-fractional stochastic Stokes–Darcy model

    Baishemirov, Z., Berdyshev, A., Baigereyev, D., Boranbek, K., 2024. Efficient numerical im- plementation of the time-fractional stochastic Stokes–Darcy model. Fractal and Fractional 8, 476. doi:10.3390/fractalfract8080476

    work page doi:10.3390/fractalfract8080476 2024

  13. [13]

    Multidimensional MHD shock tests of regularized smoothed particle hydrodynamics

    Borve, S., Omang, M., Trulsen, J., 2006. Multidimensional MHD shock tests of regularized smoothed particle hydrodynamics. The Astrophysical Journal 652, 1306–1317. doi:10. 1086/508454. 32

    2006

  14. [14]

    Attractivity of fractional functional differential equations

    Chen, F., Zhou, Y., 2011. Attractivity of fractional functional differential equations. Com- puters & Mathematics with Applications 62, 1359–1369. doi:10.1016/j.camwa.2011.03. 062

    work page doi:10.1016/j.camwa.2011.03 2011

  15. [15]

    doi:10.48550/ARXIV.2501.17179

    Chen, Z.M., Dong, B.Q., Zhang, Q., 2025.l 2 decay estimates of weak solutions to 3D fractional MHD equations in exterior domains. doi:10.48550/ARXIV.2501.17179

    work page doi:10.48550/arxiv.2501.17179 2025

  16. [16]

    Functional Fractional Calculus

    Das, S., 2011. Functional Fractional Calculus. Springer Berlin Heidelberg. doi:10.1007/ 978-3-642-20545-3

    2011

  17. [17]

    An Introduction to Magnetohydrodynamics

    Davidson, P.A., 2001. An Introduction to Magnetohydrodynamics. Cambridge University Press. doi:10.1017/cbo9780511626333

    work page doi:10.1017/cbo9780511626333 2001

  18. [18]

    Structure preserving hybrid finite volume finite element method for compressible mhd

    Fambri, F., Sonnendrücker, E., 2025. Structure preserving hybrid finite volume finite element method for compressible mhd. Journal of Computational Physics 523, 113691. doi:10.1016/j.jcp.2024.113691

    work page doi:10.1016/j.jcp.2024.113691 2025

  19. [19]

    Existence of solutions for a singular system of nonlinear fractional differential equations

    Feng, W., Sun, S., Han, Z., Zhao, Y., 2011. Existence of solutions for a singular system of nonlinear fractional differential equations. Computers & Mathematics with Applications 62, 1370–1378. doi:10.1016/j.camwa.2011.03.076

    work page doi:10.1016/j.camwa.2011.03.076 2011

  20. [20]

    Geometry of singular structures in magnetohydrodynamic flows

    Grauer, R., Marliani, C., 1998. Geometry of singular structures in magnetohydrodynamic flows. Physics of Plasmas 5, 2544–2552. doi:10.1063/1.872939

    work page doi:10.1063/1.872939 1998

  21. [21]

    A priori estimates and optimal finite element approximation of the MHD flow in smooth domains

    He, Y., Zou, J., 2018. A priori estimates and optimal finite element approximation of the MHD flow in smooth domains. ESAIM: Mathematical Modelling and Numerical Analysis 52, 181–206. doi:10.1051/m2an/2018006

    work page doi:10.1051/m2an/2018006 2018

  22. [22]

    Splitting-based structure preserving discretizations for magnetohydrodynamics

    Hiptmair, R., Pagliantini, C., 2018. Splitting-based structure preserving discretizations for magnetohydrodynamics. The SMAI Journal of computational mathematics 4, 225–257. doi:10.5802/smai-jcm.34

    work page doi:10.5802/smai-jcm.34 2018

  23. [23]

    Analysis of energy transfer through direct numerical simulations of magnetohydrodynamic Orszag–Tang vortex

    Jadhav, K., Chandy, A.J., 2021. Analysis of energy transfer through direct numerical simulations of magnetohydrodynamic Orszag–Tang vortex. Physics of Fluids 33. doi:10. 1063/5.0051476

    2021

  24. [24]

    On the divergence constraint in mixed finite element methods for incompressible flows

    John, V., Linke, A., Merdon, C., Neilan, M., Rebholz, L.G., 2017. On the divergence constraint in mixed finite element methods for incompressible flows. SIAM Review 59, 492–544. doi:10.1137/15m1047696

    work page doi:10.1137/15m1047696 2017

  25. [25]

    https://doi.org/10.48550/arXiv

    Kayanikhoo, F., Cemeljic, M., Wielgus, M., Kluzniak, W., 2023. Energy dissipation in astrophysical simulations: results of the Orszag-Tang test problem. doi:10.48550/ARXIV. 2312.06675

    work page internal anchor Pith review doi:10.48550/arxiv 2023

  26. [26]

    Numerical simulations and modeling of MHD boundary layer flow and heat transfer dynamics in Darcy-Forchheimer media with dis- tributed fractional-order derivatives

    Khan, M., Zhang, Z., Lu, D., 2023. Numerical simulations and modeling of MHD boundary layer flow and heat transfer dynamics in Darcy-Forchheimer media with dis- tributed fractional-order derivatives. Case Studies in Thermal Engineering 49, 103234. doi:10.1016/j.csite.2023.103234

    work page doi:10.1016/j.csite.2023.103234 2023

  27. [27]

    Numerical study on radiative MHD flow of viscoelas- tic fluids with distributed-order and variable-order space fractional operators

    Li, N., Wang, X., Xu, H., Qi, H., 2024. Numerical study on radiative MHD flow of viscoelas- tic fluids with distributed-order and variable-order space fractional operators. Mathematics and Computers in Simulation 215, 291–305. doi:10.1016/j.matcom.2023.07.021

    work page doi:10.1016/j.matcom.2023.07.021 2024

  28. [28]

    A discrete Grönwall inequality with applications to numerical schemes for subdiffusion problems

    Liao, H.l., McLean, W., Zhang, J., 2019. A discrete Grönwall inequality with applications to numerical schemes for subdiffusion problems. SIAM Journal on Numerical Analysis 57, 218–237. doi:10.1137/16m1175742. 33

    work page doi:10.1137/16m1175742 2019

  29. [29]

    Analysis of a multi-term variable-order time-fractional diffusion equation and its Galerkin finite element approximation

    Liu, H., Zheng, X., Fu, H., 2022. Analysis of a multi-term variable-order time-fractional diffusion equation and its Galerkin finite element approximation. Journal of Computational Mathematics 40, 814–834. doi:10.4208/jcm.2102-m2020-0211

    work page doi:10.4208/jcm.2102-m2020-0211 2022

  30. [30]

    Magnetohydrodynamic analysis and fast calculation for fractional Maxwell fluid with adjusted dynamic viscosity

    Liu, Y., Jiang, M., 2024. Magnetohydrodynamic analysis and fast calculation for fractional Maxwell fluid with adjusted dynamic viscosity. Magnetochemistry 10, 72. doi:10.3390/ magnetochemistry10100072

    2024

  31. [31]

    Numerical simulation and parameter estimation of the space-fractional magnetohydrodynamic flow and heat transfer coupled model

    Liu, Y., Jiang, X., Jia, J., 2024. Numerical simulation and parameter estimation of the space-fractional magnetohydrodynamic flow and heat transfer coupled model. Fractal and Fractional 8, 557. doi:10.3390/fractalfract8100557

    work page doi:10.3390/fractalfract8100557 2024

  32. [32]

    Nonlocal modeling and inverse parameter estimation of time-varying vehicular emissions in urban pollution dynamics

    Madiyarov, M., Alimbekova, N., Bakishev, A., Mukhamediyev, G., Yergaliyev, Y., 2025. Nonlocal modeling and inverse parameter estimation of time-varying vehicular emissions in urban pollution dynamics. Mathematics 13, 2772. doi:10.3390/math13172772

    work page doi:10.3390/math13172772 2025

  33. [33]

    Fractional calculus in bioengineering, part 1

    Magin, R.L., 2004. Fractional calculus in bioengineering, part 1. Critical Reviews in Biomedical Engineering 32, 1–104. doi:10.1615/critrevbiomedeng.v32.10

    work page doi:10.1615/critrevbiomedeng.v32.10 2004

  34. [34]

    2010.Fractional Calculus and Waves in Linear Viscoelasticity

    Mainardi, F., 2010. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduc- tion to Mathematical Models. Imperial College Press. doi:10.1142/p614

    work page doi:10.1142/p614 2010

  35. [35]

    Variable-order fractional derivatives in financial sys- tems

    Malaikah, H., Al-Abdali, J.F., 2025. Variable-order fractional derivatives in financial sys- tems. Applied Mathematics 16, 461–469. doi:10.4236/am.2025.166025

    work page doi:10.4236/am.2025.166025 2025

  36. [36]

    Metzler \ and\ author J

    Metzler, R., Klafter, J., 2000. The random walk’s guide to anomalous diffusion: a fractional dynamics approach. Physics Reports 339, 1–77. doi:10.1016/s0370-1573(00)00070-3

    work page doi:10.1016/s0370-1573(00)00070-3 2000

  37. [37]

    Modeling of turbulent buoyant flow and heat transfer in liquid metals

    Mohamad, A., Viskanta, R., 1993. Modeling of turbulent buoyant flow and heat transfer in liquid metals. International Journal of Heat and Mass Transfer 36, 2815–2826. doi:10. 1016/0017-9310(93)90101-b

    1993

  38. [38]

    Magnetohydrodynamics

    Moreau, R., 1990. Magnetohydrodynamics. Springer Netherlands. doi:10.1007/ 978-94-015-7883-7

    1990

  39. [39]

    An improved spectral homotopy analysismethod forMHDflowin asemi-porouschannel

    Motsa, S.S., Shateyi, S., Marewo, G.T., Sibanda, P., 2011. An improved spectral homotopy analysismethod forMHDflowin asemi-porouschannel. NumericalAlgorithms 60, 463–481. doi:10.1007/s11075-011-9523-0

    work page doi:10.1007/s11075-011-9523-0 2011

  40. [40]

    Small-scale structure of two-dimensional magneto- hydrodynamic turbulence

    Orszag, S.A., Tang, C.M., 1979. Small-scale structure of two-dimensional magneto- hydrodynamic turbulence. Journal of Fluid Mechanics 90, 129–143. doi:10.1017/ s002211207900210x

    1979

  41. [41]

    Approximation of the inductionless mhd problem using a stabilized finite element method

    Planas, R., Badia, S., Codina, R., 2011. Approximation of the inductionless mhd problem using a stabilized finite element method. Journal of Computational Physics 230, 2977–2996. doi:10.1016/j.jcp.2010.12.046

    work page doi:10.1016/j.jcp.2010.12.046 2011

  42. [42]

    Rotating MHD flow and heat transfer of generalized Maxwell fluid through an infinite plate with Hall effect

    Qiao, Y., Xu, H., Qi, H., 2024. Rotating MHD flow and heat transfer of generalized Maxwell fluid through an infinite plate with Hall effect. Acta Mechanica Sinica 40. doi:10.1007/ s10409-023-23274-x

    2024

  43. [43]

    A fractional study with Newtonian heating effect on heat absorbing MHD radiative flow of rate type fluid with application of novel hybrid fractional derivative operator

    Rehman, A.U., Hua, S., Riaz, M.B., Awrejcewicz, J., Xiange, S., 2023. A fractional study with Newtonian heating effect on heat absorbing MHD radiative flow of rate type fluid with application of novel hybrid fractional derivative operator. Arab Journal of Basic and Applied Sciences 30, 482–495. doi:10.1080/25765299.2023.2250063. 34

    work page doi:10.1080/25765299.2023.2250063 2023

  44. [44]

    On the genesis of the Earth’s magnetism

    Roberts, P.H., King, E.M., 2013. On the genesis of the Earth’s magnetism. Reports on Progress in Physics 76, 096801. doi:10.1088/0034-4885/76/9/096801

    work page doi:10.1088/0034-4885/76/9/096801 2013

  45. [45]

    A finite volume spectral element method for solv- ing magnetohydrodynamic (MHD) equations

    Shakeri, F., Dehghan, M., 2011. A finite volume spectral element method for solv- ing magnetohydrodynamic (MHD) equations. Applied Numerical Mathematics 61, 1–23. doi:10.1016/j.apnum.2010.07.010

    work page doi:10.1016/j.apnum.2010.07.010 2011

  46. [46]

    Unsteady MHD flow and heat transfer of fractional Maxwell viscoelastic nanofluid with Cattaneo heat flux and different particle shapes

    Shen, M., Chen, S., Liu, F., 2018. Unsteady MHD flow and heat transfer of fractional Maxwell viscoelastic nanofluid with Cattaneo heat flux and different particle shapes. Chi- nese Journal of Physics 56, 1199–1211. doi:10.1016/j.cjph.2018.04.024

    work page doi:10.1016/j.cjph.2018.04.024 2018

  47. [47]

    High-order spectral scheme with structure maintenance and fast memory algorithm for nonlocal nonlinear dif- fusion equations

    Shiyapov, K., Abdiramanov, Z., Issa, Z., Zhumaseyitova, A., 2026. High-order spectral scheme with structure maintenance and fast memory algorithm for nonlocal nonlinear dif- fusion equations. AppliedMath 6, 54. doi:10.3390/appliedmath6040054

    work page doi:10.3390/appliedmath6040054 2026

  48. [48]

    A review on variable-order frac- tional differential equations: Mathematical foundations, physical models, numerical meth- ods and applications

    Sun, H.G., Chang, A., Zhang, Y., Chen, W., 2019. A review on variable-order frac- tional differential equations: Mathematical foundations, physical models, numerical meth- ods and applications. Fractional Calculus and Applied Analysis 22, 27–59. doi:10.1515/ fca-2019-0003

    2019

  49. [49]

    Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media

    Tarasov, V.E., 2010. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media. Springer Berlin Heidelberg. doi:10.1007/ 978-3-642-14003-7

    2010

  50. [50]

    On fractional differential equations, dimensional analysis, and the double gamma function

    Vaz, J., de Oliveira, E.C., 2025. On fractional differential equations, dimensional analysis, and the double gamma function. Nonlinear Dynamics 113, 34305–34320. doi:10.1007/ s11071-025-11386-8

    2025

  51. [51]

    Pressure and convection robust finite elements for magnetohydrodynamics

    Beirão da Veiga, L., Dassi, F., Vacca, G., 2025. Pressure and convection robust finite elements for magnetohydrodynamics. Numerische Mathematik 157, 1161–1209. doi:10. 1007/s00211-025-01476-5

    2025

  52. [52]

    Nodal-based finite element methods with local projection stabilization for linearized incompressible magnetohydrodynamics

    Wacker, B., Arndt, D., Lube, G., 2016. Nodal-based finite element methods with local projection stabilization for linearized incompressible magnetohydrodynamics. Computer Methods in Applied Mechanics and Engineering 302, 170–192. doi:10.1016/j.cma.2016. 01.004

    work page doi:10.1016/j.cma.2016 2016

  53. [53]

    Convectionheatandmasstransferoffractional MHDMaxwellfluidinaporousmediumwithSoretandDufoureffects

    Zhao, J., Zheng, L., Zhang, X., Liu, F., 2016. Convectionheatandmasstransferoffractional MHDMaxwellfluidinaporousmediumwithSoretandDufoureffects. InternationalJournal of Heat and Mass Transfer 103, 203–210. doi:10.1016/j.ijheatmasstransfer.2016.07. 057. 35

    work page doi:10.1016/j.ijheatmasstransfer.2016.07 2016