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.

Contamination is becoming an important problem in great metropolitan areas. A large portion of the contaminants is emitted by the vehicle fleet. At European level, as well as in other economical areas, the regulation is becoming more and more restrictive. Euro regulations are the best example of this tendency. Specially important are the emissions of nitrogen oxide (NOx) and Particle Matter (PM). Two different strategies exist to reduce the emission of pollutants. One of them is trying to avoid their creation. Modifying the combustion process by means of different fuel injection laws or controlling the air regeneration are the typical methods. The second set of strategies is focused on the contaminant elimination. The NOx are reduced by means of catalysis and/or reducing atmosphere, usually created by injection of urea. The particle matter is eliminated using filters. This thesis is focused in this matter. Most of the strategies to reduce the emission of contaminants penalise fuel consumption. The particle filter is not an exception. Its installation, located in the exhaust duct, restricts the pass of the air. It increases the pressure along the whole exhaust line before the filter reducing the performance. Optimising the filter is then an important task. The efficiency of the filter has to be good enough to obey the contaminant normative. At the same time the pressure drop has to be as low as possible to optimise fuel consumption and performance. The objective of the thesis is to find the relation between the micro-structure and the macroscopic properties. With this knowledge the optimisation of the micro-structure is possible. The micro-structure of the filter mimics acicular mullite. It is created by procedural generation using random parameters. The relation between micro-structure and the macroscopic properties such as porosity and permeability are studied in detail. The flow field is solved using LabMoTer, a software developed during this thesis. The formulation is based on Lattice Botlzmann Methods, a new approach to simulate fluid dynamics. In addition, Walberla framework is used to solve the flow field too. This tool has been developed by Friedrich Alexander University of Erlangen Nürnberg. The second part of the thesis is focused on the particles immersed into the fluid. The properties of the particles are given as a function of the aerodynamic diameter. This is enough for macroscopic approximations. However, the discretization of the porous media has the same order of magnitude than the particle size. Consequently realistic geometry is necessary. Diesel particles are aggregates of spheres. A simulation tool is developed to create these aggregated using ballistic collision. The results are analysed in detail. The second step is to characterise their aerodynamic properties. Due to the small size of the particles, with the same order of magnitude than the separation between molecules of air, the fluid can not be approximated as a continuous medium. A new approach is needed. Direct Simulation Monte Carlo is the appropriate tool. A solver based on this formulation is developed. Unfortunately complex geometries could not be implemented on time. The thesis has been fruitful in several aspects. A new model based on procedural generation has been developed to create a micro-structure which mimics acicular mullite. A new CFD solver based on Lattice Boltzmann Methods, LabMoTer, has been implemented and validated. At the same time it is proposed a technique to optimized setup. Ballistic agglomeration process is studied in detail thanks to a new simulator developed ad hoc for this task. The results are studied in detail to find correlation between properties and the evolution in time. Uncertainty Quantification is used to include the Uncertainty in the models. A new Direct Simulation Monte Carlo solver has been developed and validated to calculate rarefied flow. / La contaminación se está volviendo un gran problema para las grandes áreas metropolitanas, en gran parte debido al tráfico. A nivel europeo, al igual que en otras áreas, la regulación es cada vez más restrictiva. Una buena prueba de ello es la normativa Euro de la Unión Europea. Especialmente importantes son las emisiones de óxidos de nitrógeno (NOx) y partículas (PM). La reducción de contaminantes se puede abordar desde dos estrategias distintas. La primera es la prevención. Modificar el proceso de combustión a través de las leyes de inyección o controlar la renovación de la carda son los métodos más comunes. La segunda estrategia es la eliminación. Se puede reducir los NOx mediante catálisis o atmósfera reductora y las partículas mediante la instalación de un filtro en el conducto de escape. La presente tesis se centra en el estudio de éste último. La mayoría de as estrategias para la reducción de emisiones penalizan el consumo. El filtro de partículas no es una excepción. Restringe el paso de aire. Como consecuencia la presión se incrementa a lo largo de toda la línea reduciendo las prestaciones del motor. La optimización del filtro es de vital importancia. Tiene que mantener su eficacia a la par que que se minimiza la caída de presión y con ella el consumo de combustible. El objetivo de la tesis es encontrar la relación entre la miscroestructura y las propiedades macroscópicas del filtro. Las conclusiones del estudio podrán utilizarse para optimizar la microestructura. La microestructura elegida imita los filtros de mulita acicular. Se genera por ordenador mediante generación procedimental utilizando parámetros aleatorios. Gracias a ello se puede estudiar la relación que existe entre la microestructura y las propiedades macroscópicas como la porosidad y la permeabilidad. El campo fluido se resuelve con LabMoTer, un software desarrollado en esta tesis. Está basado en Lattice Boltzmann, una nueva aproximación para simular fluidos. Además también se ha utilizado el framework Walberla desarrollado por la universidad Friedrich Alexander de Erlangen Nürnberg. La segunda parte de la tesis se centra en las partículas suspendidas en el fluido. Sus propiedades vienen dadas en función del diámetro aerodinámico. Es una buena aproximación desde un punto de vista macroscópico. Sin embargo éste no es el caso. El tamaño de la discretización requerida para calcular el medio poroso es similar al tamaño de las partículas. En consecuencia se necesita simular geometrías realistas. Las partículas Diesel son agregados de esferas. El proceso de aglomeración se ha simulado mediante colisión balística. Los resultados se han analizado con detalle. El segundo paso es la caracterización aerodinámica de los aglomerados. Debido a que el tamaño de las partículas precursoras es similar a la distancia entre moléculas el fluido no puede ser considerado un medio continuo. Se necesita una nueva aproximación. La herramienta apropiada es la Simulación Directa Monte Carlo (DSMC). Por ello se ha desarrollado un software basado en esta formulación. Desafortunadamente no ha habido tiempo suficiente como para implementar condiciones de contorno sobre geometrías complejas. La tesis ha sido fructífera en múltiples aspectos. Se ha desarrollado un modelo basado en generación procedimental capaz de crear una microestructura que aproxime mulita acicular. Se ha implementado y validado un nuevo solver CFD, LabMoTer. Además se ha planteado una técnica que optimiza la preparación del cálculo. El proceso de aglomeración se ha estudiado en detalle gracias a un nuevo simulador desarrollado ad hoc para esta tarea. Mediante el análisis estadístico de los resultados se han planteado modelos que reproducen la población de partículas y su evolución con el tiempo. Técnicas de Cuantificación de Incertidumbre se han empleado para modelar la dispersión de datos. Por último, un simulador basado / La contaminació s'està tornant un gran problema per a les grans àrees metropolitanes, en gran part degut al tràfic. A nivell europeu, a l'igual que en atres àrees, la regulació és cada volta més restrictiva. Una bona prova d'això és la normativa Euro de l'Unió Europea. Especialment importants són les emissions d'òxits de nitrogen (NOX) i partícules (PM). La reducció de contaminants se pot abordar des de dos estratègies distintes. La primera és la prevenció. Modificar el procés de combustió a través de les lleis d'inyecció o controlar la renovació de la càrrega són els mètodos més comuns. La segona estratègia és l'eliminació. Se pot reduir els NOX mediant catàlisis o atmòsfera reductora i les partícules mediant l'instalació d'un filtre en el vas d'escap. La present tesis se centra en l'estudi d'este últim. La majoria de les estratègies per a la reducció d'emissions penalisen el consum. El filtre de partícules no és una excepció. Restringix el pas d'aire. Com a conseqüència la pressió s'incrementa a lo llarc de tota la llínea reduint les prestacions del motor. L'optimisació del filtre és de vital importància. Ha de mantindre la seua eficàcia a la par que que es minimisa la caiguda de pressió i en ella el consum de combustible. L'objectiu de la tesis és trobar la relació entre la microescritura i les propietats macroscòpiques del filtre. Les conclusions de l'estudi podran utilisar-se per a optimisar la microestructura. La microestructura elegida imita els filtres de mulita acicular. Se genera per ordenador mediant generació procedimental utilisant paràmetros aleatoris. Gràcies ad això es pot estudiar la relació que existix entre la microestructura i les propietats macroscòpiques com la porositat i la permeabilitat. El camp fluït se resol en LabMoTer, un software desenrollat en esta tesis. Està basat en Lattice Boltzmann, una nova aproximació per a simular fluïts. Ademés també s'ha utilisat el framework Walberla, desentollat per l'Universitat Friedrich Alexander d'Erlangen Nürnberg. La segona part de la tesis se centra en les partícules suspeses en el fluït. Les seues propietats venen donades en funció del diàmetro aerodinàmic. És una bona aproximació des d'un punt de vista macroscòpic. No obstant este no és el cas. El tamany de la discretisació requerida per a calcular el mig porós és similar al tamany de les partícules. En conseqüència es necessita simular geometries realistes. Les partícules diésel són agregats d'esferes. El procés d'aglomeració s'ha simulat mediant colisió balística. Els resultats s'han analisat en detall. El segon pas és la caracterisació aerodinàmica dels aglomerats. Degut a que el tamany de les partícules precursores és similar a la distància entre molècules el fluït no pot ser considerat un mig continu. Se necessita una nova aproximació. La ferramenta apropiada és la Simulació Directa Monte Carlo (DSMC). Per això s'ha desenrollat un software basat en esta formulació. Malafortunadament no ha hagut temps suficient com per a implementar condicions de contorn sobre geometries complexes. La tesis ha segut fructífera en múltiples aspectes. S'ha desenrollat un model basat en generació procedimental capaç de crear una microestructura que aproxime mulita acicular. S'ha implementat i validat un nou solver CFD, LabMoTer. Ademés s'ha plantejat una tècnica que optimisa la preparació del càlcul. El procés d'aglomeració s'ha estudiat en detall gràcies a un nou simulador desenrollat ad hoc per ad esta tasca. Mediant l'anàlisis estadístic dels resultats s'han plantejat models que reproduixen la població de partícules i la seua evolució en el temps. Tècniques de Quantificació d'Incertea s'han empleat per a modelar la dispersió de senyes. Per últim, un simulador basat en DSMC s'ha desenrollat per a calcular fluïts rarificats. / García Galache, JP. (2017). Study of the flow field through the wall of a Diesel particulate filter using Lattice Boltzmann Methods [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/90413 / TESIS

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