A computational model is developed to simulate the impact dynamics of a droplet on solid substrates with the purpose of predicting the droplet spreading characteristics over time. Previous studies focused on finding relations between the impact parameters and outcome dynamics. A modified approach like the one used in this project revolves around modeling the moving contact lines at the interface in a multiphase flow environment. Focusing on research from an aircraft de-icing point of view, this study is considered a prerequisite in understanding the physics of droplet impact. The primary focus is on extending the application to incorporate super-cooled environments. Development of the model involved the use of the Volume-of-Fluid function coupled with the High-Resolution Interface Capturing scheme to model the moving contact line. The evolution of the moving contact line is modeled with contact angles as their inputs to understand the effect of the surface tension forces. Contact angle modeling is based on the Blended-Kistler method, which captures the contact angle evolution based on the surface tension and capillary number. Preliminary validation performed on the model proves its effectiveness in accurately simulating the impact behavior when compared to the literature, where the spread diameter and height agree well with experiments. The validated model is also compared to the in-house experiments performed at the Cavitation and Multiphase flow laboratory using different substrate materials. The substrates each show unique behavior - Impact on Glass results in the droplet depositing on the surface. Aluminum results in a full rebound and PET-G, results in a drop ejection. Based on inputs from the experiments - contact angles, spread diameter, and the maximum spread $beta$, show good agreement in comparison to the literature. / Master of Science / Computational model developed to simulate the impact dynamics of the droplet on solid surfaces, which predicts the evolution of the droplet over time in order to analyze the effect of the surface and properties of the fluid on the behavior of the droplet on impact. Focusing on research from an aircraft de-icing point of view, this study is considered a pre-requisite in understanding the physics of droplet impact, with potential scope in extending the simulation to applications at temperatures lower than $0^{circ}$ C. A model developed with the help of basic governing equations in fluid mechanics helps capture the effect of interactions between different physical states. The angle at which the droplet interacts with the surface (Contact Angle) and the diameter evolution (d/D) help us understand the effectiveness of the model to simulate droplet impact. Preliminary validation of the model is performed with respect to the literature where the droplet diameter evolution and the height variation match well with the experiments, which was the major criterion in determining the accuracy of the model. This model is compared to the in-house experiments performed at the Cavitation and Multiphase flow laboratory on different surfaces such as Glass, Aluminum, and Plastic (PET-G). The surfaces each show unique behavior with impact on Glass having the droplet deposit on the surface, aluminum resulting in the droplet bouncing after hitting the surface, and PET-G resulting in a small droplet being ejected from the bigger droplet which eventually deposits on the surface. Conclusions from the comparison between the experiments and the numerical simulation show how the model is effective in capturing the impact behavior on surfaces like glass in comparison to surfaces like Aluminum in this case that repels water.

This theses deals with flow of two immisible fluid in horizontal pipeline. First part teoretically describes immisible flow. What follows is experimental measurement in wich experimental track was adjusted for inlet of oil from the top. Water and corn germ oil were used as fluids. There were observed many modes of flow on the track. After that PIV was described and measured. PIV was done for measurement of values of velocity vectors. Simulation of one chosen mode was developed in the last chapter. This simulation was done in Ansys Fluent with help of VOF method. Simulation was done both in 3D and 2D pipeline. In the end these simulation were compared with experiment measurement and were critically evaluated.

On s'intéresse dans ce travail à la modélisation d'un écoulement diphasique gaz-liquide co-courant descendant dans les réacteurs à lit fixe de particules, procédé largement utilisé dans le domaine industriel. En raison de la complexité de l'écoulement, induite par les nombreuses configurations multiphasiques pouvant coexister au sein du lit, les modèles développés directement à l'échelle du réacteur sont généralement issus d'approches semi-empiriques, en considérant l'écoulement homogène. Or, il a été observé que des hétérogénéités locales, géométrique et hydrodynamique, telle qu'une mal-distribution de la phase liquide, entrainaient une diminution du taux de réaction et conduisait les modèles existants à surestimer la productivité d'un réacteur. La nécessité de prendre en compte les phénomènes microscopiques dans un modèle macroscopique à l'échelle du réacteur rend l'utilisation d'approches multi-échelles indispensable. L'écoulement étant cependant d'une nature complexe, le changement d'échelle ne peut se faire de façon directe et nécessite donc la mise en place d'outils de modélisation adaptés à une échelle intermédiaire. Dans une première étape, la méthode de simulation numérique directe ``Volume-Of-Fluid'' (VOF) est validée dans le cas d'un film ruisselant dans un tube capillaire. Cette méthode est ensuite utilisée, à l'échelle microscopique, afin de proposer et de valider des relations de fermeture pour un modèle de type ``réseau de pores'' pouvant être utilisé à une échelle intermédiaire, celle du Volume Elémentaire Représentatif. Ce changement d'échelle est tout d'abord effectué dans le cas d'un lit fixe en deux dimensions, c'est-à-dire un empilement de cylindres entre deux plaques. Cette configuration permet la mise en place d'un dispositif expérimental qui, couplé à des simulations VOF 2D à plus grande échelle, valide l'approche de type "réseau de pores" adoptée. Le modèle réseau est ensuite étendu au cas d'un lit fixe réel, c'est-à-dire en trois dimensions, dont la géométrie est obtenue par micro-tomographie. Les lois de comportement locales sont redéfinies à l'aide de simulations numériques directes à l'échelle microscopique. Les résultats provenant de simulations de type « réseaux de pore » sont ensuite confrontés, dans le cas d'une répartition homogène des phases, aux modèles 1D habituellement utilisés pour les écoulements diphasiques en lit fixe. Enfin, une campagne expérimentale est menée afin d'observer, par imagerie scanner, l'étalement d'un jet de liquide sur un empilement de grains. Une comparaison qualitative est ensuite effectuée entre les observations expérimentales et les simulations numériques réseaux dans le cas spécifique de l'étalement d'un jet de liquide / We study in this work the modelling of two-phase cocurrent downflows in fixed bed reactors, a process widely used in industry. Due to the flow complexity, i.e., the presence of different interface configurations and, therefore, different phase interactions, most models have been developed using empirical approaches, with the assumption of a homogeneous flow in the reactor. However, several studies showed that local heterogeneities, geometric and hydrodynamic, such as the liquid distribution, could have a great influence on the flow at the reactor-scale and, therefore, on the reactor performance. Consider the microscopic phenomena in a macroscopic model require the use of multi-scale approaches. However, due to the flow complexity, the upscaling cannot be done directly and requires the development of modelling tools suitable for an intermediate scale. In a first step, the direct numerical method \ Volume-Of-Fluid" (VOF) is validated in the case of a two-phase flow in a capillary tube with the presence of a thin film. Then, this method is used, at a microscopic level to propose and validate closure laws for a pore-network model which will be used to simulate the flow at the intermediate scale. This upscaling approach is first tested in a two-dimensional case,i.e., an array of cylinders between two walls. This configuration allows the set up of an experimental approach, coupled with 2D VOF simulations at the intermediate scale, in order to validate the pore-network approach. The pore-network approach is then extended to a real fixed bed, i.e. in three dimensions, whose geometry is obtained by micro-tomography. Local laws of the pore-network model are redefined using direct numerical simulations at a microscopic scale. Pore-network simulations are then compared, for a homogenous phase distribution, with 1D models typically used for two-phase flow in fixed beds. Finally, an experimental campaign was set up to observe, by imaging scanner, the spreading of a liquid jet on a fixed bed pilot. A qualitative comparison is then performed between experimental observations and pore-network simulations in the specific case of the spreading of a liquid jet

Abstract In ladle metallurgy, gas stirring and the behaviour of the slag layer are very important for alloying and the homogenization of the steel. When gas is injected through a nozzle located at the bottom of the ladle into the metal bath, the gas jet exiting the nozzle breaks up into gas bubbles. The rising bubbles break the slag layer and create an open-eye. The size of the open-eye is very important as the efficiency of the metal-slag reactions depend on the interaction between the slag and steel created during the stirring process, and information about the position and size of the open-eye is important for effective alloying practice. Moreover, the open-eye has an effect on the energy balance since it increases heat losses. In this study, experimental measurements and numerical simulations were performed to study the effect of different operating parameters on the formation of the open-eye and mixing time in a water model and industrial ladle. Experimental measurements were performed to study the effect of the gas flow rate, slag layer thickness, slag layer densities and number of porous plugs in a 1/5 scale water model and in a 150-ton steelmaking ladle. For numerical modelling, a multi-phase volume of fluid (VOF) model was used to simulate the system including the behaviour of the slag layer. The numerical simulation of the open-eye size and mixing time was found to be in good agreement with the experimental data obtained from the water model and data obtained from the industrial measurements. / Tiivistelmä Senkkametallurgiassa kaasuhuuhtelu ja kuonakerroksen käyttäytyminen ovat tärkeitä teräksen seostamisen ja homogenisoinnin näkökulmasta. Senkan pohjalla sijaitsevasta suuttimesta puhallettava kaasu hajoaa kupliksi, jotka rikkovat kuonakerroksen ja muodostavat avoimen silmäkkeen. Avoimen silmäkkeen koko on yhteydessä voimakkaampaan kuonan emulgoitumiseen, joka tehostaa metallisulan ja kuonan välisiä reaktioita. Tietoa avoimen silmäkkeen paikasta ja koosta tarvitaan myös tehokkaaseen seostuspraktiikkaan. Avoin silmäke vaikuttaa lisäksi prosessin energiataseeseen lisäten sen lämpöhäviöitä. Tässä tutkimuksessa tutkittiin kokeellisesti ja laskennallisesti erilaisten operointiparametrien vaikutusta avoimen silmäkkeen muodostumiseen vesimallissa ja terässenkassateollisessa senkassa. Kokeellisia mittauksia tehtiin kaasuhuuhtelun, kuonakerroksen paksuuden, ja suuttimien määrän vaikutuksen tutkimiseksi 1/5-mittakaavan vesimallissa ja 150 tonnin terässenkassa. Numeerisessa mallinnuksessa systeemin ja siihen lukeutuvan kuonakerroksen käyttäytymisen simuloimiseen käytettiin volume of fluid (VOF) –monifaasimenetelmää. Avoimen silmäkkeen kokoon ja sekoittumisaikaan liittyvien numeeristen simulointien havaittiin vastaavan hyvin vesimallista ja teollisista mittauksista saatua kokeellista aineistoa.

Cette thèse vise à améliorer le fonctionnement de la boîte à eau d’un véhicule automobile, organe destiné à évacuer l’eau de pluie et à alimenter l’habitacle en air propre et sec. Elle se concentre plus particulièrement sur le phénomène de tourbillon de vidange qui a pour effet de limiter le débit d’évacuation et d’augmenter le niveau d’eau dans la boîte. Cette hausse de niveau peut être source de problèmes allant de l’inconfort pour les occupants du véhicule jusqu’à une panne du boîtier de ventilation, chauffage et climatisation. Une étude expérimentale de l’écoulement est d’abord réalisée sur une géométrie simplifiée. La vélocimétrie par images de particules (PIV) et des mesures de niveau par sonde capacitive sont mises en oeuvre. A partir d’images obtenues par caméra rapide, un algorithme d’analyse d’images est développé pour mesurer le diamètre du noyau d’air au coeur du tourbillon. Ces techniques de mesure permettent d’étudier la structure du tourbillon et montrent qu’il existe une interaction entre le tourbillon et l’écoulement d’air destiné à la ventilation de l’habitacle. Cette interaction est liée à la dépression générée par l’écoulement d’air : elle modifie temporairement la structure du tourbillon et conduit à un nouvel état d’équilibre à un niveau inférieur au niveau initial. Un dispositif permettant de réduire efficacement le niveau d’eau en modifiant l’écoulement en amont du tourbillon est aussi étudié et a fait l’objet d’un dépôt de brevet. Des simulations numériques sont ensuite réalisées à l’aide du code de calcul OpenFOAM et reproduisent qualitativement l’écoulement observé expérimentalement. Bien que des différences existent sur la valeur du niveau d’eau, les dimensions et la position du noyau d’air sont correctement simulées. Enfin, les résultats expérimentaux et numériques sont comparés aux modèles de tourbillon de la littérature. / This thesis aims to improve the performance of the cowl box in automotive vehicles ; a system dedicated to draining rainwater and providing clean and dry air to the passenger compartment. This work places emphasis on the bathtub vortex phenomenon which leads to decreased drain rate and increased water level in the box. This increase can cause a variety of issues from passenger discomfort to the breakdown of the heating, ventilation and air conditioning unit. An experimental study of the flow was firstly conducted on a simplified geometry. Particle Image Velocimetry (PIV) and capacitive probe for water level measurement have been applied to study the flow structure. In addition, an image processing algorithm has been developped to measure the size of the air core with a fast camera. It enables to study the vortex structure and demonstrates that there is an interaction between the vortex and airflow intended for the ventilation of the passenger compartment. This interaction is caused by the depression created by the airflow which temporarily modifies the vortex structure and leads to a new equilibrium state at a lower water level. A device that allows the efficient reduction of the water level by modifying the upstream flow has also been studied. This device has been patented. Following the experimental aspect of the study, numerical simulations were generated using the OpenFOAM software package. Simulations qualitatively reproduce the flow which had been experimentally observed. Despite the difference on the water level, the size and position of the air core is correctly simulated. Simulation also demonstrates the existence of an airflow sucked through the drain in the air core. Lastly, the experimental and numerical results were compared to various vortex models existing in the literature.

Ce travail porte sur l’étude du procédé de projection plasma par voie liquide et en particulier sur l’interaction entre la phase liquide et l’écoulement de plasma. Différents phénomènes physiques ont lieu comme la fragmentation du jet, l’évaporation du liquide et la fusion des particules. La simulation numérique permet de mieux comprendre les phénomènes physiques et vient en appui des expériences qui sont parfois difficiles voire impossibles à réaliser. Un modèle compressible diphasique est développé afin de prendre en compte les effets compressibles du plasma. Des validations sont réalisées sur des cas test académiques et sur un jet d’argon pur. De bons accords sont obtenus entre les simulations et la théorie ou les expériences. Ensuite différents mélanges de gaz plasmagènes, classiquement utilisés en projection plasma, sont étudiés : l’argon/hydrogène et l’argon/hélium.Une analyse du caractère instationnaire et turbulent est ensuite menée et permet de voir l’importance des effets instationnaires dans ces écoulements. Puis des simulations de l’injection de la phase liquide sous forme de gouttes ou de jet continu dans l’écoulement de plasma sont réalisées. L’influence de la nature des gaz plasmagènes et du mode d’injection est démontrée. Les structures observées expérimentalement sont également obtenues parles simulations. / The study of the plasma spraying process is investigated in this work, in particular the interaction between a liquid phase and a plasma flow. Different physical phenomena occureas the jet fragmentation, the liquid vaporization and the particles fusion. The numerical simulation is used to better understand the physical phenomena and is an alternative to experimental measurements, which could be locally difficult to perform because of the measurement techniques limitations and the particular characteristics of the plasma. A compressible two-phase model is developped to take into account the compressible effects of the plasma flow. Validations are realized on academic test cases and on a argonjet. Good agreements are obtained between simulations and theory or experiments. Then different mixtures of plasmagen gas, classically used in plasma spraying, are studied :argon/hydrogen and argon/helium. An analysis of unsteady and turbulent character isled and allows seeing the importance of unsteady effects in these flows. Numerical simulations of the injection of a liquid phase into the plasma flow are realized. The influence of the plasmagen gas nature and the injection type is anlysed. Structures observed by experiments are also obtained by the simulations.

Cette thèse montre le développement d’un code de calcul pour les simulations numériques directes d’écoulements diphasiques compressibles avec évaporation. Un couplage entre les méthodes Level Set et VOF est réalisé pour le suivi d’interface. Afin de résoudre les équations de la mécanique des fluides, une méthode basée sur la pression est employée et, pour découpler la vitesse de la pression, une méthode de projection est effectuée. Cette méthode permet l’implicitation des termes liés à l’acoustique et donc de diminuer la contrainte sur le pas de temps. Le liquide et le gaz sont traités de manière compressible permettant des variations locales des masses volumiques grâce à l’utilisation d’équations d’état. L’évaporation est simulée de deux manières différentes ; une première, où un taux d’évaporation constant est employé et une seconde, où ce taux est calculé par la thermique. Parallèlement à ce sujet, une étude de la distribution des courbures dans une injection de liquide est réalisée. Cette étude permet d’étendre le concept de distribution des tailles de gouttes dans un spray et d’améliorer les informations disponibles dans le modèle ELSA. Enfin, une autre étude est effectuée sur la recherche d’un critère, basé sur les courbures à l’interface, pour estimer la qualité d’une simulation. / This PhD thesis shows the development of a numerical method for solving two-phase flows with vaporisation. A coupling between Level Set and VOF methods is realised for the interface capturing. In order to solve fluid mechanics equations, a pressure based method is employed and, to decouple velocity and pressure, a projection method is performed. This method allows the implicitation of the acoustic terms and the time step constraint reduction. Liquid and gas are considered as compressible allowing local density variations with equations of state. The vaporisation is computed in two different ways ; a first one where the vaporisation rate is constant and a second one, where this rate is calculated by thermodynamics. Along with this topic, a study on curvature distribution in a liquid injection configuration is realised. This study allows to extend the drop size distribution concept in a spray and to improve available informations on ELSA model. Finally, an other study is performed on thedevelopment of a criterion, based on interface curvatures, which estimates the quality of a simulation.

Water management is a serious issue that affects the performance and durability of PEM fuel cells. It is known, from previous experimental investigations, that surface wettability has influence on water behaviour and fuel cell performance. This finding has lead researchers to develop numerical tools for further investigation of the liquid water behaviour in gas channels. The Volume-of-Fluid (VOF) method has been used in a wide range of studies for its advantage of showing the multi-phase interface in a Computational Fluid Dynamics (CFD) simulation to understand liquid water behaviour in gas channels. In this thesis, numerical study has been carried out to examine the behaviour of liquid water in gas channels. The dynamic movement of the liquid water in the channel and the associated pressure drop, water saturation and water coverage of the GDL have been investigated. Firstly, flow diffusion into the GDL was examined to determine its effect on liquid droplet behaviour in a small section of a gas channel. Furthermore, the effects of the percentage of flow diffusion, GDL wettability, pore size, and water inlet velocity were investigated. Fluid diffusion into GDL found to have insignificant impact on liquid water behaviour so further investigations has been carried with a solid GDL surface. Secondly, gas channel geometry effect on liquid water behaviour was studied. Square, semicircle, triangle, trapezoid with a long base and trapezoid with a short base were compared to find suitable cross section geometry to carry wall wettability investigations. Among the examined geometries, the square cross section showed reasonable results for both scenarios of geometry design, fixed Reynolds number and fixed GDL interface. The effect of wall wettability was assessed by comparing nine different wall/GDL wettability combinations for straight and bend channels. Wall wettability found to have an impact on liquid water behaviour but not as much as GDL wettability. It affects liquid water saturation in the channel by a great deal by accumulating water in the channel edges affecting water behaviour. This was also proven in the last test case of a long channel where water accumulation was investigated by running the calculation until the percentage of water saturation is stabilized. It is also concluded that changing wall wettability from hydrophobic to hydrophilic doubles the percentage of channel occupied by liquid water and increases the time to reach steady state.

In 1999, California Air Resources Board (CARB) implemented a regulation that required all gasoline cars sold in California be fitted with an Onboard Refueling Vapor Recovery System (ORVR). The ORVR system is designed to prevent Volatile Organic Compounds (VOCs) from escaping into the atmosphere during refuelling by storing the gas vapours in a carbon canister. Due to the complex nature of the fuel system, making design changes could have large implications on the ORVR performance of the vehicle. It is therefore desirable to develop a CFD model that can predict the effects of design changes, thereby reducing the need to perform physical tests on each design iteration. This master thesis project was performed at the Fuel Systems department at Volvo Cars in order to help reduce project lead times and product development costs by incorporating CFD as a part of the fuel system development cycle. The CFD results obtained were validated through experimental tests that were also performed as part of this project. In this master thesis project, a CFD model was developed to simulate the refuelling of gasoline for a California specification Volvo XC90 with an OPW-11B pump pistol. The model was set up in STAR-CCM+ using the Eulerian Volume of Fluid model for multiphase flow, the RANS realizable k-epsilon turbulence model and the two layer all y+ wall treatment. The effects of the carbon canister were modelled as a porous baffle interface in the simulations where viscous and inertial resistances of the porous media were adjusted to obtain a desired pressure drop across the canister. This method proved to be a suitable simplification for this study. The effects of evaporation as well as a chemical adsorption model for the carbon canister have been excluded from the project due to time limitations. It was found that the CFD simulations were in good agreement with the experimental results, especially with respect to capturing the overall behaviour of the fuel system during refuelling. It was found that resolving the flow spatially (and temporally) in the filler pipe was a crucial part in ensuring solver stability. A pressure difference between experiment and simulation was also observed as a consequence of excluding evaporation from the CFD model. After the CFD model had been verified and validated, changes to different parts of the fuel system were investigated to observe their effects on ORVR performance. These included changing the recirculation line diameter, changing the carbon canister properties and changing the angle of how the pump pistol was inserted into the capless unit. It was found that the recirculation line diameter is a very sensitive design parameter and increasing the diameter would result in fuel vapour leaking back out into the atmosphere. Similarly, increasing the back pressure by swapping to a different carbon canister would result in the leakage of fuel vapour. On the other hand, insignificant changes in system behaviour were observed when the fuel pistol angle was changed.

In 1999, California Air Resources Board (CARB) implemented a regulation that required all gasoline cars sold in California be fitted with an Onboard Refueling Vapor Recovery System (ORVR). The ORVR system is designed to prevent Volatile Organic Compounds (VOCs) from escaping into the atmosphere during refuelling by storing the gas vapours in a carbon canister. Due to the complex nature of the fuel system, making design changes could have large implications on the ORVR performance of the vehicle. It is therefore desirable to develop a CFD model that can predict the effects of design changes, thereby reducing the need to perform physical tests on each design iteration. This master thesis project was performed at the Fuel Systems department at Volvo Cars in order to help reduce project lead times and product development costs by incorporating CFD as a part of the fuel system development cycle. The CFD results obtained were validated through experimental tests that were also performed as part of this project. In this master thesis project, a CFD model was developed to simulate the refuelling of gasoline for a California specification Volvo XC90 with an OPW-11B pump pistol. The model was set up in STAR-CCM+ using the Eulerian Volume of Fluid model for multiphase flow, the RANS realizable k − ε turbulence model and the two layer all y + wall treatment. The effects of the carbon canister were modelled as a porous baffle interface in the simulations where viscous and inertial resistances of the porous media were adjusted to obtain a desired pressure drop across the canister. This method proved to be a suitable simplification for this study. The effects of evaporation as well as a chemical adsorption model for the carbon canister have been excluded from the project due to time limitations. It was found that the CFD simulations were in good agreement with the experimental results, especially with respect to capturing the overall behaviour of the fuel system during refuelling. It was found that resolving the flow spatially (and temporally) in the filler pipe was a crucial part in ensuring solver stability. A pressure difference between experiment and simulation was also observed as a consequence of excluding evaporation from the CFD model. After the CFD model had been verified and validated, changes to different parts of the fuel system were investigated to observe their effects on ORVR performance. These included changing the recirculation line diameter, changing the carbon canister properties and changing the angle of how the pump pistol was inserted into the capless unit. It was found that the recirculation line diameter is a very sensitive design parameter and increasing the diameter would result in fuel vapour leaking back out into the atmosphere. Similarly, increasing the back pressure by swapping to a different carbon canister would result in the leakage of fuel vapour. On the other hand, insignificant changes in system behaviour were observed when the fuel pistol angle was changed.In 1999, California Air Resources Board (CARB) implemented a regulation that required all gasoline cars sold in California be fitted with an Onboard Refueling Vapor Recovery System (ORVR). The ORVR system is designed to prevent Volatile Organic Compounds (VOCs) from escaping into the atmosphere during refuelling by storing the gas vapours in a carbon canister. Due to the complex nature of the fuel system, making design changes could have large implications on the ORVR performance of the vehicle. It is therefore desirable to develop a CFD model that can predict the effects of design changes, thereby reducing the need to perform physical tests on each design iteration. This master thesis project was performed at the Fuel Systems department at Volvo Cars in order to help reduce project lead times and product development costs by incorporating CFD as a part of the fuel system development cycle. The CFD results obtained were validated through experimental tests that were also performed as part of this project. In this master thesis project, a CFD model was developed to simulate the refuelling of gasoline for a California specification Volvo XC90 with an OPW-11B pump pistol. The model was set up in STAR-CCM+ using the Eulerian Volume of Fluid model for multiphase flow, the RANS realizable k − ε turbulence model and the two layer all y + wall treatment. The effects of the carbon canister were modelled as a porous baffle interface in the simulations where viscous and inertial resistances of the porous media were adjusted to obtain a desired pressure drop across the canister. This method proved to be a suitable simplification for this study. The effects of evaporation as well as a chemical adsorption model for the carbon canister have been excluded from the project due to time limitations. It was found that the CFD simulations were in good agreement with the experimental results, especially with respect to capturing the overall behaviour of the fuel system during refuelling. It was found that resolving the flow spatially (and temporally) in the filler pipe was a crucial part in ensuring solver stability. A pressure difference between experiment and simulation was also observed as a consequence of excluding evaporation from the CFD model. After the CFD model had been verified and validated, changes to different parts of the fuel system were investigated to observe their effects on ORVR performance. These included changing the recirculation line diameter, changing the carbon canister properties and changing the angle of how the pump pistol was inserted into the capless unit. It was found that the recirculation line diameter is a very sensitive design parameter and increasing the diameter would result in fuel vapour leaking back out into the atmosphere. Similarly, increasing the back pressure by swapping to a different carbon canister would result in the leakage of fuel vapour. On the other hand, insignificant changes in system behaviour were observed when the fuel pistol angle was changed.