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

A mathematical resolution to log transformations and the binning effect in applied processing of data in flow cytometry

This dissertation develops a new mathematical approach that overcomes the effect of a data processing phenomenon known as "histogram binning" inherent to flow cytometry data. A real-time procedure is introduced to prove the effectiveness and fast implementation of such an approach on real-world data. The histogram binning effect is a dilemma posed by two seemingly antagonistic developments: (1) flow cytometry data in its histogram form is extended in its dynamic range to improve its analysis and interpretation, and (2) the inevitable dynamic range extension introduces an unwelcome side effect, the binning effect, which skews the statistics of the data, undermining as a consequence the accuracy of the analysis and the eventual interpretation of the data. Researchers in the field contended with such a dilemma for many years, resorting either to hardware approaches that are rather costly with inherent calibration and noise effects; or have developed software techniques based on filtering the binning effect but without successfully preserving the statistical content of the original data. The mathematical approach introduced in this dissertation is so appealing that a patent application has been filed. The contribution of this dissertation is an incremental scientific innovation based on a mathematical framework that will allow researchers in the field of flow cytometry to improve the interpretation of data knowing that its statistical meaning has been faithfully preserved for its optimized analysis. Furthermore, with the same mathematical foundation, proof of the origin of such an inherent artifact is provided. These results are unique in that new mathematical derivations are established to define and solve the critical problem of the binning effect faced at the experimental assessment level, providing a data platform that preserves its statistical content. In addition, a novel method for accumulating the log-transformed data was developed. This new method uses the properties of the transformation of statistical distributions to accumulate the output histogram in a non-integer and multi-channel fashion. Although the mathematics of this new mapping technique seem intricate, the concise nature of the derivations allow for an implementation procedure that lends itself to a real-time implementation using lookup tables, a task that is also introduced in this dissertation.

Integrated Surface-Ground Water Modeling in Wetlands With Improved Methods to Simulate Vegetative Resistance to Flow

This dissertation focused on developing an integrated surface – subsurface hydrologic simulation numerical model by programming and testing the coupling of the USGS MODFLOW-2005 Groundwater Flow Process (GWF) package (USGS, 2005) with the 2D surface water routing model: FLO-2D (O’Brien et al., 1993). The coupling included the necessary procedures to numerically integrate and verify both models as a single computational software system that will heretofore be referred to as WHIMFLO-2D (Wetlands Hydrology Integrated Model). An improved physical formulation of flow resistance through vegetation in shallow waters based on the concept of drag force was also implemented for the simulations of floodplains, while the use of the classical methods (e.g., Manning, Chezy, Darcy-Weisbach) to calculate flow resistance has been maintained for the canals and deeper waters. A preliminary demonstration exercise WHIMFLO-2D in an existing field site was developed for the Loxahatchee Impoundment Landscape Assessment (LILA), an 80 acre area, located at the Arthur R. Marshall Loxahatchee National Wild Life Refuge in Boynton Beach, Florida. After applying a number of simplifying assumptions, results have illustrated the ability of the model to simulate the hydrology of a wetland. In this illustrative case, a comparison between measured and simulated stages level showed an average error of 0.31% with a maximum error of 2.8%. Comparison of measured and simulated groundwater head levels showed an average error of 0.18% with a maximum of 2.9%.

Numerical groundwater flow modeling in the Wakal River basin, India

Increasing dependence on groundwater in the Wakal River basin, India, jeopardizes water supply sustainability. A numerical groundwater model was developed to better understand the aquifer system and to evaluate its potential in terms of quantity and replenishment. Potential artificial recharge areas were delineated using landscape and hydrogeologic parameters, Geographic Information System (GIS), and remote sensing. Groundwater models are powerful tools for recharge estimation when transmissivity is known. Proper recharge must be applied to reproduce field-measured heads. The model showed that groundwater levels could decline significantly if there are two drought years in every four years that result in reduced recharge, and groundwater withdrawal is increased by 15%. The effect of such drought is currently uncertain however, because runoff from the basin is unknown. Remote sensing and GIS revealed areas with slopes less than 5%, forest cover, and Normalized Difference Vegetative Index greater than 0.5 that are suitable recharge sites.

Efficient numerical methods for the random-field Ising model: Finite-size scaling, reweighting extrapolation, and computation of response functions

It was recently shown [Phys. Rev. Lett. 110, 227201 (2013)] that the critical behavior of the random-field Ising model in three dimensions is ruled by a single universality class. This conclusion was reached only after a proper taming of the large scaling corrections of the model by applying a combined approach of various techniques, coming from the zero-and positive-temperature toolboxes of statistical physics. In the present contribution we provide a detailed description of this combined scheme, explaining in detail the zero-temperature numerical scheme and developing the generalized fluctuation-dissipation formula that allowed us to compute connected and disconnected correlation functions of the model. We discuss the error evolution of our method and we illustrate the infinite limit-size extrapolation of several observables within phenomenological renormalization. We present an extension of the quotients method that allows us to obtain estimates of the critical exponent a of the specific heat of the model via the scaling of the bond energy and we discuss the self-averaging properties of the system and the algorithmic aspects of the maximum-flow algorithm used.

A SPH solver for simulating paramagnetic solid fluid interaction in the presence of an external magnetic field

Acknowledgment The first two authors wish to express their sincerest thanks to Iran National Science Foundation (INSF) for supporting this work under Contract Number 92021291.

Moduli Space of (0,2) Conformal Field Theories

In this thesis we study aspects of (0,2) superconformal field theories (SCFTs), which are suitable for compactification of the heterotic string. In the first part, we study a class of (2,2) SCFTs obtained by fibering a Landau-Ginzburg (LG) orbifold CFT over a compact K\"ahler base manifold. While such models are naturally obtained as phases in a gauged linear sigma model (GLSM), our construction is independent of such an embedding. We discuss the general properties of such theories and present a technique to study the massless spectrum of the associated heterotic compactification. We test the validity of our method by applying it to hybrid phases of GLSMs and comparing spectra among the phases. In the second part, we turn to the study of the role of accidental symmetries in two-dimensional (0,2) SCFTs obtained by RG flow from (0,2) LG theories. These accidental symmetries are ubiquitous, and, unlike in the case of (2,2) theories, their identification is key to correctly identifying the IR fixed point and its properties. We develop a number of tools that help to identify such accidental symmetries in the context of (0,2) LG models and provide a conjecture for a toric structure of the SCFT moduli space in a large class of models. In the final part, we study the stability of heterotic compactifications described by (0,2) GLSMs with respect to worldsheet instanton corrections to the space-time superpotential following the work of Beasley and Witten. We show that generic models elude the vanishing theorem proved there, and may not determine supersymmetric heterotic vacua. We then construct a subclass of GLSMs for which a vanishing theorem holds.

Dr. Price measuring surface water flow at SRS-1

http://digitalcommons.fiu.edu/fce_lter_photos/1338/thumbnail.jpg

Graduate student David Lagomasino downloading surface water flow data at SRS-1

http://digitalcommons.fiu.edu/fce_lter_photos/1339/thumbnail.jpg

Field investigation of contactless displacement measurement using computer vision systems for civil engineering applications

Much of the bridge stock on major transport links in North America and Europe was constructed in the 1950s and 1960s and has since deteriorated or is carrying loads far in excess of the original design loads. Structural Health Monitoring Systems (SHM) can provide valuable information on the bridge capacity but the application of such systems is currently limited by access and bridge type. This paper investigates the use of computer vision systems for SHM. A series of field tests have been carried out to test the accuracy of displacement measurements using contactless methods. A video image of each test was processed using a modified version of the optical flow tracking method to track displacement. These results have been validated with an established measurement method using linear variable differential transformers (LVDTs). The results obtained from the algorithm provided an accurate comparison with the validation measurements. The calculated displacements agree within 2% of the verified LVDT measurements, a number of post processing methods were then applied to attempt to reduce this error.

Combined Conduction-Convection-Radiation Heat Transfer of Slip Flow inside a Micro-channel Filled with a Cellular Porous Material

Combined conduction–convection–radiation heat transfer is investigated numerically in a micro-channel filled with a saturated cellular porous medium, with the channel walls held at a constant heat flux. Invoking the velocity slip and temperature jump, the thermal behaviour of the porous–fluid system are studied by considering hydrodynamically fully developed flow and applying the Darcy–Brinkman flow model. One energy equation model based on the local thermal equilibrium condition is adopted to evaluate the temperature field within the porous medium. Combined conduction and radiation heat transfer is treated as an effective conduction process with a temperature-dependent effective thermal conductivity. Results are reported in terms of the average Nusselt number and dimensionless temperature distribution, as a function of velocity slip coefficient, temperature jump coefficient, porous medium shape parameter and radiation parameters. Results show that increasing the radiation parameter (Tr)(Tr) and the temperature jump coefficient flattens the dimensionless temperature profile. The Nusselt numbers are more sensitive to the variation in the temperature jump coefficient rather than to the velocity slip coefficient. Such that for high porous medium shape parameter, the Nusselt number is found to be independent of velocity slip. Furthermore, it is found that as the temperature jump coefficient increases, the Nusselt number decrease. In addition, for high temperature jump coefficients, the Nusselt number is found to be insensitive to the radiation parameters and porous medium shape parameter. It is also concluded that compared with the conventional macro-channels, wherein using a porous material enhances the rate of heat transfer (up to about 40 % compared to the clear channel), insertion of a porous material inside a micro-channel in slip regime does not effectively enhance the rate of heat transfer that is about 2 %.

Mechanochemical synthesis of TiO2/NiFe2O4 magnetic catalysts for operation under RF field

Composite NiFe2O4–TiO2 magnetic catalysts were prepared by mechanochemical synthesis from a mixture of titania supported nickel ferrite nanoparticles and P25 titania (Evonic). The former provides fast and efficient heating under radiofrequency field, while the latter serves as an active catalyst or catalyst support. The highest heating rate was observed over a catalyst prepared for a milling time of 30 min. The catalytic activity was measured over the sulfated composite catalysts in the condensation of aniline and 3-phenylbutyric acid in a stirred tank reactor and in a continuous RF heated flow reactor in the 140–170 °C range. The product yield of 47% was obtained over the sulfated P25 titania catalyst in the flow reactor.

Ecological effects of flow and temperature regimes on fish communities in temperate rivers

Les organismes aquatiques sont adaptés à une grande variabilité hydrique et thermique des rivières. Malgré ceci, la régulation des eaux suscite des changements aux débits qui peuvent provoquer des impacts négatifs sur la biodiversité et les processus écologiques en rivière. Celle-ci peut aussi causer des modifications au niveau des régimes thermiques et des caractéristiques de l’habitat du poisson. Des données environnementales et biologiques décrivant l’habitat du poisson existent, mais elles sont incomplètes pour plusieurs rivières au Canada et de faible qualité, limitant les relations quantitatives débit-température-poissons à un petit nombre de rivières ou à une région étudiée. La recherche menée dans le cadre de mon doctorat concerne les impacts de la génération d'hydroélectricité sur les rivières; soit les changements aux régimes hydriques et thermiques reliés à la régulation des eaux sur la variation des communautés ichtyologiques qui habitent les rivières régulées et naturelles au Canada. Suite à une comparaison d’échantillonnage de pêche, une méthode constante pour obtenir des bons estimés de poisson (richesse, densité et biomasse des espèces) a été établie pour évaluer la structure de la communauté de poissons pour l’ensemble des rivières ciblées par l’étude. Afin de mieux comprendre ces changements environnementaux, les principales composantes décrivant ces régimes ont été identifiées et l’altération des régimes hydriques pour certaines rivières régulées a été quantifiée. Ces résultats ont servi à établir la relation significative entre le degré de changement biotique et le degré de changement hydrique pour illustrer les différences entre les régimes de régulation. Pour faire un complément aux indices biotiques déjà calculés pour l’ensemble des communautés de poissons (diversité, densité et biomasse des espèces par rivière), les différences au niveau des guildes de poissons ont été quantifiées pour expliquer les divers effets écologiques dus aux changements de régimes hydriques et thermiques provenant de la gestion des barrages. Ces derniers résultats servent à prédire pour quels traits écologiques ou groupes d’espèces de poissons les composantes hydriques et thermiques sont importantes. De plus, ces derniers résultats ont servi à mettre en valeur les variables décrivant les régimes thermiques qui ne sont pas toujours inclues dans les études hydro-écologiques. L’ensemble des résultats de cette thèse ont des retombées importantes sur la gestion des rivières en évaluant, de façon cohérente, l’impact de la régulation des rivières sur les communautés de poissons et en développant des outils de prévision pour la restauration des écosystèmes riverains.

Modelling the Sun's small scale global photospheric magnetic field

We present a new model for the Sun's global photospheric magnetic field during a deep minimum of activity, in which no active regions emerge. The emergence and subsequent evolution of small- scale magnetic features across the full solar surface is simulated, subject to the influence of a global supergranular flow pattern. Visually, the resulting simulated magnetograms reproduce the typical structure and scale observed in quiet Sun magnetograms. Quantitatively, the simulation quickly reaches a steady state, resulting in a mean field and flux distribution that are in good agreement with those determined from observations. A potential coronal magnetic field is extrapolated from the simulated full Sun magnetograms, to consider the implications of such a quiet photospheric magnetic field on the corona and inner heliosphere. The bulk of the coronal magnetic field closes very low down, in short connections between small-scale features in the simulated magnetic network. Just 0.1% of the photospheric magnetic flux is found to be open at 2:5 Rʘ, around 10 - 100 times less than that determined for typical HMI synoptic map observations. If such conditions were to exist on the Sun, this would lead to a significantly weaker interplanetary magnetic field than is presently observed, and hence a much higher cosmic ray flux at Earth.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.