Proper design of high performance industrial heat transfer equipment relies on accurate knowledge and prediction of the thermal boundary conditions. In order to enhance the overall gas turbine efficiency, advancements in cooling technology for gas turbines and related applications are continuously investigated to increase the turbine inlet temperature without compromising the durability of the materials used. For detailed design, local distributions are needed in addition to bulk quantities. Detailed local distributions require advanced experimental techniques whereas they are readily available using numerical tools. Numerical predictions using a computational fluid dynamics approach with popular turbulence models are benchmarked against a semi-empirical correlation for the friction in a circular channel with repeated-rib roughness to demonstrate some shortcomings of the models used. Numerical predictions varied widely depending on the turbulence modelling approach used. The need for a compatible experimental dataset to accompany numerical simulations was discussed. An exact, closed-form analytical solution to the enhanced lumped capacitance model is derived. The temperature evolution in a representative 2D turbulated surface is simulated using Fluent to validate the model and its exact solution. A case including an interface contact resistance was included as well as various rib sizes to test the validity of the model over a range of conditions. The analysis was extended to the inter-rib region to investigate the extent and magnitude of the influence of the metallic rib features on the apparent heat transfer coefficients in the inter-rib region. It was found that the thermal contamination is limited only to the regions closest to the base of the rib feature. An experimental setup was developed, capable of measuring the local heat transfer distributions on all four channel walls of a rectangular channel (with aspect ratios between 1 and 5) at Reynolds numbers up to 150,000. The setup utilizes a transient thermochromic liquid crystals technique using narrow band crystals and a four camera setup. The setup is used to test a square channel with ribs applied to one wall. Using the transient thermochromic liquid crystals technique and applying it underneath high conductivity, metallic surface features, it is possible to calculate the heat transfer coefficient using a lumped heat capacitance approach. The enhanced lumped capacitance model is used to account for heat conduction into the substrate material. Rohacell and aluminum ribs adhered to the surface were used to tandem to validate the hybrid technique against the standard technique. Local data was also used to investigate the effect of thermal contamination. Thermal contamination observed empirically was more optimistic than numerical predictions. Traditional transient thermochromic liquid crystals technique utilizes the time-to-arrival of the peak intensity of the green color signal. The technique has been extended to utilize both the red and green color signals, increasing the throughput by recovering unused data while also allowing for a reduction in the experimental uncertainty of the calculated heat transfer coefficient. The over-determined system was solved using an un-weighted least squares approach. Uncertainty analysis of the multi-color technique demonstrated its superior performance over the single-color technique. The multi-color technique has the advantage of improved experimental uncertainty while being easy to implement.

The company Kanthal produces electric heating elements that require high temperature treatment in one production step. In this process step, called sintering, the amount of heat received by the sintered material is in direct correlation to the product’s outcome.  It is therefore of interest for the company to gather information about how heat transfer happens in an electrical furnace. This study examines two different possible scenarios of how the heat transfer in the furnace could look like and which amount of heat the sintered material would receive. The relation between a gaseous ambience at a certain temperature and the temperature an object submerged into this ambience is assuming is studied in the process called "transient heat conduction".  Two models were built in Matlab, representing transient heat conduction effects on two different geometries: a plane wall and a short cylinder.  It could be shown that transient heat conduction effects turned out differently for the two models. The conclusion drawn from the results was that the wall model was susceptible to horizontal heat transfer effects, whereas the cylinder model was affected from all directions equally. Further, an analysis of the heat transfer channels within the furnace revealed that the heat leakage through the furnace muffle edges, which are in contact with air, causes a multiple in heat loss compared to the overall heat leakage.

Heat flux and heat conductance at the metal mold interface plays a key role in controlling the final metal casting strength. It is difficult to obtain these parameters through direct measurement because of the required placement of sensors, however they can be obtained through inverse heat conduction calculations. Existing inverse heat conduction methods are analyzed and classified into three categories, i.e., direct inverse methods, observer-based methods and optimization methods. The solution of the direct inverse methods is based on the linear relationship between heat flux and temperature (either in the time domain or in the frequency domain) and is calculated in batch mode. The observer-based method consists on the application of observer theory to the inverse heat conduction problem. The prominent characteristic in this category is online estimation, but the methods in this category show weak robustness. Transforming estimation problems into optimization problems forms the methods in the third category. The methods in third category show very good robustness property and can be easily extended to multidimensional and nonlinear problems. The unknown parameters in some inverse heat conduction methods can be obtained by a proposed calibration procedure. A two-index property evaluation (accuracy and robustness) is also proposed to evaluate inverse heat conduction methods and thus determine which method is suitable for a given situation. The thermocouple dynamics effect on inverse calculation is also analyzed. If the thermocouple dynamics is omitted in the inverse calculation, the time constant of thermocouple should be as small as possible. Finally, a simple model is provided simulating the temperature measurement using a thermocouple. FEA (Finite Element Analysis) is employed to simulate temperature measurement.

No description available.

No description available.

The Cone Calorimeter has been used widely for various purposes as a bench - scale apparatus. Originally the retainer frame (edge frame) was designed to reduce unrepresentative edge burning of specimens. In general, the frame has been used in most Cone tests without enough understanding of its effect. It is very important to have one - dimensional (1D) conditions in order to estimate thermal properties of materials. It has been implicitly assumed that the heat conduction in the Cone Calorimeter is 1D using the current specimen preparation. However, the assumption has not been corroborated explicitly to date. The first objective of this study was to evaluate the heat transfer behavior of a Cone specimen by examining its three - dimensional (3D) heat conduction. It is essential to understand the role of wall lining materials when they are exposed to a fire from an ignition source. Full - scale test methods permit an assessment of the performance of a wall lining material. Fire growth models have been developed due to the costly expense associated with full - scale testing. The models require heat flux maps from the ignition burner flame as input data. Work to date was impeded by a lack of detailed spatial characterization of the heat flux maps due to the use of limited instrumentation. To increase the power of fire modeling, accurate and detailed heat flux maps from the ignition burner are essential. High level spatial resolution for surface temperature can be provided from an infrared camera. The second objective of this study was to develop a heat flux mapping procedure for a room test burner flame to a wall configuration with surface temperature information taken from an infrared camera. A prototype experiment is performed using the ISO 9705 test burner to demonstrate the developed heat flux mapping procedure. The results of the experiment allow the heat flux and spatial resolutions of the method to be determined and compared to the methods currently available.

Descreve-se neste trabalho um estudo experimental da Emissão Eletrociclotrônica (EEC) no tokamak TCABR. Um radiômetro de EEC foi instalado, calibrado e utilizado para o estudo do plasma em descargas térmicas. O radiômetro é do tipo heteródino de varredura operando na faixa de 50-85 GHz, no modo extraordinário e na segunda harmônica. Determinou-se a temperatura de ruído do radiômetro e também sua estabilidade em amplitude e freqüência. Foi medida a largura de banda em freqüência do radiômetro (resolução espacial horizontal). A antena utilizada é do tipo gaussiana sendo que o raio da cintura do feixe gaussiano e a posição do foco foram experimentalmente determinados (W0 ~ 1 cm e d ~ 37 cm, respectivamente). A posição da cintura da antena foi posicionada próxima do centro geométrico do vaso. Foi feita a calibração absoluta do equipamento considerando-se os efeitos das janelas de diagnóstico (reflexão e absorção). O sistema pode operar em modo varredura, para a obtenção de perfis radiais de Te ou modo freqüência única onde se tem alta resolução temporal. As medidas da radiação EC foram feitas, na sua maioria, em descargas com densidade eletrônica média entre 1.10+19 m-3 e 1,5.1019 m-3 de forma a se ter acessibilidade da radiação EC e também minimizar-se a presença de elétrons fugitivos. Para ne > 1,5.1019 m-3 (com B0 = 1,14 T) verifica-se o corte parcial da radiação EC. Nesta condição, o corte na EEC foi utilizado na determinação do perfil radial da densidade eletrônica e aplicado em três diferentes situações: descargas com injeção adicional de gás, com a aplicação do eletrodo de polarização e em descargas com injeção de ondas de radiofreqüência na região das ondas de Alfvén para o TCABR. Usando um perfil parabólico típico para a densidade eletrônica, observou-se que, para descargas com injeção adicional de gás ou em descargas com a aplicação de ondas de radiofreqüência tem-se 0,85 < alfa < 1, e para descargas com aplicação do eletrodo de polarização obteve-se alfa ~ 0,6. Foram feitas observações simultâneas da temperatura eletrônica, a partir do sinal da EEC, e das oscilações de Mirnov (freqüência ~ 11,7 kHz) em descargas térmicas com q(r=0) > 1. Os resultados indicam a presença de um modo de ruptura dominante em rs ~ 9,5 cm com a largura da ilha magnética de W ~ 2,0-2,5 cm. Estes resultados experimentais obtidos estão em acordo com os resultados indicados por teorias de transporte na região das ilhas magnéticas. Observou-se também que a localização da ilha magnética coincide com uma região onde o perfil radial da temperatura de plasma é aproximadamente plano. Num outro cenário, com q(r=0) < 1, observaram-se oscilações dente de serra com período de ~ 0,44 ms, tempo de queda de ~ 0,12 ms, e raio de inversão em r ~ 4 cm. Neste tipo de descargas observou-se que, no perfil radial da amplitude das oscilações da Te, |DeltaTe|, devido à propagação dos pulsos dente de serra, apresentavam posições de mínimos e que estes coincidiam com as posições onde ocorrem patamares no perfil radial da temperatura eletrônica. Partindo destes resultados, juntamente com o perfil de q(r), dão-nos os modos racionais (m/n), posições (r) e larguras(W) para as ilhas magnéticas, presentes nestas descargas, a saber: m/n = 4/3 (r ~ 9 cm, W4/3 ~ 0,9 cm), m/n = 3/2 (r ~ 11,8 cm, W3/2 ~ 0,9 cm) e m/n = 2/1 (r ~ 13,7 cm, W2/1 ~ 1,4 cm). Este novo método, aqui proposto, permite a determinação direta da posição e da largura das ilhas magnéticas, em descargas onde a instabilidade dente de serra encontra-se presente.¶ / In this work, an experimental study of the Electron Cyclotron Emission (ECE) in the TCABR Tokamak, is described. An ECE sweeping heterodyne radiometer, operating in the 50-85 GHz frequency range, was installed, calibrated and used to detect ECE radiation in the second harmonic extraordinary mode from thermal plasma discharges. The noise temperature, amplitude and frequency stability of the radiometer were determined. The frequency bandwidth (horizontal space resolution) was measured. A Gaussian Antenna is used and the gaussian beam waist radius (W0 ~ 1 cm) and the focus position (d ~ 37 cm) were experimentally determined. The focus of the antenna was positioned close to the center of the vacuum vessel. The absolute calibration of the equipment was done and the effect of the diagnostic window was considered (reflections and absorptions). The system can operate in sweeping mode, in order to obtain the radial electron temperature profiles, and also in the single frequency mode to obtain temporal electron temperature profiles with high time resolution. Due to the TCABR accessibility conditions and also to avoid runaway electrons, most of the ECE measurements were obtained in discharges with line electron density between 1.0x10+19 m-3 and 1.5x1019 m-3. For ne > 1.5x1019 m-3 (with B0 = 1.14 T) the cutoff in the ECE radiation was observed. The ECE cutoff was used to determine the radial profile of the electron density and applied to three different situations: discharges with additional gas puffing, with the application of a biasing electrode and in discharges with radio-frequency wave injection. Using a parabolic profile it was observed that, 0.85 < alfa < 1.0 for discharges with additional gas injection or with RF injection, and alfa ~ 0.6 for the electrode biasing experiments. The electron temperature profiles and Mirnov oscillations (f ~ 11.7 kHz) were simultaneously measured in discharges with q (r=0) > 1. The position and the width of the magnetic island were therefore calculated. The results indicate the presence of a dominant tearing mode in rs ~ 9.5 cm and the magnetic island width W ~ 2.0-2.5 cm. These experimental results are in agreement with the theoretical results foreseen by models of heat transport in the presence of magnetic islands. It was also observed that, the position of the magnetic island coincided with the region where the temperature radial profile is approximately flat. In another plasma scenario, with q(r=0) < 1, sawteeth oscillations with a period ~ 0.44 ms, crash time ~ 0.12 ms, inversion radius r ~ 4 cm, were measured. Peculiarities in the Te radial profile associated with the amplitude of Te oscillations, |DeltaTe|, due to the sawtooth instability, were observed. It was verified that the minimum values of the |DeltaTe| profile coincided with the region where the temperature radial profile was approximately flat. These results together with the q(r) profile yielded the following mode numbers (m/n), positions (r) and widths (W) of the magnetic islands: m/n = 4/3 (r ~ 9 cm, W4/3 ~ 0.9 cm), m/n = 3/2 (r ~ 11.8 cm, W3/2 ~ 0.9 cm) and m/n = 2/1 (r ~ 13.7 cm, W2/1; ~ 1.4 cm). A novel method to determine the position and width of the magnetic islands, in discharges with the presence of sawtooth instability, is presented here.

O problema da condução de calor, envolvendo mudança de fase, foi resolvido para o caso de uma cavidade limitada por duas superfícies cilíndricas indefinidamente longas. As condições de contorno impostas consistem em manter a temperatura da superfície interna fixa e abaixo da temperatura de fusão do material que preenche a cavidade, enquanto que a temperatura da superfície externa é mantida fixa e acima da temperatura de fusão. Como condição inicial se fixou a temperatura de todo o material que preenche a cavidade no valor da temperatura da superfície externa. A solução obtida consiste em duas soluções da equação de condução de calor, uma escrita para o material solidificado e outra escrita para o material em estado líquido. As duas soluções são formalmente escritas em termos da posição da frente de mudança de fase, que é representada por uma superfície cilíndrica com raio em expansão dentro da cavidade. A posição dessa superfície é, a princípio, desconhecida e é calculada impondo o balanço de energia através da frente da mudança de fase. O balanço de energia é expresso por uma equação diferencial de primeira ordem, cuja solução numérica fornece a posição da frente como função do tempo. A substituição da posição da frente de mudança de fase em um instante particular, nas soluções da equação de condução de calor, fornece a temperatura nas duas fases naquele instante. A solução obtida é ilustrada através de exemplos numéricos. / The heat conduction problem, in the presence of a change of state, was solved for the case of an indefinitely long cylindrical layer cavity. As boundary conditions it is imposed that the internal surface of the cavity is maintained below the fusion temperature of the infilling substance and the external surface is kept above it. The solution, obtained in non-dimensional variables, consists in two closed form heat conduction equation solutions for the solidified and liquid regions, which formally depend of the, at first, unknown position of the phase change front. The energy balance through the phase change front furnishes the equation for time dependence of the front position, which is numerically solved. Substitution of the front position for a particular instant in the heat conduction equation solutions gives the temperature distribution inside the cavity at that moment. The solution is illustrated with numerical examples.

A têmpera dos aços envolve a austenitização de uma peça seguida por um resfriamento rápido para promover a formação de microestrutura martensítica. É necessário avaliar os meios de têmpera para manter o processo de têmpera sob controle. Os parâmetros mais importantes no processo de resfriamento são o coeficiente de transferência de calor e/ou o fluxo de calor entre o meio de têmpera e a peça a ser resfriada. Um dos métodos de se avaliar os meios de têmpera (meios de resfriamento) e saber o que está acontecendo dentro da peça durante o resfriamento do ponto de vista térmico é o problema inverso de condução de calor. O problema inverso de condução de calor consiste na determinação de parâmetros como fluxo de calor, taxa de resfriamento e temperatura em qualquer posição através da peça, assim como o coeficiente de transferência de calor. Esses parâmetros são obtidos a partir de medições de temperatura em um ou mais pontos dentro da peça. O escopo deste trabalho foi desenvolver um software baseado no problema inverso condução de calor para avaliar meios de resfriamento para têmpera. A validação deste código foi feita usando água, óleo de soja, óleo mineral e solução aquosa de NaNO3. / Steels quenching involves part austenitization followed by a fast cooling to promote martensitic microstructure formation. It is necessary to evaluate quenchants in order to keep the quenching process under control. The most important cooling process parameters are the heat transfer coefficient and/or the heat flux between the quenchant and the part to be cooled. One of the methods to evaluate quenchants (cooling media) and to know what is happening inside the part during the cooling in the thermal point of view is the inverse heat conduction problem. The inverse heat conduction problem consists in the determination of parameters like heat flux, cooling rate and temperature in any position across the part, as well as the heat transfer coefficient. These parameters are obtained from temperature measurements in one or more points inside the part. The scope of this work was to develop a software based in the inverse heat conduction problem in order to evaluate quenchants for quenching. The validation of this code was made using water, soybean oil, mineral oil and NaNO3 aqueous solution.

Despite advancements in the development of synthetic fibers and materials that provide better insulation, fire ground burn injuries remain a significant issue. The current test methods for fire fighting clothing were investigated to determine their adequacy in evaluating the actual performance of clothing materials. This investigation uncovered several potential problems with the current test methods. A series of new, small scale, tests were used to evaluate the shortcomings of the current test methods and develop possible improvements. A small test apparatus, designed and donated by Ktech Corporation, was used to measure the thermal properties (thermal conductivity and volumetric heat capacity) of a series of fire fighting clothing materials. The thermal properties were estimated for single fabric layers, as well as ensembles, with various levels of moisture added to simulate actual end use conditions. In addition, a skin simulant sensor was used to assess the time to 2nd degree burn for exposures similar to those required in current standards for fire fighting clothing. A one dimensional heat conduction model was developed to predict the time to 2nd degree burn for the skin simulant sensor protected with outer shell materials that may be used as wildland fire fighting clothing, using the thermal property data obtained from earlier tests. An alternative method was developed to calculate the time to 2nd degree burn for ensembles evaluated with the new skin simulant sensor. The predictions for the time to 2nd degree burn obtained from the new skin simulant sensor were compared against results obtained using the sensor specified in the current test methods. The predictions for the skin simulant sensor were consistently shorter than those from the current test sensor. The current test sensor predictions for the time to 2nd degree burn were nominally 40% to 50% higher than the predictions from the skin simulant sensor during the evaluations of outer shell materials.