Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Implementation of X-ray diffraction in the measurement of residual thermal strains

Microelectronic systems are multi-material, multi-layer structures, fabricated and exposed to environmental stresses over a wide range of temperatures. Thermal and residual stresses created by thermal mismatches in films and interconnections are a major cause of failure in microelectronic devices. Due to new device materials, increasing die size and the introduction of new materials for enhanced thermal management, differences in thermal expansions of various packaging materials have become exceedingly important and can no longer be neglected. X-ray diffraction is an analytical method using a monochromatic characteristic X-ray beam to characterize the crystal structure of various materials, by measuring the distances between planes in atomic crystalline lattice structures. As a material is strained, this interplanar spacing is correspondingly altered, and this microscopic strain is used to determine the macroscopic strain. This thesis investigates and describes the theory and implementation of X-ray diffraction in the measurement of residual thermal strains. The design of a computer controlled stress attachment stage fully compatible with an Anton Paar heat stage will be detailed. The stress determined by the diffraction method will be compared with bimetallic strip theory and finite element models.

Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Tubular micro- and nanostructures of TCO materials grown by a vapor-solid method

Microtubes and rods with nanopipes of transparent conductive oxides (TCO), such as SnO_2, TiO_2, ZnO and In_2O_3, have been fabricated following a vapor-solid method which avoids the use of catalyst or templates. The morphology of the as-grown tubular structures varies as a function of the precursor powder and the parameters employed during the thermal treatments carried out under a controlled argon flow. These materials have been also doped with different elements of technological interest (Cr, Er, Li, Zn, Sn). Energy Dispersive X-ray Spectroscopy (EDS) measurements show that the concentration of the dopants achieved by the vapor-solid method ranges from 0.5 to _3 at.%. Luminescence of the tubes has been analyzed, with special attention paid to the influence of the dopants on their optical properties. In this work, we summarize and discuss some of the processes involved not only in the anisotropic growth of these hollow micro and nanostructures, but also in their doping.

Limits to sustained energy intake. XXIII. Does heat dissipation capacity limit the energy budget of lactating bank voles?

© 2016. Published by The Company of Biologists Ltd.

Modeling of a Heat-Induced Buckling of Plates Using the Mesh-free Method

In the process of engineering design of structural shapes, the flat plate analysis results can be generalized to predict behaviors of complete structural shapes. In this case, the purpose of this project is to analyze a thin flat plate under conductive heat transfer and to simulate the temperature distribution, thermal stresses, total displacements, and buckling deformations. The current approach in these cases has been using the Finite Element Method (FEM), whose basis is the construction of a conforming mesh. In contrast, this project uses the mesh-free Scan Solve Method. This method eliminates the meshing limitation using a non-conforming mesh. I implemented this modeling process developing numerical algorithms and software tools to model thermally induced buckling. In addition, convergence analysis was achieved, and the results were compared with FEM. In conclusion, the results demonstrate that the method gives similar solutions to FEM in quality, but it is computationally less time consuming.

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.

Crack propagation in non-homogenous materials: Evaluation of mixed-mode SIFs, T-stress and kinking angle using a variant EFG method

A new variant of the Element-Free Galerkin (EFG) method, that combines the diffraction method, to characterize the crack tip solution, and the Heaviside enrichment function for representing discontinuity due to a crack, has been used to model crack propagation through non-homogenous materials. In the case of interface crack propagation, the kink angle is predicted by applying the maximum tangential principal stress (MTPS) criterion in conjunction with consideration of the energy release rate (ERR). The MTPS criterion is applied to the crack tip stress field described by both the stress intensity factor (SIF) and the T-stress, which are extracted using the interaction integral method. The proposed EFG method has been developed and applied for 2D case studies involving a crack in an orthotropic material, crack along an interface and a crack terminating at a bi-material interface, under mechanical or thermal loading; this is done to demonstrate the advantages and efficiency of the proposed methodology. The computed SIFs, T-stress and the predicted interface crack kink angles are compared with existing results in the literature and are found to be in good agreement. An example of crack growth through a particle-reinforced composite materials, which may involve crack meandering around the particle, is reported.

Improved electrochemical performance of Sr2Fe1.5Mo0.4Nb0.1O6-δ -Sm0.2Ce0.8O2-δ composite cathodes by a one-pot method for intermediate temperature solid oxide fuel cells

In this paper, Sr₂Fe_1.5Mo_0.4Nb_0.1O_6-δ (SFMNb)-xSm_0.2Ce_0.8O_2-δ (SDC) (x = 0, 20, 30, 40, 50 wt%) composite cathode materials were synthesized by a one-pot combustion method to improve the electrochemical performance of SFMNb cathode for intermediate temperature solid oxide fuel cells (IT-SOFCs). The fabrication of composite cathodes by adding SDC to SFMNb is conducive to providing extended electrochemical reaction zones for oxygen reduction reactions (ORR). X-ray diffraction (XRD) demonstrates that SFMNb is chemically compatible with SDC electrolytes at temperature up to 1100 °C. Scanning electron microscope (SEM) indicates that the SFMNb-SDC composite cathodes have a porous network nanostructure as well as the single phase SFMNb. The conductivity and thermal expansion coefficient of the composite cathodes decrease with the increased content of SDC, while the electrochemical impedance spectra (EIS) exhibits that SFMNb-40SDC composite cathode has optimal electrochemical performance with low polarization resistance (R_p) on the La_0.9Sr_0.1Ga_0.8Mg_0.2O₃ electrolyte. The R_p of the SFMNb-40SDC composite cathode is about 0.047 Ω cm² at 800 °C in air. A single cell with SFMNb-40SDC cathode also displays favorable discharge performance, whose maximum power density is 1.22 W cm^-2 at 800 °C. All results indicate that SFMNb-40SDC composite material is a promising cathode candidate for IT-SOFCs.

Synthesis of Pr0.6Sr0.4FeO 3-δ -xCe0.9Pr0.1O 2-δ cobalt-free composite cathodes by a one-pot method for intermediate-temperature solid oxide fuel cells

Cobalt-free composite cathodes consisting of Pr_0.6Sr_0.4FeO _3-δ -xCe_0.9Pr_0.1O _2-δ (PSFO-xCPO, x = 0-50 wt%) have been synthesized using a one-pot method. X-ray diffraction, scanning electron microscopy, thermal expansion coefficient, conductivity, and polarization resistance (R _P ) have been used to characterize the PSFO-xCPO cathodes. Furthermore the discharge performance of the Ni-SSZ/SSZ/GDC/PSFO-xCPO cells has been measured. The experimental results indicate that the PSFO-xCPO composite materials fully consist of PSFO and CPO phases and posses a porous microstructure. The conductivity of PSFO-xCPO decreases with the increase of CPO content, but R _P of PSFO-40CPO shows the smallest value amongst all the samples. The power density of single cells with a PSFO-40CPO composite cathode is significantly improved compared with that of the PSFO cathode, exhibiting 0.43, 0.75, 1.08 and 1.30 W cm^-2 at 650, 700, 750 and 800 °C, respectively. In addition, single cells with the PSFO-40CPO composite cathode show a stable performance with no obvious degradation over 100 h when operating at 750 °C.

Data-Driven Dynamic Modeling of Coupled Thermal and Electric Outputs of Microturbines

Microturbines are among the most successfully commercialized distributed energy resources, especially when they are used for combined heat and power generation. However, the interrelated thermal and electrical system dynamic behaviors have not been fully investigated. This is technically challenging due to the complex thermo-fluid-mechanical energy conversion processes which introduce multiple time-scale dynamics and strong nonlinearity into the analysis. To tackle this problem, this paper proposes a simplified model which can predict the coupled thermal and electric output dynamics of microturbines. Considering the time-scale difference of various dynamic processes occuring within microturbines, the electromechanical subsystem is treated as a fast quasi-linear process while the thermo-mechanical subsystem is treated as a slow process with high nonlinearity. A three-stage subspace identification method is utilized to capture the dominant dynamics and predict the electric power output. For the thermo-mechanical process, a radial basis function model trained by the particle swarm optimization method is employed to handle the strong nonlinear characteristics. Experimental tests on a Capstone C30 microturbine show that the proposed modeling method can well capture the system dynamics and produce a good prediction of the coupled thermal and electric outputs in various operating modes.

Improved Real-time State-of-Charge Estimation of LiFePO4 Battery Based on a Novel Thermoelectric Model

Li-ion batteries have been widely used in electric vehicles, and battery internal state estimation plays an important role in the battery management system. However, it is technically challenging, in particular, for the estimation of the battery internal temperature and state-ofcharge (SOC), which are two key state variables affecting the battery performance. In this paper, a novel method is proposed for realtime simultaneous estimation of these two internal states, thus leading to a significantly improved battery model for realtime SOC estimation. To achieve this, a simplified battery thermoelectric model is firstly built, which couples a thermal submodel and an electrical submodel. The interactions between the battery thermal and electrical behaviours are captured, thus offering a comprehensive description of the battery thermal and electrical behaviour. To achieve more accurate internal state estimations, the model is trained by the simulation error minimization method, and model parameters are optimized by a hybrid optimization method combining a meta-heuristic algorithm and the least square approach. Further, timevarying model parameters under different heat dissipation conditions are considered, and a joint extended Kalman filter is used to simultaneously estimate both the battery internal states and time-varying model parameters in realtime. Experimental results based on the testing data of LiFePO4 batteries confirm the efficacy of the proposed method.

Thomson scattering on non-thermal atmospheric pressure plasma jets

To characterize non-thermal atmospheric pressure plasmas experimentally, a large variety of methods and techniques is available, each having its own specific possibilities and limitations. A rewarding method to investigate these plasma sources is laser Thomson scattering. However, that is challenging. Non-thermal atmospheric pressure plasmas (gas temperatures close to room temperature and electron temperatures of a few eV) have usually small dimensions (below 1 mm) and a low degree of ionization (below 10^-4). Here an overview is presented of how Thomson scattering can be applied to such plasmas and used to measure directly spatially and temporally resolved the electron density and energy distribution. A general description of the scattering of photons and the guidelines for an experimental setup of this active diagnostic are provided. Special attention is given to the design concepts required to achieve the maximum signal photon flux with a minimum of unwanted signals. Recent results from the literature are also presented and discussed.

Development of a Rational Characterization Method for Iowa Fly Ash, HR-286, Annual Progress Report 2, 1987

The following report summarizes research activities on the project for the period December 1, 1986 to November 30, 1987. Research efforts for the second year deviated slightly from those described in the project proposal. By the end of the second year of testing, it was possible to begin evaluating how power plant operating conditions influenced the chemical and physical properties of fly ash obtained from one of the monitored power plants (Ottumwa Generating Station, OGS). Hence, several of the tasks initially assigned to the third year of the project (specifically tasks D, E, and F) were initiated during the second year of the project. Manpower constraints were balanced by delaying full scale implementation of the quantitative X-ray diffraction and differential thermal analysis tasks until the beginning of the third year of the project. Such changes should have little bearing on the outcome of the overall project.