Thesis (MTech (Mechanical Engineering))--Cape Peninsula University of Technology, 2016. / This work presents the computational uid dynamics (CFD) analysis of a light road vehicle. Simulations are conducted using the lattice Boltzmann method (LBM) with the wall adapting local eddy (WALE) turbulence model. Simulations include and compare the use of a rolling road, rotating wheels, adaptive re nement as well as showing comparison with a Reynolds-averaged Navier-Stokes (RANS) solver and the Spalart- Allmaras (SA) turbulence model. The lift coe cient of the vehicle for the most part was seen to show a much greater di erence and inconsistencies when compared to drag from the comparisons of solvers, turbulence models, re nement and the e ect of rolling road. Determining the drag of a road vehicle can be easily achieved and veri ed using multiple solvers and methods, however, the lift coe cient and its validation require a greater understanding of the vehicle ow eld as well as the solvers, turbulence models and re nement levels capable of correctly simulating the turbulent regions around a vehicle. Using the presented method, it was found that the optimisation of vehicle aerodynamics can easily be done alongside the design evolution from initial low-drag shapes to the nal detail design, ensuring aerodynamic characteristics are controlled with aesthetic change.

La modélisation des écoulements et du transport dans les milieux poreux est un domaine actif pour, notamment, progresser dans la compréhension du transfert des polluants dans les sols. Les sols présentent fréquemment des hétérogénéités comme des macropores (provoqués par la faune, la flore ou des fissures) et un certain nombre de modèles numériques utilisent les concepts de double ou de multi-perméabilité pour tenir compte de tous les types d’écoulements susceptibles de coexister dans de tels systèmes. Cependant, les modèles classiques semblent sous-estimer l’effet de la macroporosité sur l’écoulement et le transfert préférentiels et restreindre la zone d’écoulement préférentiel au seul volume occupé par la macroporosité. Diverses études expérimentales antérieures à cette thèse ont questionné cette hypothèse. Cette étude se propose de comprendre l’établissement de l’écoulement et du transport préférentiel et en particulier les mécanismes d’échange d’eau et de masse entre un macropore et une matrice poreuse environnante en condition saturée. Pour cela, des traçages de l’eau sont réalisés pour un milieu poreux modèle constitué de billes de verre, traversé par un macropore synthétique et mis en place en colonnes de laboratoire. Elution et transfert dans les colonnes sont caractérisés par suivi de la concentration en sortie et par imagerie par résonance magnétique. Un modèle numérique développé sur la base de la méthode de Boltzmann sur réseau est utilisé pour simuler numériquement des écoulements dans un système macroporé et identifier les mécanismes d’écoulements préférentiels à l’échelle de pores. Les données expérimentales montrent que le transfert du traceur est fortement dépendant du débit d’injection ainsi que du coefficient de diffusion dans l’eau. À fort débit, le transfert semble s’effectuer exclusivement dans le macropore, avec très peu d’échange avec la matrice. Pour des débits plus faibles, la percée présente une inflexion suivie d’un pic. Les images IRM montrent alors un échange significatif de traceur entre le macropore et la matrice poreuse environnante. Les simulations numériques sont utilisées pour calculer le champ de vitesse de l’écoulement dans le système en fonction du débit. Les modélisations numériques montrent que l’écoulement préférentiel est étendu dans la matrice poreuse sur une zone de même dimension que le diamètre moyen des grains indépendamment de la taille du macropore et du débit, dans la gamme de débits simulés. Ces résultats expérimentaux et numériques montrent que l’influence du macropore sur les transferts doit être étendue dans la matrice poreuse sur une zone de la taille des grains pour l’écoulement et sur une zone dépendant du coefficient de diffusion du traceur ainsi que du temps de séjour moyen de celui-ci pour le transfert des solutés / Flow and transport modeling through porous media is of primary concern nowadays, especially in order to progress in the understanding of pollutant transfers through soils. Soils present frequently heterogeneities such as macropores (caused by fauna, flora or cracks) and several numerical models use double or multi permeability concepts in order to take into account all flow types that may exist in such porous systems. Nevertheless, classical models seem underestimate the macropore effect on preferential flow and transport by restricting the preferential flow zone only to the volume occupied by the macroporosity. Various experimental studies prior to this thesis have questioned this hypothesis. This study proposes to understand the establishment of preferential flow and transport and in particular the mechanism of flow and solute exchanges between a synthetic macropore and a surrounding porous matrix in saturated condition. For this purpose, water tracing are realized for a model porous media constituted by glass beads, crossed by a synthetic macropore and implemented in laboratory columns. Breakthrough and transport in columns are characterized by monitoring the concentration at the end of the column by magnetic nuclear resonance. A numerical model developed on the basis of lattice-Boltzmann method is used to simul ate flow in macroporous system and identify preferential flow mechanisms at pore scale. Experimental data show that tracer transport is strongly dependent on injection flow rate and the diffusion coefficient in water. At high flow rate, the transport seems to occur exclusively in the macropore, with very little masse exchange with the porous matrix. At lower flow rates, the breakthrough exhibits an inflexion followed by a peak. The MRI images show a significant mass exchange of tracer between the macropore and the surrounding porous matrix. The numerical simulations are used to calculate the flow field in a porous system as a function of flow rate. They show that preferential flow is extended in porous matrix into a zone of same dimension the mean diameter of beads regardless of macropore size or injected flow rate, in the range of simulated flow rates. These experimental and numerical results show that macropore influence on transport should be extended through the surrounding porous matrix into a zone of the same size of grains diameter for flow and into a zone depending on diffusion coefficient as well as mean residence time of the studied tracer for solute transport

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.

Carbon dioxide (CO2) sequestration in the sub-surface is a potential mitigation technique for global climate change caused by greenhouse gas emissions. In order to evaluate the feasibility of this technique, understanding the behaviour of CO2 stored in geological rock formations over a range of length- and time-scales is crucial. The work presented in this dissertation contributes to the knowledge in this field by investigating the two-phase flow and entrapment processes of CO2, as well as other relevant fluids, in porous media at the pore- and centimetre-scales using a combination of lab-based nuclear magnetic resonance (NMR) experimental techniques and lattice Boltzmann (LB) numerical simulation techniques. Pulsed field gradient (PFG) NMR techniques were used to acquire displacement distributions (propagators) of brine flow through a model porous medium (100 µm glass bead packing) before and after the capillary (residual) trapping of gas-phase CO2 in the pore space. The acquired propagators were compared quantitatively with the corresponding LB simulations. In addition, magnetic resonance imaging (MRI) techniques were used to characterise the extent of CO2 trapping in the bead pack. The acquired NMR propagators were compared to LB simulations applied to various CO2 entrapment scenarios in order to investigate the pore morphology in which CO2 becomes entrapped. Subsequently, MRI drop shape analysis techniques were used to identify a pair of analogue fluids which matched certain key physical properties (specifically interfacial tension) of the supercritical CO2/water system in order to extend the work to conditions more relevant to CO2 sequestration in the sub-surface, where CO2 is likely to be present in the supercritical phase. As before, NMR propagator measurements and MRI techniques, along with LB simulations, were used to characterise the capillary trapping of the CO2 analogue phase in glass bead packs, as well as two different types of rock core plugs – relatively homogeneous Bentheimer sandstone, and heterogeneous Portland carbonate. In addition to capillary trapping, the effect of vertical permeability heterogeneity, such as is often present in underground rock formations, was investigated for the flow of miscible (water/brine) gravity currents in model porous media (glass bead packs), using MRI techniques such as 2D spin-echo imaging and phase-shift velocity imaging. Finally, a preliminary investigation was made into the effect of particle- and pore-size distributions on the gas/liquid (air/water) interface for porous media consisting of glass bead and sand packs of different average particle size using quantitative MRI techniques.

This thesis is concerned with the new formulation of a finite-volume lattice Boltzmann equation method and its implementation on unstructured meshes. The finite-volume discretization with a cell-centered tessellation is employed. The new formulation effectively adopts a total variation diminishing concept. The formulation is analyzed for the modified partial differential equation and the apparent viscosity of the model. Further, the high-order extension of the present formulation is laid out. Parallel simulations of a variety of two-dimensional benchmark flows are carried out to validate the formulation. In Chapter 1, the important notions of the kinetic theory and the most celebrated equation in the kinetic theory, ‘the Boltzmann equation’ are given. The historical developments and the theory of a discrete form of Boltzmann equation are briefly discussed. Various off-lattice schemes are introduced. Various methodologies adopted in the past for the solution of the lattice Boltzmann equation on finite-volume discretization are reviewed. The basic objectives of this thesis are stated. In Chapter2,the basic formulations of lattice Boltzmann equation method with a rational behind different boundary condition implementations are discussed. The benchmark flows are studied for various flow phenomenon with the parallel code developed in-house. In particular, the new benchmark solution is given for the flow induced inside a rectangular, deep cavity. In Chapter 3, the need for off-lattice schemes and a general introduction to the finite-volume approach and unstructured mesh technology are given. A new mathematical formulation of the off-lattice finite-volume lattice Boltzmann equation procedure on a cell-centered, arbitrary triangular tessellation is laid out. This formulation employs the total variation diminishing procedure to treat the advection terms. The implementation of the boundary condition is given with an outline of the numerical implementation. The Chapman-Enskog (CE) expansion is performed to derive the conservation equations and an expression for the apparent viscosity from the finite-volume lattice Boltzmann equation formulation in Chapter 4. Further, the numerical investigations are performed to analyze the apparent viscosity variation with respect to the grid resolution. In Chapter 5, an extensive validation of the newly formulated finite-volume scheme is presented. The benchmark flows considered are of increasing complexity and are namely (1) Posieuille flow, (2) unsteady Couette flow, (3) lid-driven cavity flow, (4) flow past a backward step and (5) steady flow past a circular cylinder. Further, a sensitivity study to the various limiter functions has also been carried out. The main objective of Chapter6is to enhance the order of accuracy of spatio-temporal calculations in the newly presented finite-volume lattice Boltzmann equation formulation. Further, efficient implementation of the formulation for parallel processing is carried out. An appropriate decomposition of the computational domain is performed using a graph partitioning tool. The order-of-accuracy has been verified by simulating a flow past a curved surface. The extended formulation is employed to study more complex unsteady flows past circular cylinders. In Chapter 7, the main conclusions of this thesis are summarized. Possible issues to be examined for further improvements in the formulation are identified. Further, the potential applications of the present formulation are discussed.

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

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.

Computer-generated imagery is ubiquitous in today's society, appearing in advertisements, video games, and computer-animated movies among other places. Much of this imagery needs to be as realistic as possible, and animators have turned to techniques such as fluid simulation to create scenes involving substances like smoke, fire, and water. The Lattice-Boltzmann Method (LBM) is one fluid simulation technique that has gained recent popularity due to its relatively simple basic algorithm and the ease with which it can be distributed across multiple processors. Unfortunately, current LBM simulations also suffer from high memory usage and restrict free surface fluids to domains of fixed size. This thesis modifies the LBM to utilize a recursive run-length-encoded (RLE) grid data structure instead of the standard fixed array of grid cells, which reduces the amount of memory required for LBM simulations as well as allowing the domain to grow and shrink as necessary to accomodate a liquid surface. The modified LBM is implemented within the open-source 3D animation package Blender and compared to Blender's current LBM simulator using the metrics of memory usage and time required to complete a given simulation. Results show that, although the RLE-based simulator can take several times longer than the current simulator to complete a given simulation, the memory usage is significantly reduced, making an RLE-based simulation preferable in a few specific circumstances.

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