Beaucoup de procédés utilisés dans l'industrie sidérurgique (coulée de lingots,coulée continue, …) peuvent générer des défauts : macro-ségrégation, mauvaises propriétés de la microstructure, défauts surfaciques. Ces problèmes peuvent être résolus par un contrôle de la température et de l’écoulement d'acier liquide. Le brassage électromagnétique (EMS) est une technique largement utilisée pour contrôler l’écoulement d'acier liquide par l’imposition d'un champ électromagnétique. Cette technique est complexe car elle couple plusieurs types de problèmes physiques:écoulement multiphasique, solidification,transfert de chaleur et induction électromagnétique à basse fréquence.En outre, l’approche expérimentale est difficile de par la dimension,l'environnement et le coût des procédés considérés. Pour ces raisons, des simulations numériques efficaces sont nécessaires pour comprendre les applications EMS et améliorer les procédés évoqués. L'objectif de cette thèse est de développer une méthodologie numérique robuste,efficace et précise pour la simulation multi-physique de l'EMS, en particulier pour le brassage dans le moule dans le cadre de la coulée continue d'acier. Cette méthodologie a été mise en oeuvre dans le code commercial THERCAST® pour être utilisé dans le cadre d’applications industrielles / Many of the processes used in thesteelmaking industry (e.g. ingot casting,continuous casting, …) can lead todefects: macro-segregation, poormicrostructure properties, surfacedefects. These issues can be solved bycontrolling the temperature and the flowof molten steel. Electromagnetic stirring(EMS) is a widely used technique to steerthe flow of liquid steel by thesuperimposition of an electro-magneticfield. This application is complex becauseit couples several physical problems:multi-phase flow, solidification, heattransfer and low frequency electromagneticinduction. In addition,experimental work is difficult because ofthe size, environment and cost of theconsidered processes. For thesereasons, efficient and effective numericalsimulations are needed to understandEMS applications and improve theaforementioned processes.The objective of this thesis is to developa robust, efficient and accurate numericalprocedure for the multi-physicssimulation of EMS, especially for in-moldstirring in the framework of continuouscasting of steel. This procedure has beenimplemented in the commercial codeTHERCAST® in order to be used forindustrial applications.

No description available.

The incompressible Navier-Stokes equations have proven formidable for nearly a century. The present difficulties are mathematical and computational in nature; the computational requirements, in particular, are exponentially exacerbated in the presence of high Reynolds number. The issues are further compounded with the introduction of markers or an immiscible fluid intended to be tracked in an ambient high Reynolds number flow; despite the overwhelming pragmatism of problems in this regime, and increasing computational efficacy, even modest problems remain outside the realm of direct approaches. Herein three approaches are presented which embody direct application to problems of this nature. An LES model based on an entropy-viscosity serves to abet the computational resolution requirements imposed by high Reynolds numbers and a one-stage compressive flux, also utilizing an entropy-viscosity, aids in accurate, efficient, conservative transport, free of low order dispersive error, of an immiscible fluid or tracer. Finally, an integral commutator and the theory of anti-dispersive spaces is introduced as a novel theoretical tool for consistency error analysis; in addition the material engenders the construction of error-correction techniques for mass lumping schemes.

[pt] Esta tese propõe um método meshfree adaptativo de advecção para problemas de fluxo bifásico de fluidos incompressíveis e imiscíveis em meios porosos heterogêneos tridimensionais. Este método se baseia principalmente na combinação do método Semi-Lagrangeano adaptativo com interpolação local meshfree usando splines poliharmônicas como funções de base radial. O método proposto é uma melhoria e uma extensão do método adaptativo meshfree AMMoC proposto por Iske e Kaser (2005) para modelagem 2D de reservatórios de petróleo. Inicialmente este trabalho propõe um modelo em duas dimensões, contribuindo com uma melhoria significativa no cálculo do Laplaciano, utilizando os métodos meshfree de Hermite e Kansa. Depois, o método é ampliado para três dimensões (3D) e para um meio poroso heterogêneo. O método proposto é testado com o problema de five spot e os resultados são comparados com os obtidos por sistemas bem conhecidos na indústria de petróleo. / [en] This thesis proposes an adaptive meshfree advection method for two-phase flow problems of incompressible and immiscible fluids through three-dimensional heterogeneous porous media. This method is based mainly on a combination of adaptive semi-Lagrangian method with local meshfree interpolation using polyharmonic splines as radial basis functions. The proposed method is an improvement and extension of the adaptive meshfree advection scheme AMMoC proposed by Iske and Kaser (2005) for 2D oil reservoir modeling. Initially this work proposes a model in two dimensions, contributing to a significant improvement in the calculation of the Laplacian, using the meshfree methods of Hermite and Kansa. Then, the method is extended to three dimensions (3D) and a heterogeneous porous medium. The proposed method is tested with the five spot problem and the results are compared with those obtained by well-known systems in the oil industry.

La reconstruction en tomographie discrète de la microstructure de l’os joue un role très important pour le diagnostic de l’ostéoporse, une maladie des os très fréquente. Le diagnostic clinique est basé sur l’absortiométrie duale de rayons X. Avec la tomographie de rayons X, une résolution spatiale élevée avec des images reconstruites in vivo requiert une dose d’irradiation élevée et un temps de balayage long, ce qui est dangereux pour le patient. Une des méthodes pour résoudre ce problème est de limiter le nombre de projections. Cependant, avec cette méthode le problème de reconstruction devient mal posé. Deux types de régularisation par Variation Totale minimisées avec la méthode Alternate Direction of Minimization Method (ADMM) et deux schémas basés sur les méthodes de régularisation Level-set sont appliquées à deux images d’os expérimentales acquises avec un synchrotron (pixel size: 15 μm). Des images de tailles variées et avec différents niveaux de bruit Gaussien additifs ajoutés aux projections sont utlisées pour étudier l’efficacité des méthodes de régularisation. Des minima locaux sont obtenus avec ces méthodes déterministes. Une approche globale d’optimisation est nécessaire pour améliorer les résultats. Des perturbations stochastiques peuvent être un moyen très utile pour échapper aux minima locaux. Dans une première approche, une équation différentielle stochastique basée sur la régularisation level-set est étudiée. Cette méthode améliore les résultats de reconstruction mais ne modifie que les frontières entre les régions 0 et 1. Ensuite une équation aux dérivées partielles stochastique est obtenue avec la régularisation TV pour améliorer la méthode stochastique level-set. A la fin de notre travail, nous avons étendu la méthode de régularisation à des images 3D avec des données réelles. Cette algorithme a été implémenté avec RTK. Nous avons aussi étendu l’approche level-set utilisée pour la tomographie binaire au cas multi-level. / Discrete tomography reconstruction of bone microstructure is important in diagnosis of osteoporosis. One way to reduce the radiation dose and scanning time in CT imaging is to limit the number of projections. This method makes the reconstruction problem highly ill-posed. A common solution is to reconstruct only a finite number of intensity levels. In this work, we investigate only binary tomography reconstruction problem. First, we consider variational regularization methods. Two types of Total Variation (TV) regularization approaches minimized with the Alternate Direction of Minimization Method (ADMM) and two schemes based on Level-set (LS) regularization methods are applied to two experimental bone cross-section images acquired with synchrotron micro-CT. The numerical experiments have shown that good reconstruction results were obtained with TV regularization methods and that level-set regularization outperforms the TV regularization for large bone image with complex structures. Yet, for both methods, some reconstruction errors are still located on the boundaries and some regions are lost when the projection number is low. Local minima were obtained with these deterministic methods. Stochastic perturbations is a useful way to escape the local minima. As a first approach, a stochastic differential equation based on level-set regularization was studied. This method improves the reconstruction results but only modifies the boundaries between the 0 and 1 regions. Then partial stochastic differential equation obtained with the TV regularization semi-norm were studied to improve the stochastic level-set method. The random change of the boundary are performed in a new way with the gradient or wavelet decomposition of the reconstructed image. Random topological changes are included to find the lost regions in the reconstructed images. At the end of our work, we extended the TV regularization method to 3D images with real data on RTK (Reconstruction Toolkit). And we also extended the level-set to the multi-level cases.

This doctoral thesis proposes a new algorithm for the reconstruction of impedance images in monitored objects. The algorithm eliminates the spatial resolution problems present in existing reconstruction methods, and, with respect to the monitored objects, it exploits both the partial knowledge of configuration and the material composition. The discussed novel method is designed to recognize certain significant fields of interest, such as material defects or blood clots and tumors in biological images. The actual reconstruction process comprises two phases; while the former stage is focused on industry-related images, with the aim to detect defects in conductive materials, the latter one concentrates on biomedical applications. The thesis also presents a description of the numerical model used to test the algorithm. The testing procedure was centred on the resulting impedivity value, influence of the regularization parameter, initial value of the numerical model impedivity, and effect exerted by noise on the voltage electrodes upon the overall reconstruction results. Another issue analyzed herein is the possibility of reconstructing impedance images from components of the magnetic flux density measured outside the investigated object. The given magnetic field is generated by a current passing through the object. The created algorithm for the reconstruction of impedance images is modeled on the proposed algorithm for EIT-based reconstruction of impedance images from voltage. The algoritm was tested for stability, influence of the regularization parameter, and initial conductivity. From the general perspective, the thesis describes the methodology for both magnetic field measurement via NMR and processing of the obtained data.

After a review of the notions of Hausdorff and Fourier dimensions from fractal geometry and Fourier analysis and the properties of local times of Brownian motion, we study the Fourier structure of Brownian level sets. We show that if δa(X) is the Dirac measure of one-dimensional Brownian motion X at the level a, that is the measure defined by the Brownian local time La at level a, and μ is its restriction to the random interval [0, L−1 a (1)], then the Fourier transform of μ is such that, with positive probability, for all 0 ≤ β < 1/2, the function u → |u|β|μ(u)|2, (u ∈ R), is bounded. This growth rate is the best possible. Consequently, each Brownian level set, reduced to a compact interval, is with positive probability, a Salem set of dimension 1/2. We also show that the zero set of X reduced to the interval [0, L−1 0 (1)] is, almost surely, a Salem set. Finally, we show that the restriction μ of δ0(X) to the deterministic interval [0, 1] is such that its Fourier transform satisfies E (|ˆμ(u)|2) ≤ C|u|−1/2, u 6= 0 and C > 0. Key words: Hausdorff dimension, Fourier dimension, Salem sets, Brownian motion, local times, level sets, Fourier transform, inverse local times. / Decision Sciences / PhD. (Operations Research)

Elastography, which is based on applying pressure and estimating the resulting deformation, involves the forward problem to obtain the strain distributions and inverse problem to construct the elastic distributions consistent with the obtained strains on observation points. This thesis focuses on the former problem whose solution is used as an input to the latter problem. The aim is to provide the inverse problem community with accurate strain estimates of a coronary artery vessel wall. In doing so, a new ultrasonic image-based elastography approach is developed. Because the accuracy and quality of the estimated strain fields depend on the resolution level of the ultrasound image and to date best resolution levels obtained in the literature are not enough to clearly see all boundaries of the artery, one of the main goals is to acquire high-resolution coronary vessel wall ultrasound images at different pressures. For this purpose, first an experimental setup is designed to collect radio frequency (RF) signals, and then image formation algorithm is developed to obtain ultrasound images from the collected signals. To segment the noisy ultrasound images formed, a geodesic active contour-based segmentation algorithm with a novel stopping function that includes local phase of the image is developed. Then, region-based information is added to make the segmentation more robust to noise. Finally, elliptical deformable template is applied so that a priori information regarding the shape of the arteries could be taken into account, resulting in more stable and accurate results. The use of this template also implicitly provides boundary point correspondences from which high-resolution, size-independent, non-rigid and local strain fields of the coronary vessel wall are obtained.

Les phénomènes d'érosion sont une cause majeure dans la rupture des ouvrages hydrauliques tels que les digues ou les barrages. Le territoire compte 10000 km de digues, d'où un enjeu important en terme de risque d'inondation. Une approche plus générale est essentielle à une compréhension fine des phénomènes. La modélisation des mécanismes de l'érosion doit permettre l'élaboration de solutions pour sécuriser les ouvrages et les rendre moins sensibles à ces phénomènes. Ce travail est consacré à la modélisation d'une interface entre un milieu de sol et un milieu d'écoulement. Le mouvement de cette interface est régi par l'érosion du sol sous l'effet de l'écoulement. L'écoulement est calculé par les équations de Stokes munies d'un terme de pénalisation permettant la représentation unifiée des comportements des deux milieux. Dans une première partie, nous aborderons la problématique de la représentation de l'érosion. La seconde partie traite plus précisément de la représentation unifiée des comportements des milieux, notamment grâce à l'utilisation des méthodes de pénalisation. La troisième partie constitue une revue des principales méthodes permettant le suivit d'une interface. Une description précise de la méthode Level Set est donnée. Dans la quatrième partie, nous abordons la résolution numérique des équations d'écoulement. La cinquième partie rassemble l'ensemble des tests de la méthode. Enfin, la dernière partie est consacrée à l'application du modèle sur trois cas : le " Hole Erosion Test ", le phénomène de suffusion et l'érosion d'un lit de rivière autour d'une pile de pont.

The main contribution of this thesis is the implementation of manufacturing constraints in shape and topology optimization. Fabrication limitations related to the casting process are formulated as mathematical constraints and introduced in the optimization algorithm. In addition, based on the same theoretical and modelization tools, we propose a novel formulation for multi-phase optimization problems, which can be extended to the optimization of structures with functionally-graded properties. A key ingredient for the mathematical formulation of most problems throughout our work is the notion of the signed distance function to a domain. This work is divided into three parts. The rst part is bibliographical and contains the necessary background material for the understanding of the thesis' main core. It includes the rst two chapters. Chapter 1 provides a synopsis of shape and topology optimization methods and emphasizes the combination of shape sensitivity analysis and the level-set method for tracking a shape's boundary. In Chapter 2 we give a short description of the casting process, from which all our manufacturing constraints derive. We explain how industrial designers account for these limitations and propose a strategy to incorporate them in shape and topology optimization algorithms. The second part is about the mathematical formulation of manufacturing constraints. It starts with Chapter 3, where the control of thickness is discussed. Based on the signed distance function, we formulate three constraints to ensure a maximum and minimm feature size, as well as a minimal distance between structural members. Then, in Chapter 4, we propose ways to handle molding direction constraints and combine them with thickness constraints. Finally, a thermal constraint coming from the solidi cation of cast parts is treated in Chapter 5 using several thermal models. Multi-phase optimization is discussed in the third part. The general problem of shape and topology optimization using multiple phases is presented in detail in Chapter 6. A "smoothed-interface" approach, based again on the signed distance function, is proposed to avoid numerical di culties related to classical "sharp-interface" problems and a shape derivative is calculated. An extension of this novel formulation to general types of material properties' gradation is shown in the Appendix A.