In this dissertation, Immersed Boundary Method (IBM) is evaluated in ribbed duct geometries to show the potential of simulating complex geometry with a simple structured grid. IBM is first investigated in well-accepted benchmark cases: channel flow and pipe flow with circular cross-section. IBM captures all the flow features with very good accuracy in these two cases. Then a two side ribbed duct geometry is test using IBM at Reynolds number of 20,000 under fully developed assumption. The IBM results agrees well with body conforming grid predictions. A one side ribbed duct geometry is also tested at a bulk Reynolds number of 1.5⨉10⁴. Three cases have been examined for this geometry: a stationary case; a case of positive rotation at a rotation number (Ro=ΩDₕ/U) of 0.3 (destabilizing); and a case of negative rotation at Ro= -0.3 (stabilizing). Time averaged mean, turbulent quantities are presented, together with heat transfer. The overall good agreement between IBM, BCG and experimental results suggests that IBM is a promising method to apply to complex blade geometries. Due to the disadvantage of IBM that it requires large amount of cells to resolve the boundary near the immersed surface, wall modeled LES (WMLES) is evaluated in the final part of this thesis. WMLES is used for simulating turbulent flow in a developing staggered ribbed U-bend duct. Three cases have been tested at a bulk Reynolds number of 10⁵: a stationary case; a positive rotation case at a rotation number Ro=0.2; and a negative rotation case at Ro=-0.2. Coriolis force effects are included in the calculation to evaluate the wall model under the influence of these effects which are known to affect shear layer turbulence production on the leading and trailing sides of the duct. Wall model LES prediction shows good agreement with experimental data. / Master of Science

Most realistic fluid flow problems are characterised by high Reynolds numbers and complex boundaries. Over the last ten years, immersed boundary methods (IBM) that are able to cope with realistic geometries have been applied to Lattice- Boltzmann methods (LBM). These methods, however, have normally been applied to low Reynolds number problems. In the present work, an iterative direct forcing IBM has been successfully coupled with a multi-domain cascaded LBM in order to investigate viscous flows around rigid, moving and wilfully deformed boundaries at a wide range of Reynolds numbers. The iterative force-correction immersed boundary method of (Zhang et al., 2016) has been selected due to the improved accuracy of the computation, while the cascaded LB formulation is used due to its superior stability at high Reynolds numbers. The coupling is shown to improve both the stability and numerical accuracy of the solution. The resulting solver has been applied to viscous flow (up to a Reynolds number of 100000) passed a NACA-0012 airfoil at a 10 degree angle of attack. Good agreement with results obtained using a body-fitted Navier-Stokes solver has been obtained. At moving or deformable boundary applications, emphasis should be given on the influence of the internal mass on the computation of the aerodynamic forces, focusing on deforming boundary motions where the rigid body approximation is no longer valid. Both the rigid body and the internal Lagrangian points approximations are examined. The resulting solver has been applied to viscous flows around an in-line oscillating cylinder, a pitching foil, a plunging SD7003 airfoil and a plunging and flapping NACA-0014 airfoil. Good agreement with experimental results and other numerical schemes has been obtained. It is shown that the internal Lagrangian points approximation accurately captures the internal mass effects in linear and angular motions, as well as in deforming motions, at Reynolds numbers up to 4 • 104. Finally, an expanded higher-order immersed boundary method which addresses two major drawbacks of the conventional IBM will be presented. First, an expanded velocity profile scheme has been developed, in order to compensate for the discontinuities caused by the gradient of the velocity across the boundary. Second, a numerical method derived from the Navier-Stokes equations in order to correct the pressure distribution across the boundary has been examined. The resulting hybrid solver has been applied to viscous flows around stationary and oscillating cylinders and examined the hovering flight of elliptical wings at low Reynolds numbers. It is shown that the proposed scheme smoothly expands the velocity profile across the boundary and increases the accuracy of the immersed boundary method. In addition, the pressure correction algorithm correctly expands the pressure profile across the boundary leading to very accurate pressure coefficient values along the boundary surface. The proposed numerical schemes are shown to be very efficient in terms of computational cost. The majority of the presented results are obtained within a few hours of CPU time on a 2.8 GHz Intel Core i7 MacBook Pro computer with a 16GB memory.

L'approche la plus courante en Mécanique des Fluides Numérique pour réaliser les simulations d'écoulement cardiovasculaire consiste à utiliser des méthodes numériques Body-fitted. Ces méthodes ont permis d'obtenir des simulations d'écoulement sanguin dans les artères qui sont précises et utiles. Toutefois, la génération du maillage body-fitted est une tâche qui demande beaucoup de temps et d'expertise à l'utilisateur.Les méthodes de Frontières Immergées sont des méthodes numériques alternatives qui ont l'avantage d'être plus simples d'emploi car elles ne requièrent aucune tâche de maillage de la part de l'utilisateur. Le travail présenté ici vise à évaluer le potentiel d'un méthode de Frontières Immergées à réaliser des simulations d'écoulement cardiovasculaire.Ce travail s'attache, dans un premier temps, à décrire les capacités de cette méthode numérique à rendre compte de l'imperméabilité et de la mobilité des parois sur des cas relativement simples mais représentatifs d'écoulements cardiovasculaires. Ensuite, des applications de la méthode à des cas d'écoulement cardiovasculaire plus complexes sont montrées. Il s'agira d'abord d'une simulation de l'écoulement dans un modèle rigide d'artère aorte. Puis, la simulation d'un écoulement à l'intérieur d'un ventricule cardiaque à paroi mobile sera montrée. / The most common approach in Computational Fluid Dynamics(CFD) for simulating blood flow into vessel is to make use of a body-fitted me-thod. This approach has lead to accurate and useful simulations of blood flowinto arteries. However, generation of the body-fitted grid is time consuming andrequires from the user an engineering knowledge.The Immersed Boundary Method has emerged as an alternate method whichdoes not require from the user any grid generation task. Simulations are done on astructured Cartesian grid which can be automatically generated. Here we addressthe question of the capability of an Immersed Boundary Method to cope withcardiovascular flow simulations.In particular, we assess the impermeable and moving properties of the wallwhen using the Immersed Boundary Method on simple but relevant vascular flowcases. Then, we show more complex and realistic cardiovascular flow simulations.The first application consists of blood flow simulation inside an aorta cross model.Then, the simulation of blood flow inside a cardiac ventricle with moving wall isshown.

We present a novel dimensionally split Cartesian cut cell method to compute inviscid, viscous and turbulent flows around rigid geometries. On a cut cell mesh, the existence of arbitrarily small boundary cells severely restricts the stable time step for an explicit numerical scheme. We solve this `small cell problem' when computing solutions for hyperbolic conservation laws by combining wave speed and geometric information to develop a novel stabilised cut cell flux. The convergence and stability of the developed technique are proved for the one-dimensional linear advection equation, while its multi-dimensional numerical performance is investigated through the computation of solutions to a number of test problems for the linear advection and Euler equations. This work was recently published in the Journal of Computational Physics (Gokhale et al., 2018). Subsequently, we develop the method further to be able to compute solutions for the compressible Navier-Stokes equations. The method is globally second order accurate in the L1 norm, fully conservative, and allows the use of time steps determined by the regular grid spacing. We provide a full description of the three-dimensional implementation of the method and evaluate its numerical performance by computing solutions to a wide range of test problems ranging from the nearly incompressible to the highly compressible flow regimes. This work was recently published in the Journal of Computational Physics (Gokhale et al., 2018). It is the first presentation of a dimensionally split cut cell method for the compressible Navier-Stokes equations in the literature. Finally, we also present an extension of the cut cell method to solve high Reynolds number turbulent automotive flows using a wall-modelled Large Eddy Simulation (WMLES) approach. A full description is provided of the coupling between the (implicit) LES solution and an equilibrium wall function on the cut cell mesh. The combined methodology is used to compute results for the turbulent flow over a square cylinder, and for flow over the SAE Notchback and DrivAer reference automotive geometries. We intend to publish the promising results as part of a future publication, which would be the first assessment of a WMLES Cartesian cut cell approach for computing automotive flows to be presented in the literature.

The fluid-structure interactions of wall-mounted slender structures, such as cilia, filaments, flaps, and flags, play an important role in a broad range of physical processes: from the coherent waving motion of vegetation, to the passive flow control capability of hair-like surface coatings. While these systems are ubiquitous, their coupled nonlinear response exhibits a wide variety of behaviours that is yet to be fully understood, especially when multiple structures are considered. The purpose of this work is to investigate, via numerical simulation, the fluid-structure interactions of arrays of slender structures over a range of input conditions. A direct modelling approach, whereby the individual structures and their dynamics are fully resolved, is realised via a lattice Boltzmann-immersed boundary model, which is coupled to two different structural solvers: an Euler-Bernoulli beam model, and a finite element model. Results are presented for three selected test cases - which build in scale from a single flap in a periodic array, to a small finite array of flaps, and finally to a large finite array - and the key behaviour modes are characterised and quantified. Results show a broad range of behaviours, which depend on the flow conditions and structural properties. In particular, the emergence of coherent waving motions are shown to be closely related to the natural frequency of the array. Furthermore, this behaviour is associated with a lock-in between the natural frequency of the array and the predicted frequency of the fluid instabilities. The original contributions of this work are: the development and application of a numerical tool for direct modelling of large arrays of slender structures; the characterisation of the behaviour of slender structures over a range of input conditions; and the exposition of key behaviour modes of slender structures and their relation to input conditions.

Aerothermodynamics is an integral component in the design and implementation of hypersonic transport systems. Accurate estimates of the aerodynamic forces and heat transfer rates are critical in trajectory analysis and for payload weight considerations. The present work seeks to investigate the ability of an unstructured Cartesian grid framework in modeling hypersonic viscous flows. The effectiveness of modeling viscous phenomena in hypersonic flows using the immersed boundary ghost cell methodology of this solver is analyzed. The capacity of this framework to predict the surface physics in a hypersonic non-reacting environment is investigated. High velocity argon gas flows past a 2-D cylinder are simulated for a set of freestream conditions (Reynolds numbers), and impact of the grid cell sizes on the quality of the solution is evaluated. Additionally, the formulation is verified over a series of hypersonic Mach numbers for the flow past a hemisphere, and compared to experimental results and empirical estimates. Next, a test case that involves flow separation and the interaction between a hypersonic shock wave and a boundary layer, and a separation bubble is investigated using various adaptive mesh refinement strategies. The immersed boundary ghost cell approach is tested with two temperature clipping strategies, and their impact on the overall solution accuracy and smoothness of the surface property predictions are compared. Finally, species diffusion terms in the conservation equations, and collision cross-section based transport coefficients are installed, and hypersonic flows in thermochemical nonequilibrium environments are studied, and comparisons of the off-surface flow properties and the surface physics predictions are evaluated. First, a 2-D cylinder in a hypersonic reacting air flow is tested with an adiabatic wall boundary condition. Next, the same geometry is tested to evaluate the viscous chemistry prediction capability of the solver with an isothermal wall boundary condition, and to identify the strengths and weaknesses of the immersed boundary ghost cell methodology in computing convective heating rates in such an environment.

Plunger pumps are needed for heavy duty sludge pumping at wastewater treatment facilities. America's leading pump manufacturer Wastecorp Inc. brought their plunger pump problem to us in late 2009. It was found that when the ow rate reaches a critical value, the plunger pump starts to generate a clicking noise. A one-dimensional model was built for studying the ow of a typical plunger pump operation. The velocities and pressures are calculated at certain interesting locations. Pressure jumps have been found while opening or closing the valves. The valve motion is then modeled with considerations to its geometry. The results show that as the plunger speed reaches a critical value, the valve moves more rapidly and more likely to hit the wall and generates a noise. We also provide a methodology to study the ow across the valve in higher resolution. A nite-di erence approach to the Navier-Stokes equations are presented with the immersed boundary method. / UOIT

Two immersed boundary methods (IBM) for the simulation of conjugate heat transfer problems with complex geometries are introduced: a finite element (IFEM) and a finite volume (IFVM) immersed boundary methods are discussed. In the IFEM a projection approach is presented for the coupled system of time-dependent incompressible Navier-Stokes equations (NSEs) and energy equation in conjunction with the immersed boundary method for solving fluid flow and heat transfer problems in the presence of rigid objects not represented by the underlying mesh. The IBM allows solving the flow for geometries with complex objects without the need of generating a body-fitted mesh. Dirichlet boundary constraints are satisfied applying a boundary force at the immersed body surface. Using projection and interpolation operators from the fluid volume mesh to the solid surface mesh (i.e., the “immersed” boundary) and vice versa, it is possible to impose the extra constraint to the NSEs as a Lagrange multiplier in a fashion very similar to the effect pressure has on the momentum equations to satisfy the divergence-free constraint. The IFEM approach presented shows third order accuracy in space and second order accuracy in time when the simulation results for the Taylor-Green decaying vortex are compared to the analytical solution. For the IFVM a ghost-cell approach with sharp interface scheme is used to enforce the boundary condition at the fluid/solid interface. The interpolation procedure at the immersed boundary preserves the overall second order accuracy of the base solver. The developed ghost-cell method is applied on a staggered configuration with the Semi-Implicit Method for Pressure-Linked Equations Revised algorithm. Second order accuracy in space and first order accuracy in time are obtained when the Taylor-Green decaying vortex test case is compared to the IFVM analytical solution. Computations were performed using the IFEM and IFVM approaches for the two-dimensional flow over a backward-facing step, two-dimensional flow past a stationary circular cylinder, three-dimensional flow past a sphere and two and three-dimensional natural convection in an enclosure with/without immersed body. The numerical results obtained with the discussed IFEM and IFVM were compared against other IBMs available in literature and simulations performed with the commercial computational fluid dynamics code STAR-CCM+/V7.04.006. The benchmark test cases showed that the numerical results obtained with the implemented immersed boundary methods are in good agreement with the predictions from STAR-CCM+ and the numerical data from the other IBMs. The immersed boundary method based of finite element approach is numerically more accurate than the IBM based on finite volume discretization. In contrast, the latter is computationally more efficient than the former.

In this dissertation, we explore direct-forcing immersed boundary methods (IBM) under the framework of the lattice Boltzmann method (LBM), which is called the direct-forcing immersed boundary-lattice Boltzmann method (IB-LBM). First, we derive the direct-forcing formula based on the split-forcing lattice Boltzmann equation, which recovers the Navier-Stokes equation with second-order accuracy and enables us to develop a simple and accurate formula due to its kinetic nature. Then, we assess the various interface schemes under the derived direct-forcing formula. We consider not only diffuse interface schemes but also a sharp interface scheme. All tested schemes show a second-order overall accuracy. In the simulation of stationary complex boundary flows, we can observe that the sharper the interface scheme is, the more accurate the results are. The interface schemes are also applied to moving boundary problems. The sharp interface scheme shows better accuracy than the diffuse interface schemes but generates spurious oscillation in the boundary forcing terms due to the discontinuous change of nodes for the interpolation. In contrast, the diffuse interface schemes show smooth change in the boundary forcing terms but less accurate results because of discrete delta functions. Hence, the diffuse interface scheme with a corrected radius can be adopted to obtain both accurate and smooth results. Finally, a direct-forcing immersed boundary method (IBM) for the thermal lattice Boltzmann method (TLBM) is proposed to simulate non-isothermal flows. The direct-forcing IBM formulas for thermal equations are derived based on two TLBM models: a double-population model with a simplified thermal lattice Boltzmann equation (Model 1) and a hybrid model with an advection-diffusion equation of temperature (Model 2). The proposed methods are validated through natural convection problems with stationary and moving boundaries. In terms of accuracy, the results obtained from the IBMs based on both models are comparable and show a good agreement with those from other numerical methods. In contrast, the IBM based on Model 2 is more numerically efficient than the IBM based on Model 1. Overall, this study serves to establish the feasibility of the direct-forcing IB-LBM as a viable tool for computing various complex and/or moving boundary flow problems.

Study of turbulent flow over a rough bed is highly important due to its numerous applications in the areas of sediment transport and pollutant discharge in streams, rivers and channels. Over the past few decades, many experimental studies have been conducted in this respect to understand the underlying phenomenon. However, there is a scarcity in the number of computational studies conducted on this topic. Therefore, a Large Eddy Simulation (LES) of turbulent flow over a rough channel bed was conducted to contribute further understanding of the influence of bed roughness on turbulent flow properties. For this purpose, an efficient, second order accurate 'immersed boundary method' was implemented into the LES code Hydro3d-GT, and validated for flow past bluff bodies. LES results from the present study showed excellent agreement with previous experimental studies on flow over rough beds. An in-depth analysis of time varying turbulent quantities (like the velocity fluctuations) revealed the presence of coherent structures in the flow. Also, a three dimensional visualization of the turbulent structures provided a good picture of the flow, especially in the near bed region, which is quite difficult to accomplish using experimental studies.