In science and engineering, many simulations are carried out over domains consisting of multiple materials separated by curves/surfaces. If partial differential equations (PDEs) are used to model these simulations, it usually leads to the so-called interface problems of PDEs whose coefficients are discontinuous. In this dissertation, we consider nonconforming immersed "nite element (IFE) methods and error analysis for interface problems. We "first consider the second order elliptic interface problem with a discontinuous diffusion coefficient. We propose new IFE spaces based on the nonconforming rotated Q1 "finite elements on Cartesian meshes. The degrees of freedom of these IFE spaces are determined by midpoint values or average integral values on edges. We investigate fundamental properties of these IFE spaces, such as unisolvency and partition of unity, and extend well-known trace inequalities and inverse inequalities to these IFE functions. Through interpolation error analysis, we prove that these IFE spaces have optimal approximation capabilities. We use these IFE spaces to develop partially penalized Galerkin (PPG) IFE schemes whose bilinear forms contain penalty terms over interface edges. Error estimation is carried out for these IFE schemes. We prove that the PPG schemes with IFE spaces based on integral-value degrees of freedom have the optimal convergence in an energy norm. Following a similar approach, we prove that the interior penalty discontinuous Galerkin schemes based on these IFE functions also have the optimal convergence. However, for the PPG schemes based on midpoint-value degrees of freedom, we prove that they have at least a sub-optimal convergence. Numerical experiments are provided to demonstrate features of these IFE methods and compare them with other related numerical schemes. We extend nonconforming IFE schemes to the planar elasticity interface problem with discontinuous Lam"e parameters. Vector-valued nonconforming rotated Q1 IFE functions with integral-value degrees of freedom are unisolvent with appropriate interface jump conditions. More importantly, the Galerkin IFE scheme using these vector-valued nonconforming rotated Q1 IFE functions are "locking-free" for nearly incompressible elastic materials. In the last part of this dissertation, we consider potential applications of IFE methods to time dependent PDEs with moving interfaces. Using IFE functions in the discretization in space enables the applicability of the method of lines. Crank-Nicolson type fully discrete schemes are also developed as alternative approaches for solving moving interface problems. / Ph. D.

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

Consider a bottom-hinged Oscillating Wave Surge Converter (OWSC): This device oscillates due to the hydrodynamic forces applied on it by the action of ocean waves. The focus of this thesis is to build upon the in-house multi-block generalized coordinate finite volume solver GenIDLEST using a collocated grid arrangement within the framework of the fractional-step method to make it compatible to simulate such systems. The first step in this process is to deploy a convection scheme which differentiates between air and water. This process is further complicated by the 1:1000 density and 1:100 viscosity ratio between the two fluids. For this purpose, a phase field method is chosen for its ease of implementation and proven boundedness and conservativeness properties. Extensive validation and verification using standard test cases, such as droplet in shear flow, Rayleigh Taylor instability, and the Dam Break Problem is carried out. This development is then coupled with the present Immersed Boundary Module which is used to simulate the presence of moving bodies and again verified against test cases, such as the Dam Break problem with a vertical obstacle and heave decay of a partially submerged buoyant cylinder. Finally, a relaxation zone technique is used to generate waves and a numerical beach technique is used to absorb them. These are then used to simulate the Oscillating Surge Wave Converter. / Master of Science / An Oscillating Wave Surge Converter can be best described as a rectangular flap, hinged at the bottom, rotating under the influence of ocean waves from which energy is harvested. The singular aim of this thesis is to model this device using Computational Fluid Dynamics (CFD). More specifically, the aim is to model this dynamic device with the full Navier Stokes Equations, which include inertial forces, arising due to the motion of the fluid, viscous forces which dissipate energy, and body forces such as gravity. This involves three key steps: 1. Modelling the air-water interface using a convection scheme. A phase field method is used to differentiate between the two fluids. This task is made more challenging because of the very large density and viscosity differences between air and water. 2. Model dynamic moving geometries in a time-dependent framework. For this, we rely on the Immersed Boundary Method. 3. Develop a numerical apparatus to generate and absorb ocean waves. For this, we rely on the Relaxation Zone and Numerical Beach Method. These developments are validated in different canonical problems and finally applied to a two-dimensional oscillating surge wave energy converter.

In this dissertation, I provide a robust, efficient inverse mapping algorithm for use in immersed simulation methods, specifically in the Flex Representation Method. I also explore a structural theory that unifies the theories of solids, shells, beams, and rigid bodies. As part of this, I preform a preliminary exploration of applying the Flex Representation Method to shells. Finally, I explore why higher order elements suffer from small critical time steps in explicit dynamics. I then propose a simple method of remedying this issue by exploiting the properties of U-splines.

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.

This thesis utilizes the Immersed Interface Finite Element Method (IIFEM) for shape optimization and material design, while also investigating the modeling and parameterization of lung tissue for Diver Underwater Explosion (UNDEX) simulations. In the first part, a shape optimization scheme utilizing a four-noded rectangular immersed-interface element is presented. This method eliminates the need for interface fitted mesh or mesh morphing, reducing computational costs while maintaining solution accuracy. Analytical design sensitivity analysis is performed to obtain gradients for the optimization formulation, and various parametrization techniques are explored. The effectiveness of the approach is demonstrated through verification and case studies. For material design, the study combines topological shape optimization with IIFEM, providing a computational approach for architecting materials with desired effective properties. Numerical homogenization evaluates effective properties, and level set-based topology optimization evolves boundaries within the unit cell to generate optimal periodic microstructures. The design space is parameterized using radial basis functions, facilitating a gradient-based optimization algorithm for optimal coefficients. The method produces geometries with smooth boundaries and distinct interfaces, demonstrated through numerical examples. The thesis then delves into modeling the mechanical response of lung tissues, particularly focusing on hyperelastic and hyperviscoelastic models. The research adopts a phased approach, emphasizing hyperelastic model parametrization while reserving hyperviscoelastic model parametrization for future studies. Alternative methods are used for parametrization, circumventing direct experimental tests on biological materials. Representative material properties are sourced from literature or refit from existing experimental data, incorporating both empirically derived data and practical values suitable for simulations. Damage parameter quantification relies on asserted quantitative relationships between injury levels and the regions or percentages of affected lung tissue. / Doctor of Philosophy / This research explores the following themes: optimizing shapes, designing materials using repetitive identical building blocks, and understanding how divers' lungs respond to underwater explosions. When computationally analyzing structures with multiple materials, the conventional method involves creating meshes that align with material interfaces, which can be intricate and time-consuming. The Immersed Interface Finite Element Method (IIFEM) is introduced as a computational approach that simplifies this process, utilizing a uniform grid for analysis regardless of interface shape. Consider a plate with a hole or other inclusions. Shape optimization seeks the optimal hole/inclusion shape for withstanding specific loading. Traditional optimization processes necessitate iterative mesh recreation, a step circumvented by employing IIFEM. This technique also extends to creating micro-building blocks of materials, enabling the architectural design of materials with desired qualities. Materials with specific properties, like strength or flexibility can be achieved. This thesis also addresses the challenge of understanding how divers' lungs respond to underwater explosions, a crucial aspect of safety. Advanced computer models are used to mimic the behavior of lung tissue under shock loads. Directly testing materials and tissues can be difficult and restricted. Techniques like gathering data from scientific papers and refitting existing experimental data are utilized to obtain the information needed. Also, it is hard to directly measure how much damage an underwater explosion does to a diver's lungs. Thus, the level of damage was quantified based on assertions about the relationship between different injury severities and how much lung tissue is affected.

No description available.

京都大学 / 0048 / 新制・課程博士 / 博士(工学) / 甲第21753号 / 工博第4570号 / 新制||工||1712(附属図書館) / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 西脇 眞二, 教授 松原 厚, 教授 黒瀬 良一 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

In this dissertation we discuss bilinear immersed finite elements (IFE) for solving interface problems. The related research works can be categorized into three aspects: (1) the construction of the bilinear immersed finite element spaces; (2) numerical methods based on these IFE spaces for solving interface problems; and (3) the corresponding error analysis. All of these together form a solid foundation for the bilinear IFEs. The research on immersed finite elements is motivated by many real world applications, in which a simulation domain is often formed by several materials separated from each other by curves or surfaces while a mesh independent of interface instead of a body-fitting mesh is preferred. The bilinear IFE spaces are nonconforming finite element spaces and the mesh can be independent of interface. The error estimates for the interpolation of a Sobolev function in a bilinear IFE space indicate that this space has the usual approximation capability expected from bilinear polynomials, which is <i>O</i>(<i>h</i>²) in <i>L</i>² norm and <i>O</i>(<i>h</i>) in <i>H</i>¹ norm. Then the immersed spaces are applied in Galerkin, finite volume element (FVE) and discontinuous Galerkin (DG) methods for solving interface problems. Numerical examples show that these methods based on the bilinear IFE spaces have the same optimal convergence rates as those based on the standard bilinear finite element for solutions with certain smoothness. For the symmetric selective immersed discontinuous Galerkin method based on bilinear IFE, we have established its optimal convergence rate. For the Galerkin method based on bilinear IFE, we have also established its convergence. One of the important advantages of the discontinuous Galerkin method is its flexibility for both <i>p</i> and <i>h</i> mesh refinement. Because IFEs can use a mesh independent of interface, such as a structured mesh, the combination of a DG method and IFEs allows a flexible adaptive mesh independent of interface to be used for solving interface problems. That is, a mesh independent of interface can be refined wherever needed, such as around the interface and the singular source. We also develop an efficient selective immersed discontinuous Galerkin method. It uses the sophisticated discontinuous Galerkin formulation only around the locations needed, but uses the simpler Galerkin formulation everywhere else. This selective formulation leads to an algebraic system with far less unknowns than the immersed DG method without scarifying the accuracy; hence it is far more efficient than the conventional discontinuous Galerkin formulations. / Ph. D.