Heat Transfer Evaluation on Curved Boundaries in Thermal Lattice Boltzmann Equation Method

2013 ◽  
Vol 136 (1) ◽  
Author(s):  
Like Li ◽  
Renwei Mei ◽  
James F. Klausner

An efficient and accurate approach for heat transfer evaluation on curved boundaries is proposed in the thermal lattice Boltzmann equation (TLBE) method. The boundary heat fluxes in the discrete velocity directions of the TLBE model are obtained using the given thermal boundary condition and the temperature distribution functions at the lattice nodes close to the boundary. Integration of the discrete boundary heat fluxes with effective surface areas gives the heat flow rate across the boundary. For lattice models with square or cubic structures and uniform lattice spacing the effective surface area is constant for each discrete heat flux, thus the heat flux integration becomes a summation of all the discrete heat fluxes with constant effective surface area. The proposed heat transfer evaluation scheme does not require a determination of the normal heat flux component or a surface area approximation on the boundary; thus, it is very efficient in curved-boundary simulations. Several numerical tests are conducted to validate the applicability and accuracy of the proposed heat transfer evaluation scheme, including: (i) two-dimensional (2D) steady-state thermal flow in a channel, (ii) one-dimensional (1D) transient heat conduction in an inclined semi-infinite solid, (iii) 2D transient heat conduction inside a circle, (iv) three-dimensional (3D) steady-state thermal flow in a circular pipe, and (v) 2D steady-state natural convection in a square enclosure with a circular cylinder at the center. Comparison between numerical results and analytical solutions in tests (i)–(iv) shows that the heat transfer is second-order accurate for straight boundaries perpendicular to one of the discrete lattice velocity vectors, and first-order accurate for curved boundaries due to the irregularly distributed lattice fractions intersected by the curved boundary. For test (v), the computed surface-averaged Nusselt numbers agree well with published results.

Author(s):  
Like Li ◽  
Renwei Mei ◽  
James F. Klausner

The evaluation of the boundary heat flux and total heat transfer in the lattice Boltzmann equation (LBE) simulations is investigated. The boundary heat fluxes in the discrete velocity directions of the thermal LBE (TLBE) model are obtained directly from the temperature distribution functions at the lattice nodes. With the rectangular lattice uniformly spaced the effective surface area for the discrete heat flux is the unit spacing distance, thus the heat flux integration becomes simply a summation of all the discrete heat fluxes with constant surface areas. The present method for the evaluation of total heat transfer is very efficient and robust for curved boundaries because it does not require the determination of the normal heat flux on the boundary and the surface area. To validate its applicability and accuracy, several numerical tests with analytical solutions are conducted, including 2-dimensional (2D) steady thermal flow in a channel, 1-D transient heat conduction in an inclined semi-infinite solid, 2-D transient conduction inside a circle, and 3-D steady thermal flow in a circular pipe. For straight boundaries perpendicular to one of the discrete velocity vectors, the total heat transfer is second-order accurate. For curved boundaries only first-order accuracy is obtained for the total heat transfer due to the irregularly distributed lattice fractions cut by the curved boundary.


Fluids ◽  
2021 ◽  
Vol 6 (7) ◽  
pp. 246
Author(s):  
Rozie Zangeneh

The Wall-modeled Large-eddy Simulation (WMLES) methods are commonly accompanied with an underprediction of the skin friction and a deviation of the velocity profile. The widely-used Improved Delayed Detached Eddy Simulation (IDDES) method is suggested to improve the prediction of the mean skin friction when it acts as WMLES, as claimed by the original authors. However, the model tested only on flow configurations with no heat transfer. This study takes a systematic approach to assess the performance of the IDDES model for separated flows with heat transfer. Separated flows on an isothermal wall and walls with mild and intense heat fluxes are considered. For the case of the wall with heat flux, the skin friction and Stanton number are underpredicted by the IDDES model however, the underprediction is less significant for the isothermal wall case. The simulations of the cases with intense wall heat transfer reveal an interesting dependence on the heat flux level supplied; as the heat flux increases, the IDDES model declines to predict the accurate skin friction.


2013 ◽  
Vol 135 (6) ◽  
Author(s):  
Paul M. Kodzwa ◽  
John K. Eaton

This paper presents isoenergetic temperature and steady-state film-cooled heat transfer coefficient measurements on the pressure surface of a modern, highly cambered transonic airfoil. A single passage model simulated the idealized two-dimensional flow path between blades in a modern transonic turbine. This set up offered a simpler construction than a linear cascade but produced an equivalent flow condition. Furthermore, this model allowed the use of steady-state, constant surface heat fluxes. We used wide-band thermochromic liquid crystals (TLCs) viewed through a novel miniature periscope system to perform high-accuracy (±0.2 °C) thermography. The peak Mach number along the pressure surface was 1.5, and maximum turbulence intensity was 30%. We used air and carbon dioxide as injectant to simulate the density ratios characteristic of the film cooling problem. We found significant differences between isoenergetic and recovery temperature distributions with a strongly accelerated mainstream and detached coolant jets. Our heat transfer data showed some general similarities with lower-speed data immediately downstream of injection; however, we also observed significant heat transfer attenuation far downstream at high blowing conditions. Our measurements suggested that the momentum ratio was the most appropriate variable to parameterize the effect of injectant density once jet lift-off occurred. We noted several nonintuitive results in our turbulence effect studies. First, we found that increased mainstream turbulence can be overwhelmed by the local augmentation of coolant injection. Second, we observed complex interactions between turbulence level, coolant density, and blowing rate with an accelerating mainstream.


Author(s):  
V. G. Razumovskiy ◽  
Eu. N. Pis’mennyy ◽  
A. Eu. Koloskov ◽  
I. L. Pioro

The results of heat transfer to supercritical water flowing upward in a vertical annular channel (1-rod channel) and tight 3-rod bundle consisting of the tubes of 5.2-mm outside diameter and 485-mm heated length are presented. The heat-transfer data were obtained at pressures of 22.5, 24.5, and 27.5 MPa, mass flux within the range from 800 to 3000 kg/m2·s, inlet temperature from 125 to 352°C, outlet temperature up to 372°C and heat flux up to 4.6 MW/m2 (heat flux rate up to 2.5 kJ/kg). Temperature regimes of the annular channel and 3-rod bundle were stable and easily reproducible within the whole range of the mass and heat fluxes, even when a deteriorated heat transfer took place. The data resulted from the study could be applicable for a reference estimation of heat transfer in future designs of fuel bundles.


Author(s):  
Jensen Hoke ◽  
Todd Bandhauer ◽  
Jack Kotovsky ◽  
Julie Hamilton ◽  
Paul Fontejon

Liquid-vapor phase change heat transfer in microchannels offers a number of significant advantages for thermal management of high heat flux laser diodes, including reduced flow rates and near constant temperature heat rejection. Modern laser diode bars can produce waste heat loads >1 kW cm−2, and prior studies show that microchannel flow boiling heat transfer at these heat fluxes is possible in very compact heat exchanger geometries. This paper describes further performance improvements through area enhancement of microchannels using a pyramid etching scheme that increases heat transfer area by ∼40% over straight walled channels, which works to promote heat spreading and suppress dry-out phenomenon when exposed to high heat fluxes. The device is constructed from a reactive ion etched silicon wafer bonded to borosilicate to allow flow visualization. The silicon layer is etched to contain an inlet and outlet manifold and a plurality of 40μm wide, 200μm deep, 2mm long channels separated by 40μm wide fins. 15μm wide 150μm long restrictions are placed at the inlet of each channel to promote uniform flow rate in each channel as well as flow stability in each channel. In the area enhanced parts either a 3μm or 6μm sawtooth pattern was etched vertically into the walls, which were also scalloped along the flow path with the a 3μm periodicity. The experimental results showed that the 6μm area-enhanced device increased the average maximum heat flux at the heater to 1.26 kW cm2 using R134a, which compares favorably to a maximum of 0.95 kw cm2 dissipated by the plain walled test section. The 3μm area enhanced test sections, which dissipated a maximum of 1.02 kW cm2 showed only a modest increase in performance over the plain walled test sections. Both area enhancement schemes delayed the onset of critical heat flux to higher heat inputs.


1968 ◽  
Vol 90 (1) ◽  
pp. 51-54 ◽  
Author(s):  
W. A. Beckman

The one-dimensional steady-state temperature distribution within an isotropic porous bed subjected to a collimated and/or diffuse radiation heat flux and a transparent flowing fluid has been determined by numerical methods. The porous bed was assumed to be nonscattering and to have a constant absorption coefficient. Part of the radiation absorbed by the porous bed is reradiated and the remainder is transferred to the fluid by convection. Due to the assumed finite volumetric heat transfer coefficient, the bed and fluid have different temperatures. A bed with an optical depth of six and with a normal incident collimated radiation heat flux was investigated in detail. The radiation incident on the bed at the fluid exit was assumed to originate from a black surface at the fluid exit temperature. The investigation covered the range of incident diffuse and collimated radiation heat fluxes expected in a nonconcentrating solar energy collector. The results are presented in terms of a bed collection efficiency from which the fluid temperature rise can be calculated.


Author(s):  
Han Wang ◽  
Qincheng Bi ◽  
Linchuan Wang ◽  
Haicai Lv ◽  
Laurence K. H. Leung

An experiment has recently been performed at Xi’an Jiaotong University to study the wall temperature and pressure drop at supercritical pressures with upward flow of water inside a 2×2 rod bundle. A fuel-assembly simulator with four heated rods was installed inside a square channel with rounded corner. The outer diameter of each heated rod is 8 mm with an effective heated length of 600 mm. Experimental parameters covered the pressure of 23–28 MPa, mass flux of 350–1000 kg/m2s and heat flux on the rod surface of 200–1000 kW/m2. According to the experimental data, it was found that the circumferential wall temperature distribution of a heated rod is not uniform. The temperature difference between the maximum and the minimum varies with heat flux and/or mass flux. Heat transfer characteristics of supercritical water in bundle were discussed with respect to various heat fluxes. The effect of heat flux on heat transfer in rod bundles is similar with that in tubes or annuli. In addition, flow resistance reflected in the form of pressure loss has also been studied. Experimental results showed that the total pressure drop increases with bulk enthalpy and mass flux. Four heat transfer correlations developed for supercritical pressures water were compared with the present test data. Predictions of Jackson correlation agrees closely with the experimental data.


2004 ◽  
Vol 27 (4) ◽  
pp. 319-339 ◽  
Author(s):  
Sutthisak Phongthanapanich ◽  
Pramote Dechaumphai

A finite element method is combined with the Delaunay triangulation and an adaptive remeshing technique to solve for solutions of both steady-state and transient heat conduction problems. The Delaunay triangulation and the adaptive remeshing technique are explained in detail. The solution accuracy and the effectiveness of the combined procedure are evaluated by heat transfer problems that have exact solutions. These problems include steady-state heat conduction in a square plate subjected to a highly localized surface heating, and a transient heat conduction in a long plate subjected to a moving heat source. The examples demonstrate that the adaptive remeshing technique with the Delaunay triangulation significantly reduce the number of the finite elements required for the problems and, at the same time, increase the analysis solution accuracy as compared to the results produced using uniform finite element meshes.


1994 ◽  
Vol 76 (5) ◽  
pp. 2084-2094 ◽  
Author(s):  
M. B. Ducharme ◽  
P. Tikuisis

The objective of the present study was to investigate the relative contribution of the convective heat transfer in the forearm and hand to 1) the total heat loss during partial immersion in cold water [water temperature (Tw) = 20 degrees C] and 2) the heat gained during partial immersion in warm water (Tw = 38 degrees C). The heat fluxes from the skin of the forearm and finger were continuously monitored during the 3.5-h immersion of the upper limb (forearm and hand) with 23 recalibrated heat flux transducers. The last 30 min of the partial immersion were conducted with an arterial occlusion of the forearm. The heat flux values decreased during the occlusion period at Tw = 20 degrees C and increased at Tw = 38 degrees C for all sites, plateauing only for the finger to the value of the tissue metabolic rate (124.8 +/- 29.0 W/m3 at Tw = 20 degrees C and 287.7 +/- 41.8 W/m3 at Tw = 38 degrees C). The present study shows that, at thermal steady state during partial immersion in water at 20 degrees C, the convective heat transfer between the blood and the forearm tissue is the major heat source of the tissue and accounts for 85% of the total heat loss to the environment. For the finger, however, the heat produced by the tissue metabolism and that liberated by the convective heat transfer are equivalent. At thermal steady state during partial immersion in water at 38 degrees C, the blood has the role of a heat sink, carrying away from the limb the heat gained from the environment and, to a lesser extent (25%), the metabolic and conductive heats. These results suggest that during local cold stress the convective heat transfer by the blood has a greater role than that suggested by previous studies for the forearm but a lesser role for the hand.


Sign in / Sign up

Export Citation Format

Share Document