The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re=600-900 and phase change material particle concentrations ¡Ü0.25 , as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (0.5-0.7) slurry flow. By using the two newly-defined parameters, named effectiveness factor and performance index, respectively, it is found that there exists an optimal relation between the channel design parameters, particle volume fraction, Reynolds number, and the wall heat flux. The influence of the particle volume fraction, particle size, and the particle viscosity, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed.

Dans un contexte mondial où l’accès à l’énergie est un problème de premier plan, l’exploitation des courants de marée avec des hydroliennes revête un intérêt certain. Les écoulements dans les zones à fort potentiel énergétique propices à l’installation d’hydroliennes sont souvent fortement turbulents. Or la turbulence ambiante impacte fortement l’hydrodynamique avoisinante et le fonctionnement de la turbine. Une prédiction fine de la turbulence et du sillage est fondamentale pour l'optimisation d'une ferme d'hydroliennes. Un modèle de simulation de l'écoulement autour de la turbine doit donc être précis et tenir compte de la turbulence ambiante. Un outil basé sur la méthode de Boltzmann sur réseau (LBM) est utilisé à ces fins, en association avec une approche de simulation des grandes échelles (LES). La LBM est une méthode instationnaire de modélisation d’écoulement fluide. Une méthode de génération de turbulence synthétique est implémentée afin de prendre en compte la turbulence ambiante des sites hydroliens. Les géométries complexes, potentiellement en mouvement, sont modélisées avec la méthode des frontières immergées (IBM). La mise en place d’un modèle de paroi est réalisée afin de réduire le cout en calcul du modèle. Ces outils sont ensuite utilisés pour modéliser en LBM-LES une hydrolienne dans un environnement turbulent. Les calculs, réalisés à deux taux de turbulence différents, sont comparés avec des résultats expérimentaux et des résultats NS-LES. Les modélisations LBM-LES sont ensuite utilisées pour analyser le sillage de l'hydrolienne. Il est notamment observé qu'un faible taux de turbulence impacte de manière significative la propagation des tourbillons de bout de pale. / In a global context where access to energy is a major problem, the exploitation of tidal currents with tidal turbines is of particular interest. Flows in areas with high energy potential suitable for the installation of tidal turbines are often highly turbulent. However, the ambient turbulence has a strong impact on the surrounding hydrodynamics and the turbine operation. A precise prediction of turbulence and wake is fundamental to the optimization of a tidal farm. A numerical model of the flow around the turbine must therefore be accurate and take into account the ambient turbulence. A tool based on the Lattice Boltzmann Method (LBM) is used for this purpose, in combination with a Large Eddy Simulation (LES) approach. The LBM is an unsteady method for modelling fluid flows. A synthetic turbulence method is implemented to take into account the ambient turbulence of tidal sites. Complex geometries, potentially in motion, are modelled using the Immersed Boundary Method (IBM). The implementation of a wall model is carried out in order to reduce the cost of the simulations. These tools are then used to model a turbine in a turbulent environment. The calculations, performed at two different turbulence rates, are compared with experimental and NS-LES results. The LBM-LES models are then used to analyze the wake of the turbine. In particular, it is observed that a low turbulence rate has a significant impact on the propagation of tip-vortices.

Les extensions bifluides des méthodes Lattice Boltzmann à frontière libre utilisent généralement des pseudopotentiels microscopiques pour modéliser l'interface. Nous avons choisi d'orienter nos recherches vers une méthode Lattice Boltzmann à capture d'interface où la fraction massique d'un des deux fluides, inconnue, est transportée. De nombreux travaux ont montré les difficultés des méthodes Lattice Boltzmann à traiter des systèmes bifluides, et ce d'autant plus que le ratio de densité est important. Nous expliquerons l'origine de ces problèmes en mettant en évidence le manque de diffusion numérique pour capturer précisément les discontinuités de contact. Pour régler cet obstacle, nous proposerons une formulation Arbitrary Lagrangian Eulerian (ALE) des méthodes Lattice Boltzmann. Cela permet de séparer le traitement des ondes matérielles de celui des ondes de pression. Une fois l'étape ALE terminée, une phase de projection ramène les variables sur la grille eulérienne de calcul initiale. Nous expliquons comment obtenir une procédure de projection ayant une précision d'ordre 2 et une interface fine et dépourvue d'oscillations. Il sera montré que la fraction massique satisfait un principe du maximum discret et qu'elle reste donc entre 0 et 1. Les simulations numériques sont en accord avec la théorie. Même si notre méthode n'est pour le moment utilisée que pour simuler des écoulements de fluides non visqueux (Equations d'Euler), nous sommes convaincus qu'elle pourra être étendue à des simulations d'écoulements bifluides visqueux. / Two-fluid extensions of Lattice Boltzmann methods with free boundaries usually consider ``microscopic'' pseudopotential interface models. In this paper, we rather propose an interface-capturing Lattice Boltzmann approach where the mass fraction variable is considered as an unknown and is advected. Several works have reported the difficulties of LBM methods to deal with such two-fluid systems especially for high-density ratio configurations. This is due to the mixing nature of LBM, as with Flux vector splitting approaches for Finite Volume methods. We here give another explanation of the lack of numerical diffusion of Lattice Boltzmann approaches to accurately capture contact discontinuities. To fix the problem, we propose an arbitrary Lagrangian-Eulerian (ALE) formulation of Lattice-Boltzmann methods. In the Lagrangian limit, it allows for a proper separated treatment of pressure waves and advection phenomenon. After the ALE solution, a remapping (advection) procedure is necessary to project the variables onto the Eulerian Lattice-Boltzmann grid.We explain how to derive this remapping procedure in order to get second-order accuracy and achieve sharp stable oscillation-free interfaces. It has been shown that mass fractions variables satisfy a local discrete maximum principle and thus stay in the range $[0,1]$. The theory is supported by numerical computations of rising bubbles (without taking into account surface tension at this current state of development).Even if our methods are currently used for inviscid flows (Euler equations) by projecting the discrete distributions onto equilibrium ones at each time step, we believe that it is possible to extend the framework formulation for multifluid viscous problems. This will be at the aim of a next work.

Indiana University-Purdue University Indianapolis (IUPUI) / The thesis includes two application studies of computational fluid dynamics. The first is new and efficient drug delivery to the posterior part of the eye, a growing health necessity worldwide. Current treatment of eye diseases, such as age-related macular degeneration (AMD), relies on repeated intravitreal injections of drug-containing solutions. Such a drug delivery has significant cant drawbacks, including short drug life, vital medical service, and high medical costs. In this study, we explore a new approach of controlled drug delivery by introducing unique porous implants. Computational modeling contains physiological and anatomical traits. We simulate the IgG1 Fab drug delivery to the posterior eye to evaluate the effectiveness of the porous implants to control the drug delivery. The computational model was validated by established computation results from independent studies and experimental data. Overall, the results indicate that therapeutic drug levels in the posterior eye are sustained for eight weeks, similar to those performed with intravitreal injection of the same drug. We evaluate the effects of the porous implant on the time evaluation of the drug concentrations in the sclera, choroid, and retina layers of the eye. Subsequent simulations were carried out with varying porosity values of a porous episcleral implant. Our computational results reveal that the time evolution of drug concentration is distinctively correlated to drug source location and pore size. The response of this porous implant for controlled drug delivery applications was examined. A correlation between porosity and fluid properties for the porous implants was revealed in this study. The second application lays in the computational modeling of the oscillating

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第18271号 / 工博第3863号 / 新制||工||1592(附属図書館) / 31129 / 京都大学大学院工学研究科航空宇宙工学専攻 / (主査)教授 稲室 隆二, 教授 泉田 啓, 教授 青木 一生 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第19681号 / 工博第4136号 / 新制||工||1638(附属図書館) / 32717 / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 西脇 眞二, 教授 稲室 隆二, 教授 松原 厚 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

Micromechanics of solidluid interactions can play a key role controlling macro-scale engineering behavior of granular media. The main objective of this study is to numerically investigate the micromechanics involved in solidluid mixtures to develop a better understanding of the macroscopic behavior of granular media for different applications. This is accomplished by developing a numerical model coupling the Discrete Element Method (DEM) and the Lattice Boltzmann Method (LBM) and employing it to study three distinct yet interrelated applications throughout the course of this research. In the first application, the DEM model is used to provide a clear relationship between energy dissipated by micro-scale mechanisms versus the traditional engineering definition based on macro-scale (continuum) parameters to develop a better understanding for the frictional behavior of granular media. Macroscopic frictional behavior of granular materials is of great importance for studying several complex problems such as fault slip and landslides. In the second application, the DEM-LBM model is employed for studying the undrained condition of dense granular media. While the majority of previous modeling approaches did not realistically represent non-uniform strain conditions that exist in geomechanical problems, including the LBM in the proposed model offers a realistic approach to simulate the undrained condition since the fluid can locally conserve the system volume. For the third application, the DEM-LBM model is used to study discontinuous shear thickening in a dense solidluid suspension. Shear thickening in a fluid occurs when the viscosity of the fluid increases with increasing applied strain rate. The DEM-LBM results for discontinuous shear thickening were compared to experimental data and proved to be an accurate approach at reproducing this phenomenon. The validated DEM-LBM model is then used to develop a physics-based constitutive model for discontinuous shear thickening-shear thinning in granular medialuid suspension. A closedorm model is then calibrated using the DEM-LBM model and validated against existing experimental test results reported in the literature. Findings of this research demonstrate how micromechanical modeling can be employed to address challenging problems in granular media involving solidluid interaction.

In this work, a two-phase lattice Boltzmann method (LBM) approach is implemented to investigate the hydrodynamic behavior of a single droplet impingement on a dry surface. LBM is a recently developed powerful technique to compute a wide range of fluid flow problems, especially in applications involving interfacial dynamics and complex geometries. Instead of solving the non-linear Navier-Stokes equations, which are complicated partial differential equations, LBM solves a set of discretized linear equations, which are easy to implement and parallelize. The fundamental idea of LBM is to recover the macroscopic properties of the fluid which obeys Navier-Stokes equations, by using simplified kinetic equations that incorporate the essential physics at the microscopic level. Considering the numerical instability induced by large density difference between two phases during the LBM simulations, the particular LBM scheme used in this study has its benefits when dealing with high density ratios. All the simulations are conducted for density ratio up to 50 in a three-dimensional Cartesian coordinate system, and three important dimensionless numbers, namely Weber number, Reynolds number and Ohnesorge number, are used for this study. To validate this multiphase LBM approach, several benchmark tests are conducted. First, the angular frequency of an oscillating droplet is calculated and compared with the corresponding theoretical value. Errors are found to be within 6.1% for all the cases. Secondly, simulations of binary droplet collisions are conducted in the range of 20

Computational fluid dynamics, abbreviated CFD, is a valuable tool for several engineering applications. Applications such as heating, cooling or drying are some examples of areas where CFD is used. The Lattice Boltzmann method, abbreviated LBM, is a popular method for CFD simulations due to its fast simulation time in comparison to traditional methods like the finite volume method. ESS is a company that has supplied LBM simulation software aimed at perfecting the baking process of car carrosserie. The aim of this thesis is to verify and validate their software against earlier works, to ensure that the solver can capture the physical reality of an impinging jet. The verification and validation are conducted over three different cases: experimental free jet, direct numerical simulation and a large eddy simulation of an impinging jet. Identical digital models were created for each case and then simulated with similar conditions for comparison. All three cases were tested with Richardson extrapolation to ensure grid convergence. The Richardson method showed that the highest error among the cases was 2.02 %. The results shows that the LBM solver can accurately predict the flow and thermal field compared to the three cases. The LBM solver is considered verified and validated for flowfield and thermal simulations. / Numeriska strömningsberäkningar är ett värdefullt verktyg för flera ingenjörsapplikationer. Applikationer som uppvärmning, kylning eller torkning är några exempel på områden där CFD används. Lattice Boltzmann-metoden, förkortat LBM, är en populär metod för CFD-simuleringar på grund av dess snabba simuleringshastighet jämfört med traditionella metoder som finita volymsmetoden. ESS är ett företag som har skapat LBM- simuleringsprogramvaran med målet att förbättra bakningsprocessen för bilkarosserier. Målet med denna avhandling är att verifiera och validera deras programvara mot tidigare arbeten för att säkerställa att lösningsmetoden kan fånga den fysiska verkligheten hos en påverkande jetstråle.Verifieringen och valideringen utförs över tre olika fall: Experimentell fri stråle, Direkt numerisk simulering och en Large eddy-simulering av en påverkande jetstråle. Identiska digitala modeller skapades för varje fall och simulerades sedan under liknande förhållanden för att kunna jämföra resultaten. Alla tre fall testades med Richardson-extrapolation för att säkerställa nätkonvergens. Richardson-metoden visade att den högsta felet bland fallen var 2,02 %. Resultaten visar att LBM-lösaren kan förutsäga flödes och termiska fältet noggrant jämfört med de tre fallen. LBM-lösaren anses vara verifierad och validerad för flödesfälts- och termiska simuleringar.

No description available.