The use of primitives in the calculation of radiative view factors

Compilations of radiative view factors (often in closed analytical form) are readily available in the open literature for commonly encountered geometries. For more complex three-dimensional (3D) scenarios, however, the effort required to solve the requisite multi-dimensional integrations needed to estimate a required view factor can be daunting to say the least. In such cases, a combination of finite element methods (where the geometry in question is sub-divided into a large number of uniform, often triangular, elements) and Monte Carlo Ray Tracing (MC-RT) has been developed, although frequently the software implementation is suitable only for a limited set of geometrical scenarios. Driven initially by a need to calculate the radiative heat transfer occurring within an operational fibre-drawing furnace, this research set out to examine options whereby MC-RT could be used to cost-effectively calculate any generic 3D radiative view factor using current vectorisation technologies

Infrared thermochromic antenna composite for self-adaptive thermoregulation

Self-adaptive thermoregulation, the mechanism living organisms use to balance their temperature, holds great promise for decarbonizing cooling and heating processes. The functionality can be effectively emulated by engineering the thermal emissivity of materials to adapt to background temperature variations. Yet, solutions that marry large emissivity switching (${\Delta}\epsilon$) with scalability, cost-effectiveness and design freedom are still lacking. Here, we fill this gap by introducing infrared dipole antennas made of tunable thermochromic materials. We demonstrate that non-spherical antennas (rods, stars and flakes) made of vanadium-dioxide can exhibit a massive (~200-fold) increase in their absorption cross-section as temperature rises. Embedding these antennas in polymer films, or simply spraying them directly, creates free-form thermoregulation composites, featuring an outstanding ${\Delta}\epsilon\sim0.6$ in spectral ranges that can be tuned at will. Our research paves the way for versatile self-adaptive heat management solutions (coatings, fibers, membranes and films) that could find application in radiative-cooling, heat-sensing, thermal-camouflage, and other

Computational Methods in Science and Engineering : Proceedings of the Workshop SimLabs@KIT, November 29 - 30, 2010, Karlsruhe, Germany

In this proceedings volume we provide a compilation of article contributions equally covering applications from different research fields and ranging from capacity up to capability computing. Besides classical computing aspects such as parallelization, the focus of these proceedings is on multi-scale approaches and methods for tackling algorithm and data complexity. Also practical aspects regarding the usage of the HPC infrastructure and available tools and software at the SCC are presented

The use of high-performance computing to solve participating media radiative heat transfer problems-results of an NSF workshop

Radiation in participating media is an important transport mechanism in many physical systems. The simulation of complex radiative transfer has not effectively exploited high-performance computing capabilities. In response to this need, a workshop attended by members active in the high-performance computing community, members active in the radiative transfer community, and members from closely related fields was held to identify how high-performance computing can be used effectively to solve the transport equation and advance the state-of-the-art in simulating radiative heat transfer. This workshop was held on March 29-30, 1994 in Albuquerque, New Mexico and was conducted by Sandia National Laboratories. The objectives of this workshop were to provide a vehicle to stimulate interest and new research directions within the two communities to exploit the advantages of high-performance computing for solving complex radiative heat transfer problems that are otherwise intractable

Efficient representations of large radiosity matrices

The radiosity equation can be expressed as a linear system, where light interactions between patches of the scene are considered. Its resolution has been one of the main subjects in computer graphics, which has lead to the development of methods focused on different goals. For instance, in inverse lighting problems, it is convenient to solve the radiosity equation thousands of times for static geometries. Also, this calculation needs to consider many (or infinite) light bounces to achieve accurate global illumination results. Several methods have been developed to solve the linear system by finding approximations or other representations of the radiosity matrix, because the full storage of this matrix is memory demanding. Some examples are hierarchical radiosity, progressive refinement approaches, or wavelet radiosity. Even though these methods are memory efficient, they may become slow for many light bounces, due to their iterative nature. Recently, efficient methods have been developed for the direct resolution of the radiosity equation. In this case, the challenge is to reduce the memory requirements of the radiosity matrix, and its inverse. The main objective of this thesis is exploiting the properties of specific problems to reduce the memory requirements of the radiosity problem. Hereby, two types of problems are analyzed. The first problem is to solve radiosity for scenes with a high spatial coherence, such as it happens to some architectural models. The second involves scenes with a high occlusion factor between patches. For the high spatial coherence case, a novel and efficient error-bounded factorization method is presented. It is based on the use of multiple singular value decompositions along with a space filling curve, which allows to exploit spatial coherence. This technique accelerates the factorization of in-core matrices, and allows to work with out-of-core matrices passing only one time over them. In the experimental analysis, the presented method is applied to scenes up to 163K patches. After a precomputation stage, it is used to solve the radiosity equation for fixed geometries and infinite bounces, at interactive times. For the high occlusion problem, city models are used. In this case, the sparsity of the radiosity matrix is exploited. An approach for radiative exchange computation is proposed, where the inverse of the radiosity matrix is approximated. In this calculation, near-zero elements are removed, leading to a highly sparse result. This technique is applied to simulate daylight in urban environments composed by up to 140k patches.La ecuación de radiosidad tiene por objetivo el cálculo de la interacción de la luz con los elementos de la escena. Esta se puede expresar como un sistema lineal, cuya resolución ha derivado en el desarrollo de diversos métodos gráficos para satisfacer propósitos específicos. Por ejemplo, en problemas inversos de iluminación para geometrías estáticas, se debe resolver la ecuación de radiosidad miles de veces. Además, este cálculo debe considerar muchos (infinitos) rebotes de luz, si se quieren obtener resultados precisos de iluminación global. Entre los métodos desarrollados, se destacan aquellos que generan aproximaciones u otras representaciones de la matriz de radiosidad, debido a que su almacenamiento requiere grandes cantidades de memoria. Algunos ejemplos de estas técnicas son la radiosidad jerárquica, el refinamiento progresivo y la radiosidad basada en wavelets. Si bien estos métodos son eficientes en cuanto a memoria, pueden ser lentos cuando se requiere el cálculo de muchos rebotes de luz, debido a su naturaleza iterativa. Recientemente se han desarrollado métodos eficientes para la resolución directa de la ecuación de radiosidad, basados en el pre-cómputo de la inversa de la matriz de radiosidad. En estos casos, el desafío consiste en reducir los requerimientos de memoria y tiempo de ejecución para el cálculo de la matriz y de su inversa. El principal objetivo de la tesis consiste en explotar propiedades específicas de ciertos problemas de iluminación para reducir los requerimientos de memoria de la ecuación de radiosidad. En este contexto, se analizan dos casos diferentes. El primero consiste en hallar la radiosidad para escenas con alta coherencia espacial, tal como ocurre en algunos modelos arquitectónicos. El segundo involucra escenas con un elevado factor de oclusión entre parches. Para el caso de alta coherencia espacial, se presenta un nuevo método de factorización de matrices que es computacionalmente eficiente y que genera aproximaciones cuyo error es configurable. Está basado en el uso de múltiples descomposiciones en valores singulares (SVD) junto a una curva de recubrimiento espacial, lo que permite explotar la coherencia espacial. Esta técnica acelera la factorización de matrices que entran en memoria, y permite trabajar con matrices que no entran en memoria, recorriéndolas una única vez. En el análisis experimental, el método presentado es aplicado a escenas de hasta 163 mil parches. Luego de una etapa de precómputo, se logra resolver la ecuación de radiosidad en tiempos interactivos, para geométricas estáticas e infinitos rebotes. Para el problema de alta oclusión, se utilizan modelos de ciudades. En este caso, se aprovecha la baja densidad de la matriz de radiosidad, y se propone una técnica para el cálculo aproximado de su inversa. En este cálculo, los elementos cercanos a cero son eliminados. La técnica es aplicada a la simulación de la luz natural en ambientes urbanos compuestos por hasta 140 mil parches

Doctor of Philosophy

dissertationRadiation is the dominant mode of heat transfer in high temperature combustion environments. Radiative heat transfer affects the gas and particle phases, including all the associated combustion chemistry. The radiative properties are in turn affected by the turbulent flow field. This bi-directional coupling of radiation turbulence interactions poses a major challenge in creating parallel-capable, high-fidelity combustion simulations. In this work, a new model was developed in which reciprocal monte carlo radiation was coupled with a turbulent, large-eddy simulation combustion model. A technique wherein domain patches are stitched together was implemented to allow for scalable parallelism. The combustion model runs in parallel on a decomposed domain. The radiation model runs in parallel on a recomposed domain. The recomposed domain is stored on each processor after information sharing of the decomposed domain is handled via the message passing interface. Verification and validation testing of the new radiation model were favorable. Strong scaling analyses were performed on the Ember cluster and the Titan cluster for the CPU-radiation model and GPU-radiation model, respectively. The model demonstrated strong scaling to over 1,700 and 16,000 processing cores on Ember and Titan, respectively

Modelling heat transfer and respiration of occupants in indoor climate

Although the terms "Human Thermal Comfort" and "Indoor Air Quality (IAQ)" can be highly subjective, they still dictate the indoor climate design (HVAC design) of a building. In order to evaluate human thermal comfort and IAQ, one of three main tools are used, a) direct questioning the subjects about their thermal and air quality sensation (voting, sampling etc.), b) measuring the human thermal comfort by recording the physical parameters such as relative humidity, air and radiation temperature, air velocities and concentration gradients of pollutants or c) by using numerical simulations either including or excluding detailed thermo-physiological models. The application of the first two approaches can only take place in post commissioning and/or testing phases of the building. Use of numerical techniques can however be employed at any stage of the building design. With the rapid development in computational hard- and software technology, the costs involved in numerical studies has reduced compared to detailed tests. Employing numerical modelling to investigate human thermal comfort and IAQ however demand thorough verification and validation studies. Such studies are used to understand the limitations and application of numerical modelling of human thermal comfort and IAQ in indoor climates. This PhD research is an endeavour to verify, validate and apply, numerical simulation for modelling heat transfer and respiration of occupants in indoor climates. Along with the investigations concerning convective and radiation heat transfer between the occupants and their surroundings, the work focuses on detailed respiration modelling of sedentary human occupants. The objectives of the work have been to: verify the convective and radiation numerical models; validate them for buoyancy-driven flows due to human occupants in indoor climates; and apply these validated models for investigating human thermal comfort and IAQ in a real classroom for which field study data was available. On the basis of the detailed verification, validation and application studies, the findings are summarized as a set of guidelines for simulating human thermal comfort and IAQ in indoor climates. This PhD research involves the use of detailed human body geometries and postures. Modelling radiation and investigating the effect of geometrical posture has shown that the effective radiation area varies significantly with posture. The simulation results have shown that by using an effective radiation area factor of 0.725, estimated previously (Fanger, 1972) for a standing person, can lead to an underestimation of effective radiation area by 13% for the postures considered. Numerical modelling of convective heat transfer and respiration processes for sedentary manikins have shown that the SST turbulence model (Menter, 1994) with appropriate resolution of near wall region can simulate the local air velocity, temperature and heat transfer coefficients to a level of detail required for prediction of thermal comfort and IAQ. The present PhD work has shown that in a convection dominated environment, the detailed seated manikins give rise to an asymmetrical thermal plume as compared to the thermal plumes generated by simplified manikins or point sources. Validated simulation results obtained during the present PhD work have shown that simplified manikins can be used without significant limitations while investigating IAQ of complete indoor spaces. The use of simplified manikins however does not seem appropriate when simulating detailed respiration effects in the immediate vicinity of seated humans because of the underestimation in the amount of re-inhaled CO2 and pollutants from the surroundings. Furthermore, the results have shown that due to the simplification in geometrical form of the nostrils, the CO2 concentration is much higher near the face region (direct jet along the nostrils) as compared to a detailed geometry (sideways jet). Simulating the complete respiration cycle has shown that a pause between exhalation and inhalation has a significant effect on the amount of re-inhaled CO2. Previous results have shown the amount of re-inhaled CO2 to range between 10 - 19%. The present study has shown that by considering the pause, this amount of re-inhaled CO2 falls down to values lower than 1%. A comparison between the simplified and detailed geometry has shown that a simplified geometry can cause an underestimation in the amount of re-inhaled CO2 by more than 37% as compared to a detailed geometry. The major contribution to knowledge delivered by this PhD work is the provision of a validated seated computational thermal manikin. This PhD work follows a structured verification and validation approach for conducting CFD simulations to predict human thermal comfort and indoor air quality. The work demonstrates the application of the validated model to a classroom case with multiple occupancy and compares the measured results with the simulation results. The comparison of CFD results with measured data advocates the use of CFD and visualizes the importance of modelling thermal manikins in indoor HVAC design rather than designing the HVAC by considering empty spaces as the occupancy has a strong influence on the indoor air flow. This PhD work enables the indoor climate researchers and building designers to employ simplified thermal manikin to correctly predict the mean flow characteristics in indoor surroundings. The present work clearly demonstrates the limitation of the PIV measurement technique, the importance of using detailed CFD manikin geometry when investigating the phenomena of respiration in detail and the effect of thermal plume around the seated manikin. This computational thermal manikin used in this work is valid for a seated adult female geometry

Modelling polarized light for computer graphics

The quality of visual realism in computer generated images is largely determined by the accuracy of the reflection model. Advances in global illumination techniques have removed to a large extent, some of the limitations on the physical correctness achievable by reflection models. While models currently used by researchers are physically based, most approaches have ignored the polarization of light. The few previous efforts addressing the polarization of light were hampered by inherently unphysical light transport algorithms. This paper, besides taking polarization of light into account in the reflection computation, also provides a basis for modelling polarization as an inherent attribute of light, using the Stokes parameters. A reflection model is developed within this framework and the implementation within a global illumination algorithm called Photon is presented

The Future's Brighter at the Edges: Modelling and Designs for Efficient Luminescent Solar Concentrators

Luminescent Solar Concentrators (LSCs) are fluorescent particle doped planar devices which collect light incident on their large collection faces, and concentrate it to the edges. There, the light can be coupled into photovoltaic devices and converted into electricity. In this thesis I present a series of studies, combining fabrication, characterisation, computational modelling and analytical solutions in order to progress the understanding of, and to increase, the efficiency of LSC devices. The work is structured around a back-bone of an in-house built experimentally verified Monte- Carlo ray tracing platform, which accurately predicts the efficiency and losses of LSC devices. I have extended the modelling platform by integrating a novel time-domain aspect to the simulations, predicting the impulse response of LSC devices. I present experimental verification of these predictions, whilst also demonstrating their application in the context of bandwidth prediction for the use of LSCs as detectors for data transmission in Visual Light Communications. The technique is not limited to this application, as it can provide a deeper understanding of the advancement of LSCs for all applications. Importantly, the modifications are simple to integrate into existing simulation platforms. I have also expanded the model to investigate the effects of curvature on LSC efficiency. Curvature and flexibility in LSCs can expand the range of available applications and fabrication techniques. Here, I demonstrate the potentially detrimental effects induced by curvature as well as the mechanisms behind them. Furthermore, I have provided a potential solution to reduce curvature induced losses by demonstrating a novel, fully flexible, LSC – distributed Bragg reflector combination. I have shown that the all-silicone based design, fabricated using scalable techniques, not only boosts LSC efficiency but also severely decouples the efficiency of such devices from the effects of curvature, paving a route to highly efficient flexible LSCs