On an indirect integral equation approach for stationary heat transfer in semi-infinite layered domains in R3 with cavities

Regions containing internal boundaries such as composite materials arise in many applications.We consider a situation of a layered domain in IR3 containing a nite number of bounded cavities. The model is stationary heat transfer given by the Laplace equation with piecewise constant conductivity. The heat ux (a Neumann condition) is imposed on the bottom of the layered region and various boundary conditions are imposed on the cavities. The usual transmission (interface) conditions are satised at the interface layer, that is continuity of the solution and its normal derivative. To eciently calculate the stationary temperature eld in the semi-innite region, we employ a Green's matrix technique and reduce the problem to boundary integral equations (weakly singular) over the bounded surfaces of the cavities. For the numerical solution of these integral equations, we use Wienert's approach [20]. Assuming that each cavity is homeomorphic with the unit sphere, a fully discrete projection method with super-algebraic convergence order is proposed. A proof of an error estimate for the approximation is given as well. Numerical examples are presented that further highlights the eciency and accuracy of the proposed method.

Analogue simulation of heat transfer in the thermal processing of food

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

Polynomial Expansions for Solutions of Higher-Order Bessel Heat Equation in Quantum Calculus

Mathematics Subject Class.: 33C10,33D60,26D15,33D05,33D15,33D90

Numerical solution of the Dirichlet initial boundary value problem for the heat equation in exterior 3-dimensional domains using integral equations

A numerical method for the Dirichlet initial boundary value problem for the heat equation in the exterior and unbounded region of a smooth closed simply connected 3-dimensional domain is proposed and investigated. This method is based on a combination of a Laguerre transformation with respect to the time variable and an integral equation approach in the spatial variables. Using the Laguerre transformation in time reduces the parabolic problem to a sequence of stationary elliptic problems which are solved by a boundary layer approach giving a sequence of boundary integral equations of the first kind to solve. Under the assumption that the boundary surface of the solution domain has a one-to-one mapping onto the unit sphere, these integral equations are transformed and rewritten over this sphere. The numerical discretisation and solution are obtained by a discrete projection method involving spherical harmonic functions. Numerical results are included.

Embedded heat pipes in cofired ceramic substrates for enhanced thermal management of electronics

A novel and new thermal management technology for advanced ceramic microelectronic packages has been developed incorporating miniature heat pipes embedded in the ceramic substrate. The heat pipes use an axially grooved wick structure and water as the working fluid. Prototype substrate/heat pipe systems were fabricated using high temperature co-fired ceramic (alumina). The heat pipes were nominally 81 mm in length, 10 mm in width, and 4 mm in height, and were charged with approximately 50–80 μL of water. Platinum thick film heaters were fabricated on the surface of the substrate to simulate heat dissipating electronic components. Several thermocouples were affixed to the substrate to monitor temperature. One end of the substrate was affixed to a heat sink maintained at constant temperature. The prototypes were tested and shown to successful and reliably operate with thermal loads over 20 Watts, with thermal input from single and multiple sources along the surface of the substrate. Temperature distributions are discussed for the various configurations and the effective thermal resistance of the substrate/heat pipe system is calculated. Finite element analysis was used to support the experimental findings and better understand the sources of the system's thermal resistance. ^

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Cancer therapy combining the modalities of hyperthermia and chemotherapy: In vitro cellular response after rapid heat accumulation in the cancer cell

Hyperthermia is usually used at a sub-lethal level in cancer treatment to potentiate the effects of chemotherapy. The purpose of this study is to investigate the role of heating rate in achieving synergistic cell killing by chemotherapy and hyperthermia. For this purpose, in vitro cell culture experiments with a uterine cancer cell line (MES-SA) and its multidrug resistant (MDR) variant MES-SA/Dx5 were conducted. The cytotoxicity, mode of cell death, induction of thermal tolerance and P-gp mediated MDR following the two different modes of heating were studied. Doxorubicin (DOX) was used as the chemotherapy drug. Indocyanine green (ICG), which absorbs near infrared light at 808nm (ideal for tissue penetration), was chosen for achieving rapid rate hyperthermia. A slow rate hyperthermia was provided by a cell culture incubator. The results show that the potentiating effect of hyperthermia to chemotherapy can be maximized by increasing the rate of heating as evident by the results from the cytotoxicity assay. When delivered at the same thermal dose, a rapid increase in temperature from 37°C to 43°C caused more cell membrane damage than gradually heating the cells from 37°C to 43°C and thus allowed for more intracellular accumulation of the chemotherapeutic agents. Different modes of cell death are observed by the two hyperthermia delivery methods. The rapid rate laser-ICG hyperthermia @ 43°C caused cell necrosis whereas the slow rate incubator hyperthermia @ 43°C induced very mild apoptosis. At 43°C a positive correlation between thermal tolerance and the length of hyperthermia exposure is identified. This study shows that by increasing the rate of heating, less thermal dose is needed in order to overcome P-gp mediated MDR.

Severe Heat Cramps in a High School Football Player: A Case Report

We present a case study of an adolescent football player who suffered from severe full body muscle cramping after supplementing with creatine for two months. A paucity of data exists regarding the safety of creatine supplementation and its side effects on dehydration, body fluid/electrolyte balance, and other heat illnesses.

Embedded Heat Pipes in Cofired Ceramic Substrates for Enhanced Thermal Management of Electronics

A novel and new thermal management technology for advanced ceramic microelectronic packages has been developed incorporating miniature heat pipes embedded in the ceramic substrate. The heat pipes use an axially grooved wick structure and water as the working fluid. Prototype substrate/heat pipe systems were fabricated using high temperature co-fired ceramic (alumina). The heat pipes were nominally 81 mm in length, 10 mm in width, and 4 mm in height, and were charged with approximately 50-80 mL of water. Platinum thick film heaters were fabricated on the surface of the substrate to simulate heat dissipating electronic components. Several thermocouples were affixed to the substrate to monitor temperature. One end of the substrate was affixed to a heat sink maintained at constant temperature. The prototypes were tested and shown to successful and reliably operate with thermal loads over 20 Watts, with thermal input from single and multiple sources along the surface of the substrate. Temperature distributions are discussed for the various configurations and the effective thermal resistance of the substrate/heat pipe system is calculated. Finite element analysis was used to support the experimental findings and better understand the sources of the system's thermal resistance.

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Heat transport in low dimensions : introduction and phenomenology

Review paper, to appear in the Springer Lecture Notes in Physics volume "Thermal transport in low dimensions: from statistical physics to nanoscale heat transfer" (S. Lepri ed.)

Heat Transfer in an Airfoil Shaped Strut

The heat transfer from a hot primary flow stream passing over the outside of an airfoil shaped strut to a cool secondary flow stream passing through the inside of that strut was studied experimentally and numerically. The results showed that the heat transfer on the inside of the strut could be reliably modeled as a developing flow and described using a power law model. The heat transfer on the outside of the strut was complicated by flow separation and stall on the suction side of the strut at high angles of attack. This separation was quite sensitive to the condition of the turbulence in the flow passing over the strut, with the size of the separated wake changing significantly as the mean magnitude and levels of anisotropy were varied. The point of first stall moved by as much as 15% of the chord, while average heat transfer levels changed by 2-5% as the inlet condition was varied. This dependence on inlet conditions meant that comparisons between experiment and steady RANS based CFD were quite poor. Differences between the CFD and experiment were attributed to anisotropic and unsteady effects. The coupling between the two flows was shown to be quite low - that is to say, heat transfer coefficients on both the inner and outer surfaces of the strut were relatively unaffected by the temperature of the strut, and it was possible to predict the temperature on the strut surface quite reliably using heat transfer data from decoupled tests, especially for CFD simulations.

SIGNIFICANCE OF HEAT DISSIPATION IN THE CATHODE CATALYST LAYER ON THE LIFETIME AND RELIABILITY OF POLYMER ELECTROLYTE FUEL CELLS

To study the dissipation of heat generated due to the formation of pinholes that cause local hotspots in the catalyst layer of the Polymer Electrolyte Fuel Cell, a two-phase non-isothermal model has been developed by coupling Darcy’s law with heat transport. The domain under consideration is a section of the membrane electrode assembly with a half-channel and a half-rib. Five potential locations where a pinhole might form were analyzed: at the midplane of the channel, midway between the channel midplane and the channel wall, at the channel or rib wall, midway between the rib midplane and the channel wall, at the midplane of the rib. In the first part of this work, a preliminary thermal model was developed. The model was then refined to account for the two-phase effects. A sensitivity study was done to evaluate the effect of the following properties on the maximum temperature in the domain: Catalyst layer thermal conductivity, the Microporous layer thermal conductivity, the anisotropy factor of the Catalyst layer thermal conductivity, the Porous transport layer porosity, the liquid water distribution and the thickness of the membrane and porous layers. Accounting for the two-phase effects, a slight cooling effect was observed across all hotspot locations. The thermal properties of the catalyst layer were shown to have a limited impact on the maximum temperature in the catalyst layer of new fuel cells without pinhole. However, as hotspots start to appear, thermal properties play a more significant role in mitigating the thermal runaway.

Comparison of DES predictions of turbulent flow and heat transfer in a 45° ribbed duct

Numerical predictions of the turbulent flow and heat transfer of a stationary duct with square ribs 45° angled to the main flow direction are presented. The rib height to channel hydraulic diameter is 0.1, the rib pitch to rib height is 10. The calculations have been carried out for a bulk Reynolds number of 50,000. The flows generated by ribs are dominated by separating and reattaching shear layers with vortex shedding and secondary flows in the cross-section. The hybrid RANS-LES approach is adopted to simulate such flows at a reasonable computation cost. The capability of the various versions of DES method, depending the RANS model, such as DES-SA, DES-RKE, DES-SST, have been compared and validated against the experiment. The significant effect of RANS model on the accuracy of the DES prediction has been shown. The DES-SST method, which was able to reproduce the correct physics of flow and heat transfer in a ribbed duct showed better performance than others.