SIGNIFICANCE OF HEAT DISSIPATION IN THE CATHODE CATALYST LAYER ON THE LIFETIME AND RELIABILITY OF POLYMER ELECTROLYTE FUEL CELLS

To study the dissipation of heat generated due to the formation of pinholes that cause local hotspots in the catalyst layer of the Polymer Electrolyte Fuel Cell, a two-phase non-isothermal model has been developed by coupling Darcy’s law with heat transport. The domain under consideration is a section of the membrane electrode assembly with a half-channel and a half-rib. Five potential locations where a pinhole might form were analyzed: at the midplane of the channel, midway between the channel midplane and the channel wall, at the channel or rib wall, midway between the rib midplane and the channel wall, at the midplane of the rib. In the first part of this work, a preliminary thermal model was developed. The model was then refined to account for the two-phase effects. A sensitivity study was done to evaluate the effect of the following properties on the maximum temperature in the domain: Catalyst layer thermal conductivity, the Microporous layer thermal conductivity, the anisotropy factor of the Catalyst layer thermal conductivity, the Porous transport layer porosity, the liquid water distribution and the thickness of the membrane and porous layers. Accounting for the two-phase effects, a slight cooling effect was observed across all hotspot locations. The thermal properties of the catalyst layer were shown to have a limited impact on the maximum temperature in the catalyst layer of new fuel cells without pinhole. However, as hotspots start to appear, thermal properties play a more significant role in mitigating the thermal runaway.

Heat Dissipation of Hybrid Iron Oxide-Gold Nanoparticles in an Agar Phantom

Hybrid iron oxide-gold nanoparticles (HNPs) have shown potential in cancer therapy as agents for tumour ablation
and thermal switches for targeted drug release. Heat generation occurs by exploitation of the surface plasmon
resonance of the gold coating, which usually occurs at the maximum UV absorption wavelength. However, lasers
at such wavelength are often expensive and highly specialised. Here, we report the heating and monitoring of heat
dissipation of HNPs suspended in agar phantoms using a relatively inexpensive Ng: YAG pulsed 1064 nm laser source.
The particles experience heating of up to 40°C with a total area of heat dissipation up to 132.73 mm2 from the 1 mm
diameter irradiation point after 60 seconds. This work reports the potential and possible drawbacks of these particles
for translation into cancer therapy based on our findings.

The storage and dissipation of magnetic energy in the quiet sun corona determined from SDO/HMI magnetograms

In recent years, higher cadence, higher resolution observations have revealed the quiet-Sun photosphere to be complex and rapidly evolving. Since magnetic fields anchored in the photosphere extend up into the solar corona, it is expected that the small-scale coronal magnetic field exhibits similar complexity. For the first time, the quiet-Sun coronal magnetic field is continuously evolved through a series of non-potential, quasi-static equilibria, deduced from magnetograms observed by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, where the photospheric boundary condition which drives the coronal evolution exactly reproduces the observed magnetograms. The build-up, storage, and dissipation of magnetic energy within the simulations is studied. We find that the free magnetic energy built up and stored within the field is sufficient to explain small-scale, impulsive events such as nanoflares. On comparing with coronal images of the same region, the energy storage and dissipation visually reproduces many of the observed features. The results indicate that the complex small-scale magnetic evolution of a large number of magnetic features is a key element in explaining the nature of the solar corona.

Heat dissipation from carbon nano-electronics

The incorporation of graphitic compounds such as carbon nanotubes (CNTs) and graphene into nano-electronic device packaging holds much promise for waste heat management given their high thermal conductivities. However, as these graphitic materials must be used in together with other semiconductor/insulator materials, it is not known how thermal transport is affected by the interaction. Using different simulation techniques, in this thesis, we evaluate the thermal transport properties - thermal boundary conductance (TBC) and thermal conductivity - of CNTs and single-layer graphene in contact with an amorphous SiO2 (a-SiO2) substrate. First, the theoretical methodologies and concepts used in our simulations are presented. In particular, two concepts are described in detail as they are necessary for the understanding of the subsequent chapters. The first is the linear response Green-Kubo (GK) theory of thermal boundary conductance (TBC), which we develop in this thesis, and the second is the spectral energy density method, which we use to directly compute the phonon lifetimes and thermal transport coefficients. After we set the conceptual foundations, the TBC of the CNT-SiO2 interface is computed using non- equilibrium molecular dynamics (MD) simulations and the new Green-Kubo method that we have developed. Its dependence on temperature, the strength of the interaction with the substrate, and tube diameter are evaluated. To gain further insight into the phonon dynamics in supported CNTs, the scattering rates are computed using the spectral energy density (SED) method. With this method, we are able to distinguish the different scattering mechanisms (boundary and CNT-substrate phonon-phonon) and rates. The phonon lifetimes in supported CNTs are found to be reduced by contact with the substrate and we use that lifetime reduction to determine the change in CNT thermal conductivity. Next, we examine thermal transport in graphene supported on SiO2. The phonon contribution to the TBC of the graphene-SiO2 interface is computed from MD simulations and found to agree well with experimentally measured values. We derive the theory of remote phonon scattering of graphene electrons and compute the heat transfer coefficient dependence on doping level and temperature. The thermal boundary conductance from remote phonon scattering is found to be an order of magnitude smaller than that of the phonon contribution. The in-plane thermal conductivity of supported graphene is calculated from MD simulations. The experimentally measured order of magnitude reduction in thermal conductivity is reproduced in our simulations. We show that this reduction is due to the damping of the flexural (ZA) modes. By varying the interaction between graphene and the substrate, the ZA modes hybridize with the substrate Rayleigh modes and the dispersion of the hybridized modes is found to linearize in the strong coupling limit, leading to an increased thermal conductance in the composite structure.

Tracking the attenuation and nonbreaking dissipation of swells using altimeters

A method for systematically tracking swells across oceanic basins is developed by taking advantage of high-quality data from space-borne altimeters and wave model output. The evolution of swells is observed over large distances based on 202 swell events with periods ranging from 12 to 18 s. An empirical attenuation rate of swell energy of about 4 × 10−7 m−1 is estimated using these observations, and the nonbreaking energy dissipation rates of swells far away from their generating areas are also estimated using a point source model. The resulting acceptance range of nonbreaking dissipation rates is −2.5 to 5.0 × 10−7 m−1, which corresponds to a dissipation e-folding scales of at least 2000 km for steep swells, to almost infinite for small-amplitude swells. These resulting rates are consistent with previous studies using in-situ and synthetic aperture radar (SAR) observations. The frequency dispersion and angular spreading effects during swell propagation are discussed by comparing the results with other studies, demonstrating that they are the two dominant processes for swell height attenuation, especially in the near field. The resulting dissipation rates from these observations can be used as a reference for ocean engineering and wave modeling, and for related studies such as air-sea and wind-wave-turbulence interactions.

Swell dissipation from 10 years of Envisat ASAR in wave mode

Swells are found in all oceans and strongly influence the wave climate and air-sea processes. The poorly known swell dissipation is the largest source of error in wave forecasts and hindcasts. We use synthetic aperture radar data to identify swell sources and trajectories, allowing a statistically significant estimation of swell dissipation. We mined the entire Envisat mission 2003–2012 to find suitable storms with swells (13 < T < 18 s) that are observed several times along their propagation. This database of swell events provides a comprehensive view of swell extending previous efforts. The analysis reveals that swell dissipation weakly correlates with the wave steepness, wind speed, orbital wave velocity, and the relative direction of wind and waves. Although several negative dissipation rates are found, there are uncertainties in the synthetic aperture radar-derived swell heights and dissipation rates. An acceptable range of the swell dissipation rate is −0.1 to 6 × 10−7 m−1 with a median of 1 × 10−7 m−1.

Ocean waves across the Arctic: attenuation due to dissipation dominates over scattering for periods longer than 19 s

The poorly understood attenuation of surface waves in sea ice is generally attributed to the combination of scattering and dissipation. Scattering and dissipation have very different effects on the directional and temporal distribution of wave energy, making it possible to better understand their relative importance by analysis of swell directional spreading and arrival times. Here we compare results of a spectral wave model – using adjustable scattering and dissipation attenuation formulations – with wave measurements far inside the ice pack. In this case, scattering plays a negligible role in the attenuation of long swells. Specifically, scattering-dominated attenuation would produce directional wave spectra much broader than the ones recorded, and swell events arriving later and lasting much longer than observed. Details of the dissipation process remain uncertain. Average dissipation rates are consistent with creep effects but are 12 times those expected for a laminar boundary layer under a smooth solid ice plate.

Inertial dissipation method applied to derive turbulent fluxed over the ocean during the Surface of the Ocean, Fluxes and Interactions with the Atmosphere Atlantic Stratocumulus Transition Experiment (SOFIA/ASTEX) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments with low to moderate wind speed

The transfer coefficients for momentum and heat have been determined for 10 m neutral wind speeds (U-10n) between 0 and 12 m/s using data from the Surface of the Ocean, Fluxes and Interactions with the Atmosphere (SOFIA) and Structure des Echanges Mer-Atmosphere, Proprietes des Heterogeneites Oceaniques: Recherche Experimentale (SEMAPHORE) experiments. The inertial dissipation method was applied to wind and pseudo virtual temperature spectra from a sonic anemometer, mounted on a platform (ship) which was moving through the turbulence held. Under unstable conditions the assumptions concerning the turbulent kinetic energy (TKE) budget appeared incorrect. Using a bulk estimate for the stability parameter, Z/L (where Z is the height and L is the Obukhov length), this resulted in anomalously low drag coefficients compared to neutral conditions. Determining Z/L iteratively, a low rate of convergence was achieved. It was concluded that the divergence of the turbulent transport of TKE was not negligible under unstable conditions. By minimizing the dependence of the calculated neutral drag coefficient on stability, this term was estimated at about -0.65Z/L. The resulting turbulent fluxes were then in close agreement with other studies at moderate wind speed. The drag and exchange coefficients for low wind speeds were found to be C-en x 10(3) = 2.79U(10n)(-1) + 0.66 (U-10n < 5.2 m/s), C-en x 10(3) = C-hn x 10(3) = 1.2 (U-10n greater than or equal to 5.2 m/s), and C-dn x 10(3) = 11.7U(10n)(-2) + 0.668 (U-10n < 5.5 m/s), which imply a rapid increase of the coefficient values as the wind decreased within the smooth flow regime. The frozen turbulence hypothesis and the assumptions of isotropy and an inertial subrange were found to remain valid at these low wind speeds for these shipboard measurements. Incorporation of a free convection parameterization had little effect.

Dissipation of the energy imparted by mid-latitude storms in the Southern Ocean

The aim of this study is to clarify the role of the Southern Ocean storms on interior mixing and meridional overturning circulation. A periodic and idealized numerical model has been designed to represent the key physical processes of a zonal portion of the Southern Ocean located between 70 and 40° S. It incorporates physical ingredients deemed essential for Southern Ocean functioning: rough topography, seasonally varying air–sea fluxes, and high-latitude storms with analytical form. The forcing strategy ensures that the time mean wind stress is the same between the different simulations, so the effect of the storms on the mean wind stress and resulting impacts on the Southern Ocean dynamics are not considered in this study. Level and distribution of mixing attributable to high-frequency winds are quantified and compared to those generated by eddy–topography interactions and dissipation of the balanced flow. Results suggest that (1) the synoptic atmospheric variability alone can generate the levels of mid-depth dissipation frequently observed in the Southern Ocean (10−10–10−9 W kg−1) and (2) the storms strengthen the overturning, primarily through enhanced mixing in the upper 300 m, whereas deeper mixing has a minor effect. The sensitivity of the results to horizontal resolution (20, 5, 2 and 1 km), vertical resolution and numerical choices is evaluated. Challenging issues concerning how numerical models are able to represent interior mixing forced by high-frequency winds are exposed and discussed, particularly in the context of the overturning circulation. Overall, submesoscale-permitting ocean modeling exhibits important delicacies owing to a lack of convergence of key components of its energetics even when reaching Δx =  1 km.

Dissipation of triadimefon on the solid/gas interface

The dissipation or triadimefon, as pure solid and in the Bayleton 5 commercial formulation, was studied under controlled and natural conditions. Volatilization and photodegradation were shown to be the main dissipation processes. The volatilization results can be described by an empirical model assuming exponential decay of the volatilization rate. The filler of the commercial formulation is determinant for the volatilization but has little effect on the photodegradation rates. The main photoproducts were identified and a reaction mechanism proposed. (C) 2001 Elsevier Science Ltd. All rights reserved.

Comparative study of the dissipation of triadimefon in greenhouse and field conditions

The dissipation of triadimefon, {1-(4-chlorophenoxy)-3,3-dimethyl-1-(1H-1,2,4-triazol-1-yl)butanone}, was studied after its application to melon leaves, glass and paper, both in greenhouse and field conditions. The dissipation rate of triadimefon in its commercial formulation Bayleton 5 was found to be lower in greenhouse than field. The results for different samples in the same conditions show that the dissipation of triadimefon was found to be biphasic. This result can be accounted by a semi-empirical model which assumes an initial fast decline of the dissipation rate, attributed to an exponential decay of the volatilization rates, followed by a second phase where the dissipation is due to a first order degradation processes.The dissipation of triadimefon, {1-(4-chlorophenoxy)-3,3-dimethyl-1-(1H- 1,2,4-triazol-1-yl)butan-one}, was studied after its application to melon leaves, glass and paper, both in greenhouse and field conditions. The dissipation rate of triadimefon in its commercial formulation Bayleton 5 was found to be lower in greenhouse than field. The results for different samples in the same conditions show that the dissipation of triadimefon was found to be biphasic. This result can be accounted by a semi-empirical model which assumes an initial fast decline of the dissipation rate, attributed to an exponential decay of the volatilization rates, followed by a second phase where the dissipation is due to a first order degradation processes.

Clustering versus evenly distributing energy dissipation in wireless sensor routing for prolong network lifetime

A novel Cluster Heads (CH) choosing algorithm based on both Minimal Spanning Tree and Maximum Energy resource on sensors, named MSTME, is provided for prolonging lifetime of wireless sensor networks. MSTME can satisfy three principles of optimal CHs: to have the most energy resource among sensors in local clusters, to group approximately the same number of closer sensors into clusters, and to distribute evenly in the networks in terms of location. Simulation shows the network lifetime in MSTME excels its counterparts in two-hop and multi-hop wireless sensor networks.

Energy dissipation and geometry effects over stepped spillways

The energy dissipation process is the major significant point in the designof hydraulic structure. The dissipation of high energy on stepped weirsprevents any damage in the weir itself and channels the energy ownstream to reduce the stilling basin size. In this study, four physical models are used to evaluate the impact of adding end sills that have a quarter circle shape at step edges. The amount of energy loss on weirs under different flow regimes is investigated by experimental work. Stepped weirs have a suitable number of steps and two different ratios of the width to height (2.22, and 2.40). The scale of the physical models is 20:1. The outcomes of the dimensional analyses refer to the critical depth for flow in weirs to the height of step yc/h, the end sill radius, and the number of steps N are more effective parameters than others inthe energy loss process. Moreover, for small values of yc/h, the energydissipation is the greatest. Any increase in yc/h leads to a decrease in theenergy dissipation, while the energy dissipation increases with the number of steps (N).

Initial CFD investigation of an MR fluid in a bicycle ergometer

A bicycle ergometer is a scientific device used by exercise physiologists which attempts to mimic on-road cycling characteristics such as foot technique, EMG activity, VO2, VCO2 and rider cardiology in a laboratory environment. Presently there are no known useful scientific ergometers that mimic these characteristics and are able to provide a satisfactory controlled resistance that is independent of speed. Previous research has suggested the use of a Magneto-Rheological (MR) Fluid as part of the ergometer design, as when used in a rotary brake application it is able to be controlled electronically to increase resistance instantly and independent of speed. In the target application, MR fluids are subject to immense tribological wear and temperature during viscous shearing, and will eventually show some degree of deterioration which is usually manifested as an increase in off-state viscosity. It is not known exactly how the fluid fails, however the amount of deterioration is related to the shear rate, temperature and duration and directly related to the power dissipation. Currently, there is very little literature that investigates the flow and thermal characteristics of MR fluid tribology using CFD. In this paper, we present initial work that aims to improve understanding of MR fluid wear via CFD modelling using Fluent, and results from the model are compared with those obtained from a experimental test rig of an MR fluid-based bicycle ergometer.

Theoretical and numerical investigation of plasmon nanofocusing in metallic tapered rods and grooves

Effective focusing of electromagnetic (EM) energy to nanoscale regions is one of the major challenges in nano-photonics and plasmonics. The strong localization of the optical energy into regions much smaller than allowed by the diffraction limit, also called nanofocusing, offers promising applications in nano-sensor technology, nanofabrication, near-field optics or spectroscopy. One of the most promising solutions to the problem of efficient nanofocusing is related to surface plasmon propagation in metallic structures. Metallic tapered rods, commonly used as probes in near field microscopy and spectroscopy, are of a particular interest. They can provide very strong EM field enhancement at the tip due to surface plasmons (SP’s) propagating towards the tip of the tapered metal rod. A large number of studies have been devoted to the manufacturing process of tapered rods or tapered fibers coated by a metal film. On the other hand, structures such as metallic V-grooves or metal wedges can also provide strong electric field enhancements but manufacturing of these structures is still a challenge. It has been shown, however, that the attainable electric field enhancement at the apex in the V-groove is higher than at the tip of a metal tapered rod when the dissipation level in the metal is strong. Metallic V-grooves also have very promising characteristics as plasmonic waveguides. This thesis will present a thorough theoretical and numerical investigation of nanofocusing during plasmon propagation along a metal tapered rod and into a metallic V-groove. Optimal structural parameters including optimal taper angle, taper length and shape of the taper are determined in order to achieve maximum field enhancement factors at the tip of the nanofocusing structure. An analytical investigation of plasmon nanofocusing by metal tapered rods is carried out by means of the geometric optics approximation (GOA), which is also called adiabatic nanofocusing. However, GOA is applicable only for analysing tapered structures with small taper angles and without considering a terminating tip structure in order to neglect reflections. Rigorous numerical methods are employed for analysing non-adiabatic nanofocusing, by tapered rod and V-grooves with larger taper angles and with a rounded tip. These structures cannot be studied by analytical methods due to the presence of reflected waves from the taper section, the tip and also from (artificial) computational boundaries. A new method is introduced to combine the advantages of GOA and rigorous numerical methods in order to reduce significantly the use of computational resources and yet achieve accurate results for the analysis of large tapered structures, within reasonable calculation time. Detailed comparison between GOA and rigorous numerical methods will be carried out in order to find the critical taper angle of the tapered structures at which GOA is still applicable. It will be demonstrated that optimal taper angles, at which maximum field enhancements occur, coincide with the critical angles, at which GOA is still applicable. It will be shown that the applicability of GOA can be substantially expanded to include structures which could be analysed previously by numerical methods only. The influence of the rounded tip, the taper angle and the role of dissipation onto the plasmon field distribution along the tapered rod and near the tip will be analysed analytically and numerically in detail. It will be demonstrated that electric field enhancement factors of up to ~ 2500 within nanoscale regions are predicted. These are sufficient, for instance, to detect single molecules using surface enhanced Raman spectroscopy (SERS) with the tip of a tapered rod, an approach also known as tip enhanced Raman spectroscopy or TERS. The results obtained in this project will be important for applications for which strong local field enhancement factors are crucial for the performance of devices such as near field microscopes or spectroscopy. The optimal design of nanofocusing structures, at which the delivery of electromagnetic energy to the nanometer region is most efficient, will lead to new applications in near field sensors, near field measuring technology, or generation of nanometer sized energy sources. This includes: applications in tip enhanced Raman spectroscopy (TERS); manipulation of nanoparticles and molecules; efficient coupling of optical energy into and out of plasmonic circuits; second harmonic generation in non-linear optics; or delivery of energy to quantum dots, for instance, for quantum computations.