La présence de fines particules de matières oxydables est rencontrée dans de nombreuses situations industrielles. Le risque d'explosion de poussières présente une menace constante pour les industries de transformation qui fabriquent, utilisent ou manipulent des poudres ou despoussières de matières combustibles. Dans le secteur nucléaire, les scénarios envisagés traitent,en particulier, le risque d'explosion de poussières de graphite liées aux opérations dedémantèlement des réacteurs Uranium Naturel Graphite Gaz. La problématique considérée, dans le cadre de ce travail de thèse, est celle de la combustion d'un mélange dilué gaz-particules.L'objectif de cette thèse est de développer un modèle Euler-Lagrange macroscopique permettantde prédire la vitesse laminaire de flamme qui est une des données essentielles pour les modèlesde vitesse de flamme turbulente utilisés dans l'évaluation des risques d'explosion de poussières.Dans un premier temps, les équations macroscopiques de transferts massique et thermique sont dérivées à partir de la méthode de prise de moyenne volumique. L'intérêt de l'approche utilisée ici est de proposer des problèmes de fermeture permettant d'estimer les coefficients de transfertseffectifs, tels que les coefficients d'échanges thermiques et le coefficient effectif de la réactionhétérogène. Dans un deuxième temps, des simulations Euler-Lagrange sont utilisées pourdéterminer la vitesse de flamme laminaire diphasique plane en fonction des caractéristiques du mélange gazeux et des poussières de graphite. Le modèle proposé dans ce travail est comparé au modèle Euler-Lagrange classique basé sur la résolution du problème de couche limite pourune particule isolée en milieu infini. Cette étude montre que les effets du taux de dilution et deséchanges indirects entre les particules ne sont pas systématiquement négligeables dans leséchanges macroscopiques entre les deux phases. D'autre part, la présente étude laisse entrevoir la potentialité de l'approche proposée pour les simulations détaillées de l'écoulement diphasique / The presence of fine particles of oxidizable materials is encountered in many industrial situations.The risk of dust explosion presents a constant threat in transformation industries that manufacture,use or manipulate powders or combustible materials dusts. In nuclear safety analysis, one of themain scenarios is the risk of graphite dust explosion that may occur during decommissioningoperations of Uranium Natural Graphite Gas reactors. The issue considered in this thesis isrelated to combustion of a dilute gas-particle mixture. This work aims at developing a macroscopicEuler-Lagrange model for predicting laminar flame velocity, which is one of the essential data forturbulent flame velocity models used to evaluate the risk of dust explosion. First, the macroscopicheat and mass transfer equations are derived using the volume averaging method. The majorinterest of the proposed approach is to provide closure problems that allow to estimate theeffective transport coefficients, such as heat exchange coefficients and the effective coefficient ofthe heterogeneous reaction. Second, Euler-Lagrange simulations are used to determine the planetwo-phase laminar flame velocity as a function of gas mixture and graphite dust characteristics.The proposed model is compared to the classical Euler-Lagrange model based on the resolutionof the boundary layer problem in the vicinity of an isolated particle in infinite medium. Results showthat the dilution rate and the indirect particle-particle exchanges are not systematically negligible inthe macroscopic exchanges between the two-phases. On the other hand, this study suggests thepotentiality of the proposed approach for detailed simulations of two-phase flow

Dropwise condensation, which collects the condensate liquid in the form of droplets, has attracted a growing interest due to much higher heat transfer coefficient. One important and challenging issue in dropwise condensation is the presence of non-condensable gas (NCG) which drastically reduces its heat transfer performance. Concerning the mechanism understanding, this thesis is aiming to investigate dropwise condensation in case of NCG by combing different methods. Firstly, convective dropwise condensation out of moist air is experimentally investigated under controllable conditions. In modeling, some crucial aspects are reasonably captured: the coupled heat and mass transfer during droplet growth by a multi-scale droplet growth model; the inter-droplet interaction defined by a distributed point sink method; the enhancement of the convective mass transfer using the droplet Sherwood number. Furthermore, a multi-component multi-phase thermal pseudopotential-based LB model is developed to advance the directly numerical simulation of dropwise condensation.

Zur Kühlung von Prozessströmen kommen aufgrund hoher Leistungsdichten häufig Verdunstungskühltürme zum Einsatz. Um die Übertragungsfläche für Wärme und Stoff zu vergrößern, werden in diesen Kühltürmen Struktureinbauten integriert. Die Weiterentwicklung von Kühlturmeinbauten und die Untersuchung der den Kühlprozess beeinflussenden Faktoren erfolgt empirisch, was eine Vielzahl von Versuchen notwendig macht. Eine numerische Simulation des Kühlprozesses kann diese Messungen unterstützen und so helfen eine Vielzahl an Versuchen einzusparen. Des Weiteren können bei versuchsbegleitender Simulation mit einem geeigneten Modell weitere Untersuchungen durchgeführt und Erkenntnisse gewonnen werden, die bei Messungen am Versuchskühlturm verborgen bleiben. In dieser Arbeit werden zwei Ansätze der numerischen Simulation eines Verdunstungskühlturms betrachtet. Es werden eine CFD-Simulation und ein vereinfachtes Modellkonzept hinsichtlich der Anwendbarkeit auf diese Problemstellung untersucht. Schwerpunkt der vorliegenden Arbeit ist die methodische Entwicklung eines solchen vereinfachten mathematischen Modells. Dieses beruht auf der physikalisch deterministischen Beschreibung der im Kühlturm ablaufenden Prozesse der Wärme- und Stoffübertragung unter Berücksichtigung des Stoffverhaltens. Aufgrund der Nichtlinearität des Stoffverhaltens und der erforderlichen Inkrementierung des Berechnungsgebiets ist ein methodisches Vorgehen erforderlich, um die Erstellung der Modellgleichungen und deren Lösung überhaupt realisieren zu können. Hierfür wird auf allgemeine Methoden der gleichungsorientierten Simulation technischer Systeme zurückgegriffen. Das entwickelte Modellkonzept wird für die Modellierung und Simulation eines Versuchskühlturms angewandt. Mit den so ermittelten Messdaten wird das Modell kalibriert und validiert. Es zeigt sich, dass mit dem erstellten Modell quantitativ und qualitativ valide Ergebnisse erzielt werden können. / Due to the high power density, the cooling of process streams is often done bei evaporative cooling towers. To enlarge the exchange area for the heat and mass transfer, these cooling towers contain integrated structural fills. The future development of cooling tower fills and the research regarding the cooling process and its influencing parameters will be carried out empirically, resulting in a large number of required experiments. A numeric simulation of the cooling process can support theses measurements and reduce the vast number of needed experiments. Furthermore, with the use of test-related simulations and adapted models, it will be possible to gain knowledge and do research in areas which are omitted during regular measurements on cooling towers. In this study it is looked to two different approaches of numeric simulation of a evaporative cooling tower. There will be an examination of a CFD-Simulation and a simplified model concept regarding their respective applicability for this problem. This work is focussed on the systematic developement of such simplified mathematical models, based on the physical deterministic description of the occurring processes of heat and mass transfer in cooling towers considering the stock behaviour. Due to the non-linearity of the stock behaviour and the required incrementation of the calculation area, a systematic approach is needed to model equations and their respective solutions. For this purpose it is necessary to access general techniques of equation-based simulations of technological systems. The developed model concept will be applied for the modelation and simulation of an experimental cooling tower. The model will be calibrated and validated with data from this experimental tower. It shows, that the results from this model are qualitatively and quantitatevily valid.

Beim Stoffübergang einer grenzflächenaktiven Substanz in einem flüssigen Zweiphasensystem kann solutale Marangoni-Instabilität einsetzen. Die weitere nichtlineare Entwicklung der Marangoni-Instabilität geht mit einer enormen Vielfalt von Strömungsmustern einher. In der Literatur wird dieser Aspekt häufig unter dem unscharfen Ausdruck „Grenzflächenturbulenz“ zusammengefasst. Diese Arbeit stellt heraus, dass drei grundlegende Strukturformen existieren: Rollzellen, Relaxationsoszillationen und Relaxationsoszillationswellen. Ein großer Teil der Komplexität der Strömungsmuster ist dadurch begründet, dass die Grundstrukturen unterschiedliche Hierarchieebenen aufweisen. Es werden die zugrunde liegenden Bedingungen für das Auftreten der jeweiligen Strukturtypen, ihre transiente Natur und die Bildung der hierarchischen Strömungsmuster untersucht. Des Weiteren betrachtet diese Arbeit die Wechselwirkungen mit Dichteeffekten, die sowohl die Charakteristik der Strukturen als auch ihre zeitliche Entwicklung beeinflussen.

Es werden die Gleichgewichtskonstanten der Zinkextraktion ermittelt. Das Wilson-Modell und das erweiterte Debye-Hückel-Gesetz werden zur Beschreibung der Aktivitäten verwendet. Die tensiometrischen Untersuchungen erfolgen am hängenden Tropfen. Die Modellauswahl zur Beschreibung der Gleichgewichtsgrenzflächenspannung erfolgt im submizellaren Konzentrationsbereich. Die pseudo-nichtionische Modellierung auf der Basis der Langmuir-Isothermen der Mehrkomponentenadsorption bei Verwendung der Stern-Isothermen für die Gegenionenanreicherung liefert die beste Datenanpassung. Durch ein einfaches Modell zur Mizell- und Aggregatbildung gelingt die Modellerweiterung. Die gemessenen dynamischen Grenzflächenspannungskurven werden sorptionskinetisch und durch diffusive Approximationen angepasst. Es zeigt sich, dass der Stofftransport diffusionsdirigiert ist. Am oszillierenden Tropfen werden die Ergebnisse bestätigt. Für das Kationenaustauscheranion wird die Gültigkeit des Maxwell-Modells zur Beschreibung der Grenzflächendilatationsrheologie nachgewiesen.

<p>With the extreme increases in power density of electronic devices, the contact thermal resistance imposed at interfaces between mating solids becomes a major challenge in thermal management. This contact thermal resistance is mainly caused by micro-scale surface asperities (roughness) and wavy profile of surface (nonflatness) which severely reduce the contact area available for heat conduction. High contact pressures (1~100 MPa) can be used to deform the surface asperities to increase contact area. Besides, a variety of conventional thermal interface materials (TIM), such as greases and pastes, are used to improve the contact thermal conductance by filling the remaining air gaps. However, there are still some applications where such TIMs are disallowed for reworkability concerns. For example, heat must be transferred across dry interfaces to a heat sink in pluggable opto-electronic transceivers which needs to repeatedly slide into / out of contact with the heat sink. Dry contact and low contact pressures are required for this sliding application.</p> <p>This dissertation presents a metallized micro-spring array as a surface coating to enhance dry contact thermal conductance under ultra-low interfacial contact pressure. The shape of the micro-springs is designed to be mechanically compliant to achieve conformal contact between nonflat surfaces. The polymer scaffolds of the micro-structured TIMs are fabricated by using a custom projection micro-stereolithography (μSL) system. By applying the projection scheme, this method is more cost-effective and high-throughput than other 3D micro-fabrication methods using a scanning scheme. The thermal conductance of polymer micro-springs is further enhanced by metallization using plating and surface polishing on their top surfaces. The measured mechanical compliance of TIMs indicates that they can deform ~10s μm under ~10s kPa contact pressures over their footprint area, which is large enough to accommodate most of surface nonflatness of electronic packages. The measured thermal resistances of the TIM at different fabrication stages confirms the enhanced thermal conductance by applying metallization and surface polishing. Thermal resistances of the TIMs are compared to direct metal-to-metal contact thermal resistance for flat and nonflat mating surfaces, which confirms that the TIM outperforms direct contact. A thin layer of soft polymer is coated on the top surfaces of the TIMs to accommodate surface roughness that has a smaller spatial period than the micro-springs. For rough surfaces, the polymer-coated TIM has reduced thermal resistance which is comparable to a benchmark case where the top surfaces of the TIM are glued to the mating surface. A polymer base is designed under the micro-spring array which can provide the advantages for handling as a standalone material or integration convenience, at the toll of an increased insertion resistance. Through-holes are designed in the base layer and coated with thermally conductive metal after metallization to enhance thermal conductance of the base layer; a thin layer of epoxy is applied between the base layer and the working surface to reduce contact thermal resistance exposed on the base layer. Cycling tests are conducted on the TIMs; the results show good early-stage reliability of the TIM under normal pressure, sliding contact, and temperature cycles. The TIM is thermally demonstrated on a pluggable application, namely, a CFP4 module, which shows enhanced thermal conductance by applying the TIM. </p> To further enhance the potential mechanical compliance of microstructured surfaces, a stable double curved beam structure with near-zero stiffness composed of intrinsic negative and positive stiffness elastic elements is designed and fabricated by introducing residual stresses. Stiffness measurements shows that the positive-stiffness single curved beam, which is the same as the top beam in the double curved beam, is stiffer than the double curved beam, which confirms the negative stiffness of the bottom beam in the double curved beam. Layered near zero-stiffness materials made of these structures are built to demonstrate the scalability of the zero-stiffness zone.

Several types of reactors are encountered in industry where reagents in a gas and a liquid phase need to be catalysed by a solid catalyst. Common reactors that are used to this end, are trickle-bed reactors, where gas and liquid flow cocurrently down a packed bed of catalyst. Apart from the catalytic process itself, several mass transfer steps can influence the rate and/or selectivity of a solid catalysed gas-liquid reaction. In trickle-bed reactors, flow morphology can have a major effect on these mass transfer steps. This study investigates the interaction between liquid flow morphology and mass transfer in trickle-bed reactors from three different angles. The primary focus is on liquid-solid mass transfer and internal diffusion as affected by the contacting between the liquid and the catalyst. First, the contacting between the liquid and the solid in trickleflow, or wetting efficiency, is characterised using colorimetry. Though this investigation is limited to the flow of nitrogen and water over a packed bed at ambient conditions, it provides useful information regarding liquid flow multiplicity behaviour and its influence on the distribution of fractional wetting on a particle scale. The colorimetric study also provides descriptions of the geometry of the liquid-solid contacting on partially wetted particles. These are used in a second investigation, for the numerical simulation of reaction and diffusion in partially wetted catalysts. This second investigation uses numerical simulations to evaluate and develop simple theoretical descriptions of liquid-solid contacting effects on catalyst particle efficiency. Special attention is given to the case where external and intraparticle mass transfer rates of both a volatile and non-volatile reagent affect the overall rate of reaction. Also, since these are not often considered in theoretical studies, some suggestions are made for the evaluation of the particle efficiency of eggshell catalyst. Finally, liquid-solid contacting is investigated in a high-pressure pilot reactor. Wetting efficiency is measured with a useful technique that does not rely on descriptions of particle kinetics or liquid-solid mass transfer rates. Liquid-solid mass transfer coefficients are also measured and results agree well with the colorimetric investigation, suggesting the existence of different types of flow within in the hydrodynamic multiplicity envelope of trickle-flow. Since it consists of different investigations of liquid-solid contacting from different angles, the study highlights several aspects of liquid-solid contacting and how it can be expected to influence trickle-bed reactor performance. / Thesis (PhD)--University of Pretoria, 2009. / Chemical Engineering / unrestricted

<div>Thermal management in next generation power electronic devices, radar applications and semiconductor packaging architectures is becoming increasingly challenging due to the need to reject localized high heat fluxes as well as large total powers. Air cooling has been considered as a simple and reliable method for thermal management compared to architectures that incorporate liquid cooling. However, air-cooled heat sinks typically require effective heat spreading to provide the requisite level of area enhancement to dissipate high heat fluxes. Compared to solid metallic heat spreaders, advanced heat sinks that incorporate two-phase heat transfer devices such as vapor chambers can significantly enhance the power dissipation capabilities in such configurations. Vapor chambers are devices that utilize evaporation/boiling processes within a sealed cavity to achieve efficient heat spreading. In high-heat-flux applications, boiling can occur within the internal wick structure of the vapor chamber at the location of the heat input (i.e., the evaporator). The maximum dryout heat flux and thermal resistance of the device is dictated by the resulting two-phase flow and heat transfer in the porous evaporator due to boiling. While various works in the literature have introduced new evaporator wick designs to improve the dryout heat flux during boiling, the enhancement is limited to small, millimeter scale hotspots or at a very high thermal resistance. In additixon, the effective design of such evaporator systems requires mechanistic models that can accurately predict the dryout limit and thermal performance. </div><div> This thesis first explores the usage of a novel ‘two-layer’ evaporator wick for passive high heat flux dissipation over large heater areas at a low thermal resistance. Moreover, a new mechanistic (first principles based) model framework is introduced for dryout limit and thermal performance prediction during boiling in capillary fed evaporators, by considering the resulting simultaneous flow of two phases (liquid and vapor) within the microscale porous media.</div><div> The novel two-layer wick concept uses a thick ‘cap’ layer of porous material to feed liquid to a thin ‘base’ layer through an array of vertical liquid-feeding ‘posts’. Vapor ‘vents’ in the cap layer allow for vapor formed during the boiling process (which is constrained to the base layer) to escape out of the wick. This two-layer structure decouples the functions of liquid resupply and capillary-fed boiling heat transfer, making the design realize high heat flux dissipation greater than 500 W/cm2 over large heat input areas of ~1 cm2. A reduced-order model is first developed to demonstrate the performance of a vapor chamber incorporating such a two-layer evaporator wick design. The model comprises simplified hydraulic and thermal resistance networks for predicting the capillary-limited maximum heat flux and the overall thermal resistance, respectively. The reduced-order model is validated against a higher fidelity numerical model and then used to analyze the performance of the vapor chamber with varying two-layer wick geometric feature sizes. The fabrication of the proposed two-layer wick is then presented. The thermal performance of the fabricated wicks is characterized using a boiling test facility that utilizes high speed visualization to identify the characteristic regimes of boiling operation in the wicks. The performance is also benchmarked to conventional single-layer wicks. </div><div> It is observed that single-layer wicks exhibit an unfavorable boiling regime where the center of the heater area dries out locally, leading to a high value of thermal resistance. The two-layer wicks avoid local dryout due to the distributed feeding provided by the posts and enhance the dryout heat flux significantly compared to single-layer wicks. A two-layer design that consists of a 10 × 10 array of liquid feeding posts provided a 400% improvement in the dryout heat flux. Following a parametric analysis of the effect of particle size, two-layer wicks composed of 180 – 212 µm particles and a 15 × 15 array of liquid feeding posts yielded a maximum heat flux dissipation of 485 W/cm2 over a 1 cm2 heat input area while also maintaining a low thermal resistance of only ~0.052 K/W. The effect of vapor venting and liquid-feeding areas is also experimentally studied. By understanding these effects, a parametrically optimized design is fabricated and shown to demonstrate an extremely high dryout limit of 512 W/cm2. We identify that the unique area-scalability of the two-layer wick design allows it to achieve an unprecedented combination of high total power and low-thermal-resistance heat dissipation over larger areas than was previously possible in the literature.</div><div> The results from the characterization of two-layer wicks revealed that the overall performance of the design was limited by the boiling process in the thin base wick layer. A fundamental model-based understanding of the resulting two-phase flow and heat transfer process in such thin capillary-fed porous media was still lacking. This lack of a mechanistic model precluded the accurate prediction of dryout heat flux and thermal performance of the two-layer wick. Moreover, such an understanding is needed for the optimal design of advanced hybrid evaporator wicks that leverage capillary-fed boiling. Despite the existence of various experimental works, there are currently no mechanistic approaches that model this behavior. To fill this unmet need, this thesis presents a new semi-empirical model for prediction of dryout and thermal resistance of capillary-fed evaporator systems. Thermal conduction across the solid and volumetric evaporation within the pores are solved to obtain the temperature distribution in the porous structure. Capillary-driven lateral liquid flow from the outer periphery of the evaporator to its center, with vapor flow across the thickness, is considered to obtain the local liquid and vapor pressures. Experiments are conducted on sintered copper particle evaporators of different particle sizes and heater areas to collect data for model calibration. To demonstrate the wider applicability of the model for other types of porous evaporators, the model is further calibrated against a variety of dryout limit and thermal resistance data collected from the literature. The model is shown to predict the experimentally observed trends in the dryout limit with mean particle/pore size, heater size, and evaporator thicknesses. This physics–based modeling approach is then implemented into a vapor chamber model to predict the thermal performance limits of air-cooled heat sinks with embedded vapor chambers. The governing energy and momentum equations of a low-cost analytical vapor chamber modeling approach is coupled with the evaporator model to capture the effect of boiling in the evaporator wick. An example case study illustrating the usage of the model is demonstrated and compared to a purely evaporation-based modeling approach, for quantifying the differences in dryout limit prediction, signifying the need to account for boiling in the evaporator wick. </div><div> The understanding gained from this thesis can be utilized for the prediction of dryout and thermal performance during boiling in capillary limited evaporator systems. The work also suggests the usage of a universal relative permeability correlation for the two-phase flow configuration studied herein for capillary-fed boiling, based on a wide calibration to experimental data. The modeling framework can also be readily leveraged to find novel and unexplored designs of advanced evaporator wicks. From an application standpoint, the new vapor chamber model developed here can be used for the improved estimation of performance limits specifically when high heat fluxes are encountered by the device. This will enable better and informed design of air-cooled heat sink architectures with embedded vapor chambers for high performance applications. </div><div><br></div>

Energy recovery ventilation (ERV) systems are designed to decrease the energy consumed by building HVAC systems. ERV’s scavenge sensible and latent energy from the exhaust air leaving a building or space and recycle this energy content to pre-condition the entering outdoor air. A few studies found in the open literature are dedicated to developing detailed numerical models to predict or simulate the performance of energy recovery wheels and desiccant wheels. However, the models are often computationally intensive, requiring a lot of time to perform parametric studies. For example, if the physical characteristics of a study target change (e.g., wheel diameter or depth) or if the system runs at different operating conditions (e.g., wheel rotation speed or airflow rate), the model parameters need to be recalculated. Hence, developing a mapping method with better computational efficiency, which will enable the opportunity to conduct extensive parametric or optimal design studies for different wheels is the goal of this research. In this work, finite difference method (FDM) numerical models of energy recovery wheels and desiccant wheels are established and validated with laboratory test results. The FDM models are then used to provide data for the development of performance mapping methods for an energy wheel or a desiccant wheel. After validating these new mapping approaches, they are employed using independent data sets from different laboratories and other sources available in the literature to identify their universality. One significant characteristic of the proposed mapping methods that makes the contribution unique is that once the models are trained, they can be used to predict performance for other wheels with different physical geometries or different operating conditions if the desiccant material is identical. The methods provide a computationally efficient performance prediction tool; therefore, they are ideal to integrate with transient building energy simulation software to conduct performance evaluations or optimizations of energy recovery/ desiccant wheel integrated HVAC systems.

Bench and pilot scale determinations of the volumetric oxygen transfer coefficient, K_La, were performed on an improved A²/O biological nutrient removal (BNR) pilot plant. Effluent from a full scale primary clarifier, used as pilot plant influent, was found to have an alpha (ratio of process to clean water K_La) of 0.71 as determined in a 21 liter bench scale reactor and an alpha of 0.332 as determined in a 0.45 m³ aeration basin of the 2.4 m³ pilot plant. Alpha of a 1:1 mixture of primary clarifier effluent with pilot plant return activated sludge was determined to be 0.94 at bench scale and 0.71 at pilot scale. An assay of alphas through the initial non aerated treatment zones of the pilot plant using the bench scale reactor indicated that alphas peaked in the effluent of the first anaerobic zone (alpha equal to 1.01) and were lower in the second anaerobic zone and first anoxic zone. An assay of alphas in the three pilot plant series sideline aeration basins indicated that alpha was maximum in the first aeration basin (alpha equal to 0.905) and were lower in the second and third aeration basins (0.716 and 0.661 respectively). A consistent increase in average surface tension was noted from the first to second to third aeration basins, however the differences were not statistically significant. A comparison of pilot plant alphas determined in the first aeration basin following anaerobic nominal hydraulic retention times of 0.0, 0.21, 0.43, and 0.64 hours yielded alpha values of 0.71, 0.94, 0.64, and 0.74 respectively. Like the assay using the bench scale reactor, the alpha values at pilot scale peaked following treatment in only one anaerobic zone (nominal HRT of 0.21 hours). The study concludes that short exposures in an initial anaerobic reactor as required for biological phosphorus removal may benefit oxygen transfer efficiency through increased alphas, however the benefits of long periods of anaerobic reaction time (over 0.43 hours) are uncertain. / Master of Science