Home Analyzing the role of length and radius of MWCNTs in a nanofluid flow influenced by variable thermal conductivity and viscosity considering Marangoni convection
Article Open Access

Analyzing the role of length and radius of MWCNTs in a nanofluid flow influenced by variable thermal conductivity and viscosity considering Marangoni convection

  • Muhammad Ramzan EMAIL logo , Ibtehal Alazman , Nazia Shahmir and Wei Sin Koh
Published/Copyright: September 20, 2025

Abstract

Carbon nanotubes (CNTs), formed by rolling graphene sheets into cylindrical shapes, are widely used in the development of next-generation lithium-ion batteries due to their outstanding thermal conductivity and large surface area, which enhance electrode performance, enable faster charging, and improve energy capacity. Inspired by such a fascinating application of CNTs, the present investigation examines the flow behavior of an aqueous-based nanofluid containing multi-walled carbon nanotubes (MWCNTs), utilizing the well-established Yamada–Ota thermal conductivity model to assess the influence of nanoparticle radius and length. The proposed flow model novelty lies in incorporating Marangoni convection and temperature-dependent thermophysical properties. The governing partial differential equations are transformed into ordinary differential equations through appropriate similarity transformations. Numerical solutions are obtained using the boundary value problem solver bvp4c, with results illustrated through graphs and tables. The findings indicate that thermal conductivity reaches its peak at specific nanoparticle length L = 25 × 10 6 and radius values R = 150 × 10 9 . Optimizing the length-to-radius ratio of MWCNTs in nanofluids can greatly improve thermal conductivity in high-performance electronic cooling systems as extended nanotube lengths enhance heat transfer, while increased radii tend to reduce it. Additionally, Marangoni convection leads to a reduction in temperature distribution while enhancing the heat transfer rate. An increase in the variable temperature parameter results in a corresponding rise in the temperature profile. Validation of the proposed model is also included in this study.

Nomenclature

A

positive dimensional constant

C 1 , C 2

constant quantities

C p

specific heat capacity

E

variable thermal conductivity parameter

f

dimensionless velocity

K NF ( T )

variable thermal conductivity of nanofluid

L

length of MWCNTs

M A

Marangoni parameter

m

exponent constant

Nu

Nusselt number

Pr

Prandtl number

q c

heat flux

R

radius of MWCNTs

s 1

nanoparticles

T

temperature

T o

ambient temperature

x , y

coordinate axis

θ R

variable viscosity parameter

Φ

particle volume fraction

σ o

surface tension at interface

ϒ T

temperature coefficient of surface tension

θ

dimensionless temperature

ρ

density

μ NF ( T )

variable viscosity of nanofluid

ψ

stream function

σ

surface tension

ζ

similarity variable

Φ s 1

particle volume fraction of nanoparticles

1 Introduction

Variable thermal conductivity and viscosity are crucial as they significantly influence fluid rheological properties and heat transfer performance. These characteristics are essential for designing efficient nanofluid systems used in applications like heat exchangers and renewable energy systems. Choosing variables over constant properties ensures greater accuracy for specific applications. While constant properties suffice for simple models or small temperature variations, temperature-dependent characteristics are essential for accurately modeling complex scenarios or high-temperature conditions. These variations help determine fluid flow behavior and predict actual thermal transport phenomena. Researchers increasingly emphasize the importance of temperature-dependent viscosity and thermal conductivity in advanced studies. Abu-Nada [1] conducted a groundbreaking theoretical study on nanofluid flow using an aluminum-water mixture based on the Tiwari and Das model, incorporating variable features with natural convection. Mahmoud and Megahed [2] then investigated the magnetohydrodynamic flow of a power-law non-Newtonian fluid film over an extended surface, considering unique variable characteristics. Subsequently, Abu-Nada et al. [3] numerically analyzed the flow of Al₂O₃–water and CuO–water nanofluid combinations in a differentially heated enclosure with variable properties. Roslan et al. [4] followed with an intriguing study using the finite difference technique to evaluate heat transfer rates in a buoyancy-driven nanofluid-filled trapezoidal enclosure under variable conditions. Sebdani et al. [5] numerically examined mixed convective nanofluid flow in a square cavity, incorporating variable viscosity and thermal conductivity. More recently, Shamshuddin and Eid [6] applied Chebyshev spectral numerical methods to study magnetized aqueous nanofluid flow with ferromagnetic nanoparticles over a stretchable spinning disk under variable features. Additional recent works exploring variable characteristics in diverse scenarios are provided in previous literature [7,8,9,10].

The Marangoni effect signifies mass transfer driven by a surface tension gradient at the interface between two liquids. When this gradient arises from temperature differences, it is identified as thermo-capillary convection. The phenomenon was first described in 1855 by physicist James Thomson, who observed it in the “tears of wine.” In this process, regions with higher surface tension pull more strongly on surrounding particles, causing the liquid to flow from areas of high to low surface tension. This gradient can result from temperature or concentration variations. The Marangoni effect has applications in crystal growth, metallurgy, soap film stabilization, Benard cells, and integrated circuit manufacturing, where it aids in drying silicon wafers. Golia and Viviani [11] introduced the concept of the Marangoni boundary layer in non-isobaric conditions, which was later explored by Montgomery and Wang [12] for a flat plate geometry. Arifin et al. [13] extended this work by studying nanoliquid flow influenced by Marangoni convection over a static flat plate using the Tiwari and Das model. Their findings showed that nanoliquids with titanium dioxide exhibited lower thermal diffusivity compared to those with copper nanoparticles. Hayat et al. [14] analytically investigated Carbon nanotube (CNT)-based nanofluid flows under Marangoni convection and thermal radiative heat flux, concluding that increasing suction reduced velocity distribution for both CNTs and nanofluids. Qayyum [15] examined the Marangoni convection flow of hybrid nanoliquids in two base fluids under gravity, finding that velocity increased with the Marangoni number, while temperature decreased. Abdullah et al. [16] analyzed the electrically conducting flow of hybrid nanoliquids along a porous extending plate under Marangoni convection, observing significant temperature increases with higher magnetic parameters. Recently, Anusha et al. [17] analyzed heat and mass transfer in Marangoni convective nanofluid flow through porous media, reporting increased axial and transverse velocities with rising Marangoni numbers. Additional recent works on Marangoni convection are cited in [18,19,20].

A nanofluid is a combination of a conventional liquid (e.g., water or oil) and suspended nanoparticles (such as oxides, carbides, or CNTs) with sizes below 100 nm. Recognized as next-generation fluids, nanofluids exhibit superior heat transfer properties, making them invaluable for various engineering and industrial applications, including electronics, food processing, medical technologies, and energy production. Recognizing the inadequate heat transfer capabilities of conventional liquids in high-tech processes, Maxwell [21] proposed the idea of mixing solid particles (particles with sizes ranging from millimeters and micrometers) with traditional liquids in 1881, to improve heat transfer capabilities. However, this concept initially faced criticism due to drawbacks such as clogging, sedimentation, and corrosion. In 1995, Choi and Eastman [22] pioneered the idea of nanofluids, updating Maxwell’s theory after nearly a century. Their idea involved adding nano-sized metallic nanoparticles to conventional liquids, rather than using millimeter or micrometer-sized particles. CNTs are considered the most effective nanoparticles. Made from graphene sheet structures, CNTs are grouped as either single-walled (SWCNTs) or multi-walled (MWCNTs). SWCNTs consist of a single layer with a diameter varying from 0.5 to 1.5 nm, while MWCNTs contain interconnected tubes with progressively larger diameters. CNTs exhibit thermal, mechanical, and physicochemical properties that are six times superior, making them crucial in chemical production, optics, microelectronic cooling, and material science. Researchers have shown keen interest in exploring nanofluid flows containing CNTs over varied geometries. Yan et al. [23] experimentally studied the trend of CNTs-based hybrid non-Newtonian nanofluid flow at different temperatures. They determined that viscosity enhances for greater estimations of particle volume fraction. Harish Babu et al. [24] analyzed the inclined magnetic flux on the hybrid nanoliquid along an exponentially stretching sheet. They came to the conclusion that compared to a fundamental liquid, the temperature distribution of a hybrid nanoliquid including CNTs boosts more significantly. Next Sneha et al. [25] assessed the role of Marangoni convection and magnetic flux over the SWCNTs-based nanofluid flow in a porous media. It was shown that the rate of heat transfer was enhanced by Marangoni convection. El Glılı and Drıouıch [26] examined the non-orthogonal stagnation flow of CNTs-nano liquid along an extending surface. Their numerical results showed that SWCNTs have a more substantial heat transmission rate compared with MWCNTs. Dero et al. [27] obtained multiple solutions of CNTs-based nanofluid flow on a permeable shrinking sheet with viscous dissipation effects. They discovered that the velocity distribution for the nanoparticles volume fraction for MWCNTs-kerosene oil is larger than that of SWCNTs-kerosene oil nanoliquid. More work on the CNTs-based nanofluid along different geometries is available in previous studies [28,29,30,31].

The reviewed literature revealed extensive studies on nanofluid flows, but only few of them have focused on Marangoni convection in such systems. Notably, none of the existing works, including those cited, provide a numerical solution for aqueous-based nanofluid flow with MWCNTs under Marangoni convection while incorporating variable thermal conductivity and viscosity. The proposed model’s uniqueness lies in considering the influence of MWCNTs length and radius, analyzed using the Yamada–Ota (Y–O) thermal conductivity model. Numerical solutions are obtained through the bvp4c technique, with results presented through illustrations and tables. Table 1 highlights the novelty of the proposed model by highlighting gap analysis.

Table 1

Gap analysis

Published works Marangoni convection Y–O thermal conductivity model Aqueous-based nanofluid Variable viscosity and thermal conductivity
[2] × × ×
[3] × ×
[14] × ×
[16] × ×
Present

√ (Effect present), × (Effect absent).

2 Mathematical formulation

The proposed mathematical model is developed based on the following assumptions:

  • We consider an incompressible, aqueous-based nanofluid flow containing MWCNTs, following the Y–O thermal conductivity model.

  • The Y–O thermal conductivity model is engaged to gauge the radius and length of nanoparticles in a nanofluid flow.

  • The assumed nanofluid combination is in thermal balance and there is no slip between customary liquid and the assumed nanoparticles.

  • Temperature-dependent thermophysical characteristics are also considered to enhance the novelty of the model.

  • The flow is subjected to Marangoni convection.

  • A cartesian coordinate system ( x , y ) , where x is determined along the interface s and y is normal to it. The whole scenario is shown in Figure 1.

  • Moreover, T o is the temperature considered far off from the surface and σ is the surface tension, with expression σ = σ o + ϒ T 2 ( T T o ) , where ϒ T is the temperature coefficient of surface tension [14].

The model equations [13] following the above assumptions are

(1) u x + v y = 0 ,

(2) u u x + v u y = 1 ρ NF μ NF ( T ) y u y + μ NF ( T ) 2 u y 2 ,

(3) u T x + v T y = 1 ( ρ C p ) NF K NF ( T ) y T y + K NF ( T ) 2 T y 2 ,

where μ NF ( T ) is the variable viscosity given as [32]

(4) μ NF ( T ) = μ NF θ R θ R θ ,

K NF ( T ) is the variable thermal conductivity and θ R is the variable viscosity parameter, given by [32]

(5) K NF ( T ) = K NF ( 1 + E θ ) ,

where E is the variable thermal conductivity parameter. The appropriate boundary constraints [3] are

(6) μ NF ( T ) u y = σ T T x , v = 0 , T = T w = T o + A x m + 1 , at y = 0 , u 0 , T T o , as y .

where m is the exponent constant. Similarity transformations are appended as follows:

(7) ζ = C 2 x m 1 3 y , ψ = C 1 x 2 + m 3 f ( ζ ) , u = C 1 C 2 x 2 m + 1 3 f ( ζ ) , T = T o + A x m + 1 θ ( ζ ) , v = C 1 2 + m 3 x m 1 3 f ( ζ ) + C 1 x m 1 3 m 1 3 ζ f ( ζ ) .

Figure 1 
               Geometry of the flow model.
Figure 1

Geometry of the flow model.

By substituting Eqs. (1))–(6) in Eq. (7), we get the subsequent system

(8) μ NF ρ NF θ R θ R θ f + θ R [ θ R θ ] 2 θ f + 2 + m 3 f f 2 m + 1 3 f 2 = 0 ,

(9) K NF ( ρ C p ) NF ( E θ 2 + ( 1 + E θ ) ) θ Pr ( m + 1 ) f θ 2 + m 3 f θ = 0 ,

The constraints on the surface are outlined below

(10) f ( 0 ) = 0 , f ( 0 ) = μ NF μ F θ R θ θ R M A , θ ( 0 ) = 1 , f ( ) 0 , θ ( ) 0 .

The above equations include the following dimensionless parameters:

(11) M A = A ϒ T ( m + 1 ) μ F C 1 C 2 2 , Pr = μ f ( C p ) f k f ,

where M A is the Marangoni parameter and Pr is the Prandtl number. The base liquid and MWCNTs thermos-physical attributes are shown in Table 2.

Table 2

Thermophysical values for base fluid and nanoparticles [33]

Base fluid/nanoparticles k ρ C p
H2O 0.613 997.1 4,179
MWCNTs 3,000 1,600 796

2.1 Thermophysical models for nanofluid

The subsequent models are considered to reflect the thermophysical characteristics [34,35]

Dynamic viscosity

μ NF = μ F / ( 1 Φ s 1 ) 2.5 .

Density

ρ NF ρ F = ( 1 Φ s 1 ) + Φ s 1 ρ s 1 ρ F .

Specific heat capacity

( ρ C p ) NF ( ρ C p ) F = ( 1 Φ s 1 ) + Φ s 1 ( ρ C p ) s 1 ( ρ C p ) F .

Thermal conductivity

(12) k NF k F = 1 + k F k s 1 L R Φ s 1 0.2 + 1 k F k s 1 Φ s 1 L R Φ s 1 0.2 + 2 Φ s 1 k s 1 k s 1 k F ln k s 1 + k F 2 k s 1 1 Φ s 1 + 2 Φ s 1 k F k s 1 k F ln k s 1 + k F 2 k F ,

where L represents the length, R is the radius of CNTs, ( k F ) is the thermal conductivity of fluid, ( k s 1 ) is the thermal conductivity of CNTs, and ( Φ s 1 ) is the particle volume fraction of CNTs.

2.2 Engineering quantity

The rate of heat flux in non-dimensional notation is given as

Nu = x ( q c ) K F ( T w T o ) ,

where

(13) q c = K NF ( T ) T y .

The Nusselt number in dimensionless form is

(14) Nu = k NF k F ( 1 + E θ ) C 2 x m + 2 3 θ ( ζ ) .

2.3 Numerical scheme

We employed the collocation technique bvp4c, an efficient computational method (a finite difference mechanism) for solving sets of coupled differential equations. Using this technique, which is elaborated below, the order of the equations governing the flow is reduced to one. A technique to transform the higher-order system to one is appended as under

y 1 f , y 2 f , y 3 f , y y 1 f , y 4 θ , y 5 θ , y y 2 θ ,

y y 1 = 1 μ NF ρ NF θ R θ R y 4 θ R [ θ R y 4 ] 2 y 5 y 3 2 + m 3 y 1 y 3 + 2 m + 1 3 ( y 2 ) 2 ,

y y 2 = 1 K NF ( ρ C p ) NF ( E y 5 2 + ( 1 + E y 4 ) ) Pr ( m + 1 ) y 2 y 4 2 + m 3 y 1 y 5 ,

(15) y 1 ( 0 ) , y 3 ( 0 ) + μ NF μ F θ R y 4 θ R M A , y 4 ( 0 ) 1 , y 2 ( ) 0 , y 4 ( ) 0 .

A step size of Δ η = 0.01 with η = 14 is considered for all profiles for convergence purposes. Another important ingredient of this numerical scheme is the consideration of the initial guess that must satisfy the assumed boundary conditions. A complete flow diagram of the numerical technique (bvp4c) is given in Figure 2.

Figure 2 
                  Flow chart of bvp4c numerical scheme.
Figure 2

Flow chart of bvp4c numerical scheme.

This technique possesses the following advantages over the contemporary analytical and numerical schemes:

  • The bvp4c automatically adapts the mesh size for enhanced accuracy. By doing so we can achieve accuracy without user interference.

  • In comparison to the traditional techniques, bvp4c efficiently handles the discontinuities.

  • Contrary to analytical techniques, it handles nonlinear boundary value problems effectively.

  • The bvp4c is simpler and straightforward to implement with less effort.

  • In comparison to finite difference, the bvp4c method engages the collocation method leading to higher accuracy.

Table 3 portrays the numerical results of heat transmission rate for different estimations of variable thermal conductivity parameter ( E ) , Marangoni parameter ( M A ) , and exponent constant ( m ) . It is noted that the heat transmission rate improves for mounting estimations of all three parameters. This is because the improved values of variable thermal conductivity parameter result in enhanced diffusion that eventually triggers the heat transfer rate. Similarly, stronger thermocapillary convection is observed for increased Marangoni parameter that improves the fluid movement at the surface. Thus, enhanced heat transfer efficacy at the fluid interface is seen. Likewise, a large value of exponent triggers the thermal conductivity in high temperatures that eventually leads to an efficient heat transfer rate.

Table 3

Numerical estimates of Nusselt number Nu against varied parameters

( E ) ( M A ) ( m ) ( Nu )
0.1 0.10 0.1 0.84463
0.2 0.90408
0.3 0.95993
0.4 1.01320
0.12 0.89717
0.13 0.92157
0.14 0.94454
0.2 0.86797
0.3 0.89031
0.4 0.91179

The variation in thermal conductivity with particle volume fractions ranging from 1 to 6% is shown in Table 4, with the MCNTs length held constant. As the particle volume fraction of MCNTs increases at a fixed length, thermal conductivity is observed to rise. This finding aligns well with the experimental results stated by Assael et al. [35].

Table 4

Numerical estimates of thermal conductivity against higher values of particle volume fraction

Φ NF k NF
0.01 1.344517154332972
0.02 1.830524151848501
0.03 2.498053161467530
0.04 3.440034976804251
0.05 4.841142856090153
0.06 7.111676268043935

3 Graphical results analysis

This section aims to inspect the consequences of emerging factors on various flow and thermal distributions.

From Figure 3, the impact of length ( L ) of MWCNTs on the thermal conductivity can be observed. It is seen that for raising the length of MWCNTs, the thermal conductivity of nanofluid flow is significantly enhanced. This enhancement in the thermal conductivity may be owing to the high intrinsic thermal conductivity of MWCNTs, their aspect ratio, and uniform dispersion of MWCNTs in water, and volume fraction.

Figure 3 
               Upshot of length of MWCNTs on the thermal conductivity 
                     
                        
                        
                           
                              
                                 K
                              
                              
                                 NF
                              
                           
                        
                        {K}_{\text{NF}}
                     
                  .
Figure 3

Upshot of length of MWCNTs on the thermal conductivity K NF .

The role of the radius ( R ) of MWCNTs over the thermal conductivity can be visualized in Figure 4. It is inferred that for improving the radius of MWCNTs, the thermal conductivity drops. Here the factors affecting thermal conductivity may include increased inter-layer scattering, thermal boundary resistance, and contact resistance in composites.

Figure 4 
               Upshot of radius of MWCNTs on the thermal conductivity 
                     
                        
                        
                           
                              
                                 K
                              
                              
                                 NF
                              
                           
                        
                        {K}_{\text{NF}}
                     
                  .
Figure 4

Upshot of radius of MWCNTs on the thermal conductivity K NF .

The role of the Marangoni parameter ( M A ) on the velocity distribution f ( ζ ) can be seen in Figure 5 . It is seen that velocity distribution escalates for mounting values of ( M A ) . Marangoni convection causes the fluid to flow at higher velocities in areas with a noticeable surface tension noted. When there are significant temperature differences close to borders or interfaces, this impact may be very noticeable.

Figure 5 
               Upshot of Marangoni parameter 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       M
                                    
                                    
                                       A
                                    
                                 
                              
                              )
                           
                        
                        ({M}_{A})
                     
                   on 
                     
                        
                        
                           
                              f
                              ′
                           
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        f^{\prime} (\zeta )
                     
                  .
Figure 5

Upshot of Marangoni parameter ( M A ) on f ( ζ ) .

The trend of velocity distribution f ( ζ ) against the particle volume fraction ( Φ ) , is shown in Figure 6. It is seen that raising counts of particle volume fraction ( Φ ) upsurges the velocity distribution. Physically, strengthening the particle volume fraction causes a rapid collision of nanoparticles, which improves the velocity distribution.

Figure 6 
               Upshot of particle volume fraction 
                     
                        
                        
                           
                              (
                              
                                 Φ
                              
                              )
                           
                        
                        (\Phi )
                     
                   on 
                     
                        
                        
                           
                              f
                              ′
                           
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        f^{\prime} (\zeta )
                     
                  .
Figure 6

Upshot of particle volume fraction ( Φ ) on f ( ζ ) .

Figure 7 reveals the impact of the variable viscosity parameter ( θ R ) on the velocity distribution f ( ζ ) . It is seen that velocity distribution escalates by greater values of ( θ R ) . This is due to the fact that rising ( θ R ) implies zero shear rate viscosity, which lowers viscous drag force and speeds up flow.

Figure 7 
               Upshot of variable viscosity parameter 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       θ
                                    
                                    
                                       R
                                    
                                 
                              
                              )
                           
                        
                        ({\theta }_{R})
                     
                   on 
                     
                        
                        
                           
                              f
                              ′
                           
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        f^{\prime} (\zeta )
                     
                  .
Figure 7

Upshot of variable viscosity parameter ( θ R ) on f ( ζ ) .

The influence of the exponent constant ( m ) on the velocity distribution f ( ζ ) is displayed in Figure 8. It is seen that velocity distribution drops significantly for greater counts of ( m ) . When m rises, the temperature variation throughout the surface reduces, resulting in a lesser velocity distribution.

Figure 8 
               Upshot of exponent constant 
                     
                        
                        
                           
                              (
                              
                                 m
                              
                              )
                           
                        
                        (m)
                     
                   on 
                     
                        
                        
                           
                              f
                              ′
                           
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        f^{\prime} (\zeta )
                     
                  .
Figure 8

Upshot of exponent constant ( m ) on f ( ζ ) .

Figure 9 is portrayed to analyze the consequences of the Marangoni parameter ( M A ) on the temperature distribution θ ( ζ ) . It is seen that temperature distribution dwindled for the greater counts of ( M A ) . The higher estimates of ( M A ) enhance the surface tension resulting in improved heat transfer, thus reducing the temperature gradient with the assumed fluid flow.

Figure 9 
               Upshot of Marangoni parameter 
                     
                        
                        
                           
                              (
                              
                                 
                                    
                                       M
                                    
                                    
                                       A
                                    
                                 
                              
                              )
                           
                        
                        ({M}_{A})
                     
                   on 
                     
                        
                        
                           θ
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        \theta (\zeta )
                     
                  .
Figure 9

Upshot of Marangoni parameter ( M A ) on θ ( ζ ) .

The role of particle volume fraction ( Φ ) on the temperature distribution is displayed in Figure 10. It is revealed that the temperature distribution of nanofluid containing MWCNTs nanoparticles escalates for mounting ( Φ ) . Adding more MWCNTs nanoparticles to the fluid enhances heat transfer capabilities owing to the high thermal conductivity of MWCNTs. This results in an improved heat transfer rate.

Figure 10 
               Upshot of particle volume fraction 
                     
                        
                        
                           
                              (
                              
                                 Φ
                              
                              )
                           
                        
                        (\Phi )
                     
                   on 
                     
                        
                        
                           θ
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        \theta (\zeta )
                     
                  .
Figure 10

Upshot of particle volume fraction ( Φ ) on θ ( ζ ) .

Finally, the role of variable thermal conductivity ( E ) on the temperature distribution can be analyzed from Figure 11. It is revealed that temperature distribution escalates for greater counts of ( E ) . This is because the temperature is increased as a consequence of the quicker heat transfer process from sheet to fluid.

Figure 11 
               Upshot of variable thermal conductivity 
                     
                        
                        
                           
                              (
                              
                                 E
                              
                              )
                           
                        
                        (E)
                     
                   on 
                     
                        
                        
                           θ
                           
                              
                                 (
                                 
                                    ζ
                                 
                                 )
                              
                           
                        
                        \theta (\zeta )
                     
                  .
Figure 11

Upshot of variable thermal conductivity ( E ) on θ ( ζ ) .

To validate the outcomes of the proposed mode, a comparison with Montgomery and Wang [12] when m = 0 is made. An excellent correlation is achieved (Figure 12).

Figure 12 
               Comparison with published work [12] when 
                     
                        
                        
                           m
                           =
                           0
                        
                        m=0
                     
                  .
Figure 12

Comparison with published work [12] when m = 0 .

4 Conclusion

This study examined the thermal properties of aqueous nanofluid flow containing immersed MWCNTs influenced by Marangoni convection. Temperature-dependent models for viscosity and thermal conductivity were utilized instead of constant values. The Y–O model was applied to analyze the impact of radius and length on nanofluid thermal conductivity. Numerical results were obtained using the bvp4c method. Key findings are summarized below:

  • Increasing the length of MWCNTs improves the nanofluid’s thermal conductivity, whereas increasing the radius reduces it.

  • The thermal conductivity is improved in order to increase the particle volume fraction while maintaining the fixed length of MWCNTs.

  • The Marangoni parameter boosts the velocity distribution while diminishing the temperature distribution.

  • Velocity distribution escalates by raising the viscosity parameter and particle volume fraction.

  • It is seen that velocity distribution drops significantly for greater counts of exponent constant.

  • By raising the variable temperature parameter, the temperature distribution escalates.

  • Marangoni convection drops the temperature distribution but enhances the heat transmission rate.

Acknowledgments

This work was supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-DDRSP2503).

  1. Funding information: This work was supported and funded by the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University (IMSIU) (grant number IMSIU-DDRSP2503).

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Appendix

u x = C 1 C 2 2 m + 1 3 x 2 m 2 3 f + C 1 C 2 ζ m 1 3 x 2 m 2 3 f ,

v y = C 1 C 2 2 m + 1 3 x 2 m 2 3 f C 1 C 2 ζ m 1 3 x 2 m 2 3 f ,

u y = C 1 C 2 2 x m f ,

2 u y 2 = C 1 C 2 3 x 4 m 1 3 f ,

T x = A ( m + 1 ) x m θ + A x m ζ θ ,

T y = A C 2 x 4 m + 2 3 θ ,

2 T y 2 = A C 2 2 x 3 m 1 3 θ ,

μ NF ( T ) y = μ NF θ R ( θ R θ ) 2 C 2 x m 1 3 θ ,

(A1) K NF ( T ) y = K NF E C 2 x m 1 3 θ .

Using the above values of partial derivative in model equations, the set of reduced differential equations will be obtained.

References

[1] Abu-Nada E. Effects of variable viscosity and thermal conductivity of Al2O3–water nanofluid on heat transfer enhancement in natural convection. Int J Heat Fluid Flow. 2009;30(4):679–90. 10.1016/j.ijheatfluidflow.2009.02.003.Search in Google Scholar

[2] Mahmoud MA, Megahed AM. MHD flow and heat transfer in a non-Newtonian liquid film over an unsteady stretching sheet with variable fluid properties. Can J Phys. 2009;87(10):1065–71. 10.1139/P09-066.Search in Google Scholar

[3] Abu-Nada E, Masoud Z, Oztop HF, Campo A. Effect of nanofluid variable properties on natural convection in enclosures. Int J Therm Sci. 2010;49(3):479–91. 10.1016/j.ijthermalsci.2009.09.002.Search in Google Scholar

[4] Roslan R, Saleh H, Hashim I. Buoyancy-driven heat transfer in nanofluid-filled trapezoidal enclosure with variable thermal conductivity and viscosity. Numer Heat Transfer, Part A: Appl. 2011;60(10):867–82. 10.1080/10407782.2011.616778.Search in Google Scholar

[5] Sebdani SM, Mahmoodi M, Hashemi SM. Effect of nanofluid variable properties on mixed convection in a square cavity. Int J Therm Sci. 2012;52:112–26. 10.1016/j.ijthermalsci.2011.09.003.Search in Google Scholar

[6] Shamshuddin MD, Eid MR. Magnetized nanofluid flow of ferromagnetic nanoparticles from parallel stretchable rotating disk with variable viscosity and thermal conductivity. Chin J Phys. 2021;74:20–37. 10.1016/j.cjph.2021.07.038.Search in Google Scholar

[7] Alharbi KAM, Shahmir N, Ramzan M, Almusawa MY, Kadry S. Bioconvective radiative unsteady Casson nanofluid flow across two concentric stretching cylinders with variable viscosity and variable thermal conductivity. Numer Heat Transfer, Part A: Appl. 2024;85(10):1653–70. 10.1080/10407782.2023.2208733.Search in Google Scholar

[8] Mandal G, Pal D. Dual solutions of radiative Ag-MoS2/water hybrid nanofluid flow with variable viscosity and variable thermal conductivity along an exponentially shrinking permeable Riga surface: Stability and entropy generation analysis. Int J Model Simul. 2023;44(6):360–85. 10.1080/02286203.2023.2171656.Search in Google Scholar

[9] Govindarajulu K, Subramanyam Reddy A, Jagadeshkumar K, Srinivas S, Kumar BR, Vajravelu K. Entropy generation on MHD pulsatile flow of third grade hybrid nanofluid in a vertical porous channel with nonuniform heat source/sink, variable viscosity, thermal conductivity, and Joule heating: a numerical study. Numer Heat Transfer, Part A: Appl. 2023;85(24):4184–203. 10.1080/10407782.2023.2255929.Search in Google Scholar

[10] Babu MS, Sankar GR, Velpula VR, Chu YM, Khan MI, Raju CSK, et al. Chemically reactive flow of viscous thermophoretic fluid over wedge with variable thermal conductivity and viscosity. Case Stud Therm Eng. 2023;45:102924. 10.1016/j.csite.2023.102924.Search in Google Scholar

[11] Golia C, Viviani A. Non isobaric boundary layers related to Marangoni flows. Meccanica. 1986;21:200–4.10.1007/BF01556486Search in Google Scholar

[12] Montgomery CD, Wang B. Prandtl number effects for Marangoni convection over a flat surface. Int J Therm Sci. 2001;40(6):564–70. 10.1016/S1290-0729(01)01244-3.Search in Google Scholar

[13] Arifin NM, Nazar R, Pop I. Similarity solution of Marangoni convection boundary layer flow over a flat surface in a nanofluid. J Appl Math. 2013;2013:634746. 10.1155/2013/634746.Search in Google Scholar

[14] Hayat T, Khan MI, Farooq M, Alsaedi A, Yasmeen T. Impact of Marangoni convection in the flow of carbon–water nanofluid with thermal radiation. Int J Heat Mass Transf. 2017;106:810–5. 10.1016/j.ijheatmasstransfer.2016.08.115.Search in Google Scholar

[15] Qayyum S. Dynamics of Marangoni convection in hybrid nanofluid flow submerged in ethylene glycol and water base fluids. Int Commun Heat Mass Transfer. 2020;119:104962. 10.1016/j.icheatmasstransfer.2020.104962.Search in Google Scholar

[16] Abdullah NN, Som ANM, Arifin NM, Ab Ghani A. The effect of MHD on Marangoni boundary layer of hybrid nanofluid flow past a permeable stretching surface. CFD Letters. 2023;15(5):65–73. 10.37934/cfdl.15.5.6573.Search in Google Scholar

[17] Anusha T, Pérez LM, Mahabaleshwar US, Zeidan D. An MHD nanofluid flow with Marangoni laminar boundary layer over a porous medium with heat and mass transfer. Int J Model Simul. 2024;1–17. 10.1080/02286203.2024.2335112.Search in Google Scholar

[18] Abbas M, Khan N, Hashmi MS, Alhefthi RK, Rezapour S, Inc M. Thermal Marangoni convection in two-phase quadratic convective flow of dusty MHD trihybrid nanofluid with non-linear heat source. Case Stud Therm Eng. 2024;57:104190. 10.1016/j.csite.2024.104190.Search in Google Scholar

[19] Abbas M, Khan N, Shehzad SA. Numerical analysis of Marangoni convected dusty second-grade nanofluid flow in a suspension of chemically reactive microorganisms. Proc Inst Mech Eng, Part C-J Mech Eng Sci. 2024;238(10):4400–17. 10.1177/09544062231209828.Search in Google Scholar

[20] Rehman A, Khun MC, Alsubaie AS, Inc M. Influence of Marangoni convection, viscous dissipation, and variable fluid viscosity of nanofluid flow on stretching surface analytical analysis. ZAMM - J Appl Math Mech. 2024;104(3):e202300413. 10.1002/zamm.202300413.Search in Google Scholar

[21] Maxwell JC. A treatise on electricity and magnetism. Vol. 1, Oxford: Clarendon Press; 1873.Search in Google Scholar

[22] Choi SU, Eastman JA. Enhancing thermal conductivity of fluids with nanoparticles. San Francisco, USA: ASME International Mechanical Engineering Congress & Exposition; 1995.Search in Google Scholar

[23] Yan SR, Toghraie D, Abdulkareem LA, Alizadeh AA, Barnoon P, Afrand M. The rheological behavior of MWCNTs–ZnO/Water–Ethylene glycol hybrid non-Newtonian nanofluid by using of an experimental investigation. J Mater Res Technol. 2020;9(4):8401–6.10.1016/j.jmrt.2020.05.018Search in Google Scholar

[24] Harish Babu D, Naidu KK, Deo S, Satya Narayana PV. Impacts of inclined Lorentz forces on hybrid CNTs over an exponentially stretching sheet with slip flow. Int J Model Simul. 2023;43(3):310–24. 10.1080/02286203.2022.2079109.Search in Google Scholar

[25] Sneha KN, Bognar G, Mahabaleshwar US, Singh DK, Singh OP. Magnetohydrodynamics effect of Marangoni nano boundary layer flow and heat transfer with CNT and radiation. J Magn Magn Mater. 2023;575:170721. 10.1016/j.jmmm.2023.170721.Search in Google Scholar

[26] El Glılı I, Drıouıch M. Impact of inclined magnetic field on non-orthogonal stagnation point flow of CNT-water through stretching surface in a porous medium. J Therm Eng. 2024;10(1):115–29. 10.18186/thermal.1429409.Search in Google Scholar

[27] Dero S, Fadhel MA, Lund LA, Shah NA. Multiple solutions of unsteady flow of CNTs nanofluid over permeable shrinking surface with effects of dissipation and slip conditions. Mod Phys Lett B. 2024;38(15):2450120. 10.1142/S0217984924501203.Search in Google Scholar

[28] Bakhtiyar NK, Esmaeili S, Javadpour R, Heris SZ. Experimental investigation of indirect heat transfer through a novel designed lab-scale setup using functionalized MWCNTs nanofluids (MWCNTs-COOH/water and MWCNTs-OH/water). Case Stud Therm Eng. 2023;45:102951. 10.1016/j.csite.2023.102951.Search in Google Scholar

[29] Pourpasha H, Heris SZ, Mousavi SB. Thermal performance of novel ZnFe2O4 and TiO2-doped MWCNT nanocomposites in transformer oil. J Mol Liq. 2024;394:123727. 10.1016/j.molliq.2023.123727.Search in Google Scholar

[30] Farooq U, Liu T, Farooq U, Majeed S. Non-similar analysis of bioconvection MHD micropolar nanofluid on a stretching sheet with the influences of Soret and Dufour effects. Appl Water Sci. 2024;14(6):116. 10.1007/s13201-024-02143-0.Search in Google Scholar

[31] Liu Z, Ali AB, Hussein RA, Singh NSS, Al-Bahrani M, Abdullaeva B, et al. Using evolutionary algorithms and group method of data handling ANN for prediction of the viscosity MWCNT-ZnO/oil SAE 50 nano-lubricant. Int Commun Heat Mass Transf. 2025;163:108749. 10.1016/j.icheatmasstransfer.2025.108749.Search in Google Scholar

[32] Swain K, Parida SK, Dash GC. MHD heat and mass transfer on stretching sheet with variable fluid properties in porous medium. AMSE J Model B. 2017;86:706–26.10.18280/mmc_b.860307Search in Google Scholar

[33] Muhammad S, Ali G, Shah Z, Islam S, Hussain SA. The rotating flow of magneto hydrodynamic carbon nanotubes over a stretching sheet with the impact of non-linear thermal radiation and heat generation/absorption. Appl Sci. 2018;8(4):482. 10.3390/app8040482.Search in Google Scholar

[34] Khashi’ie NS, Arifin NM, Hafidzuddin EH, Wahi N, Pop I. Mixed convective stagnation point flow of a thermally stratified hybrid Cu-Al2O3/water nanofluid over a permeable stretching/shrinking sheet. ASM Sci J - Akademi Sains Malays. 2019;12:17–25. 10.1038/s41598-023-34686-8. for ICoAIMS2019; Alharbi KAM, Ramzan M, Shahmir N, Ghazwani HAS, Elmasry Y, Eldin SM, et al. Comparative appraisal of mono and hybrid nanofluid flows comprising carbon nanotubes over a three-dimensional surface impacted by Cattaneo–Christov heat flux. Sci Report. 2023;13(1):7964.Search in Google Scholar PubMed PubMed Central

[35] Assael MJ, Chen CF, Metaxa I, Wakeham, WA. Thermal conductivity of suspensions of carbon nanotubes in water. Int J Thermophys. 2004;25:971–85. 10.1023/B:IJOT.0000038494.22494.04.Search in Google Scholar

Received: 2025-01-11
Revised: 2025-04-06
Accepted: 2025-05-05
Published Online: 2025-09-20

© 2025 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Research Articles
  2. Single-step fabrication of Ag2S/poly-2-mercaptoaniline nanoribbon photocathodes for green hydrogen generation from artificial and natural red-sea water
  3. Abundant new interaction solutions and nonlinear dynamics for the (3+1)-dimensional Hirota–Satsuma–Ito-like equation
  4. A novel gold and SiO2 material based planar 5-element high HPBW end-fire antenna array for 300 GHz applications
  5. Explicit exact solutions and bifurcation analysis for the mZK equation with truncated M-fractional derivatives utilizing two reliable methods
  6. Optical and laser damage resistance: Role of periodic cylindrical surfaces
  7. Numerical study of flow and heat transfer in the air-side metal foam partially filled channels of panel-type radiator under forced convection
  8. Water-based hybrid nanofluid flow containing CNT nanoparticles over an extending surface with velocity slips, thermal convective, and zero-mass flux conditions
  9. Dynamical wave structures for some diffusion--reaction equations with quadratic and quartic nonlinearities
  10. Solving an isotropic grey matter tumour model via a heat transfer equation
  11. Study on the penetration protection of a fiber-reinforced composite structure with CNTs/GFP clip STF/3DKevlar
  12. Influence of Hall current and acoustic pressure on nanostructured DPL thermoelastic plates under ramp heating in a double-temperature model
  13. Applications of the Belousov–Zhabotinsky reaction–diffusion system: Analytical and numerical approaches
  14. AC electroosmotic flow of Maxwell fluid in a pH-regulated parallel-plate silica nanochannel
  15. Interpreting optical effects with relativistic transformations adopting one-way synchronization to conserve simultaneity and space–time continuity
  16. Modeling and analysis of quantum communication channel in airborne platforms with boundary layer effects
  17. Theoretical and numerical investigation of a memristor system with a piecewise memductance under fractal–fractional derivatives
  18. Tuning the structure and electro-optical properties of α-Cr2O3 films by heat treatment/La doping for optoelectronic applications
  19. High-speed multi-spectral explosion temperature measurement using golden-section accelerated Pearson correlation algorithm
  20. Dynamic behavior and modulation instability of the generalized coupled fractional nonlinear Helmholtz equation with cubic–quintic term
  21. Study on the duration of laser-induced air plasma flash near thin film surface
  22. Exploring the dynamics of fractional-order nonlinear dispersive wave system through homotopy technique
  23. The mechanism of carbon monoxide fluorescence inside a femtosecond laser-induced plasma
  24. Numerical solution of a nonconstant coefficient advection diffusion equation in an irregular domain and analyses of numerical dispersion and dissipation
  25. Numerical examination of the chemically reactive MHD flow of hybrid nanofluids over a two-dimensional stretching surface with the Cattaneo–Christov model and slip conditions
  26. Impacts of sinusoidal heat flux and embraced heated rectangular cavity on natural convection within a square enclosure partially filled with porous medium and Casson-hybrid nanofluid
  27. Stability analysis of unsteady ternary nanofluid flow past a stretching/shrinking wedge
  28. Solitonic wave solutions of a Hamiltonian nonlinear atom chain model through the Hirota bilinear transformation method
  29. Bilinear form and soltion solutions for (3+1)-dimensional negative-order KdV-CBS equation
  30. Solitary chirp pulses and soliton control for variable coefficients cubic–quintic nonlinear Schrödinger equation in nonuniform management system
  31. Influence of decaying heat source and temperature-dependent thermal conductivity on photo-hydro-elasto semiconductor media
  32. Dissipative disorder optimization in the radiative thin film flow of partially ionized non-Newtonian hybrid nanofluid with second-order slip condition
  33. Bifurcation, chaotic behavior, and traveling wave solutions for the fractional (4+1)-dimensional Davey–Stewartson–Kadomtsev–Petviashvili model
  34. New investigation on soliton solutions of two nonlinear PDEs in mathematical physics with a dynamical property: Bifurcation analysis
  35. Mathematical analysis of nanoparticle type and volume fraction on heat transfer efficiency of nanofluids
  36. Creation of single-wing Lorenz-like attractors via a ten-ninths-degree term
  37. Optical soliton solutions, bifurcation analysis, chaotic behaviors of nonlinear Schrödinger equation and modulation instability in optical fiber
  38. Chaotic dynamics and some solutions for the (n + 1)-dimensional modified Zakharov–Kuznetsov equation in plasma physics
  39. Fractal formation and chaotic soliton phenomena in nonlinear conformable Heisenberg ferromagnetic spin chain equation
  40. Single-step fabrication of Mn(iv) oxide-Mn(ii) sulfide/poly-2-mercaptoaniline porous network nanocomposite for pseudo-supercapacitors and charge storage
  41. Novel constructed dynamical analytical solutions and conserved quantities of the new (2+1)-dimensional KdV model describing acoustic wave propagation
  42. Tavis–Cummings model in the presence of a deformed field and time-dependent coupling
  43. Spinning dynamics of stress-dependent viscosity of generalized Cross-nonlinear materials affected by gravitationally swirling disk
  44. Design and prediction of high optical density photovoltaic polymers using machine learning-DFT studies
  45. Robust control and preservation of quantum steering, nonlocality, and coherence in open atomic systems
  46. Coating thickness and process efficiency of reverse roll coating using a magnetized hybrid nanomaterial flow
  47. Dynamic analysis, circuit realization, and its synchronization of a new chaotic hyperjerk system
  48. Decoherence of steerability and coherence dynamics induced by nonlinear qubit–cavity interactions
  49. Finite element analysis of turbulent thermal enhancement in grooved channels with flat- and plus-shaped fins
  50. Modulational instability and associated ion-acoustic modulated envelope solitons in a quantum plasma having ion beams
  51. Statistical inference of constant-stress partially accelerated life tests under type II generalized hybrid censored data from Burr III distribution
  52. On solutions of the Dirac equation for 1D hydrogenic atoms or ions
  53. Entropy optimization for chemically reactive magnetized unsteady thin film hybrid nanofluid flow on inclined surface subject to nonlinear mixed convection and variable temperature
  54. Stability analysis, circuit simulation, and color image encryption of a novel four-dimensional hyperchaotic model with hidden and self-excited attractors
  55. A high-accuracy exponential time integration scheme for the Darcy–Forchheimer Williamson fluid flow with temperature-dependent conductivity
  56. Novel analysis of fractional regularized long-wave equation in plasma dynamics
  57. Development of a photoelectrode based on a bismuth(iii) oxyiodide/intercalated iodide-poly(1H-pyrrole) rough spherical nanocomposite for green hydrogen generation
  58. Investigation of solar radiation effects on the energy performance of the (Al2O3–CuO–Cu)/H2O ternary nanofluidic system through a convectively heated cylinder
  59. Quantum resources for a system of two atoms interacting with a deformed field in the presence of intensity-dependent coupling
  60. Studying bifurcations and chaotic dynamics in the generalized hyperelastic-rod wave equation through Hamiltonian mechanics
  61. A new numerical technique for the solution of time-fractional nonlinear Klein–Gordon equation involving Atangana–Baleanu derivative using cubic B-spline functions
  62. Interaction solutions of high-order breathers and lumps for a (3+1)-dimensional conformable fractional potential-YTSF-like model
  63. Hydraulic fracturing radioactive source tracing technology based on hydraulic fracturing tracing mechanics model
  64. Numerical solution and stability analysis of non-Newtonian hybrid nanofluid flow subject to exponential heat source/sink over a Riga sheet
  65. Numerical investigation of mixed convection and viscous dissipation in couple stress nanofluid flow: A merged Adomian decomposition method and Mohand transform
  66. Effectual quintic B-spline functions for solving the time fractional coupled Boussinesq–Burgers equation arising in shallow water waves
  67. Analysis of MHD hybrid nanofluid flow over cone and wedge with exponential and thermal heat source and activation energy
  68. Solitons and travelling waves structure for M-fractional Kairat-II equation using three explicit methods
  69. Impact of nanoparticle shapes on the heat transfer properties of Cu and CuO nanofluids flowing over a stretching surface with slip effects: A computational study
  70. Computational simulation of heat transfer and nanofluid flow for two-sided lid-driven square cavity under the influence of magnetic field
  71. Irreversibility analysis of a bioconvective two-phase nanofluid in a Maxwell (non-Newtonian) flow induced by a rotating disk with thermal radiation
  72. Hydrodynamic and sensitivity analysis of a polymeric calendering process for non-Newtonian fluids with temperature-dependent viscosity
  73. Exploring the peakon solitons molecules and solitary wave structure to the nonlinear damped Kortewege–de Vries equation through efficient technique
  74. Modeling and heat transfer analysis of magnetized hybrid micropolar blood-based nanofluid flow in Darcy–Forchheimer porous stenosis narrow arteries
  75. Activation energy and cross-diffusion effects on 3D rotating nanofluid flow in a Darcy–Forchheimer porous medium with radiation and convective heating
  76. Insights into chemical reactions occurring in generalized nanomaterials due to spinning surface with melting constraints
  77. Influence of a magnetic field on double-porosity photo-thermoelastic materials under Lord–Shulman theory
  78. Soliton-like solutions for a nonlinear doubly dispersive equation in an elastic Murnaghan's rod via Hirota's bilinear method
  79. Analytical and numerical investigation of exact wave patterns and chaotic dynamics in the extended improved Boussinesq equation
  80. Nonclassical correlation dynamics of Heisenberg XYZ states with (x, y)-spin--orbit interaction, x-magnetic field, and intrinsic decoherence effects
  81. Exact traveling wave and soliton solutions for chemotaxis model and (3+1)-dimensional Boiti–Leon–Manna–Pempinelli equation
  82. Unveiling the transformative role of samarium in ZnO: Exploring structural and optical modifications for advanced functional applications
  83. On the derivation of solitary wave solutions for the time-fractional Rosenau equation through two analytical techniques
  84. Analyzing the role of length and radius of MWCNTs in a nanofluid flow influenced by variable thermal conductivity and viscosity considering Marangoni convection
  85. Advanced mathematical analysis of heat and mass transfer in oscillatory micropolar bio-nanofluid flows via peristaltic waves and electroosmotic effects
  86. Review Article
  87. Examination of the gamma radiation shielding properties of different clay and sand materials in the Adrar region
  88. Special Issue on Fundamental Physics from Atoms to Cosmos - Part II
  89. Possible explanation for the neutron lifetime puzzle
  90. Special Issue on Nanomaterial utilization and structural optimization - Part III
  91. Numerical investigation on fluid-thermal-electric performance of a thermoelectric-integrated helically coiled tube heat exchanger for coal mine air cooling
  92. Special Issue on Nonlinear Dynamics and Chaos in Physical Systems
  93. Analysis of the fractional relativistic isothermal gas sphere with application to neutron stars
  94. Abundant wave symmetries in the (3+1)-dimensional Chafee–Infante equation through the Hirota bilinear transformation technique
  95. Successive midpoint method for fractional differential equations with nonlocal kernels: Error analysis, stability, and applications
Downloaded on 24.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2025-0160/html
Scroll to top button