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).

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).

Fluid dynamics and platelet deposition in a model arterial stenosis

The objective of this study was to gain further understanding and elucidation of the fluid dynamic factors and flow-induced mechanisms of the thrombogenic process of platelet deposition onto, and possible subsequent embolization from, the walls of an arterial stenosis. This has been accomplished by measurement of the axial dependence of platelet deposition within a modeled arterial stenosis for a transitional flow and a completely laminar flow field. The stenotic region of the model was collagen-coated to simulate a damaged endothelial lining of an artery. Fluid dynamics within a stenosis was studied using qualitative flow visualization, and was further compared to the in vitro platelet deposition studies. Normalized platelet density (NPD) measurements indicate decreased levels of NPD in the high shear throat region of the stenosis for a Reynolds number of 300 and a drastic increase in NPD at the throat for a Reynolds number of 175. This study provides further understanding of the flow dynamic effects on thrombus development within a stenosis. ^

The nitrate to ammonia and ceramic (NAC) process for the denitration and immobilization of low-level radioactive liquid waste (LLW)

Hazardous radioactive liquid waste is the legacy of more than 50 years of plutonium production associated with the United States' nuclear weapons program. It is estimated that more than 245,000 tons of nitrate wastes are stored at facilities such as the single-shell tanks (SST) at the Hanford Site in the state of Washington, and the Melton Valley storage tanks at Oak Ridge National Laboratory (ORNL) in Tennessee. In order to develop an innovative, new technology for the destruction and immobilization of nitrate-based radioactive liquid waste, the United State Department of Energy (DOE) initiated the research project which resulted in the technology known as the Nitrate to Ammonia and Ceramic (NAC) process. However, inasmuch as the nitrate anion is highly mobile and difficult to immobilize, especially in relatively porous cement-based grout which has been used to date as a method for the immobilization of liquid waste, it presents a major obstacle to environmental clean-up initiatives. Thus, in an effort to contribute to the existing body of knowledge and enhance the efficacy of the NAC process, this research involved the experimental measurement of the rheological and heat transfer behaviors of the NAC product slurry and the determination of the optimal operating parameters for the continuous NAC chemical reaction process. Test results indicate that the NAC product slurry exhibits a typical non-Newtonian flow behavior. Correlation equations for the slurry's rheological properties and heat transfer rate in a pipe flow have been developed; these should prove valuable in the design of a full-scale NAC processing plant. The 20-percent slurry exhibited a typical dilatant (shear thickening) behavior and was in the turbulent flow regime due to its lower viscosity. The 40-percent slurry exhibited a typical pseudoplastic (shear thinning) behavior and remained in the laminar flow regime throughout its experimental range. The reactions were found to be more efficient in the lower temperature range investigated. With respect to leachability, the experimental final NAC ceramic waste form is comparable to the final product of vitrification, the technology chosen by DOE to treat these wastes. As the NAC process has the potential of reducing the volume of nitrate-based radioactive liquid waste by as much as 70 percent, it not only promises to enhance environmental remediation efforts but also effect substantial cost savings. ^

Design and performance assessment of aortic heart valve for tissue engineering

Current artificial heart valves are classified as mechanical and bioprosthetic. An appealing pathway that promises to overcome the shortcomings of commercially available heart valves is offered by the interdisciplinary approach of cardiovascular tissue engineering. However, the mechanical properties of the Tissue Engineering Heart Valves (TEHV) are limited and generally fail in the long-term use. To meet this performance challenge novel biodegradable triblock copolymer poly(ethylene oxide)-polypropylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO or F108) crosslinked to Silk Fibroin (F108-SilkC) to be used as tri-leaflet heart valve material was investigated. ^ Synthesis of ten polymers with varying concentration and thickness (55 µm, 75 µm and 100 µm) was achieved via a covalent crosslinking scheme using bifunctional polyethylene glycol diglycidyl ether (PEGDE). Static and fatigue testing were used to assess mechanical properties of films, and hydrodynamic testing was performed to determine performance under a simulated left ventricular flow regime. The crosslinked copolymer (F108-Silk C) showed greater flexibility and resilience, but inferior ultimate tensile strength, by increasing concentration of PEGDE. Concentration molar ratio of 80:1 (F108: Silk) and thickness of 75 µm showed longer fatigue life for both tension-tension and bending fatigue tests. Four valves out of twelve designed satisfactorily complied with minimum performance requirement ISO 5840, 2005. ^ In conclusion, it was demonstrated that the applicability of a degradable polymer in conjugation with silk fibroin for tissue engineering cardiovascular use, specifically for aortic valve leaflet design, met the performance demands. Thinner thicknesses (t<75 µm) in conjunction with stiffness lower than 320 MPa (80:1, F108: Silk) are essential for the correct functionality of proposed heart valve biomaterial F108-SilkC. Fatigue tests were demonstrated to be a useful tool to characterize biomaterials that undergo cyclic loading. ^

Development of an integrated surface and subsurface model of Everglades National Park

An integrated surface-subsurface hydrological model of Everglades National Park (ENP) was developed using MIKE SHE and MIKE 11 modeling software. The model has a resolution of 400 meters, covers approximately 1050 square miles of ENP, includes 110 miles of drainage canals with a variety of hydraulic structures, and processes hydrological information, such as evapotranspiration, precipitation, groundwater levels, canal discharges and levels, and operational schedules. Calibration was based on time series and probability of exceedance for water levels and discharges in the years 1987 through 1997. Model verification was then completed for the period of 1998 through 2005. Parameter sensitivity in uncertainty analysis showed that the model was most sensitive to the hydraulic conductivity of the regional Surficial Aquifer System, the Manning's roughness coefficient, and the leakage coefficient, which defines the canal-subsurface interaction. The model offers an enhanced predictive capability, compared to other models currently available, to simulate the flow regime in ENP and to forecast the impact of topography, water flows, and modifying operation schedules.