Thermal aspects of satellite downscaling

The strong trend toward nanosatellites creates new challenges in terms of thermal balance control. The thermal balance of a satellite is determined by the heat dissipation in its subsystems and by the thermal connections between them. As satellites become smaller, heat dissipation in their subsystems tends to decrease and thermal connectivity scales down with dimension. However, these two terms do not necessarily scale in the same way, and so the thermal balance may alter and the temperature of subsystems may reach undesired levels. This paper focuses on low-Earth-orbit satellites. We constructed a generalized lumped thermal model that combines a generalized low-Earth-orbit satellite configuration with scaling trends in subsystem heat dissipation and thermal connectivity. Using satellite mass as a scaling parameter, we show that subsystems do not become thermally critical by scaling mass alone.

Temperature Dynamics Investigation At Small And Shallow Lakes Using Hydrodynamic Model

A three-dimensional time-dependent hydrodynamic and heat transport model of Lake Binaba, a shallow and small dam reservoir in Ghana, emphasizing the simulation of dynamics and thermal structure has been developed. Most numerical studies of temperature dynamics in reservoirs are based on one- or two-dimensional models. These models are not applicable for reservoirs characterized with complex flow pattern and unsteady heat exchange between the atmosphere and water surface. Continuity, momentum and temperature transport equations have been solved. Proper assignment of boundary conditions, especially surface heat fluxes, has been found crucial in simulating the lake’s hydrothermal dynamics. This model is based on the Reynolds Average Navier-Stokes equations, using a Boussinesq approach, with a standard k − ε turbulence closure to solve the flow field. The thermal model includes a heat source term, which takes into account the short wave radiation and also heat convection at the free surface, which is function of air temperatures, wind velocity and stability conditions of atmospheric boundary layer over the water surface. The governing equations of the model have been solved by OpenFOAM; an open source, freely available CFD toolbox. As its core, OpenFOAM has a set of efficient C++ modules that are used to build solvers. It uses collocated, polyhedral numerics that can be applied on unstructured meshes and can be easily extended to run in parallel. A new solver has been developed to solve the hydrothermal model of lake. The simulated temperature was compared against a 15 days field data set. Simulated and measured temperature profiles in the probe locations show reasonable agreement. The model might be able to compute total heat storage of water bodies to estimate evaporation from water surface.

Linear Approach to the Orbiting Spacecraft Thermal Problem

A linear method is developed for solving the nonlinear differential equations of a lumped-parameter thermal model of a spacecraft moving in a closed orbit. This method, based on perturbation theory, is compared with heuristic linearizations of the same equations. The essential feature of the linear approach is that it provides a decomposition in thermal modes, like the decomposition of mechanical vibrations in normal modes. The stationary periodic solution of the linear equations can be alternately expressed as an explicit integral or as a Fourier series. This method is applied to a minimal thermal model of a satellite with ten isothermal parts (nodes), and the method is compared with direct numerical integration of the nonlinear equations. The computational complexity of this method is briefly studied for general thermal models of orbiting spacecraft, and it is concluded that it is certainly useful for reduced models and conceptual design but it can also be more efficient than the direct integration of the equations for large models. The results of the Fourier series computations for the ten-node satellite model show that the periodic solution at the second perturbative order is sufficiently accurate.

Magnetic Attitude Control System for a Small Satellite. Impact on the Thermal Performance

El principal objetivo de la tesis es estudiar el acoplamiento entre los subsistemas de control de actitud y de control térmico de un pequeño satélite, con el fin de buscar la solución a los problemas relacionados con la determinación de los parámetros de diseño. Se considera la evolución de la actitud y de las temperaturas del satélite bajo la influencia de dos estrategias de orientación diferentes: 1) estabilización magnética pasiva de la orientación (PMAS, passive magnetic attitude stabilization), y 2) control de actitud magnético activo (AMAC, active magnetic attitude control). En primer lugar se presenta el modelo matemático del problema, que incluye la dinámica rotacional y el modelo térmico. En el problema térmico se considera un satélite cúbico modelizado por medio de siete nodos (seis externos y uno interno) aplicando la ecuación del balance térmico. Una vez establecido el modelo matemático del problema, se estudia la evolución que corresponde a las dos estrategias mencionadas. La estrategia PMAS se ha seleccionado por su simplicidad, fiabilidad, bajo coste, ahorrando consumo de potencia, masa coste y complejidad, comparado con otras estrategias. Se ha considerado otra estrategia de control que consigue que el satélite gire a una velocidad requerida alrededor de un eje deseado de giro, pudiendo controlar su dirección en un sistema inercial de referencia, ya que frecuentemente el subsistema térmico establece requisitos de giro alrededor de un eje del satélite orientado en una dirección perpendicular a la radiación solar incidente. En relación con el problema térmico, para estudiar la influencia de la velocidad de giro en la evolución de las temperaturas en diversos puntos del satélite, se ha empleado un modelo térmico linealizado, obtenido a partir de la formulación no lineal aplicando un método de perturbaciones. El resultado del estudio muestra que el tiempo de estabilización de la temperatura y la influencia de las cargas periódicas externas disminuye cuando aumenta la velocidad de giro. Los cambios de temperatura se reducen hasta ser muy pequeños para velocidades de rotación altas. En relación con la estrategia PMAC se ha observado que a pesar de su uso extendido entre los micro y nano satélites todavía presenta problemas que resolver. Estos problemas están relacionados con el dimensionamiento de los parámetros del sistema y la predicción del funcionamiento en órbita. Los problemas aparecen debido a la dificultad en la determinación de las características magnéticas de los cuerpos ferromagnéticos (varillas de histéresis) que se utilizan como amortiguadores de oscilaciones en los satélites. Para estudiar este problema se presenta un modelo analítico que permite estimar la eficiencia del amortiguamiento, y que se ha aplicado al estudio del comportamiento en vuelo de varios satélites, y que se ha empleado para comparar los resultados del modelo con los obtenidos en vuelo, observándose que el modelo permite explicar satisfactoriamente el comportamiento registrado. ABSTRACT The main objective of this thesis is to study the coupling between the attitude control and thermal control subsystems of a small satellite, and address the solution to some existing issues concerning the determination of their parameters. Through the thesis the attitude and temperature evolution of the satellite is studied under the influence of two independent attitude stabilization and control strategies: (1) passive magnetic attitude stabilization (PMAS), and (2) active magnetic attitude control (AMAC). In this regard the mathematical model of the problem is explained and presented. The mathematical model includes both the rotational dynamics and the thermal model. The thermal model is derived for a cubic satellite by solving the heat balance equation for 6 external and 1 internal nodes. Once established the mathematical model of the problem, the above mentioned attitude strategies were applied to the system and the temperature evolution of the 7 nodes of the satellite was studied. The PMAS technique has been selected to be studied due to its prevalent use, simplicity, reliability, and cost, as this strategy significantly saves the overall power, weight, cost, and reduces the complexity of the system compared to other attitude control strategies. In addition to that, another control law that provides the satellite with a desired spin rate along a desired axis of the satellite, whose direction can be controlled with respect to the inertial reference frame is considered, as the thermal subsystem of a satellite usually demands a spin requirement around an axis of the satellite which is positioned perpendicular to the direction of the coming solar radiation. Concerning the thermal problem, to study the influence of spin rate on temperature evolution of the satellite a linear approach of the thermal model is used, which is based on perturbation theory applied to the nonlinear differential equations of the thermal model of a spacecraft moving in a closed orbit. The results of this study showed that the temperature stabilization time and the periodic influence of the external thermal loads decreases by increasing the spin rate. However, the changes become insignificant for higher values of spin rate. Concerning the PMAS strategy, it was observed that in spite of its extended application to micro and nano satellites, still there are some issues to be solved regarding this strategy. These issues are related to the sizing of its system parameters and predicting the in-orbit performance. The problems were found to be rooted in the difficulties that exist in determining the magnetic characteristics of the ferromagnetic bodies (hysteresis rods) that are applied as damping devices on-board satellites. To address these issues an analytic model for estimating their damping efficiency is proposed and applied to several existing satellites in order to compare the results with their respective in-flight data. This model can explain the behavior showed by these satellites.

Modeling and validation of the thermal behavior of buildings for the development of demand response methods

The representation of the thermal behaviour of the building is achieved through a relatively simple dynamic model that takes into account the effects due to the thermal mass of the building components. The model of a intra-floor apartment has been built in the Matlab-Simulink environment and considers the heat transmission through the external envelope, wall and windows, the internal thermal masses, (i.e. furniture, internal wall and floor slabs) and the sun gain due to opaque and see-through surfaces of the external envelope. The simulations results for the entire year have been compared and the model validated, with the one obtained with the dynamic building simulation software Energyplus.

Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

Real-time temperature estimation for power MOSFETs considering thermal ageing effects

This paper presents a novel real-time power-device temperature estimation method that monitors the power MOSFET's junction temperature shift arising from thermal aging effects and incorporates the updated electrothermal models of power modules into digital controllers. Currently, the real-time estimator is emerging as an important tool for active control of device junction temperature as well as online health monitoring for power electronic systems, but its thermal model fails to address the device's ongoing degradation. Because of a mismatch of coefficients of thermal expansion between layers of power devices, repetitive thermal cycling will cause cracks, voids, and even delamination within the device components, particularly in the solder and thermal grease layers. Consequently, the thermal resistance of power devices will increase, making it possible to use thermal resistance (and junction temperature) as key indicators for condition monitoring and control purposes. In this paper, the predicted device temperature via threshold voltage measurements is compared with the real-time estimated ones, and the difference is attributed to the aging of the device. The thermal models in digital controllers are frequently updated to correct the shift caused by thermal aging effects. Experimental results on three power MOSFETs confirm that the proposed methodologies are effective to incorporate the thermal aging effects in the power-device temperature estimator with good accuracy. The developed adaptive technologies can be applied to other power devices such as IGBTs and SiC MOSFETs, and have significant economic implications.

Thermal Tracking and Estimation for Microprocessors

Due to increasing integration density and operating frequency of today's high performance processors, the temperature of a typical chip can easily exceed 100 degrees Celsius. However, the runtime thermal state of a chip is very hard to predict and manage due to the random nature in computing workloads, as well as the process, voltage and ambient temperature variability (together called PVT variability). The uneven nature (both in time and space) of the heat dissipation of the chip could lead to severe reliability issues and error-prone chip behavior (e.g. timing errors). Many dynamic power/thermal management techniques have been proposed to address this issue such as dynamic voltage and frequency scaling (DVFS), clock gating and etc. However, most of such techniques require accurate knowledge of the runtime thermal state of the chip to make efficient and effective control decisions. In this work we address the problem of tracking and managing the temperature of microprocessors which include the following sub-problems: (1) how to design an efficient sensor-based thermal tracking system on a given design that could provide accurate real-time temperature feedback; (2) what statistical techniques could be used to estimate the full-chip thermal profile based on very limited (and possibly noise-corrupted) sensor observations; (3) how do we adapt to changes in the underlying system's behavior, since such changes could impact the accuracy of our thermal estimation. The thermal tracking methodology proposed in this work is enabled by on-chip sensors which are already implemented in many modern processors. We first investigate the underlying relationship between heat distribution and power consumption, then we introduce an accurate thermal model for the chip system. Based on this model, we characterize the temperature correlation that exists among different chip modules and explore statistical approaches (such as those based on Kalman filter) that could utilize such correlation to estimate the accurate chip-level thermal profiles in real time. Such estimation is performed based on limited sensor information because sensors are usually resource constrained and noise-corrupted. We also took a further step to extend the standard Kalman filter approach to account for (1) nonlinear effects such as leakage-temperature interdependency and (2) varying statistical characteristics in the underlying system model. The proposed thermal tracking infrastructure and estimation algorithms could consistently generate accurate thermal estimates even when the system is switching among workloads that have very distinct characteristics. Through experiments, our approaches have demonstrated promising results with much higher accuracy compared to existing approaches. Such results can be used to ensure thermal reliability and improve the effectiveness of dynamic thermal management techniques.

Stability of planar flames as gasdynamic discontinuities

The stability of a steadily propagating planar premixed flame has been the subject of numerous studies since Darrieus and Landau showed that in their model flames are unstable to perturbations of any wavelength. Moreover, the instability was shown to persist even for very small wavelengths, i.e. there was no high-wavenumber cutoff of the instability. In addition to the Darrieus-Landau instability, which results from thermal expansion, analysis of the diffusional thermal model indicates that premixed flames may exhibit cellular and pulsating instabilities as a consequence of preferential diffusion. However, no previous theory captured all the instabilities including a high-wavenumber cutoff for each. In Class, Matkowsky & Klimenko (2003) a unified theory is proposed which, in appropriate limits and under appropriate assumptions, recovers all the relevant previous theories. It also includes additional new terms, not present in previous theories. In the present paper we consider the stability of a uniformly propagating planar flame as a solution of the unified model. The results are then compared to those based on the models of Darrieus-Landau, Sivashinsky and Matalon-Matkowsky. In particular, it is shown that the unified model is the only model to capture the Darrieus-Landau, cellular and pulsating instabilities including a high-wavenumber cutoff for each.

Caracterização experimental e numérica do desempenho térmico de um edifício "gaioleiro": identificação de oportunidades de melhoria

O sector dos edifícios é responsável por uma percentagem significativa dos consumos de energia primária e de energia eléctrica em Portugal, associada principalmente ao conforto térmico dos seus ocupantes. A União Europeia pretende uma redução de 20%, até 2020, dos consumos de energia e consequentes emissões de CO2 através da melhoria da eficiência energética dos edifícios públicos e residenciais. Em Portugal, o Plano Nacional para a Eficiência Energética (PNAEE) tem por objectivo obter uma poupança anual de energia de pelo menos 1% até ao ano de 2016, tomando como base a média de consumos de energia final, registados entre 2001 e 2005. Neste contexto, os edifícios anteriores a 1990 (primeira versão do RCCTE) podem apresentar um potencial significativo de melhoria da sua eficiência energética com base na sua reabilitação. Os edifícios “Gaioleiros” (1880 – 1930) constituindo uma parte importante do património histórico da cidade de Lisboa, para os quais a informação sobre o seu desempenho térmico é limitada, considerou-se pertinente efectuar um estudo destinado à sua caracterização experimental e numérica, face à especificidade do comportamento térmico das suas paredes caracterizadas pela elevada espessura. No presente trabalho, apresenta-se a metodologia e os resultados experimentais da medição da resistência térmica das paredes e da medição das necessidades térmicas de aquecimento da habitação. Estes resultados experimentais foram utilizados na validação do modelo de simulação térmica da habitação, que posteriormente serviu para avaliar as suas necessidades térmicas de aquecimento (Nic) e de arrefecimento (Nvc), identificar oportunidades de melhoria e avaliar o respectivo potencial de reabilitação. Neste trabalho, como contributos para uma reabilitação sustentável, apresentam-se avaliações de oportunidades de melhoria com base em estratégias de reforço do isolamento térmico. Dos resultados obtidos concluiu-se que melhorando o isolamento térmico das paredes e vãos envidraçados é possível baixar consideravelmente os consumos de energia associados à habitação, cumprindo assim as exigências estabelecidas no RCCTE ao nível dos requisitos de qualidade térmica da envolvente e consumos energéticos para edifícios novos e grandes reabilitações.

Predictive analytical and numerical modeling for orthogonal cutting

In this thesis, a predictive analytical and numerical modeling approach for the orthogonal cutting process is proposed to calculate temperature distributions and subsequently, forces and stress distributions. The models proposed include a constitutive model for the material being cut based on the work of Weber, a model for the shear plane based on Merchants model, a model describing the contribution of friction based on Zorev’s approach, a model for the effect of wear on the tool based on the work of Waldorf, and a thermal model based on the works of Komanduri and Hou, with a fraction heat partition for a non-uniform distribution of the heat in the interfaces, but extended to encompass a set of contributions to the global temperature rise of chip, tool and work piece. The models proposed in this work, try to avoid from experimental based values or expressions, and simplifying assumptions or suppositions, as much as possible. On a thermo-physical point of view, the results were affected not only by the mechanical or cutting parameters chosen, but also by their coupling effects, instead of the simplifying way of modeling which is to contemplate only the direct effect of the variation of a parameter. The implementation of these models was performed using the MATLAB environment. Since it was possible to find in the literature all the parameters for AISI 1045 and AISI O2, these materials were used to run the simulations in order to avoid arbitrary assumption.

Gas gap heat switch for CryoFree project

Cryogen-free superconducting magnet systems have become popular during the last two decades for the simple reason that with the use of liquid helium is rather cumbersome and is a scarce resource. Some available CFMS uses a mechanical cryocooler as cold source of the superconductor magnet. However, the cooling of the sample holder is still made through an open circuit of helium. A thermal management of a completely cryogen-free system is possible to be implemented by using a controlled gas gap heat switch (GGHS) between the cryocooler and the variable temperature insert (VTI). This way it would eliminate the helium open circuit. Heat switches are devices that allow to toggle between two distinct thermal states (ON and OFF state). Several cryogenic applications need good thermal contact and a good thermal insulation at different stages of operation. A versatile GGHS was designed and built with a 100 mm gap and tested with helium as exchange gas. An analytic thermal model was developed and a good agreement with the experimental data was obtained. The device was tested on a crycooler at 4 to 80 K ranges. A 285 mW/K thermal conductance was measured at ON state and 0.09 mW/K at OFF. 3000 ON/OFF thermal conductance ratio was obtained at 4 K with helium.

40 K Neon liquid energy storage unit

Cryocoolers have been progressively replacing the use of the stored cryogens in cryogenic chains used for detector cooling, thanks to their higher and higher reliability. However, the mechanical vibrations, the electromagnetic interferences and the temperature fluctuations inherent to their functioning could reduce the sensor’s sensitivity. In order to minimize this problem, compact thermal energy storage units (ESU) are studied, devices able to store thermal energy without significant temperature increase. These devices can be used as a temporary cold source making it possible to turn the cryocooler OFF providing a proper environment for the sensor. A heat switch is responsible for the thermal decoupling of the ESU from the cryocooler’s temperature that increases when turned OFF. In this work, several prototypes working around 40 K were designed, built and characterized. They consist in a low temperature cell that contains the liquid neon connected to an expansion volume at room temperature for gas storage during the liquid evaporation phase. To turn this system insensitive to the gravity direction, the liquid is retained in the low temperature cell by capillary effect in a porous material. Thanks to pressure regulation of the liquid neon bath, 900 J were stored at 40K. The higher latent heat of the liquid and the inexistence of triple point transitions at 40 K turn the pressure control during the evaporation a versatile and compact alternative to an ESU working at the triple point transitions. A quite compact second prototype ESU directly connected to the cryocooler cold finger was tested as a temperature stabilizer. This device was able to stabilize the cryocooler temperature ((≈ 40K ±1 K) despite sudden heat bursts corresponding to twice the cooling power of the cryocooler. This thesis describes the construction of these devices as well as the tests performed. It is also shown that the thermal model developed to predict the thermal behaviour of these devices, implemented as a software,describes quite well the experimental results. Solutions to improve these devices are also proposed.

40 K neon liquid energy storage unit

Cryocoolers have been progressively replacing the use of the stored cryogens in cryogenic chains used for detector cooling, thanks to their higher and higher reliability. However, the mechanical vibrations, the electromagnetic interferences and the temperature fluctuations inherent to their functioning could reduce the sensor’s sensitivity. In order to minimize this problem, compact thermal energy storage units (ESU) are studied, devices able to store thermal energy without significant temperature increase. These devices can be used as a temporary cold source making it possible to turn the cryocooler OFF providing a proper environment for the sensor. A heat switch is responsible for the thermal decoupling of the ESU from the cryocooler’s temperature that increases when turned OFF. In this work, several prototypes working around 40 K were designed, built and characterized. They consist in a low temperature cell that contains the liquid neon connected to an expansion volume at room temperature for gas storage during the liquid evaporation phase. To turn this system insensitive to the gravity direction, the liquid is retained in the low temperature cell by capillary effect in a porous material. Thanks to pressure regulation of the liquid neon bath, 900 J were stored at 40K. The higher latent heat of the liquid and the inexistence of triple point transitions at 40 K turn the pressure control during the evaporation a versatile and compact alternative to an ESU working at the triple point transitions. A quite compact second prototype ESU directly connected to the cryocooler cold finger was tested as a temperature stabilizer. This device was able to stabilize the cryocooler temperature ((≈ 40K ±1 K) despite sudden heat bursts corresponding to twice the cooling power of the cryocooler. This thesis describes the construction of these devices as well as the tests performed. It is also shown that the thermal model developed to predict the thermal behaviour of these devices,implemented as a software, describes quite well the experimental results. Solutions to improve these devices are also proposed.

Incorporation of hybrid phase change materials in plastering mortars for increased energy efficiency in buildings

Tese de Doutoramento em Engenharia Civil.