Red Blood Cell Dynamics on Non-Uniform Grids using a Lattice Boltzmann Flux Solver and a Spring-Particle Red Blood Cell Model

The Computational Haemodynamics Research Group (CHRG) in Technological University Dublin is developing a computational fluid dynamics (CFD) software package aimed specifically at physiologically-realistic modelling of blood flow. A physiologically-realistic model of blood flow involves calculating the deformation of individual red blood cells (RBCs) and the contribution of this deformation to the overall blood flow. The CHRG has developed an enhanced spring-particle RBC structural model that is capable of modelling the full stomatocyte-discocyteechinocyte (SDE) transformation. This RBC model, incorporated into a fluid dynamics solver, will provide a physiologically-realistic blood flow model. In this work the overall plasma flow is modelled using a novel technique: the lattice Boltzmann flux solver (LBFS). This is an innovative approach to solving the NavierStokes (N-S) equations for fluid flow. It involves solving the macroscopic equations using the finite volume method (FVM) and calculating the flux across the cell interfaces via a local reconstruction of the lattice Boltzmann equation (LBE). Fluidstruture interaction between the RBC and the plasma is captured by coupling the RBC solver to the LBFS via the immersed boundary method (IBM). Numerical experiments investigating RBC dynamics are performed using non-uniform grids and validated against existing experimental data in the literature. Finally all numerical solvers are developed using general purpose GPU programming (GPGPU) and this is shown to accelerate simulation runtimes significantly

Computational analysis of woodwind instruments using the lattice Boltzmann method

Through the use of the lattice Boltzmann method, a series of computational models were created to simulate air flow in woodwind instruments. Start- ing as a two-dimensional code in Matlab running on the CPU, the model went through a series of iterations before becoming a three-dimension code in Fortran that was accelerated through the use of GPU parallel computing. The accuracy and stability of the model are shown by comparison to various published benchmark tests. Thus far, the air flow in organ pipes for a two dimensional model was simulated showing oscillating flow by the labium as expected. This thesis offers the mathematical and computational background as well as a description of the implementation of the basic LBM for simulat- ing flow in musical instruments. The method described here is meant as a first step to a code that is highly flexible and can be used to study many as- pects of acoustics in musical instruments. Future applicability of the model includes observing flow at the exit of both square and round organ pipes in addition to modeling the reed-mouthpiece system of the clarinet

Kinetic Solvers with Adaptive Mesh in Phase Space

An Adaptive Mesh in Phase Space (AMPS) methodology has been developed for solving multi-dimensional kinetic equations by the discrete velocity method. A Cartesian mesh for both configuration (r) and velocity (v) spaces is produced using a tree of trees data structure. The mesh in r-space is automatically generated around embedded boundaries and dynamically adapted to local solution properties. The mesh in v-space is created on-the-fly for each cell in r-space. Mappings between neighboring v-space trees implemented for the advection operator in configuration space. We have developed new algorithms for solving the full Boltzmann and linear Boltzmann equations with AMPS. Several recent innovations were used to calculate the discrete Boltzmann collision integral with dynamically adaptive mesh in velocity space: importance sampling, multi-point projection method, and the variance reduction method. We have developed an efficient algorithm for calculating the linear Boltzmann collision integral for elastic and inelastic collisions in a Lorentz gas. New AMPS technique has been demonstrated for simulations of hypersonic rarefied gas flows, ion and electron kinetics in weakly ionized plasma, radiation and light particle transport through thin films, and electron streaming in semiconductors. We have shown that AMPS allows minimizing the number of cells in phase space to reduce computational cost and memory usage for solving challenging kinetic problems

A general pressure equation based method for incompressible two-phase flows

We present a fully-explicit, iteration-free, weakly-compressible method to simulate immiscible incompressible two-phase flows. To update pressure, we circumvent the computationally expensive Poisson equation and use the general pressure equation. In addition, the volume-of-fluid approach is used for interface capturing under the operator-split methodology. Our method is fully-explicit and stable with simple local spatial discretization, and hence, it is easy to implement. Several two- and three-dimensional canonical two-phase flows are simulated. The qualitative and quantitative results prove that our method is capable of accurately handling problems involving a range of density and viscosity ratios and surface tension effects

Real-Time Simulation of Indoor Air Flow using the Lattice Boltzmann Method on Graphics Processing Unit

This thesis investigates the usability of the lattice Boltzmann method (LBM) for the simulation of indoor air flows in real-time. It describes the work undertaken during the three years of a Ph.D. study in the School of Mechanical Engineering at the University of Leeds, England. Real-time fluid simulation, i.e. the ability to simulate a virtual system as fast as the real system would evolve, can benefit to many engineering application such as the optimisation of the ventilation system design in data centres or the simulation of pollutant transport in hospitals. And although real-time fluid simulation is an active field of research in computer graphics, these are generally focused on creating visually appealing animation rather than aiming for physical accuracy. The approach taken for this thesis is different as it starts from a physics based model, the lattice Boltzmann method, and takes advantage of the computational power of a graphics processing unit (GPU) to achieve real-time compute capability while maintaining good physical accuracy. The lattice Boltzmann method is reviewed and detailed references are given a variety of models. Particular attention is given to turbulence modelling using the Smagorinsky model in LBM for the simulation of high Reynolds number flow and the coupling of two LBM simulations to simulate thermal flows under the Boussinesq approximation. A detailed analysis of the implementation of the LBM on GPU is conducted. A special attention is given to the optimisation of the algorithm, and the program kernel is shown to achieve a performance of up to 1.5 billion lattice node updates per second, which is found to be sufficient for coarse real-time simulations. Additionally, a review of the real-time visualisation integrated within the program is presented and some of the techniques for automated code generation are introduced. The resulting software is validated against benchmark flows, using their analytical solutions whenever possible, or against other simulation results obtained using accepted method from classical computational fluid dynamics (CFD) either as published in the literature or simulated in-house. The LBM is shown to resolve the flow with similar accuracy and in less time

TURBULENT TRANSITION SIMULATION AND PARTICULATE CAPTURE MODELING WITH AN INCOMPRESSIBLE LATTICE BOLTZMANN METHOD

Derivation of an unambiguous incompressible form of the lattice Boltzmann equation is pursued in this dissertation. Further, parallelized implementation in developing application areas is researched. In order to achieve a unique incompressible form which clarifies the algorithm implementation, appropriate ansatzes are utilized. Through the Chapman-Enskog expansion, the exact incompressible Navier-Stokes equations are recovered. In initial studies, fundamental 2D and 3D canonical simulations are used to evaluate the validity and application, and test the required boundary condition modifications. Several unique advantages over the standard equation and alternative forms found in literature are found, including faster convergence, greater stability, and higher fidelity for relevant flows. Direct numerical simulation and large eddy simulation of transitional and chaotic flows are one application area explored with the derived incompressible form. A multiple relaxation time derivation is performed and implemented in a 2D cavity (direct simulation) and a 3D cavity (large eddy simulation). The Kolmogorov length scale, a function of Reynolds number, determines grid resolution in the 2D case. Comparison is made to the extensive literature on laminar flows and the Hopf bifurcation, and final transition to chaos is predicted. Steady and statistical properties in all cases are in good agreement with literature. In the 3D case the relatively new Vreman subgrid model provides eddy viscosity modeling. By comparing the center plane to the direct numerical simulation case, both steady and unsteady flows are found to be in good agreement, with a coarse grid, including prediction of the Hopf bifurcation. Multiphysics pore scale flow is the other main application researched here. In order to provide the substrate geometry, a straightforward algorithm is developed to generate random blockages producing realistic porosities and passages. Combined with advection-diffusion equations for conjugate heat transfer and soot particle transport, critical diesel particulate filtration phenomena are simulated. To introduce additional fidelity, a model is added which accounts for deposition caused by a variety of molecular and atomic forces. Detailed conclusions are presented to lay the groundwork for future extensions and improvements. Predominantly, higher lattice velocity large eddy simulation, improved parallelization, and filter regeneration