Numerical Investigation of the Electro-Thermo Convection in an Inclined Cavity Filled with a Dielectric Fluid

: The present work is a numerical analysis of electro-thermo convection, occurring in a square differentially heated cavity ﬁlled with a dielectric ﬂuid. The cavity experiences the combined effects of viscous, electrical, and thermal forces. The equations modelling the physical problem are solved via the ﬁnite volume approach. The study focuses on the effect of cavity tilt on the ﬂuid ﬂow structure and thermal performance inside the enclosure under the action of an electric ﬁeld. A parametric study was performed, where the tilt angle is getting varied between 0 ◦ and 90 ◦ , as well as the Rayleigh number (5000 ≤ Ra ≤ 250,000) and the electric ﬁeld (0 ≤ T ≤ 800). Furthermore, the electric charge injection level C, the mobility M and the Prandtl Pr numbers were all adjusted to a value of 10. The obtained results demonstrate that the hydrodynamic and thermal ﬁelds are signiﬁcantly impacted by the cavity inclination. In addition, regardless of the thermal Rayleigh’s number, high electric ﬁeld values could govern ﬂuid movement through electric forces. Electro-convection typically demonstrates an oscillating ﬂow due to the tilting of the cavity which gives rise to a bicellular regime occupying the entire cavity. A correlation has been established to estimate heat transfer by considering various system parameters such as cavity inclination, electrical Rayleigh number, and thermal Rayleigh number.


Introduction
Electric field combined to heat transfer has garnered growing consideration in the past few years due to its wide scope of benefits in many engineering fields, such as the formation of crystals in a liquid, nuclear reactors cooling, electronic equipment thermal dissipation, as well as microelectronic devices, solar technology, fuel cells [1], etc.During the last decades, several leading journals covered various issues of electrohydrodynamic improvement of heat transfer [2][3][4] within square, cylindrical, and elliptical cavities.An intensive study of EHD convection was given by Hassen et al. [2].
However, cavity inclination can also have a significant influence on heat transfer encountered in both industrial and natural fields.Existing literature includes numerous works relating to the subject, without the electric field and/or, in the presence of a magnetic field.
Numerical studies carried out by Mourabit et al. [5] examined the effect of 'T' shape cavity inclination equipped with heating blocks.The flow behavior and heat transfer rate in the case of mixed convection have been studied.The authors found that tilt promotes cooling of the left block and disadvantages cooling of the right block.
Cheong et al. [6], have examined free convection in a tilted rectangular cavity with a sine-wave temperature applied on the left wall.Streamlines, isotherms, and Nusselt numbers are illustrated and studied.It was noticed that heat exchange decreases with increasing cavity tilt for all aspect ratios of the simulated geometries.The authors established a correlation to quantify the mean Nusselt number by considering the thermal Rayleigh number, cavity tilts angle, and form factors as variables.
Yekani et al. [7] have also numerically analyzed the free convection of a nanofluid in a porous enclosure with differentially heated vertical surfaces, the numerical simulations have been carried out for various porous Grashof numbers and for different inclinations of the enclosure (from 0 • to 60 • ).The numerical findings reveal the existence of a boundary layer adjacent to the walls of the enclosure, and its thickness diminishes as the porosity decreases or the porous Grashof number increases.
Ouyahia et al. [8] investigate numerically the impact of tilt on magnetohydrodynamic, laminar-free convection within a rectangular cavity containing an Al 2 O 3 -water nanofluid.This cavity was subjected to a horizontal, constant, and uniform magnetic field.The vertical walls were differentially heated, and the inclination was varied between 0 • and 60 • .It was found that the variation of the tilt angle leads to changes in the structure of the flow.Indeed, the best thermal performances were obtained for an angle of inclination of 30 • , unlike 60 • which has the lowest Nusselt number values.
Benyoucef et al. [9] have investigated the 3D convection of a tilted rectangular cavity filled with silicone oil and compared it with experimental results.A complex behavior has been detected, the increase in the angle beyond 0 • destabilizes the convective pattern to benefit from an unstable situation for specific Rayleigh number values.The authors demonstrated that the mean Nusselt number achieves more rapid convergence to the stationary state when the angles fall within the range of 25 • to 60 • .
Hakan et al. [10] used the heatline visualization technique to visualize the heat transport in a non-isothermally heated square inclined (from 0 • to 90 • ) cavity filled with water-CuO nanofluid.The authors showed that for low Ra around 10 3 , increasing the angle of inclination leads to an improvement in the Nusselt number.For Ra = 10 4 and 10 5 , the average Nusselt number decreases as the angle of inclination are varied from 30 • to 60 • and increases again between 60 • and 90 • .
Vivek et al. [11] conducted a numerical study to analyze the effects of the interaction between surface radiation and natural convection in air-filled tilted enclosures.The study covered a wide range of tilt angles (α).The authors have shown that the presence of surface radiation results in an enhancement of the overall convective circulation rate for specific tilt angles.The Convective Nusselt number reaches its maximum value at an angle of α = +15 • and remains unaffected by the presence or absence of radiation.The significance of the interaction effects between radiation and convection is much higher in shallow enclosures compared to square enclosures.In a shallow cavity with an aspect ratio of 5, a positive inclination of 45 • reduces convection by 50%, while negative inclinations increase convection by more than 100%.
Very recently, Mandal et al. [12] studied the effect of a magnetic field on free convection in an inclined porous cavity with corrugated walls.The cavity was filled with Cu-Al 2 O 3 /water hybrid nanofluid and partially heated from the top and bottom walls.It has been found that heat transfer and air movement increase with the inclination of the cavity.The direction of the magnetic field and the volume fraction of the added nanoparticles are also parameters that control heat transfer.
Several quite recent experimental studies have also focused on natural convection in inclined cavities.Leporini et al. [13] experimentally studied the behavior of natural convection in an inclined square cavity with hot side walls and a cold bottom wall.Four different cavity inclination angles, from 0 to 90 • , and five different hot wall temperatures (10 4 ≤ Ra ≤ 10 5 ) were tested for a total of 20 different conditions.The techniques adopted for the experimental tests are particle image velocimetry (PIV) to determine velocity fields and a set of thermocouples to measure wall temperature.It has been determined that both the average and maximum velocities exhibit an increase as the Rayleigh number and tilting angle increase.The velocity field exhibits two vortices, with a symmetry axis crossing, the center of the enclosure when the tilt angle is 90 • ; as this angle increases, the symmetry tends to disappear, resulting in the presence of only one vortex for 90 • .
Abdulmunem et al. [14] published an experimental study to investigate the effects of different tilt angles on the melting of a PCM material (paraffin wax) inside a rectangular enclosure.Four different angles ranging from 0 to 90 • were tested by subjecting a single wall to a constant temperature of 373 Kelvin.Heating was carried out using a 1000 W sun simulator to achieve a uniform heat distribution across the surface of the plate.A total of 10 K-type thermocouples were strategically positioned within the paraffin wax to facilitate accurate measurements of temperature.The authors show that the overall melting time increases by more than 75% when the cavity tilt angle is reduced from 90 to 0 • .It was also noted that decreasing the angle of inclination reduces the effect of convective heat transfer during MCP melting until a purely conductive heat transfer state is reached at 0 • .This literature review unequivocally illustrates the impact of cavity tilt which is important to control heat transfer as well as flow patterning.However, as far as we are aware, no previous study has examined this effect in the presence of an electric field.Therefore, the main goal of this paper is to assess the hydrodynamic and thermal performance of a dielectric fluid enclosed within an inclined cavity under the action of an electric field.The influence of the tilt angle, the electric and thermal Rayleigh numbers, on heat transfer and flow topography was illustrated and discussed.

Physical Model and Governing Equations
The geometry in question is a layer of insulating (dielectric) liquid located in a square cavity of dimension L 2 (Figure 1).The bottom and top surfaces are thermally insulated, while the side surfaces are differentially heated and held at a constant and homogenous temperature (with θ H denoting the hot side and θ C the cold one).To achieve an enhanced heat exchange, a potential difference is created between the horizontal surfaces.The bottom wall serves as the emitting electrode and is maintained at the highest potential (V i ).The top wall is maintained at the lowest potential (V 0 ) which is considered as the receiving electrode.
Mass equation Momentum equation temperature (with θH denoting the hot side and θC the cold one).To heat exchange, a potential difference is created between the horizo tom wall serves as the emi ing electrode and is maintained at the The top wall is maintained at the lowest potential (V0) which is consi electrode.Thermal energy equation Maxwell's conservation equations ∂q ∂t + ∇.
→ j = 0 where ∇.(ε As a numerical approach, the "Vorticity-Current function" [30] (ψ-ω) formalism was adopted, using the following variable changes: To describe problems universally, it is especially useful to work with dimensionless equations.For this purpose, the following transformations are conducted, To describe problems universally, it is especially useful to work with dimensionless equations.For this purpose, the following transformations are conducted, x y ∂q ∂t In the above-mentioned equations, five dimensionless numbers are defined as follows: .
The thermal characteristic of the system is represented by the Rayleigh number Ra and the Prandtl number Pr. Ra denotes the quotient of thermal force and the viscous resistance force, whereas Pr is the momentum to thermal diffusivities ratio.The remaining three numbers are related to the electrical side.The electric Rayleigh number T describes the electric Coulomb force to viscous force.The injection level is denoted by C, and the dimensionless ionic mobility is represented by M.
The boundary and initial conditions are expressed in the following manner: For t = 0 For t > 0 To quantify the heat transfer level, we employ the following equations to calculate both the local and average Nusselt numbers:

Numerical Method
The partial differential equations (PDEs) ( 7) to ( 12) and the related boundary conditions are resolved by performing numerical integration over a control volume.The power law scheme, described in [31], has been employed to handle the convective terms.The time derivative terms are discretized using a semi-implicit Euler first-order scheme, as described in reference [30].The semi-implicit first order Euler time scheme was chosen to discretize the temporal derivatives terms.This choice was made in order to improve the precision and reliability of the solution.The scheme relies on a hybrid approach comprising explicit and implicit components, where the explicit part is used to approximate the nonlinear terms and the implicit part is used to approximate the linear terms.This combination ensures that the numerical solution is stable and accurate, and also allows for a fast convergence of the numerical solution.
After discretization, a system of algebraic equations is generated for variables such as temperature, stream function, vorticity, electric potential, and charge density.To solve this system, the Successive Over Relaxation technique (SOR) is employed, utilizing the Gauss-Seidel iterative solver [31].To guarantee the reliability of the results, grid independency checks were conducted to determine the optimal mesh size for the simulation.After a series of tests, it was found that a 101 × 101 control volume mesh was sufficient to obtain the required precision.The dimensionless time step must be chosen carefully to guarantee the stability of the calculations, for this reason a small enough time step of 10 −4 is found adequate to serving this purpose while stills being computationally efficient [32].The numerical solution can be considered as converged and the SOR algorithm can be stopped when the following convergence criterion is being verified.
where the exponent k denotes the kth step of the successive overrelaxation method.
For further details about the calculation methodology, please consult the paper by Hassen et al. [33].

Results and Discussion
In this section, the combined effect of the cavity tilt (from 0 • to 90 • ) coupled to an electric field in the presence of a temperature gradient is investigated.This refers to the rotation of the cavity in a clockwise direction, which enables the study of various geometries ranging from the position where the thermal gradients are perpendicular to gravity (α = 0) to the Rayleigh-Bénard situation (α = 90) where gravity counteracts thermal gradients.Liquid gas oil is selected as the working medium.The thermophysical properties of this fluid, detailed in [34], enable the calculation of both the Prandtl number and the ionic mobility, with values set to Pr = M = 10.The charge injection level is classified into three regimes [35]: strong (C > 5), medium (0.2 < C < 5), and weak (0 < C < 0.2).Among these regimes, previous studies have mainly opted for strong injection with C = 10, which seems to be the most widely used choice.The thermal and electrical Rayleigh numbers as well as the inclination angle were the subject of the parametric study.As shown in Table 1, the parametric study begins by investigating the impact of inclination (0 ≤ α ≤ 90) while keeping Ra and T constant, and then explores the effect of inclination under various applied electric fields (0 ≤ T ≤ 800) and different imposed thermal gradients (5000 ≤ Ra ≤ 250,000).Figure 2 presents the steady-state charge density distribution, streamlines, and isothermal lines plots for Ra = 50,000 and T = 600.In the case of 0 • tilt, it is clearly seen that the collaboration between the thermal and the electric field on the right side of the cavity generates a sudden acceleration of the cell rotational speed.The flow structure is characterized by a circular vortex with a center shifted to the hot wall.The distortions of the isothermal lines are so pronounced in the center of the cavity that a thermal gradient inversion takes place.On the other hand, on the upper right and lower left corners, a strong intensification of this gradient is observed.As for the electric charge distribution, we can also identify the inversion of the gradient from the center to more accentuated near to the isothermal walls.
Tilting the cavity up to 30 • results in a slight 10% increase in vortex velocity, but no significant alteration in the flow topography is observed.Above an inclination of 45 • , a bifurcation from a single-cell flow to a two-cell flow is noted.
This bifurcation can be explained by the fact that the cavity inclination causes a reduction in the intensity of buoyancy forces (going from ρ while the electrical forces remain unaffected.However, the left portion of the cavity is the location of collaboration between the electrical and thermal forces, generating a large and intense cell and a strong penetration of electrical charge.For isothermal lines, this also leads to a division of the cavity into two compartments.A hot and dominating left compartment and a cold one on the right side.
The same configuration persisted until 60 • ; however, when an angle of 90 • was reached, two horizontally oriented and symmetric counter-rotating cells were observed.These cells are located at the top and bottom regions of the cavity.
Figure 3 presents both the temporal evolution of the mean and local Nusselt number on the hot wall.Two different behaviors for heat transfer can be distinguished.For an angle α ∈ [0, 40 • ], an increase in the angle of tilt results in a 5% improvement in heat transfer.From 45 • , called the critical angle which separates the one-cell regime from the two-cell regime, an oscillatory regime appears, leading to a significant decrease (around 33%) in the Nusselt number.By further increasing the tilt angle, the amplitude, and frequency of the oscillations increase, leading to a subsequent increase in heat transfer.
Finally, it is essential to point out that the bottom part of the hot wall experiences the most intense heat exchange, particularly for 'y' values ranging from 0.25 to 0.3.Conversely, the top section of the heated wall exhibits inefficient heat transfer, regardless of the cavity's angle.Finally, it is essential to point out that the bo om part of the hot wall experiences the most intense heat exchange, particularly for 'y' values ranging from 0.25 to 0.3.Conversely, the top section of the heated wall exhibits inefficient heat transfer, regardless of the cavity's angle.Charge density distributions for various tilt angles and electrical Rayleigh numbers are illustrated in Figure 4.At an extremely low electrical Rayleigh number (T = 10), thermal forces dominate, and buoyancy forces entirely dictate the movement of electric charges.Indeed, electric charges exhibit an upward orientation in the vicinity of the hot wall, while they manifest a downward direction close to the cold wall.This creates a kind of ring inside the cavity.The inclination of the cavity seems to have no impact on the charge density distribution.Conversely, for a large electric Rayleigh number (T = 600), electrical forces take control.Both electric plumes along the vertical walls take place regardless of the cavity's tilt.Finally, for a moderate electrical Rayleigh number.A single plume originates along the hot wall, where the thermal and electrical forces cooperate.However, on the cold wall side, the downward thermal forces seem to be more dominant for α = 0 and 45°, preventing the development of a second plume.In this case (moderate T value), the cavity inclination affects the charge density distribution topography.As the angle α reaches 90° a shift occurs from a single-plume pa ern to a bifurcated-plume arrangement situated at both the upper and lower regions of the cavity.For α = 0°, the charge distribution aligns closely with the findings reported by Hassen et al. [36] and Traore et al.Charge density distributions for various tilt angles and electrical Rayleigh numbers are illustrated in Figure 4.At an extremely low electrical Rayleigh number (T = 10), thermal forces dominate, and buoyancy forces entirely dictate the movement of electric charges.Indeed, electric charges exhibit an upward orientation in the vicinity of the hot wall, while they manifest a downward direction close to the cold wall.This creates a kind of ring inside the cavity.The inclination of the cavity seems to have no impact on the charge density distribution.Conversely, for a large electric Rayleigh number (T = 600), electrical forces take control.Both electric plumes along the vertical walls take place regardless of the cavity's tilt.Finally, for a moderate electrical Rayleigh number.A single plume originates along the hot wall, where the thermal and electrical forces cooperate.However, on the cold wall side, the downward thermal forces seem to be more dominant for α = 0 and 45 • , preventing the development of a second plume.In this case (moderate T value), the cavity inclination affects the charge density distribution topography.As the angle α reaches 90 • a shift occurs from a single-plume pattern to a bifurcated-plume arrangement situated at both the upper and lower regions of the cavity.For α = 0 • , the charge distribution aligns closely with the findings reported by Hassen et al. [36] and Traore et al. [37].
Finally, it is essential to point out that the bottom part of the hot wall experiences the most intense heat exchange, particularly for 'y' values ranging from 0.25 to 0.3.Conversely, the top section of the heated wall exhibits inefficient heat transfer, regardless of the cavity's angle.Charge density distributions for various tilt angles and electrical Rayleigh numbers are illustrated in Figure 4.At an extremely low electrical Rayleigh number (T = 10), thermal forces dominate, and buoyancy forces entirely dictate the movement of electric charges.Indeed, electric charges exhibit an upward orientation in the vicinity of the hot wall, while they manifest a downward direction close to the cold wall.This creates a kind of ring inside the cavity.The inclination of the cavity seems to have no impact on the charge density distribution.Conversely, for a large electric Rayleigh number (T = 600), electrical forces take control.Both electric plumes along the vertical walls take place regardless of the cavity's tilt.Finally, for a moderate electrical Rayleigh number.A single plume originates along the hot wall, where the thermal and electrical forces cooperate.However, on the cold wall side, the downward thermal forces seem to be more dominant for α = 0 and 45°, preventing the development of a second plume.In this case (moderate T value), the cavity inclination affects the charge density distribution topography.As the angle α reaches 90° a shift occurs from a single-plume pattern to a bifurcated-plume arrangement situated at both the upper and lower regions of the cavity.For α = 0°, the charge distribution aligns closely with the findings reported by Hassen et al. [36] and Traore et al.The local Nusselt number variation over the hot wall according to the tilt angle is plotted in Figure 5.At T = 10, heat exchange occurs mainly in the lower part of the heated wall.For y < 0.4, the heat transfer is more significant as the inclination angle decreases.In fact, for y = 0.1, a difference of 45% in Nu number is observed between α = 0° and α = 90°, indicating the impact of the inclination angle on heat transfer.However, in the upper half of the hot wall, the heat transfer is slightly affected by the cavity inclination.By increasing the electric field strength, the impact of the tilt angle becomes less important, even in the lower part of the hot plate.For T = 600, regardless of the tilt angle, almost all the curves overlap.This was expected due to the intensification of the electric field, which reaches a level where the flow regime is electrically dominated, and it becomes the only parameter controlling the heat transfer process.The local Nusselt number variation over the hot wall according to the tilt angle is plotted in Figure 5.At T = 10, heat exchange occurs mainly in the lower part of the heated wall.For y < 0.4, the heat transfer is more significant as the inclination angle decreases.In fact, for y = 0.1, a difference of 45% in Nu number is observed between α = 0 • and α = 90 • , indicating the impact of the inclination angle on heat transfer.However, in the upper half of the hot wall, the heat transfer is slightly affected by the cavity inclination.By increasing the electric field strength, the impact of the tilt angle becomes less important, even in the lower part of the hot plate.For T = 600, regardless of the tilt angle, almost all the curves overlap.This was expected due to the intensification of the electric field, which reaches a level where the flow regime is electrically dominated, and it becomes the only parameter controlling the heat transfer process.
The mean Nusselt number according to the tilt cavity for various T values is shown in Figure 6.Two distinct behaviors can be recorded for T ≤ 200, as depicted in Figure 6a.The profile of the mean Nusselt number follows a parabolic shape with a maximum around an angle of 45 • .It is interesting to report that the outcomes for T = 10 and T = 0 are identical.This similarity arises because, for a low Rayleigh number, T < Tc below the critical Rayleigh number (commonly referred to as the linear stability criterion [37]) the electric field has an insignificant impact on natural convection.Increasing T leads to a thermal transfer enhancement of approximately 13% between T = 0 and 200.For high values of T (greater than 400), a reduction in Nu can be observed between 200 and 400 due to bifurcations from a single-cell flow to a two-cell flow.Further, increasing T results in the intensification of the average Nusselt number, reaching an improvement of 61% between T = 0 and T = 800.It should also be noted that for high electrical Rayleigh values, the profile of the average Nusselt number exhibits horizontal lines, indicating that the electrical forces are so dominant that the cavity tilt has no influence on the heat exchange process, and the only controlling parameter is the electrical force.The mean Nusselt number according to the tilt cavity for various T values is shown in Figure 6.Two distinct behaviors can be recorded for T ≤ 200, as depicted in Figure 6a.The profile of the mean Nusselt number follows a parabolic shape with a maximum around an angle of 45°.It is interesting to report that the outcomes for T = 10 and T = 0 are identical.This similarity arises because, for a low Rayleigh number, T < Tc below the critical Rayleigh number (commonly referred to as the linear stability criterion [37]) the electric field has an insignificant impact on natural convection.Increasing T leads to a thermal transfer enhancement of approximately 13% between T = 0 and 200.For high values of T (greater than 400), a reduction in Nu can be observed between 200 and 400 due to bifurcations from a single-cell flow to a two-cell flow.Further, increasing T results in the intensification of the average Nusselt number, reaching an improvement of 61% between T = 0 and T = 800.It should also be noted that for high electrical Rayleigh values, the profile of the average Nusselt number exhibits horizontal lines, indicating that the electrical forces are so dominant that the cavity tilt has no influence on the heat exchange process, and the only controlling parameter is the electrical force.The mean Nusselt number according to the tilt cavity for various T values is shown in Figure 6.Two distinct behaviors can be recorded for T ≤ 200, as depicted in Figure 6a.The profile of the mean Nusselt number follows a parabolic shape with a maximum around an angle of 45°.It is interesting to report that the outcomes for T = 10 and T = 0 are identical.This similarity arises because, for a low Rayleigh number, T < Tc below the critical Rayleigh number (commonly referred to as the linear stability criterion [37]) the electric field has an insignificant impact on natural convection.Increasing T leads to a thermal transfer enhancement of approximately 13% between T = 0 and 200.For high values of T (greater than 400), a reduction in Nu can be observed between 200 and 400 due to bifurcations from a single-cell flow to a two-cell flow.Further, increasing T results in the intensification of the average Nusselt number, reaching an improvement of 61% between T = 0 and T = 800.It should also be noted that for high electrical Rayleigh values, the profile of the average Nusselt number exhibits horizontal lines, indicating that the electrical forces are so dominant that the cavity tilt has no influence on the heat exchange process, and the only controlling parameter is the electrical force.Figure 7 plots the evolution of the local Nusselt number with respect to the cavity inclination angle for different Rayleigh numbers.It can be shown that convective heat exchange is more significant when the thermal Rayleigh number is higher.An intensification of up to 180% is achieved at an angle of 15° just by increasing Ra = 5000 to 250,000.However, we observe that the Nusselt number is only affected by inclination when the thermal Rayleigh number is relatively high (Ra > 10 5 ).Moreover, for high Ra values, the optimal heat exchanges are reached at 15-30° cavity inclination, while a reduction of 34% occurs at an angle of 90°.After conducting several numerical simulations, a multi-parameter linear regression analysis was performed to establish an empirical correlation.Although several correlations have been presented in the literature for estimating the Nusselt number in the case of classical free convection in inclined cavities, such as those suggested by Wright et al. [38], ElSherbiny et al. [39], and Khezzar et al. [40].As far as our current knowledge extends, no correlation has yet been suggested for the case of electro-thermo convection,  Figure 7 plots the evolution of the local Nusselt number with respect to the cavity inclination angle for different Rayleigh numbers.It can be shown that convective heat exchange is more significant when the thermal Rayleigh number is higher.An intensification of up to 180% is achieved at an angle of 15° just by increasing Ra = 5000 to 250,000.However, we observe that the Nusselt number is only affected by inclination when the thermal Rayleigh number is relatively high (Ra > 10 5 ).Moreover, for high Ra values, the optimal heat exchanges are reached at 15-30° cavity inclination, while a reduction of 34% occurs at an angle of 90°.After conducting several numerical simulations, a multi-parameter linear regression analysis was performed to establish an empirical correlation.Although several correlations have been presented in the literature for estimating the Nusselt number in the case of classical free convection in inclined cavities, such as those suggested by Wright et al. [38], ElSherbiny et al. [39], and Khezzar et al. [40].As far as our current knowledge extends, no correlation has yet been suggested for the case of electro-thermo convection, taking into account the inclination angle, heat transfer, and electric field.After conducting several numerical simulations, a multi-parameter linear regression analysis was performed to establish an empirical correlation.Although several correlations have been presented in the literature for estimating the Nusselt number in the case of classical free convection in inclined cavities, such as those suggested by Wright et al. [38], ElSherbiny et al. [39], and Khezzar et al. [40].As far as our current knowledge extends, no correlation has yet been suggested for the case of electro-thermo convection, taking into account the inclination angle, heat transfer, and electric field.
The correlation proposed in Equation ( 17) offers an approximate estimation of the Nusselt number, which is applicable for a wide range of thermal Rayleigh numbers, spanning from 5000 to 250,000, electrical Rayleigh numbers ranging from 10 to 600, and tilt angles varying between 0 • and 80 • .NU = 0.221 Ra 0.242 T 0.061 exp(−0.011tan(α)). ( This correlation exhibits a coefficient of determination R 2 = 0.95, a standard error σ = 0.06, and a p-value for the different parameters ranging from 0 to 0.007.All these statistical parameters, such as the p-value, standard deviation, and coefficient of determination, indicating the effectiveness and reliability of the proposed correlation. To further confirm the validity of the correlation, comparison tests were carried out in Figure 8.These tests compare the mean values of the Nusselt number (Nu) simulated and given by Equation ( 20) for several influential parameters (Ra, T, and α).The excellent agreement between simulated and correlated mean Nusselt numbers highlights the good reliability of the developed correlation.Indeed, the maximum error does not exceed 4%, confirming the accuracy of Equation (20).
Processes 2023, 11, x FOR PEER REVIEW 13 of 1 The correlation proposed in Equation ( 17) offers an approximate estimation of th Nusselt number, which is applicable for a wide range of thermal Rayleigh numbers, span ning from 5000 to 250,000, electrical Rayleigh numbers ranging from 10 to 600, and til angles varying between 0° and 80°.
This correlation exhibits a coefficient of determination R 2 = 0.95, a standard error σ 0.06, and a p-value for the different parameters ranging from 0 to 0.007.All these statistica parameters, such as the p-value, standard deviation, and coefficient of determination, in dicating the effectiveness and reliability of the proposed correlation.
To further confirm the validity of the correlation, comparison tests were carried ou in Figure 8.These tests compare the mean values of the Nusselt number (Nu) simulated and given by Equation ( 20) for several influential parameters (Ra, T, and α).The excellen agreement between simulated and correlated mean Nusselt numbers highlights the good reliability of the developed correlation.Indeed, the maximum error does not exceed 4% confirming the accuracy of Equation (20).

Conclusions
In this work, numerical simulations of elelctro-thermo convection in an inclined dif ferential-heated square cavity have been carried out.The cavity is filled with dielectri fluid and subjected to an electric field.A parametric study was conducted to investigate the effects of enclosure inclination, imposed thermal gradient, and applied electric field The key conclusions of this study can be outlined in the following:


The flow's behavior can be significantly impacted by the inclination of the cavity.By increasing this angle, a shift from a unicellular to a bicellular regime was detected.


At high electric and thermal fields, when the angle α ∈ [0.40°], increasing the tilt angl results in a 5% improvement in heat transfer.Above this angular range, an oscillation regime could be developed.


Local Nusselt numbers show that heat exchange occurs mainly in the lower section of the heated wall.A reduction in the convective heat transfer of up to 45% can b

Conclusions
In this work, numerical simulations of elelctro-thermo convection in an inclined differential-heated square cavity have been carried out.The cavity is filled with dielectric fluid and subjected to an electric field.A parametric study was conducted to investigate the effects of enclosure inclination, imposed thermal gradient, and applied electric field.
The key conclusions of this study can be outlined in the following: - The flow's behavior can be significantly impacted by the inclination of the cavity.By increasing this angle, a shift from a unicellular to a bicellular regime was detected.

Figure 3 .
Figure 3. (a) Time profile of the average Nusselt number; (b) local Nusselt number on the hot plate for Ra = 50,000 and T = 600 and for various inclinations.

Figure 3 .
Figure 3. (a) Time profile of the average Nusselt number; (b) local Nusselt number on the hot plate for Ra = 50,000 and T = 600 and for various inclinations.

Figure 3 .
Figure 3. (a) Time profile of the average Nusselt number; (b) local Nusselt number on the hot plate for Ra = 50,000 and T = 600 and for various inclinations.

Figure 4 .
Figure 4. Distributions of charge density for different tilt angles and for different values of the electric Rayleigh number, Ra = 10,000, Pr = 10, C = 10, and M = 10.

Figure 7 Figure 6 .
Figure7plots the evolution of the local Nusselt number with respect to the cavity inclination angle for different Rayleigh numbers.It can be shown that convective heat exchange is more significant when the thermal Rayleigh number is higher.An intensification of up to 180% is achieved at an angle of 15 • just by increasing Ra = 5000 to 250,000.However, we observe that the Nusselt number is only affected by inclination when the thermal Rayleigh number is relatively high (Ra > 10 5 ).Moreover, for high Ra values, the optimal heat exchanges are reached at 15-30 • cavity inclination, while a reduction of 34% occurs at an angle of 90 • .

Figure 8 .
Figure 8.Average Nusselt number according to tilt cavity angle for various Ra and T numbers lines Correlation values; Symbols: Numerical results.

Figure 8 .
Figure 8.Average Nusselt number according to tilt cavity angle for various Ra and T numbers lines: Correlation values; Symbols: Numerical results.