In this work, we develop an interface-preserving level-set method in the finite-element framework for interfacial flows with moving contact lines. In our method, the contact line is advected naturally by the flow field. Contact angle hysteresis can be easily implemented without explicit calculation of the contact angle or the contact line velocity, and meshindependent results can be obtained following a simple computational strategy. We have implemented the method in three dimensions and provide numerical studies that compare well with analytical solutions to verify our algorithm. We first develop a high-order numerical method for interface-preserving level-set reinitialization. Within the interface cells, the gradient of the level set function is determined by a weighted local projection scheme and the missing additive constant is determined such that the position of the zero level set is preserved. For the non-interface cells, we compute the gradient of the level set function by solving a Hamilton-Jacobi equation as a conservation law system using the discontinuous Galerkin method. This follows the work by Hu and Shu [SIAM J. Sci. Comput. 21 (1999) 660-690]. The missing constant for these cells is recovered using the continuity of the level set function while taking into account the characteristics. To treat highly distorted initial conditions, we develop a hybrid numerical flux that combines the Lax-Friedrichs flux and a penalty flux. Our method is accurate for non-trivial test cases and handles singularities away from the interface very well. When derivative singularities are present on the interface, a second-derivative limiter is designed to suppress the oscillations. At least (N + 1)th order accuracy in the interface cells and Nth order accuracy in the whole domain are observed for smooth solutions when Nth degree polynomials are used. Two dimensional test cases are presented to demonstrate superior properties such as accuracy, long-term stability, interface-preserving capability, and easy treatment of contact lines. We then develop a level-set method in the finite-element framework. The contact line singularity is removed by the slip boundary condition proposed by Ren and E [Phys. Fluids, vol. 19, p. 022101, 2007], which has two friction coefficients: βN that controls the slip between the bulk fluids and the solid wall and βCL that controls the deviation of the microscopic dynamic contact angle from the static one. The predicted contact line dynamics from our method matches the Cox theory very well. We further find that the same slip length in the Cox theory can be reproduced by different combinations of (βN; βCL). This combination leads to a computational strategy for mesh-independent results that can match the experiments. There is no need to impose the contact angle condition geometrically, and the dynamic contact angle automatically emerges as part of the numerical solution. With a little modification, our method can also be used to compute contact angle hysteresis, where the tendency of contact line motion is readily available from the level-set function. Different test cases, including code validation and mesh-convergence study, are provided to demonstrate the efficiency and capability of our method. Lastly, we extend our method to three-dimensional simulations, where an extension equation is solved on the wall boundary to obtain the boundary condition for level-set reinitializaiton with contact lines. Reinitialization of ellipsoidal interfaces is presented to show the accuracy and stability of our method. In addition, simulations of a drop on an inclined wall are presented that are in agreement with theoretical results. / Doctor of Philosophy / When a liquid droplet is sliding along a solid surface, a moving contact line is formed at the intersection of the three phases: liquid, air and solid. This work develops a numerical method to study problems with moving contact lines. The partial differential equations describing the problem are solved by finite element methods. Our numerical method is validated against experiments and theories. Furthermore, we have implemented our method in three-dimensional problems.

This thesis is concerned with motion planning for the classical two-phase Stefan problem in level set formulation. The interface separating the fluid phases from the solid phases is represented as the zero level set of a continuous function whose evolution is described by the level set equation. Heat conduction in the two phases is modeled by the heat equation. A quadratic tracking-type cost functional that incorporates temperature tracking terms and a control cost term that expresses the desire to have the interface follow a prescribed trajectory by adjusting the heat flux through part of the boundary of the computational domain. The formal Lagrange approach is used to establish a first-order optimality system by applying shape calculus tools. For the numerical solution, the level set equation and its adjoint are discretized in space by discontinuous Galerkin methods that are combined with suitable explicit Runge-Kutta time stepping schemes, while the temperature and its adjoint are approximated in space by the extended finite element method (which accounts for the weak discontinuity of the temperature by a dynamic local modification of the underlying finite element spaces) combined with the implicit Euler method for the temporal discretization. The curvature of the interface which arises in the adjoint system is discretized by a finite element method as well. The projected gradient method, and, in the absence of control constraints, the limited memory BFGS method are used to solve the arising optimization problems. Several numerical examples highlight the potential of the proposed optimal control approach. In particular, they show that it inherits the geometric flexibility of the level set method. Thus, in addition to unidirectional solidification, closed interfaces and changes of topology can be tracked. Finally, the Moreau-Yosida regularization is applied to transform a state constraint on the position of the interface into a penalty term that is added to the cost functional. The optimality conditions for this penalized optimal control problem and its numerical solution are discussed. An example confirms the efficacy of the state constraint. / Die vorliegende Arbeit beschäftigt sich mit einem Optimalsteuerungsproblem für das klassische Stefan-Problem in zwei Phasen. Die Phasengrenze wird als Niveaulinie einer stetigen Funktion modelliert, was die Lösung der so genannten Level-Set-Gleichung erfordert. Durch Anpassen des Wärmeflusses am Rand des betrachteten Gebiets soll ein gewünschter Verlauf der Phasengrenze angesteuert werden. Zusammen mit dem Wunsch, ein vorgegebenes Temperaturprofil zu approximieren, wird dieses Ziel in einem quadratischen Zielfunktional formuliert. Die notwendigen Optimalitätsbedingungen erster Ordnung werden formal mit Hilfe der entsprechenden Lagrange-Funktion und unter Benutzung von Techniken aus der Formoptimierung hergeleitet. Für die numerische Lösung müssen die auftretenden partiellen Differentialgleichungen diskretisiert werden. Dies geschieht im Falle der Level-Set-Gleichung und ihrer Adjungierten auf Basis von unstetigen Galerkin-Verfahren und expliziten Runge-Kutta-Methoden. Die Wärmeleitungsgleichung und die entsprechende Gleichung im adjungierten System werden mit einer erweiterten Finite-Elemente-Methode im Ort sowie dem impliziten Euler-Verfahren in der Zeit diskretisiert. Dieser Zugang umgeht die aufwändige Adaption des Gitters, die normalerweise bei der FE-Diskretisierung von Phasenübergangsproblemen unvermeidbar ist. Auch die Krümmung der Phasengrenze wird numerisch mit Hilfe der Methode der finiten Elemente angenähert. Zur Lösung der auftretenden Optimierungsprobleme werden ein Gradienten-Projektionsverfahren und, im Fall dass keine Kontrollschranken vorliegen, die BFGS-Methode mit beschränktem Speicherbedarf eingesetzt. Numerische Beispiele beleuchten die Stärken des vorgeschlagenen Zugangs. Es stellt sich insbesondere heraus, dass sich die geometrische Flexibilität der Level-Set-Methode auf den vorgeschlagenen Zugang zur optimalen Steuerung vererbt. Zusätzlich zur gerichteten Bewegung einer flachen Phasengrenze können somit auch geschlossene Phasengrenzen sowie topologische Veränderungen angesteuert werden. Exemplarisch, und zwar an Hand einer Beschränkung an die Lage der Phasengrenze, wird auch noch die Behandlung von Zustandsbeschränkungen mittels der Moreau-Yosida-Regularisierung diskutiert. Ein numerisches Beispiel demonstriert die Wirkung der Zustandsbeschränkung.

This thesis is concerned with motion planning for the classical two-phase Stefan problem in level set formulation. The interface separating the fluid phases from the solid phases is represented as the zero level set of a continuous function whose evolution is described by the level set equation. Heat conduction in the two phases is modeled by the heat equation. A quadratic tracking-type cost functional that incorporates temperature tracking terms and a control cost term that expresses the desire to have the interface follow a prescribed trajectory by adjusting the heat flux through part of the boundary of the computational domain. The formal Lagrange approach is used to establish a first-order optimality system by applying shape calculus tools. For the numerical solution, the level set equation and its adjoint are discretized in space by discontinuous Galerkin methods that are combined with suitable explicit Runge-Kutta time stepping schemes, while the temperature and its adjoint are approximated in space by the extended finite element method (which accounts for the weak discontinuity of the temperature by a dynamic local modification of the underlying finite element spaces) combined with the implicit Euler method for the temporal discretization. The curvature of the interface which arises in the adjoint system is discretized by a finite element method as well. The projected gradient method, and, in the absence of control constraints, the limited memory BFGS method are used to solve the arising optimization problems. Several numerical examples highlight the potential of the proposed optimal control approach. In particular, they show that it inherits the geometric flexibility of the level set method. Thus, in addition to unidirectional solidification, closed interfaces and changes of topology can be tracked. Finally, the Moreau-Yosida regularization is applied to transform a state constraint on the position of the interface into a penalty term that is added to the cost functional. The optimality conditions for this penalized optimal control problem and its numerical solution are discussed. An example confirms the efficacy of the state constraint. / Die vorliegende Arbeit beschäftigt sich mit einem Optimalsteuerungsproblem für das klassische Stefan-Problem in zwei Phasen. Die Phasengrenze wird als Niveaulinie einer stetigen Funktion modelliert, was die Lösung der so genannten Level-Set-Gleichung erfordert. Durch Anpassen des Wärmeflusses am Rand des betrachteten Gebiets soll ein gewünschter Verlauf der Phasengrenze angesteuert werden. Zusammen mit dem Wunsch, ein vorgegebenes Temperaturprofil zu approximieren, wird dieses Ziel in einem quadratischen Zielfunktional formuliert. Die notwendigen Optimalitätsbedingungen erster Ordnung werden formal mit Hilfe der entsprechenden Lagrange-Funktion und unter Benutzung von Techniken aus der Formoptimierung hergeleitet. Für die numerische Lösung müssen die auftretenden partiellen Differentialgleichungen diskretisiert werden. Dies geschieht im Falle der Level-Set-Gleichung und ihrer Adjungierten auf Basis von unstetigen Galerkin-Verfahren und expliziten Runge-Kutta-Methoden. Die Wärmeleitungsgleichung und die entsprechende Gleichung im adjungierten System werden mit einer erweiterten Finite-Elemente-Methode im Ort sowie dem impliziten Euler-Verfahren in der Zeit diskretisiert. Dieser Zugang umgeht die aufwändige Adaption des Gitters, die normalerweise bei der FE-Diskretisierung von Phasenübergangsproblemen unvermeidbar ist. Auch die Krümmung der Phasengrenze wird numerisch mit Hilfe der Methode der finiten Elemente angenähert. Zur Lösung der auftretenden Optimierungsprobleme werden ein Gradienten-Projektionsverfahren und, im Fall dass keine Kontrollschranken vorliegen, die BFGS-Methode mit beschränktem Speicherbedarf eingesetzt. Numerische Beispiele beleuchten die Stärken des vorgeschlagenen Zugangs. Es stellt sich insbesondere heraus, dass sich die geometrische Flexibilität der Level-Set-Methode auf den vorgeschlagenen Zugang zur optimalen Steuerung vererbt. Zusätzlich zur gerichteten Bewegung einer flachen Phasengrenze können somit auch geschlossene Phasengrenzen sowie topologische Veränderungen angesteuert werden. Exemplarisch, und zwar an Hand einer Beschränkung an die Lage der Phasengrenze, wird auch noch die Behandlung von Zustandsbeschränkungen mittels der Moreau-Yosida-Regularisierung diskutiert. Ein numerisches Beispiel demonstriert die Wirkung der Zustandsbeschränkung.

Dans cette thèse, nous calculons la forme d’un objet immergé d’aire donnée qui minimise la résistance de vague. Le corps, considéré lisse, avance à vitesse constante sous la surface libre d’un fluide qui est supposé parfait et incompressible. La résistance de vague est la traînée, c’est-à-dire la composante horizontale de la force exercée par le fluide sur l’obstacle. Nous utilisons les équations de Neumann-Kelvin 2D, qui s’obtiennent en linéarisant les équations d’Euler irrotationnelles avec surface libre. Le problème de Neumann-Kelvin est formulé comme une équation intégrale de frontière basée sur une solution fondamentale qui intègre la condition linéarisée à la surface libre. Nous utilisons une méthode de descente de gradient pour trouver un minimiseur local du problème de résistance de vague. Un gradient par rapport à la forme est calculé par la méthode de variation de frontières. Nous utilisons une approche level-set pour calculer la résistance de vague et gérer les déplacements de la frontière de l’obstacle. Nous obtenons une grande variété de formes optimales selon la profondeur de l’objet et sa vitesse. / In this thesis, we compute the shape of a fully immersed object with a given area which minimises the wave resistance. The smooth body moves at a constant speed under the free surface of a fluid which is assumed to be inviscid and incompressible. The wave resistance is the drag, i.e. the horizontal component of the force exerted by the fluid on the obstacle. We work with the 2D Neumann-Kelvin equations, which are obtained by linearising the irrotational Euler equations with a free surface. The Neumann-Kelvin problem is formulated as a boundary integral equation based on a fundamental solution which handles the linearised free surface condition. We use a gradient descent method to find a local minimiser of the wave resistance problem. A gradient with respect to the shape is calculated by a boundary variation method. We use a level-set approach to calculate the wave-making resistance and to deal with the displacements of the boundary of the obstacle. We obtain a great variety of optimal shapes depending on the depth of the object and its velocity.

Image segmentation is a process to identify regions of interest from digital images. Image segmentation plays an important role in medical image processing which enables a variety of clinical applications. It is also a tool to facilitate the detection of abnormalities such as cancerous lesions in the brain. Although numerous efforts in recent years have advanced this technique, no single approach solves the problem of segmentation for the large variety of image modalities existing today. Consequently, brain MRI segmentation remains a challenging task. The purpose of this thesis is to demonstrate brain MRI segmentation for delineation of tumors, ventricles and other anatomical structures using Insight Segmentation and Registration Toolkit (ITK) routines as the foundation. ITK is an open-source software system to support the Visible Human Project. Visible Human Project is the creation of complete, anatomically detailed, three-dimensional representations of the normal male and female human bodies. Currently under active development, ITK employs leading-edge segmentation and registration algorithms in two, three, and more dimensions. A goal of this thesis is to implement those algorithms to facilitate brain segmentation for a brain cancer research scientist.

O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de problemas de interação fluido-estrutura (IFE) em aplicações envolvendo corpos rígidos submersos, onde escoamentos incompressíveis de fluidos Newtonianos com superfície livre são considerados. Propõe-se o emprego da Análise Isogeométrica por permitir a unificação entre os procedimentos de pré-processamento e análise, melhorando assim as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação da conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, em uma formulação cinemática ALE (Arbitrary Lagrangean- Eulerian). A consideração da superfície livre no escoamento se dá tratando o fluido como um meio bifásico, através do método Level Set. O corpo rígido apresenta não linearidade na rotação e restrições representadas por vínculos elásticos e amortecedores viscosos, sendo a equação de equilíbrio dinâmico resolvida através do método de Newmark. O esquema de acoplamento sólido-fluido adotado é o particionado convencional, que impõe condições de compatibilidade cinemáticas e de equilíbrio sobre a interface sólido-fluido, analisando ambos os meios de maneira sequencial. A discretização das equações governantes é realizada através do esquema explícito de dois passos de Taylor-Galerkin, aplicado no contexto da Análise Isogeométrica. Por fim, são analisados alguns problemas da Dinâmica de Fluidos Computacional, de onde se concluiu que os resultados obtidos são bastante consistentes com os fenômenos envolvidos, com as ferramentas exclusivas da Análise Isogeométrica, como o refinamento k, melhorando a convergência dos resultados. Para escoamentos bifásicos, verificou-se que o método Level Set obteve resultados bastante promissores apresentando, entretanto, uma dissipação numérica excessiva. Propõe-se, para estudos futuros, a elaboração de esquemas numéricos que conservem melhor o volume da fase líquida do escoamento. / The present work aims to development of a numerical formulation based on Isogeometric Analysis for the study of Fluid-Structure Interaction problems in applications involving rigid bodies submerged, considering incompressible Newtonian flows with free surface. The use of the Isogeometric Analysis allows unification between the preprocessing and analysis steps, improving then the continuity of the base functions employed, both in the spatial discretization and approximation of the variables in the system of equations. The fundamental flow equations are formed by the Navier-Stokes and the mass conservation, described by de pseudo-compressibility hypothesis, in an ALE (Arbitrary Lagrangean-Eulerian) kinematic formulation. The free surface consideration of the flow is handled treating the fluid like a two- phase medium, using the Level Set method. The rigid body considers nonlinearity in rotation, and restrictions represented by elastic springs and viscous dampers, with the dynamic equilibrium equation being resolved using the Newmark’s method. The solid-fluid coupling scheme is the conventional partitioned, which imposes kinematics and equilibrium compatibility conditions on the solid-fluid interface, analyzing both mediums in a sequential manner. The governing equations are discretized using the explicit two step Taylor-Galerkin method, applied in an Isogeometric Analisys context. Finally, some Computational Fluid Dinamics problems are analysed, from which it was concluded that the results obtained are quite consistent with phenomena involved, with the unique tools of Isogeometric Analysis, such as k-refinement, improving the convergence of the results. For biphasic flows, it was verified that the Level Set method obtained very promising results, presenting, however, an excessive numerical dissipation. For future studies, it is proposed the elaboration of numerical schemes that better preserve the volume of the liquid phase of the flow.

La segmentation d'images est un domaine important dans le traitement d'images et un grand nombre d'approches différentes ent été développées pendant ces dernières décennies. L'approche des contours actifs est un des plus populaires. Dans ce cadre, cette thèse vise à développer des algorithmes robustes, qui peuvent segmenter des images avec des inhomogénéités d'intensité. Nous nous concentrons sur l'étude des énergies externes basées région dans le cadre des ensembles de niveaux. Précisément, nous abordons la difficulté de choisir l'échelle de la fenêtre spatiale qui définit la localité. Notre contribution principale est d'avoir proposé une échelle adaptative pour les méthodes de segmentation basées sur les statistiques locales. Nous utilisons l'approche d'Intersection des Intervalles de Confiance pour définir une échelle position-dépendante pour l'estimation des statistiques image. L'échelle est optimale dans le sens où elle donne le meilleur compromis entre le biais et la variance de l'approximation polynomiale locale de l'image observée conditionnellement à la segmentation actuelle. De plus, pour le model de segmentation basé sur une interprétation Bahésienne avec deux noyaux locaux, nous suggérons de considérer leurs valeurs séparément. Notre proposition donne une segmentation plus lisse avec moins de délocalisations que la méthode originale. Des expériences comparatives de notre proposition à d'autres méthodes de segmentation basées sur des statistiques locales sont effectuées. Les résultats quantitatifs réalisés sur des images ultrasonores de simulation, montrent que la méthode proposée est plus robuste au phénomène d'atténuation. Des expériences sur des images réelles montrent également l'utilité de notre approche. / Image segmentation is an important research area in image processing and a large number of different approaches have been developed over the last few decades. The active contour approach is one of the most popular among them. Within this framework, this thesis aims at developing robust algorithms, which can segment images with intensity inhomogeneities. We focus on the study of region-based external energies within the level set framework. We study the use of local image statistics for the design of external energies. Precisely, we address the difficulty of choosing the scale of the spatial window that defines locality. Our main contribution is to propose an adaptive scale for local region-based segmen¬tation methods. We use the Intersection of Confidence Intervals approach to define this pixel-dependent scale for the estimation of local image statistics. The scale is optimal in the sense that it gives the best trade-off between the bias and the variance of a Local Polynomial Approximation of the observed image conditional on the current segmenta¬tion. Additionally, for the segmentation model based on a Bayesian interpretation with two local kernels, we suggest to consider their values separately. Our proposition gives a smoother segmentation with less mis-localisations Chan the original method.Comparative experiments of our method to other local region-based segmentation me¬thods are carried out. The quantitative results, on simulated ultrasound B-mode images, show that the proposed scale selection strategy gives a robust solution to the intensity inhomogeneity artifact of this imaging modality. More general experiments on real images also demonstrate the usefulness of our approach.

Ce travail de thèse est consacré à deux thèmes de recherche en Calcul Scientifique liés par l'approximation numérique de problèmes en mécanique des fluides. Le premier thème concerne l'approximation numérique des équations de Stokes, modélisant les écoulements de fluides incompressibles à vitesse faible. Ce thème est présent dans plusieurs travaux en Calcul Scientifique. La discrétisation en temps est réalisée à l'aide de la méthode de projection. La discrétisation en espace utilise la méthode des éléments finis mixtes hybrides qui permet d'imposer de façon exacte la contrainte d'incompressibilité. Cette approche est originale: la méthode des éléments mixtes hybrides est couplée avec une méthode d'éléments finis standards. L'ordre de convergence des deux méthodes est préservé. Le second thème concerne la mise au point de méthodes numériques de type volumes finis pour la résolution de l'équation Level Set. Ces équations interviennent de manière essentielle dans la résolution des problèmes de propagation d'interfaces. Dans cette partie, nous avons développé une nouvelle méthode d'ordre 2 de type MUSCL pour résoudre le système hyperbolique résultant de l'équation Level Set. Nous illustrons ces propriétés par des applications numériques. En particulier nous avons regardé le cas du problème des deux demi-plans pour lequel notre schéma donne une approximation pour le gradient de la fonction Level Set. Par ailleurs, l'ordre de précision attendu est obtenu avec les normes $L_1$ et $L_{\infty}$ pour des fonctions régulières. Pour finir, il est à noter que notre méthode peut être facilement étendue aux problèmes d'Hamilton-Jacobi du premier et du second ordre.

L'étude des mécanismes de formation d'un spray à partir d'un écoulement liquide-gaz est un sujet de recherche fondamental, et dont les applications industrielles sont nombreuses, notamment dans le domaine de la combustion et de la propulsion. La simulation numérique de tels écoulements apparaît comme un complément indispensable aux expériences et aux théories pour aboutir à une compréhension et à une prédiction fine de tels phénomènes physiques. Durant cette étude nous avons développé une technique numérique de suivi d'interface couplée à un solveur des équations de Navier-Stokes afin d'étudier précisément la dynamique de l'interface qui sépare la phase liquide de la phase gazeuse. Nous décrivons en détail la méthode Level Set qui a été utilisée durant cette étude pour suivre le mouvement de l'interface, ainsi que les techniques numériques utilisées pour la résolution des équations de Navier-Stokes. Notamment différents schémas numériques sont testés pour que la précision du code de calcul développé soit optimal. Le couplage entre les équations du mouvement et le suivi d'interface peut être effectué par différentes techniques. L'utilisation de la Ghost Fluid Method pour réaliser ce couplage permet de préserver numériquement le caractère discontinu de l'interface liquide-gaz tout en proposant une discrétisation des équations précises et robustes. Les codes développés (2D, 3D parallélisation MPI) sont ensuite utilisés pour étudier des collisions de gouttes. Les comparaisons entre les résultats expérimentaux et les simulations numériques sont très probantes et mettent en valeur la prédictibilité des calculs effectués. D'autres configurations plus complexes sont étudiées dans le cadre d'une étude de faisabilité. L'étude de l'atomisation primaire d'un jet liquide continu en un spray de gouttelettes semble notamment être un domaine d'investigation particulièrement prometteur pour ce type de simulations. Enfin l'étude d'interfaces dites « réactives » comme la vaporisation d'un liquide et la propagation de flammes prémalangées a également été abordé en utilisant la Ghost Fluid Method pour imposer les conditions de saut spécifiques à ce type de phénomènes.

L'objet de cette thèse a été d'étudier la simulation numérique directe d'écoulements diphasiques de fluides non miscibles, incompressibles et isothermes avec la technique du raffinement adaptatif local de maillage. La résolution des équations de Navier-Stokes incompressibles est faite par le biais d'une méthode de projection explicite. La capture de l'interface est réalisée explicitement par la méthode Level Set. L'utilisation de la méthode Ghost Fluid pour le traitement des conditions de saut permet la résolution couplée des deux fluides. Le maillage adaptatif a été implémenté avec les librairies parallèles PARAMESH. Celles-ci gèrent la création d'un maillage construit par blocs qui peuvent être bissectés de façon récursive afin d'obtenir la résolution désirée. Les blocs ont tous le même nombre de cellules, ce qui permet une répartition efficace de la charge de travail en parallle. Un ensemble d'outils nécessaires à une correcte résolution des équations ont été développés, à partir d'un robuste solveur elliptique BiCG-stab préconditionné par méthode multigrille. Le code a été vérifié sur des cas tests académiques capable de maintenir la précision sur les grilles plus fines. Les performances du code en termes de réduction de charge de travail et d'efficacité de la parallelisation ont été également illustrées. Le code a été enfin testé sur la désintégration assistée d'une nappe liquide bidimensionnelle cisaillée par des courants gazeux à haute vitesse. Le code s'est avéré capable de retrouver certains phénomènes physiques comme l'oscillation longitudinale de la nappe, ainsi que de permettre une simulation multi échelles, ce qui permet de plus se rapprocher des conditions des injecteurs réels.