Abstract Multiphase Droplet Interactions with a Single Fiber By: Noor M. Farhan A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at Virginia Commonwealth University. Virginia Commonwealth University, 2019 Director: Hooman V. Tafreshi, Professor, Department of Mechanical and Nuclear Engineering Formulating the physics of droplet adhesion to a fiber is interesting intellectually and important industrially. A typical example of a droplet–fiber system in nature is the dew droplets on spider webs, where droplets first precipitate and grow on the fibers, but they eventually fall when they become too heavy. Obviously, quantifying the force of adhesion between a droplet and a fiber is crucial in designing fog harvesting devices or manufacturing filtration media for liquid–gas or liquid–liquid separation, among many other industrial applications. This study is aimed at developing a mathematical framework for the mechanical forces between a droplet and a fiber in terms of their physical and wetting properties. To this end, a series of experiments were conducted to detach ferrofluid droplets of varying volumes from fibers with different diameters and Young–Laplace contact angles (YLCAs) in a controlled magnetic field. The force of detachment was measured using a sensitive scale and used along with the results of numerical simulations to develop a semi-analytical expression for the force required to detach a droplet from a fiber. This universally-applicable expression allows one to predict the force detachment without the need to run an experiment or a computer simulation. This work also reports on the use of magnetic force to measure the force of detachment for nonmagnetic droplets for the first time. This is accomplished by adding a small amount of a ferrofluid to the original nonmagnetic droplet to create a compound droplet with the ferrofluid nesting inside or cloaking the nonmagnetic droplet. The ferrofluid is then used to induce a body force to the resulting compound droplet and thereby detach it from the fiber. The recorded detachment force is used directly (the case of nesting ferrofluid) or after scaling (the case of cloaking ferrofluid) to obtain the force of detachment for the original nonmagnetic droplet. The accuracy of these measurements was examined through comparison with numerical simulations as well as available experimental data in the literature. In addition, a simple method is developed to directly measure the intrinsic contact angle of a fiber (i.e., Young–Laplace Contact angle of the fiber material) with any arbitrary liquid. It is shown that the intrinsic contact angle of a fiber can be obtained by simply measuring the angle between the tangent to the fiber surface and the tangent to the droplet at the contact line, if the droplet possesses a clamshell conformation and is viewed from the longitudinal direction. The novelty of the proposed method is that its predictions are not affected by the volume of the droplet used for the experiment, the wettability of the fiber, the surface tension of the liquid, or the magnitude of the body force acting on the droplet during the experiment. Also, a liquid droplet interaction with granular coatings is simulated and the droplet apparent contact angle (ACA) and the transition from Cassie (fully dry) to Wenzel (fully wet) state as a function to the roughness wavelength have been studied. For a fixed droplet volume, two different granular coatings have been used, spherical and hemispherical bumps. It is demonstrated that the chemistry (YLCA) and geometrical parameters for the granular microtexture play an important effect on the droplet ACA and its transition from Cassie to Wenzel state.

L’étude des Réacteurs à Lit Bouillonnant (RLB) est à l’origine de ce projet de thèse. Ce type de réacteur chimique est très étudié en génie des procédés, en raison notamment de son utilisation pour l’hydrocraquage des charges lourdes. Des phénomènes complexes ont lieu dans un RLB, ce qui rend leur design et leur optimisation difficiles. Certains des mécanismes physiques prenant place dans les RLBs sont également des champs de recherche actifs en mécanique des fluides. Par conséquent, cette étude se concentre sur des mécanismes locaux participant à l’hydrodynamique des RLBs avec des catalyseurs cylindriques. Dans un premier temps, l’impact de l’anisotropie du catalyseur sur sa chute est étudié. Nous réalisons une étude expérimentale de la chute libre d’un cylindre en fluide au repos, afin de déterminer l’effet de l’anisotropie du corps sur sa dynamique. Les paramètres d’intérêt du problème sont le nombre d’Archimède du cylindre (Ar) et son rapport d’élongation (L/d). Les expériences sont menées avec deux caméras orthogonales, et des techniques de traitement d’images avancées sont développées pour parvenir à une mesure précise de la position et de l’orientation du corps en 3D. Pour (200 < Ar < 1100, 2 < L/d < 20), les cylindres adoptent différents types de trajectoire. Les deux principaux sont la chute rectiligne, durant laquelle l’axe du cylindre reste horizontal, et un mouvement de fluttering, qui est analysé en détail. D’autres types de mouvement plus complexes sont observés et discutés. De surcroît, le sillage du cylindre est analysé et caractérisé. De nombreuses particules sont présentes dans un RLB (40% de fraction massique environ). Les interactions entre ces corps multiples ont un impact fort sur le mouvement de chacun d’entre eux, mais sont très complexes. En première approximation, nous rendons compte de la présence de multiples particules en introduisant un milieu confiné. Nous étudions expérimentalement la chute d’un seul cylindre dans une cellule confinée verticale, dans laquelle le cylindre n’est libre de se mouvoir que dans deux directions. Le rapport d’élongation du cylindre (3<L/d<40) et son rapport de densité ( c / f = 1,16, 2,70, 4,50) sont les deux paramètres d’intérêt. Le nombre d’Archimède du cylindre se trouve entre les mêmes bornes qu’en milieu non confiné, et les deux modes principaux de mouvement sont aussi la chute rectiligne et le fluttering. Cependant, pour des paramètres (Ar,L/d) comparables, il existe des différences importantes dans le déplacement du cylindre comparé au cas non confiné. Nous avons également étudié l’interaction entre un cylindre en chute libre et un nuage de bullesascendantes. Cette étude a été menée expérimentalement dans la cellule confinée utilisée pour la seconde partie de la thèse. Des cylindres de plusieurs rapports de densité ( c / f = 1,16, 2,70, 4,50) and rapports d’élongation (3<L/d<20) ont été lâchés dans un nuage de bulles de fraction volumique de gaz comprise entre 2% et 5%. Plusieurs mécanismes d’interaction entre le cylindre et les bulles ont été identifiés (contact direct, interaction avec des perturbations du fluide. . .), et leur effet a été caractérisé. Nous avons effectué une étude statistique du mouvement du cylindre dans le nuage de bulles, et nous l’avons comparée aux résultats obtenus en milieu confiné et en fluide au repos. Les rapports de densité et d’élongation du cylindre jouent tous deux un rôle important dans son mouvement au sein du nuage de bulles. Des statistiques conditionnelles nous permettent d’approfondir notre analyse du contact entre le cylindre et les bulles, ainsi que du rôle de l’orientation du cylindre. Enfin, la dispersion du mouvement du cylindre dans le nuage est caractérisée. Un des principaux effets du nuage de bulles est d’accroître, via les contacts bulle cylindre, l’orientation du cylindre jusqu’à-ce qu’il soit presque vertical, ce qui a un effet très fort sur sa cinématique en comparaison avec le fluide au repos / The origin of this PhD thesis lies in the study of Ebullated Bed Reactors (EBRs). These chemical reactors are very active research topics in chemical processes, notably thanks to their usage in heavy oil processing. Many complex phenomena take place within EBRs, and make their design and optimization difficult. In fluid mechanics, a lot of physical mechanisms present in EBRs are active fields of study (three-phase flow, fluid-body interaction...). Hence, in the present work, a study of the mechanisms participating in the hydrodynamics of an EBR with cylindrical catalysts is performed. In a first part, the impact of the catalyst anisotropy on its fall is investigated. In order to gain insight on the effect of the body anisotropy on its fall dynamics, we investigate experimentally the free fall of a solid cylinder in a fluid at rest. The sensitivity to two dimensionless parameters, the Archimedes number (Ar) and the aspect ratio of the cylinder (L/d) is examined. Experiments are conducted with two orthogonal cameras, and advanced image processing techniques are developed in order to measure the position and orientation of the cylinder in 3D. Within the range of parameters studied (200 < Ar < 1100, 2 < L/d < 20), the cylinders adopt different types of falling motion. Two main types of paths are observed, the first one is a rectilinear fall of the cylinder that keeps its axis horizontal, and the second one is a fluttering oscillatory motion. Other more complex types of motion are observed and discussed. The fluttering motion of the cylinder is analyzed in details. On top of the study of the body motion, the cylinder wake is also visualized and characterized. A large number of particles are present at the same time inside an EBRs (about 40% of the mass). Interactions between multiple objects have a strong impact on the motion of each individual particle, but are very complex. In a first approximation, we take into account the presence of numerous particles by introducing a confined medium. We study experimentally the fall of a single cylinder in a confined vertical thin-gap cell, where the cylinders are free to move in only two directions. The cylinder elongation ratio (3<L/d<40) and density ratio ( c / f = 1.16, 2.70, 4.50) are the two parameters of interest. The Archimedes number of the cylinder lies within the same range as in the unconfined medium, and the two main modes of motion of the cylinder are a rectilinear motion, and a fluttering one. However, for the same parameters (Ar,L/d), the motion of the cylinder in the confined cell is strongly different in form to that in the unconfined medium. We also studied the interaction between a freely falling cylinder and a rising swarm of bubbles. This investigation was performed experimentally, in the confined cell used in the second part. Cylinders of various density ratio ( c / f = 1.16, 2.70, 4.50) and elongation ratio (3<L/d<20) are released in a bubble swarm of gas volume fraction between 2% and 5%. The cylinder motion is greatly modified by the bubble swarm. Several mechanisms of interaction between the cylinder and the bubbles are identified (direct contact, interactions with fluid perturbations...), and their effect is characterized. We perform a statistical analysis of the cylinder motion in the swarm, and compare it to results in the confined fluid at rest. The cylinder density ratio and elongation ratio both play an important role in its motion in the bubble swarm. Conditional statistics allow us to further investigate the effect of the contact between the cylinder and a bubble, and of the cylinder orientation in the swarm. Finally, the dispersion of the cylinder motion in the swarm is characterized. A major effect of the bubble swarm is to increase, through bubble-cylinder contacts, the probability of the cylinder to be in nearly vertical orientations. This drastically changes the kinematics of the cylinder as compared to its motion in the fluid at rest

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.

L’accrétion de givre sur les surfaces rencontrées en aéronautique (ailes, entrées d’airmoteur, sonde) est considérée comme un risque majeur pour la sécurité aérienne. Lesconséquences observées sont la dégradation des performances aérodynamiques pouvantconduire au décrochage, des perturbations dans les moteurs pouvant aller jusqu’`a sonextinction ou le colmatage des sondes. C’est pourquoi les avionneurs développent dessystèmes de protection thermiques contre le givre. L’eau accumulée sur les surfaces resteainsi à l’état liquide et forme un film mince. Les propriétés dynamiques (hauteur, vitesseet étalement) et thermiques (température, taux d’évaporation) du film en présence d’unécoulement d’air cisaillé permettent de prédire un éventuel regel du film d’eau en dehorsdes zones protégées ("runback ice "). Comme les essais en vols ou en soufflerie sontsouvent complexes à mettre en ouvre et onéreux, la simulation numérique est devenueun outil efficace et complémentaire pour dimensionner ces systèmes. L’objet principal decette thèse est le développement de modèles intégrés dans un outil numérique permettantde prédire le transport d’eau liquide sur une surface sous forme de film ou de ruisseletsou de gouttes. Une approche intégrale de type Saint Venant est adoptée ce qui permet dedécrire la dynamique macroscopique d’un film 3D pour des configurations et des temps decalcul raisonnables par rapport `a un calcul DNS. Une formulation augmentée du secondordre en espace pour le traitement des termes de courbure est proposée, ce qui autorisel’utilisation de maillages surfaciques non structurés généraux. Contrairement aux modèlesdisponibles dans la littérature, celui propos´e dans ce manuscrit présente l’avantage de tenircompte des phénomènes capillaires et de mouillage sans limite de validité en termesd’angle de contact statique. Une équation de conservation de l’énergie garantissant laconsistance thermodynamique des solutions calculées est dérivée du système augmentérégissant la dynamique du fluide. Une discrétisation de type Volumes Finis du systèmed’équation est proposée. Des simulations numériques valident le modèle pour des configurationsacadémiques de mouillage statiques et dynamiques. La transition d’un film continuen ruisselets est également simulée. / The ice accretion on surfaces encountered in aeronautics (wings, nacelle lips, sensors)is considered as a major risk for aviation safety. The consequences observed are the degradationof the aerodynamic performances that can lead to the aircraft stall, disturbancesin the engines that can lead to flame out or clogging of the sensors. That is why aircraftmanufacturers are developing thermal protection systems against icing. The wateraccumulated on the surfaces thus remains in the liquid state and forms a thin film. Thedynamic properties (thickness, velocity and spreading) and thermal properties (temperature,evaporation rate) of the film in the presence of a sheared air flow make it possible apotential refreezing of the water film on unprotected surfaces (”runback ice” phenomena).Since flight or wind tunnel tests are generally expensive and difficult to set up, numericalsimulation has become an effective and complementary tool to design these systems. Themain purpose of this thesis is to develop a model integrated in a numerical tool to predictthe transport of liquid water on a surface which might take the form of a film, a rivulet ora droplet. An integral approach based on a shallow water type model is adopted. It makesit possible to describe the macroscopic dynamics of a three-dimensional liquid film onrealistic configurations and within reasonable computing times compared to a full Navier-Stokes computation. An extended formulation is proposed, it corresponds to a second orderdifferential system and thus allows to use arbitrary surface meshes. Contrary to modelsavailable in the literature, the one proposed in this manuscript has the advantage of takinginto account capillary and wetting phenomena without validity limit in term of staticcontact angle. An energy conservation equation ensuring the thermodynamic consistencyof the calculated solutions is derived from the extended model governing fluid dynamics.A Finite Volume discretization of the system is proposed. Numerical simulations validatethe model for both static and dynamic academical wetting configurations. The transitionof a continuous film into rivulets is also simulated.

In single phase, steady, turbulent axisymmetric jets, the time-averaged velocity field can be characterised by the decay in centreline velocity and increased spread with increasing distance from the jet orifice. In a two-phase or ‘particle-laden’ jet, the particles will modulate the jet turbulence and exchange momentum with the gas phase. Consequently, these effects reduce both the centreline velocity decay and spreading rates with respect to the single-phase jet. Empirical exponential scaling factors were found by previous authors to describe the reduced centreline decay and spreading rates well for low Stokes numbers. In this thesis, power-law scaling factors are found to scale well a wide range of centreline velocity decay and spreading rate data published over the past 40 years, for a wide range of Stokes numbers. The power-law scaling is composed of three different regimes. For low Stokes numbers St₀ ≲20, it is found that the gas phase centreline velocity, u₀/uc collapses if plotted as a function of x/D(1 + Ø₀)⁻¹, and the velocity profile half widths r₁/ ₂ collapse if plotted as a function of x/D(1+Ø₀)⁻¹. Here, u₀ is the exit velocity, Ø₀ is the exit mass loading, x is the axial coordinate and D is the pipe diameter. For intermediate Stokes numbers, u₀/uc collapses if plotted as a function of x/D(1 + Ø₀)⁻¹ and r₁/ ₂ collapses if plotted as a function of x/D(1 + Ø₀)⁻¹/². For high Stokes numbers St₀ ≳ 200, u₀/uc collapses if plotted as a function of x/D(1 + Ø₀)⁻¹/² and the half width is approximately independent of Ø₀. In addition to the velocity of the gas phase, other aspects of particle- laden jets are found to be amenable to scaling by power-law functions. It is found that reported solid phase mass flux data scales similarly to gas phase measurements. Limited solid phase concentration and entrainment measurements reported in the literature are also found to scale by power-law functions. Whereas that limited data was obtained from the literature, measurements of the distribution of particles in particle-laden jets were conducted to further assess the validity of the scaling regimes to the solid phase. A planar light scattering technique is conducted to measure the distribution of particles in an axisymmetric jet and their subsequent scaling (or lack thereof) are reported for a variation in Ø₀, Stokes number and gas phase jet exit density. For Stokes numbers based on the pipe friction velocity St* ₀ ∼ 1, half widths of particle distributions were found to scale with x/D(1+Ø₀)⁻¹/² . The apparent centreline concentration was found to be independent of Ø₀ at this same St* ₀ . For Stokes numbers based on the pipe friction velocity St*₀ < 1, half widths are independent of Ø₀. The effect of the other parameters, i.e. Stokes number and density ratio, on centreline distributions and half widths are also investigated. Measurements of particle distributions, delivered via an annular channel, in a triangular oscillating jet (OJ) flow are also reported for a variation in momentum ratio, the ratio of OJ momentum to channel momentum and mass loading. The results of the variation in momentum ratio on particle distributions are compared with an existing precessing jet (PJ) study. It is the aim of this study to determine the experimental conditions for which the OJ nozzle is superior to the PJ nozzle. The use of an OJ nozzle is preferable at an industrial scale by virtue of its lower driving pressure compared with a PJ nozzle. It is found that particle distributions in a PJ flow spread at a greater rate with increasing momentum ratio compared with the spread of particles in an OJ flow. However, at momentum ratios approximately less than unity, the absolute spread from an OJ is greater. This also corresponds to nozzle driving pressure less than approximately 10kPA. For an increase in mass loading, the spread of particle distribution in the OJ decreases and recirculation increases. / http://proxy.library.adelaide.edu.au/login?url= http://library.adelaide.edu.au/cgi-bin/Pwebrecon.cgi?BBID=1337352 / Thesis (M.Eng.Sc.) -- University of Adelaide, School of Mechanical Engineering, 2008

In single phase, steady, turbulent axisymmetric jets, the time-averaged velocity field can be characterised by the decay in centreline velocity and increased spread with increasing distance from the jet orifice. In a two-phase or ‘particle-laden’ jet, the particles will modulate the jet turbulence and exchange momentum with the gas phase. Consequently, these effects reduce both the centreline velocity decay and spreading rates with respect to the single-phase jet. Empirical exponential scaling factors were found by previous authors to describe the reduced centreline decay and spreading rates well for low Stokes numbers. In this thesis, power-law scaling factors are found to scale well a wide range of centreline velocity decay and spreading rate data published over the past 40 years, for a wide range of Stokes numbers. The power-law scaling is composed of three different regimes. For low Stokes numbers St₀ ≲20, it is found that the gas phase centreline velocity, u₀/uc collapses if plotted as a function of x/D(1 + Ø₀)⁻¹, and the velocity profile half widths r₁/ ₂ collapse if plotted as a function of x/D(1+Ø₀)⁻¹. Here, u₀ is the exit velocity, Ø₀ is the exit mass loading, x is the axial coordinate and D is the pipe diameter. For intermediate Stokes numbers, u₀/uc collapses if plotted as a function of x/D(1 + Ø₀)⁻¹ and r₁/ ₂ collapses if plotted as a function of x/D(1 + Ø₀)⁻¹/². For high Stokes numbers St₀ ≳ 200, u₀/uc collapses if plotted as a function of x/D(1 + Ø₀)⁻¹/² and the half width is approximately independent of Ø₀. In addition to the velocity of the gas phase, other aspects of particle- laden jets are found to be amenable to scaling by power-law functions. It is found that reported solid phase mass flux data scales similarly to gas phase measurements. Limited solid phase concentration and entrainment measurements reported in the literature are also found to scale by power-law functions. Whereas that limited data was obtained from the literature, measurements of the distribution of particles in particle-laden jets were conducted to further assess the validity of the scaling regimes to the solid phase. A planar light scattering technique is conducted to measure the distribution of particles in an axisymmetric jet and their subsequent scaling (or lack thereof) are reported for a variation in Ø₀, Stokes number and gas phase jet exit density. For Stokes numbers based on the pipe friction velocity St* ₀ ∼ 1, half widths of particle distributions were found to scale with x/D(1+Ø₀)⁻¹/² . The apparent centreline concentration was found to be independent of Ø₀ at this same St* ₀ . For Stokes numbers based on the pipe friction velocity St*₀ < 1, half widths are independent of Ø₀. The effect of the other parameters, i.e. Stokes number and density ratio, on centreline distributions and half widths are also investigated. Measurements of particle distributions, delivered via an annular channel, in a triangular oscillating jet (OJ) flow are also reported for a variation in momentum ratio, the ratio of OJ momentum to channel momentum and mass loading. The results of the variation in momentum ratio on particle distributions are compared with an existing precessing jet (PJ) study. It is the aim of this study to determine the experimental conditions for which the OJ nozzle is superior to the PJ nozzle. The use of an OJ nozzle is preferable at an industrial scale by virtue of its lower driving pressure compared with a PJ nozzle. It is found that particle distributions in a PJ flow spread at a greater rate with increasing momentum ratio compared with the spread of particles in an OJ flow. However, at momentum ratios approximately less than unity, the absolute spread from an OJ is greater. This also corresponds to nozzle driving pressure less than approximately 10kPA. For an increase in mass loading, the spread of particle distribution in the OJ decreases and recirculation increases. / http://proxy.library.adelaide.edu.au/login?url= http://library.adelaide.edu.au/cgi-bin/Pwebrecon.cgi?BBID=1337352 / Thesis (M.Eng.Sc.) -- University of Adelaide, School of Mechanical Engineering, 2008

La présence de nuages influence la capacité oxydante de la troposphère, mais aussi le bilan radiatif de la planète notamment à travers les processus liés à la chimie multiphase qui restent mal connus. Le but de ce travail était d'améliorer notre compréhension des processus physico-chimiques au sein des nuages. Pour cela, nous avons développé un modèle numérique de chimie multiphase couplé à un modèle de microphysique (Berry et Reinhardt, 1974), fondé sur le modèle de chimie gazeuse de Madronich et Calvert (1990), le mécanisme réactionnel développé par Jacob (1986), une compilation exhaustive des données de la littérature et des collaborations avec des cinéticiens. Afin d'effectuer le lien nécessaire à la compréhension de la chimie multiphase entre les données de laboratoire et les données de terrain, le modèle dans sa version non couplée a été appliqué dans le cadre d'un événement nuageux échantillonné au Puy de Dôme (Voisin et coll., 2000). Ces résultats montrent en général la capacité du modèle à reproduire les comportements observés et sa capacité d'analyse de la réactivité du système chimique nuageux (Leriche et coll., 2000a). Notamment, une nouvelle voie réactionnelle d'oxydation du S(IV) par l'acide pernitrique en acide fort a été mise en évidence. Finalement, le modèle couplé chimie/microphysique a été appliqué au même événement nuageux sur la base d'un scénario académique afin de quantifier l'influence de la formation de la pluie sur le régime chimique précédemment étudié. Les résultats principaux obtenus (Leriche et coll., 2000b) montrent que la présence de nuages exerce deux effets différents sur la chimie troposphérique : un effet direct de lessivage des espèces chimiques par transfert de masse, solubilité et réactivité, et un effet indirect lié aux transferts microphysiques de l'eau nuageuse en eau précipitante et à la redistribution d'espèces réactives entre l'air interstitiel, l'eau nuageuse et l'eau précipitante.

<p>In this thesis a set of boundary conditions for multiphase flow is suggested.</p><p>Characteristic-based boundary conditions are reviewed for single-phase flow. The problem of open-boundary conditions is investigated, and to avoid drifting values, the use of control functions is proposed.</p><p>The use of control functions is also verified with a new test which assesses the quality of the boundary conditions. Particularly, P- and PI-control functions are examined. PI-controllers have the ability to specify a given variable exactly at the outlet as well as at the inlet, without causing spurious reflections which are amplified.</p><p>Averaged multiphase flow equations are reviewed, and a simplified model is established. This model is used for the boundary analysis and the computations. Due to the averaging procedure, signal speeds are reduced to the order of the flow speed. This leads to numerical challenges. For a horizontal channel flow, a splitting of the interface pressure model is suggested. This bypasses the numerical problems associated with separation by gravity, and a physical realistic model is used. In this case, the inviscid model is shown to possess complex eigenvalues, and still the characteristic boundary conditions give reasonable results.</p><p>The governing equations are solved with a Runge-Kutta scheme for the time integration. For the spatial discretisation, a finite-volume and a finite-difference method are used. Both implementations give equivalent results. In single-phase flow, the results improve significantly when a numerical filter is applied. For two-dimensional two-phase flow, the computations are unstable without a numerical filter.</p>

<p>This Master Thesis presents TMVOC simulations of a NAPL-spill (non-aqueous phase liquid) and following pumping remediation. TMVOC is a simulation program for three-phase non-isothermal multicomponent flow in saturated-unsaturated heterogeneous media. The models presented are based on an actual remediation project. The aim of the thesis was to study if the historical development of the NAPL-spill could be simulated and how long time the pumping remediation would take. A 3D-model and a radially symmetric cylindrical model were created.</p><p>A large effort of the work done was in taking the complex TMVOC model in use and modifying it for the problem at hand. Therefore, the numerical results of the simulations should be considered as preliminary and as forming basis for future studies.</p><p>The results from the spill simulation and historical pumping simulation indicated that the spill volume could be less than the estimated 1400 m³, perhaps around 700 m³, assuming a leakage time of 30 years.</p><p>The historical pumping simulation of a 700 m³ diesel spill showed good agreement with measured values for some wells, but overestimated the recovery in other wells. The overestimation could be due to the fact that the 3D-model did not take seasonal changes in the groundwater level into consideration. Also, the model did not account for any heterogeneity or compartmentalization in soil material properties that could explain the differences between the wells. </p><p>Assuming the same spill of 700 m³, future pumping was simulated. The results from these simulations indicated the remediation time to be long due to fast decreasing mobility of the NAPL phase. The NAPL flow rate to the wells was halved in a couple of years. Much of the NAPL was distributed over a large area at near residual saturation with the highest NAPL saturation found at the opposite side of the pumping wells in the model.</p><p> </p><p>Future simulation studies should address the effect of discretization as well as the effect of uncertainties in material properties e.g. conductivity, residual NAPL saturation and soil heterogeneity.</p>

Contamination of ground water by non-aqueous phase liquids (NAPLs) has received increasing attention. The most common approach to numerical modeling of NAPL movement through the unsaturated zone is the use of the finite difference or finite element methods to solve a set of partial differential equations derived from Darcy's law and the continuity equations (Abriola and Pinder, 1985; Kaluarachchi and Parker, 1989). These methods work well in many settings, but have given little insights as to why certain non-ideal flow phenomena will occur. The network modeling method, which considers flow at the pore-scale, was used in this study to better understand macroscopic flow phenomena in porous media. Pore-scale network models approximate porous medium as a connected network of pores and channels. Two and three-dimensional pore-scale network models were constructed in this study. A uniform statistical distribution was assumed to represent the random arrangement of pore and tube sizes. Both hysteristic scanning curves and intermediate fluid distribution are studied. The simulation results suggested that network models may be used to predict the characteristic curves of three-phase systems. The results also suggested that three-dimensional models are necessary to study the three-phase problems. Two-dimensional models do not provide realistic results as evidenced by their inability to provide scale-invariant representation of flow processes. The network sizes used in this study ranged from 10 x 5 (50) to 156 x 78 (12168) pores for two-dimensional and from 10 x 5 x 5 (250) to 100 x 50 x 5 (25000) pores for three-dimensional domains. The domain size of 100 x 50 x 5 pores was large enough to provide descriptions independent of the domain scale. The one important limitation of network models is the considerable computational requirements. The use of very high speed computers is essential. Except for this limitation, the network model provides an invaluable technique to study fluid transport mechanisms in the vadose zone. / Graduation date: 1996