Numerical Simulation Of Heat Transfer Research Articles

Numerical study of blood-based MHD tangent hyperbolic hybrid nanofluid flow over a permeable stretching sheet with variable thermal conductivity and cross-diffusion

Abstract Modern heat transport processes such as fuel cells, hybrid engines, microelectronics, refrigerators, heat exchangers, grinding, coolers, machining, and pharmaceutical operations may benefit from the unique properties of nanoliquids. By considering Al 2 O 3 {{\rm{Al}}}_{2}{{\rm{O}}}_{3} nanoparticles as a solo model and Al 2 O 3 – Cu {{\rm{Al}}}_{2}{{\rm{O}}}_{3}{\rm{\mbox{--}}}{\rm{Cu}} as hybrid nanocomposites in a hyperbolic tangent fluid, numerical simulations for heat and mass transfer have been established. To compare the thermal acts of the nanofluid and hybrid nanofluid, bvp4c computes the solution for the created mathematical equations with the help of MATLAB software. The impacts of thermal radiation, such as altering thermal conductivity and cross-diffusion, as well as flow and thermal facts, including a stretchy surface with hydromagnetic, and Joule heating, were also included. Furthermore, the hybrid nanofluid generates heat faster than a nanofluid. The temperature and concentration profiles increase with the Dufour and the Soret numbers, respectively. The upsurge permeability and Weissenberg parameter decline to the velocity. An upsurge variable of the thermal conductivity grows to the temperature profile. Compared to the nanofluids, the hybrid nanofluids have higher thermal efficiency, making them a more effective working fluid. The magnetic field strength significantly reduces the movement and has a striking effect on the width of the momentum boundary layer.

Numerical simulation of heat transfer and nanoparticle transport in a nanofluid within a tumour surrounding a blood vessel

A Newtonian nanofluid containing suspended nanoparticles can substantially improve heat transfer due to enhanced energy transport mechanisms. This theoretical study investigates heat and mass transfer in biological tissues using such a nanofluid under a magnetic field. These properties have promising medical and engineering applications. The nonlinear governing equations were transformed into ordinary differential equations using similarity variables and numerically solved with MATLAB boundary value problem solver bvp4c, subject to appropriate boundary conditions. Results demonstrated increasing the heat source parameter dramatically raised tumor interstitial temperature. This heating, along with improved nanoparticle accumulation within the tumor due to the thermal effects, are together essential for effective hyperthermia treatment. The model provides new insights into tuning heat and mass transport mechanisms in biological tissues via nanofluids for therapeutic applications. Therefore, the findings of this study may improve the efficacy of thermal therapy in treating cancer.

Liquid-liquid Taylor flow heat transfer: Correlations for predicting Nusselt number and pressure drop based on four independent variables

The numerical simulation of heat transfer in Taylor flows involving water-light mineral oil, water-hexadecane, and water-dodecane configurations was conducted, considering water as the primary phase in three scenarios and oils as the primary phases in three others. Additionally, an investigation into heat transfer in Water-AR20 silicon oil with water as the primary phase was undertaken. In each Taylor flow scenario, two correlations were presented for predicting the thermally developed Nusselt number and pressure drop per unit length, obviating the necessity for experimental observation or numerical simulation. In cases where water served as the secondary phase, fully developed Nusselt number of the two-phase flow increased as the dimensionless effective diameter of the water droplets increased. Anomalously, when water was the primary phase, an optimal dimensionless effective diameter, equivalent to 1 for the secondary phase droplets, resulted in maximum Nusselt. The findings indicated that utilizing water as the primary phase leads to superior heat transfer capacity compared to employing oils as the primary phase. Furthermore, it was deduced that, under the condition where the effective diameter of the secondary oil droplets is set at 1 and three other variables are the same, alterations in oil type yields negligible variations in the Nusselt number. The findings of this paper can be a guide in applying liquid–liquid Taylor flow as a coolant in microchannel heatsinks.

Numerical simulation of combined convection and radiation heat transfer in a hybrid nanofluid inside an open fins cavity under a magnetic field

Numerical simulation of combined convection and radiation heat transfer in a hybrid nanofluid inside an open fins cavity under a magnetic field

Numerical modeling of mixed convective nanofluid flow with fractal stochastic heat and mass transfer using finite differences

This study presents the first comprehensive numerical simulation of heat and mass transfer in fractal-like mixed convective nanofluid flows. The flow of non-Newtonian nanofluids over flat and oscillating sheets is modelled mathematically, and a finite difference scheme is used to solve this model. The two-stage scheme can tackle fractal and fractal stochastic mathematical models of partial differential equations. The consistency in the mean square is proved, and Fourier series stability analysis is adopted to find stability conditions for fractal stochastic partial differential equation. The scheme is applied to solve the unsteady Casson nanofluid flow over the flat and oscillatory sheet, which affects thermal radiation, heat source, and chemical reaction. The existence of the solution is also provided for the Navier-Stokes equation of the considered flow model using fractal time derivative. The graph illustrates that the proposed fractal technique achieves faster convergence than the Crank-Nicolson approach. Applications in energy systems, materials science, and environmental engineering are just a few of the domains that could benefit from a better understanding of mixed convective nanofluid flows with fractal features, and that is what this research study hopes to accomplish. Scientists and engineers may better develop efficient and environmentally friendly systems by simulating and analyzing these complicated processes with the suggested finite difference technique.

3D numerical simulations of mixed convective heat transfer and correlation development for a thermal manikin head

The head represents 10% of the body's total surface area. Unprotected, it accounts for a significant portion of overall heat loss when exposed to cold conditions. This study was motivated by a need to clarify how the human head interacts with its environment in terms of heat exchange. Accurate estimations of heat transfer coefficients on the human head are essential for conducting thermal comfort and safety analyses in buildings. In this study, a thermal head resembling a real male human head is utilized to investigate heat transfer between the body and the surrounding environment. A three-dimensional computational fluid dynamics (CFD) model is proposed to simulate steady-state dry heat loss from the human head within a chamber. This model provides predictions for heat flux, temperature, and velocity distribution surrounding the head. A straightforward correlation, derived from numerical and experimental findings, is introduced to forecast the average Nusselt number for the head under combined natural and forced convection. This correlation, relying on dimensionless parameters (Grashof, Reynolds, and Prandtl numbers), offers enhanced accuracy, simplicity, and fewer terms. The predicted average Nusselt numbers from the proposed correlation for mixed convection closely match CFD and experimental results, with relative percentage differences within ±2%, signifying excellent accuracy across a broader range of flow conditions, including temperature differences and air velocities. Additionally, the study explores the impact of head diameter on overall heat transfer.

Numerical simulation of mass and heat transfer for water extraction from icy lunar regolith

The confirmation of water ice’s existence in the permanent shadow area led to extensive research in the field of water ice mining from icy lunar regolith. A priori numerical simulation of water ice mining is necessary for guiding the development of lunar water ice mining schemes more reasonably. A 2D axisymmetric numerical simulation model capable of simulating the thermal extraction process about mining water ice from icy lunar regolith is constructed, which is executed in the COMSOL Multiphysics. The thermal extraction cases of lunar regolith with different initial water ice content and heating fluxes are simulated. The EER (energy efficiency ratio) is used to evaluate the efficiency of thermal extraction. The results show that the EER is higher as the initial water ice content is increased, which means more power is used to heat water and less power is used to heat the regolith. The icy lunar regolith with initial water ice content higher than 5.0 wt% is found to be more valuable, over which the EER at the end of thermal extraction will not increase much as the initial water ice content increases. However, the higher heating flux leads to the lower EER at the end of thermal extraction. The speed and economics of thermal extraction are suggested to be weighted before the mission’s implementation. The status of thermal diffusion (thermal transpiration) is studied, and the results indicate that thermal diffusion and advection both can be ignored in thermal extraction modeling, unless the average magnitude of temperature gradient and pressure gradient exceed the maximum of 75764 K/m and 11053 Pa/m in our calculation cases, respectively.

Numerical Simulation of Transient Flow and Heat Transfer in a Steel Ladle During Holding Period

Numerical Simulation of Transient Flow and Heat Transfer in a Steel Ladle During Holding Period

Numerical simulation of heat and mass transfer in laser directed energy deposition

In order to accurately predict the heat transfer and flow process of laser directed energy deposition, the numerical simulation of whole stage after the powder leaves the nozzle was carried out. First, the powder flow mode with laser body heat source in the first stage was established, and the average temperature and velocity of the powder at the convergence plane were obtained as 1300 K and 13 m/s, respectively. Based on this, the coupling model of heat flow between powder and substrate in the second stage is improved, and the influence of multiple source terms such as powder temperature, powder velocity, gas velocity, surface tension and thermal buoyancy is considered comprehensively. The evolution of the temperature and flow field of the cladding layer are analyzed in detail under the spatiotemporal variation, and it is clear that there are two convection currents with different clockwise directions in the molten pool, and the Marangoni convection occupies the main position. Finally, FeCrNiMo alloy powder was deposited on 42CrMo substrate, and its morphology, hardness and tensile strength in different directions were analyzed. The average height, width and depth deviations of the experiment and simulation are 17.3 %, 15.5 % and 22.9 %, respectively. The average strength of the tensile sample V2 prepared by the cladding layer is 11.8 % higher than that of the substrate.

Study on the turbulent Prandtl number model for liquid metal flow and heat transfer in a narrow rectangular channel

In order to improve the core power density of liquid metal reactors, plate-type fuel assemblies are typically used, which have obvious advantages in heat transfer capacity. For the low Prandtl number fluid, the RANS turbulence model-based CFD method needs to be modified by inserting the turbulent Prandtl number model. This paper investigated the turbulent Prandtl number model for the turbulent heat transfer of liquid metal flow in a narrow rectangular channel. First, the existing turbulent Prandtl number models were summarized and classified as the local turbulence Prandtl number model and the global turbulence Prandtl number model. They were evaluated based on the experimental data. The results showed that the prediction results of each model are not ideal. Then, a typical numerical example was used to analyze the temperature field and flow field in the narrow gap of a narrow rectangular channel. It was shown that the entire flow region only includes the near-wall region and the transition region, and the transition region accounts for the main proportion. A semi-empirical local turbulence Prandtl number model was obtained by solving the turbulent transport equation, and the unknown coefficients were determined iteratively with the experimental data. The error was controlled within ±20%. Finally, the applicability of the newly proposed turbulent Prandtl number model was explored, and the results showed that it has good applicability for the different runner sizes and different liquid metal types. This study can provide a reference for the numerical simulation of turbulent heat transfer of metal liquid in the narrow rectangular channel.

Flow and heat transfer on the pressure side of C3X vane with various cooling hole blockages

Performance of film cooling under blockage has been extensively studied on flat plates. However, the ideal film cooling patterns obtained from flat plates may not be directly applicable to realistic engineering applications involving curved vanes or surfaces. The pressure gradient on a curved surface can significantly alter the behavior and trajectory of the coolant jet. This can result in variations in film cooling effectiveness and heat transfer characteristics compared to flat plate configurations. Therefore, in this study, a conjugate heat transfer numerical simulation was performed on the C3X vane under three sorts of cooling hole blockages including all-round blockage (RB), leading-edge blockage (LB), and trailing edge blockage (TB) with blockage ratio of 0.2, 0.35 and 0.5, respectively. Results show that RB changes the cooling effectiveness in the mid-span section nonlinearly, the cooling effectiveness was decreased by (0.019, 0.023, 0.038). LB is found beneficial to local film cooling downstream the cooling hole near the mid-span and tip area of a vane at blockage ratio of 0.2 with an increment of 0.021 in the mid-span section and 0.013 in 86% span section. As the blockage ratio increases to 0.35 and 0.5, LB produces adverse influences and the mid-span cooling effectiveness was decreased by 0.013 and 0.043, respectively. TB riseed the cooling hole outlet temperature marginally. The cooling effectiveness near the hole outlet in the mid-span section was decreased by (0.020, 0.034, 0.024). However, under blockage ratio of 0.2 and 0.35, the coolant jet is found to reattach the vane surface, resulting in improved film cooling effectiveness.

Enhanced parallel computation for time-fractional fluid dynamics: A fast time-stepping method with Newton-Krylov-Schwarz solver

This paper presents a sum-of-exponentials domain decomposition method for the numerical simulation of two-dimensional unsteady fluid flow and heat transfer using a time-fractional fluid model. We employ a fast time-stepping approach to discretize the time-fractional derivatives, followed by the application of a parallel Newton-Krylov-Schwarz algorithm to solve the resulting discrete nonlinear system. The numerical experiments demonstrate a superior performance of the fast time-stepping method compared to the traditional L1 scheme, particularly in the context of parallel computation, as it substantially reduces the computational complexity and memory usage. The algorithm exhibits robust performance across a wide range of model parameters and solver settings. Notably, it achieves 60% parallel efficiency with the use of 768 processor cores. This study underscores the efficacy of parallel processing as a potent computational strategy for addressing the challenges in solving time-fractional partial differential equations.

Effects of arcade design on wind and thermal environment inside an idealized 2D urban canyon with realistic solar heating

The arcade design has been widely used in the architecture design since it can provide space for commercial retails and improve local thermal comfort. However, the influences of different arcade design parameters on wind and thermal environment under realistic solar heating conditions are not fully understood, which hinders further application of the arcade design. This paper aims to investigate the effects of different arcade design parameters on wind and thermal environment under different solar heating conditions. Numerical wind flow, solar radiation and heat transfer simulations of a two-dimensional deep street canyon with aspect ratio of three are performed. The steady Reynolds-Averaged Navier-Stokes model with standard k−ε turbulence model is employed in this investigation, while considering solar radiation model. Moreover, the realistic solar heating simulations are achieved by setting different local solar times, i.e., LST0900 (leeward wall heating), LST1200 (street ground heating) and LST1500, (windward wall heating). Under realistic solar heating conditions, the building wall heating and the street ground heating are shown to evidently influence the strength and location of the predominant recirculation inside the street canyon, especially for the windward wall heating. Moreover, analysis results have shown that the arcade width has larger influences on wind and thermal environment than that of the arcade height. Particularly, the wind velocity increases 2.1 to 3.8 times, the air temperature decreases 0.42–2.24 °C and the Physiological Equivalent Temperature value drops between 2.2 °C and 7.0 °C when the arcade width increases to half of building width.

Numerical Heat Transfer Simulation of Oil Shale Large-Size Downhole Heater

Downhole heaters are critical for effectively achieving in situ oil shale cracking. In this study, we simulate the heat transfer performance of a large-scale helical baffle downhole heater under various operational conditions. The findings indicate that at 160 m3/h and 6 kW the outlet temperature can reach 280 °C. Controlling heating power or increasing the injected gas flow effectively mitigates heat accumulation on the heating rod’s surface. The outlet temperature curve exhibits two phases. Simultaneously, a balance in energy exchange between the injected gas and heating power occurs, mitigating high-temperature hotspots. Consequently, the outlet temperature cannot attain the theoretical maximum temperature, referred to as the actual maximum temperature. Employing h/∆p13 as the indicator to evaluate heat transfer performance, optimal performance occurs at 100 m3/h. Heat transfer performance at 200 m3/h is significantly impacted by heating power, with the former being approximately 6% superior to the latter. Additionally, heat transfer performance is most stable below 160 m3/h. The gas heating process is categorized into three stages based on temperature distribution characteristics within the heater: rapid warming, stable warming, and excessive heating. The simulation findings suggest that the large-size heater can inject a higher flow rate of heat-carrying gas into the subsurface, enabling efficient oil shale in situ cracking.

Numerical simulation of heat transfer during spark plasma sintering of porous SiC

The strong covalent bonding of SiC makes it difficult to be fully sintered using conventional sintering techniques. Therefore, electrical current-assisted sintering techniques such as spark plasma sintering are needed to sinter SiC. The finite numerical simulations have been generally used to predict and guide sintering processes of ceramics (i.e.TiC,AlO3). In this work, the numerical simulations of the temperature distribution, electrical field and thermal field on the porous SiC cylindrical samples by Comsol were investigated in details. The influence of porosity on the temperature distribution and electrical field distribution of the SiC samples were analyzed. Both the maximum current density (1.19 × 107A/m2) and the maximum temperature (2030 °C) of the SPS device occurred at the punch/spacer interface. The temperature values at the centre of the samples were proportional to the porosity, and the temperatures at the center of the samples increased from ∼1960 °C to ∼2010 °C at t = 2000s, with the porosity increasing from 1% to 60%. The uniformity of the temperature distribution deteriorated as the porosity increases, and the temperature gradients between the center and edge of the samples with the porosity of 1% and 60%, were ∼80 °C and ∼130 °C at t = 2000s, respectively.

Flow Reversal of Unsteady Oscillatory Flow in Combined Convection Fully Developed Vertical Channel

Numerical simulations of flow and heat transfer in an unsteady oscillatory fully developed combined convection are conducted. The temperature on the right wall varies with sinusoidal rhythm over time within a locked mean temperature, whereas the temperature on the left wall remains constant. Both walls are also considered to be stationary. A numerical solution for the momentum and energy equations obtained using the Runga- Kutta method is used to analyze the effect of periodic oscillation on the velocity and temperature configurations, as well as the progression of heat generations. The combined convection process of heat transfer and its flow pattern in a vertical channel is important in many applications, including environmental, industrial, and engineering. This study is prompted by a problem with the heat transfer process, which is difficult, expensive, and time-consuming. The purpose of this study is to develop a mathematical model, find a numerical solution, and conduct parametric studies on the double-diffusive fully developed vertical channel with unsteady oscillatory flow, as well as to investigate the effect of various parameters on the velocity and temperature profiles. The effect of different dimensional parameters is tested and the flow reversal phenomenon is discussed. The dimensional equations are transformed into non-dimensional equations using the similarity technique. Then, the built-in program dsolve in Maple is used to numerically solve the boundary value problem (BVP). A validation study is conducted on a previously reported problem to confirm the validity of the present calculation. The flow and temperature numerical findings are represented visually. It is found that a large development of flow is caused by the plate temperature oscillation with high amplitude, therefore, the nature of the flow differs from the non-oscillating case.

Numerical simulation of heat transfer and entropy generation due to the nanofluid natural convection with viscous dissipation in an inclined square cavity

In this work, the effects of viscous dissipation on heat transfer and entropy generation of the nanofluid natural convection in an inclined square cavity are numerically investigated. A fractional-step semiimplicit finite difference algorithm based on the projection method on a staggered grid is proposed to solve the laminar natural convection problem, which has the advantage of larger time steps. The square cavity is filled with a nanofluid composed of (Cu) copper nanoparticles and water. and the viscous dissipative behavior of the mixture flow is not negligible. The studied parameters are Rayleigh number Ra ( 10 4 and 10 5 ), Eckert number Ec ( 0 − 2 ), the volume fraction of solid particles ϕ ( 0 − 0.06 ), and inclination angle of square cavity α ( 0 − π / 2 ). The Prandtl number is fixed to Pr = 6.2 . The results show that at any inclination angle, the increase in viscous dissipation leads to weakened heat transfer on the hot wall, enhanced heat transfer on the cold wall, and weakened flow in the square cavity. For the base solution without nanoparticles, the Eckert number has the greatest impact, with its effects on the average Nusselt number on the hot wall, maximum streamfunction, and average Nusselt number on the cold wall being ∼13, 0.15, 18.69%, and at Ra = 10 4 , and 12, 0.5, 28.87% at Ra = 10 5 , respectively. Research on entropy generation shows that as Eckert number increases, the entropy generation due to heat transfer increases, while due to fluid friction decreases. As the Rayleigh number increases, the effect of viscous dissipation increases. As the inclination angle increases, the effects of volume fraction and Eckert number weaken. Adding solid particles can effectively weaken the effect of viscous dissipation.

Influence of the Shrinkage of the Inner Layer of Steel Tubes on Permissible Thermal Load

This paper presents the results of numerical simulations of heat transfer in a tube with a protective chromium layer on the inner surface made of steel, with different shrinkage temperatures. Shrinkage in the steel is an unfavorable phenomenon because it causes cracks in the chrome coating. The cracks contribute to the peeling of the protective material on the inner surface of the tube. Wear and damage to the chrome layer significantly shortens the tube’s service life. The influence of the type of steel with medium carbon content on heat transfer for a sequence of ten, thirty, and sixty heat impulses was examined. Simulations were carried out for two selected steels with clearly different shrinkage temperatures, i.e., 30HN2MFA and X37CrMoV5-1 (1.2343). In 30HN2MFA steel, the shrinkage effect occurred at a temperature of approx. 750 °C, while in X37CrMoV5-1 steel it occurred at a temperature of approx. 870 °C. In the computational model, 30 cross-sections of the three-meter-long tube were analyzed. A time-dependent heat transfer coefficient was calculated in each zone. Heat transfer simulations were carried out using COMSOL version 6.1 software. This paper shows that for X37CrMoV5-1 steel, the shrinkage temperature on the inner surface of the tube was reached after approx. 60 heat impulses, while for 30HN2MFA steel, it was reached after approx. 30. The greatest differences in the number of impulses occurred for a pipe with a 200 µm thick chrome layer.

A comprehensive fluid–solid coupling dynamic simulation for spatiotemporal distribution of regression rate in hybrid rocket motors

The spatiotemporal distribution characteristics of the regression rate are crucial aspects of the research on Hybrid Rocket Motor (HRM). This study presents a pioneering effort in achieving a comprehensive numerical simulation of fluid dynamics and heat transfer in both the fluid and solid regions throughout the entire operation of an HRM. To accomplish this, a dynamic grid technique that incorporates fluid–solid coupling is utilized. To validate the precision of the numerical simulations, a firing test is conducted, with embedded thermocouple probes being used to measure the inner temperature of the fuel grain. The temperature variations in the solid fuel obtained from both experiment and simulations show good agreement. The maximum combustion temperature and average thrust obtained from the simulations are found to deviate from the experimental results by only 3.3% and 2.4%, respectively. Thus, it can be demonstrated that transient numerical simulations accurately capture the fluid–solid coupling characteristics and transient regression rate. The dynamic simulation results of inner flow field and solid region throughout the entire working stage reveal that the presence of vortices enhances the blending of combustion gases and improves the regression rate at both the front and rear ends of the fuel grain. In addition, oscillations of the regression rate obtained in the simulation can also be well corresponded with the corrugated surface observed in the experiment. Furthermore, the zero-dimension regression rate formula and the formula describing the axial location dependence of the regression rate are fitted from the simulation results, with the corresponding coefficients of determination (R2) of 0.9765 and 0.9298, respectively. This research serves as a reference for predicting the performance of HRM with gas oxygen and polyethylene, and presents a credible way for investigating the spatiotemporal distribution of the regression rate.

Numerical simulation of hybrid nanofluid flow and heat transfer across parallel surfaces with suction/injection and magnetic effect

The heat and mass transfer through the hybrid nanofluid (Hnf) flow with the significance of magnetic field across two spinning parallel plates has been described. The Hnf is prepared by the dispersion of SiO2 (Silicon dioxide) and MoS2 (Molybdenum Disulfide) nanoparticles (NPs) in the ethylene glycol (EG). The MoS2 is anti-friction compound used in automobiles and other types of heavy machineries in industry. It reduces the engine noise, reduces fuel consumption and enhances engine life. Similarly, SiO2 is used in structural materials, food processing, pharmaceutical industries and as an electrical insulator in microelectronic apparatus. Based on the remarkable applications of the hybrid nanofluid (MoS2-SiO2/EG), the flow has been modeled in form of PDEs, which are numerically handled through the parametric method (PCM). The results are compared for velocity, energy and concentration profiles with another numerical technique bvp4c (Matlab code). It has been detected that the results derived from the PCM are reliable and accurate. Furthermore, the velocity field declines with the upshot of Reynold number and suction/injection factor. The heat dissemination rate enhances from 2.85% to 9.89%, whereas the mass diffusion rate enriches form 2.32% to 9.56% as the values of nanoparticles varies from 0.01 to 0.03.