Implementation of a phenomenological evaporation model into a porous network simulation for water management in low temperature fuel cells

A phenomenological transition film evaporation model was introduced to a pore network model with the consideration of pore radius, contact angle, non-isothermal interface temperature, microscale fluid flows and heat and mass transfers. This was achieved by modeling the transition film region of the menisci in each pore throughout the porous transport layer of a half-cell polymer electrolyte membrane (PEM) fuel cell. The model presented in this research is compared with the standard diffusive fuel cell modeling approach to evaporation and shown to surpass the conventional modeling approach in terms of predicting the evaporation rates in porous media. The current diffusive evaporation models used in many fuel cell transport models assumes a constant evaporation rate across the entire liquid-air interface. The transition film model was implemented into the pore network model to address this issue and create a pore size dependency on the evaporation rates. This is accomplished by evaluating the transition film evaporation rates determined by the kinetic model for every pore containing liquid water in the porous transport layer (PTL). The comparison of a transition film and diffusive evaporation model shows an increase in predicted evaporation rates for smaller pore sizes with the transition film model. This is an important parameter when considering the micro-scaled pore sizes seen in the PTL and becomes even more substantial when considering transport in fuel cells containing an MPL, or a large variance in pore size. Experimentation was performed to validate the transition film model by monitoring evaporation rates from a non-zero contact angle water droplet on a heated substrate. The substrate was a glass plate with a hydrophobic coating to reduce wettability. The tests were performed at a constant substrate temperature and relative humidity. The transition film model was able to accurately predict the drop volume as time elapsed. By implementing the transition film model to a pore network model the evaporation rates present in the PTL can be more accurately modeled. This improves the ability of a pore network model to predict the distribution of liquid water and ultimately the level of flooding exhibited in a PTL for various operating conditions.

Development of a direct evaporation bismuth Hall thruster

Hall thrusters have been under active development around the world since the 1960’s. Thrusters using traditional propellants such as xenon have been flown on a variety of satellite orbit raising and maintenance missions with an excellent record. To expand the mission envelope, it is necessary to lower the specific impulse of the thrusters but xenon and krypton are poor performers at specific impulses below 1,200 seconds. To enhance low specific impulse performance, this dissertation examines the development of a Hall-effect thruster which uses bismuth as a propellant. Bismuth, the heaviest non-radioactive element, holds many advantages over noble gas propellants from an energetics as well as a practical economic standpoint. Low ionization energy, large electron-impact crosssection and high atomic mass make bismuth ideal for low-specific impulse applications. The primary disadvantage lies in the high temperatures which are required to generate the bismuth vapors. Previous efforts carried out in the Soviet Union relied upon the complete bismuth vaporization and gas phase delivery to the anode. While this proved successful, the power required to vaporize and maintain gas phase throughout the mass flow system quickly removed many of the efficiency gains expected from using bismuth. To solve these problems, a unique method of delivering liquid bismuth to the anode has been developed. Bismuth is contained within a hollow anode reservoir that is capped by a porous metallic disc. By utilizing the inherent waste heat generated in a Hall thruster, liquid bismuth is evaporated and the vapors pass through the porous disc into the discharge chamber. Due to the high temperatures and material compatibility requirements, the anode was fabricated out of pure molybdenum. The porous vaporizer was not available commercially so a method of creating a refractory porous plate with 40-50% open porosity was developed. Molybdenum also does not respond well to most forms of welding so a diffusion bonding process was also developed to join the molybdenum porous disc to the molybdenum anode. Operation of the direct evaporation bismuth Hall thruster revealed interesting phenomenon. By utilizing constant current mode on a discharge power supply, the discharge voltage settles out to a stable operating point which is a function of discharge current, anode face area and average pore size on the vaporizer. Oscillations with a 40 second period were also observed. Preliminary performance data suggests that the direct evaporation bismuth Hall thruster performs similar to xenon and krypton Hall thrusters. Plume interrogation with a Retarding Potential Analyzer confirmed that bismuth ions were being efficiently accelerated while Faraday probe data gave a view of the ion density in the exhausted plume.

Investigations of Cloud Microphysical Response to Mixing Using Digital Holography

Cloud edge mixing plays an important role in the life cycle and development of clouds. Entrainment of subsaturated air affects the cloud at the microscale, altering the number density and size distribution of its droplets. The resulting effect is determined by two timescales: the time required for the mixing event to complete, and the time required for the droplets to adjust to their new environment. If mixing is rapid, evaporation of droplets is uniform and said to be homogeneous in nature. In contrast, slow mixing (compared to the adjustment timescale) results in the droplets adjusting to the transient state of the mixture, producing an inhomogeneous result. Studying this process in real clouds involves the use of airborne optical instruments capable of measuring clouds at the `single particle' level. Single particle resolution allows for direct measurement of the droplet size distribution. This is in contrast to other `bulk' methods (i.e. hot-wire probes, lidar, radar) which measure a higher order moment of the distribution and require assumptions about the distribution shape to compute a size distribution. The sampling strategy of current optical instruments requires them to integrate over a path tens to hundreds of meters to form a single size distribution. This is much larger than typical mixing scales (which can extend down to the order of centimeters), resulting in difficulties resolving mixing signatures. The Holodec is an optical particle instrument that uses digital holography to record discrete, local volumes of droplets. This method allows for statistically significant size distributions to be calculated for centimeter scale volumes, allowing for full resolution at the scales important to the mixing process. The hologram also records the three dimensional position of all particles within the volume, allowing for the spatial structure of the cloud volume to be studied. Both of these features represent a new and unique view into the mixing problem. In this dissertation, holographic data recorded during two different field projects is analyzed to study the mixing structure of cumulus clouds. Using Holodec data, it is shown that mixing at cloud top can produce regions of clear but humid air that can subside down along the edge of the cloud as a narrow shell, or advect down shear as a `humid halo'. This air is then entrained into the cloud at lower levels, producing mixing that appears to be very inhomogeneous. This inhomogeneous-like mixing is shown to be well correlated with regions containing elevated concentrations of large droplets. This is used to argue in favor of the hypothesis that dilution can lead to enhanced droplet growth rates. I also make observations on the microscale spatial structure of observed cloud volumes recorded by the Holodec.