We give a concise explication of the theory of level set methods for modeling motion of an interface as well as the numerical implementation of these methods. We then introduce the geometry of a camera and the mathematical models for 3D reconstruction with a few examples both simulated and from a real camera. We finally describe the model for 3D surface reconstruction from n-camera views using level set methods.

Under the right conditions, a moving sphere may successfully enter, and leave, a soap bubble without rupturing that bubble. The physics behind this phenomena are not well understood, nor the limiting factors (such as sphere size, speed, etc.). This work, investigating this phenomenon using high speed photography, has produced several results which are presented. First, several distinct regimes, noted while photographing the interactions between the spheres and the bubbles, are classified and discussed. Next a probabilistic examination of the soap bubbles rupture by the moving spheres is presented. Then a conjecture for the limiting sphere sizes and speeds is presented. And finally some interesting phenomena, noted in the course of this investigation, are presented and discussed.

京都大学 / 新制・課程博士 / 博士(工学) / 甲第24226号 / 工博第5054号 / 新制||工||1789(附属図書館) / 京都大学大学院工学研究科機械理工学専攻 / (主査)教授 平山 朋子, 教授 岩井 裕, 教授 松原 厚 / 学位規則第4条第1項該当 / Doctor of Philosophy (Engineering) / Kyoto University / DFAM

No description available.

This dissertation addresses path planning for an autonomous vehicle navigating in a two dimensional environment for which an a priori map is inaccurate and for which the environment is sensed in real-time. For this class of application, planning decisions must be made in real-time. This work is motivated by the need for fast autonomous vehicles that require planning algorithms to operate as quickly as possible. In this dissertation, we first study the case in which there are only static obstacles in the environment. We propose a hybrid receding horizon control path planning algorithm that is based on level-set methods. The hybrid method uses global or local level sets in the formulation of the receding horizon control problem. The decision to select a new level set is made based on certain matching conditions that guarantee the optimality of the path. We rigorously prove sufficient conditions that guarantee that the vehicle will converge to the goal as long as a path to the goal exists. We then extend the proposed receding horizon formulation to the case when the environment possesses moving obstacles. Since all of the results in this dissertation are based on level-set methods, we rigorously investigate how level sets change in response to new information locally sensed by a vehicle. The result is a dynamic fast marching algorithm that usually requires significantly less computation that would otherwise be the case. We demonstrate the proposed dynamic fast marching method in a successful field trial for which an autonomous surface vehicle navigated four kilometers through a riverine environment. / Ph. D.

Computational micromechanics analysis of carbon nanotube-epoxy nanocomposites, containing aligned nanotubes, is performed using the mesh independent intrinsic extended finite element method (IXFEM). The IXFEM employs a localized intrinsic enrichment strategy to treat arbitrary discontinuities defined through the level-set method separate from the problem domain discretization, i.e. the finite element (FE) mesh. A global domain decomposition identifies local subdomains for building distinct partition of unities that appropriately suit the approximation. Specialized inherently enriched shape functions, constructed using the moving least square method, enhance the approximation space in the vicinity of discontinuity interfaces, maintaining accuracy of the solution, while standard FE shape functions are used elsewhere. Comparison of the IXFEM in solving validation problems with strong and weak discontinuities against a standard finite element method (FEM) and analytic solutions validates the enriched intrinsic bases, and shows anticipated trends in the error convergence rates. Applying the IXFEM to model composite materials, through a representative volume element (RVE), the filler agents are defined as individual weak bimaterial interfaces. Though a series of RVE studies, calculating the effective elastic material properties of carbon nanotube-epoxy nanocomposite systems, the benefits in substituting the conventional mesh dependent FEM with the mesh independent IXFEM when completing micromechanics analysis, investigating effects of high filler count or an evolving microstructure, are demonstrated. / Master of Science

A otimização topológica de estruturas está relacionada à concepção de projetos que executem suas funções com nível de segurança adequado empregando a quantidade mínima de material. Neste trabalho, determina-se a geometria ótima de estruturas planas por meio do acoplamento do Método dos Elementos de Contorno (MEC) ao Método Level Set (MLS). O algoritmo é composto por 3 etapas: problema mecânico, otimização topológica e reconstrução da estrutura. O problema mecânico é resolvido pelas equações algébricas do MEC. A otimização topológica é determinada pelo MLS, este representa a geometria do corpo e suas evoluções por meio da função Level Set (LS) avaliada em seu nível zero. Na reconstrução realiza-se o remalhamento, pois a cada iteração a estrutura é modificada. O acoplamento proposto resulta na geometria ótima da estrutura sem a necessidade da aplicação de filtros. Os exemplos analisados mostram que algoritmo desenvolvido capta adequadamente a geometria ótima das estruturas. Com esse trabalho, avança-se no campo das aplicações do acoplamento MEC-MLS e no desenvolvimento de soluções inovadoras para problemas complexos de engenharia. / In general, the topological optimization of structures is related to design projects that perform their functions with appropriate security levels using the minimum amount of material. This research determines the optimal geometry of 2D structures by coupling the Boundary Blement Method (BEM) to Level Set Method (LSM). The algorithm consists of 3 steps: mechanical model, topology optimization and structure reconstruction. The mechanical model is solved by BEM algebraic equations. The topology optimization is determined using the MLS, the geometry of the body is determined by the Level Set (LS) function evaluated at the zero level. The reconstruction achieves the remeshing, because for each iteration of the structure is modified. The proposed coupling results in the optimal geometry of the structure without the filters application. The examples show that the algorithm developed captures adequately the optimal geometry of the structures. With this dissertation, it is possible advance in the field of applications of the BEM - LSM and develop innovative solutions to complex engineering problems.

A otimização topológica de estruturas está relacionada à concepção de projetos que executem suas funções com nível de segurança adequado empregando a quantidade mínima de material. Neste trabalho, determina-se a geometria ótima de estruturas planas por meio do acoplamento do Método dos Elementos de Contorno (MEC) ao Método Level Set (MLS). O algoritmo é composto por 3 etapas: problema mecânico, otimização topológica e reconstrução da estrutura. O problema mecânico é resolvido pelas equações algébricas do MEC. A otimização topológica é determinada pelo MLS, este representa a geometria do corpo e suas evoluções por meio da função Level Set (LS) avaliada em seu nível zero. Na reconstrução realiza-se o remalhamento, pois a cada iteração a estrutura é modificada. O acoplamento proposto resulta na geometria ótima da estrutura sem a necessidade da aplicação de filtros. Os exemplos analisados mostram que algoritmo desenvolvido capta adequadamente a geometria ótima das estruturas. Com esse trabalho, avança-se no campo das aplicações do acoplamento MEC-MLS e no desenvolvimento de soluções inovadoras para problemas complexos de engenharia. / In general, the topological optimization of structures is related to design projects that perform their functions with appropriate security levels using the minimum amount of material. This research determines the optimal geometry of 2D structures by coupling the Boundary Blement Method (BEM) to Level Set Method (LSM). The algorithm consists of 3 steps: mechanical model, topology optimization and structure reconstruction. The mechanical model is solved by BEM algebraic equations. The topology optimization is determined using the MLS, the geometry of the body is determined by the Level Set (LS) function evaluated at the zero level. The reconstruction achieves the remeshing, because for each iteration of the structure is modified. The proposed coupling results in the optimal geometry of the structure without the filters application. The examples show that the algorithm developed captures adequately the optimal geometry of the structures. With this dissertation, it is possible advance in the field of applications of the BEM - LSM and develop innovative solutions to complex engineering problems.

Patient-specific multi-physics simulations have the potential to improve the diagnosis, treatment, and scientific inquiry of heart valve dynamics. It has been shown that the flow characteristics within the left ventricle are important to correctly capture the aortic and mitral valve motion and corresponding fluid dynamics, motivating the use of patient-specific imaging to describe the aortic and mitral valve geometries as well as the motion of the left ventricle (LV). The LV position can be captured at several time points in the cardiac cycle, such that its motion can be prescribed a priori as a Dirichlet boundary condition during a simulation. Valve leaflet motion, however, should be computed from soft-tissue models and incorporated using fully-coupled Fluid Structure Interaction (FSI) algorithms. While FSI simulations have in part or wholly been achieved by multiple groups, to date, no high-throughput models have been developed, which are needed for use in a clinical environment. This project seeks to enable patient-derived moving LV boundary conditions, and has been developed for use with a previously developed immersed boundary, fixed Cartesian grid FSI framework. One challenge in specifying LV motion from medical images stems from the low temporal resolution available. Typical imaging modalities contain only tens of images during the cardiac cycle to describe the change in position of the left ventricle. This temporal resolution is significantly lower than the time resolution needed to capture fluid dynamics of a highly deforming heart valve, and thus an approach to describe intermediate positions of the LV is necessary. Here, we propose a primarily Eulerian means of representing LV displacement. This is a natural extension, since an Eulerian framework is employed in the CFD model to describe the large displacement of the heart valve leaflets. This approach to using Eulerian interface representation is accomplished by applying “morphing” techniques commonly used in the field of computer graphics. For the approach developed in the current work, morphing is adapted to the unique characteristics of a Cartesian grid flow solver which presents challenges of adaptive mesh refinement, narrow band approach, parallel domain decomposition, and the need to supply a local surface velocity to the flow solver that describes both normal and tangential motion. This is accomplished by first generating a skeleton from the Eulerian interface representation, and deforming the skeleton between image frames to determine bulk displacement. After supplying bulk displacement, local displacement is determined using the Eulerian fields. The skeletons are also utilized to automate the simulation setup to track the locations upstream and downstream where the system inflow/outflow boundary conditions are to be applied, which in the current approach, are not limited to Cartesian domain boundaries.

This thesis presents numerical techniques for solving problems of incompressible flow coupled to scalar transport equations using finite element discretizations in space. The two applications considered in this thesis are multi-phase flow, modeled by level set or phase field methods, and planetary mantle convection based on the Boussinesq approximation. A systematic numerical study of approximation errors in evaluating the surface tension in finite element models for two-phase flow is presented. Forces constructed from a gradient in the same discrete function space as used for the pressure are shown to give the best performance. Moreover, two approaches for introducing contact line dynamics into level set methods are proposed. Firstly, a multiscale approach extracts a slip velocity from a micro simulation based on the phase field method and imposes it as a boundary condition in the macro model. This multiscale method is shown to provide an efficient model for the simulation of contact-line driven flow. The second approach combines a level set method based on a smoothed color function with a the phase field method in different parts of the domain. Away from contact lines, the additional information in phase field models is not necessary and it is disabled from the equations by a switch function. An in-depth convergence study is performed in order to quantify the benefits from this combination. Also, the resulting hybrid method is shown to satisfy an a priori energy estimate. For the simulation of mantle convection, an implementation framework based on modern finite element and solver packages is presented. The framework is capable of running on today's large computing clusters with thousands of processors. All parts in the solution chain, from mesh adaptation over assembly to the solution of linear systems, are done in a fully distributed way. These tools are used for a parallel solver that combines higher order time and space discretizations. For treating the convection-dominated temperature equation, an advanced stabilization technique based on an artificial viscosity is used. For more efficient evaluation of finite element operators in iterative methods, a matrix-free implementation built on cell-based quadrature is proposed. We obtain remarkable speedups over sparse matrix-vector products for many finite elements which are of practical interest. Our approach is particularly efficient for systems of differential equations.