Abstract: This study presents the numerical simulation of three-dimensional incompressible steady and laminar fluid flow and conjugate heat transfer of a trapezoidal microchannel heat sink using water as a cooling fluid in a silicon substrate. Navier-Stokes equations with conjugate energy equation are discretized by finite-volume method. We perform numerical computations for a range of 50 ≦ Re ≦ 600, 0.05W ≦ P ≦ 0.8W, 20W/cm2 ≦q"≦ 40W/cm2. The present study demonstrates the numerical optimization of a trapezoidal microchannel heat sink design using the response surface methodology (RSM) and the genetic algorithm method (GA). The results show that the average Nusselt number increases with an increase in the Reynolds number or pumping power, and the thermal resistance decreases as the pumping power increases. The thermal resistance of a trapezoidal microchannel is minimized for a constant heat flux and constant pumping power.
Abstract: A 3D-conjugate numerical investigation was conducted to predict heat transfer characteristics in a rectangular cross-sectional micro-channel employing simultaneously developing Tow-phase flows. The sole purpose for analyzing two phase flow heat transfer in rectangular micro channel is to pin point what are the different factors affecting this phenomenon. Different methods and techniques have been undertaken to analyze the equations arising constituting the flow of heat from gas phase to liquid phase and vice versa.Different models of micro channels have been identified and analyzed. How the geometry of micro channels affects their activity i.e. of circular and non-circular geometry has also been reviewed. To the study the results average Nusselt no plotted against the Reynolds no has been taken into consideration to study average heat exchange in micro channels against applied heat flux. High heat fluxes up to 140 W/cm2 were applied to investigate micro-channel thermal characteristics.
Abstract: The ice rink floor is the largest heat exchanger in an ice rink. The important part of the floor consists of concrete, and the thermophysical properties of this concrete have strong influence on the energy usage of the ice rink. The thermal conductivity of concrete can be increased by using iron ore as ballast. In this study, the Transient Plane Source (TPS) method showed an increase up to 58.2% of thermal conductivity comparing the improved concrete to standard concrete. Moreover, two alternative ice rink floor designs are suggested to incorporate the improved concrete. A 2D simulation was developed to investigate the temperature distribution in the conventional and the suggested designs. The results show that the suggested designs reduce the temperature difference between the ice surface and the brine by 1-4˚C, when comparing with convectional designs at equal heat flux. This primarily leads to an increased coefficient of performance (COP) in the primary refrigeration cycle and secondly to a decrease in the secondary refrigerant pumping power. The suggested designs have great potential to reduce the energy usage of ice rinks. Depending on the load scenario in the ice rink, the saving potential lies in the range of 3-10% of the refrigeration system energy usage. This calculation is based on steady state conditions and the potential with improved dynamic behavior is expected to increase the potential saving.
Abstract: This paper presents a study on the effect of
second-order slip on forced convection through a long isoflux heated
or cooled planar microchannel. The fully developed solutions of flow
and thermal fields are analytically obtained on the basis of the
second-order Maxwell-Burnett slip and local heat flux boundary
conditions. Results reveal that when the average flow velocity
increases or the wall heat flux amount decreases, the role of thermal
creep becomes more insignificant, while the effect of second-order slip
becomes larger. The second-order term in the Deissler slip boundary
condition is found to contribute a positive velocity slip and then to lead
to a lower pressure drop as well as a lower temperature rise for the
heated-wall case or to a higher temperature rise for the cooled-wall
case. These findings are contrary to predictions made by the
Karniadakis slip model.
Abstract: A mathematical heat transfer model for the prediction of transient heating of the slab in a direct-fired walking beam type reheating furnace has been developed by considering the nongray thermal radiation with given furnace environments. The furnace is modeled as radiating nongray medium with carbon dioxide and water with five-zoned gas temperature and the furnace wall is considered as a constant temperature lower than furnace gas one. The slabs are moving with constant velocity depending on the residence time through the non-firing, charging, preheating, heating, and final soaking zones. Radiative heat flux obtained by considering the radiative heat exchange inside the furnace as well as convective one from the surrounding hot gases are introduced as boundary condition of the transient heat conduction within the slab. After validating thermal radiation model adopted in this work, thermal fields in both model and real reheating furnace are investigated in terms of radiative heat flux in the furnace and temperature inside the slab. The results show that the slab in the furnace can be more heated with higher slab emissivity and residence time.
Abstract: Implicit in most large-scale numerical analyses of the crystal growth from the melt is the assumption that the shape and position of the phase boundary are determined by the transport phenomena coupled strongly to the melt hydrodynamics. In the present numerical study, the interface shape-effect on the convective interactions in a Czochralski oxide melt is described. It was demonstrated that thermocapillary flow affects inversely the phase boundaries of distinct shapes. The inhomogenity of heat flux and the location of the stagnation point at the crystallization front were investigated. The forced convection effect on the point displacement at the boundary found to be much stronger for the flat plate interface compared to the cone-shaped one with and without the Marangoni flow.
Abstract: The radiative heat transfer problem is investigated numerically for 2D complex geometry biomass pyrolysis reactor composed of two pyrolysis chambers and a heat recuperator. The fumes are a mixture of carbon dioxide and water vapor charged with absorbing and scattering particles and soot. In order to increase gases residence time and heat transfer, the heat recuperator is provided with many inclined, vertical, horizontal, diffuse and grey baffles of finite thickness and has a complex geometry. The Finite Volume Method (FVM) is applied to study radiative heat transfer. The blocked-off region procedure is used to treat the geometrical irregularities. Eight cases are considered in order to demonstrate the effect of adding baffles on the walls of the heat recuperator and on the walls of the pyrolysis rooms then choose the best case giving the maximum heat flux transferred to the biomass in the pyrolysis chambers. Ray effect due to the presence of baffles is studied and demonstrated to have a crucial effect on radiative heat flux on the walls of the pyrolysis rooms. Shadow effect caused by the presence of the baffles is also studied. The non grey radiative heat transfer is studied for the real existent configuration. The Weighted Sum of The Grey Gases (WSGG) Model of Kim and Song is used as non grey model. The effect of soot volumetric fraction on the non grey radiative heat flux is investigated and discussed.
Abstract: This work deals with heat and mass transfer by steady laminar boundary layer flow of a Newtonian, viscous fluid over a vertical flat plate with variable surface heat flux embedded in a fluid saturated porous medium in the presence of thermophoresis particle deposition effect. The governing partial differential equations are transformed into no-similar form by using special transformation and solved numerically by using an implicit finite difference method. Many results are obtained and a representative set is displaced graphically to illustrate the influence of the various physical parameters on the wall thermophoresis deposition velocity and concentration profiles. It is found that the increasing of thermophoresis constant or temperature differences enhances heat transfer rates from vertical surfaces and increase wall thermophoresis velocities; this is due to favorable temperature gradients or buoyancy forces. It is also found that the effect of thermophoresis phenomena is more pronounced near pure natural convection heat transfer limit; because this phenomenon is directly a temperature gradient or buoyancy forces dependent. Comparisons with previously published work in the limits are performed and the results are found to be in excellent agreement.
Abstract: Laminar mixed Convection heat transfer of a nanofluid with prescribed constant heat flux on the inner wall of horizontal annular tube has been studied numerically based on two-phase mixture model in different Rayleigh Numbers and Azimuth angles. Effects of applying of different volume fractions of Al2O3 nanoparticles in water as a base fluid on hydrodynamic and thermal behaviors of the fluid flow such as axial velocity, secondary flow, temperature, heat transfer coefficient and friction coefficient at the inner and outer wall region, has been investigated. Conservation equations in elliptical form has been utilized and solved in three dimensions for a steady flow. It is observed that, there is a good agreement between results in this work and previously published experimental and numerical works on mixed convection in horizontal annulus. These particles cause to increase convection heat transfer coefficient of the fluid, meanwhile there is no considerable effect on friction coefficient.
Abstract: The unsteady transient free convection flow of an incompressible dissipative viscous fluid between parallel plates at different distances have been investigated under porous medium. Due to presence of heat flux under the influence of uniform transverse magnetic field the velocity distribution and the temperature distribution, is shown graphically. Since exact solution is not possible so we find parametrical solution by perturbation technique. The result is shown in graph for different parameters. We notice that heat generation effects fluid velocity keeping in which of free convection which cools.
Abstract: Heat transfer and laminar fluid flow over backward facing step with and without obstacle numerically studied in this paper. The finite volume method adopted to solve continuity, momentum and energy equations in two dimensions. Backward facing step without obstacle and with different dimension of obstacle were presented. The step height and expansion ratio of channel were 4.8mm and 2 respectively, the range of Reynolds number varied from 75 to 225, constant heat flux subjected on downstream of wall was 2000W/m2, and length of obstacle was 1.5, 3, and 4.5mm with width 1.5mm. The separation length noticed increase with increase Reynolds number and height of obstacle. The result shows increase of heat transfer coefficient for backward facing step with obstacle in compared to those without obstacle. The maximum enhancement of heat transfer observed at 4.5mm of height obstacle due to increase recirculation flow after the obstacle in addition that at backward. Streamline of velocity showing the increase of recirculation region with used obstacle in compared without obstacle and highest recirculation region observed at obstacle height 4.5mm. The amount of enhancement heat transfer was varied between 3-5% compared to backward without obstacle.
Abstract: In this study the augmentation of heat transfer in a rectangular channel with triangular vortex generators is evaluated. The span wise averaged Nusselt number, mean temperature and total heat flux are compared with and without vortex generators in the channel at a blade angle of 30° for Reynolds numbers 800, 1200, 1600, and 2000. The use of vortex generators increases the span wise averaged Nusselt number compared to the case without vortex generators considerably. At a particular blade angle, increasing the Reynolds number results in an enhancement in the overall performance and span wise averaged Nusselt number was found to be greater at particular location for larger Reynolds number. The total heat flux from the bottom wall with vortex generators was found to be greater than that without vortex generators and the difference increases with increase in Reynolds number.
Abstract: The study focuses to investigate the thermal response of delaminations and develop mathematical models using numerical results to obtain the optimum heat requirement and time to identify delaminations in GLARE type of Fibre Metal Laminates (FML) in both reflection mode and through-transmission (TT) mode of step pulsed active thermography (SPAT) method in the type of nondestructive testing and evaluation (NDTE) technique. The influence of applied heat flux and time on various sizes and depth of delaminations in FML is analyzed to investigate the thermal response through numerical simulations. A finite element method (FEM) is applied to simulate SPAT through ANSYS software based on 3D transient heat transfer principle with the assumption of reflection mode and TT mode of observation individually.
The results conclude that the numerical approach based on SPAT in reflection mode is more suitable for analysing smaller size of near-surface delaminations located at the thermal stimulator side and TT mode is more suitable for analysing smaller size of deeper delaminations located far from thermal stimulator side or near thermal detector/Infrared camera side. The mathematical models provide the optimum q and T at the required MRTD to identify unidentified delamination 7 with 25015.0022W/m2 at 2.531sec and delamination 8 with 16663.3356 W/m2 at 1.37857sec in reflection mode. In TT mode, the delamination 1 with 34954W/m2 at 13.0399sec, delamination 2 with 20002.67W/m2 at 1.998sec and delamination 7 with 20010.87 W/m2 at 0.6171sec could be identified.
Abstract: In this research, the capability of neural networks in
modeling and learning complicated and nonlinear relations has been
used to develop a model for the prediction of changes in the diameter
of bubbles in pool boiling distilled water. The input parameters used
in the development of this network include element temperature, heat
flux, and retention time of bubbles. The test data obtained from the
experiment of the pool boiling of distilled water, and the
measurement of the bubbles form on the cylindrical element. The
model was developed based on training algorithm, which is
typologically of back-propagation type. Considering the correlation
coefficient obtained from this model is 0.9633. This shows that this
model can be trusted for the simulation and modeling of the size of
bubble and thermal transfer of boiling.
Abstract: In this research, the changes in bubbles diameter and
number that may occur due to the change in heat flux of pure water
during pool boiling process. For this purpose, test equipment was
designed and developed to collect test data. The bubbles were graded
using Caliper Screen software. To calculate the growth and
nucleation rates of bubbles under different fluxes, population balance
model was employed. The results show that the increase in heat flux
from q=20 kw/m2 to q= 102 kw/m2 raised the growth and nucleation
rates of bubbles.
Abstract: In this paper the influence of a vertical plate’s thermal capacity is numerically investigated in order to evaluate the evolution of the thermal boundary layer structure, as well as the convective heat transfer coefficient and the velocity and temperature profiles. Whereas the heat flux of the heated vertical plate is evaluated under time depending boundary conditions. The main important feature of this problem is the unsteadiness of the physical phenomena. A 2D CFD model is developed with the Ansys Fluent 14.0 environment and is validated using unsteady data obtained for plasterboard studied under a dynamic temperature evolution. All the phenomena produced in the vicinity of the thermal conductive vertical plate (plasterboard) are analyzed and discussed. This work is the first stage of a holistic research on transient free convection that aims, in the future, to study the natural convection in the vicinity of a vertical plate containing Phase Change Materials (PCM).
Abstract: Heat transfer and laminar air flow over a double backward-facing step numerically studied in this paper. The simulations was performed by using ANSYS ICEM for meshing process and using ANSYS fluent 14 (CFD) for solving. The k-ɛ standard model adopted with Reynolds number varied between 98.5 to 512 and three step height at constant heat flux (q=2000 W/m2). The top of wall and bottom of upstream are insulated with bottom of downstream is heated. The results show increase in Nusselt number with increases of Reynolds number for all cases and the maximum of Nusselt number happens at the first step in compared to the second step. Due to increase of cross section area of downstream to generate sudden expansion then Nusselt number decrease but the profile of Nusselt number keep same trend for all cases where increase after the first and second steps. Recirculation region after the first and second steps are denoted by contour of streamline velocity. The higher augmentation of heat transfer rate observed for case 1 at Reynolds number of 512 and heat flux q=2000 W/m2.
Abstract: An experimental investigation was conducted to study the effect of surface roughness on friction factor and heat transfer characteristics in single-phase fluid flow in a stainless steel micro-tube having diameter of 0.85 mm and average internal surface roughness of 1.7 μm with relative surface roughness of 0.002. Distilled water and R134a liquids were used as the working fluids and testing was conducted with Reynolds numbers ranging from 100 to 10,000 covering laminar, transition and turbulent flow conditions. The experiments were conducted with the micro-tube oriented horizontally with uniform heat fluxes applied at the test section. The results indicated that the friction factor of both water and R134a can be predicted by the Hagen-Poiseuille equation for laminar flow and the modified Miller correlation for turbulent flow and early transition from laminar to turbulent flows. The heat transfer results of water and R134a were in good agreement with the conventional theory in the laminar flow region and lower than the Adam’s correlation for turbulent flow region which deviates from conventional theory.
Abstract: With the long-term objective of Critical Heat Flux (CHF) prediction, a Direct Numerical Simulation (DNS) framework for simulation of subcooled and saturated nucleate pool boiling is developed. One case of saturated, and three cases of subcooled boiling at different subcooling levels are simulated. Grid refinement study is also reported. Both boiling and condensation phenomena can be computed simultaneously in the proposed numerical framework. Computed bubble detachment diameters of the saturated nucleate pool boiling cases agree well with the experiment. The flow structures around the growing bubble are presented and the accompanying physics is described. The relation between heat flux evolution from the heated wall and the bubble growth is studied, along with investigations of temperature distribution and flow field evolutions.
Abstract: Combined conduction-free convection heat transfer in
vertical eccentric annuli is numerically investigated using a finitedifference
technique. Numerical results, representing the heat transfer
parameters such as annulus walls temperature, heat flux, and heat
absorbed in the developing region of the annulus, are presented for a
Newtonian fluid of Prandtl number 0.7, fluid-annulus radius ratio 0.5,
solid-fluid thermal conductivity ratio 10, inner and outer wall
dimensionless thicknesses 0.1 and 0.2, respectively, and
dimensionless eccentricities 0.1, 0.3, 0.5, and 0.7. The annulus walls
are subjected to thermal boundary conditions, which are obtained by
heating one wall isothermally whereas keeping the other wall at inlet
fluid temperature. In the present paper, the annulus heights required
to achieve thermal full development for prescribed eccentricities are
obtained. Furthermore, the variation in the height of thermal full
development as function of the geometrical parameter, i.e.,
eccentricity is also investigated.