Home Impact of permeability and fluid parameters in couple stress media on rotating eccentric spheres
Article Open Access

Impact of permeability and fluid parameters in couple stress media on rotating eccentric spheres

  • Amal Al-Hanaya and Shreen El-Sapa EMAIL logo
Published/Copyright: December 31, 2024

Abstract

This study elucidates the axisymmetric rotation of two eccentric spheres immersed in an incompressible couple stress fluid within a porous annular region. The investigation incorporates boundary conditions applied to the surfaces of both spheres, which rotate axially at differing angular velocities. Utilizing a collocation scheme alongside a semi-analytical approach under the assumption of low Reynolds numbers, this research delves into the hydrodynamic couple force (torque) exerted by the Brinkman couple stress fluid on the interior sphere. The analysis reveals that the dimensionless torque becomes increasingly significant as permeability enhances the size ratio, couple stress fluid parameters, and separation distance. The computed torque values due to stress and couple stress are presented graphically, providing visual clarity to the findings. Additionally, all the results align closely with the numerical analysis of two eccentric rigid spheres in conventional couple stress fluids, devoid of permeability effects or slip conditions, as demonstrated in previous studies (Al-Hanaya and El-Sapa (Effects of slippage and permeability of couple stress fluid squeezed between two concentric rotating spheres. Phys Fluids. 2023;35:103112); Al-Hanaya et al. (Axisymmetric motion of an incompressible couple stress fluid between two eccentric rotating spheres. J Appl Mech Tech Phys. 2022;63(5):1–9)).

Nomenclature

( r , θ , ϕ )

spherical polar coordinates

q

velocity vector of fluid (m/s)

U

characteristic velocity

ξ

couple stress parameter a 2 μ η

α

permeability parameter a 2 K

μ , λ

viscosity coefficient (Pa s)

η , η

couple stress viscosity coefficients (Pa s)

ρ

density of the fluid (kg/m 3 )

ω

spin vector

f

body force per unit mass (m/s2)

c

body couple per unit mass (m2/s2)

M

couple stress diadic

p

pressure (Pa (N/m2))

E

rate of deformation tensor

m

trace of couple stress tensor

Ω

characteristic angular velocity (1/s)

1 Introduction

For the past 50 years, researchers have concentrated on understanding fluid flow through porous materials. Numerous applications, including heat exchange, separation, filtration, increased soil contamination, shale gas extraction, oil recovery, and catalytic assistance in industry, require an understanding of flow behavior in porous media. Fluid flow in a porous material that is isotropic and homogenous on a large scale is governed by either Darcy’s law or Brinkman’s equations. Porosity has a negligible impact on Darcy’s law. However, because Darcy’s law cannot resolve an issue with high porosity and a high tangential rate, one runs into difficulties. This problem is solved by applying Brinkman’s equation [1]. Leontev [2] studied the Newtonian fluid flow past spheres and cylinders in a porous media under Navier boundary conditions. El-Sapa [3] proposed an analytical method for thermophoresis of particles in the center of a spherical cavity that was filled with Brinkman’s porous medium. The mobility of a spherical particle was investigated by Faltas et al. in the context of a semi-infinite Brinkman flow [4]. Madasu and Bucha [5] studied magnetohydrodynamics’ effect on micropolar fluid flow around a sphere in Brinkman’s porous medium.

Couple stress fluids are useful in real life. Bio-fluids, colloidal fluids, liquid crystals, and others are couple stress fluids. Polar effects in couple stress fluid flow were initially studied by Stokes [6]. Stokes [7] proposed a basic theory to explain fluid polar phenomena and generalize classical theory. Stokes discovered that the polar effect in fluids is a result of couple stresses and body couplings. He also observed that the mechanical action is analogous to both forces and momentum distributions. Couple stresses lead to a size-dependent effect that is not accounted for in classical hydrodynamics. However, it is essential in specific cases, such as fluid flows with suspended particles, colloidal fluids, or liquid crystals. Many studies have examined the utilization of these fluids. Jangili et al. [8] examined how thermal conductivity and viscosity affect couple stress fluid flow irreversibility between asymmetrically heated parallel plates. Al-Hanaya et al. [9] researched the axisymmetric movement of fluid couple stress between two eccentric spinning spheres. They found that the solid sphere’s normalized torque rotates inside the outer sphere without slipping. The normalized torque rises with the size ratio and separation distance. Couple stress fluids have been utilized in many porous media, such as in the studies conducted by El-Sapa and Almoneef [10], focusing on the axisymmetric flow of aerosol particles under a slippage regime.

Many natural applications exist for the flow motion theory of two eccentric spheres in micropolar fluids and rheological fluids. Fluid microstructure can affect limited and asymmetric flow inflow performance. Majumdar investigated the asymmetrical homogenous incompressible viscous liquid flow in the annulus between two rigid spheres [11]. Kamel et al. [12] studied micropolar liquid flow around two rolling rigid spheres. Munson [13] studied eccentric Newtonian fluid flow. Recently, Sherief et al. [14] used semi-analytical methods to solve the no-slip problem of pushing a sphere into a non-concentric spherical hollow that was filled with creeping micropolar fluid. El-Sapa [15] examined the relationship between a spherical envelope with a slip regime and a non-concentric rigid sphere submerged in a micropolar fluid. Al-Hanaya et al. [9] considered a slow, steady motion of a couple stress fluid between two non-concentric, rotating spherical boundaries.

The torque generated by the rotation of solid particles in incompressible Newtonian, non-Newtonian, or porous media has applications in the domains of chemistry, biomedicine, environmental science, and other scientific disciplines. Scientists have analyzed the occurrence of torque in the movement of spheres, cylinders, and other geometric shapes in a couple of stress fluids. Davis [16] studied the force and torque on a rotating sphere close to and within a fluid-filled rotating sphere. Madasu and Sarkar [17] conducted an assessment of the torque exerted on the slow rotation of a slightly deformed slip sphere in a Brinkman medium. In the context of the Theory of Porous Media, this work introduces a two-scale Mesh-in-Element (MIEL) approach that allows for finite element discretization and significant scale coupling. It investigates the theoretical underpinnings, algorithmic implementation, and application of MIEL for modeling heterogeneous, inelastic materials, building on prior research and verifying its efficacy using numerical examples by Maike et al. [18]. Also, El-Sapa and Al-Hanaya [19] examined how variations in slip conditions at the fluid interfaces and the permeability of the medium affect the flow characteristics by employing mathematical modeling and analytical techniques in two concentric spheres with couple stress fluid in the annulus. Further studies [20,21] explore fluid dynamics within annular geometries, focusing on peristaltic motion and magnetohydrodynamic effects. Their initial study formulates a mathematical model for blood flow in a porous annular space between deformable tubes, providing valuable insights for peristaltic endoscopes. Then, they develop a mathematical model for creeping electro-magnetohydrodynamic peristaltic propulsion in an annular gap created by sinusoidally deforming tubes. Together, these studies contribute to a deeper theoretical understanding of fluid dynamics in both medical and engineering applications.

According to the aforementioned literature review, this study looks into the Brinkman couple stress fluid between two spinning eccentric spheres, which has many potential applications, such as tribology, which is understanding the lubrication behavior of non-Newtonian fluids in eccentric geometries; also biomechanics such as modeling the flow of blood or other biological fluids in vessels with irregular shapes; and finally, chemical engineering, such as analyzing the mixing and transport of non-Newtonian fluids in industrial processes. The solutions are obtained with analytical and numerical approaches, and there are other numerical methods, such as finite element or finite difference methods, that can be employed to solve the governing equations.

2 Field equations and constitutive relations

In 1984, Stokes [6] extended the well-known Navier–Stokes theory to develop the micro-continuum theory of couple stress fluid. The Stokes theory takes into account the impact of particle size as well as the presence of polar effects in the fluid medium, such as couple stress, body couple stress, and asymmetric stress. Previous studies [18,19] provide the revised equations, which regulate the continuity and momentum of an incompressible fluid with couple stress in Brinkman’s medium as

(2.1) q = 0 ,

(2.2) p + μ q + η q + μ K q = 0 ,

where q is the volume-averaged velocity, p is the pore average pressure, μ represents the viscosity of the fluid, η is the first couple stress viscosity coefficient, and η is the second couple stress viscosity coefficient. If the couple stress coefficient η is taken to be zero, then the equation of motion (2.2) reduces to the classical Navier–Stokes’s equation. The fluid velocity via a porous medium is governed by the following Darcy–Brinkman equation, but we assume that μ e is the effective viscosity when the fluid flow takes place in the porous medium, which is equal to μ that is suggested by Ochoa-Tapia and Whitaker [24,25]. In addition, K is the Darcy permeability of the permeable medium and also is a scalar for isotropic porous medium. Otherwise, K is a second-order tensor [26] and the Brinkman equation of couple stress fluid reduces to the Darcy equation when K 0 and to the Stokes equation of couple stress fluid when K .

The constitutive equations of the stress dyadic Π and couple stress dyadic M of the couple stress fluids are given by [7,9]

(2.3) Π = p I + λ q + 2 μ E + 1 2 I M ,

(2.4) M = m I + 4 η ω + 4 η ( ω ) T ,

where E = 1 2 ( q + ( q ) T ) is the deformation dyadic, I is the unit dyadic, and ω is the vorticity vector. Here, ( ) T denotes the transpose of a dyadic. The viscosity coefficients in couple stress fluid equations satisfy these inequalities:

(2.5) μ 0 , η 0 , η η , 3 λ + 2 μ 0 .

Introducing the dimensionless quantities to the governing Eq. (2.2), we obtain the following:

(2.6) q = q a Ω , ω = ω Ω , p = p μ Ω , = a , τ r r ˆ = τ r r μ Ω , α 2 = a 2 K , ξ 2 = a 2 μ η , r = r a , m = a m μ Ω .

Substituting Eq. (2.8) in Eq. (2.7) and then dropping the primes, we obtain

(2.7) p + q + 1 ξ 2 q + α 2 q = 0 ,

where α is the permeability parameter, and ξ is the length-dependent parameter of the first couple stress fluid coefficient. If ξ = ( η = 0 ) , Eq. (2.7) denotes the improved Stokes equation for non-polar fluids.

Furthermore, the relations of non-dimensional stress and couple stress are obtained by

(2.8) τ r ϕ = q ϕ r q ϕ r + 1 2 m r θ r + 1 r ( 2 m r θ + m θ r ) + cot θ r ( m θ θ m ϕ ϕ ) + 1 r m θ θ θ ,

where

(2.9) m r r = m ¯ + 4 1 ξ 2 + 1 ξ 2 ω r r ,

(2.10) m r θ = 4 ξ 2 ω θ r + 4 ξ 2 1 r ω r θ ω θ ,

(2.11) m θ θ = m ¯ + 4 1 r 1 ξ 2 + 1 ξ 2 ω θ θ + ω r ,

(2.12) m θ r = 4 ξ 2 1 r ω r θ ω θ + 4 ξ 2 ω θ r ,

(2.13) m ϕ ϕ = m ¯ + 4 1 r 1 ξ 2 + 1 ξ 2 ( ω r + cot θ ω θ ) ,

where ξ 2 = a 2 μ η is the length-dependent parameter on the second couple stress fluid coefficient.

3 Axisymmetric solutions to the problem

Suppose a solid sphere rotates through unbounded Brinkman couple stress fluids. Therefore, the origin of the spherical coordinate system ( r , θ , ϕ ) is positioned in the center of the sphere. The flow field functions are independent of ϕ due to symmetry. Additionally, the vectors of velocity and vorticity have the form:

(3.1) q = q ϕ ( r , θ ) e ϕ , ω = ω r e r + ω θ e θ ,

where e r , e θ , and e ϕ are the unit vectors along the coordinate lines r , θ , and ϕ . Applying Eq. (3.1) into (2.7), we obtain the fourth-order partial differential equation as follows:

(3.2) ( L 1 λ 1 2 ) ( L 1 λ 2 2 ) ( r sin θ q ϕ ) = 0 ,

where the axisymmetric Stokesian differential operator L 1 and the roots are expressed as

(3.3) L 1 = 2 r 2 + 1 ζ 2 r 2 2 ξ 2 , ζ = cos θ , λ 1 2 + λ 2 2 = ξ 2 , λ 1 2 λ 2 2 = ξ 2 α 2 λ i j = ± ξ 2 ± ξ ξ 2 4 α 2 2 , i , j = 1 , 2 .

The non-zero vorticity components ω r and ω θ are expressed as

(3.4) ω r = 1 2 ( q ) e r = 1 2 1 r sin θ θ ( sin θ q ϕ ) , ω θ = 1 2 ( q ) e θ = 1 2 1 r r ( r q ϕ ) .

On the rotating solid sphere, the following boundary conditions have two possible situations that were made assumptions by Stokes [9]: boundary conditions (A) declare that there are no couple stresses at the border of the region and (B) the second states that the angular velocity of the boundary is equivalent to the vorticity along the boundary. Physically, according to the first boundary condition, the tangential velocity at the sphere’s surface depends on the polar angle θ and is proportionate to the angular velocity. This condition emphasizes that there is no surface slip for the tangential component of the velocity and illustrates how the fluid or surface interacts with the rotating sphere. The second boundary condition indicates that the sphere is not under any extra couple stresses. This indicates that there are no moments created by the applied forces that would result in rotation around an axis. By demonstrating that the stresses are entirely normal or shear with no twisting effects, this condition streamlines the analysis:

  1. The no-slip condition:

    (3.5) q ϕ = a Ω sin θ ,

  2. No couple stresses (Type A):

    (3.6) m i j n i = 0 ,

where Ω is the characteristic angular velocity, and the range of the slip coefficient, β , is zero to infinity. Only the kind of fluid flow and the boundary’s substance affect this parameter. From this study, the perfect slipcase may be reached when the slip coefficient disappears, and the classical situation of no slip can be inferred as a special case when the slip parameter approaches infinity. In addition, n is the unit normal to the solid sphere’s surface. The aforementioned boundary conditioning was proposed by Stokes [6]. According to physical theory, a force distribution can only result from mechanical interactions at the boundaries in this particular situation.

Consequently, the differential Eq. (3.2) has the following generic solution:

(3.7) q ϕ ( r , θ ) = j = 1 2 n = 1 A n ( j ) r 1 2 K n + 1 2 ( λ 1 j r ) + B n ( j ) r 1 2 I n + 1 2 ( λ 2 j r ) P n 1 ( cos θ ) ,

where the functions K n ( ) and I n ( ) are the modified Bessel functions of the first and second kinds of order n , respectively. where A n ( j ) , B n ( j ) , j = 1 , 2 are the constants. Also, P n 1 ( ) represents the associated Legendre polynomials of degree n and order one. The components of vorticity are obtained by

(3.8) ω r ( r , θ ) = 1 2 j = 1 2 n = 1 n ( n + 1 ) r 3 2 P n ( cos θ ) × A n ( j ) K n + 1 2 ( λ 1 j r ) + B n ( j ) I n + 1 2 ( λ 2 j r ) ,

(3.9) ω θ ( r , θ ) = 1 2 j = 1 2 n = 1 r 3 2 P n 1 ( cos θ ) × A n ( j ) n K n + 1 2 ( λ 1 j r ) + λ 1 j r K n 1 2 ( λ 1 j r ) + B n ( j ) n I n + 1 2 ( λ 2 j r ) λ 2 j r I n 1 2 ( λ 2 j r ) .

Furthermore, the couple stress functions are obtained from (2.9)–(2.13) as follows:

(3.10) m r θ ( r , θ ) = 2 j = 1 2 n = 1 r 5 2 P n 1 ( cos θ ) × A n ( j ) ( ξ 2 + ξ 2 ) λ 1 j r K n 1 2 ( λ 1 j r ) + ( ξ 2 λ 1 j 2 r 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) K n + 1 2 ( λ 1 j r ) + B n ( j ) ( ( ξ 2 + ξ 2 ) × λ 2 j r I n 1 2 ( λ 2 j r ) ( ξ 2 λ 2 j 2 r 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) I n + 1 2 ( λ 2 j r ) ,

(3.11) m θ r ( r , θ ) = 2 j = 1 2 n = 1 r 5 2 P n 1 ( cos θ ) × A n ( j ) ( ξ 2 + ξ 2 ) λ 1 j r K n 1 2 ( λ 1 j r ) + ( ξ 2 λ 1 j 2 r 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) K n + 1 2 ( λ 1 j r ) + B n ( j ) ( ( ξ 2 + ξ 2 ) × λ 2 j r I n 1 2 ( λ 2 j r ) ( ξ 2 λ 2 j 2 r 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) I n + 1 2 ( λ 2 j r ) ,

(3.12) m θ θ ( r , θ ) = m + 2 ( ξ 2 + ξ 2 ) j = 1 2 n = 1 ( n + 1 ) r 5 2 P n ( cos θ ) × [ A n ( j ) csc θ λ 1 j r K n 1 2 ( λ 1 j r ) + n ( csc θ 1 ) K n + 1 2 ( λ 1 j r ) + B n ( j ) csc θ λ 2 j r I n 1 2 ( λ 2 j r ) n ( csc θ 1 ) I n + 1 2 ( λ 2 j r ) + r 5 2 cot θ P n 1 ( cos θ ) × A n ( j ) λ 1 j r K n 1 2 ( λ 1 j r ) + n K n + 1 2 ( λ 1 j r ) B n ( j ) λ 2 j r I n 1 2 ( λ 2 j r ) n I n + 1 2 ( λ 2 j r ) ,

(3.13) m ϕ ϕ ( r , θ ) = m + 2 ( ξ 2 + ξ 2 ) r 5 2 j = 1 2 n = 1 n ( n + 1 ) P n ( cos θ ) × A n ( j ) K n + 1 2 ( λ 1 j r ) + B n ( j ) I n + 1 2 ( λ 2 j r ) + cot θ P n 1 ( cos θ ) A n ( j ) λ 1 j r K n 1 2 ( λ 1 j r ) + n K n + 1 2 ( λ 1 j r ) B n ( j ) λ 2 j r I n 1 2 ( λ 2 j r ) n I n + 1 2 ( λ 2 j r ) .

Hence, from Eq. (3.6), we obtain m r r = 0 on the surface of the object.

(3.14) m r r ( r , θ ) = 2 ( ξ 2 + ξ 2 ) j = 1 2 n = 1 r 3 2 n ( n + 1 ) P n 1 ( cos θ ) × A n ( j ) λ 1 j r K n 1 2 ( λ 1 j r ) + ( n + 2 ) K n + 1 2 ( λ 1 j r ) B n ( j ) λ 2 j r I n 1 2 ( λ 2 j r ) ( n + 2 ) I n + 1 2 ( λ 2 j r ) .

Inserting Eqs (3.7) and (3.10)–(3.14) into (2.8), we obtain

(3.15) τ r ϕ ( r , θ ) = r 7 2 j = 1 2 n = 1 2 σ cos θ csc 2 θ P n ( cos θ ) × A n ( j ) λ 1 j r K n 1 2 ( λ 1 j r ) + n K n + 1 2 ( λ 1 j r ) + B n ( j ) λ 2 j r I n 1 2 ( λ 2 j r ) n I n + 1 2 ( λ 2 j r ) + csc 2 θ P n 1 ( cos θ ) [ A n ( j ) ( ( cos 2 θ ( σ 1 + σ csc θ ) + σ 3 + σ csc θ ) λ 1 j r K n 1 2 ( λ 1 j r ) + ( cos 2 θ ( σ 7 + n σ csc θ ) + σ 5 + n σ csc θ ) K n + 1 2 ( λ 1 j r ) + B n ( j ) ( cos θ 2 ( σ 2 + σ csc θ ) σ 4 σ csc θ ) λ 2 j r I n 1 2 ( λ 2 j r ) + ( cos 2 θ ( σ 8 + n σ csc θ ) + σ 6 + n σ csc θ ) I n + 1 2 ( λ 2 j r ) ,

where

σ = ( n + 1 ) ( ξ 2 + ξ 2 ) , σ 1 = λ 1 j 2 r 2 ξ 2 r 2 2 ξ 2 2 ξ 2 , σ 2 = λ 2 j 2 r 2 ξ 2 r 2 2 ξ 2 2 ξ 2 , σ 3 = λ 1 j 2 r 2 ξ 2 r 2 ξ 2 ξ 2 , σ 4 = λ 2 j 2 r 2 ξ 2 r 2 ξ 2 ξ 2 , σ 5 = λ 1 j 2 r 2 ξ 2 n ( n + 2 ) r 2 n ( 2 n + 3 ) ( ξ 2 + ξ 2 ) , σ 6 = λ 2 j 2 r 2 ξ 2 n ( n + 2 ) r 2 n ( 2 n + 3 ) ( ξ 2 + ξ 2 ) , σ 7 = λ 1 j 2 r 2 ξ 2 n ( n + 2 ) r 2 2 n σ , σ 8 = λ 2 j 2 r 2 ξ 2 n ( n + 2 ) r 2 2 n σ .

4 Torque due to stress and couple stress tensors

Torque in fluids arises from the forces due to stress and, in couple stress fluids, additional microstructural effects. In classical fluids, torque results from pressure and viscous forces acting on surfaces, calculated using the stress tensor. In couple stress fluids, which account for internal body and surface couples, torque includes contributions from both the stress tensor and the couple stress tensor, reflecting the fluid’s microstructure. These effects are significant in systems such as suspensions, polymer solutions, and microfluidic devices, influencing rotational dynamics, lubrication, and energy dissipation. This enhances the fluid’s ability to model behaviors in complex systems.

The torque T z e z of the fluid flow affecting the axisymmetric particle is controlled on the solid of revolution and calculated according to the following [27]:

(4.1) T z = 2 π μ a 3 0 π τ r ϕ r = a sin 2 θ d θ .

Substituting Eq. (3.15) into formula (4.1), we obtain

(4.2) T z = 2 π μ a 1 2 j = 1 2 π A 1 ( j ) λ 1 j a K 1 2 ( λ 1 j a ) + K 3 2 ( λ 1 j a ) + B 1 ( j ) λ 2 j a I 1 2 ( λ 2 j a ) I 3 2 ( λ 2 j a ) + I A 1 ( j ) 2 3 σ 1 + 1 π 2 σ + 2 σ 3 λ 1 j a K 1 2 ( λ 1 j a ) + 2 3 σ 7 + 1 π 2 σ + 2 σ 5 K 3 2 ( λ 1 j a ) + B 1 ( j ) 2 3 σ 2 π 2 σ 2 σ 4 λ 2 j a I 1 2 ( λ 2 j a ) + 2 3 σ 8 + π 2 σ + 2 σ 6 I 3 2 ( λ 2 j a ) .

The torque T m due to the couple stress tensor exerted on a spherical particle, a 1 is

(4.3) T m = 2 π μ a 2 0 π ( m r r cos θ m r θ sin θ ) r = a sin θ d θ .

Substituting (3.10) and (3.14) into (4.3), we obtain

(4.4) T m = 2 I π 2 μ a 1 2 j = 1 2 × A 1 ( j ) ( ξ 2 + ξ 2 ) λ 1 j a K 1 2 ( λ 1 j a ) + ( ξ 2 λ 1 j 2 a 2 + ( ξ 2 + 3 ξ 2 ) ) K 3 2 ( λ 1 j a ) + B 1 ( j ) ( ξ 2 + ξ 2 ) λ 2 j a I 1 2 ( λ 2 j a ) ( ξ 2 λ 2 j 2 a 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) I 3 2 ( λ 2 j a ) .

5 Interaction between two eccentric spheres

We examine the steady couple stress incompressible fluid flow produced by two rotating eccentric spheres of radii a 1 and a 2 with two angular velocities Ω 1 and Ω 2 , respectively, around their common diameter ( z -axis), as illustrated in Figure 1. The spherical coordinates, in this case, are ( r 2 , θ 2 , ϕ 2 ) with the origin at the outside sphere center and ( r 1 , θ 1 , ϕ 1 ) based on the inside solid sphere center, which is situated at a distance of h from the cavity center. Suppose that couple stress fluids rotate consistently in an axisymmetric direction around a sphere of radius a 1 . The flow field functions are invariant with respect to ϕ due to symmetry. Furthermore, the velocity and vorticity vectors are expressed as

(5.1) q = q ϕ ( r , θ ) e ϕ and ω = ω r e r + ω θ e θ ,

where e r , e θ , and e ϕ represent the unit vectors along the coordinate lines r , θ , and ϕ . Furthermore, the coordinate systems ( r 1 , θ 1 , ϕ 1 ) and ( r 2 , θ 2 , ϕ 2 ) are connected by the following relations ( r 1 , θ 1 ) and ( r 2 , θ 2 ) is given by r 1 2 = r 2 2 + h 2 2 r 2 h cos θ 2 or by r 2 2 = r 1 2 + h 2 + 2 r 1 h cos θ 1 . The equations’ linearity enables the application of the superposition principle. However, we choose the velocity q ϕ ( j ) and the tangential couple stress functions m r θ ( j ) in the form:

(5.2) q ϕ ( j ) ( r j , θ j ) = q ϕ ( 1 ) ( r 1 , θ 1 ) + q ϕ ( 2 ) ( r 2 , θ 2 ) , m r j θ j ( j ) ( r j , θ j ) = m r 1 θ 1 ( 1 ) ( r 1 , θ 1 ) + m r 2 θ 2 ( 2 ) ( r 2 , θ 2 ) .

From the differential Eq. (3.2), four boundary conditions are required to describe the problem. The boundary conditions on the surfaces of the spheres are established as follows:

(5.3) r 1 = a 1 : q ϕ = a 1 Ω 1 sin θ 1 , m r 1 θ 1 = 0 ,

(5.4) r 2 = a 2 : q ϕ = a 2 Ω 2 sin θ 2 , m r 2 θ 2 = 0 ,

where Ω j , j = 1 , 2 are the angular velocities of the two spheres. Accordingly, the mechanical interaction at the wall is equivalent to the no-slip condition at the boundary, which states that all solid borders have zero fluid velocity. In addition, boundary walls can prevent fluid elements from rotating, so there is no zero couple stress when it comes to boundary walls. The field functions, the velocity, and the couple stress are given by

(5.5) q ϕ ( r , θ ) = j = 1 2 n = 1 A n ( j ) r 1 1 2 K n + 1 2 ( λ 1 j r 1 ) P n 1 ( cos θ 1 ) + j = 1 2 n = 1 B n ( j ) r 2 1 2 I n + 1 2 ( λ 2 j r 2 ) P n 1 ( cos θ 2 ) ,

(5.6) m θ r ( r , θ ) = 2 j = 1 2 n = 1 r 1 5 2 P n 1 ( cos θ 1 ) A n ( j ) × ( ξ 2 + ξ 2 ) λ 1 j r 1 K n 1 2 ( λ 1 j r 1 ) + ( ξ 2 λ 1 j 2 r 1 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) K n + 1 2 ( λ 1 j r 1 ) + 2 j = 1 2 n = 1 r 2 5 2 P n 1 ( cos θ 2 ) B n ( j ) ( ξ 2 + ξ 2 ) λ 2 j r 2 I n 1 2 ( λ 2 j r 2 ) ( ξ 2 λ 2 j 2 r 2 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) I n + 1 2 ( λ 2 j r 2 ) .

Figure 1 
               Sketch of Brinkman couple stress fluid between two rotating eccentric spheres.
Figure 1

Sketch of Brinkman couple stress fluid between two rotating eccentric spheres.

Substituting the field functions from (5.5) and (5.6) into the boundary conditions from Eqs (5.3) and (5.4), we have the linear system of equations in the four unknowns A n ( j ) , B n ( j ) , j = 1 , 2 as

(5.7) a 1 Ω 1 sin θ 1 = j = 1 2 n = 1 A n ( j ) a 1 1 2 K n + 1 2 ( λ 1 j a 1 ) P n 1 ( cos θ 1 ) + j = 1 2 n = 1 B n ( j ) r 2 1 2 I n + 1 2 ( λ 2 j r 2 ) r 1 = a 1 P n 1 ( cos θ 2 ) ,

(5.8) a 2 Ω 2 sin θ 2 = j = 1 2 n = 1 A n ( j ) r 1 1 2 K n + 1 2 ( λ 1 j r 1 ) r 2 = a 2 P n 1 ( cos θ 1 ) + j = 1 2 n = 1 B n ( j ) a 2 1 2 I n + 1 2 ( λ 2 j a 2 ) P n 1 ( cos θ 2 ) ,

(5.9) 0 = 2 j = 1 2 n = 1 a 1 5 2 P n 1 ( cos θ 1 ) A n ( j ) × ( ξ 2 + ξ 2 ) λ 1 j a 1 K n 1 2 ( λ 1 j a 1 ) + ( ξ 2 λ 1 j 2 a 1 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) × K n + 1 2 ( λ 1 j a 1 ) + 2 j = 1 2 n = 1 r 2 5 2 P n 1 ( cos θ 2 ) × B n ( j ) ( ξ 2 + ξ 2 ) λ 2 j r 2 I n 1 2 ( λ 2 j r 2 ) ( ξ 2 λ 2 j 2 r 2 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) I n + 1 2 ( λ 2 j r 2 ) r 1 = a 1 ,

(5.10) 0 = 2 j = 1 2 n = 1 a 1 5 2 P n 1 ( cos θ 1 ) A n ( j ) × ( ξ 2 + ξ 2 ) λ 1 j a 1 K n 1 2 ( λ 1 j a 1 ) + ( ξ 2 λ 1 j 2 a 1 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) K n + 1 2 ( λ 1 j a 1 ) r 2 = a 2 + 2 j = 1 2 n = 1 r 2 5 2 P n 1 ( cos θ 2 ) B n ( j ) × ( ξ 2 + ξ 2 ) λ 2 j r 2 I n 1 2 ( λ 2 j r 2 ) ( ξ 2 λ 2 j 2 r 2 2 + ( ξ 2 + ξ 2 ) n ( n + 2 ) ) I n + 1 2 ( λ 2 j r 2 ) .

Then, from Eqs (4.2) and (4.4), the torque due to stress and couple stress forces T z and T m acting on the inside particle a 1 = 1 is given by

(5.11) T z = 2 π μ a 1 1 2 j = 1 2 π × A 1 ( j ) λ 1 j a 1 K 1 2 ( λ 1 j a 1 ) + K 3 2 ( λ 1 j a 1 ) + B 1 ( j ) λ 2 j a 1 I 1 2 ( λ 2 j a 1 ) I 3 2 ( λ 2 j a 1 ) + I A 1 ( j ) 2 3 σ 1 + 1 π 2 σ + 2 σ 3 λ 1 j a 1 K 1 2 ( λ 1 j a 1 ) + 2 3 σ 7 + 1 π 2 σ + 2 σ 5 K 3 2 ( λ 1 j a 1 ) + B 1 ( j ) 2 3 σ 2 π 2 σ 2 σ 4 λ 2 j a 1 I 1 2 ( λ 2 j a 1 ) + 2 3 σ 8 + π 2 σ + 2 σ 6 I 3 2 ( λ 2 j a 1 ) ,

(5.12) T m = 2 I π 2 μ a 1 1 2 j = 1 2 A 1 ( j ) ( ξ 2 + ξ 2 ) λ 1 j a 1 K 1 2 ( λ 1 j a 1 ) + ( ξ 2 λ 1 j 2 a 1 2 + ( ξ 2 + 3 ξ 2 ) ) K 3 2 ( λ 1 j a 1 ) + B 1 ( j ) ( ξ 2 + ξ 2 ) λ 2 j a 1 I 1 2 ( λ 2 j a 1 ) ( ξ 2 λ 2 j 2 a 1 2 + ( ξ 2 + 3 ξ 2 ) ) I 3 2 ( λ 2 j a 1 ) .

6 Results and discussions

In order to satisfy the boundary condition over the entire semicircular arc of a sphere, it is necessary to resolve an infinite number of unknown constants. As a result of boundary collocation, the infinite series can be truncated into manageable finite representations by applying the boundary condition at a finite number of discrete points along the arc. It is possible to truncate the series after N terms by imposing a boundary condition at N discrete points, which leads to a system of 2 N simultaneous linear equations. In order to clarify the unknown constants, the equations can be solved using a matrix reduction technique. By increasing N , the accuracy of this truncation can be improved, with the magnitude of the truncation error diminishing as N approaches infinity. This combination of methods provides a robust framework for analyzing the axisymmetric creeping flow generated by the particles in a couple stress fluid. For the application of boundary conditions in fluid dynamics, this method emphasizes the strategic selection of points along the semicircular arc of a sphere.

Overall, our initial point of selection is θ = π 2 , since this defines the area projected from the particle normal to the direction of motion. Furthermore, the points θ = 0 and θ = π are also significant, although their use may lead to the formation of a singular coefficient matrix in the equations for the unknown constants A n and B n . There are four fundamental collocation points that are chosen to address this issue while maintaining symmetry: θ = ε , θ = π 2 ε , and θ = π a . The optimal value of ε is 0.01 , which allows the numerical results for drag force to converge to at least four significant figures. Using a sufficient number of well-distributed collocation points, the solution for the torque can be approximated accurately, regardless of the particle’s actual form or boundary conditions in previous studies [2830].

The consequences of these studies are then explored in detail, along with how the presence of couple stresses affect the fluid’s behavior under different flow circumstances. When a fluid is able to withstand internal couple stresses, for example, its flow resistance changes, or it produces unique stress distributions compared with other fluids. In addition to discussing flow across porous media, it would also be interesting to discuss the impact of couple stress fluids on the microstructure of the porous matrix, which could alter permeability or other important parameters. The relevant parameters are the length of the first and second couple stress parameters, ξ and ξ , which take the range from 1.0 to , the permeability parameter, α , which value from 1 to 10 , the separation distance δ = ( a 2 a 1 ) h = 1.001 10 , the size ratio a 1 a 2 = 0.1 0.99 , and the angular speeds Ω 2 Ω 1 . This work discussed the normalized torque with respect to stress tensor T z T and couple stress tensor T m T that appears graphically from Figures 211 and Table 1.

Figure 2 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 z
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{z}/{T}_{\infty }
                     
                   
                  versus the separation distance for various values of angular velocity ratio 
                     
                        
                        
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                        
                        {\Omega }_{2}/{\Omega }_{1}
                     
                   at 
                     
                        
                        
                           ξ
                           =
                           4.0
                           ,
                           
                              
                                 a
                              
                              
                                 1
                              
                           
                           ⁄
                           
                              
                                 a
                              
                              
                                 2
                              
                           
                           =
                           0.2
                           ,
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           =
                           2.0
                        
                        \xi =4.0,{a}_{1}/{a}_{2}=0.2,{\xi }^{^{\prime} }=2.0
                     
                  , and 
                     
                        
                        
                           α
                           =
                           1.5
                        
                        \alpha =1.5
                     
                  .
Figure 2

Normalized torque T z T versus the separation distance for various values of angular velocity ratio Ω 2 Ω 1 at ξ = 4.0 , a 1 a 2 = 0.2 , ξ = 2.0 , and α = 1.5 .

Table 1

Dimensionless torque T z T on the rigid sphere a 1 with ξ = ξ = 30 and the two spheres rotating in opposite directions

a 2 a 1 h a 1 a 2 α = 0.0 [9] α = 1.0 α = 2.0 α = 3.0 α = 4.0 α = 10.0
1.05 0.1 2.122441 0.009290 0.013862 0.018227 0.021941 0.024624
0.2 2.269303 0.014601 0.014031 0.018232 0.021941 0.024624
0.3 2.430441 0.026206 0.016550 0.018576 0.021983 0.024624
0.4 2.631674 0.038509 0.023930 0.021166 0.022742 0.024624
0.5 2.910755 0.047317 0.034613 0.028292 0.026541 0.024649
0.6 3.335021 0.066856 0.051530 0.041589 0.035944 0.025059
0.7 4.056916 0.091365 0.074623 0.062736 0.054287 0.028426
0.8 5.530868 0.141426 0.119717 0.103982 0.092054 0.044943
0.9 10.023727 0.296392 0.259082 0.231912 0.210972 0.119431
0.99 91.579170 3.090827 2.735826 2.466769 2.249198 1.119105
0.999 907.585754 30.995344 27.398773 24.645763 22.394852 10.498260
4.0 0.1 2.010692 0.009295 0.013862 0.018227 0.021941 0.024624
0.2 2.026726 0.015348 0.014058 0.018232 0.021941 0.024624
0.3 2.069914 0.028412 0.016953 0.018630 0.021989 0.024624
0.4 2.156706 0.040973 0.024844 0.021446 0.022817 0.024624
0.5 2.312791 0.047696 0.034379 0.027981 0.026323 0.024646
0.6 2.587229 0.066756 0.051112 0.041254 0.035741 0.025053
0.7 3.093276 0.089896 0.073039 0.061348 0.053186 0.028289
0.8 4.170103 0.137255 0.115798 0.100496 0.089002 0.043807
0.9 7.512265 0.284453 0.248203 0.222316 0.202500 0.115038
0.99 68.504509 2.943968 2.599263 2.342640 2.135705 1.033162
0.999 678.908569 29.496065 25.994425 23.357748 21.205435 9.537851
10.0 0.1 2.010401 0.009307 0.013862 0.018227 0.021941 0.024624
0.2 2.024912 0.015545 0.014066 0.018233 0.021941 0.024624
0.3 2.065131 0.028598 0.016992 0.018636 0.021990 0.024624
0.4 2.147774 0.041094 0.024896 0.021464 0.022822 0.024624
0.5 2.298784 0.047732 0.034379 0.027976 0.026321 0.024646
0.6 2.567110 0.066757 0.051095 0.041242 0.035736 0.025054
0.7 3.065037 0.089818 0.072955 0.061275 0.053130 0.028285
0.8 4.128208 0.137021 0.115577 0.100300 0.088830 0.043746
0.9 7.433049 0.283770 0.247572 0.221747 0.201987 0.114743
0.99 67.771484 2.935644 2.591443 2.335432 2.129021 1.027992
0.999 671.643616 29.411222 25.914204 23.283249 21.135777 9.480927
0.1 2.010351 0.009310 0.013862 0.018227 0.021941 0.024624
0.2 2.024594 0.015580 0.014068 0.018233 0.021941 0.024624
0.3 2.064278 0.028634 0.016999 0.018637 0.021990 0.024624
0.4 2.146158 0.041117 0.024906 0.021467 0.022823 0.024624
0.5 2.296226 0.047739 0.034379 0.027975 0.026321 0.024646
0.6 2.563412 0.066757 0.051092 0.041240 0.035735 0.025054
0.7 3.059825 0.089803 0.072939 0.061262 0.053120 0.028284
0.8 4.120455 0.136976 0.115536 0.100263 0.088798 0.043734
0.9 7.418374 0.283641 0.247452 0.221639 0.201890 0.114686
0.99 67.635643 2.934221 2.590061 2.334135 2.127805 1.027024
0.999 670.301086 29.415516 25.912762 23.278940 21.129986 9.472232

Figure 2 illustrates the normalized torque T z T as a function of the separation distance ( a 2 a 1 ) h for different values of the angular velocity ratio Ω 2 Ω 1 . The data are presented for a fixed set of parameters: ξ = 4.0 , a 1 = 0.2 , a 2 = 2.0 , and α = 1.5 . The curves show how normalized torque varies with separation distance for various angular velocity ratios, ranging from Ω 2 Ω 1 = 4.0 to other values. The torque remains relatively constant for most values of ( a 2 a 1 ) h , suggesting a stable behavior in the system across different conditions. The variations in the curves indicate how changes in the angular velocity ratio influence the torque for the given configuration. This figure provides insights into the relationship between torque and separation distance, highlighting the impact of angular velocity ratios in the system dynamics.

Figure 3 presents the normalized torque T z T as a function of the separation distance ( a 2 a 1 ) h for various values of the ratio a 1 a 2 . The analysis is based on a fixed set of parameters: ξ = 10.0 , Ω 2 Ω 1 = 10.0 , and α = 1.0 . The curves indicate that as the separation distance approaches zero, the normalized torque experiences a significant drop, particularly for lower values of a 1 a 2 (e.g., 0.1 ). This behavior highlights the strong dependency of torque on the geometric configuration and the angular velocity ratio. The marked differences in the torque values suggest that the system dynamics are highly sensitive to the separation distance, with larger values of a 1 a 2 leading to more stable torque levels at increased distances.

Figure 3 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 z
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{z}/{T}_{\infty }
                     
                   
                  versus the separation distance for various values of angular velocity ratio 
                     
                        
                        
                           
                              
                                 a
                              
                              
                                 1
                              
                           
                           ⁄
                           
                              
                                 a
                              
                              
                                 2
                              
                           
                        
                        {a}_{1}/{a}_{2}
                     
                   at 
                     
                        
                        
                           ξ
                           =
                           10.0
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           10.0
                           ,
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           →
                           ∞
                        
                        \xi =10.0,{\Omega }_{2}/{\Omega }_{1}=10.0,{\xi }^{^{\prime} }\to \infty 
                     
                  , and 
                     
                        
                        
                           α
                           =
                           1.0
                        
                        \alpha =1.0
                     
                  .
Figure 3

Normalized torque T z T versus the separation distance for various values of angular velocity ratio a 1 a 2 at ξ = 10.0 , Ω 2 Ω 1 = 10.0 , ξ , and α = 1.0 .

Figure 4 depicts the normalized torque T z T as a function of the separation distance ( a 2 a 1 ) h for various values of the parameter ξ while holding a 1 a 2 = 0.4 , Ω 2 Ω 1 = 4.0 , and α = 1.0 . The analysis reveals a pronounced sensitivity of the normalized torque to variations in the separation distance, particularly as it approaches zero. At lower values of ξ (e.g., 2.5), the torque exhibits a dramatic increase, highlighting the enhanced coupling effects within the system under tighter spatial constraints. Conversely, higher values of ξ reflect a more gradual torque response, indicating a transition toward a more stable dynamical regime. This behavior can be attributed to the interplay between angular momentum transfer and the geometric configuration, which collectively dictate the system’s mechanical stability. The steep gradients observed in the curves illustrate the non-linear characteristics of torque as influenced by both the separation distance and the angular velocity ratio.

Figure 4 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 z
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{z}/{T}_{\infty }
                     
                   
                  versus the separation distance for various values of angular velocity ratio 
                     
                        
                        
                           ξ
                        
                        \xi 
                     
                   at 
                     
                        
                        
                           
                              
                                 a
                              
                              
                                 1
                              
                           
                           ⁄
                           
                              
                                 a
                              
                              
                                 2
                              
                           
                           =
                           0.4
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           4.0
                           ,
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           =
                           1.0
                        
                        {a}_{1}/{a}_{2}=0.4,{\Omega }_{2}/{\Omega }_{1}=4.0,{\xi }^{^{\prime} }=1.0
                     
                  , and 
                     
                        
                        
                           α
                           =
                           1.0
                        
                        \alpha =1.0
                     
                  .
Figure 4

Normalized torque T z T versus the separation distance for various values of angular velocity ratio ξ at a 1 a 2 = 0.4 , Ω 2 Ω 1 = 4.0 , ξ = 1.0 , and α = 1.0 .

Figure 5 illustrates the normalized torque T z T as a function of the separation distance ( a 2 a 1 ) h for various values of the angular velocity ratio ξ , with parameters fixed at a 1 a 2 = 0.4 , Ω 2 Ω 1 = 4.0 , ξ = 10.0 , and α = 1.0 . The data reveal a striking characteristic: as a s h increases, the normalized torque stabilizes at values near unity, irrespective of the specific value of ξ . This observation underscores a significant decoupling of the torque from the separation distance for larger ξ , indicating that the system approaches a state of dynamic equilibrium. At lower values of ξ (e.g., 5.0), there is a more pronounced response to changes in separation, suggesting that the system remains sensitive to spatial configurations and angular momentum interactions. However, as ξ increases, the curves converge, indicating a reduction in the influence of angular velocity on the torque dynamics. This behavior highlights the transition from a regime dominated by angular momentum transfer dynamics to one where the geometric factors play a more pivotal role. Overall, Figure 5 provides valuable insights into the torque characteristics of the system, illustrating how varying angular velocities can modulate the torque response over different spatial configurations, ultimately leading to a regime of stability at higher separations.

Figure 5 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 z
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{z}/{T}_{\infty }
                     
                   
                  versus the separation distance for various values of angular velocity ratio 
                     
                        
                        
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                        
                        {\xi }^{^{\prime} }
                     
                   at 
                     
                        
                        
                           
                              
                                 a
                              
                              
                                 1
                              
                           
                           ⁄
                           
                              
                                 a
                              
                              
                                 2
                              
                           
                           =
                           0.4
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           4.0
                           ,
                           ξ
                           =
                           10.0
                        
                        {a}_{1}/{a}_{2}=0.4,{\Omega }_{2}/{\Omega }_{1}=4.0,\xi =10.0
                     
                  , and 
                     
                        
                        
                           α
                           =
                           1.0
                        
                        \alpha =1.0
                     
                  .
Figure 5

Normalized torque T z T versus the separation distance for various values of angular velocity ratio ξ at a 1 a 2 = 0.4 , Ω 2 Ω 1 = 4.0 , ξ = 10.0 , and α = 1.0 .

Figure 6 presents the normalized torque T z T as a function of the separation distance ( a 2 a 1 ) h for various values of the angular velocity ratio α , while maintaining a 1 a 2 = 0.4 , Ω 2 Ω 1 = 4.0 , and ξ = . The trends observed in this figure reveal critical insights into the torque behavior as influenced by the angular velocity ratio. At lower values of α (e.g., 1.0 ), the normalized torque exhibits a relatively steep rise with decreasing separation distance, indicating a robust coupling between the components of the system. This behavior suggests that as the components come closer, the interactions become increasingly significant, leading to enhanced torque generation. The underlying physics can be attributed to the increased angular momentum transfer and the proximity-induced coupling effects, which dominate the torque dynamics in this regime. As α increases, particularly beyond 2.0 , a noticeable shift in the torque response is observed. The curves begin to flatten out, suggesting that the system approaches a more stable regime where the torque becomes less sensitive to changes in separation distance.

Figure 6 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 z
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{z}/{T}_{\infty }
                     
                   
                  versus the separation distance for various values of angular velocity ratio 
                     
                        
                        
                           α
                        
                        \alpha 
                     
                   at 
                     
                        
                        
                           
                              
                                 a
                              
                              
                                 1
                              
                           
                           ⁄
                           
                              
                                 a
                              
                              
                                 2
                              
                           
                           =
                           0.9
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           4.0
                        
                        {a}_{1}/{a}_{2}=0.9,{\Omega }_{2}/{\Omega }_{1}=4.0
                     
                  , and 
                     
                        
                        
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           =
                           ξ
                           →
                           ∞
                        
                        {\xi }^{^{\prime} }=\xi \to \infty 
                     
                  .
Figure 6

Normalized torque T z T versus the separation distance for various values of angular velocity ratio α at a 1 a 2 = 0.9 , Ω 2 Ω 1 = 4.0 , and ξ = ξ .

This indicates that higher angular velocities contribute to a more uniform distribution of forces within the system, thereby reducing the torque fluctuations as the separation varies. For larger values of α (e.g., 10.0 ), the normalized torque approaches a near-constant value across a wide range of separation distances. This behavior highlights a significant decoupling of torque from spatial configurations, suggesting that the dynamics are increasingly governed by the balance of forces rather than the geometric arrangement of the components. Overall, Figure 6 illustrates how the angular velocity ratio significantly influences the torque response in the system. The transition from a sensitivity to separation distance at lower α values to a stable, nearly constant torque at higher α values emphasizes the complex interplay between angular momentum dynamics and geometric configurations, ultimately shedding light on the stability characteristics of the system under varying operational conditions.

Figure 7 illustrates the normalized torque T 1 T as a function of the size ratio ξ 1 ξ 2 for various values of the angular velocity ratio Ω 2 Ω 1 , with other parameters fixed: α = 1.5 , ξ = 1.0 , and ξ = 4.0 . The graph clearly demonstrates how the normalized torque behaves under different conditions of size and velocity ratios. As the size ratio increases, the torque initially shows a positive trend for certain angular velocity ratios, particularly for the higher values of Ω 2 Ω 1 (represented by the magenta and cyan curves). In contrast, for lower velocity ratios (shown by the black and blue curves), the normalized torque decreases and eventually becomes negative. The trend indicates a significant influence of the angular velocity ratio on the torque’s direction and magnitude, highlighting how the relative sizes of the objects and their rotational velocities impact the overall system dynamics.

Figure 7 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 z
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{z}/{T}_{\infty }
                     
                   
                  versus the size ratio for various values of angular velocity ratio 
                     
                        
                        
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                        
                        {\Omega }_{2}/{\Omega }_{1}
                     
                   at 
                     
                        
                        
                           
                              
                                 
                                    
                                       a
                                    
                                    
                                       2
                                    
                                 
                                 ‒
                                 
                                    
                                       a
                                    
                                    
                                       1
                                    
                                 
                              
                              
                                 h
                              
                           
                           =
                           1.5
                           ,
                           α
                           =
                           1.0
                           ,
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           =
                           1.0
                        
                        \frac{{a}_{2}‒{a}_{1}}{h}=1.5,\alpha =1.0,{\xi }^{^{\prime} }=1.0
                     
                  , and 
                     
                        
                        
                           ξ
                           =
                           4.0
                        
                        \xi =4.0
                     
                  .
Figure 7

Normalized torque T z T versus the size ratio for various values of angular velocity ratio Ω 2 Ω 1 at a 2 a 1 h = 1.5 , α = 1.0 , ξ = 1.0 , and ξ = 4.0 .

Figure 8 presents the normalized torque T m T as a function of the size ratio ξ 1 ξ 2 for various values of permeability parameter α , while keeping the other parameters constant: ξ = 10.0 , ξ = 20.0 , and Ω 2 Ω 1 = 2.0 . The graph shows a clear increase in normalized torque with the increasing size ratio, and this effect becomes more pronounced as the permeability parameter α increases. For smaller values of α (represented by the black and blue curves), the torque increases gradually with the size ratio, but the slope remains relatively mild. As α increases, particularly for α = 4.0 and α = 5.0 (indicated by the pink and green curves), the torque increases at a much sharper rate. This suggests that higher permeability values significantly amplify the effect of the size ratio on the normalized torque, leading to a much steeper rise in torque. This behavior is important for systems where permeability plays a key role, such as in porous materials or fluid–structure interaction problems, where changes in permeability can have a substantial impact on the overall mechanical or fluid dynamic performance.

Figure 8 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 m
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{m}/{T}_{\infty }
                     
                   
                  versus the size ratio for various values of permeability 
                     
                        
                        
                           α
                        
                        \alpha 
                     
                   at 
                     
                        
                        
                           
                              
                                 
                                    
                                       a
                                    
                                    
                                       2
                                    
                                 
                                 ‒
                                 
                                    
                                       a
                                    
                                    
                                       1
                                    
                                 
                              
                              
                                 h
                              
                           
                           =
                           10.0
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           ‒
                           2.0
                        
                        \frac{{a}_{2}‒{a}_{1}}{h}=10.0,{\Omega }_{2}/{\Omega }_{1}=‒2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           =
                           10.0
                        
                        {\xi }^{^{\prime} }=10.0
                     
                  , and 
                     
                        
                        
                           ξ
                           =
                           20.0
                        
                        \xi =20.0
                     
                  .
Figure 8

Normalized torque T m T versus the size ratio for various values of permeability α at a 2 a 1 h = 10.0 , Ω 2 Ω 1 = 2.0 , ξ = 10.0 , and ξ = 20.0 .

Figure 9 illustrates the normalized torque T m T as a function of the size ratio ξ 1 ξ 2 for various values of angular velocity ratio Ω 2 Ω 1 , while keeping other parameters fixed: α = 1.0 , ξ = 1.0 , and ξ = 4.0 . The results show that the torque behavior is strongly influenced by the angular velocity ratio. For negative values of Ω 2 Ω 1 , represented by the curves for Ω 2 Ω 1 = 4.0 (black) and Ω 2 Ω 1 = 2.0 (blue), the normalized torque decreases initially as the size ratio increases and eventually rises steeply for larger size ratios. The other curves, representing lower angular velocity ratios (e.g., Ω 2 Ω 1 = 0.0 , Ω 2 Ω 1 = 2.0 , and Ω 2 Ω 1 = 4.0 ), exhibit a more gradual increase in torque with the size ratio. The trends indicate that for negative angular velocity ratios, the torque tends to decrease initially before increasing, whereas for positive angular velocity ratios, the torque steadily increases with the size ratio. This behavior underscores the significant role of angular velocity differences in controlling the torque generated in such systems, particularly in configurations where the angular velocities of the two objects differ.

Figure 9 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 m
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{m}/{T}_{\infty }
                     
                   
                  versus the size ratio for various values of angular velocity ratio 
                     
                        
                        
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                        
                        {\Omega }_{2}/{\Omega }_{1}
                     
                   at 
                     
                        
                        
                           
                              
                                 
                                    
                                       a
                                    
                                    
                                       2
                                    
                                 
                                 ‒
                                 
                                    
                                       a
                                    
                                    
                                       1
                                    
                                 
                              
                              
                                 h
                              
                           
                           =
                           1.5
                        
                        \frac{{a}_{2}‒{a}_{1}}{h}=1.5
                     
                  , 
                     
                        
                        
                           α
                           =
                           1.0
                        
                        \alpha =1.0
                     
                  , 
                     
                        
                        
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           =
                           1.0
                        
                        {\xi }^{^{\prime} }=1.0
                     
                  , and 
                     
                        
                        
                           ξ
                           =
                           4.0
                        
                        \xi =4.0
                     
                  .
Figure 9

Normalized torque T m T versus the size ratio for various values of angular velocity ratio Ω 2 Ω 1 at a 2 a 1 h = 1.5 , α = 1.0 , ξ = 1.0 , and ξ = 4.0 .

Figure 10 illustrates the normalized torque T m T 0 as a function of the size ratio for different values of the length of the first couple stress parameter ξ . The graph plot values for ξ = 2.5 , 4.0 , 6.0 , and 8.0 using various line styles. The curves demonstrate how the normalized torque behaves as the size ratio approaches 0.9, revealing distinct trends based on the chosen parameter values. Notably, as ξ increases, the torque approaches a limiting value, indicating a convergence in behavior among the different parameter settings.

Figure 10 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 m
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{m}/{T}_{\infty }
                     
                   
                  versus the size ratio for various values of the length of the first couple stress parameter 
                     
                        
                        
                           ξ
                        
                        \xi 
                     
                   at 
                     
                        
                        
                           
                              
                                 
                                    
                                       a
                                    
                                    
                                       2
                                    
                                 
                                 ‒
                                 
                                    
                                       a
                                    
                                    
                                       1
                                    
                                 
                              
                              
                                 h
                              
                           
                           →
                           ∞
                           ,
                           α
                           =
                           1.0
                           ,
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                           →
                           ∞
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           ‒
                           1.0
                        
                        \frac{{a}_{2}‒{a}_{1}}{h}\to \infty ,\alpha =1.0,{\xi }^{^{\prime} }\to \infty ,{\Omega }_{2}/{\Omega }_{1}=‒1.0
                     
                  .
Figure 10

Normalized torque T m T versus the size ratio for various values of the length of the first couple stress parameter ξ at a 2 a 1 h , α = 1.0 , ξ , Ω 2 Ω 1 = 1.0 .

Figure 11 displays the normalized torque T m T 0 as a function of the size ratio Ω 2 Ω 1 for various values of the second couple stress parameter ξ . The plotted lines correspond to different parameter settings: ξ = 2.5 , 4.0 , 6.0 , 8.0 , and 10.0 . Notably, as the size ratio increases, the normalized torque exhibits an upward trend, indicating that larger size ratios lead to higher torque values. The divergence in the curves suggests that the choice of the second couple stress parameter significantly influences the relationship between torque and size ratio. As ξ increases, the torque response becomes more pronounced, particularly at larger size ratios. This behavior highlights the importance of couple stress effects in analyzing torque in systems with varying geometrical configurations.

Figure 11 
               Normalized torque 
                     
                        
                        
                           
                              
                                 T
                              
                              
                                 m
                              
                           
                           ⁄
                           
                              
                                 T
                              
                              
                                 ∞
                              
                           
                        
                        {T}_{m}/{T}_{\infty }
                     
                   
                  versus the size ratio for various values of the length of the second couple stress parameter 
                     
                        
                        
                           
                              
                                 ξ
                              
                              
                                 ′
                              
                           
                        
                        {\xi }^{^{\prime} }
                     
                   at 
                     
                        
                        
                           
                              
                                 
                                    
                                       a
                                    
                                    
                                       2
                                    
                                 
                                 ‒
                                 
                                    
                                       a
                                    
                                    
                                       1
                                    
                                 
                              
                              
                                 h
                              
                           
                           =
                           3.0
                           ,
                           α
                           =
                           1.0
                           ,
                           ξ
                           =
                           20
                           ,
                           
                              
                                 Ω
                              
                              
                                 2
                              
                           
                           ⁄
                           
                              
                                 Ω
                              
                              
                                 1
                              
                           
                           =
                           ‒
                           0.1
                        
                        \frac{{a}_{2}‒{a}_{1}}{h}=3.0,\alpha =1.0,\xi =20,{\Omega }_{2}/{\Omega }_{1}=‒0.1
                     
                  .
Figure 11

Normalized torque T m T versus the size ratio for various values of the length of the second couple stress parameter ξ at a 2 a 1 h = 3.0 , α = 1.0 , ξ = 20 , Ω 2 Ω 1 = 0.1 .

Our discussions can be summarized as the following main points:

  • The variation of normalized torque with the separation distance and the size ratio between spheres for different angular velocity ratios shows that the decreasing torque trend at higher angular velocity ratios underscores the reduced influence of couple stress effects at elevated speeds, which is critical for dynamic system applications.

  • The normalized torque as a function of parameters such as couple stress properties, size ratios, and permeability indicates that increasing couple stress parameters or adjusting size ratios significantly influences torque, offering pathways to optimize fluid behavior in engineering systems. The larger size ratios can enhance torque generation in systems with positive angular velocity ratios.

  • This diagram illustrates how torque affects with couple stress parameter within the fluid medium. For small separation distances, it arrives at the maximum value, and at a certain point, it is near to ( a 2 a 1 ) h 3 , reversing its direction to reach its minimum value. This aids in the development of predictive models and improving fluid–structure interaction analysis.

  • For negative angular velocity ratios, the torque initially decreases before increasing, indicating a complex relationship influenced by the size ratio. But, for positive angular velocity ratios, the torque consistently increases with the size ratio, highlighting the direct impact of angular velocity differences on torque generation.

  • As the first couple stress parameter ( ξ ) increases, the torque values stabilize, suggesting that higher couple stress characteristics lead to more consistent torque responses across varying size ratios.

  • The normalized torque from couple stress increases with the size ratio Ω 2 Ω 1 , indicating a stronger torque response at larger size ratios. This behavior is influenced by the second couple stress parameter ξ , with higher values leading to more pronounced torque variations.

7 Conclusions

The study investigates the torque acting on a solid sphere positioned within a spherical container filled with a porous medium saturated by incompressible couple stress fluids. It is hypothesized that the inner sphere moves along the axis of symmetry. A semi-analytical method is employed to derive the field functions and compute the torque associated with both stress and couple stress forces, with graphical representations provided. Additionally, numerical calculations of the torque, including limit cases, are presented in tabular form. The study of couple stress fluids in porous media has important applications. Furthermore, the future directions of this work could include (a) investigating the effects of varying permeability and fluid properties on stress distribution and flow behavior in couple stress media under different rotational speeds, (b) exploring the impact of eccentricity on stability and dynamic response in rotating systems, potentially incorporating non-linear fluid dynamics, and (c) developing predictive models that integrate experimental data to enhance the understanding of fluid–structure interactions in such media, informing applications in engineering and materials science. Also, the studies by Kumar and Jangili [31] and Siva et al. [32] delve into the complexities of electroosmotic flow in couple stress fluids within microchannels. While the former emphasizes heat transfer analysis in rotating systems, the latter investigates the impact of electromagnetohydrodynamics and slip-dependent zeta potentials on rotational flow dynamics. Together, they enhance the understanding of fluid behavior in microfluidic systems under coupled physical and thermal effects. On the other side [3335], these studies explore various effects on magnetohydrodynamic flows, focusing on the Soret, Joule, and Hall effects in rotating mixed convective scenarios. They investigated the behavior of elastico-viscous fluids through porous media and the implications of ion slip in unsteady conditions. These studies contribute to a deeper understanding of fluid dynamics in engineering applications involving magnetic fields. These can be taken into account in the future work.

Acknowledgments

This study was supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project Number (PNURSP2024R215), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

  1. Funding information: This study was supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project Number (PNURSP2024R215), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

  2. Author contributions: All authors were involved in the overall planning and design of the study. Shreen El-Sapa was responsible for writing, collecting data, and analyzing the results. Dr. Amal Al-Hanaya wrote the initial draft of the manuscript, which was then reviewed and improved upon by all authors. 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: All data generated or analysed during this study are included in this published article.

References

[1] Deo S, Gupta BR. Drag on a porous sphere embedded in another porous medium. J Porous Media. 2010;13(11):1009–16. 10.1615/JPorMedia.v13.i11.70. Search in Google Scholar

[2] Leontev NE. Flow past a cylinder and a sphere in a porous medium within the framework of the Brinkman equation with the Navier boundary condition. Fluid Dynam. 2014;49(2):232–7. 10.1134/S0015462814020112. Search in Google Scholar

[3] El-Sapa S. Effect of permeability of Brinkman flow on thermophoresis of a particle in a spherical cavity. Europ J Mech-B/Fluids. 2020;79:315–23. 10.1016/j.euromechflu.2019.09.017. Search in Google Scholar

[4] Faltas MS, Sherief HH, Allam AA, Ahmed BA. Mobilities of a spherical particle straddling the interface of a semi-infinite Brinkman flow. J Fluids Eng. 2021;143(7):071402. 10.1115/1.4049931. Search in Google Scholar

[5] Madasu KP, Bucha T. Influence of MHD on micropolar fluid flow past a sphere implanted in porous media. Indian J Phys. 2021;95(6):1175–83. 10.1007/s12648-020-01759-7. Search in Google Scholar

[6] Stokes VK. Couple stresses in fluids. In: Theories of fluids with microstructure: an introduction. Springer; 1984. p. 34–80. 10.1007/978-3-642-82351-0_4Search in Google Scholar

[7] Stokes VK. Theories of fluids with microstructure: an introduction. Berlin, Germany: Springer Science and Business Media; 2012. Search in Google Scholar

[8] Jangili S, Adesanya SO, Ogunseye HA, Lebelo R. Couple stress fluid flow with variable properties: a second law analysis. Math Meth Appl Sci. 2019;42(1):85–98. 10.1002/mma.5325Search in Google Scholar

[9] Al-Hanaya A, El-Sapa S, Ashmawy EA. Axisymmetric motion of an incompressible couple stress fluid between two eccentric rotating spheres. J Appl Mech Tech Phys. 2022;63(5):1–9. 10.1134/S0021894422050078Search in Google Scholar

[10] El-Sapa S, Almoneef A. Axisymmetric migration of an aerosol particle embedded in a Brinkman medium of a couple stress fluid with slip regime. Europ J Pure Appl Math. 2022;15(4):1566–92. 10.29020/nybg.ejpam.v15i4.4549Search in Google Scholar

[11] Majumdar SR. On the slow motion of viscous liquid in space between two eccentric spheres. J Phys Soc Japan. 1969;26(3):827–40. 10.1143/JPSJ.26.827Search in Google Scholar

[12] Kamel MT, Saint John NB, Chan Man Fong CF. Micropolar fluid flow between two eccentric coaxially rotating spheres. Acta Mech. 1993;99:155–71. 10.1007/BF01177242Search in Google Scholar

[13] Munson BR. Viscous incompressible flow between eccentric coaxially rotating spheres. Phys Fluids. 1974;17(3):528. 10.1063/1.1694750. Search in Google Scholar

[14] Sherief HH, Faltas MS, El-Sapa S. A general formula for the drag on a solid of revolution body at low Reynolds numbers in a microstretch fluid. Meccanica. 2017;52(11–12):2655–64. 10.1007/s11012-017-0617-1Search in Google Scholar

[15] El-Sapa S. Interaction between a non-concentric rigid sphere immersed in a micropolar fluid and a spherical envelope with slip regime. J Mol Liquids. 2022;351:118611. 10.1016/j.molliq.2022.118611Search in Google Scholar

[16] Davis AMJ. Force and torque on a rotating sphere close to and within a fluid-filled rotating sphere. In: Proceedings of the 59th Annual Meeting of the APS Division of Fluid Dynamics. American Physical Society; November 2006. Search in Google Scholar

[17] Madasu KP, Sarkar P. Torque on the slow rotation of a slightly deformed slip sphere in a Brinkman medium. Pramana. 2023;97(1):44. 10.1007/s12043-022-02505-wSearch in Google Scholar

[18] Maike S, Schröder J, Bluhm J, Ricken T. A mesh-in-element method for the theory of porous media. Z Angew Math Mech. 2024;64. 10.1002/nme.7565.Search in Google Scholar

[19] El-Sapa S, Al-Hanaya A. Effects of slippage and permeability of couple stress fluid squeezed between two concentric rotating spheres. Phys Fluids. 2023;35:103112. 10.1063/5.0171851Search in Google Scholar

[20] Yadav PK, Roshan M. Mathematical modeling of blood flow in an annulus porous region between two coaxial deformable tubes: An advancement to peristaltic endoscope. Chin J Phys. 2024;88:89–109. 10.1016/j.cjph.2024.01.017Search in Google Scholar

[21] Yadav PK, Roshan M. Mathematical modeling of creeping electromagnetohydrodynamic peristaltic propulsion in an annular gap between sinusoidally deforming permeable and impermeable curved tubes. Phys Fluids. 2024;36(7):071907. 10.1063/5.0217370. Search in Google Scholar

[22] Sharma R, Thakur K. On couple-stress fluid heated from below in porous medium in hydromagnetics. Czechoslovak J Phys. 2000;50:753–8. 10.1023/A:1022886903213Search in Google Scholar

[23] Chand R, Rana G, Yadav D. Thermal instability in a layer of couple stress nanofluid saturated porous medium. J Theoretic Appl Mech. 2017;47:69. 10.1515/jtam-2017-0005Search in Google Scholar

[24] Ochoa-Tapia JA, Whitaker S. Momentum transfer at the boundary between a porous medium and a homogeneous fluid- II. Comparison with experiment. Int J Heat Mass Transfer. 1995;38(14):2647–55. 10.1016/0017-9310(94)00347-XSearch in Google Scholar

[25] Ochoa-Tapia JA, Whitaker S. Momentum transfer at the boundary between a porous medium and a homogeneous fluid-I. Theoretical development. Int J Heat Mass Transfer. 1995;38(14):2635–46. 10.1016/0017-9310(94)00346-WSearch in Google Scholar

[26] Madasu KP, Sarkar P. Couple stress fluid past a sphere embedded in a porous medium. Archive Mech Eng. 2022;5–19. 10.24425/ame.2021.139314Search in Google Scholar

[27] Ramkissoon H. On the rotation of an axially symmetric body in couple stress fluids. Z Angew Math Mech. 1984;64:297–302. 10.1002/zamm.19840640712Search in Google Scholar

[28] Faltas MS, Nashwan MG, Ahmed BA, Hamdy P. Oscillations of a spherical particle in the presence of a flat interface separating two fluid phases. Chinese J Phys. 2024;89:1377–99. 10.1016/j.cjph.2024.04.014Search in Google Scholar

[29] Faltas MS, Sherief HH, El-Maghraby NM, Wanas EF. The electrokinetic flow of a micropolar fluid in a microtube with velocity and spin velocity slippage. Chinese J Phys. 2024;89:504–27. 10.1016/j.cjph.2023.10.034Search in Google Scholar

[30] Faltas MS, Ashmawy EA, Sherief HH, Othman HA. Oscillations of coaxial hydrophobic spherical colloidal particles in a micropolar fluid. Phys Fluids. 2024;36(5):052017. 10.1063/5.0208069Search in Google Scholar

[31] Kumar B, Jangili S. Heat transfer analysis of electroosmotic flow of couple stress fluid through a rotating circular microchannel. Int J Thermal Sci. February 2025;208:109359. 10.1016/j.ijthermalsci.2024.109359Search in Google Scholar

[32] Siva T, Dubey D, Jangili S. Rotational flow dynamics of electroosmotic transport of couple stress fluid in a microfluidic channel under electromagnetohydrodynamic and slip-dependent zeta potential effects. Phys Fluids. 2024;36(9):092006. 10.1063/5.0215302. Search in Google Scholar

[33] Veera Krishna M, Swarnalathamma BV, Chamkha AJ. Investigations of Soret. Joule and Hall effects on MHD rotating mixed convective flow past an infinite vertical porous plate. J Ocean Eng Sci. 2019;4(3):263–75. 10.1016/j.joes.2019.05.002Search in Google Scholar

[34] Veera Krishna M, Chamkha AJ. Hall and ion slip effects on MHD rotating flow of elastico-viscous fluid through porous medium. Int Commun Heat Mass Transfer. 2020;113:104494. 10.1016/j.icheatmasstransfer.2020.104494Search in Google Scholar

[35] Veera Krishna M, Ameer Ahamad N, Chamkha AJ. Hall and ion slip effects on unsteady MHD free convective rotating flow through a saturated porous medium over an exponential accelerated plate. Alexandr Eng J. 2020;59(2):565–77. 10.1016/j.aej.2020.01.043Search in Google Scholar

Received: 2024-09-30
Revised: 2024-12-08
Accepted: 2024-12-25
Published Online: 2024-12-31

© 2024 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. Regular Articles
  2. Numerical study of flow and heat transfer in the channel of panel-type radiator with semi-detached inclined trapezoidal wing vortex generators
  3. Homogeneous–heterogeneous reactions in the colloidal investigation of Casson fluid
  4. High-speed mid-infrared Mach–Zehnder electro-optical modulators in lithium niobate thin film on sapphire
  5. Numerical analysis of dengue transmission model using Caputo–Fabrizio fractional derivative
  6. Mononuclear nanofluids undergoing convective heating across a stretching sheet and undergoing MHD flow in three dimensions: Potential industrial applications
  7. Heat transfer characteristics of cobalt ferrite nanoparticles scattered in sodium alginate-based non-Newtonian nanofluid over a stretching/shrinking horizontal plane surface
  8. The electrically conducting water-based nanofluid flow containing titanium and aluminum alloys over a rotating disk surface with nonlinear thermal radiation: A numerical analysis
  9. Growth, characterization, and anti-bacterial activity of l-methionine supplemented with sulphamic acid single crystals
  10. A numerical analysis of the blood-based Casson hybrid nanofluid flow past a convectively heated surface embedded in a porous medium
  11. Optoelectronic–thermomagnetic effect of a microelongated non-local rotating semiconductor heated by pulsed laser with varying thermal conductivity
  12. Thermal proficiency of magnetized and radiative cross-ternary hybrid nanofluid flow induced by a vertical cylinder
  13. Enhanced heat transfer and fluid motion in 3D nanofluid with anisotropic slip and magnetic field
  14. Numerical analysis of thermophoretic particle deposition on 3D Casson nanofluid: Artificial neural networks-based Levenberg–Marquardt algorithm
  15. Analyzing fuzzy fractional Degasperis–Procesi and Camassa–Holm equations with the Atangana–Baleanu operator
  16. Bayesian estimation of equipment reliability with normal-type life distribution based on multiple batch tests
  17. Chaotic control problem of BEC system based on Hartree–Fock mean field theory
  18. Optimized framework numerical solution for swirling hybrid nanofluid flow with silver/gold nanoparticles on a stretching cylinder with heat source/sink and reactive agents
  19. Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations
  20. Convective flow of a magnetohydrodynamic second-grade fluid past a stretching surface with Cattaneo–Christov heat and mass flux model
  21. Analysis of the heat transfer enhancement in water-based micropolar hybrid nanofluid flow over a vertical flat surface
  22. Microscopic seepage simulation of gas and water in shale pores and slits based on VOF
  23. Model of conversion of flow from confined to unconfined aquifers with stochastic approach
  24. Study of fractional variable-order lymphatic filariasis infection model
  25. Soliton, quasi-soliton, and their interaction solutions of a nonlinear (2 + 1)-dimensional ZK–mZK–BBM equation for gravity waves
  26. Application of conserved quantities using the formal Lagrangian of a nonlinear integro partial differential equation through optimal system of one-dimensional subalgebras in physics and engineering
  27. Nonlinear fractional-order differential equations: New closed-form traveling-wave solutions
  28. Sixth-kind Chebyshev polynomials technique to numerically treat the dissipative viscoelastic fluid flow in the rheology of Cattaneo–Christov model
  29. Some transforms, Riemann–Liouville fractional operators, and applications of newly extended M–L (p, s, k) function
  30. Magnetohydrodynamic water-based hybrid nanofluid flow comprising diamond and copper nanoparticles on a stretching sheet with slips constraints
  31. Super-resolution reconstruction method of the optical synthetic aperture image using generative adversarial network
  32. A two-stage framework for predicting the remaining useful life of bearings
  33. Influence of variable fluid properties on mixed convective Darcy–Forchheimer flow relation over a surface with Soret and Dufour spectacle
  34. Inclined surface mixed convection flow of viscous fluid with porous medium and Soret effects
  35. Exact solutions to vorticity of the fractional nonuniform Poiseuille flows
  36. In silico modified UV spectrophotometric approaches to resolve overlapped spectra for quality control of rosuvastatin and teneligliptin formulation
  37. Numerical simulations for fractional Hirota–Satsuma coupled Korteweg–de Vries systems
  38. Substituent effect on the electronic and optical properties of newly designed pyrrole derivatives using density functional theory
  39. A comparative analysis of shielding effectiveness in glass and concrete containers
  40. Numerical analysis of the MHD Williamson nanofluid flow over a nonlinear stretching sheet through a Darcy porous medium: Modeling and simulation
  41. Analytical and numerical investigation for viscoelastic fluid with heat transfer analysis during rollover-web coating phenomena
  42. Influence of variable viscosity on existing sheet thickness in the calendering of non-isothermal viscoelastic materials
  43. Analysis of nonlinear fractional-order Fisher equation using two reliable techniques
  44. Comparison of plan quality and robustness using VMAT and IMRT for breast cancer
  45. Radiative nanofluid flow over a slender stretching Riga plate under the impact of exponential heat source/sink
  46. Numerical investigation of acoustic streaming vortices in cylindrical tube arrays
  47. Numerical study of blood-based MHD tangent hyperbolic hybrid nanofluid flow over a permeable stretching sheet with variable thermal conductivity and cross-diffusion
  48. Fractional view analytical analysis of generalized regularized long wave equation
  49. Dynamic simulation of non-Newtonian boundary layer flow: An enhanced exponential time integrator approach with spatially and temporally variable heat sources
  50. Inclined magnetized infinite shear rate viscosity of non-Newtonian tetra hybrid nanofluid in stenosed artery with non-uniform heat sink/source
  51. Estimation of monotone α-quantile of past lifetime function with application
  52. Numerical simulation for the slip impacts on the radiative nanofluid flow over a stretched surface with nonuniform heat generation and viscous dissipation
  53. Study of fractional telegraph equation via Shehu homotopy perturbation method
  54. An investigation into the impact of thermal radiation and chemical reactions on the flow through porous media of a Casson hybrid nanofluid including unstable mixed convection with stretched sheet in the presence of thermophoresis and Brownian motion
  55. Establishing breather and N-soliton solutions for conformable Klein–Gordon equation
  56. An electro-optic half subtractor from a silicon-based hybrid surface plasmon polariton waveguide
  57. CFD analysis of particle shape and Reynolds number on heat transfer characteristics of nanofluid in heated tube
  58. Abundant exact traveling wave solutions and modulation instability analysis to the generalized Hirota–Satsuma–Ito equation
  59. A short report on a probability-based interpretation of quantum mechanics
  60. Study on cavitation and pulsation characteristics of a novel rotor-radial groove hydrodynamic cavitation reactor
  61. Optimizing heat transport in a permeable cavity with an isothermal solid block: Influence of nanoparticles volume fraction and wall velocity ratio
  62. Linear instability of the vertical throughflow in a porous layer saturated by a power-law fluid with variable gravity effect
  63. Thermal analysis of generalized Cattaneo–Christov theories in Burgers nanofluid in the presence of thermo-diffusion effects and variable thermal conductivity
  64. A new benchmark for camouflaged object detection: RGB-D camouflaged object detection dataset
  65. Effect of electron temperature and concentration on production of hydroxyl radical and nitric oxide in atmospheric pressure low-temperature helium plasma jet: Swarm analysis and global model investigation
  66. Double diffusion convection of Maxwell–Cattaneo fluids in a vertical slot
  67. Thermal analysis of extended surfaces using deep neural networks
  68. Steady-state thermodynamic process in multilayered heterogeneous cylinder
  69. Multiresponse optimisation and process capability analysis of chemical vapour jet machining for the acrylonitrile butadiene styrene polymer: Unveiling the morphology
  70. Modeling monkeypox virus transmission: Stability analysis and comparison of analytical techniques
  71. Fourier spectral method for the fractional-in-space coupled Whitham–Broer–Kaup equations on unbounded domain
  72. The chaotic behavior and traveling wave solutions of the conformable extended Korteweg–de-Vries model
  73. Research on optimization of combustor liner structure based on arc-shaped slot hole
  74. Construction of M-shaped solitons for a modified regularized long-wave equation via Hirota's bilinear method
  75. Effectiveness of microwave ablation using two simultaneous antennas for liver malignancy treatment
  76. Discussion on optical solitons, sensitivity and qualitative analysis to a fractional model of ion sound and Langmuir waves with Atangana Baleanu derivatives
  77. Reliability of two-dimensional steady magnetized Jeffery fluid over shrinking sheet with chemical effect
  78. Generalized model of thermoelasticity associated with fractional time-derivative operators and its applications to non-simple elastic materials
  79. Migration of two rigid spheres translating within an infinite couple stress fluid under the impact of magnetic field
  80. A comparative investigation of neutron and gamma radiation interaction properties of zircaloy-2 and zircaloy-4 with consideration of mechanical properties
  81. New optical stochastic solutions for the Schrödinger equation with multiplicative Wiener process/random variable coefficients using two different methods
  82. Physical aspects of quantile residual lifetime sequence
  83. Synthesis, structure, IV characteristics, and optical properties of chromium oxide thin films for optoelectronic applications
  84. Smart mathematically filtered UV spectroscopic methods for quality assurance of rosuvastatin and valsartan from formulation
  85. A novel investigation into time-fractional multi-dimensional Navier–Stokes equations within Aboodh transform
  86. Homotopic dynamic solution of hydrodynamic nonlinear natural convection containing superhydrophobicity and isothermally heated parallel plate with hybrid nanoparticles
  87. A novel tetra hybrid bio-nanofluid model with stenosed artery
  88. Propagation of traveling wave solution of the strain wave equation in microcrystalline materials
  89. Innovative analysis to the time-fractional q-deformed tanh-Gordon equation via modified double Laplace transform method
  90. A new investigation of the extended Sakovich equation for abundant soliton solution in industrial engineering via two efficient techniques
  91. New soliton solutions of the conformable time fractional Drinfel'd–Sokolov–Wilson equation based on the complete discriminant system method
  92. Irradiation of hydrophilic acrylic intraocular lenses by a 365 nm UV lamp
  93. Inflation and the principle of equivalence
  94. The use of a supercontinuum light source for the characterization of passive fiber optic components
  95. Optical solitons to the fractional Kundu–Mukherjee–Naskar equation with time-dependent coefficients
  96. A promising photocathode for green hydrogen generation from sanitation water without external sacrificing agent: silver-silver oxide/poly(1H-pyrrole) dendritic nanocomposite seeded on poly-1H pyrrole film
  97. Photon balance in the fiber laser model
  98. Propagation of optical spatial solitons in nematic liquid crystals with quadruple power law of nonlinearity appears in fluid mechanics
  99. Theoretical investigation and sensitivity analysis of non-Newtonian fluid during roll coating process by response surface methodology
  100. Utilizing slip conditions on transport phenomena of heat energy with dust and tiny nanoparticles over a wedge
  101. Bismuthyl chloride/poly(m-toluidine) nanocomposite seeded on poly-1H pyrrole: Photocathode for green hydrogen generation
  102. Infrared thermography based fault diagnosis of diesel engines using convolutional neural network and image enhancement
  103. On some solitary wave solutions of the Estevez--Mansfield--Clarkson equation with conformable fractional derivatives in time
  104. Impact of permeability and fluid parameters in couple stress media on rotating eccentric spheres
  105. Review Article
  106. Transformer-based intelligent fault diagnosis methods of mechanical equipment: A survey
  107. Special Issue on Predicting pattern alterations in nature - Part II
  108. A comparative study of Bagley–Torvik equation under nonsingular kernel derivatives using Weeks method
  109. On the existence and numerical simulation of Cholera epidemic model
  110. Numerical solutions of generalized Atangana–Baleanu time-fractional FitzHugh–Nagumo equation using cubic B-spline functions
  111. Dynamic properties of the multimalware attacks in wireless sensor networks: Fractional derivative analysis of wireless sensor networks
  112. Prediction of COVID-19 spread with models in different patterns: A case study of Russia
  113. Study of chronic myeloid leukemia with T-cell under fractal-fractional order model
  114. Accumulation process in the environment for a generalized mass transport system
  115. Analysis of a generalized proportional fractional stochastic differential equation incorporating Carathéodory's approximation and applications
  116. Special Issue on Nanomaterial utilization and structural optimization - Part II
  117. Numerical study on flow and heat transfer performance of a spiral-wound heat exchanger for natural gas
  118. Study of ultrasonic influence on heat transfer and resistance performance of round tube with twisted belt
  119. Numerical study on bionic airfoil fins used in printed circuit plate heat exchanger
  120. Improving heat transfer efficiency via optimization and sensitivity assessment in hybrid nanofluid flow with variable magnetism using the Yamada–Ota model
  121. Special Issue on Nanofluids: Synthesis, Characterization, and Applications
  122. Exact solutions of a class of generalized nanofluidic models
  123. Stability enhancement of Al2O3, ZnO, and TiO2 binary nanofluids for heat transfer applications
  124. Thermal transport energy performance on tangent hyperbolic hybrid nanofluids and their implementation in concentrated solar aircraft wings
  125. Studying nonlinear vibration analysis of nanoelectro-mechanical resonators via analytical computational method
  126. Numerical analysis of non-linear radiative Casson fluids containing CNTs having length and radius over permeable moving plate
  127. Two-phase numerical simulation of thermal and solutal transport exploration of a non-Newtonian nanomaterial flow past a stretching surface with chemical reaction
  128. Natural convection and flow patterns of Cu–water nanofluids in hexagonal cavity: A novel thermal case study
  129. Solitonic solutions and study of nonlinear wave dynamics in a Murnaghan hyperelastic circular pipe
  130. Comparative study of couple stress fluid flow using OHAM and NIM
  131. Utilization of OHAM to investigate entropy generation with a temperature-dependent thermal conductivity model in hybrid nanofluid using the radiation phenomenon
  132. Slip effects on magnetized radiatively hybridized ferrofluid flow with acute magnetic force over shrinking/stretching surface
  133. Significance of 3D rectangular closed domain filled with charged particles and nanoparticles engaging finite element methodology
  134. Robustness and dynamical features of fractional difference spacecraft model with Mittag–Leffler stability
  135. Characterizing magnetohydrodynamic effects on developed nanofluid flow in an obstructed vertical duct under constant pressure gradient
  136. Study on dynamic and static tensile and puncture-resistant mechanical properties of impregnated STF multi-dimensional structure Kevlar fiber reinforced composites
  137. Thermosolutal Marangoni convective flow of MHD tangent hyperbolic hybrid nanofluids with elastic deformation and heat source
  138. Investigation of convective heat transport in a Carreau hybrid nanofluid between two stretchable rotatory disks
  139. Single-channel cooling system design by using perforated porous insert and modeling with POD for double conductive panel
  140. Special Issue on Fundamental Physics from Atoms to Cosmos - Part I
  141. Pulsed excitation of a quantum oscillator: A model accounting for damping
  142. Review of recent analytical advances in the spectroscopy of hydrogenic lines in plasmas
  143. Heavy mesons mass spectroscopy under a spin-dependent Cornell potential within the framework of the spinless Salpeter equation
  144. Coherent manipulation of bright and dark solitons of reflection and transmission pulses through sodium atomic medium
  145. Effect of the gravitational field strength on the rate of chemical reactions
  146. The kinetic relativity theory – hiding in plain sight
  147. Special Issue on Advanced Energy Materials - Part III
  148. Eco-friendly graphitic carbon nitride–poly(1H pyrrole) nanocomposite: A photocathode for green hydrogen production, paving the way for commercial applications
Downloaded on 10.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2024-0112/html
Scroll to top button