Fully-developed heat transfer in annuli for viscoelastic fluids with viscous dissipation

Analytical solutions are obtained for heat transfer in concentric annular flows of viscoelastic fluids modeled by the simplified Phan-Thien–Tanner constitutive equation. Solutions for thermal and dynamic fully developed flow are presented for both imposed constant wall heat fluxes and imposed constant wall temperatures, always taking into account viscous dissipation.Equations are presented for the normalized temperature profile, the bulk temperature, the inner and outer wall temperatures and, through their definitions for the inner and outer Nusselt numbers as a function of all relevant non-dimensional parameters. Some special results are discussed in detail. Given the complexity of the derived equations, for ease of use compact exact expressions are presented for the Nusselt numbers and programmes to calculate all quantities are made accessible on the internet. Generally speaking, fluid elasticity is found to increase the heat transfer for imposed heating at the wall, especially in combination with internal heat generation by viscous dissipation, whereas for imposed wall temperatures it reduces heat transfer when viscous dissipation is weak.     

DNS of turbulent channel flow with conjugate heat transfer: Effect of thermal boundary conditions on the second moments and budgets

Budgets of turbulent heat fluxes and temperature variance obtained from the Direct Numerical Simulation of an incompressible periodic channel flow with a Reynolds number of 150 (based on friction velocity) and a Prandtl number of 0.71 are presented and analysed for four cases: locally imposed temperature at the wall (constant Dirichlet), locally imposed heat flux (constant Neumann), heat exchange coefficient (Robin) and 3D conjugate heat transfer. The dissipation rate associated with the temperature variance is strongly impacted by the thermal boundary condition. For non-conjugate cases, a straightforward analytical analysis establishes the connection between the boundary condition, the temperature variance and the wall-normal part of the thermal dissipation rate at the wall. For the conjugate case, the two-point correlations of the thermal field in the solid domain confirms the existence of very large scale thermal structures.     

Flow and heat transfer of an electrically conducting fluid of second grade in a porous medium over a stretching sheet subject to a transverse magnetic field

An analysis is performed for flow and heat transfer of a steady laminar boundary layer flow of an electrically conducting fluid of second grade in a porous medium subject to a transverse uniform magnetic field past a semi-infinite stretching sheet with power-law surface temperature or power-law surface heat flux. The effects of viscous dissipation, internal heat generation of absorption and work done due to deformation are considered in the energy equation. The variations of surface temperature gradient for the prescribed surface temperature case (PST) and surface temperature for the prescribed heat flux case (PHF) with various parameters are tabulated. The asymptotic expansions of the solutions for large Prandtl number are also given for the two heating conditions. It is shown that, when the Eckert number is large enough, the heat flow may transfer from the fluid to the wall rather than from the wall to the fluid when Eckert number is small. A physical explanation is given for this phenomenon.     

A Local Thermal Non-Equilibrium Analysis of Fully Developed Forced Convective Flow in a Tube Filled with a Porous Medium

A local thermal non-equilibrium model has been considered for the case of thermally fully developed flow within a constant heat flux tube filled with a porous medium. Exact temperature profiles for the fluid and solid phases are found after combining the two individual energy equations and then transforming them into a single ordinary differential equation with respect to the temperature difference between the solid phase and the wall subject to constant heat flux. The exact solutions for the case of metal-foam and air combination reveal that the local thermal equilibrium assumption may fail for the case of constant heat flux wall. The Nusselt number is presented as a function of the Peclet number, which shows a significant increase due to both high stagnant thermal conductivity and thermal dispersion resulting from the presence of the metal-foam.     

Conjugate forced convection and heat conduction in circular microchannels

The effects of axial heat conduction in the solid walls of microchannels of circular cross-sections are analyzed here. A systematic approach is adopted, with the aim of pointing out the influence of geometrical parameters and of solid wall thermal conductivity on microchannel heat transfer. The reliability of a commonly adopted criterium, based on the so-called axial conduction number, to assess the relevance of axial heat conduction is also discussed. Numerical simulations concern the simultaneously developing laminar flow of a constant property fluid in microchannels of different length, wall thickness and wall material, heated with a uniform heat flux at the outer surface, for different values of the Reynolds number. Moreover, since often in experimental tests the two end sections of the microchannel wall are not perfectly insulated, the effects of heat losses through these sections are also considered. A hybrid finite element procedure, which implies the step-by-step solution of the parabolized momentum equations in the fluid domain, followed by the solution of the energy equation in the entire domain, corresponding to both the solid and the fluid parts, is used for the numerical simulations.     

A high-order discontinuous Galerkin method for simulating incompressible fluid-thermal-structural interaction problems

This paper presents a framework for the application of the discontinuous Galerkin(DG) finite element method to the multi-physics simulation of the solid thermal deformation interacting with incompressible flow problems in two-dimensions. Recent applications of the DG method are primarily for thermoelastic problems in a solid domain or fluid-structure interaction problems without heat transfer. Based on a recently published conjugate heat transfer solver, the incompressible Navier-Stokes equation, the fluid advection-diffusion equation, the Boussinesq term, the solid heat equation and the solid linear elastic equation are solved using an explicit DG formulation. A Dirichlet-Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux computed at interface quadrature points in the fluid-solid interface. Formal hp convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved for each solver. The algorithm is then further validated against several existing benchmark cases including the in-plane loaded square, the Timoshenko Beam, the laminated beam subject to thermal-loads and the lid-driven cavity with a flexible bottom wall. The computational effort demonstrates that for all cases the highest order accurate algorithm has several magnitudes lower error than the second-order schemes for a given computational effort. It is a strong justification for the development of such high order discretisations. The solver can be employed to predict thermal deformation of a structure due to convective and conductive heat transfer at low Mach, such as chip deformation on a printed circuit board, wave-guide structure optimization, thermoelectric cooler simulation, and optics mounting method verification.     

DNS,LES and RANS of turbulent heat transfer in boundary layer with suddenly changing wall thermal conditions

The objectives of this study are to investigate a thermal field in a turbulent boundary layer with suddenly changing wall thermal conditions by means of direct numerical simulation (DNS), and to evaluate predictions of a turbulence model in such a thermal field, in which DNS of spatially developing boundary layers with heat transfer can be conducted using the generation of turbulent inflow data as a method. In this study, two types of wall thermal condition are investigated using DNS and predicted by large eddy simulation (LES) and Reynolds-averaged Navier–Stokes equation simulation (RANS). In the first case, the velocity boundary layer only develops in the entrance of simulation, and the flat plate is heated from the halfway point, i.e., the adiabatic wall condition is adopted in the entrance, and the entrance region of thermal field in turbulence is simulated. Then, the thermal boundary layer develops along a constant temperature wall followed by adiabatic wall. In the second case, velocity and thermal boundary layers simultaneously develop, and the wall thermal condition is changed from a constant temperature to an adiabatic wall in the downstream region. DNS results clearly show the statistics and structure of turbulent heat transfer in a constant temperature wall followed by an adiabatic wall. In the first case, the entrance region of thermal field in turbulence can be also observed. Thus, both the development and the entrance regions in thermal fields can be explored, and the effects upstream of the thermal field on the adiabatic region are investigated. On the other hand, evaluations of predictions by LES and RANS are conducted using DNS results. The predictions of both LES and RANS almost agree with the DNS results in both cases, but the predicted temperature variances near the wall by RANS give different results as compared with DNS. This is because the dissipation rate of temperature variance is difficult to predict by the present RANS, which is found by the evaluation using DNS results.     

Numerical analysis of transient laminar forced convection of nanofluids in circular ducts

In this study, forced convection heat transfer characteristics of nanofluids are investigated by numerical analysis of incompressible transient laminar flow in a circular duct under step change in wall temperature and wall heat flux. The thermal responses of the system are obtained by solving energy equation under both transient and steady-state conditions for hydro-dynamically fully-developed flow. In the analyses, temperature dependent thermo-physical properties are also considered. In the numerical analysis, Al₂O₃/water nanofluid is assumed as a homogenous single-phase fluid. For the effective thermal conductivity of nanofluids, Hamilton–Crosser model is used together with a model for Brownian motion in the analysis which takes the effects of temperature and the particle diameter into account. Temperature distributions across the tube for a step jump of wall temperature and also wall heat flux are obtained for various times during the transient calculations at a given location for a constant value of Peclet number and a particle diameter. Variations of thermal conductivity in turn, heat transfer enhancement is obtained at various times as a function of nanoparticle volume fractions, at a given nanoparticle diameter and Peclet number. The results are given under transient and steady-state conditions; steady-state conditions are obtained at larger times and enhancements are found by comparison to the base fluid heat transfer coefficient under the same conditions.     

TEMPERATURE PROFILES OF LOCAL THERMAL NONEQUILIBRIUM FOR THERMAL DEVELOPING FORCED CONVECTION IN POROUS MEDIUM PARALLEL PLATE CHANNEL

Based on the two-energy equation model, taking into account viscous dissipation due to the interaction between solid skeleton and pore fluid flow, temperature expressions of the solid skeleton and pore fluid flow are obtained analytically for the thermally developing forced convection in a saturated porous medium parallel plate channel, with walls being at constant temperature. It is proved that the temperatures of the two phases for the local thermal nonequilibrium approach to the temperature derived from the one-energy equation model for the local thermal equilibrium when the heat exchange coefficient goes to infinite. The temperature profiles are shown in figures for different dimensionless parameters and the effects of the parameters on the local thermal nonequilibrium are revealed by parameter study.     

Effects of pressure work and viscous dissipation on Graetz problem for gas flows in parallel-plate channels

The analytical solutions are obtained for the Graetz problem with pressure work and viscous dissipation in the thermal entrance region of the parallel-plate channels for two basic boundary conditions of uniform wall temperature and uniform wall heat flux involving fully developed laminar gas flows. The asymptotic Nusselt number is found to be zero instead of the conventionally accepted value of 7.54 for the uniform wall temperature case and (140/17)/ [1+(27/17) PrEc] for uniform wall heat flux case. The effects of pressure work and viscous dissipation contribute significantly to the asymptotic results for heat transfer and cannot be neglected under any circumstances in the case of uniform wall temperature. Sample results are presented to illustrate the effects of pressure work and viscous dissipation on heat transfer characteristics in the thermal entrance region.     

Modeling the turbulent heat and momentum transfer in flows under different thermal conditions

Two-equation turbulence models for velocity and temperature (scalar) fields are developed to calculate wall shear flows under various flow conditions and related turbulent heat transfer under various wall thermal conditions. In the present models, we make the modified dissipation rates of both turbulent energy and temperature variance zero at a wall, though the wall limiting behavior of velocity and temperature fluctuations is reproduced exactly. Thus, the models assure computational expediency and convergence. Also, the present k- model is construted using a new type of expression for the Reynolds stress proposed by Abe et al. [Trans. JSME B 61 (1995) 1714–1721], whose essential feature lies in introducing the explicit algebraic stress model concept into the nonlinear k- formulation, and the present two-equation heat transfer model is constructed to properly take into account the effects of wall thermal conditions on the eddy diffusivity for heat. The models are tested with five typical velocity fields and four typical thermal fields. Agreement with experiment and direct simulation data is quite satisfactory.     

Heat Transfer Effects on a Fully Premixed Methane Impinging Flame

A numerical assessment of different thermal conditions for an impinging flame configuration is investigated using large-eddy simulation. The cases of study correspond to a turbulent methane flame at equivalence ratio ER = 0.8 and temperature T = 298 K exiting at 30 m/s with a nozzle-to-plate distance over diameter of H/D = 2. Computational cases based on different thermal conditions are compared to a conjugate case, in which fluid and solid domains are solved simultaneously. A solid material defined with enhanced conductivity properties is used to represent the wall in the conjugate case, so that the characteristic time scales of the solid are reduced. The results indicate that the heat transfer condition influences not only the mean temperature and gradients, but also the temperature fluctuations in the near-wall region. It is found that adiabatic, isothermal and more sophisticated Robin-type boundary conditions contribute to underpredict/overpredict the temperature field near the wall. As the time scales of fluid and solid are of the same order, the use of conjugate approaches is required to predict the correct flow fields near the wall and the skin temperature.     

Laminar forced convection with viscous dissipation in a concentric annular ductConvection laminaire forcée avec dissipation visqueuse dans un conduit annulaire concentrique

Forced convection heat transfer in fully developed flows of viscous dissipating fluids in concentric annular ducts is analyzed analytically. Special attention has been paid to the effect of the viscous dissipation. Two different cases of the thermal boundary conditions are considered: uniform heat flux at the outer wall and adiabatic inner wall (Case A) and uniform heat flux at the inner wall and adiabatic outer wall (Case B). Solutions for the velocity and temperature distributions and the Nusselt number are obtained for different values of the aspect ratio and the Brinkman number. The present analytical results for the case without the viscous dissipation effect are compared with those available in the literature and an excellent agreement is observed. To cite this article: M. Avc?, O. Ayd?n, C. R. Mecanique 334 (2006).     

On selection of reference temperature of heat transfer coefficient for complicated flows

Convective heat transfer coefficient is closely related with flow and thermal conditions. To define heat transfer coefficient, a reference temperature needs to be properly selected, which can be the fluid bulk mean temperature for internal flows or the temperature at the far field for external flows. For complicated flows, the adiabatic wall temperature is commonly adopted as the reference temperature, while other options can also be applied. This paper analyzed some of the potential selections of the reference temperature for different flow settings, including film cooling, jet impingement with cross flows, and a mixing flow in a straight duct with or without internal heat source. It is observed that heat transfer coefficient changes dramatically with selection of reference temperatures. In case of constant wall temperature, using adiabatic wall temperature as reference temperature can result in negative heat transfer coefficient, which means the heat flux has a different direction with the defined driving temperature difference. To avoid the inconsistency due to the reference temperature, an innovative method is proposed to calculate the heat transfer coefficient of complicated flows.     

Two-phase annular flow and evaporative heat transfer in a microchannel

A physical and mathematical model has been developed to predict the two-phase flow and heat transfer in a microchannel with evaporative heat transfer. Sample solutions to the model were obtained for both constant wall temperature and constant wall heat flux conditions. Results are provided for evaporation rate, liquid film thickness, liquid and vapor phase pressure and temperature distributions. In addition to the sample calculations that were used to illustrate the transport characteristics, computations based on the current model were performed to generate results for comparisons with the experimental results of Qu and Mudawar (2004) where two different mass flow rates of the working fluid were used in the experiment. The comparisons of total pressure drops with the experimental data of Qu and Mudawar (2004) cover the wall heat flux range of 142.71-240 W/cm² with a total channel mass flux of 400.1 kg/m² s and also the wall heat flu range of 99.54-204.39 W/cm² with total channel mass flux of 401.9 kg/m² s. The calculated results from the current model match closely with those of Qu and Mudawar (2004).     

Multiphysics Simulation Combining Large-Eddy Simulation,Wall Heat Conduction and Radiative Energy Transfer to Predict Wall Temperature Induced by a Confined Premixed Swirling Flame

A multi-physics simulation combining large-eddy simulation, conjugate heat transfer and radiative heat transfer is used to predict the wall temperature field of a confined premixed swirling flame operating under atmospheric pressure. The combustion model accounts for the effect of enthalpy defect on the flame structure whose stabilization is here sensitive to the wall heat losses. The conjugate heat transfer is accounted for by solving the heat conduction within the combustor walls and with the Hybrid-Cell Neumann-Dirichlet coupling method, enabling to dynamically adapt the coupling period. The latter coupling procedure is enhanced to determine statistics (mean, RMS, \(\ldots \)) in a permanent regime accurately and efficiently thanks to an acceleration technique which is derived and validated. The exact radiative heat transfer equation is solved with an advanced Monte Carlo method with a local control of the statistical error. The coupled simulation is carried out with or without accounting for radiation. Excellent results for the wall temperature are achieved by the fully coupled simulation which are then further analyzed in terms of radiative effects, global energy budget and fluctuations of wall heat flux and temperature.     

单层微通道热沉拓扑结构对芯片中心热点的降温效能研究

通过热流固耦合模拟分析得到了不同微通道结构热沉基底的温度场及微通道内速度场,研究了相同入流功率下不同单层微通道拓扑结构对中心有高热流密度热点芯片的散热能力。结果表明:相同入流功率(0.05W)下,不同结构的散热能力排序由高到低为Y分形、弯曲散射、直散射(双侧出流)、直螺旋、直散射(单侧出流)、圆螺旋、树状分形、直槽结构;采用中心入流可有效降低芯片中心热点附近的温度,对于中心入流的散射结构,采用对称出流结构可提升其流动传热性能;Y分形结构具有良好的流动传热特性,对于热源面和中心热点均具有良好的散热效果。     

Effects of viscous dissipation on laminar forced convection with axially periodic wall heat flux

The laminar forced convection in a circular duct is investigated in the case of a sinusoidal axial variation of the wall heat flux. The axial heat conduction in the fluid is neglected, while the effect of viscous dissipation is taken into account. The heat transfer in the thermally developed region, where the temperature is the sum of a linear function and a periodic function of the axial coordinate, is analysed. Both the temperature field and the local Nusselt number are evaluated analytically. Comparisons with the solution in the absence of viscous heating are performed. It is shown that the effect of viscous dissipation on the temperature field may be relevant especially in the case of a sinusoidal wall heat flux distribution with a vanishing mean value. Received on 24 July 1998     

Forced convection heat transfer of Giesekus viscoelastic fluid in pipes and channels

A theoretical solution is presented for the convective heat transfer of Giesekus viscoelastic fluid in pipes and channels, under fully developed thermal and hydrodynamic flow conditions, for an imposed constant heat flux at the wall. The fluid properties are taken as constant and axial conduction is negligible. The effect of Weissenberg number (We), mobility parameter (α) and Brinkman number (Br) on the temperature profile and Nusselt number are investigated. The results emphasize the significant effect of viscous dissipation and fluid elasticity on the Nusselt number in all circumstances. For wall cooling and the Brinkman number exceeds a critical value (Br ₁), the heat generated by viscous dissipation overcomes the heat removed at the wall and fluid heats up longitudinally. Fluid elasticity shifts this critical Brinkman number to higher values.     

Thermo-mechanical modeling of turbulent heat transfer in gas–solid flows including particle collisions

A thermo-mechanical turbulence model is developed and used for predicting heat transfer in a gas–solid flow through a vertical pipe with constant wall heat flux. The new four-way interaction model makes use of the thermal k_θ–τ_θ equations, in addition to the hydrodynamic k–τ transport, and accounts for the particle–particle and particle–wall collisions through a Eulerian/Lagrangian formulation. The simulation results indicate that the level of thermal turbulence intensity and the heat transfer are strongly affected by the particle collisions. Inter-particle collisions attenuate the thermal turbulence intensity near the wall but somewhat amplify the temperature fluctuations in the pipe core region. The hydrodynamic-to-thermal times-scale ratio and the turbulent Prandtl number in the region near the wall increase due to the inter-particle collisions. The results also show that the use of a constant or the single-phase gas turbulent Prandtl number produces error in the thermal eddy diffusivity and thermal turbulent intensity fields. Simulation results also indicate that the inter-particle contact heat conduction during collision has no significant effect in the range of Reynolds number and particle diameter studied.