Design of Latent Thermal Energy Storage Systems

This dissertation documents the results of a theoretical and numerical study of time dependent storage of energy by melting a phase change material. The heating is provided along invading lines, which change from single-line invasion to tree-shaped invasion. Chapter 2 identifies the special design feature of distributing energy storage in time-dependent fashion on a territory, when the energy flows by fluid flow from a concentrated source to points (users) distributed equidistantly on the area. The challenge in this chapter is to determine the architecture of distributed energy storage. The chief conclusion is that the finite amount of storage material should be distributed proportionally with the distribution of the flow rate of heating agent arriving on the area. The total time needed by the source stream to ‘invade’ the area is cumulative (the sum of the storage times required at each storage site), and depends on the energy distribution paths and the sequence in which the users are served by the source stream. Chapter 3 shows theoretically that the melting process consists of two phases: “invasion” thermal diffusion along the invading line, which is followed by “consolidation” as heat diffuses perpendicularly to the invading line. This chapter also reports the duration of both phases and the evolution of the melt layer around the invading line during the two-dimensional and three-dimensional invasion. It also shows that the amount of melted material increases in time according to a curve shaped as an S. These theoretical predictions are validated by means of numerical simulations in chapter 4. This chapter also shows that the heat transfer rate density increases (i.e., the S curve becomes steeper) as the complexity and number of degrees of freedom of the structure are increased, in accord with the constructal law. The optimal geometric features of the tree structure are detailed in this chapter. Chapter 5 documents a numerical study of time-dependent melting where the heat transfer is convection dominated, unlike in chapter 3 and 4 where the melting is ruled by pure conduction. In accord with constructal design, the search is for effective heat-flow architectures. The volume-constrained improvement of the designs for heat flow begins with assuming the simplest structure, where a single line serves as heat source. Next, the heat source is endowed with freedom to change its shape as it grows. The objective of the numerical simulations is to discover the geometric features that lead to the fastest melting process. The results show that the heat transfer rate density increases as the complexity and number of degrees of freedom of the structure are increased. Furthermore, the angles between heat invasion lines have a minor effect on the global performance compared to other degrees of freedom: number of branching levels, stem length, and branch lengths. The effect of natural convection in the melt zone is documented.

Properties of InAlN/GaN Heterostructures Prepared by Molecular Beam Epitaxy

InAlN thin films and InAlN/GaN heterostructures have been intensively studied over recent years due to their applications in a variety of devices, including high electron mobility transistors (HEMTs). However, the quality of InAlN remains relatively poor with basic material and structural characteristics remain unclear.

Molecular beam epitaxy (MBE) is used to synthesize the materials for this research, as MBE is a widely used tool for semiconductor growth but has rarely been explored for InAlN growth. X-ray photoelectron spectroscopy (XPS) is used to determine the electronic and chemical characteristics of InAlN surfaces. This tool is used for the first time in application to MBE-grown InAlN and heterostructures for the characterization of surface oxides, the bare surface barrier height (BSBH), and valence band offsets (VBOs).

The surface properties of InAlN are studied in relation to surface oxide characteristics and formation. First, the native oxide compositions are studied. Then, methods enabling the effective removal of the native oxides are found. Finally, annealing is explored for the reliable growth of surface thermal oxides.

The bulk properties of InAlN films are studied. The unintentional compositional grading in InAlN during MBE growth is discovered and found to be affected by strain and relaxation. The optical characterization of InAlN using spectroscopy ellipsometry (SE) is also developed and reveals that a two-phase InAlN model applies to MBE-grown InAlN due to its natural formation of a nanocolumnar microstructure. The insertion of an AlN interlayer is found to mitigate the formation of this microstructure and increases mobility of whole structure by fivefold.

Finally, the synthesis and characterization of InAlN/GaN HEMT device structures are explored. The density and energy distribution of surface states are studied with relationships to surface chemical composition and surface oxide. The determination of the VBOs of InAlN/GaN structures with different In compositions are discussed at last.

Preservation of beam emittance in the presence of ion motion in future high-energy plasma-wakefield-based colliders.

The preservation of beam quality in a plasma wakefield accelerator driven by ultrahigh intensity and ultralow emittance beams, characteristic of future particle colliders, is a challenge. The electric field of these beams leads to plasma ions motion, resulting in a nonlinear focusing force and emittance growth of the beam. We propose to use an adiabatic matching section consisting of a short plasma section with a decreasing ion mass to allow for the beam to remain matched to the focusing force. We use analytical models and numerical simulations to show that the emittance growth can be significantly reduced.

Application of Numerical Methods to Study Arrangement and Fracture of Lithium-Ion Microstructure

The focus of this work is to develop and employ numerical methods that provide characterization of granular microstructures, dynamic fragmentation of brittle materials, and dynamic fracture of three-dimensional bodies.

We first propose the fabric tensor formalism to describe the structure and evolution of lithium-ion electrode microstructure during the calendaring process. Fabric tensors are directional measures of particulate assemblies based on inter-particle connectivity, relating to the structural and transport properties of the electrode. Applying this technique to X-ray computed tomography of cathode microstructure, we show that fabric tensors capture the evolution of the inter-particle contact distribution and are therefore good measures for the internal state of and electronic transport within the electrode.

We then shift focus to the development and analysis of fracture models within finite element simulations. A difficult problem to characterize in the realm of fracture modeling is that of fragmentation, wherein brittle materials subjected to a uniform tensile loading break apart into a large number of smaller pieces. We explore the effect of numerical precision in the results of dynamic fragmentation simulations using the cohesive element approach on a one-dimensional domain. By introducing random and non-random field variations, we discern that round-off error plays a significant role in establishing a mesh-convergent solution for uniform fragmentation problems. Further, by using differing magnitudes of randomized material properties and mesh discretizations, we find that employing randomness can improve convergence behavior and provide a computational savings.

The Thick Level-Set model is implemented to describe brittle media undergoing dynamic fragmentation as an alternative to the cohesive element approach. This non-local damage model features a level-set function that defines the extent and severity of degradation and uses a length scale to limit the damage gradient. In terms of energy dissipated by fracture and mean fragment size, we find that the proposed model reproduces the rate-dependent observations of analytical approaches, cohesive element simulations, and experimental studies.

Lastly, the Thick Level-Set model is implemented in three dimensions to describe the dynamic failure of brittle media, such as the active material particles in the battery cathode during manufacturing. The proposed model matches expected behavior from physical experiments, analytical approaches, and numerical models, and mesh convergence is established. We find that the use of an asymmetrical damage model to represent tensile damage is important to producing the expected results for brittle fracture problems.

The impact of this work is that designers of lithium-ion battery components can employ the numerical methods presented herein to analyze the evolving electrode microstructure during manufacturing, operational, and extraordinary loadings. This allows for enhanced designs and manufacturing methods that advance the state of battery technology. Further, these numerical tools have applicability in a broad range of fields, from geotechnical analysis to ice-sheet modeling to armor design to hydraulic fracturing.