Theoretical Analysis of Squeeze Lubrication Using Double ZZ Transform: Application of Non-Newtonian Fluids

ABSTRACT


INTRODUCTION
The idea of "integral transforms give us power and ease to solve such fundamental problems as initial value ones of a linear and of an integral differential equation, and numerous mathematicians across the world became interested in this problem during the last decades" of their applications of fluid mechanics was of great interest to the authors such as Sumudu, ARA, Saban, SEA, and SEJI [1][2][3][4][5][6][7].
Most problems in applied science and engineering fall under double integral equations or partial differential equations that are used to describe the nature and behavior of the physical outcome.Finally, the researcher concluded that, two years ago, 2016, Ranjit and Waghmare [8] hard work led to developing single transform to double transform partial equations that help in most mathematics problems to date.Ranjit and Waghmare 2016 applied the double Laplace double transform method to get the linear/nonlinear space-time exact solutions of the time-fractional telegraph when they considered t double Laplace finds the solution.Conclusion Ranjit and Waghmare t double Laplace transform method reduces the computational work volume since it is easier in solving partial equations as compared to single Laplace because the solution to the problems does not require an A domain polynomial, Lagrange multiplier value, He's polynomials, and small parameters.
Especially, compared to the numerical ones, the double Laplace transform has quite common use to solve the partial differential equations with unknown functions of two variables, which is general, and less easy to check in a numerical way.Therefore, it is highly advisable to check the performance of numerical methods to find the model error with analytical ones, which have developed and become more available.One of highly effective such analytical methods is the integral transformations, which are used to find analytical solutions of partial differential equations.Apart from the discussed double Laplace transform, there are other applications in the literature as the double extensions of the double Laplace transform and double Shuhu and double ARA-Sumudu transforms [9][10][11].
In 2016, Zafar [12] introduced a new integral transform 'ZZ Transform Method'' and it was used for the new integral transform technique application to solve first order linear differential equation, for example, Law of Natural Growth or Law of Natural Decay.ZZ transform Method is used for differential equation with variable coefficients.In 2022 a new transform, the double ZZ transform, was developed by researchers to extend the one-dimensional ZZ transform to two dimensions [13,14].We have also derived a few important theorems and properties.this method the Partial Differential equation is resolved without converting it to ordinary differential equation i.e., it is not necessary to determine the complete solution of Ordinary differential equation.This is the most important advantage of this method.Thus, it is very convenient & effective method.To demonstrate the effectiveness and great accuracy of the suggested transform, the transform is employed to find a solution to the popular squeeze action issue for unstable and incompressible fluid flow through the porous elastic of layers.By pressing the synovial fluid, the squeeze lubrication generates the hydrodynamic pressure.Both human articular mobility and synovial joints are required.The interstitial fluid must be removed from the solid organic matrix when the cartilages are squeezed under physiologically normal conditions (see Figure 1).The synovial fluid causes the interstitial fluid of the solid organic matrix to move relatively and to get drained by convection and diffusion [15,16].The amount of hyaluronic acid in the synovial fluid rises with lubrication and therefore imbibition and excretion of the solid mass increases its concentration [17,18].As cited by Wegamir increasing the concentration of Hyaluronic acid increases the viscosity of synovial fluid.Squeeze lubrication occurs when the bearing surfaces move at right angles to each other.The pressure in the fluid film is generated by the movement of the articular surfaces parallel to each other.As the approaching surfaces come closer, the competing surfaces force the air out of the area of contact.The pressure is caused by the attraction force that exists between the surfaces in the field of the lubricant, which persists as the fluid film is compressed.The squeeze lubrication is useful in short, high-pressure utilization.The demand of this device comes from a frequent pressure generation [20].
The pressure yielded in the femur-humerus and between the head of the femur and the tibia on the intervention of the synovial fluid is an example.It is significant because it is more directly related to the gait phase while it can only sustain twice the weight of the individual due to movement.The lubricating fluid is flowing through the pores of the articular cartilage to the lubricated when the movements are taking place.The amount of fluid remains constant, but the pressure emerged from the action of one layer of fluid compressing one another varies centuries.

BASIC DEFINITION OF FRACTIONAL CALCULUS
In this section, the study provides some basic definitions and theories of the double ZZ transform that will be used in this article:
Table 1 shows double ZZ transform for some special functions [13].

ASSUMPTIONS OF SQUEEZE LUBRICATION
Some of assumptions have been adopted in the current study: (1) Non-Newtonian fluid (2) Incompressible squeeze lubrication fluid (3) Unsteady (4) Neglect bodily forces like gravity and magnetic field (5) In fluid dynamics, the no-slip condition for viscous fluids assumes that at a solid boundary, the fluid will have zero velocity relative to the boundary.The fluid velocity at all fluid-solid boundaries is equal to that of the solid boundary.

Basic equation
The main equations that describe the theoretical analysis of the characteristics of the squeeze pressure region of the synovial fluid for the human joint can be derived and be a function of the phases of movement and gait of the synovial joint and the mathematical model and result analysis of the squeeze lubrication when the above-mentioned boundary conditions are addressed.This is a deterministic analysis and the effects of squeeze lubrication and the pressure and the friction force between the synovial fluid particles and the synovial fluid particles and layers against themselves are considered simultaneously.

Mathematical formulation of the problem
The continuity equation simply expresses the law of conservation of mass (mass per unit time entering the tube must flow out at same rate).The equation of continuity in cylindrical coordinates when the fluid is non-compressible, is written as follows: The equations of motion of a real fluid can be developed from consideration of the force action on a small element of the fluid including the shear stresses generated by fluid motion and viscosity.These equations are called Navier-Stokes is written as follows: After applying the assumptions of squeeze lubrication, the equations responsible the flow with (turbulentregular) through porous elastic layers become follows: Now, modifying the film region Eq. ( 3) of the Navier-Stokes equation to include the tangential component of fluid velocity: Then: , we can achieve the porosity non-Newtonian: (1 +  +  +  0 ) (

Squeeze film pressure
For theoretical and computational reasons, it is crucial to include the non-dimensional parameters in the pressure's governing equations.The different lubricating system parameters should be presented in non-dimensional form as well.
After that:

Friction force
The friction force between layers the particles of synovial fluid is very important factor.We will study set of internal factors that affect the force of friction, assume that non-Newtonian fluid and Newton low of viscosity be: where, μ is dynamic viscosity and the term (

Squeeze pressure
The variation of the dimensionless squeeze film pressure (p * ) generated by the squeeze film action as a function of dimensionless distance (r * ) for different values of fluid velocity ( ⃗ ⃗ ) is shown in Figure 2 and Table 2 with using Eq.(20).It is observed increase in the velocity of the synovial fluid particles leads to the movement of the particles, which leads to an increase in the rate of collision of the particles and the viscosity of the synovial fluid increases, and thus the pressure increases.

Table 2. Relationship between velocity of particles and pressure
After 20-Cycle Time 1.9  The different dimensionless pressure distribution for different value weight of human is seen in Figure 3 and Table 3.It is found that the pressure increases with the increase in weight of the human and time of cycle, and the explanation is that the longer the cycle, which leads to an increase in the rate of flow of lubricant from the porosity and with the different activities performed by the person, the weight doubles and the production of lubricants decreases.
The variation of the dimensionless squeeze film pressure (p * ) generated by the squeeze film action as a function of dimensionless distance (r * ) for different values of fluid concentration (C0) is seen in Figure 4 and Table 4.It is cleared the dimensionless pressure increasing and become more with decreasing value of synovial fluid concentration (C0=0.2) inversely, increasing fluid initial concentration (C0=0.6)that lead to decreasing in dimensionless pressure.The relationship between dimensionless hydrodynamic pressure (p * ) with different porosity of articular cartilage (β) is shown in Figure 5.As a result, it can be seen that reducing the porosity causes additional outflow of synovial fluid from synovial cells, which, in turn, increases the hydrodynamic pressure between layers.At the brushes in non-Newtonian lubrication for film thickness involves the cycle time becomes a minimum film thickness, and again, time is necessary.The effect of a film thickness parameter (H) of gab between two articulates on the variation (p * ) is demonstrated in Figure 6.

Kinematic friction force
Kinetic friction force is the force that occurs when resisting motion due to contact between a layer and a particle that is moving against it.The current study has dealt with a set of internal factors that affect the force of friction and after applying Eq. ( 29) to these factors using a 12 mathematical program.Dimensional friction force (ξ * ) as a function of the gait phase (ω) for different values of different viscosity (), high viscosity means an increase in the cohesion force between the molecules, which leads to a high kinetic friction force between molecules.The friction force is related to the phases of movement in a direct relationship as shown in the Table 5 and Figure 8. Dimensional friction force (ξ * ) as a function of the gait phase (ω) for different initial concentration (c0).The force of adhesion between particles when it is high, the concentration of the liquid is high, especially in the phase of heel contact as shown Figure 9.   Dimensional friction force (ξ * ) as a function of the gait phase (ω) for different nature of the contacting surfaces (δ),the nature of the surface: Friction depends on the nature of the contacting surfaces, as rough surfaces will need more force to move them than if those surfaces are smooth, as the force of friction decreases with the increase in the smoothness of the body to a certain degree, but if it exceeds that degree, the friction actually increases between two surfaces They are very smooth due to the increased electrostatic forces between their molecules as shown in the Table 6 and Figure 10 .
Dimensional friction force ( * ) as a function of the gait phase (ω) for regular velocity (  ).The lubricated liquid contains particles.The cohesion force between the particles increases with the increase in the flow velocity between the layers (superficial zone-middle zone-deep zone) and the walking stages, resulting in a high friction force of up to 80%.as shown in Figures 11 and 12.The lubricated fluid contains particles that have a high percentage of viscosity depending on the pattern of movement.The flow of the fluid loaded with particles increases as the speed increases.The friction between the particles increases with it.The rate of friction is higher in (superficial zone).
The force of friction does not depend on the area of contacting surfaces in moving objects, or on the relative velocity of those objects, but rather on the interconnecting forces arising between them, as shown in Figure 13.The force of friction is closely related to the phase of movement.In the phase (loading response), the flow of fluid increases through the porosity of the cartilage, which generates pressure between the layers that reduces friction.In the phase (early mid stance), the weight doubles and the flow of lubricants decreases, which makes the force of friction between the layers higher.In the phase (terminal), it increases in strength.Friction greatly affects the performance of the joint, and to reduce damage, we raise the fitness parameters, which are the focus of our study.

CONCLUSIONS
In this paper, we investigate the flow of non-Newtonian fluids through the porous elastic layers of human articular cartilage and describe the impact of lubricating fluid flow on performance.In this research, we highlight (squeeze lubrication) the lubricated liquid and identify the most important parameters that affect the flow of the liquid and which are reflected in the movement mechanism, which represents an important element for performing daily activities.To accomplish these goals, the characteristics of lubrication through the cohesive force of particles, initial concentration, surface smoothness, vertical roughness, and human body weight are studied.These expressions are computed using the double ZZ transformation and are employed to describe numerous scientific problems.The current findings indicate that using single transforms to solve these equations is more difficult than using the double transform.In addition, both the theory and practise of lubrication are evident when hydrodynamic pressure is generated between distinct layers.This pressure contributes significantly to the kinematic force of friction between strata and particles.Kinetic friction between the layers is greatly affected by pressure, as high pressure generates high stress peaks that protect the layers and provide safety for the cartilage, and the opposite happens when pressure decreases.According to the analysis and computation of the results, the pressure and friction forces increase as the cycle time increases.Dynamic friction and pressure play a major role in the stages of walking.When the activities are recording, the applied pressure is increased.It is appropriate to learn the patterns of movement to reduce the harm, and the developed laws must define what combination of parameters influences on directions.It is possible to create doubles of ZZ applications and resolve coupled the differential equations and systems and PDEs with variable coefficients.The mentioned new double transform presents additional novelties comparable with the existing alternative methods as follows: The numerous advantages of our double ZZ sequence are that it can be used with a numerical iterative solution of nonlinear problems of PDEs.

Figure 10 .
Figure 10.Variation of dimensionless Kinematic friction force (ξ * ) with gait phase (ω) for different the nature of the contacting surfaces: (a) surface roughness; (b) surface smooth

Figure 16 .
Figure 16.Variation of dimensionless Kinematic friction force (ξ * ) with gait phase (ω) for different weight of human

Figure 17 .
Figure 17.Variation of dimensionless Kinematic friction force (ξ * ) with gait phase (ω) for different flow force Dimensional friction force (ξ * ) as a function of the gait phase (ω) for Kinematic friction force.The relationship between the kinetic friction force and the fluid flow force is a direct relationship, the greater the value of the fluid flow force, the greater the kinetic friction force, and the greater the body's resistance to movement, as shown in Figure 17.

Table 1 .
Double ZZ transform for some special functions [

Table 3 .
Relationship between weight of human and pressure

Table 4 .
Relationship between initial constriction and pressure

Table 6 .
Relationship between parametric of the nature of the contacting surfaces: (a) surface roughness; (b) surface smooth (ζ) and Kinematic friction force (ξ * )