Thermophysical Properties of Nanofluid in Two-Phase Fluid Flow through a Porous Rectangular Medium for Enhanced Oil Recovery

It is necessary to sustain energy from an external reservoir or employ advanced technologies to enhance oil recovery. A greater volume of oil may be recovered by employing nanofluid flooding. In this study, we investigated oil extraction in a two-phase incompressible fluid in a two-dimensional rectangular porous homogenous area filled with oil and having no capillary pressure. The governing equations that were derived from Darcy’s law and the mass conservation law were solved using the finite element method. Compared to earlier research, a more efficient numerical model is proposed here. The proposed model allows for the cost-effective study of heating-based inlet fluid in enhanced oil recovery (EOR) and uses the empirical correlations of the nanofluid thermophysical properties on the relative permeability equations of the nanofluid and oil, so it is more accurate than other models to determine the higher recovery factor of one nanoparticle compared to other nanoparticles. Next, the effect of nanoparticle volume fraction on flooding was evaluated. EOR via nanofluid flooding processes and the effect of the intake temperatures (300 and 350 K) were also simulated by comparing three nanoparticles: SiO2, Al2O3, and CuO. The results show that adding nanoparticles (<5 v%) to a base fluid enhanced the oil recovery by more than 20%. Increasing the inlet temperature enhanced the oil recovery due to changes in viscosity and density of oil. Increasing the relative permeability of nanofluid while simultaneously reducing the relative permeability of oil due to the presence of nanoparticles was the primary reason for EOR.


Introduction
New technologies are being developed to boost oil production to comply with rising energy demand [1], while continuously depleting resources. Oil recovery processes are utilized to obtain maximum oil from hydrocarbon reservoirs [2] under economic feasibility constraints. Primary recovery refers to extracting oil from a reservoir using natural energy and pressure. Because the inherent energy of the reservoir declines throughout the production process, only 5% to 15% of hydrocarbons may be recovered by using this recovery method.
Enhanced oil recovery (EOR) is implemented to improve residual oil extraction from the homogenous medium in the reservoir [3]. The secondary and the tertiary recoveries are included in the EOR process, typically constituting sophisticated technologies. The secondary EOR process is achieved through water injection or gas injection [4]. After primary recovery, fluid injection is used to restore the reservoir pressure. In contrast, the fluid injection in the tertiary recovery alters the interaction between the reservoir rock and the fluid [5]. The displacement of oil from porous media using water or solutions including surfactants or polymer additives is a popular method of extracting oil from terrigenous sources [6].
Silicon SiO 2 , aluminum oxide Al 2 O 3 , and copper oxide CuO nanoparticles are the most frequently utilized nanoparticles in EOR. As a result of their characteristics, nanoparticles can address a few issues with conventional approaches and improve the displacing fluid's sweep efficiency, because the base fluids' thermal conductivity and interfacial characteristics are both improved by nanoparticles. Furthermore, they can withstand high temperatures and pressures, are ecologically friendly and are cost-effective. Almost all EOR approaches are expensive, but the advantages of nanofluid flooding can be more economical. In addition, nanoparticles can optimize the properties of displacing fluids such as density, viscosity, the interfacial tension between oil and water, thermal conductivity, and specific heat [2].
Several papers have recently been published in which nanofluid flooding is proposed as an agent for EOR [27][28][29][30][31]. The effect of nanofluid flooding on enhancing oil extraction from porous media has also been reported. Despite several experimental and modeling investigations, there is currently no exact explanation for the impacts of nanoparticles on the rise in oil recovery during nanofluid flooding, which limits the adoption of the method in the industry. Some suggested theories for understanding the influence of nanofluids on oil extraction include surface tension-lowering techniques [32][33][34][35], the alteration of the porous medium's wettability [36][37][38], and the change in viscosity [39][40][41][42][43][44].
Researchers have looked at the transfer of microscopic particles in a two-phase flow through an open porous medium to simulate the EOR process [6,[45][46][47][48]. Polysilicon nanoparticles were studied by Ju and Fan [49]. Their results showed that the environment shifts from oily to watery, promoting oil extraction. They developed a mathematical model for small particle transport in a two-phase flow in porous media [50].
The presence of nanoparticles also has a considerable impact on the thermophysical properties of the base fluid [51][52][53][54][55]. However, a limited number of research studies are available on thermophysical properties of nanofluid in porous media for EOR [54]. Furthermore, nanoparticles seem to enhance the effective characteristics of the injected fluid, particularly in the process of washing the oil from a porous medium [56,57].
In the present study, a numerical method of finite element is used to solve the governing equations that were derived from Darcy's law and the mass conservation law of a two-phase incompressible fluid in a two-dimensional porous homogenous area. The thermophysical properties of nanoparticles were studied for their volume fraction and their kind and influence on the extraction of oil. Compared to earlier research, a more efficient numerical model is proposed. The proposed model allows for the cost-effective study of heating-based inlet fluid in EOR and uses the empirical correlations of the nanofluid thermophysical properties on the relative permeability equations of the nanofluid and oil, so it is more accurate than other models to determine the higher recovery factor of one nanoparticle compared to other nanoparticles. Moreover, the EOR via nanofluid flooding processes and the effect of the intake temperatures (300 and 350 K) were also simulated by comparing three nanoparticles, SiO 2 , Al 2 O 3 and CuO, with different nanoparticle volume fractions. Finally, we studied how increasing the inlet temperature enhanced the oil recovery due to changes in viscosity and density of oil. This work aims to model the EOR procedure by flooding, taking into consideration the efficient characteristics of nanofluids conveyed through a model porous medium.

Numerical Implementation
Computational fluid dynamics (CFDs) is a methodology used to examine systems, such as fluid flow and the heat transport system, through computer simulations. When compared to experimental investigations, the CFD technique has the potential to study critical and unique circumstances in a process, reduce response time and research costs, and gain a complete and detailed understanding of the process [58][59][60][61][62][63][64].
FEM techniques, which are one of the most used approaches for CFD research, were employed here. The primary framework of the employed simulation software COMSOL Multiphysics is fully based on FEM as a numerical method. The numerical analysis community frequently uses finite element methods to investigate numerical approaches for fluid flow.
There is a lot of work on CFD and finite element techniques in literature, particularly in the EOR process and porous medium modeling [54,[65][66][67][68]. Four phenomena occur when nanofluid is introduced into porous media: adsorption, desorption, blockage, and migration with flowing fluid. This section introduces one proposed model that others have recently adopted for nanofluid injection in oil reservoirs for EOR [69][70][71][72].

Model Assumptions
The following assumptions were used to build the two-phase flow mathematical model with distributed nanoparticles [73]:

•
Fluids were both incompressible. • Fluid flow in pore medium obeyed Darcy's law.

•
The flow was in an isothermal state.

•
There was no capillary pressure such that permeability and porosity were independent of time and space and that the impact of gravity was neglected. • Darcy's law and the mass conservation law provided the essential equations. Mass conservation was comparable to volume conservation since the two phases are incompressible.

Geometry Creation
The geometry of the model was produced using COMSOL Multiphysics software with the following attributes: a 2D porous medium with porosity of 0.3, absolute permeability of 1e − 9 m 2 , and size of 1 × 2 m 2 . Figure 1 indicates that the porous medium was homogenous for its distribution.   The boundary layer mesh is important because the default triangular mesh in COM-SOL software does not work well for fluid flow problems, and there is a boundary layer along the contact between the nanofluid and the solid, as shown in Figure 2a. A no-slip condition from our assumption generated the boundary layer of the nanofluid flow through a stationary porous medium, and the layer close to the solid has zero velocity. Figure 2b below depicts the development of a boundary layer mesh for achieving a refined solution to the fluid flow problems by using extremely small, thin components near the solid wall.

Grid Independency
Varying grid types with different mesh sizes were developed for the investigation of grid independence. To compare grids, fluid pressure drops in the porous medium were considered. The total pressure at the entrance and exit were estimated. The mesh size, number of elements generated, as well as the pressure drop between the intake and exit are all shown in Table 1. As a result, mesh number 3 is chosen as the ideal mesh number and is utilized in all the simulations that follow.
Although COMSOL Multiphysics meshing algorithms attempt to eliminate lowquality parts, this is not always achievable for all geometries. A quality of 1 is the greatest feasible value for any quality measure and signifies an ideal element within the specified quality measure. At the interval's opposite end, 0 denotes a degraded element. As a result of the element measurements, the meshing of the present research geometry is of excellent quality, with a quality of 1, as shown in Figure 3a.
In addition, Figure 3b depicts the plot of convergence, which is time-domain convergence for time-dependent problems (variables), and finds a solution to the system of equations with error less than 1 − 5 at a specific time, 250 s.

Grid Independency
Varying grid types with different mesh sizes were developed for the investigation of grid independence. To compare grids, fluid pressure drops in the porous medium were considered. The total pressure at the entrance and exit were estimated. The mesh size, number of elements generated, as well as the pressure drop between the intake and exit are all shown in Table 1. As a result, mesh number 3 is chosen as the ideal mesh number and is utilized in all the simulations that follow. Although COMSOL Multiphysics meshing algorithms attempt to eliminate lowquality parts, this is not always achievable for all geometries. A quality of 1 is the greatest feasible value for any quality measure and signifies an ideal element within the specified quality measure. At the interval's opposite end, 0 denotes a degraded element. As a result of the element measurements, the meshing of the present research geometry is of excellent quality, with a quality of 1, as shown in Figure 3a.

Governing Equations
The equations associated with the mathematical modeling of flooding problems using nanofluids accounted for the heat transfer. The nanofluid flooding equations were first studied as water flooding equations by assuming that the nanoparticle's concentration equals 0 for examining the mathematical and physical equations of the problems in the oil reservoir. The porous medium's heat equation was then supplied to accomplish the model.
Because the initial pressure difference between the well and the tank was extremely large, the oil exited the tank automatically until the well and tank pressures equalized, at which time recovery procedures were required. Because it was considered that flow cannot exit boundaries 2 and 3, and that flow can only move between boundaries 1 and 4, the zero-flow boundary condition was considered.
In recent years, nanofluids have been used instead of water to increase oil recovery and to improve the flooding function [26]. To determine the rate of oil recovery, the equations of saturated oil ( o ) and water-based nanofluid saturation ( nf ) in the environment should be determined. The porous medium was first filled with oil. Waterbased nanofluid entered the porous medium from boundary 1 with a velocity nf0 of 0.001 m/s, and displaced oil, which was permitted to flow out of the porous domain via boundary 4. At the inlet (boundary 1), the velocity o of oil was 0: The pressure at the outlet (boundary 4) was set to 0 Pa: The total velocity = nf + o was constant in time and space because of the assumptions and boundary conditions stated above, and Darcy's law and the mass conservation law provide the essential equations. Mass conservation is comparable to volume conservation since the two phases were assumed to be incompressible.
For the two phases of water-based nanofluid and oil, Darcy's law yields [74]: In addition, Figure 3b depicts the plot of convergence, which is time-domain convergence for time-dependent problems (variables), and finds a solution to the system of equations with error less than 1e − 5 at a specific time, 250 s.

Governing Equations
The equations associated with the mathematical modeling of flooding problems using nanofluids accounted for the heat transfer. The nanofluid flooding equations were first studied as water flooding equations by assuming that the nanoparticle's concentration equals 0 for examining the mathematical and physical equations of the problems in the oil reservoir. The porous medium's heat equation was then supplied to accomplish the model.
Because the initial pressure difference between the well and the tank was extremely large, the oil exited the tank automatically until the well and tank pressures equalized, at which time recovery procedures were required. Because it was considered that flow cannot exit boundaries 2 and 3, and that flow can only move between boundaries 1 and 4, the zero-flow boundary condition was considered.
In recent years, nanofluids have been used instead of water to increase oil recovery and to improve the flooding function [26]. To determine the rate of oil recovery, the equations of saturated oil (S o ) and water-based nanofluid saturation (S nf ) in the environment should be determined. The porous medium was first filled with oil. Water-based nanofluid entered the porous medium from boundary 1 with a velocity u nf0 of 0.001 m/s, and displaced oil, which was permitted to flow out of the porous domain via boundary 4. At the inlet (boundary 1), the velocity u o of oil was 0: The pressure at the outlet (boundary 4) was set to 0 Pa: The total velocity u = u nf + u o was constant in time and space because of the assumptions and boundary conditions stated above, and Darcy's law and the mass conservation law provide the essential equations. Mass conservation is comparable to volume conservation since the two phases were assumed to be incompressible.
For the two phases of water-based nanofluid and oil, Darcy's law yields [74]: and where the velocities of phase 1 water-based nanofluid and phase 2 oil are u nf and u o , respectively; the mobilities of phases 1 and 2 are m nf and m o , respectively; m nf = k rnf /µ nf ; and m o = k ro /µ o . The functions k nf and k o are dependent on the saturation, and K is the absolute permeability. The laws of volume conservation are: where the saturation of phase 1 (nanofluid) is S nf and of phase 2 (oil) is S o , while the porosity is ε. The initial condition was taken as S nf = 0 to solve saturation equations. In the reservoir, there was often not just oil, but also some water. In this study, for first saturation, a quantity of 0 was considered. Because the input fluid contains no oil, the inlet nanofluid saturation (S nf ) was 1. It should also be noted that there is a correlation between the amount of oil saturation and nanofluid flooding [75]: Equation (9) was obtained by adding Equations (6) and (7) and by considering the fluids incompressible, and by using Equation (8). and The entire volumetric flow (velocity) is denoted by u. In this study, we simplify by assuming that u is constant.
Since water-based nanofluid is liquid, the values of effective density and viscosity of nanofluid to solve Equation (13) should be utilized. Equations (14) and (15) were used to derive the effective density and viscosity, respectively [76]: where ρ nf and µ nf are the density and viscosity of nanofluid, respectively, and ∅ is the nanoparticle volume fraction of the nanofluid. Relative permeability was required for water-based nanofluid and oil for solving the saturation equations. The effective saturation (S e ) was utilized in the relative permeabilities: where (S e ) is the effective saturation, (S rnf ) is the residual nanofluid saturation, and (S ro ) is the residual oil saturation. Equation (17) [77] defines the relative permeability of waterbased nanofluid and oil in the Brooks-Corey model: where λ is the distribution index for the pore size, and its value is 1. The capillary pressure or the difference in pressure between nanofluid and oil equals 0 since it was assumed that there is no capillary pressure in this model. It can then be expressed as follows. then The oil temperature in the reservoir increased at a higher heat transfer rate between the two phases when the intake fluid was high. When nanoparticles were added, the thermal conductivity coefficient k e f f increased, allowing for more efficient heat transmission involving nanofluids and oil. Under this situation, the flooding approach would benefit from two factors as temperatures rise. Lowering the viscosity of the oil eased the nanofluid for replacing the oil from the porous medium since its stickiness and heavyweight were reduced, and it drained faster than the same water flow. Another effect of rising temperatures was that oil density decreased, making it lighter in volume and needing a smaller amount of energy to be evacuated from reservoirs.
The energy equation in porous media is introduced in this section [78]: ρC In the environment, the initial temperature is 300 K. The border condition is specified as Equation (22) for the input of the heat equation: The target variable for this study is T in , denoting the temperature of the input fluid and influencing the performance of the recovery process. The condition of the foreign border 4 is: Table 2 illustrates the properties of water at a constant temperature. The four fluid parameters used in the two-phase and energy equations are density, viscosity, specific heat capacity, and conduction heat transfer coefficient. It also summarizes the results of past studies on oil characteristics. With rising temperatures, the density and viscosity of water and oil change. Transferring water heat to oil decreases the viscosity and density of oil. Therefore, the oil could be transported more readily and recovered. This approach must be investigated in this study. According to Equation (24), the density of water varies with temperature (for temperatures ranging from 300 K to 350 K) [79]. Table 3 also shows how the viscosity of water changes with temperature.  Table 3. Effects of temperature on viscosity of water [79].

273
.15 < T < 413.15 After assessing how temperature affects the density and viscosity of water, the mechanism of temperature in affecting the density and viscosity of the oil was then evaluated. There are data about oil density at three temperatures. Oil characteristics, like water, are temperature-dependent. The change in oil density with temperature is shown in Table 4. Similarly, oil viscosity operates in the same way. Table 5 [ 82] shows the oil's viscosity at three distinct temperatures. Table 4. Temperature-related changes in oil density [82].  As the oil temperature inside the reservoir rises, the viscosity and the density of the oil gradually reduce, as indicated in Tables 4 and 5. The oil becomes smoother and easier to evacuate from the porous medium. Therefore, lowering these parameters would also aid in EOR.

The Nano-Particles' Effect on Thermophysical Properties of Nanofluid
The nanoparticles SiO 2 , Al 2 O 3 , and CuO were evaluated in this study. These nanoparticles were combined with water at specific compositions to generate a nanofluid that flows into the porous medium.
The thermophysical parameters of SiO 2 , Al 2 O 3 , and CuO particles are shown in Table 6. By using the numerical values of these thermophysical parameters and the nanofluid empirical equations in Table 7, the nanofluid characteristics were extracted, and the effects of adding various nanoparticles on the physical characteristics were examined. Furthermore, for all three nanofluids, the influence of a volume fraction ranging from 1% to 5% on the thermophysical characteristics of the base fluid was investigated.  Table 7. Properties of nanofluids.

Nanofluid Properties Equation (Combining Nanoparticle and Water) Unit
Density (ρ nf ) [86] Specific heat capacity C P nf ) [87] The impact of the volume fraction on the density of water nanofluids SiO 2 , Al 2 O 3 , and CuO is shown in Figure 4a. The density of all three types of nanofluids increased as the volume fraction increased, with the increase being greater for CuO water nanofluids than SiO 2 water and Al 2 O 3 nanofluids. Due to the lowest density of SiO 2 nanoparticles compared to the others, SiO 2 water showed the lowest nanofluid density. Furthermore, for all three nanofluids, the influence of a volume fraction ranging from 1% to 5% on the thermophysical characteristics of the base fluid was investigated.  Table 7. Properties of nanofluids.

Nanofluid Properties Equation (Combining Nanoparticle and Water) Unit
Density (ρ nf ) [86] ρ nf = ∅ + (1 − ∅) kg m 3 Specific heat capacity (C P nf ) [87] C P nf = ∅ + (1 − ∅) J kg. K Thermal conductivity coefficient (K nf ) [88,89] K nf = + 2 − 2∅( − ) [90,91] μ nf = (1 + 39.11∅ + 533.9∅ 2 ) Pa. s The impact of the volume fraction on the density of water nanofluids SiO 2 , Al 2 O 3 , and CuO is shown in Figure 4a. The density of all three types of nanofluids increased as the volume fraction increased, with the increase being greater for CuO water nanofluids than SiO 2 water and Al 2 O 3 nanofluids. Due to the lowest density of SiO 2 nanoparticles compared to the others, SiO 2 water showed the lowest nanofluid density.    Figure 4b depicts the effect of increasing the amount of nanoparticles in the base fluid on thermal capacity. It demonstrates a reduction in thermal capacity as the ∅ value increases. The analysis of the three nanofluid combinations revealed that SiO 2 in the base water had the maximum thermal capacity. The thermal capacity varies depending on the nanofluid density and the nanoparticles' thermal capacity. Figure 5a demonstrates how increasing the concentration of nanoparticles ∅ influences the nanofluid's viscosity. Using the viscosity equation (see Table 7), it is clear that the viscosity was only impacted by the parameter ∅ and viscosity of the base fluid and was unaffected by the kind of nanoparticle. The viscosity increased with increasing concentration for all three nanoparticles. Figure 5b shows the impact of increasing the concentration of nanoparticles on the nanofluid's thermal conductivity. As the concentration of nanoparticles increases, so does the heat conductivity. The greatest thermal conductivity was found in an Al 2 O 3 water nanofluid, followed by SiO 2 and CuO. and was unaffected by the kind of nanoparticle. The viscosity increased with increasing concentration for all three nanoparticles. Figure 5b shows the impact of increasing the concentration of nanoparticles on the nanofluid's thermal conductivity. As the concentration of nanoparticles increases, so does the heat conductivity. The greatest thermal conductivity was found in an Al 2 O 3 water nanofluid, followed by SiO 2 and CuO.
The fluid characteristics were a mixture of nanofluids and oils, and their saturation must also be considered in the case of heat transmission. Table 8 shows the thermophysical properties of the energy equation.
The reservoir rock characteristics were required in the energy calculation for a porous medium. The Rafaat et al. [73] model was adopted because the relative permeability requires the diameter of the granules inside the reservoir; these characteristics are provided in Table 9.   The fluid characteristics were a mixture of nanofluids and oils, and their saturation must also be considered in the case of heat transmission. Table 8 shows the thermophysical properties of the energy equation. Table 8. Total characteristics that may be used in the energy equation.

Nanofluid Properties Equation (Combining Nanoparticle and Water) Unit
Specific heat capacity Cp tot The reservoir rock characteristics were required in the energy calculation for a porous medium. The Rafaat et al. [73] model was adopted because the relative permeability requires the diameter of the granules inside the reservoir; these characteristics are provided in Table 9. Table 9. Properties of reservoir rocks [92].

Verification
The simulation was compared to the findings of an experimental study provided by Maghzi et al. [78] to validate the numerical solver. In their experiment, the silica nanoparticles in distilled water were employed to increase the amount of oil recovered from the porous medium. The geometry of their model featured a 2D porous medium with porosity of 0.33, absolute permeability of 200 mD, and size of 6 × 6 cm 2 . In this work, the results are provided utilizing the COMSOL solver. In Figure 6, the experimental data are compared to the model findings for changes in the amount of oil retrieved from the medium versus pore volume and exhibited a good agreement with experimental findings [93].
Specific heat capacity (C P ) 851 kg. K Thermal conductivity coefficient ( ) 2.2 W m. K Diameter of particle 3 nm

Verification
The simulation was compared to the findings of an experimental study provided by Maghzi et al. [78] to validate the numerical solver. In their experiment, the silica nanoparticles in distilled water were employed to increase the amount of oil recovered from the porous medium. The geometry of their model featured a 2D porous medium with porosity of 0.33, absolute permeability of 200 mD, and size of 6 × 6 cm 2 . In this work, the results are provided utilizing the COMSOL solver. In Figure 6, the experimental data are compared to the model findings for changes in the amount of oil retrieved from the medium versus pore volume and exhibited a good agreement with experimental findings [93]. Figure 6. Evaluation of curve changes in the amount of oil recovered from porous medium computed from the model and experimental data [93].

Results
The solution domain and porous medium were considered homogenous in this study. Our research focus was the impacts of nanofluids in water flooding applications. The simulations were run twice, once for water and again for nanofluid inlet flow. As previously stated, empirical research has shown that this technique significantly impacts petroleum outflow. The purpose of this study was to validate and analyze its relevance and impact on EOR. The problem was initially addressed with a concentration of 0 (purified water) and then by contemplating the 5% nanoparticle concentration with water, assuming silicone nanoparticles and water fluid usage.
As demonstrated in Figure 7a-c, the quantity of nanofluid saturation in the environment was significantly larger and the extracted oil from porous media was

Results
The solution domain and porous medium were considered homogenous in this study. Our research focus was the impacts of nanofluids in water flooding applications. The simulations were run twice, once for water and again for nanofluid inlet flow. As previously stated, empirical research has shown that this technique significantly impacts petroleum outflow. The purpose of this study was to validate and analyze its relevance and impact on EOR. The problem was initially addressed with a concentration of 0 (purified water) and then by contemplating the 5% nanoparticle concentration with water, assuming silicone nanoparticles and water fluid usage.
As demonstrated in Figure 7a-c, the quantity of nanofluid saturation in the environment was significantly larger and the extracted oil from porous media was boosted. In Figure 7a-c, the red contour depicts the porous medium region where the trapped oil was replaced by the nanofluid, and the oil was extracted. The blue area of the contour refers to a porous medium region that held oil that had not yet been released. boosted. In Figure 7a-c, the red contour depicts the porous medium region where the trapped oil was replaced by the nanofluid, and the oil was extracted. The blue area of the contour refers to a porous medium region that held oil that had not yet been released. As demonstrated in Figure 8a-c, the quantity of oil saturation in the environment was significantly lower and the nanofluid inserted into porous media was boosted. In Figure 8a-c, the blue contour depicts the porous medium region where the trapped oil was replaced by the nanofluid, and the oil was extracted. The red area of the contour refers to a porous medium region that has held oil that has not yet been released.
As seen in Figure 9a-c, the quantity of pressure in the porous medium did not change significantly as a result of pressure varying with velocity. Since Darcy's law is only applicable for slow flow, most groundwater flow cases fall into this category. In this investigation, the velocity nf0 was 0.001 m/s, implying a Reynolds number smaller than unity, indicating that the flow was laminar, and that Darcy's equation was used. Figure 10 presents the oil recovery factors for flooding of water and nanofluids. This figure shows that flooding by nanofluids was more effective than flooding by water. As demonstrated in Figure 10, adding nanoparticles to the base fluid improved the recovery factor. The addition of nanoparticles to the base fluid increased the relative weight of the input fluid, increasing energy and allowing for more efficient oil production. Furthermore, because of the nanoparticles' strong thermal transfer characteristics, the heat was transferred to the oil inside the reservoir more quickly. Oil recovery was improved due to the decreased viscosity and density of reservoir oil.
Another critical parameter in EOR and the solution of two-phase flows is the relative permeability of the two fluids. In the current study, the relative permeability of water and oil fluid was altered when nanoparticles were added to the water. As demonstrated in Figure 11, raising the relative permeability of nanofluid while simultaneously reducing the permeability of oil due to the presence of nanoparticles were the primary reasons for EOR. Because the relative permeability of the water was enhanced as a result of this procedure, it drove the oil flow better. In conclusion, employing nanofluids instead of water increased the performance of the EOR process.  As demonstrated in Figure 8a-c, the quantity of oil saturation in the environment was significantly lower and the nanofluid inserted into porous media was boosted. In Figure 8a-c, the blue contour depicts the porous medium region where the trapped oil was replaced by the nanofluid, and the oil was extracted. The red area of the contour refers to a porous medium region that has held oil that has not yet been released.   As seen in Figure 9a-c, the quantity of pressure in the porous medium did not change significantly as a result of pressure varying with velocity. Since Darcy's law is only applicable for slow flow, most groundwater flow cases fall into this category. In this investigation, the velocity u nf0 was 0.001 m/s, implying a Reynolds number smaller than unity, indicating that the flow was laminar, and that Darcy's equation was used.    Figure 10 presents the oil recovery factors for flooding of water and nanofluids. This figure shows that flooding by nanofluids was more effective than flooding by water. As demonstrated in Figure 10, adding nanoparticles to the base fluid improved the recovery factor. The addition of nanoparticles to the base fluid increased the relative weight of the input fluid, increasing energy and allowing for more efficient oil production. Furthermore, because of the nanoparticles' strong thermal transfer characteristics, the heat was transferred to the oil inside the reservoir more quickly. Oil recovery was improved due to the decreased viscosity and density of reservoir oil.  Another critical parameter in EOR and the solution of two-phase flows is the relative permeability of the two fluids. In the current study, the relative permeability of water and oil fluid was altered when nanoparticles were added to the water. As demonstrated in Figure 11, raising the relative permeability of nanofluid while simultaneously reducing the permeability of oil due to the presence of nanoparticles were the primary reasons for EOR. Because the relative permeability of the water was enhanced as a result of this procedure, it drove the oil flow better. In conclusion, employing nanofluids instead of water increased the performance of the EOR process. There are several ways nanofluids may be used with nanoparticles, and they are dependent on the purposes. Nanoparticles with smaller diameters and higher densities are preferred over other particles in the EOR process. When the input fluid reached the reservoir at a high temperature, the thermal characteristics of the nanoparticles were also considered. According to the characteristics indicated, the three widely utilized nanoparticles in the procedure were SiO 2 , Al 2 O 3 , and CuO. As seen in Figure 12, these three nanoparticles have only a different impact on the recovery factor of EOR. However, the recovery rate of silicon nanoparticles is approximately 2% higher than that of the other There are several ways nanofluids may be used with nanoparticles, and they are dependent on the purposes. Nanoparticles with smaller diameters and higher densities are preferred over other particles in the EOR process. When the input fluid reached the reservoir at a high temperature, the thermal characteristics of the nanoparticles were also considered.
According to the characteristics indicated, the three widely utilized nanoparticles in the procedure were SiO 2 , Al 2 O 3 , and CuO. As seen in Figure 12, these three nanoparticles have only a different impact on the recovery factor of EOR. However, the recovery rate of silicon nanoparticles is approximately 2% higher than that of the other nanoparticles.

Conclusions
To simulate nanofluid and oil reservoir flooding, a two-dimensional rectangular porous homogenous area was filled with oil to simulate nanofluid and oil reservoir flooding. The effect of nanoparticle volume fraction on flooding and how it varied from water flooding was investigated. The EOR process was also simulated and compared with three types of nanoparticles: SiO 2 , Al 2 O 3 , and CuO. Studies were conducted to determine how the fluid intake temperatures affect the recovery factor, leading to the following conclusions: • Adding 5% silicon nanoparticles to a base fluid in a porous homogenous medium enhanced the oil recovery by more than 20%. • Even though the difference in the effects of the three nanoparticles described was less than 2%, silicon had a 2% better oil recovery rate than Al 2 O 3 and CuO. • Using two different temperatures to calculate the EOR recovery coefficient for heterogeneous porous geometry of 300 and 350 K revealed that temperature was By using nanoparticles in the fluid, the specific heat capacity of the fluid decreased, and the heat transfer coefficient of the fluid was exceptionally high. Due to the ease and speed of heat transfer input temperature, in this case, rapid changes and increased process performance were possible. The recovery coefficient was determined for two different fluid input temperatures in Figure 13. The amount of oil produced from the reservoir increased dramatically when the fluid intake temperature rose. Furthermore, because of the nanoparticles' strong thermal transfer characteristics, the heat was transferred to the oil inside the reservoir more quickly. As a result of the decreased viscosity and density of reservoir oil, the recovery factor was improved.

Conclusions
To simulate nanofluid and oil reservoir flooding, a two-dimensional rectangular porous homogenous area was filled with oil to simulate nanofluid and oil reservoir flooding. The effect of nanoparticle volume fraction on flooding and how it varied from water flooding was investigated. The EOR process was also simulated and compared with three types of nanoparticles: SiO 2 , Al 2 O 3 , and CuO. Studies were conducted to determine how the fluid intake temperatures affect the recovery factor, leading to the following conclusions: • Adding 5% silicon nanoparticles to a base fluid in a porous homogenous medium enhanced the oil recovery by more than 20%. • Even though the difference in the effects of the three nanoparticles described was less than 2%, silicon had a 2% better oil recovery rate than Al 2 O 3 and CuO. • Using two different temperatures to calculate the EOR recovery coefficient for heterogeneous porous geometry of 300 and 350 K revealed that temperature was crucial in the EOR process due to its influence on oil viscosity and density decrease. This impact was amplified when utilizing nanofluids for heat transmission due to the efficient and rapid heat transfer.

•
Raising the relative permeability of nanofluid while simultaneously reducing the relative permeability of oil due to the presence of nanoparticles was the primary reason for EOR.

•
The recovery factor of silicon nanoparticles is higher than the recovery factors of aluminum oxide and copper oxide nanoparticles. • Using the empirical correlations of the nanofluid thermophysical properties on the relative permeability equations of the nanofluid and oil in the proposed model, it is more accurate than other models to determine the higher recovery factor of one nanoparticle compared to other nanoparticles. Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.