Dans une série de travaux Zied Ammari et Francis Nier ont développé des méthodes pour étudier la dynamique de champ moyen bosonique pour des états quantiques généraux pouvant présenter des corrélations. Ils ont obtenu des formules pour décrire la dynamique des corrélations, ou plus généralement des matrices densité réduites d'ordre arbitraire. Cette thématique a été largement développée ces dernières années. Norbert Mauser en a été un des contributeurs, ainsi que sur la notion de mesure de Wigner qui est la clé de l'analyse développée par Z. Ammari et F. Nier. En général, il est admis que l'asymptotique de champ moyen est une bonne approximation du problème à N particules quand N dépasse la dizaine. Cela concerne l'asymptotique de la matrice densité réduite à une particule qui ne décrit pas la dynamique des corrélations. Un objectif est de tester la validité de la dynamique de champ moyen pour les matrices densité réduites à 2-particules. Pour des tests numériques, les modèles discrets qui n'ont pas été vraiment traités en détail dans les travaux précédents de Z. Ammari et F. Nier semblent bien adaptés. La thèse comprendra donc plusieurs étapes: adapter les résultats précédents de Z. Ammari et F. Nier à des modèles discrets , développer des méthodes numériques pour des systèmes simples mais pertinents, permettant de valider l'approximation de champ moyen et les formules pour la dynamique des corrélations. Au niveau numérique, on utilise des schémas numériques symplectiques, développés spécifiquement ces dernières années pour la discrétisation des équations hamiltoniennes. Une dernière étape concerne la combinaison des deux asymptotiques, champ moyen et approximation des modèles continus par les modèles discrets. / In a serie of works Z. Ammari and F. Nier developed methods to study the dynamics of bosonic mean field for general quantum states which can present correlations. They obtained formulas to describe the dynamics of the correlations, or more generally reduced density matrices with an arbitrary order. This topic was widely developed these last years. N.J. Mauser was one of contributors, as well as on the notion of Wigner measure which is the key of the analysis developed by Z. Ammari and F. Nier. Generally, the mean field asymptotic is admitted is a good approximation of the N-body problem when N exceed about ten. It concerns the asymptotics of the reduced density matrices for one particle which does not describe the dynamics of the correlations. An objective is to test the validity of the mean field dynamics for reduced density matrices for 2 particles. For numerical tests, the discrete models which were not really handled in detail in the previous works of Z. Ammari and F. Nier seem adapted well. The thesis will thus include several steps: adapt the previous results from Z. Ammari and F. Nier to discrete models , develop numerical methods, for simple but relevant systems, allowing to validate the approximation of mean field and the formulas for the dynamics of the correlations. About numerics, symplectic numerical scheme are used, developed specifically these last years for the discretization of the hamiltonian equations. A last possible step concerns the combination of both asymptotics, that is mean field and approximation of the continuous models by the discrete models.

Ce travail a pour but de proposer une technique de mesure permettant d’identifier les bas nombres d’onde de pressions pariétales dues à des écoulements turbulents à faible nombre de Mach. En effet, les caractéristiques de ces excitations vibroacoustiques dans ces nombres d’onde sont mal connues alors qu’ils peuvent être des sources principales de vibration et de rayonnement acoustique dans des problèmes d’ingénierie qui concernent essentiellement le secteur industriel des transports. Cette méconnaissance vient du fait que la composante acoustique de l’excitation qui est située dans les bas nombre d’onde,a une amplitude très petite devant celle de la composante aérodynamique et que l’énergie d’origine acoustique se trouve noyée dans le bruit de mesure lorsqu’on utilise des capteurs de pression.L’objectif de cette thèse est d’étudier comment la méthode inverse vibratoire de Résolution Inverse (RI) et ses variantes RIFF (Filtrée Fenêtrée) ou RIC (Corrigée)permettent d’identifier ces bas nombres d’onde. L’intérêt d’utiliser de telles méthodes inverses vibratoires est que la structure est utilisée comme capteur. Ainsi, les composantes responsables des vibrations, se situant dans les bas nombres d’onde, peuvent être mieux identifiées ou extraites puisqu’elles sont naturellement filtrées par la dynamique de la structure.Dans un premier temps, la méthode RIFF est testée dans le cas de simulations d’une couche limite turbulente. Les résultats obtenus permettent de comprendre comment cette méthode identifie l’excitation dans les bas nombres d’onde. La méthode peut ainsi paraître particulièrement intéressante, puisqu’elle permet d’extraire une composante très faible, mais fortement responsable des vibrations et du bruit rayonné par la plaque et pratiquement impossible à observer en utilisant des microphones affleurants.Dans un deuxième temps, la mise en oeuvre de ces techniques est étudiée en appliquant la méthode RIC. Cette approche permet de diminuer drastiquement le nombre de capteurs et de rendre la méthode inverse utilisable avec les moyens actuels. La méthode est ainsi testée sur la simulation d’une couche limite turbulente et sur une expérimentation en soufflerie où l’écoulement turbulent est généré par une marche montante. La proposition d’un indicateur permet alors d’identifier si les composantes acoustiques et aérodynamiques sont fortement séparées et si les résultats obtenus par la méthode RIC correspondent uniquement à la composante acoustique ou non. / The aim of this work is to propose a measurement technique for the identification of the low wavenumbers of wall pressures due to turbulent flows with low Mach number. Indeed, the characteristics of these vibroacoustic excitations in these wavenumbers are not well known whereas they can be the principal sources of vibration and acoustic radiation in a lot of engineering problems in the transport industries. This lack of knowledge is due to the low amplitude of the acoustic component of the excitation that is localized in the low wavenumbers and the acoustic energy is below the level of measurement noise when pressure sensors are used.This thesis aims to study how the inverse methods of vibration called the Force AnalysisTechnique (FAT) and its Corrected version CFAT identify these low wavenumbers.The advantage of using these inverse methods of vibration is that the structure is used as a sensor. Thus, the components responsible for the vibrations, and localized in the low wavenumbers, can be identified and extracted as they are naturally filtered by the dynamics of the structure.At first, the FAT method is tested in a simulation of a turbulent boundary layer.The results are used to demonstrate how this method identifies the excitation in the low wavenumbers. The method can appear particularly interesting, since it allows to extract a very small component that is highly responsible for vibration and noise radiated by the plate and that is impossible to measure using flush-mounted sensors like microphones.In a second study, the implementation of these techniques is studied by applying the CFAT method. This approach allows to drastically reduce the number of sensors and make the inverse method usable with existing techniques of measurement. The method is tested on the simulation of a turbulent boundary layer and on an experimentation in a wind tunnel where a turbulent flow is generated by a forward-facing step. A new indicator is introduced and identifies if the acoustic and aerodynamic components arestrongly separated and if the results obtained by the CFAT method correspond only to the acoustic component or not.

Trois essais cliniques ARISTOTLE, RE-LY et ROCKET-AF, sont réalisés sur les quatre nouveaux anticoagulants oraux qui sont l’Apixaban, le Dabigatran 110mg, Dabigatran 150mg et le Rivaroxaban, afin d’améliorer la prise en charge de la fibrillation auriculaire (FA). Ils ne permettent pas de déterminer la meilleure option thérapeutique. L’objectif de cette thèse est de développer une approche d’aide multicritère à la décision pour la problématique du choix thérapeutique appliquée au cas de la FA afin d’évaluer ces quatre options thérapeutiques médicamenteuses. La méthodologie d’aide multicritère développée, PROMETHEE-GAIA appliquée à la FA permet de comparer ces quatre nouveaux anticoagulants (l’Apixaban, le Dabigatran 110mg, Dabigatran 150mg et le Rivaroxaban) en se basant sur les cinq critères essentiels élaborés selon un processus de concertation cadrée, qui sont l’efficacité, la sécurité, la fonction rénale, l’observance et le prix. Les résultats de l’évaluation de ces quatre nouveaux anticoagulants est un classement de ces options thérapeutiques par leur ordre de performance dans la prise en charge des patients de la FA selon ces cinq critères en considération.

Une meilleure compréhension des relations entre la structure complexe de l'os cortical et ses propriétés mécaniques est nécessaire à l'évaluation de la qualité osseuse. Les méthodes conventionnelles ex vivo de mesure de l'élasticité à l'échelle du millimètre ont des limitations liées à l'anisotropie du tissu, à son inhomogénéité et à la petite taille des échantillons. Au contraire, la spectroscopie par résonance ultrasonore (RUS) est bien adaptée à la mesure de petits échantillons anisotropes. Cette méthode estime l'élasticité à partir des fréquences de résonance de l'échantillon, et l'amortissement à partir de la largeur des pics de résonance. Son application à l'os était considérée difficile, du fait de l'amortissement important des modes de vibration, qui induit un recouvrement des pics de résonance et complique la mesure des fréquences. Pour surmonter cette difficulté, des adaptions de la méthode - dans la mesure, le traitement du signal et l'estimation des propriétés du matériau - ont été proposées. Elles ont été validées sur de l'os cortical et sur des échantillons de polymère et de matériau composite imitant l'os. La précision de la méthode a été démontrée, ainsi que sa capacité à mesurer tous les termes du tenseur d'élasticité à partir d'un seul échantillon. De plus, une nouvelle formulation Bayésienne de l'inversion apporte une solution automatique à un problème qui nécessitait une stratégie fastidieuse d'essai-erreur ou de complexes modifications du dispositif expérimental. Finalement, l'application à une grande collection d'échantillons de tibias humains démontre que la méthode RUS pourrait être utilisée en routine pour la mesure de la viscoélasticité de l'os. / Deep understanding of the structure-function relationships of cortical bone in the context of bone quality assessment is still missing. Currently available methods to measure millimeter-scale elasticity ex vivo have limitations arising from theanisotropy of the tissue, its heterogeneity, and the small size of the specimens.Resonant ultrasound spectroscopy is particularly suitable for the measurement of small anisotropic specimens. This method estimates elasticity from the free resonant frequencies of a specimen, and damping from the width of the resonant peaks. Its application to cortical bone was considered challenging because of the high damping of the vibrations modes, which causes overlapping of the resonant peaks and prevents a direct measurement of the resonant frequencies. To overcome the difficulty,adaptations of all the steps of RUS – measurement, signal processing and inverse estimation of the material properties – have been introduced. Validation of each step of the procedure has been achieved by application to several test samples, including a cortical bone specimen and bone-mimicking composite and polymer specimens.RUS was shown to be precise and accurate, with the advantage of providing the complete stiffness tensor from the measurement of a single specimen. Additionally, an original Bayesian formulation of the inversion provides an automated solution toa problem that was previously solved by tedious trial-and-error procedures or complex additions to the basic experimental setup. Finally, the application to a large collection of human tibiae specimens demonstrates that RUS can be considered a routine method to characterize the viscoelasticity of bone.

Dans cette thèse, nous étudions des problèmes de minimisation d'énergies discrètes et nous cherchons à comprendre pourquoi une structure périodique peut être un minimiseur pour une énergie d'interaction, c'est ce que l'on appelle un problème de cristallisation. Après avoir montré qu'un réseau de R^d soumis à un certain potentiel paramétré peut être vu comme un minimum local, nous démontrons des résultats d'optimalité du réseau triangulaire parmi les réseaux de Bravais du plan pour certaines énergies par point, avec ou sans densité fixée. Finalement, nous démontrons, à partir des travaux de Sandier et Serfaty sur les gaz de Coulomb bidimensionnels, la conjecture de Rakhmanov-Saff-Zhou, c'est-à-dire l'existence d'un terme d'ordre n dans le développement asymptotique de l'énergie logarithmique optimale pour n points sur la sphère unité de R^3. De plus, nous montrons l'équivalence entre la conjecture de Brauchart-Hardin-Saff portant sur la valeur de ce terme d'ordre n et celle de Sandier-Serfaty sur l'optimalité du réseau triangulaire pour une énergie coulombienne renormalisée / In this thesis, we study minimization problems for discrete energies and we search to understand why a periodic structure can be a minimizer for an interaction energy, that is called a crystallization problem. After showing that a given Bravais lattice of R^d submitted to some parametrized potential can be viewed as a local minimum, we prove that the triangular lattice is optimal, among Bravais lattices of R^2, for some energies per point, with or without a fixed density. Finally, we prove, from Sandier and Serfaty works about 2D Coulomb gases, Rakhmanov-Saff-Zhou conjecture, that is to say the existence of a term of order n in the asymptotic expansion of the optimal logarithmic energy for n points on the 2-sphere. Furthermore, we show the equivalence between Brauchart-Hardin-Saff conjecture about the value of this term of order n and Sandier-Serfaty conjecture about the optimality of triangular lattice for a coulombian renormalized energy

L'utilisation des microondes pour le diagnostic est en plein essor dans le domaine médical. Une des toutes dernières applications concerne la détection d'AVC (Accident vasculaire Cérébral) par imagerie microonde. La Société EMTensor GmbH basée à Vienne en Autriche étudie actuellement un tel système en collaboration avec le LEAT, le LJAD de l’Université Côte d’Azur et le LJLL de Sorbonne Université, pour le diagnostic et le contrôle de l'efficacité de traitement. Le but de ce travail était de modéliser le système de mesure de l'imagerie du cerveau, développé par la société EMTensor GmbH. Il s'agit d'un système d'émission/réception composé de 160 antennes disposées en 5 couronnes de 32 antennes réparties sur une cuve métallique cylindrique de section circulaire semi-ouverte. Un des enjeux majeurs de ce travail consiste en la modélisation et la simulation électromagnétique (EM) du système complet incluant un modèle réaliste de cerveau. La difficulté réside à la fois dans la taille du problème EM à simuler en raison du rapport entre la taille considérable du système et la taille très faible de certaines inhomogénéités à l’intérieur du cerveau, et dans la grande hétérogénéité des permittivités diélectriques présentes à l’intérieur du cerveau. Nous avons décidé d’utiliser un code open source, FreeFem++ pour cette modélisation car il permet de déployer du calcul hautement parallèle et la décomposition de domaines, qui sont bien adaptés à la complexité du problème EM. Dans un premier temps, nous avons comparé les résultats de simulation du système de mesure à vide (sans le cerveau) aux mesures et aux résultats obtenus par le logiciel de simulation EM HFSS basé sur la FEM comme FreeFem++. Nous avons ensuite simulé un modèle de tête tridimensionnel virtuel, à partir de coupe d’image du cerveau (scanner et IRM), en partenariat avec EMTensor en recherchant la position et le type d'AVC (ischémique et hémorragique). L'influence du bruit de mesure, la valeur du gel d'adaptation utilisé, le couplage entre les capteurs et le couplage entre la tête et les capteurs sont également étudiés. Afin de valider ces modèles, deux cas simples ont été étudiés. Un grand tube et un petit tube en plastique sont remplis de liquide d'adaptation symbolisant les caractéristiques diélectriques d'un cerveau afin de retrouver la forme du tube utilisé. Nous avons montré qu’il est possible de développer des algorithmes de reconstruction pour montrer permettant de retrouver la forme des objets par imagerie qualitative. Enfin, avec les partenaires et l'entreprise d'EMTensor nous avons appliqué une méthode quantitative à la détection d’un AVC ischémique par la tomographie microonde. Le problème direct repose sur l’utilisation de FreeFem++, en utilisant des éléments d'ordre supérieur et des préconditionneurs parallèles pour la méthode de décomposition de domaine. Nous avons résolu le problème inverse par un algorithme de minimisation, afin de reconstruire des images tomographiques du cerveau dans des temps compatibles avec les impératifs médicaux définis par les cliniciens. / The use of microwaves for diagnosis is booming in the medical field. One of the latest applications is the detection of strokes by microwave imaging. The company EMTensor GmbH based in Vienna, Austria is currently studying such a system in collaboration with LEAT, the LJAD of the Côte d’Azur University and the LJLL of Sarbonne University, for the diagnosis and control of the treatement efficiency. The purpose of this work is to model the brain imaging measurement system developed by EMTensor GmbH. It is a transmission/ reception system consisting of 160 antennas arranged in 5 rings of 32 antennas distributed on a cylinder metal tank of semi-open circular section. One of the major issues of this work is the modeling and electromagnetic simulation (EM) of the complete system including a realistic brain model. The difficulty lies both in the size of the EM problem to be simulated beacause of the relationship between the considerable size of the system and the the very small size of certain inhomogeneities within the brain, and the great heterogeneity of the dielectric permittivities present inside the brain. We decided to use an open source software, FreeFem++ for this modelling because it is well adapted to high performance computing through domain decomposition methods, which is mandatory for the complexity of the EM problem. First, we compared the simulation results of the vacuum matching measurement system (without the brain) to the measurements and the results obtained by the FEM-based EM HFSS simulation software to those obtained by FreeFem++. We then simulated a virtual threedimensional head model, from brain imaging system cuts (CT scan and MRI), in partnership with EMTensor, looking for the position and type of stroke (ischemic and hemorragic). The influence of the measurement noise, the value of the adaptation gel used, the coupling between the sensors and the coupling between the head and the sensors are also studied. In order to validate these models, two simple cases have been studied. A large tube and a small plastic tube are fielld with adaptation liquid with the dielectric characteristic of a brain to find the shape of the tubes used by qualitative imaging. Finally, with the MEDIMAX project partners and the EMTensor company we applied a quantitative method to the detection of ischemic stroke by the microwave tomography. The direct problem has been solved with the help of FreeFem++, using hight order elements and parallel preconditioners for the domain decomposition method. We solved the inverse problem by a minimization algorithm, in order to reconstruct tomographic images of the brain in times compatible with medical imperatives defined by clinicians.”

Cette thèse se compose de deux parties principales. Dans la première partie, le Chapitre 3 est consacré à l'étude du comportement limite d'un système contrôlé singulièrement perturbé avec deux variables d'état qui sont faiblement couplées. Afin de prouver notre résultat d'approximation, nous utilisons la méthode de moyennisation et un résultat récent sur le contrôle nonexpansif. La principale nouveauté de notre approche est de permettre la dynamique limite de dépendre de l'état initial du système rapide. Notons que dans la littérature, le comportement limite d'un tel système a été généralement traité dans des conditions qui garantissent que la limite est indépendante de l'état initial du système rapide. Dans le Chapitre 4, nous généralisons les résultats du Chapitre 3 supposant une condition de nonexpansivité plus générale. De plus, nous considérons un exemple ou la nouvelle condition de nonexpansivité est satisfaite, mais pas la condition de nonexpansivité du Chapitre 3. Dans la deuxième partie de la thèse, le Chapitre 5 porte sur les représentations stables des Hamiltoniens convexes associant à un Hamiltonien donné des fonctions correspondant au problème de Bolza en controle optimal. Dans le Chapitre 6 nous étudions également la stabilité des solutions des équations d'Hamilton-Jacobi-Bellman sous contraintes d'état en exploitant la stabilité des fonctions valeur d'une famille de problèmes de contrôle optimal de Bolza sous contraintes d'état. Nous montrons que sous des hypothèses appropriées, la fonction valeur est la solution unique d'équation d'Hamilton-Jacobi-Bellman et que les solutions sont stables par rapport à l'Hamiltonien et les contraintes d'état. / This thesis consists of two main parts. In the first part, Chapter 3 is devoted to the investigation of the limit behavior of a singularly perturbed control system with two state variables which are weakly coupled. In order to prove our approximation result we use the so called averaging method and a recent result on nonexpansive control. The main novelty of our averaging approach lies in the fact that the limit dynamic may depend on the initial condition of the fast system. In the literature, the investigation of the limit behavior of such systems has been usually addressed under conditions that ensure that the limit dynamic is independent from the initial condition of the fast system. In Chapter 4, we generalise the results of Chapter 3 by considering a more general nonexpansivity condition. Moreover, we consider an example where the new nonexpansity condition is satisfied but the nonexpansivity condition of Chapter 3 does not hold true. The second part deals with Hamilton-Jacobi equations under state constraints. Chapter 5 focuses on the stable representation of convex Hamiltonians by functions describing a Bolza optimal control problem. In Chapter 6 we investigate stability of solutions of Hamilton-Jacobi-Bellman equations under state constraints by studying stability of value functions of a suitable family of Bolza optimal control problems under state constraints. We show that under suitable assumptions, the value function is a unique viscosity solution to Hamilton-Jacobi-Bellman equation and that solutions are stable with respect to Hamiltonians and state constraints.

Cette thèse porte sur un problème couplé des lanceurs, à savoir une manœuvre de l'attitude couplée avec la trajectoire minimisant le temps de manœuvre. La difficulté de ce problème vient essentiellement du phénomène de chattering et du couplage des dynamiques n'ayant pas la même échelle de temps. Avec une analyse géométrique des extrémales venant de l'application du principe du maximum de pontryagin, nous donnons des conditions suffisantes sous lesquelles le phénomène de chattering se produit, pour des systèmes affines bi-entrée. Nons appliquons ensuite ce résultat à notre problème, et montrons que le phénomène de chattering arrive pour les trajectoires optimales, pour certaines données terminales. A l'aide de cette analyse théorique préliminaire, nous mettons en œuvre une méthode de résolution indirecte efficace, combinée à une méthode de continuation prédicteur-correcteur. En cas de chattering, deux stratégies sous-optimales sont proposées: soit une méthode directe dont le contrôle est approché par un contrôle constant par morceaux, soit en stoppant la continuation avant l'échec dû au chattering. Avec le tir multiple et plusieurs paramètres de continuations supplémentaires, cette méthode de résolution est appliquée à chercher une manœuvre de pull-up avec des contraintes sur l'état en minimisant le temps-énergie pour des lanceurs aéroportés. Les résultats numériques permettent de mettre en évidence l'efficacité et la robustesse de notre méthode de résolution. / In this thesis, we investigate the minimum time control problem for the control and guidance of a launch vehicle, whose motion is described by its attitude kinematics and dynamics but also by its trajectory dynamics. The difficulty of this problem is essentially due to the chattering phenomenon and to the coupling of dynamics of different time scales. With a refined geometric study of the extremals coming from the application of the pontryagin maximum principle, we establish a general result for bi-input control-affine systems, providing sufficient conditions under which the chattering phenomenon occurs. We show how this result can be applied to our problem. Based on this preliminary theoretical analysis, we implement an efficient indirect numerical method, combined with numerical predictor-corrector continuation, in order to compute numerically the optimal solutions of the problem. In case of chattering, two sub-optimal strategies are designed: one is a direct method in which the control is approximated by a piecewise constant control, and the other consists of stopping the continuation procedure before its failure due to chattering. With several additional numerical continuation steps, we apply finally the developed indirect approach to the minimum time-energy pull-up maneuver problem, in which state constraints are also considered, for airborne launchers. Numerical simulations illustrate the efficiency and robustness of our method.

On identifie simultanément les coefficients d'emmagasinement et de transmissivité hydraulique dans un écoulement souterrain gouvernent par une équation parabolique linéaire. Ces deux paramètres sont supposés être des fonctions constantes par morceaux en espace. Les inconnues du problème sont non seulement les valeurs de ces coefficients mais aussi la géométrie des zones dans lesquelles ces coefficients sont constants. Le problème est formule comme la minimisation d'une fonction de moindres carres calculant la différence entre les mesures et les quantités correspondantes évaluées avec la valeur courante des paramètres. L'objectif principal de ce travail est la construction d'une technique de paramétrisation adaptative guidée par des indicateurs de raffinement. L'utilisation d'indicateurs de raffinement, nous permet de construisons la paramétrisation de façon itérative, on allant d'une paramétrisation à une seule zone à une paramétrisation avec m zones où m est une valeur optimale à identifier. Nous distinguons les cas ou les deux paramètres ont la même paramétrisation et le cas où les deux paramètres ont des paramétrisations différentes. Pour améliorer la résolution du problème inverse d'estimation de paramètres, nous incorporons des estimateurs d'erreurs a posteriori. / We identify simultaneously storage and hydraulic transmissivity coefficients in groundwater flow governed by a linear parabolic equation. Both parameters are assumed to be functions piecewise constant in space. The unknowns are the coefficient values as well as the geometry of the zones where these coefficients are constant. This problem is formulated as minimizing a least-square function calculating the difference between measurements and the corresponding quantities computed with the current parameters values. The main point of this work is to construct an adaptative parameterization technique guided by refinement indicators. Using refinement indicators, we build the parameterization iteratively, going from a one zone parametrization to a parametrization with $m$ zones where $m$ is an optimal value to identify. We distinguish the cases where the two parameters have the same parameterization and different parameterizations.\\ To improve the resolution of the inverse problem, we incorporate a posteriori error estimations.

La tomographie en 2D et 3D sont largement utilisée dans l’imagerie médicale ainsi que dans le Contrôle Non Destructif (CND) pour l’industrie. Dans toutes les deux applications, il est nécessaire de réduire le nombre de projections. Dans certains cas, la reconstruction doit être faite avec un nombre d’angle de projections limité. Les données mesurées sont toujours avec des erreurs (erreurs de mesure et de modélisation). Nous sommes donc presque toujours dans la situation de problèmes inversés mal posés. Le rôle des méthodes probabilistes et de la modélisation a priori devient crucial. Pour la modélisation a priori, en particulier dans les applications NDT, l’objet à l’examen est composé de plusieurs matériaux homogènes, avec plusieurs blocs continus séparés par des discontinuités et des contours. Ce type d’objet est dit continu par morceaux. L’objet de cette thèse est sur la reconstruction des objets continu ou constante par morceaux, ou plus généralement homogène par morceaux. En résumé, deux méthodes principales sont proposées dans le contexte de l’inférence bayésienne. La première méthode consiste à reconstruire l’objet en imposant que sa transformée de Haar soit parcimonieuse. Un modèle bayésien hiérarchique est proposé. Dans cette méthode, les variables et les paramètres sont estimés et les hyper-paramètres sont initialisés selon la définition des modèles antérieurs. La deuxième méthode reconstruit les objets en estimant simultanément les contours. L’objet continu par morceaux est modélisé par un modèle markovien non-homogène, qui dépend du gradient de l’objet, et le gradient dépend aussi de l’estimation de l’objet. Cette méthode est également semi-supervisé, avec les paramètres estimés automatiquement. Ces méthodes sont adaptées aux reconstructions de grande taille de données 3D, dans lesquelles le processeur GPU est utilisé pour accélérer les calculs. Les méthodes sont validées avec des données simulées et des données réelles, et sont comparées avec plusieurs méthodes classiques. / 2D and 3D X-ray Computed Tomography (CT) is widely used in medical imaging as well as in Non Destructive Testing (NDT) for industrial applications. In both domains, there is a need to reduce the number of projections. In some cases we may also be limited in angles. The measured data are always with errors (measurement and modelling errors). We are consequently almost always in the situation of ill-posed inverse problems. The role of the probabilistic methods and the prior modelling become crucial. For prior modelling, in particular in NDT applications, the object under examination is composed with several homogeneous materials, with several continuous blocs separated by some discontinuities and contours. This type of object is called the piecewise-continuous object. The focus of this thesis on the reconstruction of the picewise continuous or constant, or more generally piecewise homogeneous objects. In summary two main methods are proposed in the context of the Bayesian inference. The first method consists in reconstructing the object while enforcing the sparsity of the discrete Haar transformation coefficients of the object. A hierarchical Bayesian model is proposed. In this method, the unknown variables and parameters are estimated and the hyper-parameters are initialized according to the definition of prior models. The second method reconstruct objects while the contours are estimated simultaneously. The piecewise continuous object is modeled by a non-homogeneous Markovian model, which depends on the gradient of the object, while the gradient also depends on the estimation of the object. In this methods, the semi-supervised system model is also achieved, with the parameters estimated automatically. Both methods are adapted to the 3D big data size reconstructions, in which the GPU processor is used to accelerate the computation. The methods are validated with both simulated and real data, and are compared with several conventional state-of-the-art methods.