Simulação do processo de transferência de calor não linear condução-radiação por meio de um esquema linear em diferenças finitas.

Neste trabalho o processo não linear de transmissão de calor condução-radiação é abordado num contexto bidimensional plano e simulado com o uso de um esquema linear em diferenças finitas. O problema original é tratado como o limite de uma sequencia de problemas lineares, do tipo condução-convecção. Este limite, cuja existência é comprovada, é facilmente obtido a partir de procedimentos básicos, accessíveis a qualquer estudante de engenharia, permitindo assim o emprego de hipóteses mais realistas, já que não se tem o limitante matemático para a abordagem numérica de uma equação diferencial parcial elíptica. Neste trabalho foi resolvido o problema de condução de calor em regime permanente em uma placa com condições de contorno convectivas e radioativas utilizando-se o software MatLab, vale ressaltar, que a mesma metodologia é aplicável para geometrias mais complexas.

Transient electrothermal simulation of power semiconductor devices

In this paper, a new thermal model based on the Fourier series solution of heat conduction equation has been introduced in detail. 1-D and 2-D Fourier series thermal models have been programmed in MATLAB/Simulink. Compared with the traditional finite-difference thermal model and equivalent RC thermal network, the new thermal model can provide high simulation speed with high accuracy, which has been proved to be more favorable in dynamic thermal characterization on power semiconductor switches. The complete electrothermal simulation models of insulated gate bipolar transistor (IGBT) and power diodes under inductive load switching condition have been successfully implemented in MATLAB/Simulink. The experimental results on IGBT and power diodes with clamped inductive load switching tests have verified the new electrothermal simulation model. The advantage of Fourier series thermal model over widely used equivalent RC thermal network in dynamic thermal characterization has also been validated by the measured junction temperature.© 2010 IEEE.

Thermal transport in out-of-equilibrium quantum harmonic chains

We address the problem of heat transport in a chain of coupled quantum harmonic oscillators, exposed to the influences of local environments of various nature, stressing the effects that the specific nature of the environment has on the phenomenology of the transport process. We study in detail the behavior of thermodynamically relevant quantities such as heat currents and mean energies of the oscillators, establishing rigorous analytical conditions for the existence of a steady state, whose features we analyze carefully. In particular, we assess the conditions that should be faced to recover trends reminiscent of the classical Fourier law of heat conduction and highlight how such a possibility depends on the environment linked to our system.

Nanofluids development and characterization for heat exchanging intensification

A desmaterialização da economia é um dos caminhos para a promoção do desenvolvimento sustentável na medida em que elimina ou reduz a utilização de recursos naturais, fazendo mais com menos. A intensificação dos processos tecnológicos é uma forma de desmaterializar a economia. Sistemas mais compactos e mais eficientes consomem menos recursos. No caso concreto dos sistemas envolvendo processo de troca de calor, a intensificação resulta na redução da área de permuta e da quantidade de fluido de trabalho, o que para além de outra vantagem que possa apresentar decorrentes da miniaturização, é um contributo inegável para a sustentabilidade da sociedade através do desenvolvimento científico e tecnológico. O desenvolvimento de nanofluidos surge no sentido de dar resposta a estes tipo de desafios da sociedade moderna, contribuindo para a inovação de produtos e sistemas, dando resposta a problemas colocados ao nível das ciências de base. A literatura é unânime na identificação do seu potencial como fluidos de permuta, dada a sua elevada condutividade, no entanto a falta de rigor subjacente às técnicas de preparação dos mesmos, assim como de um conhecimento sistemático das suas propriedades físicas suportado por modelos físico-matemáticos devidamente validados levam a que a operacionalização industrial esteja longe de ser concretizável. Neste trabalho, estudou-se de forma sistemática a condutividade térmica de nanofluidos de base aquosa aditivados com nanotubos de carbono, tendo em vista a identificação dos mecanismos físicos responsáveis pela condução de calor no fluido e o desenvolvimento de um modelo geral que permita com segurança determinar esta propriedade com o rigor requerido ao nível da engenharia. Para o efeito apresentam-se métodos para uma preparação rigorosa e reprodutível deste tipo de nanofluido assim como das metodologias consideradas mais importantes para a aferição da sua estabilidade, assegurando deste modo o rigor da técnica da sua produção. A estabilidade coloidal é estabelecida de forma rigorosa tendo em conta parâmetros quantificáveis como a ausência de aglomeração, a separação de fases e a deterioração da morfologia das nanopartículas. Uma vez assegurado o método de preparação dos nanofluídos, realizou-se uma análise paramétrica conducente a uma base de dados obtidos experimentalmente que inclui a visão central e globalizante da influência relativa dos diferentes fatores de controlo com impacto nas propriedades termofísicas. De entre as propriedades termofísicas, este estudo deu particular ênfase à condutividade térmica, sendo os fatores de controlo selecionados os seguintes: fluido base, temperatura, tamanho da partícula e concentração de nanopartículas. Experimentalmente, verificou-se que de entre os fatores de controlo estudados, os que maior influência detêm sobre a condutividade térmica do nanofluido, são o tamanho e concentração das nanopartículas. Com a segurança conferida por uma base de dados sólida e com o conhecimento acerca da contribuição relativa de cada fator de controlo no processo de transferência de calor, desenvolveu-se e validou-se um modelo físico-matemático com um caracter generalista, que permitirá determinar com segurança a condutividade térmica de nanofluidos.

Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

LOSS of EXPONENTIAL STABILITY FOR A THERMOELASTIC SYSTEM WITH MEMORY

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Sensible and latent heat loss from the body surface of Holstein cows in a tropical environment

The general principles of the mechanisms of heat transfer are well known, but knowledge of the transition between evaporative and non-evaporative heat loss by Holstein cows in field conditions must be improved, especially for low-latitude environments. With this aim 15 Holstein cows managed in open pasture were observed in a tropical region. The latent heat loss from the body surface of the animals was measured by means of a ventilated capsule, while convective heat transfer was estimated by the theory of convection from a horizontal cylinder and by the long-wave radiation exchange based on the Stefan-Boltzmann law. When the air temperature was between 10 and 36 degrees C the sensible heat transfer varied from 160 to -30 W m(-2), while the latent heat loss by cutaneous evaporation increased from 30 to 350 W m(-2). Heat loss by cutaneous evaporation accounted for 20-30% of the total heat loss when air temperatures ranged from 10 to 20 degrees C. At air temperatures > 30 degrees C cutaneous evaporation becomes the main avenue of heat loss, accounting for approximately 85% of the total heat loss, while the rest is lost by respiratory evaporation.

A two phase flow nonlinear parameter estimation for capillary tube-suction line heat exchanger

Here we present two-phase flow nonlinear parameter estimation for HFC's flow through capillary tube-suction line heat exchangers, commonly used as expansion devices in small refrigeration systems. The simplifying assumptions adopted are: steady state, pure refrigerant, one-dimensional flow, negligible axial heat conduction in the fluid, capillary tube and suction line walls. Additionally, it is considered that the refrigerant is free from oil and both phases are assumed to be at the same pressure, that is, surface tension effects are neglected. Metastable flow effects are also disregarded, and the vapor is assumed to be saturated at the local pressure. The so-called homogeneous model, involving three, first order, ordinary differential equations is applied to analyze the two-phase flow region. Comparison is done with experimental measurements of the mass flow rate and temperature distribution along capillary tubes working with refrigerant HFC-134a in different operating conditions.

Termodinâmica de um conjunto de partículas em um bilhar bidimensional dependente do tempo: um gás bidimensional simplificado

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Experimental investigations of heat transfer coefficients of cutting fluids in metal cutting processes: analysis of workpiece phenomena in a given case study

This paper reports an experimental method to estimate the convective heat transfer of cutting fluids in a laminar flow regime applied on a thin steel plate. The heat source provided by the metal cutting was simulated by electrical heating of the plate. Three different cooling conditions were evaluated: a dry cooling system, a flooded cooling system and a minimum quantity of lubrication cooling system, as well as two different cutting fluids for the last two systems. The results showed considerable enhancement of convective heat transfer using the flooded system. For the dry and minimum quantity of lubrication systems, the heat conduction inside the body was much faster than the heat convection away from its surface. In addition, using the Biot number, the possible models were analyzed for conduction heat problems for each experimental condition tested.

Simulation of Temperatures and Stresses during Construction of a RCC Dam

Se describen las variaciones de temperaturas y de tensiones durante la construcción de presas de hormigón compactado. The curing of concrete is an exothermic process. The heat of hydration generated induces temperature increases in the concrete, which will disappear in the long term by heat conduction in the concrete mass and thermal exchanges with the environment. The problem is of particularly interest for large concrete masses, as is the case of dams, because the time involved in the heat diffusion process grows with the square of the dimensions and a hotter dam interior implies the possibility of cracking the exposed surfaces of the dam. The Cuira dam, currently being built in Venezuela using roller compacted concrete, is a 134 m high, arch-gravity dam. In support of the design, different strategies were analysed, including various combinations of cooling of the water and the aggregate in order to achieve acceptable results. The calculations were conducted with Abaqus, taking into account all the necessary mechanical and thermal characteristics, as well as the relevant non-linearities. The analyses led to the conclusion that no cooling was required, even taking into account the stress state imposed by an early and rapid filling of the reservoir.

Transition from isentropic to isothermal expansion in laser produced plasma

The transition that the expansion flow of laser-produced plasmas experiences when one moves from long, low intensity pulses (temperature vanishing at the isentropic plasma-vacuum front,lying at finite distance) to short, intense ones (non-zero, uniform temperature at the plasma-vacuum front, lying at infinity) is studied. For plznar geometry and lqge ion number Z, the transition occurs for dq5/dt=0.14(27/8)k712Z’1zn$/m4f, 12nK,,; mi, and K are laser intensity, critical density,ion mass, and Spitzer’s heat conduction coefficient. This result remains valid for finite Zit,h ough the numerical factor in d$/dt is different. Shorter wavelength lasers and higher 4 plasmas allow faster rising pulses below transition.

Effects of radiation losses on very lean methane/air flames propagating upward in a vertical tube

The stationary upward propagation of a very lean methane/air flame in a long vertical tube open at the bottom and closed at the top is simulated numerically using a single overall chemical reaction to model combustion and assuming an optically thin gas and a transparent or non-reflecting tube wall to approximately account for radiation losses from CO2CO2 and H2OH2O. Buoyancy plays a dominant role in the propagation of these flames and causes a large region of low velocity of the burnt gas relative to the flame to appear below the flame front when the equivalence ratio is decreased. The size of this region scales with the radius of the tube, and its presence enhances the effect of radiation losses, which would be otherwise negligible for a standard flammability tube, given the small concentration of radiating species. Heat conduction is found to be important in the low velocity region and to lead to a conduction flux from the flame to the burnt gas that causes extinction at the flame tip for a value of the equivalence ratio near the flammability limit experimentally measured in the standard tube. The effect of radiation losses decreases with the radius of the tube. Numerical results and order-of-magnitude estimates show that, in the absence of radiation, a very lean flame front fails to propagate only after recirculation of the burnt gas extends to its reaction region and drastically changes its structure. This condition is not realized for the standard flammability tube, but it seems to account for the flammability limit measured in a tube of about half the radius of the standard tube.

Modelling of transient heat transfer in annealing furnaces

This is a study of heat transfer in a lift-off furnace which is employed in the batch annealing of a stack of coils of steel strip. The objective of the project is to investigate the various factors which govern the furnace design and the heat transfer resistances, so as to reduce the time of the annealing cycle, and hence minimize the operating costs. The work involved mathematical modelling of patterns of gas flow and modes of heat transfer. These models are: Heat conduction and its conjectures in the steel coils;Convective heat transfer in the plates separating the coils in the stack and in other parts of the furnace; and Radiative and convective heat transfer in the furnace by using the long furnace model. An important part of the project is the development of numerical methods and computations to solve the transient models. A limited number of temperature measurements was available from experiments on a test coil in an industrial furnace. The mathematical model agreed well with these data. The model has been used to show the following characteristics of annealing furnaces, and to suggest further developments which would lead to significant savings: - The location of the limiting temperature in a coil is nearer to the hollow core than to the outer periphery. - Thermal expansion of the steel tends to open the coils, reduces their thermal conductivity in the radial direction, and hence prolongs the annealing cycle. Increasing the tension in the coils and/or heating from the core would overcome this heat transfer resistance. - The shape and dimensions of the convective channels in the plates have significant effect on heat convection in the stack. An optimal design of a channel is shown to be of a width-to-height ratio equal to 9. - Increasing the cooling rate, by using a fluidized bed instead of the normal shell and tube exchanger, would shorten the cooling time by about 15%, but increase the temperature differential in the stack. - For a specific charge weight, a stack of different-sized coils will have a shorter annealing cycle than one of equally-sized coils, provided that production constraints allow the stacking order to be optimal. - Recycle of hot flue gases to the firing zone of the furnace would produce a. decrease in the thermal efficiency up to 30% but decreases the heating time by about 26%.

Mechanisms of drying of particulate slurries

The literature on the evaporation of pure liquid drops and the drying of drops of solutions and slurries has been reviewed with particular reference to spray drying. A 0.1-0.2 mm glass filament-thermocouple was constructed and used to study simultaneously, heat and mass transfer from a single suspended drop placed in a humidity and temperature controlled, 28 mm OD vertical wind tunnel. Heat conduction through the filament was minimised eg at 100¦C it accounted for only 9.3% of the total heat transferred to a drop. Evaporation of single water drops was also studied in a 101 mm OD vertical wind tunnel. The Nusselt number was found to be a function of the Reynolds, Prandtl and Transfer number over an air temperature range of 17¦C to 107¦C. The proposed correlation is: Nu = 2+(-12.96B+0.76)Re¦-5Pr0-33 Experimental drying studies were carried out on single suspended 1 to 2.5 mm diameter drops of aqueous sodium sulphate decahydrate, sodium chloride, potassium sulphate, copper sulphate and sodium acetate solutions and slurries at temperatures of 20¦C to 124¦C. Dried crusts were examined by Scanning Electron Microscopy. The drying history of any material depended upon the nature of the crust formed. Sodium acetate formed a non-rigid skin prior to the formation of a rigid crust. A modified receding evaporation interface model was proposed for the drying of solutions and slurries. This covered both the constant rate period prior to crust formation and the subsequent falling rate period. The model was solved numerically for the variation in core temperature, drop weight and crust thickness. Good agreement was obtained between model predictions and experimental results for materials forming rigid crusts i.e. sodium sulphate decahydrate, sodium chloride, potassium sulphate and copper sulphate. However, the drying histories of drops of 10-20% weight initial concentration sodium acetate were unpredictable since formation of a non-rigid skin deviated from the model assumption of a rigid outer surface. At higher initial concentrations (40% weight) where a rigid crust was formed for sodium acetate, good agreement was obtained between experimental results and model predictions. Single suspended drop studies are concluded to provide a valuable insight into the drying mechanisms of specific solutions and slurries.