Dissertationen zum Thema „Evaporative flux“

Thesis: S.M., Massachusetts Institute of Technology, Department of Mechanical Engineering, 2017.
Cataloged from PDF version of thesis.
Includes bibliographical references (pages 60-63).
We investigated the experimental performance of a nanoporous membrane for ultra-high heat flux dissipation from high performance integrated circuits. The biporous evaporation device utilizes thermally-connected, mechanically-supported, high capillarity membranes that maximize thin film evaporation and high permeability liquid supply channels that allow for lower viscous pressure losses. The 600 nm thick membrane was fabricated on a silicon on insulator (SOI) wafer, fusion-bonded to a separate wafer with larger liquid channels. Spreading effects and overall device performance arising from non-uniform heating and evaporation of methanol was captured experimentally. Heat fluxes up to 412 W/cm2, over an area of O.4x 5 mm, and with a temperature rise of 24.1 K from the heated substrate to ambient vapor, were obtained. These results are in good agreement with a high-fidelity, coupled fluid convection and solid conduction compact model, which was necessitated by computational feasibility, which incorporates non-equilibrium and sub-continuum effects at the liquid-vapor interface. This work provides a proof-of-concept demonstration of our biporous evaporation device. Simulations from the validated model, at optimized operating conditions and with improved working fluids, predict heat dissipation in excess of 1 kW/cm2 with a device temperature rise below 30 K, for this scalable cooling approach.
by Jay D. Sircar.
S.M.

Le séchage assisté par champ électrique est une technique prometteuse pour la manipulation des produits thermosensibles, mais son étude à l'échelle de la goutte reste limitée en raison du manque de dispositifs expérimentaux adaptés. Cette étude présente un dispositif permettant d'analyser la déformation et le séchage de gouttes pendantes électriquement chargées de CuSO4.5H4O sous un champ électrique (0-800 kV/m) et à différentes concentrations en solides (0-25 % wt/wt). Les résultats montrent une réduction de la tension interfaciale pouvant atteindre 42 % avec l'intensité du champ. L'expansion de surface et la déformation sont également influencées. En séchage diffusif, le champ réduit initialement le taux d'évaporation et le flux évaporatoire des gouttes déformées avant qu'il n'augmente et se stabilise, tandis qu'en séchage convectif, l'expansion de surface réduit temporairement le flux avant stabilisation. Ces résultats apportent de nouvelles perspectives sur les effets électrohydrodynamiques lors du séchage à l'échelle de la goutte
Electrically assisted drying is a promising technique for thermosensitive products, however its study at the droplet scale remains limited due to the lack of suitable experimental setups. This study introduces a novel device to analyze the deformation and drying of pendant electrically charged droplets of CuSO4.5H4O under an electric field (0-800 kV/m) and varying solid concentrations (0 - 25 wt/wt %). Results show up to a 42 % reduction in interfacial tension with increasing field intensity. Surface expansion and deformation are influenced by both field intensity and solid content. In terms of drying, results show that under diffusive drying, the field initially reduces the evaporative flux of highly deformed droplets before it increases and stabilizes, while in convective drying, surface expansion temporarily reduces the flux before stabilization. These findings provide new insights into electrohydrodynamic effects on droplet scale drying dynamics

A novel cooling mechanism based on evaporation of thin liquid films is presented for thermal management of confined heat sources, such as microprocessor hotspots. The underlying idea involves utilization of thin nanoporous membranes for maintaining microscopically thin liquid films by capillary action, while providing a pathway for the vapor generated due to evaporation at the liquid-vapor interface. The vapor generated by evaporation is continuously removed by using a dry sweeping gas keeping the membrane outlet dry. This thesis presents a detailed theoretical, computational and experimental investigation of the heat and mass transfer mechanisms that result in dissipating heat. Performance analysis of this cooling mechanism demonstrates heat fluxes over 600W/cm2 for sufficiently thin membrane and film thicknesses (~1-5µm) and by using air jet impingement for advection of vapor from the membrane surface. Based on the results from this performance analysis, a monolithic micro-fluidic device is designed and fabricated incorporating micro and nanoscale features. This MEMS/NEMS device serves multiple functionalities of hotspot simulation, temperature sensing, and evaporative cooling. Subsequent experimental investigations using this microfluidic device demonstrate heat fluxes in excess of 600W/cm2 at 90 C using water as the evaporating coolant. In order to further enhance the device performance, a comprehensive theoretical and computational analysis of heat and mass transfer at micro and nanoscales is carried out. Since the coolant is confined using a nanoporous membrane, a detailed study of evaporation inside a nanoscale cylindrical pore is performed. The continuum analysis of water confined within a cylindrical nanopore determines the effect of electrostatic interaction and Van der Waals forces in addition to capillarity on the interfacial transport characteristics during evaporation. The detailed analysis demonstrates that the effective thermal resistance offered by the interface is negligible in comparison to the thermal resistance due to the thin film and vapor advection. In order to determine the factors limiting the performance of the MEMS device on a micro-scale, a device-level detailed computational analysis of heat and mass transfer is carried out, which is supported by experimental investigation. Identifying the contribution of various simultaneously occurring cooling mechanisms at different operating conditions, this analysis proposes utilization of hydrophilic membranes for maintaining very thin liquid films and further enhancement in vapor advection at the membrane outlet to achieve higher heat fluxes.

The boreal forests represent 8 % of all forested areas on the earth and have a significant role in the control of greenhouse gases and an impact on global climate change. The main objective of this thesis is to increase the understanding of how evaporation and heat-flux processes in the boreal forest zone are affecting the regional and global climate.

A meteorological mesoscale model with an advanced land-surface parameterization has been utilized to study aggregation of fluxes of water vapour and heat. The model has been compared against four other methods for flux estimation in a southern boreal landscape. The results show that the mesoscale model is successfully reproducing 24-hour averages of fractionally weighted mast measurements of sensible and latent heat flux.

The model was also evaluated against in-situ observations of surface fluxes and other meteorological variables. The results reveal that a correct initialization of soil moisture is crucial to simulate a realistic partitioning of the sensible and latent heat fluxes. Significant differences in surface fluxes and friction velocities between two apparently similar forest sites indicate the need for careful assessment of areal representativity when comparing mesoscale model results with in-situ observations.

A parameterization for the absorption of solar radiation of high-latitude sparse forests was implemented and tested in the model that significantly improved the simulation of high wintertime midday sensible heat fluxes. A scheme for heat storage in vegetation was also implemented which improved the results, but the scheme needs further evaluation for high latitude forests.

Two commonly used strategies for the description of land-surface heterogeneity, the effective parameter approach and the mosaic approach, were tested in the mesoscale model against airborne observations of sensible and latent heat fluxes. The results show that the mosaic approach produces better results especially when small lakes are present in model grid-squares.

Norra halvklotets barrskogsbälte representerar 8 % av all skogsbeklädd mark på jorden och har stor betydelse för kontrollen av växthusgaser och påverkan på globala klimatförändringar. Syftet med denna avhandling är att öka förståelsen av hur avdunstning och värmeflöden i den boreala skogszonen påverkar klimatet regionalt och globalt.

En meteorologisk mesoskalemodell med en avancerad landyteparameterisering har använts för att studera aggregering av avdunstning och värmeflöden. Modellen jämfördes med fyra andra metoder för uppskattning av värmeflöden i den boreala skogszonens södra delar. Resultaten visade att mesoskalemodellen reproducerar 24-timmarsmedelvärden av sensibelt och latent värmeflöde från areellt viktade mastmätningar med bra resultat.

Modellen utvärderades även mot markbaserade mätningar av sensibelt och latent värme och andra meteorologiska variabler. Resultaten visar att en korrekt initialisering av markvatteninnehållet är avgörande för att simulera en realistisk uppdelning av de sensibla och latenta värmeflödena. Markanta skillnader i markyteflöden och friktionshastigheter mellan två liknande skogsmätstationer påvisar nödvändigheten av en noggrann bedömning av den areella representativiteten när man jämför resultat från mesoskalemodellen med markbaserade mätningar.

En parameterisering för absorption av solstrålning i glesa skogsbestånd på höga breddgrader infördes och testades i modellen vilket markant förbättrade simuleringen av de höga sensibla värmeflöden som observerats vid middagstid på vintern. Ett uttryck för att beskriva värmelagring i vegetationen infördes också vilket förbättrade resultaten, men uttrycket behöver vidare utvärdering för skogsbestånd på höga breddgrader.

Två ofta använda strategier för att beskriva markytans heterogenitet, effektiva parametermetoden och mosaikmetoden, testades i mesoskalemodellen mot flygburna observationer av sensibla och latenta värmeflöden. Resultaten visar att mosaikmetoden ger bättre resultat särskilt när mindre sjöar förekommer i modellrutorna.

Cette thèse étudie théoriquement les processus de transport au voisinage de la ligne triple de contact liquide-gaz-solide et leur impact sur l'évaporation macroscopique. Dans la première partie de la thèse, l'hydrodynamique au voisinage de la ligne de contact est étudiée sous les conditions de mouillage partiel. L'évaporation induite par le chauffage du substrat dans l'atmosphère de vapeur du même fluide est considérée. La relaxation de la singularité hydrodynamique de la ligne triple est considérée. La principale conclusion de la thèse est que l'effet Kelvin (dépendance de la température de saturation de la pression) est suffisant en soi, pour faire disparaitre la singularité des variables hydrodynamiques. La microrégion (le voisinage de la ligne de contact) est résolue numériquement et analytiquement pour de faibles pentes de l'interface liquide-vapeur. Les échelles de longueur caractéristiques du problème sont identifiées et la nature multi-échelle du phénomène est prise en compte. Les études paramétriques effectuées révèlent le rôle de la résistance thermique de l'interface vapeur-liquide, de la longueur de glissement, du terme thermocapillaire, du recul de vapeur et ainsi que des forces de surface. Une extension de l'approximation de lubrification pour de pentes élevées de l'interface gaz-liquide à l'évaporation est discutée. Dans la seconde partie de la thèse, le modèle précédemment établi pour la microrégion est couplé à des simulations numériques de la croissance d'une bulle de vapeur. Le départ de la bulle de vapeur de la paroi chauffante pendant l'ébullition a également été étudiée. Il a été proposé dans la thèse, que sous des charges thermiques élevées, l'augmentation de l'angle de contact apparent provoque l'étalement de la bulle de vapeur sur la paroi chauffante. Ce phénomène peut conduire, au séchage de la paroi observé pendant la crise d'ébullition.

Cette thèse étudie théoriquement les processus de transport au voisinage de la ligne triple de contact liquide-gaz-solide et leur impact sur l'évaporation macroscopique. Dans la première partie de la thèse, l'hydrodynamique au voisinage de la ligne de contact est étudiée sous les conditions de mouillage partiel. L'évaporation induite par le chauffage du substrat dans l'atmosphère de vapeur du même fluide est considérée. La relaxation de la singularité hydrodynamique de la ligne triple est considérée. La principale conclusion de la thèse est que l'effet Kelvin (dépendance de la température de saturation de la pression) est suffisant en soi, pour faire disparaitre la singularité des variables hydrodynamiques. La microrégion (le voisinage de la ligne de contact) est résolue numériquement et analytiquement pour de faibles pentes de l'interface liquide-vapeur. Les échelles de longueur caractéristiques du problème sont identifiées et la nature multi-échelle du phénomène est prise en compte. Les études paramétriques effectuées révèlent le rôle de la résistance thermique de l'interface vapeur-liquide, de la longueur de glissement, du terme thermocapillaire, du recul de vapeur et ainsi que des forces de surface. Une extension de l'approximation de lubrification pour de pentes élevées de l'interface gaz-liquide à l'évaporation est discutée. Dans la seconde partie de la thèse, le modèle précédemment établi pour la microrégion est couplé à des simulations numériques de la croissance d'une bulle de vapeur. Le départ de la bulle de vapeur de la paroi chauffante pendant l'ébullition a également été étudiée. Il a été proposé dans la thèse, que sous des charges thermiques élevées, l'augmentation de l'angle de contact apparent provoque l'étalement de la bulle de vapeur sur la paroi chauffante. Ce phénomène peut conduire, au séchage de la paroi observé pendant la crise d'ébullition
This thesis theoretically investigates the transport processes in the vicinity of the triple gas-liquid-solid contact line and its impact on macroscopic evaporation. In the first part of the thesis, the hydrodynamics close to the contact line at partial wetting is studied. Specifically, evaporation into the atmosphere of pure vapor driven by heating of the substrate is considered. The question of singularity relaxation is addressed. The main finding of the thesis is that the Kelvin effect (dependence of saturation temperature on pressure) is sufficient by itself to relax the hydrodynamic contact line singularity. The proposed microregion (the contact line vicinity) model for small interface slopes is solved numerically. Asymptotic solutions are found for some specific cases. The governing length scales of the problem are identified and the multiscale nature of the phenomenon is addressed. Parametric studies revealing the role of the thermal resistance of vapor-liquid interface, slip length, thermocapillary term, the vapor recoil and surface forces are also performed. An extension of the lubrication approximation for high slopes of the gas-liquid interface at evaporation is discussed. In the second part of the thesis, the previously established microregion model is coupled to a simplified single vapor bubble growth numerical simulation. The bubble departure from the heater at boiling is also studied. It was proposed in the thesis, that under high heat loads, the increase of the apparent contact angle causes the vapor bubble to spread over the heated substrate. Such a behavior may cause the heater dry-out that occurs during the boiling crisis

Thesis: Ph. D., Massachusetts Institute of Technology, Department of Mechanical Engineering, 2016.
Cataloged from PDF version of thesis.
Includes bibliographical references (pages 119-123).
Heat dissipation is a critical limitation in a range of electronic devices including microprocessors, solar cells, laser diodes and power amplifiers. The most demanding devices require dissipation of heat fluxes in excess of 1 kW/cm2 with heat transfer coefficients more than 30 W/cm 2K. Advanced thermal management solutions using phase change heat transfer are the most promising approach to address these challenges, yet current solutions are limited due to the combination of heat flux, thermal resistance, size and flow stability. This thesis reports the design, fabrication and experimental characterization for an evaporation device with a nanoporous membrane for high heat flux dissipation. Evaporation in the thin film regime is achieved using nanopores with reduced liquid film thicknesses while liquid pumping is enhanced using the capillary pressure of the 120 nm pores. The membrane is mechanically supported by ridges that form liquid supply channels and also serve as a heat conduction path to the evaporating meniscus at the surface of the membrane. The combination of high capillarity pores with high permeability channels facilitates theoretical critical heat fluxes over 2 kW/cm2 and heat transfer coefficients over 100 W/cm2K. Proof-of-concept devices were fabricated using a two-wafer stack consisting of a bonded silicon-on-insulator (SOI) wafer to a silicon wafer. Pores with diameters 110 - 130 nm were defined with interference lithography and etched in the SOI. Liquid supply microchannels were etched on a silicon wafer and the two wafers were fusion bonded together to form a monolithic evaporator. Once bonded, the membrane was released by etching through the backside of the SOI. Finally, platinum heaters and Resistive Temperature Detectors (RTDs) were deposited by e-beam evaporation and liftoff to heat the sample and measure the device temperature during experiments, respectively. Samples were experimentally characterized in a custom environmental chamber for comparison to the model using R245fa, methanol, pentane, water and isopropyl alcohol as working fluids. A comparison of the results with different working fluids demonstrates that transport at the liquid-vapor interface is the dominant thermal resistance in the system, suggesting a figure of merit: ... The highest heat flux recorded was with pentane at ... and the highest heat transfer coefficient recorded was with ... not including the substrate resistance. However, the samples were observed to clog with soluble, nonvolatile contaminants which limited operation to several minutes. The clogging behavior was captured in a mass diffusion model and a new configuration was suggested which is resistant to clogging. Evaporation from nanopores represents a new paradigm in phase change cooling with a figure of merit that favors high volatility, low surface tension fluids rather than water. The models and experimental results validate the functionality and understanding of the proposed approach and provide recommendations for enhancements in performance and understanding as well as strategies for resistance to clogging. This work demonstrates that nanoporous membranes have the potential for ultra-high heat flux dissipation to address next generation thermal management needs.
by Daniel Frank Hanks.
Ph. D.

La methode de co-evaporation dans des conditions d'ultravide s'avere la mieux adaptee au depot de couches tres minces. Elles est employee ici pour la fabrication d'oxydes supraconducteurs a haute temperature critique du type (bio)#2sr#2ca#n##1cu#no#2#(#n#+#1#) qui ont des specificites structurales comme par exemple leur caractere lamellaire et la versatilite de leurs structures. Les films minces de ces oxydes ont des applications en electronique et en instrumentation (interconnections, squid, composants hyperfrequence, circuits logiques. . . ). Apres une presentation des oxydes supraconducteurs (en particulier les composes a base de bismuth) et des differentes methodes de depot, le systeme de co-evaporation adapte a l'utilisation d'oxygene en ultravide et les techniques d'analyses sont decrits. Nous montrons l'efficacite d'une source a plasma d'oxygene que nous avons mise au point sur l'oxydation du cuivre. Nous mettons en evidence et nous etudions en detail les problemes d'incorporation du bismuth lors des depots. Nous discutons ensuite les proprietes des films supraconducteurs obtenus. L'utilisation d'oxygene dans le systeme ultravide, les perspectives ouvertes par les evaporations en flux simultanes et en flux modules sont egalement discutees

Le processus d'étalement et d'évaporation d'une goutte de suspensions de particules sur une surface solide est très intéressant permettant la formation de motifs. Une étude expérimentale à été effectuée avec du sang total humain et avec des fluides purs dans une chambre sous atmosphère contrôlée en humidité relative. Pour des angles de contact faible, le processus d'étalement/évaporation peut être divisé en deux régimes. Un premier régime rapide gouverné par un équilibre entre les forces visqueuses et les forces capillaires et un deuxième régime plus lent dominé par la cinétique d'évaporation. Nous montrons que les bio-colloïdes jouent un rôle significatif sur la dynamique de la ligne de contact. La vitesse moyenne de la ligne de contact suit la même dynamique d'étalement que le modèle de Tanner, où le temps d'étalement et les paramètres géométriques de la goutte sont fonctions de l'humidité. Dans cette étude, nous montrons que l'humidité relative influence les paramètres géométriques de la goutte et par conséquent le motif a la fin du processus d'évaporation. Un modèle purement diffusif pur a été obtenu dont le diamètre de mouillage et l'angle de contact sont fonction de l'humidité. Pour l'analyse morphologique des motifs de craquelures, une méthode de segmentation manuelle a été utilisée comme une méthode de référence pour la validation de la méthode de segmentation automatique développée dans "iBlood". Par cette méthode, nous montrons que la cinétique d'évaporation influence la distribution structurelle et morphologique des cellules de forme trapézoïdale, et par conséquent, l'espacement des fractures moyenne finale
Spreading/evaporation process of droplets over solid surfaces is a fundamental process and a wide research field because of number of applications in printing, micro-electronics, DNA analysis and even in biomedical. This experimental work aims to investigate the effect of relative humidity on the contact line dynamics, on the evaporation dynamics and on the final pattern of a drop of whole human blood. The spreading of a pure fluid model that has the same physical properties as human blood was studied and compared to the blood. We showed that bio-colloids play significant effect on the dynamics of contact line and the pinning effect of the drop. For low contact angles, we showed that the spreading/evaporation process could be divided into two regimes. A fast first regime determined by a balance between viscous forces and capillary forces and a second slower regime dominated by the evaporation rate. Physical mechanisms that are responsible for the spreading enhancement are proposed and discussed. The average velocity of the contact line was found to follow the same behaviour as Tanner's model, where the spreading dynamics and geometrical parameters of the droplet are function of relative humidity. The experimental measurements are in a good agreement with the purely diffusive model where the equilibrium wetting radius and contact angle are function of relative humidity. For the morphological analysis of crack patterns, a manual segmentation method was used as a reference for the validation of the automatic developed segmentation method. We showed that the evaporation rate influences structural distribution of plaques in the corona region and hence, the mean crack spacing

Thesis: Ph. D., Massachusetts Institute of Technology, Department of Mechanical Engineering, 2016.
Cataloged from PDF version of thesis.
Includes bibliographical references (pages 129-133).
The generation of concentrated heat loads in advanced microprocessors, power amplifiers, and concentrated photovoltaics present significant thermal management challenge for defense, space and commercial applications. Liquid to vapor phase-change strategies are promising due to the high latent heat of vaporization of the working fluid. In particular, capillary pumped thin-film evaporation from micropillar wicks has received significant attention owing to advances in micro/nano-fabrication and the potential to dissipate high heat fluxes by increasing the evaporative area. Yet, predictive tools to design various wicking structures are not available due to limited understanding of the thermal-fluidic transport. This thesis reports experimental characterization and modeling of capillary-limited thin-film evaporation from micropillar wicks. We fabricated test devices and experimentally characterized the thermal performance of well-defined silicon micropillar wicks. The experiments were designed to investigate the capillary-limited dryout heat flux by ensuring pure thin-film evaporation in the absence of nucleate boiling. The tests were performed in a temperature controlled saturated vapor environment to accurately control the operating conditions. We also developed a unified semi-analytical thermal-fluidic model that incorporates the capillary pressure, permeability, and thermal resistance to help explain the experimental results. We then extended this work to study capillary-limited thin-film evaporation for dissipating extreme heat fluxes. We experimentally dissipated =6 kW/cm2 from a 640x620 [mu]m2 footprint, the largest heat flux reported to date when compared to past thin-film evaporation studies with similar size hotspots. We also demonstrated the potential of our devices to cool concurrent hotspots as well as when moderate uniform background heat flux was superposed with a hotspot. Our thermal management strategy is self-regulating and provides on-demand cooling unlike existing thermal management solutions. To gain insight into the fundamental physics of fluidic and thermal transport within the micropillar wick and explain the ultra-high heat fluxes demonstrated in our experiments, we developed a semi-analytical thermal-fluidic model that can predict the capillary-limited dryout heat flux via thin-film evaporation. The model compares well with our experiments. The results of this investigation will assist to better understand the fluidic and thermal transport of thin liquid films in microstructured wicks during thin-film evaporation. These studies suggest that capillary-pumped thin-film evaporation is a promising thermal management strategy for the next generation of high performance electronics. The insights gained from this thesis can be used as guidelines to improve the design and optimize the heat transfer performance of wicking structures which are commonly used in phase-change based thermal management devices such as heat pipes, vapor chambers, and other closed-loop configurations.
by Solomon Adera.
Ph. D.

Thesis (M.S.) -- University of Maryland, College Park, 2007.
Thesis research directed by: Dept. of Mechanical Engineering. Title from t.p. of PDF. Includes bibliographical references. Published by UMI Dissertation Services, Ann Arbor, Mich. Also available in paper.

Study of liquid film dynamics is of significant importance to the understanding and control of various industrial processes that involve spray cooling (condensation), heating (boiling), coating, cleaning and lubrication. For instance, the critical heat flux (CHF) of boiling heat transfer is one of the key parameters ensuring the efficiency and safety of nuclear power plants under both operational and accident conditions, which occurs as the liquid layers (microlayer and macrolayer) near the heater wall lose their integrity. However, an experimental quantification of thin liquid film dynamics is not straightforward, since the measurement at micro-scale is a challenge, and further complicated by the chaotic nature of boiling process. The object of present study is to develop experimental methods for the diagnosis of liquid film dynamics, and to obtain data for the film instability under various conditions. A dedicated test facility was designed and constructed where micro conductive probes and confocal optical sensors were used to measure the thickness and dynamic characteristics of a thin liquid film on various heater surfaces, while a high speed camera was used to get visual observation. Extensive tests were performed to calibrate and verify the two thickness measuring systems. The micro conductive measuring system was proven to have a high reliability and repeatability with maximum system error less than 5µm, while the optical measuring system is capable of recording the film dynamics with spatial resolution of less than 1 mm. The simultaneous measurement on the same liquid film shows that the two techniques are in a good agreement with respect to accuracy, but the optical sensors have a much higher acquisition rate up to 30 kHz, which are more suitable for rapid process. The confocal optical sensors were therefore employed to measure the dynamic thickness of liquid films (ethanol, hexane and water) evaporating on various horizontal heater surfaces (aluminum, copper, silicon, stainless steel and titanium) to investigate the influences of heat flux, the surface and liquid properties on the film instability and the critical thickness. The critical thickness of water film evaporating on various surfaces was measured in the range of 60-150 mm, increasing with the increased contact angle or increased heat flux (evaporating rate) and also with the decreased thermal conductivity of the heater material. The data suggest the conjugate heat transfer nature of the evaporating liquid film dynamics at higher heat fluxes of interest to boiling and burnout. In the case of hexane on the aged titanium surface with contact angle of ~3o, the liquid film is found resilient to rupture, with film oscillations at relatively large amplitude ensuing as the averaged film thickness decreases below 15 µm. To interpret our experimental findings on liquid film evolution and its critical thickness at rupture, a theoretical analysis is also performed to analyze the dynamics of liquid films evaporating on heater surfaces. While the influences of liquid properties, heat flux, and thermal conductivity of heater surface are captured by the simulation of the lubrication theory, influence of the wettability is considered via a minimum free energy criterion. The thinning processes of the liquid films are generally captured by the simulation of the lubrication theory. For the case with ideally uniform heat flux over the heater surface, the instability of the liquid film occurs at the thickness level of tens micro meters, while for the case of non-uniform heating, the critical thicknesses for the film rupture are closer to  the experimental data but still underestimated by the lubrication theory simulation. By introducing the minimum free energy criterion to considering the influence of surface wettability, the obtained critical thicknesses have a good agreement with the experimental ones for both titanium and copper surfaces, with a maximum deviation less than ±10%. The simulations also explain why the critical thickness on a copper surface is thinner than that on a titanium surface. It is because the good thermal conductivity of copper surface leads to uniform temperature distribution on the heat surface, which is responsible for the resilience of the liquid film to rupture. A silicon wafer with an artificial cavity fabricated by Micro Electronic Mechanical System (MEMS) technology was used as a heater to investigate the dynamics of a single bubble in both a thick and thin liquid layer under low heat flux (<60 kW/m2). The maximum departure diameter of an isolated bubble in a thick liquid film was measured to be 3.2 mm which is well predicted by the Fritz equation. However, in a thin liquid layer with its thickness less than the bubble departure diameter, the bubble was stuck on the heater surface with a dry spot beneath. A threshold thickness of the liquid film which enables the dry spot rewettable was obtained, and its value linearly increases with increasing heat flux. In addition, another test section was designed to achieve a constant liquid film flow on a titanium nano-heater surface which helps to successfully carry boiling in the liquid film from low heat flux until CHF. Again, the confocal optical sensor was employed to measure the dynamics of the liquid film on the heater surface under varied heat flux conditions.  A statistical analysis of the measured thickness signals that emerge in a certain period indicates three distinct liquid film thickness ranges: 0~50 µm as microlayer, 50~500 µm as macrolayer, 500~2500 µm as bulk layer. With increasing heat flux, the bulk layer disappears, and then the macrolayer gradually decreases to ~105 µm, beyond which instability of the liquid film may lose its integrity and CHF occurs. In addition, the high-speed camera was applied to directly visualize and record the bubbles dynamics and liquid film evolution. Dry spots were observed under some bubbles occasionally from 313 kW/m2 until CHF with the maximum occupation fraction within 5%.  A dry spot was rewetted either by liquid receding after the rupture of a bubble or by the liquid spreading from bubbles’ growth in the vicinity. This implies that the bubbles’ behavior (growth and rupture) and their interactions in particular are of paramount importance to the integrity of liquid film under nucleate boiling regime.
QC 20111205
VR-2005-5729, MSWI

The evaporation phenomenon of a liquid droplet was investigated by using microfabricated heaters. All 32 microheaters were designed to have the same resistance. Gold microheaters worked both as temperature indicators and as heaters. The first experiment was performed under a constant voltage mode to investigate the temperature and heat flux variation of the heated surface by the evaporating droplet. The second experiment was performed under constant temperature mode to investigate the spatial and temporal heat flux variation of the constant temperature heater surface by the evaporating droplet heater. Droplet evaporation was recorded with a CCD camera. Experimental data showed temperature and heat flux variations inside and outside of the droplet with respect to time and radial position from the center of the droplet by tomographic deconvolution.

Thesis: S.M., Massachusetts Institute of Technology, Department of Mechanical Engineering, 2014.
Cataloged from PDF version of thesis.
Includes bibliographical references (pages 44-46).
Heat dissipation is a limiting factor in the performance of integrated circuits, power electronics and laser diodes. State-of-the-art solutions typically use air-cooled heat sinks, which have limited performance owing to the use of air. One of the promising approaches to address these thermal management needs is liquid vapor phase-change. In this thesis, we present a study into the design and modeling of a cooling device based on thin film evaporation from a nanoporous membrane supported on microchannels. The concept utilizes the capillary pressure generated by the small pores to drive the liquid flow and largely reduces the viscous loss due to the thinness of the membrane. The interfacial transport has been re-investigated where we use the moment method to solve the Boltzmann Transport Equation. The pore-level transport has been modeled coupling liquid transport, vapor transport and the interfacial balance. The interfacial transport inside the pore also serves as a boundary condition for the device-level model. The heat transfer and pressure drop performance have been modeled and design guidelines are provided for the membrane-based cooling system. The optimized cooling device is able to dissipate 1 kW/cm² heat flux with a temperature rise less than 30 K from the vapor side. Future work will focus on more fundamental understanding of the mass and energy accommodation at the liquid vapor interface.
by Zhengmao Lu.
S.M.

This work was written as a cumulative doctoral thesis based on reviewed publications. Climate projections are mainly based on the results of numeric simulations from global or regional climate models. Up to now processes between atmosphere and land surface are only rudimentarily known. This causes one of the major uncertainties in existing models. In order to reduce parameterisation uncertainties and to find a reasonable description of sub grid heterogeneities, the determination and evaluation of parameterisation schemes for modelling require as many datasets from different spatial scales as possible. This work contributes to this topic by implying different datasets from different platforms. Its objective was to analyse the spatial heterogeneity of land surface parameters and energy flux densities obtained from both satellite observations with different spatial and temporal resolutions and in-situ measurements. The investigations were carried out for two target areas in Germany. First, satellite data for the years 2002 and 2003 were analysed and validated from the LITFASS-area (Lindenberg Inhomogeneous Terrain - Fluxes between Atmosphere and Surface: a longterm Study). Second, the data from the experimental field sites of the FLUXNET cluster around Tharandt from the years 2006 and 2007 were used to determine the NDVI (Normalised Difference Vegetation Index for identifying vegetated areas and their "condition"). The core of the study was the determination of land surface characteristics and hence radiant and energy flux densities (net radiation, soil heat flux, sensible and latent heat flux) using the three optical satellite sensors ETM+ (Enhanced Thematic Mapper), MODIS (Moderate Resolution Imaging Spektroradiometer) and AVHRR 3 (Advanced Very High Resolution Radiometer) with different spatial (30 m – 1 km) and temporal (1 day – 16 days) resolution. Different sensor characteristics and different data sets for land use classifications can both lead to deviations of the resultant energy fluxes between the sensors. Thus, sensor differences were quantified, sensor adaptation methods were implemented and a quality analysis for land use classifications was performed. The result is then a single parameterisation scheme that allows for the determination of the energy fluxes from all three different sensors. The main focus was the derivation of the latent heat flux (L.E) using the Penman-Monteith (P-M) approach. Satellite data provide measurements of spectral reflectance and surface temperatures. The P-M approach requires further surface parameters not offered by satellite data. These parameters include the NDVI, Leaf Area Index (LAI), wind speed, relative humidity, vegetation height and roughness length, for example. They were derived indirectly from the given satellite- or in-situ measurements. If no data were available so called default values from literature were taken. The quality of these parameters strongly influenced the exactness of the radiant- and energy fluxes. Sensitivity studies showed that NDVI is one of the most important parameters for determination of evaporation. In contrast it could be shown, that the parameters as vegetation height and measurement height have only minor influence on L.E, which justifies the use of default values for these parameters. Due to the key role of NDVI a field study was carried out investigating the spatial variability and sensitivity of NDVI above five different land use types (winter wheat, corn, grass, beech and spruce). Methods to determine this parameter not only from space (spectral), but also from in-situ tower measurements (broadband) and spectrometer data (spectral) were compared. The best agreement between the methods was found for winter wheat and grass measurements in 2006. For these land use types the results differed by less than 10 % and 15 %, respectively. Larger differences were obtained for the forest measurements. The correlation between the daily MODIS-NDVI data and the in-situ NDVI inferred from the spectrometer and the broadband measurements were r=0.67 and r=0.51, respectively. Subsequently, spatial variability of land surface parameters and fluxes were analysed. The several spatial resolutions of the satellite sensors can be used to describe subscale heterogeneity from one scale to the other and to study the effects of spatial averaging. Therefore land use dependent parameters and fluxes were investigated to find typical distribution patterns of land surface properties and energy fluxes. Implying the distribution patterns found here for albedo and NDVI from ETM+ data in models has high potential to calculate representative energy flux distributions on a coarser scale. The distribution patterns were expressed as probability density functions (PDFs). First results of applying PDFs of albedo, NDVI, relative humidity, and wind speed to the L.E computation are encouraging, and they show the high potential of this method. Summing up, the method of satellite based surface parameter- and energy flux determination has been shown to work reliably on different temporal and spatial scales. The data are useful for detailed analyses of spatial variability of a landscape and for the description of sub grid heterogeneity, as it is needed in model applications. Their usability as input parameters for modelling on different scales is the second important result of this work. The derived vegetation parameters, e.g. LAI and plant cover, possess realistic values and were used as model input for the Lokalmodell of the German Weather Service. This significantly improved the model results for L.E. Additionally, thermal parameter fields, e.g. surface temperature from ETM+ with 30 m spatial resolution, were used as input for SVAT-modelling (Soil-Vegetation-Atmosphere-Transfer scheme). Thus, more realistic L.E results were obtained, providing highly resolved areal information
Die vorliegende Arbeit wurde auf der Grundlage begutachteter Publikationen als kumulative Dissertation verfasst. Klimaprognosen basieren im Allgemeinen auf den Ergebnissen numerischer Simulationen mit globalen oder regionalen Klimamodellen. Eine der entscheidenden Unsicherheiten bestehender Modelle liegt in dem noch unzureichenden Verständnis von Wechselwirkungsprozessen zwischen der Atmosphäre und Landoberflächen und dem daraus folgenden Fehlen entsprechender Parametrisierungen. Um das Problem einer unsicheren Modell-Parametrisierung aufzugreifen und zum Beispiel subskalige Heterogenität in einer Art und Weise zu beschreiben, dass sie für Modelle nutzbar wird, werden für die Bestimmung und Evaluierung von Modell-Parametrisierungsansätzen so viele Datensätze wie möglich benötigt. Die Arbeit trägt zu diesem Thema durch die Verwendung verschiedener Datensätze unterschiedlicher Plattformen bei. Ziel der Studie war es, aus Satellitendaten verschiedener räumlicher und zeitlicher Auflösung sowie aus in-situ Daten die räumliche Heterogenität von Landoberflächenparametern und Energieflussdichten zu bestimmen. Die Untersuchungen wurden für zwei Zielgebiete in Deutschland durchgeführt. Für das LITFASS-Gebiet (Lindenberg Inhomogeneous Terrain - Fluxes between Atmosphere and Surface: a longterm Study) wurden Satellitendaten der Jahre 2002 und 2003 untersucht und validiert. Zusätzlich wurde im Rahmen dieser Arbeit eine NDVI-Studie (Normalisierter Differenzen Vegetations Index: Maß zur Detektierung von Vegetationflächen, deren Vitalität und Dichte) auf den Testflächen des FLUXNET Clusters um Tharandt in den Jahren 2006 und 2007 realisiert. Die Grundlage der Arbeit bildete die Bestimmung von Landoberflächeneigenschaften und daraus resultierenden Energieflüssen, auf Basis dreier optischer Sensoren (ETM+ (Enhanced Thematic Mapper), MODIS (Moderate Resolution Imaging Spectroradiometer) und AVHRR 3 (Advanced Very High Resolution Radiometer)) mit unterschiedlichen räumlichen (30 m – 1 km) und zeitlichen (1 – 16 Tage) Auflösungen. Unterschiedliche Sensorcharakteristiken, sowie die Verwendung verschiedener, zum Teil ungenauer Datensätze zur Landnutzungsklassifikation führen zu Abweichungen in den Ergebnissen der einzelnen Sensoren. Durch die Quantifizierung der Sensorunterschiede, die Anpassung der Ergebnisse der Sensoren aneinander und eine Qualitätsanalyse von verschiedenen Landnutzungsklassifikationen, wurde eine Basis für eine vergleichbare Parametrisierung der Oberflächenparameter und damit auch für die daraus berechneten Energieflüsse geschaffen. Der Schwerpunkt lag dabei auf der Bestimmung des latenten Wärmestromes (L.E) mit Hilfe des Penman-Monteith Ansatzes (P-M). Satellitendaten liefern Messwerte der spektralen Reflexion und der Oberflächentemperatur. Die P-M Gleichung erfordert weitere Oberflächenparameter wie zum Beispiel den NDVI, den Blattflächenindex (LAI), die Windgeschwindigkeit, die relative Luftfeuchte, die Vegetationshöhe oder die Rauhigkeitslänge, die jedoch aus den Satellitendaten nicht bestimmt werden können. Sie müssen indirekt aus den oben genannten Messgrößen der Satelliten oder aus in-situ Messungen abgeleitet werden. Stehen auch aus diesen Quellen keine Daten zur Verfügung, können sogenannte Standard- (Default-) Werte aus der Literatur verwendet werden. Die Qualität dieser Parameter hat einen großen Einfluss auf die Bestimmung der Strahlungs- und Energieflüsse. Sensitivitätsstudien im Rahmen der Arbeit zeigen die Bedeutung des NDVI als einen der wichtigsten Parameter in der Verdunstungsbestimmung nach P-M. Im Gegensatz dazu wurde deutlich, dass z. B. die Vegetationshöhe und die Messhöhe einen relativ kleinen Einfluss auf L.E haben, so dass für diese Parameter die Verwendung von Standardwerten gerechtfertigt ist. Aufgrund der Schlüsselrolle, welche der NDVI in der Bestimmung der Verdunstung einnimmt, wurden im Rahmen einer Feldstudie Untersuchungen des NDVI über fünf verschiedenen Landnutzungstypen (Winterweizen, Mais, Gras, Buche und Fichte) hinsichtlich seiner räumlichen Variabilität und Sensitivität, unternommen. Dabei wurden verschiedene Bestimmungsmethoden getestet, in welchen der NDVI nicht nur aus Satellitendaten (spektral), sondern auch aus in-situ Turmmessungen (breitbandig) und Spekrometermessungen (spektral) ermittelt wird. Die besten Übereinstimmungen der Ergebnisse wurden dabei für Winterweizen und Gras für das Jahr 2006 gefunden. Für diese Landnutzungstypen betrugen die Maximaldifferenzen aus den drei Methoden jeweils 10 beziehungsweise 15 %. Deutlichere Differenzen ließen sich für die Forstflächen verzeichnen. Die Korrelation zwischen Satelliten- und Spektrometermessung betrug r=0.67. Für Satelliten- und Turmmessungen ergab sich ein Wert von r=0.5. Basierend auf den beschriebenen Vorarbeiten wurde die räumliche Variabilität von Landoberflächenparametern und Flüssen untersucht. Die unterschiedlichen räumlichen Auflösungen der Satelliten können genutzt werden, um zum einen die subskalige Heterogenität zu beschreiben, aber auch, um den Effekt räumlicher Mittelungsverfahren zu testen. Dafür wurden Parameter und Energieflüsse in Abhängigkeit der Landnutzungsklasse untersucht, um typische Verteilungsmuster dieser Größen zu finden. Die Verwendung der Verteilungsmuster (in Form von Wahrscheinlichkeitsdichteverteilungen – PDFs), die für die Albedo und den NDVI aus ETM+ Daten gefunden wurden, bietet ein hohes Potential als Modellinput, um repräsentative PDFs der Energieflüsse auf gröberen Skalen zu erhalten. Die ersten Ergebnisse in der Verwendung der PDFs von Albedo, NDVI, relativer Luftfeuchtigkeit und Windgeschwindigkeit für die Bestimmung von L.E waren sehr ermutigend und zeigten das hohe Potential der Methode. Zusammenfassend lässt sich feststellen, dass die Methode der Ableitung von Oberflächenparametern und Energieflüssen aus Satellitendaten zuverlässige Daten auf verschiedenen zeitlichen und räumlichen Skalen liefert. Die Daten sind für eine detaillierte Analyse der räumlichen Variabilität der Landschaft und für die Beschreibung der subskaligen Heterogenität, wie sie oft in Modellanwendungen benötigt wird, geeignet. Ihre Nutzbarkeit als Inputparameter in Modellen auf verschiedenen Skalen ist das zweite wichtige Ergebnis der Arbeit. Aus Satellitendaten abgeleitete Vegetationsparameter wie der LAI oder die Pflanzenbedeckung liefern realistische Ergebnisse, die zum Beispiel als Modellinput in das Lokalmodell des Deutschen Wetterdienstes implementiert werden konnten und die Modellergebnisse von L.E signifikant verbessert haben. Aber auch thermale Parameter, wie beispielsweise die Oberflächentemperatur aus ETM+ Daten in 30 m Auflösung, wurden als Eingabeparameter eines Soil-Vegetation-Atmosphere-Transfer-Modells (SVAT) verwendet. Dadurch erhält man realistischere Ergebnisse für L.E, die hochaufgelöste Flächeninformationen bieten

This work was written as a cumulative doctoral thesis based on reviewed publications. Climate projections are mainly based on the results of numeric simulations from global or regional climate models. Up to now processes between atmosphere and land surface are only rudimentarily known. This causes one of the major uncertainties in existing models. In order to reduce parameterisation uncertainties and to find a reasonable description of sub grid heterogeneities, the determination and evaluation of parameterisation schemes for modelling require as many datasets from different spatial scales as possible. This work contributes to this topic by implying different datasets from different platforms. Its objective was to analyse the spatial heterogeneity of land surface parameters and energy flux densities obtained from both satellite observations with different spatial and temporal resolutions and in-situ measurements. The investigations were carried out for two target areas in Germany. First, satellite data for the years 2002 and 2003 were analysed and validated from the LITFASS-area (Lindenberg Inhomogeneous Terrain - Fluxes between Atmosphere and Surface: a longterm Study). Second, the data from the experimental field sites of the FLUXNET cluster around Tharandt from the years 2006 and 2007 were used to determine the NDVI (Normalised Difference Vegetation Index for identifying vegetated areas and their "condition"). The core of the study was the determination of land surface characteristics and hence radiant and energy flux densities (net radiation, soil heat flux, sensible and latent heat flux) using the three optical satellite sensors ETM+ (Enhanced Thematic Mapper), MODIS (Moderate Resolution Imaging Spektroradiometer) and AVHRR 3 (Advanced Very High Resolution Radiometer) with different spatial (30 m – 1 km) and temporal (1 day – 16 days) resolution. Different sensor characteristics and different data sets for land use classifications can both lead to deviations of the resultant energy fluxes between the sensors. Thus, sensor differences were quantified, sensor adaptation methods were implemented and a quality analysis for land use classifications was performed. The result is then a single parameterisation scheme that allows for the determination of the energy fluxes from all three different sensors. The main focus was the derivation of the latent heat flux (L.E) using the Penman-Monteith (P-M) approach. Satellite data provide measurements of spectral reflectance and surface temperatures. The P-M approach requires further surface parameters not offered by satellite data. These parameters include the NDVI, Leaf Area Index (LAI), wind speed, relative humidity, vegetation height and roughness length, for example. They were derived indirectly from the given satellite- or in-situ measurements. If no data were available so called default values from literature were taken. The quality of these parameters strongly influenced the exactness of the radiant- and energy fluxes. Sensitivity studies showed that NDVI is one of the most important parameters for determination of evaporation. In contrast it could be shown, that the parameters as vegetation height and measurement height have only minor influence on L.E, which justifies the use of default values for these parameters. Due to the key role of NDVI a field study was carried out investigating the spatial variability and sensitivity of NDVI above five different land use types (winter wheat, corn, grass, beech and spruce). Methods to determine this parameter not only from space (spectral), but also from in-situ tower measurements (broadband) and spectrometer data (spectral) were compared. The best agreement between the methods was found for winter wheat and grass measurements in 2006. For these land use types the results differed by less than 10 % and 15 %, respectively. Larger differences were obtained for the forest measurements. The correlation between the daily MODIS-NDVI data and the in-situ NDVI inferred from the spectrometer and the broadband measurements were r=0.67 and r=0.51, respectively. Subsequently, spatial variability of land surface parameters and fluxes were analysed. The several spatial resolutions of the satellite sensors can be used to describe subscale heterogeneity from one scale to the other and to study the effects of spatial averaging. Therefore land use dependent parameters and fluxes were investigated to find typical distribution patterns of land surface properties and energy fluxes. Implying the distribution patterns found here for albedo and NDVI from ETM+ data in models has high potential to calculate representative energy flux distributions on a coarser scale. The distribution patterns were expressed as probability density functions (PDFs). First results of applying PDFs of albedo, NDVI, relative humidity, and wind speed to the L.E computation are encouraging, and they show the high potential of this method. Summing up, the method of satellite based surface parameter- and energy flux determination has been shown to work reliably on different temporal and spatial scales. The data are useful for detailed analyses of spatial variability of a landscape and for the description of sub grid heterogeneity, as it is needed in model applications. Their usability as input parameters for modelling on different scales is the second important result of this work. The derived vegetation parameters, e.g. LAI and plant cover, possess realistic values and were used as model input for the Lokalmodell of the German Weather Service. This significantly improved the model results for L.E. Additionally, thermal parameter fields, e.g. surface temperature from ETM+ with 30 m spatial resolution, were used as input for SVAT-modelling (Soil-Vegetation-Atmosphere-Transfer scheme). Thus, more realistic L.E results were obtained, providing highly resolved areal information
Die vorliegende Arbeit wurde auf der Grundlage begutachteter Publikationen als kumulative Dissertation verfasst. Klimaprognosen basieren im Allgemeinen auf den Ergebnissen numerischer Simulationen mit globalen oder regionalen Klimamodellen. Eine der entscheidenden Unsicherheiten bestehender Modelle liegt in dem noch unzureichenden Verständnis von Wechselwirkungsprozessen zwischen der Atmosphäre und Landoberflächen und dem daraus folgenden Fehlen entsprechender Parametrisierungen. Um das Problem einer unsicheren Modell-Parametrisierung aufzugreifen und zum Beispiel subskalige Heterogenität in einer Art und Weise zu beschreiben, dass sie für Modelle nutzbar wird, werden für die Bestimmung und Evaluierung von Modell-Parametrisierungsansätzen so viele Datensätze wie möglich benötigt. Die Arbeit trägt zu diesem Thema durch die Verwendung verschiedener Datensätze unterschiedlicher Plattformen bei. Ziel der Studie war es, aus Satellitendaten verschiedener räumlicher und zeitlicher Auflösung sowie aus in-situ Daten die räumliche Heterogenität von Landoberflächenparametern und Energieflussdichten zu bestimmen. Die Untersuchungen wurden für zwei Zielgebiete in Deutschland durchgeführt. Für das LITFASS-Gebiet (Lindenberg Inhomogeneous Terrain - Fluxes between Atmosphere and Surface: a longterm Study) wurden Satellitendaten der Jahre 2002 und 2003 untersucht und validiert. Zusätzlich wurde im Rahmen dieser Arbeit eine NDVI-Studie (Normalisierter Differenzen Vegetations Index: Maß zur Detektierung von Vegetationflächen, deren Vitalität und Dichte) auf den Testflächen des FLUXNET Clusters um Tharandt in den Jahren 2006 und 2007 realisiert. Die Grundlage der Arbeit bildete die Bestimmung von Landoberflächeneigenschaften und daraus resultierenden Energieflüssen, auf Basis dreier optischer Sensoren (ETM+ (Enhanced Thematic Mapper), MODIS (Moderate Resolution Imaging Spectroradiometer) und AVHRR 3 (Advanced Very High Resolution Radiometer)) mit unterschiedlichen räumlichen (30 m – 1 km) und zeitlichen (1 – 16 Tage) Auflösungen. Unterschiedliche Sensorcharakteristiken, sowie die Verwendung verschiedener, zum Teil ungenauer Datensätze zur Landnutzungsklassifikation führen zu Abweichungen in den Ergebnissen der einzelnen Sensoren. Durch die Quantifizierung der Sensorunterschiede, die Anpassung der Ergebnisse der Sensoren aneinander und eine Qualitätsanalyse von verschiedenen Landnutzungsklassifikationen, wurde eine Basis für eine vergleichbare Parametrisierung der Oberflächenparameter und damit auch für die daraus berechneten Energieflüsse geschaffen. Der Schwerpunkt lag dabei auf der Bestimmung des latenten Wärmestromes (L.E) mit Hilfe des Penman-Monteith Ansatzes (P-M). Satellitendaten liefern Messwerte der spektralen Reflexion und der Oberflächentemperatur. Die P-M Gleichung erfordert weitere Oberflächenparameter wie zum Beispiel den NDVI, den Blattflächenindex (LAI), die Windgeschwindigkeit, die relative Luftfeuchte, die Vegetationshöhe oder die Rauhigkeitslänge, die jedoch aus den Satellitendaten nicht bestimmt werden können. Sie müssen indirekt aus den oben genannten Messgrößen der Satelliten oder aus in-situ Messungen abgeleitet werden. Stehen auch aus diesen Quellen keine Daten zur Verfügung, können sogenannte Standard- (Default-) Werte aus der Literatur verwendet werden. Die Qualität dieser Parameter hat einen großen Einfluss auf die Bestimmung der Strahlungs- und Energieflüsse. Sensitivitätsstudien im Rahmen der Arbeit zeigen die Bedeutung des NDVI als einen der wichtigsten Parameter in der Verdunstungsbestimmung nach P-M. Im Gegensatz dazu wurde deutlich, dass z. B. die Vegetationshöhe und die Messhöhe einen relativ kleinen Einfluss auf L.E haben, so dass für diese Parameter die Verwendung von Standardwerten gerechtfertigt ist. Aufgrund der Schlüsselrolle, welche der NDVI in der Bestimmung der Verdunstung einnimmt, wurden im Rahmen einer Feldstudie Untersuchungen des NDVI über fünf verschiedenen Landnutzungstypen (Winterweizen, Mais, Gras, Buche und Fichte) hinsichtlich seiner räumlichen Variabilität und Sensitivität, unternommen. Dabei wurden verschiedene Bestimmungsmethoden getestet, in welchen der NDVI nicht nur aus Satellitendaten (spektral), sondern auch aus in-situ Turmmessungen (breitbandig) und Spekrometermessungen (spektral) ermittelt wird. Die besten Übereinstimmungen der Ergebnisse wurden dabei für Winterweizen und Gras für das Jahr 2006 gefunden. Für diese Landnutzungstypen betrugen die Maximaldifferenzen aus den drei Methoden jeweils 10 beziehungsweise 15 %. Deutlichere Differenzen ließen sich für die Forstflächen verzeichnen. Die Korrelation zwischen Satelliten- und Spektrometermessung betrug r=0.67. Für Satelliten- und Turmmessungen ergab sich ein Wert von r=0.5. Basierend auf den beschriebenen Vorarbeiten wurde die räumliche Variabilität von Landoberflächenparametern und Flüssen untersucht. Die unterschiedlichen räumlichen Auflösungen der Satelliten können genutzt werden, um zum einen die subskalige Heterogenität zu beschreiben, aber auch, um den Effekt räumlicher Mittelungsverfahren zu testen. Dafür wurden Parameter und Energieflüsse in Abhängigkeit der Landnutzungsklasse untersucht, um typische Verteilungsmuster dieser Größen zu finden. Die Verwendung der Verteilungsmuster (in Form von Wahrscheinlichkeitsdichteverteilungen – PDFs), die für die Albedo und den NDVI aus ETM+ Daten gefunden wurden, bietet ein hohes Potential als Modellinput, um repräsentative PDFs der Energieflüsse auf gröberen Skalen zu erhalten. Die ersten Ergebnisse in der Verwendung der PDFs von Albedo, NDVI, relativer Luftfeuchtigkeit und Windgeschwindigkeit für die Bestimmung von L.E waren sehr ermutigend und zeigten das hohe Potential der Methode. Zusammenfassend lässt sich feststellen, dass die Methode der Ableitung von Oberflächenparametern und Energieflüssen aus Satellitendaten zuverlässige Daten auf verschiedenen zeitlichen und räumlichen Skalen liefert. Die Daten sind für eine detaillierte Analyse der räumlichen Variabilität der Landschaft und für die Beschreibung der subskaligen Heterogenität, wie sie oft in Modellanwendungen benötigt wird, geeignet. Ihre Nutzbarkeit als Inputparameter in Modellen auf verschiedenen Skalen ist das zweite wichtige Ergebnis der Arbeit. Aus Satellitendaten abgeleitete Vegetationsparameter wie der LAI oder die Pflanzenbedeckung liefern realistische Ergebnisse, die zum Beispiel als Modellinput in das Lokalmodell des Deutschen Wetterdienstes implementiert werden konnten und die Modellergebnisse von L.E signifikant verbessert haben. Aber auch thermale Parameter, wie beispielsweise die Oberflächentemperatur aus ETM+ Daten in 30 m Auflösung, wurden als Eingabeparameter eines Soil-Vegetation-Atmosphere-Transfer-Modells (SVAT) verwendet. Dadurch erhält man realistischere Ergebnisse für L.E, die hochaufgelöste Flächeninformationen bieten

Thesis (M. S.)--Mechanical Engineering, Georgia Institute of Technology, 2009.
Committee Chair: Abdel-Khalik, Said; Committee Member: Ammar, Mostafa H.; Committee Member: Ghiaasiaan, S. Mostafa; Committee Member: Hertel, Nolan E.; Committee Member: Liu, Yingjie.

To mitigate the global greenhouse gases (GHGs) emissions, carbon dioxide (CO2) capture and storage (CCS) has the potential to play a significant role for reaching mitigation target. Oxy-fuel combustion is a promising technology for CO2 capture in power plants. Advantages compared to CCS with the conventional combustion technology are: high combustion efficiency, flue gas volume reduction, low fuel consumption, near zero CO2 emission, and less nitrogen oxides (NOx) formation can be reached simultaneously by using the oxy-fuel combustion technology. However, knowledge gaps relating to large scale coal based and natural gas based power plants with CO2 capture still exist, such as combustors and boilers operating at higher temperatures and design of CO2 turbines and compressors. To apply the oxy-fuel combustion technology on power plants, much work is focused on the fundamental and feasibility study regarding combustion characterization, process and system analysis, and economic evaluation etc. Further studies from system perspective point of view are highlighted, such as the impact of operating conditions on system performance and on advanced cycle integrated with oxy-fuel combustion for CO2 capture. In this thesis, the characterization for flue gas recycle (FGR) was theoretically derived based on mass balance of combustion reactions, and system modeling was conducted by using a process simulator, Aspen Plus. Important parameters such as FGR rate and ratio, flue gas composition, and electrical efficiency etc. were analyzed and discussed based on different operational conditions. An advanced evaporative gas turbine (EvGT) cycle with oxy-fuel combustion for CO2 capture was also studied. Based on economic indicators such as specific investment cost (SIC), cost of electricity (COE), and cost of CO2avoidance (COA), economic performance was evaluated and compared among various system configurations. The system configurations include an EvGT cycle power plant without CO2 capture, an EvGT cycle power plant with chemical absorption for CO2 capture, and a combined cycle power plant. The study shows that FGR ratio is of importance, which has impact not only on heat transfer but also on mass transfer in the oxy-coal combustion process. Significant reduction in the amount of flue gas can be achieved due to the flue gas recycling, particularly for the system with more prior upstream recycle options. Although the recycle options have almost no effect on FGR ratio, flue gas flow rate, and system electrical efficiency, FGR options have significant effects on flue gas compositions, especially the concentrations of CO2 and H2O, and heat exchanger duties. In addition, oxygen purity and water/gas ratio, respectively, have an optimum value for an EvGT cycle power plant with oxy-fuel combustion. Oxygen purity of 97 mol% and water/gas ratio of 0.133 can be considered as the optimum values for the studied system. For optional operating conditions of flue gas recycling, the exhaust gas recycled after condensing (dry recycle) results in about 5 percentage points higher electrical efficiency and about 45 % more cooling water consumption comparing with the exhaust gas recycled before condensing (wet recycle). The direct costs of EvGT cycle with oxy-fuel combustion are a little higher than the direct costs of EvGT cycle with chemical absorption. However, as plant size is larger than 60 MW, even though the EvGT cycle with oxy-fuel combustion has a higher COE than the EvGT cycle with chemical absorption, the EvGT cycle with oxy-fuel combustion has a lower COA. Further, compared with others studies of natural gas combined cycle (NGCC), the EvGT system has a lower COE and COA than the NGCC system no matter which CO2 capture technology is integrated.
QC 20111123

Le travail que nous presentons dans ici, est consacre a deux objectifs de recherche principaux sur les supraconducteurs: tout d'abord l'elaboration des films minces supraconducteurs de yba#2cu#3o#7#-#x par evaporation reactive sequentielle, ensuite la definition et la caracterisation de deux dispositifs electroniques: un commutateur rapide non lineaire supraconducteur, fonctionnant a la temperature de l'azote liquide, et un transistor a flux de vortex

A cultura da cana é importante para o agronegócio brasileiro. O manejo adequado da palhada da cana-de-açúcar no solo pode influenciar no desenvolvimento e no crescimento da cultura, podendo resultar em ganhos ambientais e econômicos. Neste trabalho, objetivou-se avaliar o efeito do manejo com palha (CP) e sem palha (SP) na evapotranspiração de um canavial irrigado. Um experimento foi conduzido em canavial irrigado com dois tratamentos (mantendo-se 100% da palhada e sem cobertura de palha sobre o solo). Neste experimento, a evapotranspiração foi medida a partir do balanço hídrico \"in situ\" do solo (ETc-BHS) e pelo método da razão de Bowen (ETc-MRB) e a transpiração utilizando-se sensores de fluxo de seiva (FS). A evaporação do solo foi estimada com base em dados da cultura e meteorológicos. A modelagem da transpiração e de seus componentes baseou-se na abordagem no \"modelo da grande folha\" (\"big leaf model\") proposto por Monteith. Para tanto, foram definidos e avaliados submodelos para a simulação do saldo de radiação no dossel, da resistência aerodinâmica, da resistência de cobertura no ambiente da cultura. Ao longo do ciclo da cultura, análises biométricas foram realizadas a cada 30 dias e esses dados utilizados na parametrização do modelo. O desempenho dos submodelos e do modelo parametrizado de Penman-Monteith (PM) foi estatisticamente avaliados. Avalia-se um modelo para estimativa da evaporação em função da cobertura da palhada. O cultivo de cana-de-açúcar utilizando os resíduos vegetais favorece o incremento da transpiração, da evapotranspiração da cultura (ETc) e reduz a evaporação de água no solo. A palha influenciou as variáveis biométricas da cultura, como perfilhamento, massa fresca da parte aérea (MFPA) e massa seca da parte aérea (MSPA). Ao relacionar a transpiração estimada com o modelo de PM com a medida com FS, verifica-se um coeficiente de determinação (R2) de 0,72 e de 0,48 e um coeficiente angular (a) de 1,22 e 0,92, para a área manejada CP e SP, respectivamente. A evaporação de água no solo é maior na área SP, principalmente nos períodos com precipitações acima de 40 mm. Quando se relaciona a ETc-BHS e a ETc-MRB com a ETc estimada com base no método de PM, nota-se um R2 de 0,67 e 0,54 na área CP e 0,59 e 0,54 para a área SP, respectivamente. Independente do manejo adotado, aos 302 dias após o corte, o canavial apresentou ATR, POL e % de fibra ideal para a comercialização. O teor de clorofila na área manejada CP é superior ao manejado SP. O cultivo de cana-de-açúcar utilizando os resíduos vegetais favorece o incremento da transpiração, da ETc e reduz a evaporação de água no solo.
Cane culture is important for Brazilian agribusiness. Proper management of sugarcane trash in the soil can influence the development and growth of the crop, which can result in environmental and economic gains. The objective of this study was to evaluate the effect of trash management (CP) and without trash (SP) on the evapotranspiration of an irrigated sugar cane. An experiment was conducted in irrigated sugar cane fields with two treatments (100% of the trash and without trash cover on the soil). In this experiment, evapotranspiration was measured from soil water balance (ETc-BHS) and the Bowen ratio method (ETc-MRB) and transpiration using sap flow sensors (FS). Soil evaporation was estimated based on crop and meteorological data. The modeling of the transpiration and its components was based on the approach in the \"big leaf model\" proposed by Monteith. For that, submodels were defined and evaluated for the simulation of canopy radiation balance, aerodynamic resistance, and cover resistance in the culture environment. Throughout the crop cycle, biometric analyzes were performed every 30 days and these data were used in the parameterization of the model. The performance of the submodels and the parameterized model of Penman-Monteith (PM) was statistically evaluated. A model for estimation of evaporation is evaluated as a function of the trash cover. The cultivation of sugarcane using the vegetal residues favors the increase of transpiration, evapotranspiration of the crop (ETc) and reduces the evaporation of water in the soil. The trash influenced the biometric variables of the crop, such as tillering, fresh mass of aerial part (MFPA) and dry mass of aerial part (MSPA). When correlating the estimated transpiration with the PM model with the FS measurement, it is verified a determination coefficient (R2) of 0.72 and 0.48 and an angular coefficient (a) of 1.22 and 0.92, for the managed area CP and SP, respectively. The water evaporation in the soil is higher in the SP area, especially in periods with rainfall above 40 mm. When ETc-BHS and ETc-MRB are related to ETc estimated on the basis of the PM method, an R2 of 0.67 and 0.54 is observed in the CP area and 0.59 and 0.54 for the area SP, respectively. Regardless of the management adopted, at 302 days after cutting, the sugarcane plantation presented ATR, POL and% of fiber ideal for commercialization. The chlorophyll content in the managed area CP is superior to the managed SP. The cultivation of sugarcane using the vegetal residues favors the increase of the transpiration, ETc and reduces the evaporation of water in the soil.

The focus of the research is on investigations of various methods for obtaining sensible heat flux (H) for estimating evaporation. The key for this approach is the application of the shortened energy balance equation, and in the case of methods based on the Monin-Obukhov similarity theory (MOST), such as surface-layer scintillometrv and temperature variance with adjusted for stability using air temperature skewness, and surface renewal (SR), die iterative procedures. The application of the shortened energy balance requires that errors associated with measurement of net irradiance (Rnet) and soil heat flux (S) are kept to a minimum To this end. methodology for the calibration of net radiometers for both the infrared and short wave irradiances receive attention. A field study attempts to quantify the error in soil heat flux measurement for a mesic grassland. A standard, convenient and accurate method for calibrating net radiometers would assist in unravelling reasons for the perplexing lack of surface energy balance closure when employing the eddy covariance (EC) flux estimation method as well as improve on the accuracy of the energy balance residual method for estimating evaporation. A relatively inexpensive, accurate and quick laboratory method, based on physical theory, for non-steady radiative conditions above a large water-heated or water-cooled radiator containing circulated water, with surface-embedded thermocouples is used to obtain reproducible net radiometer calibration factor's for the infrared waveband for a wide range in net irradiance. When applied, the method would reduce error m the most important term of the shortened energy balance and assist in energy balance closure aspects of EC measurements. The SLS method, reliant on MOST, is used for estimating a really-averaged H for a mesic grassland for a 30-month period. Comparisons with EC measurements feature prominently in this unique study. These comparisons include using different MOST procedures and the influence of the Bowen ratio on SLS measurement: of if is investigated. Furthermore, since there are reports in the literature that the EC method may underestimate H and or latent energy flux (LE), resulting in the shortened energy balance not being closed, effort is devoted to this aspect. Other methods used for comparison purposes are the traditional Bowen ratio energy balance (BREB), SR, TV and ETo (grass reference) methods. The TV and SLS and/or EC measurements of H are compared above three contrasting canopy surfaces. It is shown that other high frequency air temperature-based methods, for example, for the first time the TV method with adjustment for skewness, may pave the way for evaporation stations from which real-time and sub-hourly estimates may be obtained relatively inexpensively. Another area of research that receives attention is the placement height of EC instruments above short-canopy surfaces and a spectral analysis of the vertical wind speed and some temperature measurement: for close-canopv placement heights. The SR method is used to estimate, for the first time, open-water evaporation. The ideal SR method applied above canopies is the most inexpensive micrometeorological method for estimating H, but the SR weighting factor a needs to be determined using EC and for this reason, the TV method with adjustment for skewness was investigated. Finally, a unique implementation of SR uses an iterative method for calculating H. A similar iterative procedure is applied for MOST and ETo calculations.
Thesis (DScAgric)-University of KwaZulu-Natal, Pietermaritzburg, 2010.

Evapotransipration by vegetation cover is an important component of the water budget and energy balance in any ecosystem. A key to more improved water management therefore lies in improving our understanding of evapotranspiration, the process that drives water use by plants. Estimations of the turbulent fluxes are required for various applications in micrometeorology, hydrology, environmental studies and agriculture. Numerous methods for estimation of turbulent fluxes have been developed and tested. Direct measurements of fluxes are usually achieved by the eddy covariance (EC) method, which is considered as the most reliable. However, the application of the EC method is often problematic. The necessary sensors for wind, temperature and humidity must respond very fast (resolution of 10 Hz or better) and at the same time must not show noticeable drift. This makes them delicate, expensive and difficult to calibrate among other problems associated with the method. Due to their ability to integrate atmospheric processes along a path length that may range between a few hundred metres to a few kilometres, optical methods based on the analysis of scintillation appear to be an alternative and possible supplement to classical micrometeorological methods such as the EC method, which may provide local fluxes typically at the scale of 100 m. The use of the scintillometry technique in surface flux measurements is therefore gaining in popularity. The accuracy of the measurements obtained by one method is judged by comparison of the measurements obtained by those of another method considered as the standard. For turbulent flux measurements, the EC method is taken as the standard method for the determination of sensible heat fluxes. This research presents the measurement of sensible heat fluxes using the surface layer scintillometer (SLS). The SLS system used has a dual-beam and a recommended path length of between 50 and 250 m. The method was tested against the EC method for different Bowen ratio (f3) values, as required by the theory, under different atmospheric stability conditions, as well as for different wind directions relative to the SLS beam path and slanting beam path orientation. Also presented is an analysis of the different forms of the Monin-Obukhov Similarity (MOST) functions used in micrometeorology and suggested by various authors, done by comparing the resulting sensible heat flux measured by the SLS method with the ones calculated through an iterative determination of the Monin-Obukhov parameters. A comparison of the structure function parameter of temperature (Ci ) corrected for fJ and those measured (using SLS) was carried out, with the results showing very good correspondence between the corrected and uncorrected ci values, indicating that not correcting for fJ for SLS measured ci does not result in significant error in the resulting ci values, and hence sensible heat flux estimates. A comparison of the sensible heat flux Fh obtained using EC and SLS methods for fJ < 0.6 and fJ > 0.6 followed and the results also show good correspondence between the values obtained using the EC and SLS methods, although the agreement is slightly improved for cases when fJ > 0.6. A sensitivity analysis indicates that both the ECand SLS-measurements of Fh are influenced by fJ values. A sensitivity analysis on the influence of fJ on Fh measurements by both the EC and SLS methods further indicates that the influence of fJ on Fh measurements is not large enough to warrant correcting Fh measurements for fJ . The F" measurements by the EC method appears to be influenced more by fJ especially for fJ values less than 0.74. A comparison of the various methods for computing the empirical similarity functions used by MOST was also carried out and the results show a significant difference in the Fh computed following the various methods suggested by different researchers. As for the agreement between the EC and SLS methods determination of Fh for the different atmospheric stability conditions, there seems to be a better agreement in the Fh measurements as noted by correlation coefficients closer to 1 and greater tvalues obtained during unstable atmospheric conditions in the colder months of June and August while reduced agreement in the values is recorded in the warmer summer period from November to December. Also noted is a slight difference in the EC measurements compared to the SLS measurement of F". The difference in the measurements is noticed for unstable atmospheric conditions. Also noted is that EC and SLS measurements of Fh differ slightly when the atmospheric condition is nearneutral. However the agreement between the Fh values measured by the two measurement methods is still good. was set up in an inclined position, with the receiver set at 0.68 m above the ground level and transmitter at 1.68 m, resulting in an effective height difference of 1.00 m. There was generally good agreement in the 2-min measurements of F" by the two methods for the SLS set up in inclined position, with the 30-min data resulting in even better agreements. The findings confirm that the SLS set up does not impair its performance in measuring sensible heat fluxes. This also shows that the SLS would also work well in non-ideal (heterogeneous) conditions which the inclined optical beam path mimics. For those days when wind direction was mainly approximately perpendicular to the beam, the F" values obtained by SLS and EC methods are more in agreement than when the wind direction was either irregular or parallel to the SLS beam path. Wind speed also seems to influence the F" estimates by the two methods since the agreement in the Fh values obtained by the two methods is greater when wind speed is higher compared to times of the day when the wind speed is reduced. The atmospheric stability influences the peak position of footprint with the peak footprint position being further from the measurement point when the atmospheric stability condition is closer to stable as denoted by the Obukhov length of -5 and closer to the measurement point for convectively unstable atmospheric conditions as shown by the Obukhov length of -30. Also shown is that a larger fetch is required when the atmosphere is convectively unstable as indicated by the contours plotted on top of the footprint plots. In general, there seems to be very good agreement in the sensible heat flux values obtained by the two methods, especially since SLS offers areal-averaged sensible heat flux measurements compared to the EC method which basically provides a point measurement. The SLS method therefore offers a better alternative for obtaining sensible heat flux from larger and heterogeneous area - although to a limit of250 m since beyond 250 m, the method suffers from a saturation problem.
Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2007.

The high population growth rate and rapid urbanization that the world is experiencing today has aggravated the competition for the already scarce resource ¡V water ¡V between the agricultural sector and the other economic sectors. Moreover, within the agricultural sector, water is increasingly being used for commercial plantations as opposed to growing food crops, threatening food security. Therefore, it is very important that this scarce resource is managed in an efficient and sustainable manner, for now and future use. This requires understanding the process of evaporation for accurate determination of water-use from agricultural lands. In the past, direct measurements of evaporation have proven difficult because of the cost and complexity of the available equipments, and level of expertise involved. This justifies a quest for relatively simple, accurate and inexpensive methods of determining evaporation for routine field applications. Estimation of sensible heat flux (H) from high frequency air temperature measurements and then calculating latent energy flux (ƒÜE) and hence evaporation as a residual of the shortened surface energy balance equation, assuming that closure is met, is appealing in this sense. Concurrent net irradiance (Rn) and soil heat flux (G) measurements can be conducted with relative ease for use in the energy balance equation. Alternately, evaporation can also be mathematically modelled, using single- or multi-layer models depending on vegetation cover, from less expensive routine meteorological observations. Therefore, the ultimate objective of this study is to estimate and model H and ƒÜE, and thereby evaporation, accurately over sparsely vegetated agricultural lands at low cost and effort. Temperature-variance (TV) and surface renewal (SR) methods, which use high-frequency (typically 2 to 10 Hz) air temperature measurements, are employed for estimation of H. The TV method is based on the Monin and Obukhov Similarity Theory (MOST) and uses statistical measures of the high frequency air temperature to estimate H, including adjustments for stability. The SR method is based on the principle that an air parcel near the surface is renewed by an air parcel from above and, to determine H, it uses higher order air temperature differences between two consecutive sample measurements lagged by a certain time interval. Single- and double-layer models that are based on energy and resistance combination theory were also used to estimate evaporation and H from sparse vegetation. Single- and double-layer models that were extended to include inputs of radiometric temperature in order to estimate H were also used. The transmission of solar irradiance to the soil beneath in sparse canopies is variable and depends on the vegetation density, cover and apparent position of the sun. A three-dimensional radiation interception model was developed to estimate this transmission of solar irradiance and was used as a sub-module in the double-layer models. Estimations of H from the TV (HTV), SR (HSR) and double-layer models were compared against H obtained from eddy covariance (HEC), and the modelled ƒÜE (single- and double-layer) were compared with that obtained from the shortened energy balance involving HEC. Besides, long-term ƒÜE calculated from the shortened energy balance using HTV and HSR were compared with those calculated using HEC. Unshielded and naturally-ventilated fine-wire chromel-constantan thermocouples (TCs), 75 ƒÝm in diameter, at different heights above the ground over sparse Jatropha curcas trees, mixed grassland community and bare fallow land were used to measure air temperature. A three-dimensional sonic anemometer mounted at a certain height above the ground surface was also used to measure virtual temperature and wind speed at all three sites. All measurements were done differentially at 10-Hz frequency. Additional measurements of Rn, G and soil water content (upper 60 mm) were also made. The Jatropha trees were planted in a 3-m plant and inter-row spacing in a 50 m ¡Ñ 60 m plot with the surrounding plots planted to a mixture of Jatropha trees and Kikuyu grass. Average tree height and leaf area index measurements were taken on monthly and bimonthly basis respectively. An automatic weather station about 10 m away from the edge of the Jatropha plot was also used to obtain solar irradiance, air temperature and relative humidity, wind speed and direction and precipitation data. Soil water content was measured to a depth of 1000 mm from the surface at 200 mm intervals. Soil and foliage surface temperatures were measured using two nadir-looking infrared thermometers with one mounted directly above bare soil and the other above the trees. The three-dimensional solar irradiance interception model was validated using measurements conducted on different trees and planting patterns. Solar irradiance above and below tree canopies was measured using LI-200 pyranometer and tube solarimeters respectively. Leaf area density (LAD) was estimated from LAI, canopy shape and volume measurements. It was also determined by scanning leaves using either destructive sampling or tracing method. The performance of the TV method over sparse vegetation of J. curcas, mixed grassland community and fallow land was evaluated against HEC. Atmospheric stability conditions were identified using (i) sensor height (z) and Obukhov length (L) obtained from EC and (ii) air temperature difference between two thermocouple measurement heights. The HTV estimations, adjusted and not adjusted for skewness (actual and estimated) of air temperature (sk), for unstable conditions only and for all stability conditions were used. An improved agreement in terms of slope, coefficient of determination (r2) and root mean square error (RMSE), almost over all surfaces, was obtained when the temperature difference rather than the z/L means of identifying stability conditions was used. The agreement between the HTV and HEC was improved for estimations adjusted for actual sk than not adjusted for sk. Improved agreement was also noted when HTV was adjusted using estimated sk compared to not adjusting for sk over J. curcas. The TV method could be used to estimate H for surfaces with varying homogeneity with reasonable accuracy. Long-term water-use of a fetch-limited sparse vegetation of J. curcas was determined as a residual of the shortened surface energy balance involving HTV and HSR and compared with those estimated using HEC. Concurrent measurements of Rn and G were also performed. The long-term water-use of J. curcas trees calculated from the shortened surface energy balance involving HTV and HSR agreed very well when compared with those obtained from HEC. The seasonal HTV and HSR also agreed very well when compared with HEC. Changes in structure of the canopy and environmental conditions appeared to influence partitioning of the available energy into H and ƒÜE. The seasonal total evaporation for the EC, TV and SR methods amounted to 626, 640 and 674 mm respectively with a total rainfall of 690 mm. Footprint analysis also revealed that greater than 80% of the measured flux during the day originates from within the surface of interest. The TV and SR methods, therefore, offer a relatively low-cost means for long-term estimation of H, and ƒÜE, hence the total evaporation, using the shortened surface energy balance along with measurements of Rn and G. Evaporation and biomass production estimations from tree crops requires accurate representation of solar irradiance transmission through the canopy. A relatively simple three-dimensional, hourly time-step tree-canopy radiation interception model was developed and validated using measurements conducted on isolated trees, hedgerows and tree canopies arranged in tramline mode. Measurements were obtained using tube solarimeters placed 0.5 m from each other starting from the base of a tree trunk in four directions, along and perpendicular to the row up to mid-way between trees and rows. Model-simulations of hourly radiant transmittance were in good agreement with measurements with an overall r2 of 0.91; Willmott.s index of agreement of 0.96; and general absolute standard deviation of 17.66%. Agreement between model-estimations and measurements, however, was influenced by distance and direction of the node from the tree trunk, sky conditions, symmetry of the canopy, and uniformity of the stand and leaf distribution of the canopy. The model could be useful in planning and management applications for a wide range of tree crops. Penman-Monteith (PM) equation and the Shuttleworth and Wallace (SW) model, representing single- and dual-source models respectively, were used to determine the total evaporation over a sparse vegetation of J. curcas from routine automatic weather station observations, resistance parameters and vegetation indices. The three-dimensional solar irradiance interception model was used as a sub-module in the SW model. The total evaporation from the sparse vegetation was also determined as a residual of the shortened surface energy balance using measurements of Rn, G and HEC. The PM equation failed to reproduce the .measured. daily total evaporation during periods of low LAI, with improved agreement with increased LAI. The SW model, however, produced total evaporation estimates that agreed very well with the .measured. with a slope of 0.96, r2 of 0.91 and RMSE of 0.45 mm for a LAI ranging from 0 (no leaves) to 1.83 m2 m-2. The SW model also estimated soil evaporation and plant transpiration separately, and about 66 % of the cumulative evaporation was attributed to soil evaporation. These findings suggest that the PM equation should be replaced by the SW model for surfaces that assume a range of LAI values during the growing season. The H was estimated using (i) SW model that was further developed to include surface radiometric temperature measurements; (ii) one-layer model, but linked with a two-layer model for estimation of excess resistance, that uses surface radiometric temperature; and (iii) the SW model (unmodified). The agreement between modelled and measured H, using 10-min data, was in general reasonably good with RMSE (W m-2) of 45.11, 43.77 and 39.86 for the three models respectively. The comparative results that were achieved from (iii) were not translated into the daily data as all models appeared to have a tendency to underestimate H. The resulting RMSEs for the daily H data for the three models were (MJ m-2) 1.16, 1.17 and 1.18 respectively. It appears that similar or better agreement between measured and estimated H can be forged without the need for surface radiometric temperature measurements. The study showed, in general, that high frequency air temperature measurements can be used to estimate H with reasonable accuracy using the simple and relatively low-cost TV and SR methods. Moreover, these methods can be used to calculate ƒÜE, hence ET, as a residual of the shortened surface energy balance equation along with measurements of Rn and G assuming that energy balance closure is met. The simple and low-cost nature of these methods makes replication of measurements easier and their robust nature allows long-term measurements of energy fluxes. The study also showed that H and ƒÜE can be modeled using energy and resistance combination equations with reasonable accuracy. It also reiterated that the SW-type models, which treat the plant canopy and soil components separately, are more appropriate for estimation of H and ƒÜE over sparse vegetation as opposed to the PM-type models.
Thesis (Ph.D.)-University of KwaZulu-Natal, Pietermaritzburg, 2010.

This work was written as a cumulative doctoral thesis based on reviewed publications. Climate projections are mainly based on the results of numeric simulations from global or regional climate models. Up to now processes between atmosphere and land surface are only rudimentarily known. This causes one of the major uncertainties in existing models. In order to reduce parameterisation uncertainties and to find a reasonable description of sub grid heterogeneities, the determination and evaluation of parameterisation schemes for modelling require as many datasets from different spatial scales as possible. This work contributes to this topic by implying different datasets from different platforms. Its objective was to analyse the spatial heterogeneity of land surface parameters and energy flux densities obtained from both satellite observations with different spatial and temporal resolutions and in-situ measurements. The investigations were carried out for two target areas in Germany. First, satellite data for the years 2002 and 2003 were analysed and validated from the LITFASS-area (Lindenberg Inhomogeneous Terrain - Fluxes between Atmosphere and Surface: a longterm Study). Second, the data from the experimental field sites of the FLUXNET cluster around Tharandt from the years 2006 and 2007 were used to determine the NDVI (Normalised Difference Vegetation Index for identifying vegetated areas and their "condition"). The core of the study was the determination of land surface characteristics and hence radiant and energy flux densities (net radiation, soil heat flux, sensible and latent heat flux) using the three optical satellite sensors ETM+ (Enhanced Thematic Mapper), MODIS (Moderate Resolution Imaging Spektroradiometer) and AVHRR 3 (Advanced Very High Resolution Radiometer) with different spatial (30 m – 1 km) and temporal (1 day – 16 days) resolution. Different sensor characteristics and different data sets for land use classifications can both lead to deviations of the resultant energy fluxes between the sensors. Thus, sensor differences were quantified, sensor adaptation methods were implemented and a quality analysis for land use classifications was performed. The result is then a single parameterisation scheme that allows for the determination of the energy fluxes from all three different sensors. The main focus was the derivation of the latent heat flux (L.E) using the Penman-Monteith (P-M) approach. Satellite data provide measurements of spectral reflectance and surface temperatures. The P-M approach requires further surface parameters not offered by satellite data. These parameters include the NDVI, Leaf Area Index (LAI), wind speed, relative humidity, vegetation height and roughness length, for example. They were derived indirectly from the given satellite- or in-situ measurements. If no data were available so called default values from literature were taken. The quality of these parameters strongly influenced the exactness of the radiant- and energy fluxes. Sensitivity studies showed that NDVI is one of the most important parameters for determination of evaporation. In contrast it could be shown, that the parameters as vegetation height and measurement height have only minor influence on L.E, which justifies the use of default values for these parameters. Due to the key role of NDVI a field study was carried out investigating the spatial variability and sensitivity of NDVI above five different land use types (winter wheat, corn, grass, beech and spruce). Methods to determine this parameter not only from space (spectral), but also from in-situ tower measurements (broadband) and spectrometer data (spectral) were compared. The best agreement between the methods was found for winter wheat and grass measurements in 2006. For these land use types the results differed by less than 10 % and 15 %, respectively. Larger differences were obtained for the forest measurements. The correlation between the daily MODIS-NDVI data and the in-situ NDVI inferred from the spectrometer and the broadband measurements were r=0.67 and r=0.51, respectively. Subsequently, spatial variability of land surface parameters and fluxes were analysed. The several spatial resolutions of the satellite sensors can be used to describe subscale heterogeneity from one scale to the other and to study the effects of spatial averaging. Therefore land use dependent parameters and fluxes were investigated to find typical distribution patterns of land surface properties and energy fluxes. Implying the distribution patterns found here for albedo and NDVI from ETM+ data in models has high potential to calculate representative energy flux distributions on a coarser scale. The distribution patterns were expressed as probability density functions (PDFs). First results of applying PDFs of albedo, NDVI, relative humidity, and wind speed to the L.E computation are encouraging, and they show the high potential of this method. Summing up, the method of satellite based surface parameter- and energy flux determination has been shown to work reliably on different temporal and spatial scales. The data are useful for detailed analyses of spatial variability of a landscape and for the description of sub grid heterogeneity, as it is needed in model applications. Their usability as input parameters for modelling on different scales is the second important result of this work. The derived vegetation parameters, e.g. LAI and plant cover, possess realistic values and were used as model input for the Lokalmodell of the German Weather Service. This significantly improved the model results for L.E. Additionally, thermal parameter fields, e.g. surface temperature from ETM+ with 30 m spatial resolution, were used as input for SVAT-modelling (Soil-Vegetation-Atmosphere-Transfer scheme). Thus, more realistic L.E results were obtained, providing highly resolved areal information.
Die vorliegende Arbeit wurde auf der Grundlage begutachteter Publikationen als kumulative Dissertation verfasst. Klimaprognosen basieren im Allgemeinen auf den Ergebnissen numerischer Simulationen mit globalen oder regionalen Klimamodellen. Eine der entscheidenden Unsicherheiten bestehender Modelle liegt in dem noch unzureichenden Verständnis von Wechselwirkungsprozessen zwischen der Atmosphäre und Landoberflächen und dem daraus folgenden Fehlen entsprechender Parametrisierungen. Um das Problem einer unsicheren Modell-Parametrisierung aufzugreifen und zum Beispiel subskalige Heterogenität in einer Art und Weise zu beschreiben, dass sie für Modelle nutzbar wird, werden für die Bestimmung und Evaluierung von Modell-Parametrisierungsansätzen so viele Datensätze wie möglich benötigt. Die Arbeit trägt zu diesem Thema durch die Verwendung verschiedener Datensätze unterschiedlicher Plattformen bei. Ziel der Studie war es, aus Satellitendaten verschiedener räumlicher und zeitlicher Auflösung sowie aus in-situ Daten die räumliche Heterogenität von Landoberflächenparametern und Energieflussdichten zu bestimmen. Die Untersuchungen wurden für zwei Zielgebiete in Deutschland durchgeführt. Für das LITFASS-Gebiet (Lindenberg Inhomogeneous Terrain - Fluxes between Atmosphere and Surface: a longterm Study) wurden Satellitendaten der Jahre 2002 und 2003 untersucht und validiert. Zusätzlich wurde im Rahmen dieser Arbeit eine NDVI-Studie (Normalisierter Differenzen Vegetations Index: Maß zur Detektierung von Vegetationflächen, deren Vitalität und Dichte) auf den Testflächen des FLUXNET Clusters um Tharandt in den Jahren 2006 und 2007 realisiert. Die Grundlage der Arbeit bildete die Bestimmung von Landoberflächeneigenschaften und daraus resultierenden Energieflüssen, auf Basis dreier optischer Sensoren (ETM+ (Enhanced Thematic Mapper), MODIS (Moderate Resolution Imaging Spectroradiometer) und AVHRR 3 (Advanced Very High Resolution Radiometer)) mit unterschiedlichen räumlichen (30 m – 1 km) und zeitlichen (1 – 16 Tage) Auflösungen. Unterschiedliche Sensorcharakteristiken, sowie die Verwendung verschiedener, zum Teil ungenauer Datensätze zur Landnutzungsklassifikation führen zu Abweichungen in den Ergebnissen der einzelnen Sensoren. Durch die Quantifizierung der Sensorunterschiede, die Anpassung der Ergebnisse der Sensoren aneinander und eine Qualitätsanalyse von verschiedenen Landnutzungsklassifikationen, wurde eine Basis für eine vergleichbare Parametrisierung der Oberflächenparameter und damit auch für die daraus berechneten Energieflüsse geschaffen. Der Schwerpunkt lag dabei auf der Bestimmung des latenten Wärmestromes (L.E) mit Hilfe des Penman-Monteith Ansatzes (P-M). Satellitendaten liefern Messwerte der spektralen Reflexion und der Oberflächentemperatur. Die P-M Gleichung erfordert weitere Oberflächenparameter wie zum Beispiel den NDVI, den Blattflächenindex (LAI), die Windgeschwindigkeit, die relative Luftfeuchte, die Vegetationshöhe oder die Rauhigkeitslänge, die jedoch aus den Satellitendaten nicht bestimmt werden können. Sie müssen indirekt aus den oben genannten Messgrößen der Satelliten oder aus in-situ Messungen abgeleitet werden. Stehen auch aus diesen Quellen keine Daten zur Verfügung, können sogenannte Standard- (Default-) Werte aus der Literatur verwendet werden. Die Qualität dieser Parameter hat einen großen Einfluss auf die Bestimmung der Strahlungs- und Energieflüsse. Sensitivitätsstudien im Rahmen der Arbeit zeigen die Bedeutung des NDVI als einen der wichtigsten Parameter in der Verdunstungsbestimmung nach P-M. Im Gegensatz dazu wurde deutlich, dass z. B. die Vegetationshöhe und die Messhöhe einen relativ kleinen Einfluss auf L.E haben, so dass für diese Parameter die Verwendung von Standardwerten gerechtfertigt ist. Aufgrund der Schlüsselrolle, welche der NDVI in der Bestimmung der Verdunstung einnimmt, wurden im Rahmen einer Feldstudie Untersuchungen des NDVI über fünf verschiedenen Landnutzungstypen (Winterweizen, Mais, Gras, Buche und Fichte) hinsichtlich seiner räumlichen Variabilität und Sensitivität, unternommen. Dabei wurden verschiedene Bestimmungsmethoden getestet, in welchen der NDVI nicht nur aus Satellitendaten (spektral), sondern auch aus in-situ Turmmessungen (breitbandig) und Spekrometermessungen (spektral) ermittelt wird. Die besten Übereinstimmungen der Ergebnisse wurden dabei für Winterweizen und Gras für das Jahr 2006 gefunden. Für diese Landnutzungstypen betrugen die Maximaldifferenzen aus den drei Methoden jeweils 10 beziehungsweise 15 %. Deutlichere Differenzen ließen sich für die Forstflächen verzeichnen. Die Korrelation zwischen Satelliten- und Spektrometermessung betrug r=0.67. Für Satelliten- und Turmmessungen ergab sich ein Wert von r=0.5. Basierend auf den beschriebenen Vorarbeiten wurde die räumliche Variabilität von Landoberflächenparametern und Flüssen untersucht. Die unterschiedlichen räumlichen Auflösungen der Satelliten können genutzt werden, um zum einen die subskalige Heterogenität zu beschreiben, aber auch, um den Effekt räumlicher Mittelungsverfahren zu testen. Dafür wurden Parameter und Energieflüsse in Abhängigkeit der Landnutzungsklasse untersucht, um typische Verteilungsmuster dieser Größen zu finden. Die Verwendung der Verteilungsmuster (in Form von Wahrscheinlichkeitsdichteverteilungen – PDFs), die für die Albedo und den NDVI aus ETM+ Daten gefunden wurden, bietet ein hohes Potential als Modellinput, um repräsentative PDFs der Energieflüsse auf gröberen Skalen zu erhalten. Die ersten Ergebnisse in der Verwendung der PDFs von Albedo, NDVI, relativer Luftfeuchtigkeit und Windgeschwindigkeit für die Bestimmung von L.E waren sehr ermutigend und zeigten das hohe Potential der Methode. Zusammenfassend lässt sich feststellen, dass die Methode der Ableitung von Oberflächenparametern und Energieflüssen aus Satellitendaten zuverlässige Daten auf verschiedenen zeitlichen und räumlichen Skalen liefert. Die Daten sind für eine detaillierte Analyse der räumlichen Variabilität der Landschaft und für die Beschreibung der subskaligen Heterogenität, wie sie oft in Modellanwendungen benötigt wird, geeignet. Ihre Nutzbarkeit als Inputparameter in Modellen auf verschiedenen Skalen ist das zweite wichtige Ergebnis der Arbeit. Aus Satellitendaten abgeleitete Vegetationsparameter wie der LAI oder die Pflanzenbedeckung liefern realistische Ergebnisse, die zum Beispiel als Modellinput in das Lokalmodell des Deutschen Wetterdienstes implementiert werden konnten und die Modellergebnisse von L.E signifikant verbessert haben. Aber auch thermale Parameter, wie beispielsweise die Oberflächentemperatur aus ETM+ Daten in 30 m Auflösung, wurden als Eingabeparameter eines Soil-Vegetation-Atmosphere-Transfer-Modells (SVAT) verwendet. Dadurch erhält man realistischere Ergebnisse für L.E, die hochaufgelöste Flächeninformationen bieten.

Terrestrial evapotranspiration (ET) is a significant component of the energy and water balances at the land surface. However, direct, continuous measurements of ET are spatially limited and only available since the 1990s. Due to this lack of observations, detecting and attributing long-term regional trends in ET remains difficult. This dissertation aims to alleviate the data limitation and detect long-term trends by developing a method to infer ET from data collected at common weather stations, which are spatially and temporally abundant. The methodology used to infer ET from historical meteorological data is based on an emergent relation between the land surface and atmospheric boundary layer. We refer to this methodology as the Evapotranspiration from Relative Humidity at Equilibrium method, or the “ETRHEQ method”. In the first section of this dissertation, we develop the ETRHEQ method for use at common weather stations and demonstrate the utility of the method at twenty eddy covariance sites spanning a wide range of climate and plant functional types. Next, we apply the ETRHEQ method at historical weather stations across the continental U.S. and show that ET estimates obtained via the ETRHEQ method compare well with watershed scale ET, as well as ET estimates from land surface models. From 1961 to 1997, we find negligible or increasing trends in summertime ET over the central U.S. and the west coast and negative trends in the eastern and western U.S. From 1998 to 2014, we find a sharp decline in summertime ET across the entire U.S. We show that this decline is consistent with decreasing transpiration associated with declines in humidity. Lastly, we assess the sensitivity of ET to perturbations in soil moisture and humidity anticipated with climate change. We demonstrate that the response of ET to changing humidity and soil moisture is strongly dependent on the biological and hydrological state of the surface, particularly the degree of water stress and vegetation fraction. In total, this dissertation demonstrates the utility of the ETRHEQ method as a means to estimate ET from weather station data and highlights the critical role of vegetation in modulating ET variability.

碩士
國立交通大學
機械工程學系
98
An experiment is carried out in the present study to investigate the characteristics of time periodic evaporation heat transfer for refrigerant R-410A flowing in a horizontal narrow annular duct subject to an imposed time periodic mass flux oscillation or heat flux oscillation. The mass flux oscillation and heat flux oscillation are both in the form of triangular waves. The gap of the duct is fixed at 2.0 mm. In the study, the effects of the refrigerant mass flux oscillation, heat flux oscillation, R-410A saturation temperature, imposed heat flux and vapor quality of the refrigerant on the temporal evaporating flow heat transfer and the photos of the evaporating flow will be examined in detail. The present experiment is conducted for the mean refrigerant mass flux varied from 300 to 500 kg/m2s, the amplitude of the mass flux oscillation is fixed at 10, 20 and 30% with the period of the mass flux oscillation tp fixed at 20, 60 and 120 seconds. Besides, the mean imposed heat flux is varied from 0 KW/ m2 to 15 KW/ m2, the amplitude of the heat flux oscillation is chosen to vary from 10% to 50% of the mean heat flux , and the period of the q oscillation tp is also fixed at 20, 60, 120 seconds. The mean refrigerant saturation temperature is set at 5, 10 and 15 ℃ for the mean refrigerant vapor quality varied from 0.05 to 0.95. The measured evaporation heat transfer data are expressed in terms of the variations of the heated wall temperature and evaporation heat transfer coefficient with time. In the first part of the study the results for the R-410A evaporation subject to the mass flux oscillation are presented. The measured heat transfer data for the R-410A evaporating flow for a constant coolant mass flux are first compared with the time-average data for a time periodic mass flux oscillation. This comparison shows that the mass flux oscillation exerts negligible influences on the time-average evaporation heat transfer. Then, we present the data to elucidate the effects of the experimental parameters on the amplitude of Tw oscillation over wide ranges of the experimental parameters. The results indicate that the Tw oscillation is stronger for higher amplitude and a longer period of the mass flux oscillation. However, a small time lag in the Tw oscillation is also noted. Moreover, at the intermediate vapor quality changes in the evaporating flow patterns between that dominated by the nucleation bubbles and liquid film take place cyclically. Furthermore, after the time lag the heated pipe wall temperature decreases and the evaporation heat transfer gets better as the mass flux decreases in the first half of the periodic cycle. In the second half of the cycle in which the mass flux increases the opposite processes occur. These unusual changes of the heating surface temperature and heat transfer coefficient with the mass flux oscillation are attributed to the strong effects of the mass flux oscillation on the state of the refrigerant at the duct inlet and hence on the changes of the vapor quality and liquid film thickness in the evaporating flow. In the second part of the study results for the R-410A evaporation subject to the heat flux oscillation are also presented. Effects of the mean level and oscillation amplitude and period of the heat flux on the time periodic R-410A evaporation heat transfer have been investigated in detail. We first note that the time-average heat transfer coefficients for the time periodic evaporation of R-410A are not affected to a noticeable degree by the amplitude and period of the imposed heat flux oscillation. Then, the heated pipe wall temperature and evaporating flow pattern also oscillate periodically in time and at the same frequency as the heat flux oscillation. Experiment also shows that the resulting oscillation amplitudes of the wall temperature get longer for a longer period and a larger amplitude of the imposed heat flux oscillation and for a higher mean imposed heat flux. A significant time lag in the heated surface temperature oscillation is also noted, which apparently results from the thermal inertia of the copper inner pipe. Moreover, at the intermediate vapor quality changes in the evaporating flow pattern between that dominated by the nucleation bubbles and liquid film take place cyclically. Furthermore, after the time lag the heated pipe wall temperature decreases and the evaporation heat transfer gets worse as the heat flux decreases in the first half of the periodic cycle. In the second half of the cycle in which the heat flux increases the opposite processes occur. These changes of the heating surface temperature and heat transfer coefficient with the heat flux oscillation are attributed to the strong effects of the heat flux oscillation and hence on the changes of the vapor quality and liquid film thickness in the evaporating flow. The effects of heat flux oscillation at extremely short and long periods have been explored. Due to the existence of the thermal inertia of the heated copper duct, the resulting heated surface temperature does not oscillate with time at an extremely short period of the imposed heat flux oscillation. But the oscillation amplitude of the heated surface temperature gets noticeably stronger for an extremely long period of the imposed heat flux oscillation.

As warships can carry weapons on board, the unlikely occurrence of fire is one of the most feared events on board. Until recently, halon 1301 (bromo-tri-fluoro methane, CF3Br) has been the primary fire-fighting agent for protecting the machinery spaces of ships. Halon 1301 is not only harmful to humans, but it also depletes the ozone layer. Water-mist fire suppression systems (WMFSS) have been considered as a potential candidate for the replacement of halon-based fire suppression systems by fire protection industries. WMFSS is already being used in commercial buildings, passenger and naval ships, etc. However, it is essential to examine the efficacy of water-mist droplets in suppressing fires. The efficacy of a water-mist system can be investigated in two ways: (i) experimental investigation; and (ii) numerical analysis. This study is a combination of an experimental study (water mist spray without fire) and two types of numerical studies using (a) semi-empirical equations based model developed in this study; and (b) a state of the art computational fluids dynamics (CFD) based fire model.