Computations were run to study heat transfer coefficient augmentation with film cooling for a simulated gas turbine blade leading edge. The realizable k-[epsilon] turbulence model (RKE) and Shear Stress Transport k-[omega] turbulence model (SST) were used for the computational simulations. RKE computations completed at a unity density ratio were confirmed to be consistent with experimental measurements conducted by Yuki et al.(1998) and Johnston et al. (1999) whereas SST computations exhibited significant discrepancies. Moreover the effect of the density ratio on heat transfer coefficient augmentation was studied because experimental measurements of heat transfer coefficient augmentation with film cooling are generally constrained to unity density ratio tests. It was shown that heat transfer coefficient augmentation can be simulated using unity density ratio jets, but only when scaled with the momentum flux ratio of the coolant jets. / text

The one dimensional transient heat conduction equation was numerically modeled through matrix diagonalization and three time-discretization schemes. The discrete methods were first-order backward, second-order backward, and implicit finite difference schemes. All simulations used the central difference formula in the space dimension. Two relevant physical systems were considered: a uniformly conducting slab and a melting block of ice. The latter lead to a moving boundary system, or Stefan problem. The multiphysics of melting was numerically modeled through alternating updates of temperature and melt front profiles. Iterative simulations were run with regularly refined discretization meshes in both systems. In the case of the conducting slab, temperature at a fixed point in space and time was considered. For the Stefan problem, the melt front movement after a set time was the physical solution of interest. The accuracy of the convergent results was increased using Richardson acceleration and the Wynn's epsilon algorithm. Accuracy was improved for the moving boundary problem as well, but to a significantly lesser degree. The relative errors improved by five and two orders of magnitude for the conducting block and melting ice simulations, respectively. These relative errors were used to determine that matrix diagonalization is the most accurate numerical solution among the four considered. In both simulation convergence and acceleration potential, matrix diagonalization was superior to the implicit and explicit discretization solutions. However, matrix diagonalization required significantly more computational time. With the enhancement of convergence acceleration, the finite difference schemes obtained similar relative errors to the diagonalization model. This demonstrated the value of convergence acceleration in the classic dilemma for every programmer. There is always a balance struck between model sophistication, accuracy, and computational time. Convergence acceleration allows for a simpler numerical model to achieve comparable accuracy, and in less time than that required for sophisticated numerical models. The numerical models were also compared for stability through parameters that arose in each simulation. These parameters were the Courant-Friederichs-Lewy (CFL) condition and diagonalized eigenvalues. Though diagonalization was found to be the most accurate, it was determined that the backwards finite difference solutions are the easiest to evaluate for stability. In these solution methods, the CFL value allows the stability to be determined prior to running the simulation.

This thesis is a study of the heat transfer correlations between a wire and liquid helium cooled to either 1.9 or 4.3 K. The wire resembles a part of a superconducting magnet used in the Large Hadron Collider (LHC) particle accelerator currently being built at CERN. The magnets are cooled to 1.9 K and using helium as a coolant is very efficient, especially at extremely low temperatures since it then becomes a superfluid with an apparent infinite thermal conductivity. The cooling of the magnet is very important, since the superconducting wires need to be thermally stable. Thermal stability means that a superconductive magnet can remain superconducting, even if a part of the magnet becomes normal conductive due to a temperature increase. This means that if heat is generated in a wire, it must be transferred to the helium by some sort of heat transfer mechanism, or along the wire or to the neighbouring wires by conduction. Since the magnets need to be superconductive for the operation of the particle accelerator, it is crucial to keep the wires cold. Therefore, it is necessary to understand the heat transfer mechanisms from the wires to the liquid helium. The scope of this thesis was to describe the heat transfer mechanisms from a heater immersed in liquid and superfluid helium. By performing both experiments and simulations, it was possible to determine properties like heat transfer correlations, critical heat flux limits, and the differences between transient and steady-state heat flow. The measured values were in good agreement with values found in literature with a few exceptions. These differences could be due to measurement errors. A numerical program was written in Matlab and it was able to simulate the experimental temperature and heat flux response with good accuracy for a given heat generation.

In the future, new challenges will occur during the product development in the vehicular industry when emission legislations getting tighter. This will also affect the truck cooling system and therefore increase needs for analysing the system at different levels of the product development. Volvo 3P wishes for these reasons to examine the possibility to use AMESim as a future 1D analysis tool. This tool can be used as a complement to existing analysis methods at Volvo 3P. It should be possible to simulate pressure, flow and heat transfer both steady state and transient. In this thesis work a cooling system of a FH31 MD13 520hp truck with an engine driven coolant pump is studied. Further a model of the cooling system is built in AMESim together with necessary auxiliary system such as oil circuits. The model is validated using experimental data that have been produced by Volvo 3P at the Gothenburg facility. The results from validation and other simulations show that the model gives a good picture of the cooling system. It also gives information about pressure, flow and heat transfer in steady state conditions. Further a design modification is done, showing how a change affects the flow in the cooling system. The conclusion is that a truck cooling system can be built and simulated in AMESim. Further, it shows that AMESim meets the requirements Volvo 3P in Gothenburg has set up for the future 1D analysis tool and thereby AMESim is a good complement to the already existing analysis method.

Mist or aerosol explosions present a serious hazard to process industries. Heat transfer fluids are widely used in the chemical process industry, are flammable above their flash points, and can cause aerosol explosions. Though the possibility of aerosol explosions has been widely documented, knowledge about their explosive potential is limited. Studying the formation of such aerosols by emulating leaks in process equipment will help define a source term for aerosol dispersions and aid in characterizing their explosion hazards. Analysis of the problem of aerosol explosions reveals three major steps: source term calculations, dispersion modeling, and explosion analysis. The explosion analysis, consisting of ignition and combustion, is largely affected by the droplet size distribution of the dispersed aerosol. The droplet size distribution of the dispersed aerosol is a function of the droplet size distribution of the aerosol formed from the leak. Existing methods of dealing with the problem of aerosol explosions are limited to enhancing the dispersion to prevent flammable concentrations and use of explosion suppression mechanisms. Insufficient data and theory on the flammability limits of aerosols renders such method speculative at best. Preventing the formation of aerosol upon leaking will provide an inherently safer solution to the problem. The research involves the non-intrusive measurement of heat transfer fluid aerosol sprays using a Malvern Diffraction Particle Analyzer. The aerosol is generated by plain orifice atomization to simulate the formation and dispersion of heat transfer fluid aerosols through leaks in process equipment. Predictive correlations relating aerosol droplet sizes to bulk liquid pressures, temperatures, thermal and fluid properties, leak sizes, and ambient conditions are presented. These correlations will be used to predict the conditions under which leaks will result in the formation of aerosols and will ultimately help in estimating the explosion hazards of heat transfer fluid aerosols. Heat transfer fluid selection can be based on liquids that are less likely to form aerosols. Design criteria also can incorporate the data to arrive at operating conditions that are less likely to produce aerosols. The goal is to provide information that will reduce the hazards of aerosol explosions thereby improving safety in process industries.

The aim of the work is to investigate heat transfer of the staggered tube bank to downward statically stable foam flow.

Sensible energy storage for solar domestic hot water (SDHW) systems is space consuming and heavy. Latent heat energy storage systems (LHESSs) offer a solution to this problem. However, the functionality of a LHESS during simultaneous charging/discharging, an operating mode encountered when used with a SDHW, had not been studied experimentally. A small scale vertical cylindrical LHESS, with dodecanoic acid as the phase change material (PCM), was studied during separate and simultaneous charging/discharging. Natural convection was found to have a strong influence during melting, but not during solidification. During simultaneous operation heat transfer was limited by the high thermal resistance of the solid PCM. However, when the PCM was melted, direct heat transfer occurred between the hot and cold heat transfer fluids, indicating the significance of the PCM phase on heat transfer in the system. The results of this research will lead to more optimally designed LHESS for use with SDHW. ?

This thesis work presents a numerical study of the effects of fins and thermal fluid velocities on the storage characteristics of a cylindrical latent heat energy storage system (LHESS). The work consists of two main components: 1. The development of a numerical method to study and solve the phase change heat transfer problems encountered in a LHESS during charging of the system, which results in melting of the phase change material (PCM). The numerical model is based on the finite element method. The commercial software COMSOL Multiphysics was used to implement it. The effective heat capacity method was applied in order to account for the large amount of latent energy stored during melting of a PCM, and the moving interface between the solid and liquid phases. The fluid flow, heat transfer and phase change processes were all validated using known analytical solutions or correlations. 2. Due to the low thermal conductivity of PCMs, the heat transfer characteristics of an enhanced LHESS was studied numerically. The effects of fins and the thermal fluid velocity on the melting rate of the PCM in the LHESS were analyzed. Results obtained for configurations having between 0 and 27 fins show that the heat transfer rate increases with addition of fins and thermal fluid velocity. The effect of the HTF velocity was observed to be small with few fin configurations since the thermal resistance offered by the LHESS system, mostly PCM, is vastly more important under these conditions; while its effect becomes more pronounced with addition of fins, since the overall thermal resistance decreases greatly with the addition of fins. The total energy stored after 12 hours for 0 and 27 fins configurations range between 3.6 MJ and 39.7 MJ for a thermal fluid velocity of 0.05 m/s and between 3.7 MJ and 57 MJ for a thermal fluid velocity of 0.5 m/s. The highest system efficiencies for the 0.05 m/s and 0.5 m/s, obtained with 27 fins configuration are 68.9% and 97.9% respectively.

This work presents an experimental and numerical study of phase change behaviour in a horizontal cylindrical latent heat energy storage system (LHESS). Fins with two orientations, straight fins and angled fins, are added into the PCM to enhance heat transfer. The PCM used in this study is lauric acid which has desirable thermal properties for LHESS. The experimental work concentrates on studying the heat transfer mechanism during phase change, impacts of the HTF inlet temperature and HTF flow rates. Moreover, heat transfer enhancement effectiveness of straight fins and angles fins is compared. Numerical model is simulated using COMSOL Multiphysics software package. It is observed that conduction is the dominant heat transfer mechanism during the initial stage of charging, and natural convection plays a more important role afterwards. Conduction plays a major role during solidification. Complete melting time is affected by the HTF inlet temperature and HTF flow rates.

Bubble nucleation at subcooled flow boiling condition in a horizontal annular channel with a square cross section by the use of high-speed camera is investigated. The channel represents a scaled-down version of a single rod of CANDU reactor core. The experiments were performed by the use of water at pressures between 1-3 atm, constant heat flux of 0.124 MW/m2, liquid bulk subcooling of 32-1oC and mean flow velocities of 0.3-0.4 m/s. Bubble lift-off diameters were obtained from direct high speed videography. The developed model for the bubble lift-off diameter was obtained by analyzing the forces acting on a bubble. Furthermore, a model for the bubble growth rate constant was suggested. The proposed model was then compared to experimental data and it has shown a good agreement with the experimental data. Additionally, the effects of liquid bulk subcooling, liquid pressure and mean flow velocity on bubble lift-off diameter were investigated.