The purpose of the thesis was to explore the possibilities of using a Level-Set method to design a time-optimal path planar of a subject to direction-dependent maximum velocities. A promising application for such a planning approach lies in sailboat navigation planning, because of the dynamic ocean waves, current, wind and the characteristics of a sailboat. In the thesis, we developed an IOS application to simulate such scenario as environment properties with wind, static obstacles and the sailboat mapped into direction-dependent velocities in different locations of the environment.  Considering the wind is the main power for the sailboat, a wind speed generation function was created, based on different locations. The Level-Set method is widely used in image processing because of its various advantages, for instance, the ability to deal with topology change and stability. It also can be applied in path planning, in which the process of the Level-Set method can be considered as a continuous wave front propagating with a speed from the start location. A grid-based map was used to represent the environment. While the wave front was crossing the cell on the grid, a time was recorded for every cell, following the negative gradient direction of such crossing time, and then an optimal path could be found. In addition, we used the Narrow Band method to speed up the calculation of processing the level set equation. Finally, this report gives the results of the experiments of static obstacle avoidance, wind effects and smooth path planning.

This thesis discusses how to use topology and shape optimization, specifically the level set method, for innovative design. The level set method is a numerical algorithm that simulates the expansion of dynamic implicit surfaces. In this research, the equations for manufacturability are generated and solved through use of the level set method joined with the COMSOL multi-physics package. Specific constraints are added to make the optimization practical for engineering design. The resulting method was applied to design the best underlying support structure, conforming to both curvature and manufacturability constraints, for the longerons used with the International Space Station solar panels.

An interface-preserving level set method is incorporated into the Reynolds-Averaged Navier-Stokes (RANS) numerical method to simulate the application of the green water phenomena around a platform and the breaking wave above the deck. In the present study, this method is used to evaluate the laminar in two dimension plane with fixed orthogonal grids. In this method, it is assumed that the free surface is modeled as immiscible two-phase flow (air and water). A level set function can present the individual fluids, and the interface between two-phase is represented by the zero level set. In addition, the level set evolution equation is coupled with the conservation equations for mass and momentum, which will be solved in the transformed plane. For different purposes, there are several block domains in the application grid. Chimera domain decomposition technique is employed to handle such embedding, overlapping, or matching grids. Several simple test cases were performed to demonstrate the feasibility of this method. The comparisons between the ENO scheme and the WENO scheme will be illustrated in the Zalesak's disk case and will further prove that the WENO scheme is superior to the ENO scheme. The propagation of continuous wave case will validate some properties of wave and determine the importance of some parameters in code. Moreover, the method will be applied in simulation of green water around a two dimensional platform. By configuring different deck heights, some distinct phenomena can be represented. Lastly, it is crucial to observe the green water phenomena around the platform deck by applying the velocity-extrapolation routine.

Pore-scale simulation has become increasingly important in recent years as a tool to understand multiphase flow behavior. Wettability affects aspects of flow such as capillary-pressure saturation curves, residual saturation of each phase, and relative permeability. Simulation of wettability at the pore-scale is still a non-trivial problem, and many different approaches exist to model it. In this work, we implement a variational level set formulation to impose different contact angles at the solid-fluid-fluid contact line for two-phase flow in simple rhomboidal pore geometries, and calculate the maximum mean curvature (equivalently capillary pressure) for each case. We compare our results with a detailed set of analytical and experimental results in a range of pore geometries of varying wettability from Mason and Morrow (1994), and demonstrate the accuracy of this method. While the simulations shown are for relatively simple geometries, the method has the ability to handle arbitrarily complex geometry (such as input from X-ray microtomography imaging). / text

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第18781号 / 理博第4039号 / 新制||理||1581(附属図書館) / 31732 / 京都大学大学院理学研究科物理学・宇宙物理学専攻 / (主査)准教授 藤 定義, 教授 佐々 真一, 教授 早川 尚男 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM

Numerical simulation of cell motility is one of the difficult problems in computational science. It belongs to a class of problems which involve moving interfaces between flowing or deforming media. Different numerical techniques are being developed for different application areas and in this work an attempt is made to apply two popular numerical techniques used in the field of computational multiphase flows to a cell motility problem. An unsteady cell motility problem is considered to simulate numerically based on a two-dimensional mathematical model. Two important numerical methods, the Level set method and the Front tracking methods are applied to the cell motility problem to study several cases and to verify the convergence of the solution. With the assumption of no mechanical or physical obstructions to the cell, the results of the numerical simulations show that the domain shapes converge to a circular shape as they reach the steady state condition. The final steady state velocities with which the domains move and the final steady state area to which they converge are observed to be independent of domain shapes. Moreover all shapes converge to exactly same radius of circle and move with same velocity after reaching steady state condition.

In bubble dynamics, an underwater bubble may evolve from being singly-connected to being toroidal. Furthermore, two or more individual bubbles may merge to form a single large bubble. These dynamics involve significant topological changes such as merging and breaking, which may not be handled well by front-tracking boundary element methods. In the level set method, topological changes are handled naturally through a higher-dimensional level set function. This makes it an attractive method for bubble simulation. In this paper, we present a method that combines the level set method and the boundary element method for the simulation of bubble dynamics. We propose a formulation for the update of a potential function in the level set context. This potential function is non-physical off the bubble surface but consistent with the physics on the bubble surface. We consider only axisymmetric cavitation bubbles in this paper. Included in the paper are some preliminary results and findings. / Singapore-MIT Alliance (SMA)

In this dissertation, novel variational optical-flow and active-contour methods are investigated to address challenging problems in computer vision and shape-based imaging. Starting from traditional applications of these methods in computer vision, such as object segmentation, tracking, and detection, this research subsequently applies similar active contour techniques to the realm of shape-based imaging, which is an image reconstruction technique estimating object shapes directly from physical wave measurements. In particular, the first and second part of this thesis deal with the following two physically inspired computer vision applications. Optical Flow for Vision-Based Flame Detection: Fire motion is estimated using optimal mass transport optical flow, whose motion model is inspired by the physical law of mass conservation, a governing equation for fire dynamics. The estimated motion fields are used to first detect candidate regions characterized by high motion activity, which are then tracked over time using active contours. To classify candidate regions, a neural net is trained on a set of novel motion features, which are extracted from optical flow fields of candidate regions. Coupled Photo-Geometric Object Features: Active contour models for segmentation in thermal videos are presented, which generalize the well-known Mumford-Shah functional. The diffusive nature of heat processes in thermal imagery motivates the use of Mumford-Shah-type smooth approximations for the image radiance. Mumford-Shah's isotropic smoothness constraint is generalized to anisotropic diffusion in this dissertation, where the image gradient is decomposed into components parallel and perpendicular to level set curves describing the object's boundary contour. In a limiting case, this anisotropic Mumford-Shah segmentation energy yields a one-dimensional ``photo-geometric'' representation of an object which is invariant to translation, rotation and scale. These properties allow the photo-geometric object representation to be efficiently used as a radiance feature; a recognition-segmentation active contour energy, whose shape and radiance follow a training model obtained by principal component analysis of a training set's shape and radiance features, is finally applied to tracking problems in thermal imagery. The third part of this thesis investigates a physics-driven active contour approach for shape-based imaging. Adjoint Active Contours for Shape-Based Imaging: The goal of this research is to estimate both location and shape of buried objects from surface measurements of waves scattered from the object. These objects' shapes are described by active contours: A misfit energy quantifying the discrepancy between measured and simulated wave amplitudes is minimized with respect to object shape using the adjoint state method. The minimizing active contour evolution requires numerical forward scattering solutions, which are obtained by way of the method of fundamental solutions, a meshfree collocation method. In combination with active contours being implemented as level sets, one obtains a completely meshfree algorithm; a considerable advantage over previous work in this field. With future applications in medical and geophysical imaging in mind, the method is formulated for acoustic and elastodynamic wave processes in the frequency domain.

Pas de résumé / The level set method was introduced by Osher & Sethian (1988) as a general technique to capture moving interfaces. It has been used to study crystal growth, to simulate water and fire for computer graphics applications, to study two-phase flows and in many other fields. The wellknown problem of the level set method is the following: if the flow velocity is not constant, the level set scalar may become strongly distorted. Thus, the numerical integration may suffer from loss of accuracy. In level set methods, this problem is remedied by the reinitialization procedure, i.e. by reconstruction of the level set function in a way to satisfy the eikonal equation. We propose an alternative approach. We modify directly the level set equation by embedding a source term. The exact expression of this term is such that the eikonal equation is automatically satisfied. Furthermore on the interface, this term is equal to zero. In the meantime, the advantage of our approach is this: the exact expression of the source term allows for the possibility of derivation of its local approximate forms, of first-and-higher order accuracy. Compared to the extension velocity method, this may open the simplifications in realization of level set methods. Compared to the standard approach with the reinitialization procedure, this may give the economies in the number of level set re-initializations, and also, due to reduced number of reinitializations, one may expect an improvement in resolution of zero-set level. Hence, the objective of the present dissertation is to describe and to assess this approach in different test cases.

No description available.