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.

Optimal control (motion planning) of the free interface in classical two-phase Stefan problems is considered. The evolution of the free interface is modeled by a level set function. The first-order optimality system is derived on a formal basis. It provides gradient information based on the adjoint temperature and adjoint level set function. Suitable discretization schemes for the forward and adjoint systems are described. Numerical examples verify the correctness and flexibility of the proposed scheme.

A multi-layered periodic structure is investigated for optimal shape design in diffraction gratings. A periodic dielectric material is used as the scattering profile for a planar incident wave. Designing optimal profiles for scattering is a type of inverse problem. The ability to fabricate such materials on the order of the wavelength of the incoming light is key for design strategies. We compute a finite element approximation on a variational setup of the forward problem. On the inverse and optimal design problem, we discuss the stability of the designs and develop computational strategies based on a level-set evolutionary approach.

Resin infusion-based processes are good candidates for manufacturing thin composite materials parts such as those used in aeronautics for instance. These processes consist in infusing a liquid resin into a stacking of fibrous preforms under the action of a mechanical pressure field applied onto this stacking where a stiff- distribution medium is also placed to create a resin feeding. Both physical and mechanical properties of the final pieces are rather difficult to predict and control. Numerical simulation are perfectly suited to master these processes. In the present work a numerical finite element modeling framework is proposed to simulate infusion processes. The flow of the assumed Newtonian resin is described in the distribution medium, a highly porous medium, through Stokes' equations and through Darcy's equations in the fibrous preforms, very low permeability media. This coupled Stokes-Darcy flow is modeled in a monolithic approach which consists in using a single mesh for both media. The mixed velocity- pressure formulation is then discretized by linear-linear finite elements, stabilized by a so-called ASGS multi-scale approach. Both Stokes-Darcy interface and fluid front are represented individually thanks to "Level-Set" functions, and some specific coupling conditions are prescribed on the interface separating both fluid and porous media. During the process, orthotropic preforms undergo finite strains, either during the compaction stage when resin is not yet present, or during resin infusion. Resin pressure then tends to make the preforms swell. Preforms deformations are represented through an updated Lagrangian formulation for finite deformations. Dry preforms possess a non-linear elastic behaviour in their transverse direction - across their thickness- given by existing experimental measurements. The effect of the presence of resin in the wet preforms is accounted for using a Terzaghi's equivalent model. Also, when preforms deform their porosity will change, and so will their permeability, modifying the resin flow. A Carman-Kozeny model is then used to relate porosity and permeability. After the Stokes-Darcy coupling is validated both on numerous tests cases and using the method of manufactured solutions, various 2D and 3D simulations of injection and infusion-based processes are analyzed.The latter includ- ing preform deformations along with resin flow. Comparisons with existing experimental measurements permit to validate the approach on a simple geometry. Last, some ex- tensions to more complex 3D cases are proposed as outlooks, including curvatures and thickness variations.

Piracy is still a significant threat to ships in a maritime environment. Areas such as the coast of Somalia and the Strait of Malacca are still plagued by pirates, and the total international cost of piracy numbers in the billions of dollars. The first line of defence against these threats is early detection and thus maritime surveillance has become an increasingly important task over the years. While surveillance has traditionally been a manual task using crew members in lookout positions on parts of the ship, much work is being done to automate this task using digital cameras equipped with computer vision software. While these systems are beneficial in that they do not grow tired like their human counterparts, the maritime environment is a challenging task for computer vision systems. This dissertation aims to address some of these challenges by presenting a system that is able to use prior knowledge of an object’s shape to aid in detection and tracking of the object. Additionally, it aims to test this system under various environmental conditions (such as weather). The system is based around the segmentation technique known as the level set method, which uses a contour in the image that is evolved to separate regions of interest. The system is split into two parts, comprising of an object detection stage that initially finds objects in a scene, and an object tracking stage that tracks detected objects for the rest of the sequence. The object detection stage uses a kernel density estimation-based background subtraction and a binary image level set filter, while the object tracker makes use of a tracking level set algorithm for its functionality. The object detector was tested using a group of 4 sequences, of which it was able to find a prior-known object in 3. The object tracker was tested on a group of 10 sequences for 300 frames a sequence. In 6 of these sequences the object tracker was able to successfully track the object in every single frame. It is shown that the developed video tracking system outperforms level set–based systems that don’t use prior shape knowledge, working well even where these systems fail. / Thesis (M.Sc.Eng.)-University of KwaZulu-Natal, Durban, 2012.

Le soudage hybride arc / laser représente une solution adéquate à l'assemblage de tôles d'acier de forte épaisseur. La présence d'une source laser en amont de la torche MIG permet d'accroître la productivité du procédé tout en assurant une excellente qualité de la soudure. Cependant la phénoménologie complexe de ce procédé multiphysique n'est pas encore totalement maîtrisée, ce qui motive le développement d'outils de simulation numérique. La présente étude s'est déroulée dans le cadre d'un projet multipartenaires "SISHYFE" (ANR Matériaux et Procédés).Dans ce but un modèle éléments finis 3D non stationnaire a été développé. Construit à partir de l'approche level set, il est en mesure de simuler un procédé multi-passes de soudage hybride arc / laser. Il a été développé autour de quatre principaux axes. (1) Une modélisation plus réaliste de l'apport de chaleur a conduit à la définition d'un nouveau modèle basé sur le rayonnement thermique pour décrire la source MIG. (2) Grâce au couplage entre l'intégration de termes source de matière et de chaleur dans les équations de conservation et l'actualisation de l'interface gaz / métal, le modèle proposé est capable de simuler le développement d'un cordon de soudure. (3) L'intégration d'une modélisation des écoulements du bain de fusion dans une approche level set a été évaluée et son impact sur la distribution thermique dans le métal analysé. (4) La simulation de la formation des contraintes pendant et après soudage a été possible grâce à l'adaptation d'un solveur thermomécanique.La modélisation finale a permis de simuler une configuration industrielle de soudage hybride arc / laser multi-passes. Des mesures expérimentales effectuées par les partenaires du projet "SISHYFE" ont été utilisées afin d'évaluer le modèle et d'éprouver sa capacité à reproduire l'expérience.

Reservoir characterization is one of the most important problems in petroleum engineering. It involves forward reservoir modeling that predicts the fluid behavior in the reservoir and inverse problem that calibrates created reservoir models with given data. In this dissertation, we focus on two problems in the field of reservoir characterization: depth of investigation in heterogeneous reservoirs, and history matching and uncertainty quantification of channelized reservoirs. The concept of depth of investigation is fundamental to well test analysis. Much of the current well test analysis relies on analytical solutions based on homogeneous or layered reservoirs. However, such analytic solutions are severely limited for heterogeneous and fractured reservoirs, particularly for unconventional reservoirs with multistage hydraulic fractures. We first generalize the concept to heterogeneous reservoirs and provide an efficient tool to calculate drainage volume using fast marching methods and estimate pressure depletion based on geometric pressure approximation. The applicability of proposed method is illustrated using two applications in unconventional reservoirs including flow regime visualization and stimulated reservoir volume estimation. Due to high permeability contrast and non-Gaussianity of channelized permeability field, it is difficult to history match and quantify uncertainty of channelized reservoirs using traditional approaches. We treat facies boundaries as level set functions and solve the moving boundary problem (history matching) with the level set equation. In addition to level set methods, we also exploit the problem using pixel based approach. The reversible jump Markov Chain Monte Carlo approach is utilized to search the parameter space with flexible dimensions. Both proposed approaches are demonstrated with two and three dimensional examples.

Essen, University, Diss., 2005--Duisburg.

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