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.

Radio-frequency microelectromechanical systems (RF MEMS) are widely used for contact actuators and capacitive switches. These devices typically consist of a metallic membrane which is activated by a time-periodic electrostatic force and makes periodic contact with a contact pad. The increase in switch capacitance at contact causes the RF signal to be deflected and the switch thus closes. Membrane motion is damped by the surrounding gas, typically air or nitrogen. As the switch opens and closes, the flow transitions between the continuum and rarefied regimes. Furthermore, creep is a critical physical mechanism responsible for the failure in these devices, especially those operating at high RF power. Simultaneous and accurate modeling of all these different physics is required to understand the dynamical membrane response in these devices and to estimate device lifetime and to improve MEMS reliability. It is advantageous to model fluid and structural mechanics and electrostatics within a single comprehensive numerical framework to facilitate coupling between them. In this work, we develop a single unified finite volume method based numerical framework to study this multi-physics problem in RF MEMS. Our objective required us to develop structural solvers, fluid flow solvers, and electrostatic solvers using the finite volume method, and efficient mechanisms to couple these different solvers. A particular focus is the development of flow solvers which work efficiently across continuum and rarefied regimes. A number of novel contributions have been made in this process. Structural solvers based on a fully implicit finite volume method have been developed for the first time. Furthermore, strongly implicit fluid flow solvers have also been developed that are valid for both continuum and rarefied flow regimes and which show an order of magnitude speed-up over conventional algorithms on serial platforms. On parallel platforms, the solution techniques developed in this thesis are shown to be significantly more scalable than existing algorithms. The numerical methods developed are used to compute the static and dynamic response of MEMS. Our results indicate that our numerical framework can become a computationally efficient tool to model the dynamics of RF MEMS switches under electrostatic actuation and gas damping. / text

Cette thèse a pour objectif de comprendre la dynamique d’une sphère traversant une interface liquide-liquide. Cette situation, se rencontre dans de nombreuses applications, allant du cycle du carbone dans l’océan (sédimentation de neige marine), aux procédés d’enrobage, en passant par la détection de phase dans l’industrie pétrolière. Pour étudier cette configuration, trois approches sont privilégiées. Un dispositif expérimental muni d’une caméra haute fréquence est utilisé de manière à explorer la dynamique conjointe de la sphère et de l’interface sur une large gamme de paramètres. Le couplage entre une méthode Volume of Fluid (VoF) et une méthode de frontières immergées (IBM) est réalisé et validé dans le but de simuler numériquement ce problème. Enfin des modèles théoriques sont mis en place de manière à interpréter physiquement les différents comportements observés. Ces trois démarches complémentaires permettent de caractériser le passage d’une configuration de flottaison à l’entraînement colonnaire notamment en fonction du rapport entre effets gravitationnels et capillaires. La dynamique de la colonne emportée est très riche (instabilité capillaire, visqueuse, fragmentation, ...). Le bon accord entre les expériences et les simulations numériques permet d’évaluer avec confiance l’influence de chaque paramètre sans dimension (au nombre de 5) à l’aide d’une étude paramétrique numérique. / The goal of this work is to understand the dynamics of a sphere passing through a liquid-liquid interface. Such a configuration is met in different applications, such as oceanic carbon cycle (sedimentation of marine snow), coating processes and phase detection in oil industry. To this aim, three different aproaches are employed. An experimental device, in which various sets of fluids and spheres are used, has been designed to analyze different types of configuration. A combination of an Immersed Boundary Method (IBM) with a Volume of Fluid (VoF) method is used to compute the flow field. Finally theoretical models are derived to better understand the observed behaviours. These three approaches give insights to understand whether a sphere can float or sink. The behaviour of the tail of light fluid towed by the sphere appears to be extremely rich (capillary and viscous instabilities, fragmentation, ...). The agreement between experimental and numerical results allows us to perform an extensive numerical study of the influence of all dimensionless parameters

Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / Mingjun Wei / In this study, a global Proper Orthogonal Decomposition (POD)-Galerkin based Reduced Order model (ROM) is proposed. It is extended from usual fixed-domain problems to more general fluid-solid systems with moving boundaries/interfaces. The idea of the extension is similar to the immersed boundary method in numerical simulations which uses embedded forcing terms to represent boundary motions and domain changes. This immersed boundary method allows a globally defined fixed domain including both fluid and solid, where POD-Galerkin projection can be directly applied. However, such a modified approach cannot get away with the unsteadiness of boundary terms which appear as time-dependent coefficients in the new Galerkin model. These coefficients need to be pre-computed for prescribed periodic motion, or worse, to be computed at each time step for non-prescribed (e.g. with fluid-structure interaction) or non-periodic situations. Though computational time for each unsteady coefficient is smaller than the coefficients in a typical Galerkin model, because the associated integration is only in the close neighborhood of moving boundaries. The time cost is still much higher than a typical Galerkin model with constant coefficients. This extra expense for moving-boundary treatment eventually undermines the value of using ROMs. An aggressive approach is to decompose the moving boundary/domain to orthogonal modes and derive another low-order model with fixed coefficients for boundary motion. With this domain decomposition, an approach including two coupled low-order models both with fixed coefficients is proposed. Therefore, the new global ROM with decomposed approach is more efficient. Though the model with the domain decomposition is less accurate at the boundary, it is a fair trade-off for the benefit on saving computational cost. The study further shows, however, that the most time-consuming integration in both approaches, which come from the unsteady motion, has almost negligible impact on the overall dynamics. Dropping these time-consuming terms reduces the computation cost by at least one order while having no obvious effect on model accuracy. Based on this global POD-Galerkin based ROM with forcing term, an improved ROM which can handle the parametric variation of body motions in a certain range is also presented. This study shows that these forcing terms not only represent the moving of the boundary, but also decouple the moving parameters from the computation of model coefficients. The decoupling of control parameters provides the convenience to adapt the model for the prediction on states under variation of control parameters. An improved ROM including a shit mode seems promising in model adaptation for typical problems in a fixed domain. However, the benefit from adding a shit mode to model diminishes when the method is applied to moving-boundary problems. Instead, a combined model, which integrates data from a different set of parameters to generate the POD modes, provides a stable and accurate ROM in a certain range of parametric space for moving-boundary problems. By introducing more data from a different set of parameters, the error of the new model can be further reduced. This shows that the combined model can be trained by introducing more and more information. With the idea of the combined model, the improved global ROM with forcing terms shows impressive capability to predict problems with different unknown moving parameters, and can be used in future parametric control and optimization problems.

Un large spectre d’applications en ingénierie est concerné par les écoulements de fluides en interaction avec des structures poreuses, allant de problèmes à petite échelle jusqu’à des problématiques de plus grande échelle. Ces structures poreuses, souvent à géométries complexes, peuvent se déplacer ou se déformer en réponse au forçage exercé par l’écoulement environnant.Le but de ce travail est de proposer un modèle numérique pour la simulation macroscopique d’écoulements de fluide interagissant avec des milieux poreux mobiles à géométries complexes, qui soit facile d’implémentation et pouvant être utilisé dans une large gamme d’applications. Pour atteindre cet objectif, la méthode de Lattice Boltzmann est utilisée pour résoudre l’écoulement dans des milieux poreux à l’échelle d’un volume représentatif élémentaire. Pour l’implémentation du mouvement désiré, le concept de frontières immergées est adopté. Dans ce contexte, un nouveau modèle est proposé pour traiter des milieux poreux en volume, dont la résistance à l’écoulement environnant est modélisé par la loi de Brinkman-Forchheimer-Darcy étendue.L’algorithme est d’abord testé sur l’écoulement à travers un cylindre fixe. La simplicité de ce cas test académique permet de caractériser finement la précision de la méthode. Le modèle est ensuite utilisé pour simuler des écoulements de fluide autour et à travers des corps poreux mobiles, à la fois pour des géométries confinées et pour des écoulements ouverts. L’invariance Galiléenne des équations moyennées macroscopiques gouvernant la dynamique du fluide est démontrée. D’excellents accords avec les résultats de référence sont obtenus pour les différents cas testés. / A wide spectrum of engineering problems is concerned with fluid flows in interaction with porous structures, ranging from small length-scale problems to large ones. These structures, often of complex geometry, may move/deform in response to the forces exerted by the surrounding flow. Despite the advancements in computational fluid dynamics, the numerical simulation of such configurations - a valuable tool for the study of the flow physics involved - remains a challenging task.The aim of the present work is to propose a numerical model for the macroscopic simulation of fluid flows interacting with moving porous media of complex geometry, that is easy to implement and can be used in a range of applications. To achieve this, the Lattice Boltzmann method is employed for solving the flow in porous media at the representative elementary volume scale. For the implementation of the desired body motion, the concept of the Immersed Boundary method is adopted. In this context, a novel model is proposed for dealing with moving volumetric porous media, whose resistance to the surrounding flow obeys the Brinkman-Forchheimer-extended Darcy law. The algorithm is initially tested for flow past a static cylinder. The simplicity of this academic test case allows us to assess in detail the accuracy of the proposed method. The model is later used to simulate fluid flows around and through moving porous bodies, both in a confined geometry and in open space. We are able to demonstrate the Galilean invariance of the macroscopic volume-averaged flow governing equations. Excellent agreement with reference results is obtained in all cases.

The objective of this work is to develop a novel 3-D biological particulate dynamics framework to simulate blood flow in the micro circulation. This entails the amalgamation of concepts from various fields namely blood flow dynamics, solid mechanics, fluid-structure interaction and computational data structures. It is envisioned that this project will serve as a harbinger for implementing a multi-scale simulation model with applications in a vast array of situations from blood flows in heart valves to studying cancer metastasis. The primary motivation for this work stems from the need for establishing a simple, effective and holistic framework for performing blood flow simulations, taking into account the extremely 3-D nature of flow, the particle interactions and fluid structure interaction between blood and its constituent elements. Many current models to simulate blood cells rely on finite element methods which render large scale simulations extremely computationally intensive. The development of a framework for simulating blood flow is tied together with achieving a framework for performing an investigation of cancer metastasis. Cancer initially develops at a primary site and spreads through the body to secondary sites using the circulatory systems of the body - the blood circulatory system and the lymphatic system. It is known that all the cancer cells that enter into the circulation do not survive the harsh environment, though the exact cause of this is still undetermined. Moreover, the mechanical properties of cancer cells are not well documented and appropriate computational models require that experiments be conducted to determine the same. Thus the end goal of this work is to establish a system to analyze and simulate 3-D blood particulate dynamics, including cancer cells, from a holistic standpoint in order to understand more about the phenomenon of blood flow as a whole, and cancer metastasis in particular.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第18271号 / 工博第3863号 / 新制||工||1592(附属図書館) / 31129 / 京都大学大学院工学研究科航空宇宙工学専攻 / (主査)教授 稲室 隆二, 教授 泉田 啓, 教授 青木 一生 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

No description available.

Knowledge of the detailed flow dynamics at the interior of branching corals is critical for a full understanding of nutrient uptake, mass transport, wave dissipation, and other essential processes. These physiological processes depend on the local velocity field, local concentration gradients of nutrients and waste, and the turbulent stresses developed on and above the coral surface. Though the large-scale hydrodynamics over coral reefs are well studied, the interior hydrodynamics, between the branches, remains uncharacterized due to limited optical and acoustic access to the interior. In the current thesis, a three-dimensional immersed boundary method in the large eddy simulation framework was used to compute the flow inside several branching coral colony geometries in order to study the effects of branch density and surface structure on the flow fields in the coral interiors. Two different Pocillopora colony species were studied at different Reynolds numbers. A ray-tracing algorithm was used for capturing the arbitrary branches of these complex geometries to obtain the three-dimensional flow fields within these colonies for the first time. The analysis showed the formation of vortices at the colony interior that stir the water column and thus passively enhance mass transport, compensating for the reduced mean velocity magnitude compared to the free stream value, within the densely branched Pocillopora meandrina colony. Further analysis showed that the mean streamwise velocity profile changes shape along the streamwise direction inside P. meandrina, whereas the mean velocity profile did not change shape from the front to the back for the loosely branched Pocillopora colony, Pocillopora eydouxi. Moreover, turbulent flow field quantities were computed for both these structures, and for two almost identical Montipora capitata colony geometries, one with, and one without roughness elements called verrucae. The analyses demonstrated significant differences in the mean velocity profiles, Reynolds stress, and other flow quantities with changes in colony branch density and surface structure. / Doctor of Philosophy / Coral reefs are the largest marine ecosystem, and play a critical role in protecting coastal areas against flooding and erosion. The majority of the world's corals are currently under threat from rising ocean temperatures, which disrupt the symbiotic relationship between the coral polyp and its symbiont algae causing coral bleaching. Bleaching involves processes mediated by the flow at the coral surface, but relatively little is known regarding the local flow dynamics between the branches of coral reefs. The current research seeks to characterize internal coral hydrodynamics, leading to insights about many critical physiological and other processes in corals, like drag formation, mixing, and mass or nutrient transport to and from the coral. In the current study, the influence of the coral branch density and surface structure on the resulting interbranch flow field were investigated by simulating the flow resulting from uniform oncoming ocean flow conditions using three-dimensional immersed boundary large eddy simulations. The detailed velocity and pressure fields were found throughout the interior of the colonies studied. A distinct mass transport mechanism was found inside one densely branched colony studied. For this coral, Pocillopora meandrina, the flow speed reduces substantially inside the coral because of the high branch density. But corals depend on the ocean flow to bring nutrients to the polyps on their surface. We found that P. meandrina sheds hundreds of small vortices from its branches, which stir the overlying water column, increasing the mass transport rate, and compensating almost exactly for the reduced flow in the interior. The study also included computing the flow through three other coral colony geometries, and comparisons of their mean velocity profiles and turbulent flow statistics in order to examine the impact of the colony branch density and surface structures on the resulting hydrodynamic flow field. The current investigation of coral hydrodynamics may lead to an increased understanding of coral health and physiological activity, and may help in designing effective interventions for the challenges facing corals, which could have impacts in the fields of coral restoration, coastal protection, and public policy in the United States and abroad.