Impact of Moving Walls and Entropy Generation on Doubly Diffusive Mixed Convection of Casson Fluid in Two-Sided Driven Enclosure

In this study, numerical simulations are conducted with the goal of exploring the impact of the direction of the moving wall, solute and thermal transport, and entropy production on doubly diffusive convection in a chamber occupied by a Casson liquid. Wall movement has a significant impact on convective flow, which, in turn, affects the rate of mass and heat transfer; this sparked our interest in conducting further analysis. The left and right (upright) walls are preserved with constant (but different) thermal and solutal distributions, while the horizontal boundaries are impermeable to mass transfer and insulated from heat transfer. Numerical solutions are acquired using the control volume technique. Outcomes under a variety of Casson fluid parameters, including Ri, Gr, buoyancy ratio, and direction of the moving wall(s), are explored, and the influences of entropy generation are comprehensively investigated. While the flow field consists of a single cell in case I, it is dual-cellular in case III for all values of the considered parameters. Comparing the three cases, the average heat and mass transport presented lower values in case III due to the movement of an isothermal (left) wall against the buoyant force, while these values are enhanced in case I. The obtained results are expected to be useful in thermal engineering, material, food, and chemical processing applications.


Introduction
The convective movement of non-Newtonian-type fluids in various geometries is important due to their applications in daily life and scientific requirements [1][2][3][4].A Casson fluid (CF) is a type of non-Newtonian liquid; in particular, CF presents a shear-thinning type of non-Newtonian behavior [5,6].The thermal mixing of non-Newtonian liquids is essential in food processing and chemical processes [1,7].Additionally, external forces, such as wall movement, dynamically influence the flow field in various fields [8][9][10][11][12].Hence, in this study, both mechanisms are incorporated to observe the flow and the thermal and solutal fields of a Casson liquid in a closed chamber.This survey initially describes the convective flow of Casson fluids; second, single or double lid-driven cavity (LDC) flows are detailed; and finally, the entropy generation (EG) of convection flows in the enclosure is discussed.
The convective flow of non-Newtonian liquids in closed chambers has been investigated in several studies performed throughout the past decade [13][14][15][16][17][18][19][20].Pop and Sheremet [21] studied the radiative and buoyant convection of CF in a bounded domain with the existence of viscous dissipation.They proved that the CF parameter favors stream intensification and heat transfer augmentation.Sheremet and Pop [22] studied the (thermal) radiation and convective current of a CF in a square chamber, and achieved thermal energy drops after raising the elastic number.Hamid et al. [23] studied the free convective stream of CF in a partly heated trapezoidal chamber using the Galerkin-type finite element method (GFEM), and detected that the CF parameter decreased the velocity.Hussain et al. [24] have determined the effects of entropy generation and magnetic strength on the doubly diffusive convection of CF fluid in a staggered cavity using GFEM.The magnetic convection energy transport of CF driven by combined solutal and thermal diffusion in a trapezoidal enclosure has been explored by Shahzad et al. [25], where the impact of a tilted magnetic field ruled by the Lorentz field law force was considered.Using FEM, Aghighi et al. [26] carried out numerical research on doubly diffusive convective movement of a CF based on a viscoplastic stress model in a square-shaped chamber.They deduced that the buoyancy ratio increases the extreme yield stress by up to 400%, and that the Lewis number provides a maximum value of only 20%.Hussain et al. [27] have studied EG for the free convection of CF in a partially bottom-heated (square-shaped) porous domain with horizontal magnetic force strength, inclined angle, and viscous dissipation.GFEM was used to obtain the solution.They found that the entire entropy production rises when increasing the inclination of the domain.Hirpho and Ibrahim [28] examined the combined convection of a hybrid Casson nanofluid (CNF) in a partly heated cavity using FEM.They found that the thermal transport behavior strongly varied with an increase in the Richardson numbers.Using GFEM, Vishnu Ganesh et al. [29] investigated convective flow CNF in a complex wavy chamber with an inner solid body.Patel et al. [30] investigated thermal radiation and Dufour and Soret effects on Casson nanoliquid flow in a channel with entropy generation.
The flow due to a moving lid or wall of a cavity together with heat and solute transfer is an important problem in the area of fluid dynamics and thermal sciences, due to its extensive range of applications in metal casting, mixing, drying dynamics of lakes, cooling areas, and so on [31][32][33][34].Some studies on a lid-driven cavity (LDC) have been performed by Sivasankaran et al. [8][9][10][11][12][31][32][33][34]; however, these studies were limited to only single-wall (top or side) motion of the cavities.Two-sided lid-driven cavity (2SLDC) flow has raised interest due to its applications in various fields, such as material processing, coating, and cooling.Some notable investigations focused on 2SLDC can be found in [35].Multiple solutions of convective flow in a 2SLDC with an aspect ratio of 1.96 were examined by Luo and Yang [36].Chen et al. [37] explored the effect of parallel and antiparallel motions of the (horizontal) walls of 2SLDC on steady convective flow.Asiaei et al. [38] explored the impact of multi-layered porous materials on the combined convective flow of nanoliquid in a 2SLDC, and also examined the effect of interior heat generation and EG.They found that greater thermal transport yields higher values of the entropy generation number.Maghsoudi and Siavashi [39] optimized the pore size of the porous structure in a mixed convection nanoliquid in a 2SLDC to enhance energy transport.Bhopalam et al. [40] studied the incompressible flow in an oscillating 2SLDC, and established that the dynamics of the vortex in oscillatory LDC flow is more complex than that in steady LDC flow.Azzouz and Houat [41] explored the complex topological flow pattern of a three-sided lid-driven cavity (3SLDC), showing that three types of direction support multiple alterations in the flow pattern.The LDC flow in a cavity with a top-lid, divided into two parts that could act as stirrers during a chemical mixing operation, was investigated by Turkyilmazoglu [42].The right part moved at a uniform speed, while the left part either moved or remained stationary.Kashyap and Dass [43] examined the tilting angle effect of the 2SLDC with a hot porous block placed at the center of the cavity on the mixed convective flow under seven distinct cavity inclinations.They deduced that the titling angle of the 2SLDC has a substantial impact on thermal transport as well as the generation of entropy.Chowdhury and Rathish Kumar [44] examined the impact of the aspect ratio of the 2SLDC on the convective current of a non-Newtonian-type liquid in a 2SLDC with moving horizontal (top/bottom) boundary walls.The amount of entropy created by any type of irreversible process-including mass and heat transfer, chemical reactions, viscous flow, and Joule heating-is known as entropy generation (EG) [45][46][47][48][49]. Every process in nature yields either zero (0) or positive values for the EG rate.This is a crucial component of the second law of thermodynamics.The analysis of EG in closed domains has been considered in several studies [45][46][47][48][49].
A considerable number of studies in the literature have focused on the convective current of non-Newtonian fluids in a closed domain.However, few studies have dealt with the convective flow of CF in a closed domain while considering the mass (solute) transfer effect.Furthermore, the influence of wall movement on Casson liquid streams has not yet been examined in the existing literature.This concept provides the main motivation for the current study: to numerically survey the interaction between wall movement and the doubly diffusive convection of CF inside a rectangular chamber, while also taking EG into consideration.

Mathematical Model
We consider the time-dependent incompressible laminar Casson liquid flow in a twodimensional wall-driven square enclosed chamber of length L D , as specified in Figure 1.The thermal conditions are constant on the vertical sides and thermal insulation is executed on the horizontal sides; furthermore, constant solute conditions are applied to the vertical sides, and the mass impermeable condition is applied to the horizontal sides.The (vertical) walls of the enclosed chamber move at a constant speed (V o ) in different or matching directions.Three cases were analyzed according to the movement of the wall(s).The direction of gravity's influence is negative in the y-axis.The liquid's thermal properties remain unchanged, with the exception of density (Boussinesq approximation), which is varied as follows: where β s and β T are solutal and thermal expansion coefficients, respectively.It is supposed that the liquid and foreign mass do not react chemically.
the 2SLDC has a substantial impact on thermal transport as well as the generation of entropy.Chowdhury and Rathish Kumar [44] examined the impact of the aspect ratio of the 2SLDC on the convective current of a non-Newtonian-type liquid in a 2SLDC with moving horizontal (top/bottom) boundary walls.
The amount of entropy created by any type of irreversible process-including mass and heat transfer, chemical reactions, viscous flow, and Joule heating-is known as entropy generation (EG) [45][46][47][48][49]. Every process in nature yields either zero (0) or positive values for the EG rate.This is a crucial component of the second law of thermodynamics.The analysis of EG in closed domains has been considered in several studies [45][46][47][48][49].
A considerable number of studies in the literature have focused on the convective current of non-Newtonian fluids in a closed domain.However, few studies have dealt with the convective flow of CF in a closed domain while considering the mass (solute) transfer effect.Furthermore, the influence of wall movement on Casson liquid streams has not yet been examined in the existing literature.This concept provides the main motivation for the current study: to numerically survey the interaction between wall movement and the doubly diffusive convection of CF inside a rectangular chamber, while also taking EG into consideration.

Mathematical Model
We consider the time-dependent incompressible laminar Casson liquid flow in a twodimensional wall-driven square enclosed chamber of length LD, as specified in Figure 1.The thermal conditions are constant on the vertical sides and thermal insulation is executed on the horizontal sides; furthermore, constant solute conditions are applied to the vertical sides, and the mass impermeable condition is applied to the horizontal sides.The (vertical) walls of the enclosed chamber move at a constant speed (Vo) in different or matching directions.Three cases were analyzed according to the movement of the wall(s).The direction of gravity's influence is negative in the y-axis.The liquid's thermal properties remain unchanged, with the exception of density (Boussinesq approximation), which is varied as follows: where βs and βT are solutal and thermal expansion coefficients, respectively.It is supposed that the liquid and foreign mass do not react chemically.The CF rheological model is presented below [21,22]: The CF rheological model is presented below [21,22]: where π = e ij e ij and e ij are the (i, j)th elements for the rate of deformation, π c denotes a critical rate for product of deformation rate, the plastic's dynamic viscosity is denoted by µ pd , and the liquid yield stress is indicated by P ys .The Boussinesq approximation holds for the CF.Based on the aforementioned hypotheses, the main governing model equations are ∂u ∂x ∂u ∂t where β is the CF parameter.The relevant boundary conditions and initial values for the present problem are as follows: The following variables are used to perform non-dimensionalization of the governing model equations: , t * = tV 0 L D , and P = p The dimensionless modeled equations are as follows: ∂U ∂X where N = Gr s Gr T is the buoyancy ratio, is the solutal Grashof number, Re = V 0 L D ν is the Reynolds number, Pr = ν α (= 10) is the Prandtl number, Ri = Gr T Re 2 is the Richardson number, and the Schmidt number is Sc = ν D .The initial and boundary values, expressed in dimensionless form, are as follows: Entropy 2024, 26, 245 5 of 26 The rate of heat transport across the chamber is a crucial constraint in thermal science applications.At the left (hot) barrier, the local Nusselt number is deduced as Nu = − ∂θ ∂X X =0 . The averaged Nusselt number, which may be stated as follows, is used to calculate the total heat transport across the chamber:

NudY
Additionally, the local Sherwood number alongside the cold (right) barrier of the domain is derived as Sh = − ∂C ∂X X=1 . The whole solute transfer across the chamber is represented by the averaged Sherwood number, which is expressed as follows:

Analysis of Entropy Generation
In thermal and industrial applications, buoyancy-induced convection in an enclosed chamber is of great importance.Nonetheless, the concept of entropy generation allows one to identify ideal parameters for real-world use.Entropy production may be computed from the previously known temperature, velocity, and solute concentration fields, as it is associated with the irreversible nature of heat and mass transport as well as viscosity effects.These methods can be used to calculate the (local) entropy generation per unit area of an incompressible fluid.The generation of local entropy is a result of liquid friction, heat, and mass transfer.
The three quantities mentioned above are added to obtain the total (local) entropy production: Applying (10) in a conventional manner, the non-dimensional formula for entropy production is achieved.
Integrating the (dimensionless) local entropy throughout the system, one can obtain the total entropy production: The ratio of total irreversibility to heat transmission irreversibility is described by the local Bejan number: When Be loc > 1 2 , the irreversibility of heat transmission predominates; when Be loc < 1 2 , the irreversibility of fluid resistance predominates.Heat transport and viscous irreversibilities are equivalent at Be loc = 1 2 .The average Bejan number is employed in this study to establish the relative significance of the irreversibility of heat transmission across the domain of interest. .
The ForTran-90 computer software, which was developed for temperature and stream fields, was expanded to compute the formation of entropy within the chamber.

Numerical Technique and Validation
As the FVM is a very famous technique for solving fluid flow problems, it was used to solve the governing model (non-dimensional) Equations ( 10)-( 14), subject to the boundary constraints (15).Convective and diffusive terms were discretized using the QUICK and central difference techniques, respectively.A set number of non-overlapping control volumes were used to sub-divide the computational domain of interest.A test for grid independency was executed with different grid sizes from 42 × 42 to 162 × 162.It was found that a grid of 122 × 122 provided a grid-independent solution.The tri-diagonal matrix method was used to carry out the iterative procedure.The iteration's convergence condition, which is 10 −6 , was maintained.Validation of the present computational code was conducted against existing outcomes for natural convection in a closed cavity [50], as detailed in Table 1.Consistency across the outcomes ensured that the current code for studying the topic is numerically accurate.All numerical simulations were performed using ForTran 90 in-house code.

Results and Discussion
The purpose of the conducted numerical study was to determine how moving-wall orientation, heat, and mass transfer affect the doubly diffusive convection current in a square box occupied with CF.The CF parameter, moving wall direction, solutal Grashof number (Gr s = 10 4 &10 6 ), thermal Grashof number (Gr T = 10 4 &10 6 ), Reynolds number (Re = 10 to 10 3 ), Richardson number (0.01 ≤ Ri ≤ 100), and CF parameter (0.01 ≤ β ≤ 1) were the governing parameters in the experiment.To examine the different effects, simulations were performed under a variety of combinations of the aforementioned parameters.The influence of the moving direction of the vertical isothermal walls was explored in mixed convective flow with heat and mass distribution of Casson liquids in a container.Three cases were considered, which were designed according to different directions of movement; namely, in case I, the hot (left-side) wall moves up while the cold wall (right-side) moves down, while the hot wall moves down while the cold wall moves up in case II.Both walls move up in case III, as shown in Figure 1.
Figure 2 shows the streamlines obtained with different values of Ri under the three cases of wall movement with β = 0.01.There is just one clockwise rotating eddy for all values of Ri in case I, where one eddy fills the entire domain as a result of the buoyancy force and the heated barrier moving upward.The eddy's core region elongates horizontally, and two inner eddies form when Ri = 1 and Ri = 100.The opposite trend was observed in case II; that is, a single anti-clockwise rotating eddy was observed inside the cavity.Here, the hot (left) wall moves down and against the buoyancy force near the isothermal wall.The cell stretched horizontally for all Ri numbers in both cases I and II.However, the stream pattern was totally different in case III, where both barriers moved upwards.This resulted in a dual cell structure for all Ri values.When both walls moved in the same direction (i.e., upwards), the shear force drove the liquid near the wall to move up and form a dual cell pattern.The mechanism of heat transmission changed as the Richardson number increased, but the flow structure was not significantly impacted.
Figure 3 displays the flow patterns under the three cases and three modes of heat transfer with β = 1.Although a single cell pattern was observed for all Ri values in case I, the location of the eddy core region depended on the value of Ri.In the forced convection regime (Ri = 0.01), the cell stretched diagonally and the core region of the eddy was observed at the lower left portion of the cavity.When Ri = 1 and 100 (i.e., in mixed and natural convection regimes), the cell stretched horizontally and two inner cells formed in the core area of the eddy.When considering case II, an anti-clockwise rotating single cell formed under forced and mixed convection regimes.Nevertheless, the flow pattern completely changed in the free convection regime, where a tri-cellular pattern was observed.This flow structure was not observed for lower values of β.The dual eddy arrangement formed for all values of Ri in case III.However, the size of the cells depended on the value of R; that is, the mode of heat transfer.In the forced convective flow regime, the two eddies were of the same size and power, and the cell generated by the hot wall became strengthened.Further increasing Ri (Ri = 100), the hot cell gained more strength and occupied the majority of the container.Consequently, the eddy generated by the cold wall shrank in size.Comparing Figures 2 and 3, the stream pattern was affected much more in the case β = 1, compared to that under β = 0.01, when changing the values of Ri.
studying the topic is numerically accurate.All numerical simulations were performed using ForTran 90 in-house code.8.825 8.935

Results and Discussion
The purpose of the conducted numerical study was to determine how moving-wall orientation, heat, and mass transfer affect the doubly diffusive convection current in a square box occupied with CF.The CF parameter, moving wall direction, solutal Grashof number ( were the governing parameters in the experiment.To examine the different effects, simulations were performed under a variety of combinations of the aforementioned parameters.The influence of the moving direction of the vertical isothermal walls was explored in mixed convective flow with heat and mass distribution of Casson liquids in a container.Three cases were considered, which were designed according to different directions of movement; namely, in case I, the hot (left-side) wall moves up while the cold wall (right-side) moves down, while the hot wall moves down while the cold wall moves up in case II.Both walls move up in case III, as shown in Figure 1.
Figure 2 shows the streamlines obtained with different values of Ri under the three cases of wall movement with  = 0.01.There is just one clockwise rotating eddy for all values of Ri in case I, where one eddy fills the entire domain as a result of the buoyancy force and the heated barrier moving upward.The eddy's core region elongates horizontally, and two inner eddies form when Ri = 1 and Ri = 100.The opposite trend was observed in case II; that is, a single anti-clockwise rotating eddy was observed inside the cavity.Here, the hot (left) wall moves down and against the buoyancy force near the isothermal wall.The cell stretched horizontally for all Ri numbers in both cases I and II.However, the stream pattern was totally different in case III, where both barriers moved upwards.This resulted in a dual cell structure for all Ri values.When both walls moved in the same direction (i.e., upwards), the shear force drove the liquid near the wall to move up and form a dual cell pattern.The  Figure 4 demonstrates the thermal distribution inside the container for different Ri and wall movement directions (three cases) with  = 0.01.In case I, thermal boundary layers formed along the moving isothermal walls.Higher thermal gradients were found at the top portion of the container under Ri = 0.01.In forced convection (Ri = 100), the isotherms were distributed equally inside the container, and vertical thermal stratification was observed.However, higher thermal gradients were found at the bottom of the cavity for Ri = 0.01 and 1 in case II.When the walls moved in the same direction (case III), the thermal distribution totally varied from case I and case II.At Ri = 0.01, a high thermal gradient existed at the central part of the container for  = 0.01 and 1.When Ri = 100, the isotherms were distributed over almost the whole container.The convection mode of heat transmission was predominant in every case examined here.Figure 5   Figure 4 demonstrates the thermal distribution inside the container for different Ri and wall movement directions (three cases) with β = 0.01.In case I, thermal boundary layers formed along the moving isothermal walls.Higher thermal gradients were found at the top portion of the container under Ri = 0.01.In forced convection (Ri = 100), the isotherms were distributed equally inside the container, and vertical thermal stratification was observed.However, higher thermal gradients were found at the bottom of the cavity for Ri = 0.01 and 1 in case II.When the walls moved in the same direction (case III), the thermal distribution totally varied from case I and case II.At Ri = 0.01, a high thermal gradient existed at the central part of the container for Ri = 0.01 and1.When Ri = 100, the isotherms were distributed over almost the whole container.The convection mode of heat transmission was predominant in every case examined here.Figure 5 displays the thermal pattern inside the container for the three cases of moving walls and under different values of Ri with β = 1.Every figure demonstrates the convection mechanism of heat transmission in detail.Strong thermal boundary layers developed alongside the isothermal walls at Ri = 0.01 and 1 in cases I and II.In the free convection regime, the isotherms dispersed effectively within the container.No thermal boundary layer existed in case III.However, there were significant temperature variations in the cavity's center (X = 0.5).When increasing Ri = 1, the thermal gradient moved slightly towards the right side.At Ri = 100, these thermal gradients appeared close to the cold wall.at  = 0.01 and 1 in cases I and II.In the free convection regime, the isotherms dispersed effectively within the container.No thermal boundary layer existed in case III.However, there were significant temperature variations in the cavity's center ( = 0.5).When increasing  = 1, the thermal gradient moved slightly towards the right side.At Ri = 100, these thermal gradients appeared close to the cold wall.Figure 6 demonstrates the solute distributions inside the container for diverse values of Ri and the three cases of moving wall directions with  = 0.01.The solutal mixing inside the container was strong in case I, and solutal gradient layers formed along the moving vertical barriers at  = 0.01 and 1 for cases I and II.We observed a similar trend of solutal mixing while both walls were moving oppositely (Cases I and II).As the walls in cases I and II were moving in the same direction, the solute mixing in case III completely differed from that in cases I and II.Here, a very strong solutal gradient existed in the middle of the container, due to the dual cell stream pattern due to moving the wall in the same direction.The isomasslines were distributed equally inside the container at Ri = 100 for all cases of the moving walls. Figure 7 demonstrates the solutal mixing inside the container for various values of Ri and the three cases of moving walls with  = 1.It can be seen, from these graphs, that the mass transmission is clearly influenced by the wall's movement and its direction.The liquid was well-mixed inside the container in all cases, Figure 6 demonstrates the solute distributions inside the container for diverse values of Ri and the three cases of moving wall directions with β = 0.01.The solutal mixing inside the container was strong in case I, and solutal gradient layers formed along the moving vertical barriers at Ri = 0.01and 1 for cases I and II.We observed a similar trend of solutal mixing while both walls were moving oppositely (Cases I and II).As the walls in cases I and II were moving in the same direction, the solute mixing in case III completely differed from that in cases I and II.Here, a very strong solutal gradient existed in the middle of the container, due to the dual cell stream pattern due to moving the wall in the same direction.The isomasslines were distributed equally inside the container at Ri = 100 for all cases of the moving walls. Figure 7 demonstrates the solutal mixing inside the container for various values of Ri and the three cases of moving walls with β = 1.It can be seen, from these graphs, that the mass transmission is clearly influenced by the wall's movement and its direction.The liquid was well-mixed inside the container in all cases, boosting the strong solutal mixing.The mass distribution showed a similar trend as discussed for Figure 6 in the considered cases.Comparing Figures 6 and 7, the influence of the CF parameter on mass transfer can be clearly understood.The kinetic energy (KE) obtained under different values of the Casson fluid parameter, Ri, and the three cases are displayed in Figure 8.The KE behaved non-linearly with increasing values of Ri, and declined in most of the considered cases when increasing the Casson fluid parameter.Case III provided a lower KE than the other two cases in most of the simulations.Figure 9 shows the skin friction (SF) along the moving right wall under numerous values of , Ri, and the three cases of the moving wall directions with N = 0.5.The skin friction increased when increasing the value of  at  ≤ 1, was almost constant when Ri = 10, and declined when Ri = 100 in case I.As the hot (left) wall moves up in cases I and III, a similar trend of variation in SF was observed when changing the Ri and  values.The SF declined when increasing the Ri values for all  in cases I and III.The deviation of the SF value between the two Ris increased when intensifying the CF parameter; that is, the SF increased (decreased) when increasing the value of the CF parameter under forced and mixed (free) convection regimes.When considering case II, SF declined The kinetic energy (KE) obtained under different values of the Casson fluid parameter, Ri, and the three cases are displayed in Figure 8.The KE behaved non-linearly with increasing values of Ri, and declined in most of the considered cases when increasing the Casson fluid parameter.Case III provided a lower KE than the other two cases in most of the simulations.Figure 9 shows the skin friction (SF) along the moving right wall under numerous values of β, Ri, and the three cases of the moving wall directions with N = 0.5.The skin friction increased when increasing the value of β at Ri ≤ 1, was almost constant when Ri = 10, and declined when Ri = 100 in case I.As the hot (left) wall moves up in cases I and III, a similar trend of variation in SF was observed when changing the Ri and β values.The SF declined when increasing the Ri values for all β in cases I and III.The deviation of the SF value between the two Ris increased when intensifying the CF parameter; that is, the SF increased (decreased) when increasing the value of the CF parameter under forced and mixed (free) convection regimes.When considering case II, SF declined when increasing the β value; in particular, the skin friction declined with increasing β for all Ri except Ri = 1 in case II.When Ri = 1, SF behaved constantly for all values of β (Case II).Comparing the three cases, cases I and III presented higher drag along the moving wall.As the wall moved downwards and opposite to the buoyancy force in Case II, the skin friction value took a negative sign, and was lower than in the other two cases.
To illustrate the impact of the averaged heat transmission across the container, the averaged Nusselt number is displayed, against different CF parameter and Ri values, for the three cases in Figure 10.The averaged heat transfer rate remained constant with a change in the Casson fluid parameter when Ri ≥ 10 in case I.When Ri ≤ 1, the averaged heat transfer rate declined slightly first, then became constant in case I.The averaged Nusselt number declined when increasing the CF parameter in case II for all values of Ri; in particular, the averaged heat transfer experienced a substantial reduction when β changed from 0.01 to 0.1.In case III, a similar decline also occurred when β increased from 0.01 to 0.1.From the three cases, it was found that the averaged Nusselt number deteriorated when increasing the value of Ri; that is, the mode of forced convection (Ri << 1) provides higher heat transport than the free and mixed convection regimes (Ri ≥ 1).Comparing the three cases, the averaged Nusselt number obtained lower values in case III; that is, the averaged heat transport was low in this case.This is due to the movement of an isothermal (left-side) wall against the buoyancy force.However, the left hot wall moves up, supporting the buoyancy force to produce a higher heat transmission rate in cases I and III.The natural convection mode was not influenced significantly by the heat transmission rate when varying the CF parameter values in all three cases.
when increasing the  value; in particular, the skin friction declined with increasing  for all Ri except Ri = 1 in case II.When Ri = 1, SF behaved constantly for all values of  (Case II).Comparing the three cases, cases I and III presented higher drag along the moving wall.As the wall moved downwards and opposite to the buoyancy force in Case II, the skin friction value took a negative sign, and was lower than in the other two cases.paring the three cases, the averaged Nusselt number obtained lower values in case III; that is, the averaged heat transport was low in this case.This is due to the movement of an isothermal (left-side) wall against the buoyancy force.However, the left hot wall moves up, supporting the buoyancy force to produce a higher heat transmission rate in cases I and III.The natural convection mode was not influenced significantly by the heat transmission rate when varying the CF parameter values in all three cases.Figure 11 displays the averaged thermal transfer rate across the container against the CF parameter and buoyancy ratio parameter (N) for the three cases.It can be observed, from Figure 11, that case I offered the higher heat transport rate among the three cases, whereas case III provided a lower heat transport rate.When raising the buoyancy ratio, the averaged Nusselt number declined in cases I and III, while the opposite trend was Figure 11 displays the averaged thermal transfer rate across the container against the CF parameter and buoyancy ratio parameter (N) for the three cases.It can be observed, from Figure 11, that case I offered the higher heat transport rate among the three cases, whereas case III provided a lower heat transport rate.When raising the buoyancy ratio, the averaged Nusselt number declined in cases I and III, while the opposite trend was observed in case II.When changing the Casson liquid parameter from 0.01 to 0.1, a sharp decrease in the averaged Nusselt number was found in Cases I and II, after which the averaged Nusselt number declined gradually with increasing β.The average heat transmission decreased as the CF parameter increased, as this results in higher friction between the liquid particles.
To explore the mass transport rate, the averaged Sherwood number is illustrated against the Casson fluid parameter and Ri in Figure 12 for the three cases of wall movement.The figure clearly shows that the mass transport rate deteriorated when increasing the Ri value in all three cases; that is, the forced convective regime provided a higher mass transport rate in all cases.Comparing the three cases, Case II offered a higher mass transport rate, whereas case III provided a lower mass transfer rate under all modes of heat mass transfer.At Ri = 100, the mass transfer rate remained constant when changing the value of the Casson parameter.It was also observed that the averaged Sherwood number decreased with increasing β.As a dual cell structure is created due to the same-direction movement of both walls, the fluid was not well-mixed, resulting in a deterioration of mass transfer in case III. Figure 13 shows the effect of the averaged Sherwood number for several values of the CF parameter and buoyancy ratio with Ri = 1.A sudden decline in the averaged Sherwood number was observed when altering the value of β from 0.01 to 0.1, after which it dropped gradually with an increase in the value of the CF parameter.The mass transport rate deteriorated as the buoyancy ratio (N) value increased in case I.However, cases II and III showed the opposite trend; for example, the mass transfer rate was enhanced under case I, whereas a lower rate of transfer was observed under case III.This observation is due to the direction(s) of the moving wall(s).
Figure 14 displays the Bejan number for numerous values of β and Ri, under the three cases of wall movement.The Bejan number was higher than 0.85 (i.e., Be > ½), indicating that convection heat transmission was dominant in the system.When β = 1, the Bejan number tended to unity ( Be → 1 ) for all values of Ri in the three cases.It is evident that heat transport irreversibility is dominant.For lower values of β (β = 0.01, 0.1), the Bejan number increased when decreasing the value of Ri.In all the cases of wall movement, a uniform trend of the Bejan number was observed when increasing the value of β. Figure 15 shows the averaged temperature for several values of Ri, β, and the three cases of wall movement.It can be seen, from these figures, that T Avg increased strongly with increasing Ri and β values for all three cases.The averaged temperature was high (minimum) in the free (forced) convective regime for all three cases.Figure 16 demonstrates the (RMSD-T Avg ) RMSD averaged temperature inside the container for numerous values of β, Ri, and the three cases of wall movement.The value of RMSD-T Avg increased when increasing the value of Ri in cases I and II; however, it behaved oppositely to Ri in case III.The value of RMSD-T Avg was always below ½.  was enhanced under case I, whereas a lower rate of transfer was observed under case III.This observation is due to the direction(s) of the moving wall(s).was enhanced under case I, whereas a lower rate of transfer was observed under case III.This observation is due to the direction(s) of the moving wall(s).increasing Ri and  values for all three cases.The averaged temperature was high (minimum) in the free (forced) convective regime for all three cases.Figure 16 demonstrates the (RMSD-TAvg) RMSD averaged temperature inside the container for numerous values of , Ri, and the three cases of wall movement.The value of RMSD-TAvg increased when increasing the value of Ri in cases I and II; however, it behaved oppositely to Ri in case III.The value of RMSD-TAvg was always below ½.

Conclusions
The impacts of the movement direction of the wall(s) and entropy generation on doubly diffusive mixed convective stream, as well as the heat and mass transfer of a Casson liquid in a closed container, were numerically explored in this study, for which three cases

Conclusions
The impacts of the movement direction of the wall(s) and entropy generation on doubly diffusive mixed convective stream, as well as the heat and mass transfer of a Casson liquid in a closed container, were numerically explored in this study, for which three cases of wall movement were considered.The following main observations were detected in this study: ➢ Single-or multi-cellular flows were observed, according to the different directions of the moving wall(s).A single stream pattern was seen in case I, whereas a dual cellular stream existed in case III for all values of the considered parameters.➢ Kinetic energy varied non-linearly with the CF parameter under all convective regimes.➢ Comparing the three cases, the averaged heat transport was low in case III due to the movement of an isothermal (left) wall against the buoyant force.➢ The resistance of various liquid layers to their relative motion increased with an increase in the CF parameter, resulting in a decline in the averaged heat transport when increasing the CF parameter.➢ While strengthening the buoyancy ratio did not support the mass transfer rate in case I, cases II and III presented the opposite trend.➢ The mass transfer rate was enhanced in case I and diminished in case III due to the direction of the moving walls.➢ The averaged heat transfer was high under the forced convection mode in all three cases, and was minimal under the free convection mode in all three cases.➢ The Bejan number increased when increasing the value of the Casson fluid parameter.

Figure 1 .
Figure 1.Physical configuration with three cases.

Figure 1 .
Figure 1.Physical configuration with three cases.
Figure4demonstrates the thermal distribution inside the container for different Ri and wall movement directions (three cases) with  = 0.01.In case I, thermal boundary layers formed along the moving isothermal walls.Higher thermal gradients were found at the top portion of the container under Ri = 0.01.In forced convection (Ri = 100), the isotherms were distributed equally inside the container, and vertical thermal stratification was observed.However, higher thermal gradients were found at the bottom of the cavity for Ri = 0.01 and 1 in case II.When the walls moved in the same direction (case III), the thermal distribution totally varied from case I and case II.At Ri = 0.01, a high thermal gradient existed at the central part of the container for  = 0.01 and 1.When Ri = 100, the isotherms were distributed over almost the whole container.The convection mode of heat transmission was predominant in every case examined here.Figure5displays the thermal pattern inside the container for the three cases of moving walls and under different values of Ri with  = 1.Every figure demonstrates the convection mechanism of heat transmission in detail.Strong thermal boundary layers developed alongside the isothermal walls

Figure 3 .
Figure 3. Streamlines for different Ri values under the three cases with β = 1.

Figure 5 .
Figure 5. Isotherms for different Ri values under the three cases with β = 1.

9 Figure 7 .
Figure 7. Isomasslines for different Ri values under the three cases with β = 1.

Figure 8 .
Figure 8. Kinetic energy for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).Figure 8. Kinetic energy for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 8 .
Figure 8. Kinetic energy for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).Figure 8. Kinetic energy for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 9 .
Figure 9. Skin friction for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).Figure 9. Skin friction for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 9 .
Figure 9. Skin friction for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).Figure 9. Skin friction for different values of the Casson fluid parameter and Ri in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 10 .
Figure 10.Averaged Nusselt number for different values of Ri and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 10 .
Figure 10.Averaged Nusselt number for different values of Ri and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 11 .
Figure 11.Averaged Nusselt number for different values of N and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).Figure 11.Averaged Nusselt number for different values of N and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 11 .
Figure 11.Averaged Nusselt number for different values of N and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).Figure 11.Averaged Nusselt number for different values of N and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 12 .
Figure 12.Averaged Sherwood number for different values of Ri and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 12 .Figure 13 .
Figure 12.Averaged Sherwood number for different values of Ri and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Figure 14 .
Figure 14.Bejan number for different values of Ri and Casson fluid parameters in Case 1 (a), Case 2 (b), and Case 3 (c).

Author Contributions:
Conceptualization, S.S.; methodology, S.S.; software, S.S.; validation, S.S.; investigation, S.S.; data curation, S.S. and M.B.; writing-original draft preparation, S.S.; writingreview and editing, S.S., M.B. and A.K.A.; visualization, S.S. and M.B.; project administration, S.S.; funding acquisition, S.S.All authors have read and agreed to the published version of the manuscript.Funding: This work was supported by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia, under Grant No. (D-216-130-1439).The authors, therefore, gratefully acknowledge the DSR technical and financial support.Institutional Review Board Statement: Not applicable.Data Availability Statement: Data are unavailable for privacy reasons.