Entropy Analysis and Melting Heat Transfer in the Carreau Thin Hybrid Nanoﬂuid Film Flow

: Melting heat transfer has a vital role in forming energy storage devices such as ﬂexible thin ﬁlm supercapacitors. This idea should be welcomed in the thin ﬁlm theoretical models to sustain technological advancement, which could later beneﬁt humankind. Hence, the present work endeavors to incorporate the melting heat transfer effect on the Carreau thin hybrid nanoﬂuid ﬁlm ﬂow over an unsteady accelerating sheet. The mathematical model that obeyed the boundary layer theory has been transformed into a solvable form via an apt similarity transformation. Furthermore, the collocation method, communicated through the MATLAB built-in bvp4c function, solved the model numerically. Non-uniqueness solutions have been identiﬁed, and solutions with negative ﬁlm thickness are unreliable. The melting heat transfer effect lowers the heat transfer rate without affecting the liquid ﬁlm thickness, while the Carreau hybrid nanoﬂuid contributes more entropy than the Carreau nanoﬂuid in the ﬂow regime.


Introduction
The thin film flow research has a great prospect in technological advancement due to its significance in producing electronic devices such as integrated circuits and microscopic fluidic devices [1]. For example, one can see how a solid surface is coated by a thin liquid film in those manufacturing operations. Another potential manufacturing process subset to the thin film flow application is the cast film extrusion that produces polymer sheets and films [2]. Realising the strength of the thin film flow research as an enzyme to attain the next stage of technological development, Wang [3] pioneered the problem of the thin film flow past an accelerating sheet and attested the unavailability of similarity solutions when the flow unsteadiness' rate exceeds 2. Then, Usha and Sridharan [4] revisited the flow problem in [3] asymmetrically and proved that the similarity solutions are absent when the flow unsteadiness' rate exceeds 4. As time went by, the researchers learned that heat transfer analysis is crucial in the thin film flow problem, aligned with the initiative to comprehend the heat exchangers and chemical processing equipment's design. Thus, Andersson et al. [5] solved the thin film flow and heat transfer past an accelerating sheet; they also devised a novel similarity solution for the temperature field. This contribution of [5] is remarkable and highly assists in investigating the heat transfer aspect in the present work. After that, Wang [6] presented the analytic solutions for the thin film flow and heat transfer problem over an unsteady accelerating sheet. The strong contributions of [3][4][5][6] are the impetus for the thin film flow and heat transfer research under the following effects: thermocapilarity [7,8], general surface temperature [9,10], thermal radiation [11], magnetohydrodynamics (MHD) [12], viscous dissipation [13], slip effects [14,15].
The strength of the non-Newtonian fluid in illuminating varying fluid viscosity under the applied force has vast industrial applications and managed to attract the researchers' attention to be considered under various settings, such as in [16,17]. On the other hand, the researchers' consideration of the non-Newtonian fluid in the thin film flow problem is raised because the protective coating applied on an extrudate is a non-Newtonian fluid. Therefore, Andersson et al. [18] investigated the power-law thin liquid film flow past an accelerating sheet and found a contradict trend in the fluid velocity when the power-law fluid adapts to the pseudoplastic and dilatant features respectively. Furthermore, Chen [19] enhanced the work of [18] by incorporating the heat transfer characteristic as it is an essential factor to decide the final product's quality. Subsequently, the researchers' consideration of other generalised non-Newtonian fluid models such as the Carreau fluid model increased due to its validity for high and low shear rates. Myers [20] critically analysed the generalised non-Newtonian fluid's potential in the thin film flow and suggested that the Carreau fluid model is a better choice due to its accuracy rate. Accordingly, there is a number of significant works reported within the scope of the Carreau thin film flow; see [21][22][23][24].
Besides that, hybrid nanofluid is an incredible invention by humankind to uplift technological advancement to the next level. Choi and Eastman [25] introduced the brilliant idea of suspending the nanosized metal element in the fluid to boost its heat transfer rate. Although nanofluid hits the peak of the researchers' interest due to its applications in the heat transfer equipment, nanofluid is incompatible with some specified real-world applications that require substitution between some nanofluids' properties [26]. Thus, the hybrid nanofluid is proposed to encounter this issue through the experimental works [27,28], and the hybrid nanofluid managed to gain vast interest from the researchers. Shortly, the theoretical works, such as [29] gained momentum in the hybrid nanofluid after the valuable works of Devi and Devi [30,31]. The hybrid nanofluid also succeeded in the thin film flow over an accelerating sheet owing to industrial applications such as microfluidics [32]. For instance, Sadiq et al. [33] explored the Maxwell thin hybrid nanofluid film flow across the Darcy-Forchheimer porous media and inferred that the augmentation of the heat transfer rate in the hybrid Maxwell nanofluid is better than the single-typed Maxwell nanofluid.
Melting heat transfer is a common phase experienced by industrial processes such as casting. Epstein and Cho [34] is one of the earliest works involved in the melting heat transfer effect in the laminar boundary layer flow. Epstein and Cho [34] solved the boundary layer flow problem past the horizontal positioned static flat surface along with the melting heat transfer effect. Ishak et al. [35] extended the work in [29] by considering the flow past the moving sheet and reported that the melting heat transfer effect is the decreasing function of the convective heat transfer rate. Then, Khashi'ie et al. [36] reconsidered the problem solved in [35] by incorporating the presence of hybrid nanoparticles and corroborated the finding as mentioned above in [35]. Even though the melting heat transfer effect has been probed under several settings; see [37,38]; yet, the inspection of the melting heat transfer effect in the thin film flow is scarce. Thus, this motivates the present work to scruntise the melting heat transfer effect in the thin film flow.
Overall, the present work attempts to solve the problem of the Carreau thin hybrid nanofluid film flow past an accelerating sheet under the influence of the melting heat transfer. The present model is relatively new, and entropy generation analysis is performed. The present work also adapts the similarity transformation suggested by Andersson et al. [5], and the formulated mathematical model has been solved in the built-in collocation method, the bvp4c, to produce the approximate solutions. Furthermore, non-uniqueness solutions have been reported for every case of governing parameters' variations. The findings of the present work may serve as a reference for improving the material processing industry.

Mathematical Model
Ruminate the Carreau fluid flow bounded by a thin liquid film and a horizontally placed accelerating sheet from a narrow opening at the Cartesian coordinate system origin. The two-dimensional flow is assumed to be incompressible, unsteady while the thin liquid film has an unvarying thickness, H(t). Figure 1 depicts the flow setup, and y−coordinate is located normal to the x−coordinate. The sheet's accelerated act, which portrays the stretching sheet situation, brings about the fluid motion delimited by the thin film and the accelerating sheet. The sheet is accelerated with speed U w (x, t) = bx (1−σt) , where b and σ are positive constants with dimension time -1 , σt = 1, while b > 0 conveys the stretching rate. The sheet surface is impermeable and melts. The melting surface temperature is denoted by T m , whereas T is the fluid temperature. The wall temperature, T w is defined [39]. Here, the slit temperature and reference temperature are denoted by T s and T 0 , respectively. Besides, the end effects and gravity are assumed to be very small and thus omitted. The formulated boundary layer model in the present work is only sensible if the liquid film thickness does not overlap with the boundary layer thickness. Otherwise, the present formulated model becomes irrational [40]. Also, the planar thin liquid film is assumed to be smooth and free of any surface waves [5].
uniqueness solutions have been reported for every case of governing parameters' variations. The findings of the present work may serve as a reference for improving the material processing industry.

Mathematical Model
Ruminate the Carreau fluid flow bounded by a thin liquid film and a horizontally placed accelerating sheet from a narrow opening at the Cartesian coordinate system origin. The two-dimensional flow is assumed to be incompressible, unsteady while the thin liquid film has an unvarying thickness, ( ).
H t Figure 1 depicts the flow setup, and y − coordinate is located normal to the x − coordinate. The sheet's accelerated act, which portrays the stretching sheet situation, brings about the fluid motion delimited by the thin film and the accelerating sheet. The sheet is accelerated with speed where b and σ are positive constants with dimension time -1 , conveys the stretching rate. The sheet surface is impermeable and melts. The melting surface temperature is denoted by  [39]. Here, the slit temperature and reference temperature are denoted by s T and 0 , T respectively. Besides, the end effects and gravity are assumed to be very small and thus omitted. The formulated boundary layer model in the present work is only sensible if the liquid film thickness does not overlap with the boundary layer thickness. Otherwise, the present formulated model becomes irrational [40]. Also, the planar thin liquid film is assumed to be smooth and free of any surface waves [5]. The Carreau fluid's Cauchy stress tensor is given as [20] , p τ η where ( ) ( ) The Carreau fluid's Cauchy stress tensor is given as [20] where Here, τ is the Cauchy stress tensor, p is the pressure, I denotes the identity tensor, η 0 signifies the zero-shear-rate viscosity, η ∞ is the infinite-shear-rate viscosity, λ implies the material time constant, and n represents the power-law index. The shear rate, . γ can be elaborated as In Equation (3), Π is the second invariant strain rate tensor and A 1 is the Rivlin-Ericksen tensor expressed further as The most practical cases where η 0 η ∞ is considered. Normally, the value of η ∞ is determined by the extrapolation procedure or chosen to be zero (suggested theoretical value) [41]. Thus, in the present work, the value of η ∞ is set to zero, and affect Equation (1) to become The Carreau fluid model shows pseudoplastic, dilatant, and Newtonian features when 0 < n < 1, n > 1 and n = 0, respectively, where n is the power-law index. Under these assumptions, the governing liquid film flow of the Carreau fluid can be written as [42].

∂u ∂x
∂u ∂T ∂t where u and v are the velocity components along the xand ydirections, respectively, λ is a material time constant, n signifies the power-law index. Meanwhile, ρC p hn f , µ hn f , ρ hn f , and k hn f are the hybrid nanofluid's heat capacity, dynamic viscosity, density and thermal conductivity, respectively. The further definition of µ hn f , ρ hn f , k hn f , and ρC p hn f are expressed in Table 1. Table 1. The hybrid nanofluid's correlation properties' definitions (see [43]).

Au-Cu/(CMC/H 2 O) Mathematical Relation
Density Thermal capacity According to Table 1, the nanoparticle volume fraction is φ and φ = 0 reduces the model into a regular fluid. Next, φ 1 and φ 2 signify the Au's and Cu's nanoparticle volume fraction, respectively. The total volume concentration of two types of nanoparticles suspended in the hybrid nanofluid is determined as φ hn f = φ 1 + φ 2 . Meanwhile, ρ f and ρ hn f are the densities of the base fluid and the hybrid nanoparticle, respectively, k f and k hn f are the thermal conductivities of the base fluid and the hybrid nanoparticles, respectively, ρC p f and ρC p hn f are the heat capacitance of the base fluid and the hybrid nanoparticle, respectively. These correlations are based on physical assumptions and agree with the conservation of mass and energy. Thus, the physical properties of the base fluid (water), gold (Au) and copper (Cu) hybrid nanofluids are given in Table 2. Table 2. The thermophysical properties of selected nanoparticles and base fluid (sodium carboxymethyl cellulose (CMC)/water) (see [44,45] The Equations (6)-(8) are getting along with the boundary conditions At y = h, the kinematic constraint is enforced in the fluid motion through v = dh/dt. The wall shear stress and heat flux disappear entirely at the adiabatic free surface and thus ∂u/∂y = ∂T/∂y = 0 at y = h. Next, we introduce the similarity transformations as follows (Andersson et al. [5]): where prime infers the derivative concerning ζ. Employing the similarity conversion as in (10) and (11) into the governing model (6)-(9) satisfies the continuity equation, and the remaining equations are transformed as follows: with the boundary conditions wherein We = 3 , is the local Weissenberg number [46], ω = σ b is the dimensionless measure of unsteadiness, the Prandtl number is defined as the melting heat transfer parameter is signified by . Moreover, β is an unknown constant that conveying the dimensionless film thickness. β also implies the similarity variable (ζ) value at the free surface, and hence the expression in (11) can take the following form: This unknown constant β must be calculated as an integral part of the boundary-value problem. Thus, elucidates the film thickness's rate of change. On the other hand, when n = 1 and We = 0, the Carreau fluid model in Equations (12) and (13) reveals the Newtonian characteristics.
The physical quantities of interest in the present work are the local skin friction coefficient C f x and the local Nusselt number (Nu x ), which can be defined as follows: Here, the wall shear stress (τ w ) and the heat flux from the surface of the sheet (q w ) are given by [47].
The local Reynolds number is defined as Re

Entropy Analysis
The entropy generation analysis is ideal for calculating the dissipated heat energy and measuring any flow systems' performance deterioration. This nonconserved property can be communicated in the following dimensional form [48]: where terms on heat irreversibility is given followed by the fluid friction irreversibility. By using Equations (10) and (11) into (20), the nondimensional form of the volumetric entropy generation can be formed as follows: Here, the local entropy generation rate N s , the Brinkman number (Br), and the temperature difference parameter (ε) are further defined as follows: Mathematics 2021, 9, 3092 7 of 19 In the entropy generation analysis, it is necessary to compute the Bejan number to identify which entropy generation appears to dominate the flow system, either frictional or thermal entropy generation. Thus, the Bejan number is calculated in the present work, and it takes the following form: is the irreversibility ratio. Expression in (23) signifies the heat transfer irreversibility and total entropy ratio in the flow system. Contact melting, lubrication, and electronic cooling are some heat transfer applications where the Bejan number is vital. Equation (23) also can be utilised to determine the following effects in the flow system (see Table 3): Table 3. Physical significance of the Bejan number [49].

2.
When Φ = 1, Be = 0.5 Both heat transfer irreversibility and fluid friction irreversibility are equally influencing the flow system.

3.
When Be = 0 Fluid friction irreversibility is highly influencing the flow system.

Results and Discussion
The transport phenomena in the thin film flow regime can be learned by plotting the velocity and temperature profiles. In addition, calculating the local skin friction coefficient, the local Nusselt number and the dimensionless film thickness when the pertinent parameters vary is also necessary to inspect the present model's performance. Thus, all numerical outputs were generated by setting the governing parameters' values within the following fixed range: 0 ≤ We ≤ 0.3, 0 ≤ χ ≤ 2.5, 0.8 ≤ ω ≤ 1.4, and 0.6 ≤ n ≤ 1.6. The Prandtl number is fixed to 8 throughout the computation process. Also, φ 1 and φ 2 represents the gold (Au) and copper (Cu) nanoparticle volume fractions, respectively. The numerical outputs are compared between the Carreau hybrid nanofluid case, where φ 1 = 0.02, φ 2 = 0.03, and the single-typed Carreau nanofluid is considered with φ 1 = 0, φ 2 = 0.03. These parameter values are chosen based on the availability of the numerical solutions. However, those parameter values lie within the acceptable range established in previously published works. Equations (12)- (14), which convey the simplified form of the present thin film flow problem, are solved using the bvp4c function found in the MATLAB 2019a software. This built-in collocation code eases the solving process even though the present work has dimensionless film thickness as the unknown parameter [50]. Besides that, all computed numerical results are accurate within 1 × 10 −10 . In order to test the precision of the present method, the thin film flow problem studied by Wang [6] have been resolved via the bvp4c function, and the comparison of the results is given in Table 4. Table 4 proves that the built-in collocation method agrees well with the numerical results produced via the homotopy analysis method in [6]. Meanwhile, the CPU time for calculating the non-uniqueness solutions is presented in Table 5. In this sample, it is apparent that the CPU time increases from the first to the second solution in every case of n. Before the presentation and discussion of the results go further, it is appropriate to confer about the non-uniqueness numerical solutions. It is undeniable that more than one numerical solutions are obtainable by providing a good set of guess values since that is the built-in bvp4c routine's requirement. The present work found that the second solutions always yield negative film thickness. The negative film thickness implies the thin liquid film's distortion, and hence an ideal thin liquid film cannot be formed [23]. Therefore, the trends showed by the second solution are disregarded. Now, Table 6 demonstrates the trend of β for the Carreau hybrid nanofluid and Carreau nanofluid when We increases. The increment of We from 0 to 0.3 affects the dimensionless film thickness to decrease by 0.25% for Carreau hybrid nanofluid and decrement by 0.18% for the Carreau nanofluid. The reason for this occurrence can be collected from Figure 2 and Table 7. Figure 2 shows that the fluid velocity increases insignificantly across the flow regime when We increases. The increment in We elucidates a longer relaxation rate; thus, the Carreau fluid takes more time to react with the external forces. Since the accelerating sheet imposes drag force towards the Carreau fluid (this can be evident by the negative values of C f x Re 1/2 x in Table 7), the fluid velocity increases slightly past the unsteady accelerating sheet, which elevates the wall shear stress and increases the values of C f x Re 1/2 x in a minimal amount.  Meanwhile, Figure 3 displays the temperature profiles across the thin film flow vicinity. For the Carreau fluid associated with the hybrid nanoparticles, the temperature increases when We increases. This behavior is evident when .
. . ζ = 0 72727 Such an interesting difference in the trend might be due to the suspended nanoparticles' thermal conductivity in the base fluid. Besides that, Table 8 tabulates changes in the heat transfer rate at the accelerating impermeable surface. The flow with the copper nanoparticles spectacle a gradual increase in ( ), θ ′ 0 and this is acceptable because copper has better thermal conductivity and results in an increased rate in heat exchange. However, the hybrid nanofluid does not exhibit a gradual degree of improvement. The heat transfer rate decreases when We 's value increases from 0 to 0.3. The increment in the fluid relaxation time affects the Carreau hybrid nanofluid to become warm, which lowers heat flux from the accelerating surface and diminishes ( ). θ ′ 0  Table 7. Numerical outputs of C f x Re 1/2 x for the hybrid nanofluid (φ 1 = 0.02, φ 2 = 0.03) and nanofluid when ω = 0.8, n = 0.8, Pr = 8, and χ = 1.5. Meanwhile, Figure 3 displays the temperature profiles across the thin film flow vicinity. For the Carreau fluid associated with the hybrid nanoparticles, the temperature increases when We increases. This behavior is evident when 0 ≤ ζ < 0.72727. It should be noted that the fluid is under the influence of shear thinning effect (n = 0.8), and hence the fluid viscosity may decrease with the act of the accelerating sheet and the melting heat transfer's effect. At this moment, when the effect of We is amplified, the fluid temperature augments since the relaxation time is prolonged. A similar result has been reported by Hayat et al. [51]. However, after ζ = 0.72727, the fluid far from the accelerating surface is less affected with the melting heat transfer and thus, the fluid temperature decreases when We increases. On the other hand, in the case of the Carreau nanofluid, fluid temperature increases with declining values of We when 0 ≤ ζ < 0.72727, and the opposite trend is observed after ζ = 0.72727. Such an interesting difference in the trend might be due to the suspended nanoparticles' thermal conductivity in the base fluid. Besides that, Table 8 tabulates changes in the heat transfer rate at the accelerating impermeable surface. The flow with the copper nanoparticles spectacle a gradual increase in θ (0), and this is acceptable because copper has better thermal conductivity and results in an increased rate in heat exchange. However, the hybrid nanofluid does not exhibit a gradual degree of improvement. The heat transfer rate decreases when We's value increases from 0 to 0.3. The increment in the fluid relaxation time affects the Carreau hybrid nanofluid to become warm, which lowers heat flux from the accelerating surface and diminishes θ (0).  Table 7. Numerical outputs of 1 2 Re fx x C for the hybrid nanofluid ( ) The results' discussion is further by examining the behavior of the temperature profiles and ( ) θ ′ 0 when χ increases. Even though the temperature profiles in Figure 4 reveal an unusual degree of dissonance, Table 9 informs that the convective heat transfer rate at the accelerating sheet increases in the single-typed nanofluid but deteriorates in the Carreau fluid with hybrid suspensions. The increment in χ indicates more cold fluid molecules exist from the melting accelerating sheet towards the warm fluid. Therefore, the nanofluid temperature declines at the moving surface and augments the heat exchange rate. The work of Khan et al. [52] also conveyed such similar result. Conversely, the hybrid nanofluid temperature may have retained a low thermal conductivity, ensuring the low heat flux rate and reducing ( ) θ ′ 0 at the accelerating surface. Table 10 con-  Table 8. Numerical outputs of θ (0) for the hybrid nanofluid (φ 1 = 0.02, φ 2 = 0.03) and nanofluid when ω = 0.8, n = 0.8, Pr = 8, and χ = 1.5. The results' discussion is further by examining the behavior of the temperature profiles and θ (0) when χ increases. Even though the temperature profiles in Figure 4 reveal an unusual degree of dissonance, Table 9 informs that the convective heat transfer rate at the accelerating sheet increases in the single-typed nanofluid but deteriorates in the Carreau fluid with hybrid suspensions. The increment in χ indicates more cold fluid molecules exist from the melting accelerating sheet towards the warm fluid. Therefore, the nanofluid temperature declines at the moving surface and augments the heat exchange rate. The work of Khan et al. [52] also conveyed such similar result. Conversely, the hybrid nanofluid temperature may have retained a low thermal conductivity, ensuring the low heat flux rate and reducing θ (0) at the accelerating surface. Table 10 confirms that when the Carreau fluid with the presence of the nanoparticles changes its character from the shear thinning to the shear thickening feature, the dimensionless film thickness slightly increases, in increments of about 0.035% and 0.025% for the Carreau hybrid nanofluid and Carreau single-typed nanofluid, respectively.

We
Mathematics 2021, 9,3092 12 of 21 acter from the shear thinning to the shear thickening feature, the dimensionless film thickness slightly increases, in increments of about 0.035% and 0.025% for the Carreau hybrid nanofluid and Carreau single-typed nanofluid, respectively.    Table 9. Numerical outputs of ( )   Meanwhile, the temperature profiles in Figure 5 display that when ζ ≤ 0.444, the Carreau hybrid nanofluid temperature decreases while n increases from 0.6 to 1.6. This is because the dilatant feature retards the heat energy transmission. However, when the liquid film vicinity travel at ζ > 0.444, the fluid temperature becomes an increasing function of n. The area far from the accelerating sheet is possibly less affected by the shear force from the accelerating sheet, so heat transmission is reduced. Also, from Figure 5, it is observed that the fluid temperature is low at the area far from the accelerating sheet compared to the fluid temperature near the sheet's surface. Table 11 identifies that the heat transfer rate at the accelerating sheet gradually decreases for the Carreau hybrid nanofluid when n increases. This is true because Metzner et al. [53] corroborated that dilatant fluid has a lower heat transfer rate than the shear-thinning fluid. On the other hand, the opposite trend is perceived for the single-typed hybrid nanofluid. In the single-typed nanofluid considered in the present work, the nanoparticle volume fraction is less than the hybrid nanofluid, which may increase the fluid's thermal conductivity as the low-temperature molecules enter the flow stream. Thus, a moderate increase in θ (0)'s values can be noticed along with the increment of n. Table 12 delivers the decrement of C f x Re 1/2 x along with the increment of n. The strengthening effect of n reduces the wall shear stress at the accelerating sheet. Hence, the values of C f x Re 1/2 x decline. Meanwhile, the temperature profiles in Figure 5 display that when . , ζ ≤ 0 444 the Carreau hybrid nanofluid temperature decreases while n increases from 0.6 to 1.6. This is because the dilatant feature retards the heat energy transmission. However, when the liquid film vicinity travel at . , ζ > 0 444 the fluid temperature becomes an increasing function of .
n The area far from the accelerating sheet is possibly less affected by the shear force from the accelerating sheet, so heat transmission is reduced. Also, from Figure 5, it is observed that the fluid temperature is low at the area far from the accelerating sheet compared to the fluid temperature near the sheet's surface. Table 11 identifies that the heat transfer rate at the accelerating sheet gradually decreases for the Carreau hybrid nanofluid when n increases. This is true because Metzner et al. [53] corroborated that dilatant fluid has a lower heat transfer rate than the shear-thinning fluid. On the other hand, the opposite trend is perceived for the single-typed hybrid nanofluid. In the singletyped nanofluid considered in the present work, the nanoparticle volume fraction is less than the hybrid nanofluid, which may increase the fluid's thermal conductivity as the lowtemperature molecules enter the flow stream. Thus, a moderate increase in ( ) θ ′ 0 's values can be noticed along with the increment of .
n Table 12 delivers the decrement of 1 2 Re fx x C along with the increment of . n The strengthening effect of n reduces the wall shear stress at the accelerating sheet. Hence, the values of 1 2 Re fx x C decline.       Figure 8, the state of the Carreau fluid from portraying the shear-thinning trait and then to the shear-thickening feature yields more heat energy incorporated to the system but minimal. However, the heat transfer irreversibility highly influences the fluid flow system, although Figure 9 shows the minor decrement in Be when the values of n increases.
from Figure 8, the state of the Carreau fluid from portraying the shear-thi then to the shear-thickening feature yields more heat energy incorporated but minimal. However, the heat transfer irreversibility highly influences system, although Figure 9 shows the minor decrement in Be when the v creases.

Conclusions
The present work endeavored to investigate the performance of the Carreau thin hybrid nanofluid film flow and heat transfer while the melting heat transfer effect imposed on the accelerating sheet. Interestingly, this is the original work in the thin film Carreau hybrid nanofluid theoretical model, considering the impact of the melting heat transfer. Two approximate solutions were identified for every variation case. The Carreau hybrid nanofluid's heat transfer rate decreases when the fluid adapts to the dilatant feature. Besides that, the presence of the hybrid nanoparticles promotes entropy in the flow system compared to the mono-typed nanoparticles. The numerical solutions with negative film thickness indicated defective thin film flow and unreliable. However, the melting heat

Conclusions
The present work endeavored to investigate the performance of the Carreau thin hybrid nanofluid film flow and heat transfer while the melting heat transfer effect imposed on the accelerating sheet. Interestingly, this is the original work in the thin film Carreau hybrid nanofluid theoretical model, considering the impact of the melting heat transfer. Two approximate solutions were identified for every variation case. The Carreau hybrid nanofluid's heat transfer rate decreases when the fluid adapts to the dilatant feature.
Besides that, the presence of the hybrid nanoparticles promotes entropy in the flow system compared to the mono-typed nanoparticles. The numerical solutions with negative film thickness indicated defective thin film flow and unreliable. However, the melting heat transfer effect does not affect the liquid film thickness.