In the present study, a fifth-order weighted essentially non-oscillatory (WENO) scheme was built for solving the surface-capturing level-set equation. Combined with the level-set equation, the three-dimensional Reynolds averaged Navier-Stokes (RANS) equations were employed for the prediction of nonlinear wave-interaction and wave-breaking phenomena over sloping beaches. In the level-set finite-analytic Navier-Stokes (FANS) method, the free surface is represented by the zero level-set function, and the flows are modeled as immiscible air-water two phase flows. The Navier-Stokes equations for air-water two phase flows are formulated in a moving curvilinear coordinate system and discretized by a 12-point finite-analytical scheme using the finite-analytic method on a multi-block over-set grid system. The Pressure Implicit with Splitting of Operators / Semi-Implicit Method for Pressure-Linked Equation Revised (PISO/SIMPLER) algorithm was used to determine the coupled velocity and pressure fields. The evolution of the level-set method was solved using the third-order total variation diminishing (TVD) Runge-Kutta method and fifth-order WENO scheme. The accuracy was confirmed by solving the Zalesak's problem. Two major subjects are discussed in the present study. First, to identify the WENO scheme as a more accurate scheme than the essentially non-oscillatory scheme (ENO), the characteristics of a nonlinear monochromatic wave were studied systematically and comparisons of wave profiles using the two schemes were conducted. To eliminate other factors that might produce wave profile fluctuation, different damping functions and grid densities were studied. To damp the reflection waves efficiently, we compared five damping functions. The free-surface elevation data collected from gauges distributed evenly in a numerical wave tank are analyzed to demonstrate the damping effect of the beach. Second, as a surface-tracking numerical method built on curvilinear coordinates, the level-set RANS model was tested for nonlinear bichromatic wave trains and breaking waves on a sloping beach with a complex free surface. As the wave breaks, the velocity of the fluid flow surface became more complex. Numerical modeling was performed to simulate the two-phase flow velocity and its corresponding surface and evolution when the wave passed over different sloping beaches. The breaking wave test showed that it is an efficient technique for accurately capturing the breaking wave free surface. To predict the breaking points, different wave heights and beach slopes are simulated. The results show that the dependency of wave shape and breaking characteristics to wave height and beach slope match the results provided by experiments.

<p>Consider a mixture of two immiscible, incompressible fluids e.g. oil and water. Since the fluids do not mix, an interface between the two fluids will form and move in time. The motion of the two fluids can be modelled by the incompressible Navier-Stokes equations for two phase flow with surface tension together with a representation of the moving interface. The parameters in the Navier-Stokes equations will depend on the position and other properties of the interface. The interface should move with the velocity of the flow at the interface. Since the fluids are incompressible, the density of each fluid is constant. Mass conservation then implies that the volume occupied by each of the two fluids should not change with time. The object of this thesis has been to develop a new numerical method to simulate incompressible two phase flow accurately that conserves mass and volume of each fluid correctly.</p><p>Numerical simulations of incompressible two phase flow with surface tension have been a challenge for many years. Several methods have been developed and used prior to the work presented in this thesis. The two most commonly used methods are volume of fluid methods and level set methods. There are advantages and disadvantages of both of the methods.</p><p>In volume of fluid methods the interface is represented by a discontinuity of a globally defined function. Because of the discontinuity it is hard both to move the interface as well as to calculate properties of the interface such as curvature. Specially designed methods have to be used, and all these methods are low order accurate. Volume of fluid methods do however conserve the volumes of the two fluids correctly.</p><p>In level set methods the interface is represented by the zero contour of the globally defined signed distance function. This function is smooth across the interface. Since the function is smooth, standard methods for partial differential equations can be used to advect the interface accurately. A reinitialization is however needed to make sure that the level set function remains a signed distance function. During this process the zero contour might move slightly. Because of this, the volume conservation of the method becomes poor.</p><p>In this thesis we present a new level set method. The method is designed such that the volume of each fluid is conserved, at least approximately. The interface is represented by the 0.5 contour of a regularized characteristic function. As for standard level set methods, the interface is moved first by an advective step, and then reinitialized. Unlike traditional level set methods, we can formulate the reinitialization as a conservation law. Conservative methods can then be used to move and to reinitialize the level set function numerically. Since the level set function is a regularized characteristic function, we can expect good conservation of the volume bounded by the interface.</p><p>The method is discretized using both finite differences and finite elements. Uniform and adaptive grids are used in both two and three space dimensions. Good convergence as well as volume conservation is observed. Theoretical studies are performed to investigate the conservation and the computational time needed for reinitialization.</p>

Extreme hydrological events associated with global warming are likely to produce an increasing number of flooding scenarios resulting in significant bridge inundation and associated damages. During large floods, the presence of a bridge in an open channel triggers a highly turbulent flow field including 3D complex coherent structures around bridge structures. These turbulence structures are highly energetic and possess high sediment entrainment capacity which increases scouring around the bridge foundation and consequently lead to structural stability problems or even failure of the structure. Hence, understanding the complex turbulent flow field for these extreme flow conditions is crucial to estimate the failure risks for existing bridges and better design of future bridges. This research employs the method Large Eddy Simulation (LES) to predict accurately the 3D turbulent flow around bridge structures. The LES code is refined with a novel free surface algorithm based on the Level Set Method (LSM) to determine the complex water surface profiles. The code is used to analyze the hydrodynamics of compound channel flow with deep and shallow overbanks, free flow around a bridge abutment, pressure flow with a partially submerged bridge deck and bridge overtopping flow. All simulations are validated with data from complementary physical model tests under analogous geometrical and flow conditions. Primary velocity, bed shear stress, turbulence characteristics and 3D coherent flow structures are examined thoroughly for each of the flow cases to explain the hydrodynamics of these complex turbulent flows.

Vibration control is inevitable in many fields, including mechanical and civil engineering. This matter becomes more crucial for lightweight systems, like those made of magnesium. One of the most commonly practiced methods in vibration control is to apply constrained layer damping (CLD) patches to the surface of a structure. In order to consider the weight efficiency of the structure, the best shape and locations of the patches should be determined to achieve the optimum vibration suppression with the lowest amount of damping patch. In most research work done so far, the shape of patches are assumed to be known and only their optimum locations are found. However, the shape of the patches plays an important role in vibration suppression that should be included in the overall optimization procedure. In this research, a novel topology optimization approach is proposed. This approach is capable of finding the optimum shape and locations of the patches simultaneously for a given surface area. In other words, the damping optimization will be formulated in the context of the level set technique, which is a numerical method used to track shapes and locations concurrently. Although level set technique offers several key benefits, its application especially in time-varying problems is somewhat cumbersome. To overcome this issue, a unique programming technique is suggested that utilizes MATLAB© and COMSOL© simultaneously. Different 2D structures will be considered and CLD patches will be optimally located on them to achieve the highest modal loss factor. Optimization will be performed while having different amount of damping patches to check the effectiveness of the technique. In all cases, certain constraints are imposed in order to make sure that the amount of damping material remains constant and equal to the starting value. Furthermore, different natural frequencies will be targeted in the damping optimization, and their effects will also be explained. The level set optimization technique will then be expanded to 3D structures, and a novel approach will be presented for defining an efficient 4D level set function to initialize the optimization process. Vibrations of a satellite dish will be optimally suppressed using CLD patches. Dependency of the optimum shape and location of patches to different parameters of the models such as natural frequencies and initial starting point will be examined. In another practical example, excessive vibrations of an automotive dash panel will be minimized by adding damping materials and their optimal distribution will be found. Finally, the accuracy of the proposed method will be experimentally confirmed through lab tests on a rectangular plate with nonsymmetrical boundary conditions. Different damping configurations, including the optimum one, will be tested. It will be shown that the optimum damping configuration found via level set technique possesses the highest loss factor and reveals the best vibration attenuation. The proposed level set topology optimization method shows high capability of determining the optimum damping set in structures. The effective coding method presented in this research will make it possible to easily extend this method to other physical problems such as image processing, heat transfer, magnetic fields, etc. Being interconnected, the physical part will be modeled in a finite element package like COMSOL and the optimization advances by means of Hamilton-Jacobi partial differential equation. Thus, the application of the proposed method is not confined to damping optimization and can be expanded to many engineering problems. In summary, this research: - offers general solution to 2D and 3D CLD applications and simultaneously finds the best shape and location of the patches for a given surface area (damping material); - extends the level set technique to concurrent shape and location optimization; - proposes a new numerical implementation to handle level set optimization problems in any complicated structure; - makes it possible to perform level set optimization in time dependent problems; - extends level set approach to higher order problems.

Changes in retinal structure are related to different eye diseases. Various retinal imaging techniques, such as fundus imaging and optical coherence tomography (OCT) imaging modalities, have been developed for non-intrusive ophthalmology diagnoses according to the vasculature changes. However, it is time consuming or even impossible for ophthalmologists to manually label all the retinal structures from fundus images and OCT images. Therefore, computer aided diagnosis system for retinal imaging plays an important role in the assessment of ophthalmologic diseases and cardiovascular disorders. The aim of this PhD thesis is to develop segmentation methods to extract clinically useful information from these retinal images, which are acquired from different imaging modalities. In other words, we built the segmentation methods to extract important structures from both 2D fundus images and 3D OCT images. In the first part of my PhD project, two novel level set based methods were proposed for detecting the blood vessels and optic discs from fundus images. The first one integrates Chan-Vese's energy minimizing active contour method with the edge constraint term and Gaussian Mixture Model based term for blood vessels segmentation, while the second method combines the edge constraint term, the distance regularisation term and the shape-prior term for locating the optic disc. Both methods include the pre-processing stage, used for removing noise and enhancing the contrast between the object and the background. Three automated layer segmentation methods were built for segmenting intra-retinal layers from 3D OCT macular and optic nerve head images in the second part of my PhD project. The first two methods combine different methods according to the data characteristics. First, eight boundaries of the intra-retinal layers were detected from the 3D OCT macular images and the thickness maps of the seven layers were produced. Second, four boundaries of the intra-retinal layers were located from 3D optic nerve head images and the thickness maps of the Retinal Nerve Fiber Layer (RNFL) were plotted. Finally, the choroidal layer segmentation method based on the Level Set framework was designed, which embedded with the distance regularisation term, edge constraint term and Markov Random Field modelled region term. The thickness map of the choroidal layer was calculated and shown.

abstract: Many physical phenomena and industrial applications involve multiphase fluid flows and hence it is of high importance to be able to simulate various aspects of these flows accurately. The Dynamic Contact Angles (DCA) and the contact lines at the wall boundaries are a couple of such important aspects. In the past few decades, many mathematical models were developed for predicting the contact angles of the inter-face with the wall boundary under various flow conditions. These models are used to incorporate the physics of DCA and contact line motion in numerical simulations using various interface capturing/tracking techniques. In the current thesis, a simple approach to incorporate the static and dynamic contact angle boundary conditions using the level set method is developed and implemented in multiphase CFD codes, LIT (Level set Interface Tracking) (Herrmann (2008)) and NGA (flow solver) (Desjardins et al (2008)). Various DCA models and associated boundary conditions are reviewed. In addition, numerical aspects such as the occurrence of a stress singularity at the contact lines and grid convergence of macroscopic interface shape are dealt with in the context of the level set approach. / Dissertation/Thesis / Masters Thesis Mechanical Engineering 2015

Distributed modal filters based on piezoelectric polymer have especially become popular in the field of active vibration control to reduce the problem of spillover. While distributed modal filters for one-dimensional structures can be found analytically based on the orthogonality between the mode shapes, the design for two-dimensional structures is not straightforward. It requires a continuous gain variation in two dimensions, which is not realizable from the current manufacturing point of view. In this thesis, a structural optimization problem is considered to approximate distributed modal sensors for two-dimensional plate structures, where the thickness is constant but the polarization can switch between positive and negative. The problem is solved through an explicit parametric level set method. In this framework, the boundary of a domain is represented implicitly by the zero isoline of a level set function. This allows simultaneous shape and topology changes. The level set function is approximated by a linear combination of Gaussian radial basis functions. As a result, the structural optimization problem can be directly posed in terms of the parameters of the approximation. This allows to apply standard optimization methods and bypasses the numerical drawbacks, such as reinitialization, velocity extension and regularization, which are associated with the numerical solution of the Hamilton-Jacobi equation in conventional methods.Since the level set method based on the shape derivative formally only allows shape but not topology transformation, the optimization problem is firstly tackled with a derivative-free optimization algorithm. It is shown that the approach is able to find approximate modal sensor designs with only few design variables. However, this approach becomes unsuitable as soon as the number of optimization variables is growing. Therefore, a sensitivity-based optimization approach is being applied, based on the parametric shape derivative which is with respect to the parameters of the radial basis functions. Although the shape derivatives does not exist at points where the topology changes, it is demonstrated that an optimization routine based on a SQP solver is able to perform topological changes during the optimization and finds optimal designs even from poor initial designs. In order to include the sensors' distribution as design variable, the parametric level set approach is extended to multiple level sets. It turns out that, despite the increased design space, optimal solutions always converge to full-material polarization designs. Numerical examples are provided for a simply supported as well as a cantilever square plate. / Doctorat en Sciences de l'ingénieur et technologie / info:eu-repo/semantics/nonPublished

Consider a mixture of two immiscible, incompressible fluids e.g. oil and water. Since the fluids do not mix, an interface between the two fluids will form and move in time. The motion of the two fluids can be modelled by the incompressible Navier-Stokes equations for two phase flow with surface tension together with a representation of the moving interface. The parameters in the Navier-Stokes equations will depend on the position and other properties of the interface. The interface should move with the velocity of the flow at the interface. Since the fluids are incompressible, the density of each fluid is constant. Mass conservation then implies that the volume occupied by each of the two fluids should not change with time. The object of this thesis has been to develop a new numerical method to simulate incompressible two phase flow accurately that conserves mass and volume of each fluid correctly. Numerical simulations of incompressible two phase flow with surface tension have been a challenge for many years. Several methods have been developed and used prior to the work presented in this thesis. The two most commonly used methods are volume of fluid methods and level set methods. There are advantages and disadvantages of both of the methods. In volume of fluid methods the interface is represented by a discontinuity of a globally defined function. Because of the discontinuity it is hard both to move the interface as well as to calculate properties of the interface such as curvature. Specially designed methods have to be used, and all these methods are low order accurate. Volume of fluid methods do however conserve the volumes of the two fluids correctly. In level set methods the interface is represented by the zero contour of the globally defined signed distance function. This function is smooth across the interface. Since the function is smooth, standard methods for partial differential equations can be used to advect the interface accurately. A reinitialization is however needed to make sure that the level set function remains a signed distance function. During this process the zero contour might move slightly. Because of this, the volume conservation of the method becomes poor. In this thesis we present a new level set method. The method is designed such that the volume of each fluid is conserved, at least approximately. The interface is represented by the 0.5 contour of a regularized characteristic function. As for standard level set methods, the interface is moved first by an advective step, and then reinitialized. Unlike traditional level set methods, we can formulate the reinitialization as a conservation law. Conservative methods can then be used to move and to reinitialize the level set function numerically. Since the level set function is a regularized characteristic function, we can expect good conservation of the volume bounded by the interface. The method is discretized using both finite differences and finite elements. Uniform and adaptive grids are used in both two and three space dimensions. Good convergence as well as volume conservation is observed. Theoretical studies are performed to investigate the conservation and the computational time needed for reinitialization. / QC 20101122

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

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.