Sand transport in hilly terrain geometry is different and complex to understand compared to horizontal pipeline, due to the influence of the geometry that greatly affect multiphase flow and sand behaviour at the dip. The overall aim of this research work is to use experimental method to investigate the effects of multiphase flow behaviour on sand transport in a dip configuration. Experimental work was carried out to understand the complex dynamic mechanisms that exist during sand multiphase flow using 2” inch dip test facility with different inclination angles of ±24° and ±12° configurations. In order determine the multiphase flow parameters and sand flow regimes, extensive data were collected and analysed from continuous local measurement of instantaneous liquid hold up and sand hold up using conductivity rings. Significant observations were made during this study from which several conclusions were made. In the air–water test, three slug behaviours were observed at the dip: complete stratified flow downhill with slug initiation at dip; stratified flow (with energetic ripple) downhill with slug initiation and slug growth upward dip; and aerated slug downhill and slug growth at the dip. These behaviours are different from published work on this subject with low angle of inclination. The data analysis revealed that the two types of slug initiation mechanisms (wave growth and wave coalescence) observed are geometry specifics. The slug translational velocities (at the dip and uphill section) were used as criterion to determine the flow condition for each slug initiation mechanism at the dip. Significant observations were made during this study from which several conclusions were made. In the air–water test, three slug behaviours were observed at the dip: complete stratified flow downhill with slug initiation at dip; stratified flow (with energetic ripple) downhill with slug initiation and slug growth upward dip; and aerated slug downhill and slug growth at the dip. These behaviours are different from published work on this subject with low angle of inclination. The data analysis revealed that the two types of slug initiation mechanisms (wave growth and wave coalescence) observed are geometry specifics. The slug translational velocities (at the dip and uphill section) were used as criterion to determine the flow condition for each slug initiation mechanism at the dip. Five sand-water flow regimes (full suspension, streak, saltation, sand dune, and sand bed) were established by physical observation and data analysis. It was also observed that sand streaks were denser towards the central line of pipe bottom in the downhill pipe than that in uphill pipe. At downhill pipe section, there were sand gathering toward the central line of the pipe bottom. The characteristics of sand transportation at the dip section were found slightly different from downhill and uphill pipe for higher sand concentrations. When dense streak occurred at the downhill, the sand particles become dispersed at the dip. The minimum transport conditions (MTC) were determined at different sand concentration. The sand minimum transport condition in the dip section was found to be slightly lower than those in the downhill and uphill section. The minimum transport condition for a single phase water flow for the 24˚ dip. test section was slightly higher (with difference of about 0.1m/s) than that of the 12˚ at the downward and upward of the dip section at low sand concentration. In addition, local sand measurements using conductivity time series results at the downhill and uphill section showed the influence of sand concentration and flow condition on sand flow patterns. The air-water-sand results showed that sand deposits occurred in uphill section after sand transport at the downhill and dip sections. This was due to different flow regimes exhibited at the different pipe sections. The stratified (wavy) flow was the dominant flow in downhill pipe; therefore sand was observed transporting within the liquid film as thin streak for most of test conditions. The slug initiation at the dip section was observed to be a major factor that influences the sand behaviour. Sand particles in the slug unit (at the dip and uphill of the pipe) were observed to be entrained in the slug body once slug is initiated, thereby changing the force vector generating turbulence flow at the front of slug body. Once the sand particles entered the film zone of the slug unit, they immediately stopped moving forward due to the film velocity significantly lower than the slug body coupled with gravity effect. . Sand particles were found to be falling back while travelling with the water film at some conditions, until they were picked up by the next slug body. The results of this work provide a better understanding to the study of multiphase flow for pipeline design and most especially sand behaviour at the dip. The sand dune regime is identified distinctively using conductivity ring technique which would assist in determining the operating conditions that allow sand dune formation. The knowledge of flow condition at full suspension of sand is an important parameter to determine the erosion rate over the life span of the pipeline. Also, the quantity of sand bed and flow condition of sand settling at the dip is useful information for production chemist in order to determine the effectiveness of corrosion inhibitor at the bottom of the pipe. In conclusion, sand transport characteristics and MTC were strongly dependent on the gas-liquid flow regime and pipe geometry; and cannot be generalised on the superficial liquid and gas velocities of the transport fluid.

Membrane filtration is used in a variety of industries, including water treatment and the food industry. Membrane systems include microfiltration and reverse osmosis processes. Membranes used in reverse osmosis are nonporous or pores at 0.2-2 A. This work will focus on mechanical microfiltration. These filtration systems suffer from an accumulation of the rejected material near the membrane surface. This causes additional resistance to the flow through the membrane (flux), resulting in a decline in the performance of the system. Sparging gas bubbles into the mixture has been shown to improve performance. The flow field promotes the transport of material away from the membrane surface and into the bulk. The goal is to predict the sparging that will achieve the maximum flux. Existing flux prediction models often assume steady shear at the membrane surface but in bubbling regimes the shear stresses are unsteady. In this thesis a model is developed to calculate the flux based not solely on shear but on the behaviour and resistance of suspended particles in a gas-liquid flow field. The bubble shape and flow field is calculated using computation fluid dynamics (CFD). The flow around a bubble in gap between two parallel flat sheet membranes is investigated. The calculated bubble shape correlates well with the results seen in experiments. The bubble rise velocity with respect to gap width is shown to transition between that expected in the literature for extended flow for large gap widths and that for a two dimensional case for smaller gap widths. The transitional region however, does not behave as may be expected. The rise velocity does not monotonically decrease as the gap width is reduced. The particle concentration is found by the solution of the convection-diffusion equation, where the convection velocity terms are given by the results of the CFD calculation. The permeate flux is then calculated using a resistance model giving the enhancement due to the bubble. The model is also applied to single phase crossflow. As the shear stresses are steady in this single-phase flow regime, established membrane shear linked mass-transfer coefficient methods can be employed. Good agreement is found between the model and theory. The flux results obtained when the model is applied to the flow around the bubble show a peak in performance with respect to the gap between the membranes for a given bubble volume. The optimal flux enhancement is found to correlate well with the bubble size compared to the flow area. The results show a bubble width of around 60% of the flow width provides the best flux performance.

<div> Underlying physics in bubble coalescence is critical for understanding bubble transportation. It is one of the major mechanisms of microfluidics. Understanding the mechanism has benefits in the design, development, and optimization of microfluidics for various applications. The underlying physics in bubble coalescence is investigated numerically using the free energy-based lattice Boltzmann method by massive parametrization and classification.</div><div><br></div><div> Firstly, comprehensive GPU (Graphics Processing Unit) parallelization, convergence check, and validation are carried out to ensure the computational efficiency and physical accuracy for the numerical simulations.</div><div><br></div><div> Then, the liquid-gas system is characterized by an Ohnesorge number (Oh). Two distinct coalescence phenomena with and without oscillation, are separated by a critical Oh (~0.477)number. For the oscillation cases(Oh<0.477), the mechanism of damped oscillation in microbubble coalescence is explored in terms of the competition between driving and resisting forces. Through an analogy to the conventional damped harmonic oscillator, the saddle-point trajectory over the entire oscillation can be well predicted analytically. Without oscillation in the range of 0.50r^-n </div><div><br></div><div> After that, the liquid-gas-solid interface is taken into consideration in the liquid-gas system. Six cases based on the experiment set-ups are simulated first for validation of the computational results. Based on these, a hypothesis is established about critical factors to determine if coalescence-induced microbubble detachment (CIMD) will occur. From the eighteen experimental and computational cases, we conclude that when the radius ratio is close to 1 and the father bubble is larger, then it will lead to CIMD.</div><div><br></div><div> Lastly, the effects of initial conditions on the coalescence of two equal-sized air microbubbles (R₀) in water are investigated. In both initial scenarios, the neck bridge evolution exhibits a half power-law scaling, r/R₀=A₀(t/t_i)^1/2 after development time. The development time is caused by the significant bias between the capillary forces contributed by the meniscus curvature and the neck bridge curvature. Meanwhile, the physical mechanism behind each behavior has been explored.</div>

Tandis que l'emploi des techniques de renforcement des structures s'est largement généralisé et diversifié, les méthodes de calcul et de simulation du comportement de telles structures, par nature composites, exigent encore de nombreux développements, tant sur le plan théorique (recours aux techniques d'homogénéisation), que numérique. Ainsi, dans le domaine du génie civil, une modélisation qualifiée de multiphasique a été récemment proposée pour les ouvrages en sols renforcés par inclusions linéaires continues souples (terre armée, géotextiles, etc.) ou raides (inclusions "rigides", pieux, etc.).Ce présent travail a pour but de développer une méthode de calcul rapide et fiable, à travers cette modélisation multiphasique, pour le dimensionnement vis-à-vis de sollicitations dynamiques dans le cas où l'ouvrage est renforcé par un groupe de pieux ou d'inclusions rigides, en se restreignant au cas de l'élastodynamique, c'est à dire d'un comportement élastique linéaire des différents constituants du sol renforcé. Il consiste d'une part à analyser la propagation d"ondes sismiques au sein du massif renforcé et d'autre part à mettre en oeuvre un outil numérique basé sur la méthode des éléments finis pour déterminer les fonctions d'impédance d'un sol renforcé par un réseau régulière d'inclusions vertical. D'où les effets de l'interaction sol-inclusions, ainsi de flexion et de cisaillement des inclusions sont pris en compte. / The reinforcement of structures in becoming an increasingly used technique, although the simulation and design of such structures still require many developpements both in theory (use of homogenization techniques) and numeric. Thus, in the civil engineering domains, a qualified multiphase model has been recently proposed for soil reinforced by continuous linear inclusions flexible (reinforced earth structures, geotextiles, etc. ) or stiff (rigid inclusions, piles, etc.).The present work aims to develop a fast and reliable method of calculation, through such a multiphase model, for the design with dynamic load applied on reinforced soil by a group of piles or rigid inclusions, by restricting the elastodynamic case, that is a linear elastic behavior for both constituents. It consists firstly to analyze the propagation of seismic waves in the solid reinforced and secondly to implement a fem.-based numerical code for determining the impedance functions of a reinforced soil by a regular vertical network inclusions. In this model, the interacton soil-inclusions and also the shear and flexural effects of inclusions are both taken into account.

Adverse health effects attributable to both short- and long-term exposure to air pollution have turned the focus on different microenvironments. The interior of vehicles is of relevance as road traffic emissions and re-suspension of road dust are major sources of pollutants associated with adverse health effects. Hence, the air quality inside vehicles deserves attention regarding human health. This thesis presents a new virtual methodology, using CFD, to study the distribution of fine particulate matter, PM2.5, inside a car cabin. In the CFD model, unsteady RANS and Lagrangian particle tracking were used to simulate particles entering from the exterior. In this study, a practical measurement of interior particle concentrations was also carried out as a first attempt to validate the CFD model. The objective was to find positions inside the cabin where elevated concentrations of PM2.5 are present. The results from the CFD simulations showed that significantly higher concentrations are present at head height in the front row. Due to a discrepancy in the investigated positions in the CFD model and the practical measurement, the simulation results could not be validated. Nevertheless, the simulation results in this study have provided guidelines for future measurements of interior particle concentrations.

Numerical simulation of cell motility is one of the difficult problems in computational science. It belongs to a class of problems which involve moving interfaces between flowing or deforming media. Different numerical techniques are being developed for different application areas and in this work an attempt is made to apply two popular numerical techniques used in the field of computational multiphase flows to a cell motility problem. An unsteady cell motility problem is considered to simulate numerically based on a two-dimensional mathematical model. Two important numerical methods, the Level set method and the Front tracking methods are applied to the cell motility problem to study several cases and to verify the convergence of the solution. With the assumption of no mechanical or physical obstructions to the cell, the results of the numerical simulations show that the domain shapes converge to a circular shape as they reach the steady state condition. The final steady state velocities with which the domains move and the final steady state area to which they converge are observed to be independent of domain shapes. Moreover all shapes converge to exactly same radius of circle and move with same velocity after reaching steady state condition.

Modelling of laminar multiphase flow is extremely important in a wide range of engineering and scientific applications. The particle phases are often difficult to model, especially when particles display a range of sizes and velocities at each location in space. Lagrangian methods can be too expensive and many Eulerian methods, though often computationally more affordable, suffer from model deficiencies and mathematical artifacts that lead to non-physical results. For example, efficient Eulerian models that can accurately predict the crossing of multiple streams of non-interacting particles in laminar flow have traditionally been lacking. The predictive capabilities of modern techniques from the kinetic theory of gases to the treatment of disperse multiphase flows are investigated. In particular, several moment-methods, including a recently proposed fourteen-moment approximation to the underlying kinetic equation describing particle motion, are considered and their abilities to predict particle-stream crossing are assessed. Furthermore, a new polydisperse model has been proposed for treatment of flows that display a range of particles sizes. The proposed model is an extension of the well-known maximum-entropy ten-moment model from rarefied gas dynamics with an addition for the treatment of a range of particle diameters. This model allows for anisotropic variance of particle velocities in phase space and directly treats correlations between particle diameter and velocity. The derivation and mathematical structure, of the proposed models are presented. A fine-volume discretization solution procedure for the resulting moment equations is described and used for performing numerical experiments. Results for flow problems that are designed to demonstrate the fundamental behaviour of each model are presented. It is shown that the new models offer clear advantages in terms of accuracy as compared to traditional Eulerian models for multiphase flows.

Physics of immiscible two-phase flow in porous media is relevant for various industrial and environmental applications. Wettability defined as the relative affinity of fluids with the solid surface has a significant impact on the dynamics of immiscible displacement. Although wettability effects on the macroscopic fluid flow behaviour are well known, there is a lack of pore-scale understanding. Considering the crucial role of wettability in a diverse range of applications; this research aims to provide a pore-scale picture of interface configuration induced by variations in the wetting characteristics of porous media. Besides, this study also relates the pore-scale interfacial phenomena with the macroscopic response of fluids. High-resolution direct numerical simulations (DNS) at multiscale (single capillary and a highly heterogeneous porous media) were performed using computational fluid dynamics (CFD) approach in which the Navier-Stokes equation coupled with the volume of fluid method is solved to represent immiscible displacement. Numerical results demonstrate that at pore scale as the wettability of porous media changes from strong to intermediate wet the effects of pore geometry (that includes corner angle and orientation angle) on the interfacial dynamics also enhances. This was demonstrated by the non-monotonic behaviour of entry capillary pressure at the junction of pore, curvature reversal in the converging-diverging capillary and the co-existence of concave and convex interfaces in heterogeneous porous media with uniform contact angle distribution. In addition to simulations, theoretical argument is also presented that rationalize the underlying physics of complex, yet intriguing interfacial phenomena shown by DNS. Overall this research extends the fundamental understanding of multiphase flow in porous media and paves the way for future studies on porous media.

The effects of bubbles on a jet issued below and parallel to a free surface are experimentally studied. The jet under study is isothermal and in fresh water, with air injectors that allow variation of the inlet air volume fraction for 0% to 13%. Measurements of the jet exit conditions, water velocity, water entrainment, Reynolds stresses and surface currents have been performed using LDV, PIV and surface PIV. Air volume fraction, bubble velocity, chord length and free surface elevation and RMS have been obtained using local phase detection probes. Visualization was performed using laser-induced fluorescence. Measurements show that water entrainment decreases up to 22% with the presence of bubbles, but surface current strength increases up to 60% with 0.4 l/min of air injection. The mean free surface elevation and turbulent fluctuation significantly increase with the injection of air. The water normal Reynolds stresses are damped by the presence of bubbles in the bulk of the liquid, but very close to the free surface the effect is reversed and the normal Reynolds stresses increase slightly for the bubbly flow. Flow visualizations show that the two-phase jet is lifted with the presence of bubbles and attaches to the free surface sooner. Significant bubble coalescence is observed, leading to an increase of 20% in mean bubble size as the jet develops. The coalescence near the free surface is particularly strong, due to the time it takes the bubbles to pierce the free surface, resulting in a considerable increase in the local air volume fraction.

Drinking water distribution networks such as South East Water Ltd. (SEWL), Melbourne Water, Sydney Water, etc. are supposed to transport only dissolved matter rather than a few visible particles. However, it is almost impossible to make the drinking water free from suspended solid particles. The ability to determine the origins of these particles varies between different water supply systems, with possible sources being from catchment, treatment processes, biofilm growth within the water supply pipes, and corrosion products. Improvement of our understanding of the complex hydrodynamic behavior of suspended and/or deposited particles involved in these distribution pipe networks requires mathematical and physical models. Computational Fluid Dynamics (CFD) along with analytical turbulent model is one of the most popular mathematical techniques, which has the ability to predict the behavior of complex flows for such multiphase flow applications. This study has been completed mainly in two steps. A CFD investigation was carried out to predict the hydrodynamic behavior of turbid particle flowing through a horizontal pipe networks including loop consist of bends and straight pipes. Furthermore, an extended analytical model was re-developed for the liquid-solid system to look at the similar behavior of the solid particles flowing in a turbulent field. These two parallel studies will provide better understandings about the turbidity spikes movements in the distribution networks. A comprehensive CFD investigation was carried out for particle deposition in a horizontal pipe loop consisting of four 900 bends in a turbulent flow field. A satisfactory agreement was established with the experimental data as validation. This was a steady state multi-particle problem, which helped to understand the deposition characteristics for different particle sizes and densities at upstream and downstream sides of the bends as well as its circumference. Particle concentration was seen high at the bottom wall in the pipe flow before entering the bends, but for the downstream of bend the deposition was not seen high at the bottom as seen in upstream of bend rather inner side of the bend wall (600 skewed from bottom). The larger particles clearly showed deposition near the bottom of the wall except downstream. As expected, the smaller particles showed less tendency of deposition and this was more pronounced at higher velocity. Due to the high stream line curvature and associated centrifugal force acting on the fluid at different depths the particles became well mixed and resulted in homogeneous distribution near the bend regions. The hydrodynamic behavior of particles flowing in a turbulent unsteady state flowing through a horizontal pipe was also studied to compare with the drinking water distribution networks data. In this numerical simulation six different flow-profiles and particle-load profiles were used to compute particles deposition and re-entrainment into the systems and to identify the conditions of the deposition and suspension mechanisms. Results showed that after a certain length of pipe and period of time after downward velocity gradient, when the velocity was constants over time, the shear stress was sufficiently high enough to cause the particle deposition on and roll along the bottom wall of pipe wall and created a secondary group of particle peak (called kink). Finally, an extended analytical Turbulent Diffusion Model for liquid-solid phase was developed following an existing gas-liquid turbulence model. This turbulent diffusion model was then compared with the results of the CFD investigation making use of the same boundary conditions. The comparison established good agreement between these two models. The influence of velocity on the particle size distribution was dominant over the influence of the superficial liquid velocity, which was also explained by using the new parameter, velocity ratio. This velocity ratio was defined as the ratio between the free flight and gravitational velocity. Due to some inevitable assumptions used in the analytical model, the results showed typically less deposition as compared with the CFD investigation.