Relevant Thesis-Based Degree Programs
Affiliations to Research Centres, Institutes & Clusters
- Development of massively parallel codes for granular and particle-laden flow simulation
- Multi-scale modeling of particle-laden flows
- Shape effect in granular and particle-laden flows
- Heat and mass transfer in fluidized beds
- Blood flow and application to targeted drug delivery
Complete these steps before you reach out to a faculty member!
- Familiarize yourself with program requirements. You want to learn as much as possible from the information available to you before you reach out to a faculty member. Be sure to visit the graduate degree program listing and program-specific websites.
- Check whether the program requires you to seek commitment from a supervisor prior to submitting an application. For some programs this is an essential step while others match successful applicants with faculty members within the first year of study. This is either indicated in the program profile under "Admission Information & Requirements" - "Prepare Application" - "Supervision" or on the program website.
- Identify specific faculty members who are conducting research in your specific area of interest.
- Establish that your research interests align with the faculty member’s research interests.
- Read up on the faculty members in the program and the research being conducted in the department.
- Familiarize yourself with their work, read their recent publications and past theses/dissertations that they supervised. Be certain that their research is indeed what you are hoping to study.
- Compose an error-free and grammatically correct email addressed to your specifically targeted faculty member, and remember to use their correct titles.
- Do not send non-specific, mass emails to everyone in the department hoping for a match.
- Address the faculty members by name. Your contact should be genuine rather than generic.
- Include a brief outline of your academic background, why you are interested in working with the faculty member, and what experience you could bring to the department. The supervision enquiry form guides you with targeted questions. Ensure to craft compelling answers to these questions.
- Highlight your achievements and why you are a top student. Faculty members receive dozens of requests from prospective students and you may have less than 30 seconds to pique someone’s interest.
- Demonstrate that you are familiar with their research:
- Convey the specific ways you are a good fit for the program.
- Convey the specific ways the program/lab/faculty member is a good fit for the research you are interested in/already conducting.
- Be enthusiastic, but don’t overdo it.
G+PS regularly provides virtual sessions that focus on admission requirements and procedures and tips how to improve your application.
ADVICE AND INSIGHTS FROM UBC FACULTY ON REACHING OUT TO SUPERVISORS
These videos contain some general advice from faculty across UBC on finding and reaching out to a potential thesis supervisor.
Graduate Student Supervision
Doctoral Student Supervision
Dissertations completed in 2010 or later are listed below. Please note that there is a 6-12 month delay to add the latest dissertations.
Biological fluids such as blood accomplish many vital tasks in the human body, including carrying oxygen and nutrients to tissues, regulating internal temperature and pH, or transporting white blood cells to infected areas. A better understanding of these fluids can provide insight into many pathologies such as the formation of aneurysms and the effect of sickle cell disease on the flow of red blood cells, as well as help design efficient diagnosis tools on microfluidic devices. Such fluids are composed of a continuous viscous phase and suspended bodies, including rigid particles and deformable membranes enclosing an inner fluid, referred to as capsules. In this thesis, we develop numerical tools aiming to simulate cell-resolved biological fluids such as blood. In a first part, we focus on the dispersed solid phase, a field known as granular mechanics. In this context, we implement a contact force able to accurately model static assemblies of granular media. After extensive validation, we use this contact model in a purely granular setting to study avalanches of entangled particles. Our numerical results are compared to experiments and show very good qualitative and quantitative agreement. Moreover, we present a variety of novel avalanching behaviors, as well as an intermittent regime in which reproducibility is lost. After analyzing the microstructure of granular assemblies in this regime, we conclude that it likely arises from mesoscale clusters of particles. In a second part, we concentrate on flowing biological capsules. We develop an adaptive front-tracking method which enables simulations of capsules in very large geometries for a wide range of Reynolds number. We validate our solver extensively and we show excellent qualitative and quantitative agreement with the literature. We then study the dynamics of capsules flowing through a sharp corner, a commonly encountered geometry in microfluidic devices. We analyze the trajectory, normalized velocity and area variations of the capsules and we show that in our case of strong confinement, the capsules interact weakly unless they are located very close to each other. Finally, we present and implement a fully Eulerian alternative method to simulate flowing capsules, and we highlight its advantages and limits.
Polydisperse particle-laden flows are ubiquitous in both nature and industry across diverse disciplines. An in-depth understanding of the intricate interactions between the dispersed solid phase and carrier fluid phase is essential in terms of designing and optimizing the key components of various industrial processes. Due to the large separation of temporal and spatial scales across different phases in real-world problems, previous studies encountered difficulty establishing a proper two-way coupling with both satisfactory accuracy and efficiency: either the various interactions are coarse-grained and oversimplified by the average drag closures in large-scale models, or the computational costs may soar beyond affordability while fully resolving the boundary of smallest particles in industrial-level simulations. To address this challenge, we develop two deterministic models based on the particle-resolved simulation data to estimate the force and torque fluctuations within polydisperse sphere assembly, which provides a channel for bottom-up knowledge transfer to evaluate the macroscale interphase interactions via the microscale hydrodynamic behaviors. Accordingly, the primary contributions of this dissertation feature the following aspects: First, we develop a highly scalable Immersed Boundary-Lattice Boltzmann method (IB-LBM) that is implemented on the adaptive quadtree/octree grids to model the complex solid-fluid interactions in different contexts such as flows laden with fixed or moving particles. Second, we perform a series of particle-resolved simulations of flow past random arrays of moderately to strongly bidisperse and polydisperse spheres which seemed computationally impossible on uniform grids in the past. We explore the statistical distributions of hydrodynamic forces and torques exerted on individual spheres over a wide range of simulation parameters. Finally, we extend two data-driven methods to predict such force and torque distributions, namely a Microstructure-informed Probability-driven Point-particle (MPP) model and a Physics-Informed Neural Network (PINN). We demonstrate their applicability from different perspectives including prediction and generalization performance, model complexity (i.e., computational efficiency), and interpretability in the form of binary and trinary interactions. We highlight the potential of our PINN model, which appropriately balances accuracy and efficiency, to substitute the conventional average drag closures for solid-fluid coupling in Eulerian-Lagrangian simulations.
Particle-laden flows where a dispersion of a solid phase is carried by a fluid phase are at the core of numerous industrial and natural processes, such as fluvial sediment transport and fluidized-bed reactors. The dynamics of each phase is intimately coupled with that of the other phase, leading to the emergence of complex, nontrivial interactions that can span wide ranges of spatial and temporal scales. The focus of this thesis is two-fold; namely, analysis of particle shape effects, and modeling hydrodynamic forces and torques in particle-laden flows. To this end, direct numerical simulations are performed for the generation of high-fidelity data, based on which all analyses of this thesis are carried out. In the first part, we scrutinize the dynamics of an isolated polyhedron, i.e. a cube, in highly inertial regimes and various density ratios. Robust helical motions and wake patterns are found for Reynolds numbers at which a sphere moves rectilinearly. An isolated cube exhibits remarkably larger rotational and lateral motions compared to a sphere, by which the effective drag on the particle is greatly affected. We then extend the analysis to inertial suspensions of cubes, where detailed comparisons are made with their counterpart sphere suspensions for various solid volume fractions. While strong clustering occurs in sphere suspensions, cube suspensions are found to be remarkably more homogeneous, as evident from their microstructure and momentum transfer properties. As demonstrated by their intensive transverse velocity fluctuations, cubes are more likely to break up and escape clusters, thus resisting local accumulation and making suspensions better mixed. In the second part, we develop a novel probability-driven point-particle model for the prediction of hydrodynamic forces and torques based on local microstructure in dense arrays of spheres. Following probabilistic arguments, necessary statistical information is extracted from particle-resolved simulations to construct force/torque-conditioned probability distribution maps, which are in turn used as basis functions for a regressive-type model. We subsequently show that our model is capable of predicting a substantial part of the observed force and torque variations, and is thus conceived to be highly promising for accurate interphase coupling in Euler-Lagrange simulations.
A theoretical and numerical study of yield-stress fluid creeping flow about a particle is presented motivated by theoretical aspects and industrial applications. Yield stress fluids can hold rigid particles statically buoyant if the yield stress is large enough. In addressing sedimentation of rigid particles in viscoplastic fluids, we should know this critical `yield number' beyond which there is no motion. As we get close to this limit, the role of viscosity becomes negligible in comparison to the plastic contribution in the leading order, since we are approaching the zero-shear-rate limit. Admissible stress fields in this limit can be found by using the characteristics of the governing equations of perfect plasticity (i.e., the sliplines). This approach yields a lower bound of the critical plastic drag force or equivalently the critical yield number. Admissible velocity fields also can be postulated to calculate the upper bound. This analysis methodology is examined for different families of particle shapes. Numerical experiments of either resistance or mobility problems in a viscoplastic fluid validate the predictions of slipline theory and reveal interesting aspects of the flow in the yield limit. For instance, the critical limit is not unique and here we show that for the same critical limit we may have different shaped particles that are cloaked inside the same unyielded envelope. The critical limit (or critical plastic drag coefficient) is related to the unyielded envelope rather than the particle shape. We show how to calculate the unyielded envelope directly. Here we also address the case of having multiple particles, which introduces interesting new phenomena. Firstly, plug regions can appear between the particles and connect them together, depending on the proximity and yield number. This can change the yielding behaviour since the combination forms a larger (and heavier) "particle". Moreover, small particles (that cannot move alone) can be pulled/pushed by larger particles or assembly of particles. Increasing the number of particles leads to interesting chain dynamics, including breaking and reforming.
Master's Student Supervision
Theses completed in 2010 or later are listed below. Please note that there is a 6-12 month delay to add the latest theses.
In the case of an oil spill in a marine environment, an important aspect of an early response is to confine the oil spill and prevent it from getting dispersed. At the laboratory level, chemists are capable of designing new solidifiers known as "gellants", that when applied by spraying or injection in the marine oil-spill, are capable of making the oil regroup and create a gel-like material (gelled-oil). Even though gellants make the removing of oil-spill from a marine environment easier, it is not clear whether gelled-oil is more stable in terms of break-up and how gellants affect the dispersion rate of the oil-spill. Gellants create a gelled-oil emulsion in the marine environment, which can be modelled as a dispersion of non-Newtonian droplets (gelled-oil) in a Newtonian matrix (sea water). In this study we develop a better understanding of such a liquid/liquid system subjected to imposed shear. Both the rheological nature of the oil and the oil/water surface tension are the key parameters. We address this problem from a computational viewpoint using an advanced open source academic code to perform parallel simulations. We analyze two main problems: 1) the deformation of a single droplet in a simple shear flow, and 2) the rheological behaviour of an emulsion in a simple shear flow. We show that, since applying gellants increase the surface tension and viscosity of the dispersed oil phase, the oil-water emulsion is more stable (i.e. the dispersed oil droplets deform less) and the dispersion rate of the oil in the marine environment is reduced. Moreover, we obtain results showing that the elasticity of the dispersed gelled-oil phase has a non-monotonic impact on the flow features and a very limited influence on the general behaviour of the emulsion. The conducted analysis in this project and its outcome can help to provide recommendations on how these gelled materials behave.