Estudo do consumo energético associado ao processo de cozedura na indústria cerâmica

O presente trabalho realizou-se no âmbito do Projeto Galp 20-2020 implementado na empresa Sanindusa, em parceira com a Universidade de Aveiro. Como principal objetivo pretendia-se estudar o consumo de energia durante o processo de cozedura das peças cerâmicas com o intuito de avaliar a viabilidade da substituição do refratário existente na empresa. Uma vez que esta medida implicaria um tempo de retorno do investimento demasiado elevado, optou-se por estudar o consumo de energia associado a uma placa horizontal, de massa unitária, de qualquer tipo de material refratário. Para tal, foi desenvolvido um Modelo Teórico capaz de calcular o calor absorvido por uma placa de refratário ao longo de todo o percurso dentro do Forno Túnel 2, existente na empresa e utilizado no processo de produção das peças cerâmicas. Através deste estudo foi possível concluir que o Modelo Teórico é bastante útil na medida em que permite a criação de vários cenários, através da alteração de diversas variáveis, permitindo conhecer qual o impacto de cada uma no consumo de energia neste equipamento. Como trabalho futuro propõe-se o desenvolvimento do Modelo Teórico apresentado para o estudo da transferência de calor bidimensional. Esta melhoria permitiria analisar o consumo de energia associado a diferentes formas de material refratário utilizado no processo de cozedura de peças cerâmicas (o que não é possível quando se aplica o conceito de transferência de calor unidimensional).

Moisture Desorption Studies on Polymer Hydrated and Vacuum Extruded Bentonite Clay Mat

Moisture desorption observations from two bentonite clay mats subjected to ten environmental zones with individually different combinations of laboratory-controlled constant temperatures (between 20 °C and 40 °C) and relative humidity (between 15% and 70%) are presented. These laboratory observations are compared with predictions from mathematical models, such as thin-layer drying equations and kinetic drying models proposed by Page, Wang and Singh, and Henderson and Pabis. The quality of fit of these models is assessed using standard error (SE) of estimate, relative percent of error, and coefficient of correlation. The Page model was found to better predict the drying kinetics of the bentonite clay mats for the simulated tropical climates. Critical study on the drying constant and moisture diffusion coefficient helps to assess the efficacy of a polymer to retain moisture and control desorption through water molecule bonding. This is further substantiated with the Guggenheim–Anderson–De Boer (GAB) desorption isotherm model which is presented.

A three dimensional finite volume approach to the thermo-mechanical modelling of the shape casting of metals

This paper presents a three dimensional, thermos-mechanical modelling approach to the cooling and solidification phases associated with the shape casting of metals ei. Die, sand and investment casting. Novel vortex-based Finite Volume (FV) methods are described and employed with regard to the small strain, non-linear Computational Solid Mechanics (CSM) capabilities required to model shape casting. The CSM capabilities include the non-linear material phenomena of creep and thermo-elasto-visco-plasticity at high temperatures and thermo-elasto-visco-plasticity at low temperatures and also multi body deformable contact with which can occur between the metal casting of the mould. The vortex-based FV methods, which can be readily applied to unstructured meshes, are included within a comprehensive FV modelling framework, PHYSICA. The additional heat transfer, by conduction and convection, filling, porosity and solidification algorithms existing within PHYSICA for the complete modelling of all shape casting process employ cell-centred FV methods. The termo-mechanical coupling is performed in a staggered incremental fashion, which addresses the possible gap formation between the component and the mould, and is ultimately validated against a variety of shape casting benchmarks.

Multi-physics modelling - A vital component of virtual manufacturing

One of the core tasks of the virtual-manufacturing environment is to characterise the transformation of the state of material during each of the unit processes. This transformation in shape, material properties, etc. can only be reliably achieved through the use of models in a simulation context. Unfortunately, many manufacturing processes involve the material being treated in both the liquid and solid state, the trans-formation of which may be achieved by heat transfer and/or electro-magnetic fields. The computational modelling of such processes, involving the interactions amongst various interacting phenomena, is a consider-able challenge. However, it must be addressed effectively if Virtual Manufacturing Environments are to become a reality! This contribution focuses upon one attempt to develop such a multi-physics computational toolkit. The approach uses a single discretisation procedure and provides for direct interaction amongst the component phenomena. The need to exploit parallel high performance hardware is addressed so that simulation elapsed times can be brought within the realms of practicality. Examples of Multiphysics modelling in relation to shape casting, and solder joint formation reinforce the motivation for this work.

Multi-physics modeling of aluminium reduction cells

Aluminium cells involve a range of complex physical processes which act simultaneously to provide a narrow satisfactory operating range. These processes involve electromagnetic fields, coupled with heat transfer and phase change, two phase fluid flow with a range of complexities plus the development of stress in the cell structure. All of these phenomena are coupled in some significant sense and so to provide a comprehensive model of these processes involves their representation simultaneously. Conventionally, aspects of the process have been modeled separately using uncoupled estimates of the effects of the other phenomena; this has enabled the use of standard commercial CFD and FEA tools. In this paper we will describe an approach to the modeling of aluminium cells which describes all the physics simultaneously. This approach uses a finite volume approximation for each of the phenomena and facilitates their interactions directly in the modeling-the complex geometries involved are addressed by using unstructured meshes. The very challenging issues to be overcome in this venture will be outlined and some preliminary results will be shown.

The shallow water approximation applied to the aluminium electrolysis process

The industrial production of aluminium is an electrolysis process where two superposed horizontal liquid layers are subjected to a mainly vertical electric current supplied by carbon electrodes. The lower layer consists of molten aluminium and lies on the cathode. The upper layer is the electrolyte and is covered by the anode. The interface between the two layers is often perturbed, leading to oscillations, or waves, similar to the waves on the surface of seas or lakes. The presence of electric currents and the resulting magnetic field are responsible for electromagnetic (Lorentz) forces within the fluid, which can amplify these oscillations and have an adverse influence on the process. The electrolytic bath vertical to horizontal aspect ratio is such, that it is advantageous to use the shallow water equations to model the interface motion. These are the depth-averaging the Navier-Stokes equations so that nonlinear and dispersion terms may be taken into account. Although these terms are essential to the prediction of wave dynamics, they are neglected in most of the literature on interface instabilities in aluminium reduction cells where only the linear theory is usually considered. The unknown variables are the two horizontal components of the fluid velocity, the height of the interface and the electric potential. In this application, a finite volume resolution of the double-layer shallow water equations including the electromagnetic sources has been developed, for incorporation into a generic three-dimensional computational fluid dynamics code that also deals with heat transfer within the cell.

A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Estudo do Sujamento da Superfície de Transferência de Calor do Evaporador de uma Bomba de Calor aquecido pelo efluente de uma fábrica de pasta para papel

O presente trabalho tem como principal objectivo a avaliação do sujamento da superfície de transferência de calor do evaporador de uma bomba de calor, quando se utiliza o efluente da fábrica de pasta e papel da Portucel Soporcel em Cacia como fluido de aquecimento. Para a realização deste estudo montou-se uma instalação, composta por uma bomba de calor, um circuito de água de arrefecimento do condensador da bomba de calor e um circuito de água residual quente para o aquecimento do evaporador da bomba de calor. O ensaio decorreu durante um período de 84 dias, durante o qual foram registadas as temperaturas em vários pontos dos circuitos e os caudais de circulação. A evolução temporal da quantidade de calor transferida no evaporador, assim como, a comparação do coeficiente global de transferência de calor em condições ideias e experimentais, permitiu verificar que a razão U/U0 para o evaporador da bomba de calor não teve um comportamento decrescente ao longo do período de ensaio, permitindo concluir que, nas condições operatórias da experiência, não ocorreu sujamento na superfície do evaporador. No estudo de corrosão do cobre quando exposto à água residual em questão, obteve-se, para a velocidade de corrosão, o valor de 1,56 mg/(dm2 .dia), indicando assim que a corrosão do cobre naquele meio é praticamente inexistente.

An algorithm for the inclusion of RC compact models of packages into board level thermal simulation tools

International audience

Quantifying uplift and denudation of a thermally heterogeneous crust: a detailed multi-thermochronometeric study of central west Britain

Topography is often thought as exclusively linked to mountain ranges formed by plates collision. It is now, however, known that apart from compression, uplift and denudation of rocks may be triggered by rifting, like it happens at elevated passive margins, and away from plate boundaries by both intra-plate stress causing reactivation of older structures, and by epeirogenic movements driven by mantle dynamics and initiating long-wavelength uplift. In the Cenozoic, central west Britain and other parts of the North Atlantic margins experienced multiple episodes of rock uplift and denudation that have been variable both at spatial and temporal scales. The origin of topography in central west Britain is enigmatic, and because of its location, it may be related to any of the processes mentioned above. In this study, three low temperature thermochronometers, the apatite fission track (AFT) and apatite and zircon (U-Th-Sm)/He (AHe and ZHe, respectively) methods were used to establish the rock cooling history from 200◦C to 30◦C. The samples were collected from the intrusive rocks in the high elevation, high relief regions of the Lake District (NW England), southern Scotland and northern Wales. AFT ages from the region are youngest (55–70 Ma) in the Lake District and increase northwards into southern Scotland and southwards in north Wales (>200 Ma). AHe and ZHe ages show no systematic pattern; the former range from 50 to 80 Ma and the latter tend to record the post-emplacement cooling of the intrusions (200–400 Ma). The complex, multi-thermochronometric inverse modelling suggests a ubiquitous, rapid Late Cretaceous/early Palaeogene cooling event that is particularly marked in Lake District and Criffell. The timing and rate of cooling in southern Scotland and in northern Wales is poorly resolved as the amount of cooling was less than 60◦C. The Lake District plutons were at >110◦C prior to the early Palaeogene; cooling due to a combined effect of high heat flow, from the heat producing granite batholith, and the blanketing effect of the overlying low conductivity Late Mesozoic limestones and mudstones. Modelling of the heat transfer suggests that this combination produced an elevated geothermal gradient within the sedimentary rocks (50–70◦C/km) that was about two times higher than at the present day. Inverse modelling of the AFT and AHe data taking the crustal structure into consideration suggests that denudation was the highest, 2.0–2.5 km, in the coastal areas of the Lake District and southern Scotland, gradually decreasing to less than 1 km in the northern Southern Uplands and northern Wales. Both the rift-related uplift and the intra-plate compression poorly correlate with the timing, location and spatial distribution of the early Palaeogene denudation. The pattern of early Palaeogene denudation correlates with the thickness of magmatic underplating, if the changes of mean topography, Late Cretaceous water depth and eroded rock density are taken into consideration. However, the uplift due to underplating alone cannot fully justify the total early Palaeogene denudation. The amount that is not ex- plained by underplating is, however, roughly spatially constant across the study area and can be referred to the transient thermal uplift induced by the mantle plume arrival. No other mechanisms are required to explain the observed pattern of denudation. The onset of denudation across the region is not uniform. Denudation started at 70–75 Ma in the central part of the Lake District whereas the coastal areas the rapid erosion appears to have initiated later (65–60 Ma). This is ~10 Ma earlier than the first vol- canic manifestation of the proto-Iceland plume and favours the hypothesis of the short period of plume incubation below the lithosphere before the volcanism. In most of the localities, the rocks had cooled to temperatures lower than 30◦C by the end of the Palaeogene, suggesting that the total Neogene denudation was, at a maximum, several hundreds of metres. Rapid cooling in the last 3 million years is resolved in some places in southern Scotland, where it could be explained by glacial erosion and post-glacial isostatic uplift.

COMPREHENSIVE ELECTRICAL/OPTICAL/THERMAL CHARACTERIZATIONS OF HIGH POWER LIGHT EMITTING DIODES AND LASER DIODES

Thermal characterizations of high power light emitting diodes (LEDs) and laser diodes (LDs) are one of the most critical issues to achieve optimal performance such as center wavelength, spectrum, power efficiency, and reliability. Unique electrical/optical/thermal characterizations are proposed to analyze the complex thermal issues of high power LEDs and LDs. First, an advanced inverse approach, based on the transient junction temperature behavior, is proposed and implemented to quantify the resistance of the die-attach thermal interface (DTI) in high power LEDs. A hybrid analytical/numerical model is utilized to determine an approximate transient junction temperature behavior, which is governed predominantly by the resistance of the DTI. Then, an accurate value of the resistance of the DTI is determined inversely from the experimental data over the predetermined transient time domain using numerical modeling. Secondly, the effect of junction temperature on heat dissipation of high power LEDs is investigated. The theoretical aspect of junction temperature dependency of two major parameters – the forward voltage and the radiant flux – on heat dissipation is reviewed. Actual measurements of the heat dissipation over a wide range of junction temperatures are followed to quantify the effect of the parameters using commercially available LEDs. An empirical model of heat dissipation is proposed for applications in practice. Finally, a hybrid experimental/numerical method is proposed to predict the junction temperature distribution of a high power LD bar. A commercial water-cooled LD bar is used to present the proposed method. A unique experimental setup is developed and implemented to measure the average junction temperatures of the LD bar. After measuring the heat dissipation of the LD bar, the effective heat transfer coefficient of the cooling system is determined inversely. The characterized properties are used to predict the junction temperature distribution over the LD bar under high operating currents. The results are presented in conjunction with the wall-plug efficiency and the center wavelength shift.

A geometric study of liquid retention in open-cell metal foams

Open-cell metal foams show promise as an emerging novel material for heat exchanger applications. The high surface-area-to-volume ratio suggests increased compactness and decrease in weight of heat exchanger designs. However, the metal foam structure appears conducive to condensate retention, which would degenerate heat transfer performance. This research investigates the condensate retention behavior of aluminum open-cell metal foams through the use of static dip tests and geometrical classification via X-ray Micro-Computed Tomography. Aluminum open-cell metal foam samples of 5, 10, 20, and 40 pores per inch (PPI), all having a void fraction greater than 90%, were included in this investigation. In order to model the condensate retention behavior of metal foams, a clearer understanding of the geometry was required. After exploring the ideal geometries presented in the open literature, X-ray Micro-Computed Tomography was employed to classify the actual geometry of the metal foam samples. The images obtained were analyzed using specialized software from which geometric information including strut length and pore shapes were extracted. The results discerned a high variability in ligament length, as well as features supporting the ideal geometry known as the Weaire-Phelan unit cell. The static dip tests consisted of submerging the metal foam samples in a liquid, then allowing gravity-induced drainage until steady-state was reached and the liquid remaining in the metal foam sample was measured. Three different liquids, water, ethylene glycol, and 91% isopropyl alcohol, were employed. The behaviors of untreated samples were compared to samples subjected to a Beomite surface treatment process, and no significant differences in retention behavior were discovered. The dip test results revealed two distinct regions of condensate retention, each holding approximately half of the total liquid retained by the sample. As expected, condensate retention increased as the pores sizes decreased. A model based on surface tension was developed to predict the condensate retention in the metal foam samples and verified using a regular mesh. Applying the model to both the ideal and actual metal foam geometries showed good agreement with the dip test results in this study.