Ce travail de thèse de doctorat propose et évalue différentes solutions pour caractériser, avec des outils optiques et électromagnétiques non intrusifs, les nanoparticules et agrégats observés dans différents systèmes physiques : suspensions colloïdales, aérosols et plasma poussiéreux. Deux types de modèles sont utilisés pour décrire la morphologie d'agrégats fractals (p. ex. : suies issues de la combustion) et agrégats compacts (qualifiés de « Buckyballs » et observés dans des aérosols produits par séchage de nano suspensions). Nous utilisons différentes théories et modèles électromagnétiques (T-Matrice et approximations du type dipôles discrets) pour calculer les diagrammes de diffusion (ou facteur de structure optique) de ces agrégats, de même que leurs spectres d'extinction sur une large gamme spectrale. Ceci, dans le but d'inverser différentes données expérimentales. Différents outils numériques originaux ont également été mis au point pour parvenir à une analyse morphologique quantitative de clichés obtenus par microscopie électronique. La validation expérimentale des outils théoriques et numériques développés au cours de ce travail est focalisée sur la spectrométrie d'extinction appliquée à des nano agrégats de silice, tungstène et silicium. / This Ph.D. work provides and evaluates various solutions to characterize, with optical/electromagnetic methods nanoparticles and aggregates of nanoparticles found in suspensions, aerosols and dusty plasmas. Two main models are introduced to describe the morphology of particle aggregates with fractal-like (for particles in plasmas and combustion systems) and Buckyballs-like (aerosols, suspensions) shapes. In addition, the author proposes various solutions and methods (T-Matrix, Rayleigh type approximations) to calculate the scattering diagrams (optical structure factors) of fractal aggregates as well as algorithms to inverse extinction spectra. As a reference case for the performed analysis, several tools to describe the morphology of fractal aggregates from electron microscopy images have been also developed. The experimental validation carried out with the Light Extinction Spectrometry (LES) technique (for nano silica beads, tungsten, dusty plasma and silicon aggregates) clearly proves the validity of the algorithms developed as well as the potential of the LES technique.

La profilométrie optique est une technique de métrologie de surface rapide et non destructive. Dans ce mémoire, nous avons abordé cette problématique par des méthodes inverses de diffraction électromagnétique et dans une configuration de type Microscopie Tomographique Optique par Diffraction (ODTM). La surface est sondée par un éclairement sous plusieurs angles d'incidences ; la mesure en amplitude et en phase du champ lointain diffracté constitue les données du problème. Des profils de surfaces ont été reconstruits en considérant différents modèles de diffraction, parmi lesquelles une méthode approchée fondée sur les approximations de diffusion simple et de paraxialité. La résolution latérale de cette méthode et des techniques classiques de profilométrie est limitée par le critère d'Abbe-Rayleigh, défini sur la base de l'ouverture numérique pour l'éclairement et la détection du champ. Afin de dépasser cette limite de résolution, nous avons développé une méthode itérative de Newton-Kantorovitch régularisée. L'opérateur de diffraction y est rigoureusement modélisé par une méthode des moments, résolution numérique des équations du formalisme intégral de frontière, et l'expression de la dérivée de Fréchet de cet opérateur est obtenue par la méthode des états adjoints, à partir du théorème de réciprocité. Pour les surfaces unidimensionnelles métalliques, notre technique permet d'inverser à partir de données synthétiques des surfaces très rugueuses avec une résolution au delà du critère d'Abbe-Rayleigh. / Optical profilometry is a nondestructive and fast noncontact surface metrology technique. In this thesis, we have tackled this issue with inverse scattering electromagnetic methods and in an Optical Digital Tomographic Microscopy (ODTM) configuration. The surface is probed with illuminations under several incidence angles; the measure of far scattered field amplitude and phase constitutes the problem data. Surface profiles have been reconstructed using different scattering models among which an approximate theory based on single scattering and paraxiality. The lateral resolution of this technique and classical profilometric approaches is limited by the so-called Abbe-Rayleigh's criterion defined out of the numerical aperture for illumination and field detection. In order to overpass this resolution limit, we have developed a regularized iterative Newton-Kantorovitch's method. The scattering operator is rigorously modelized with the method of moments, that is a numerical solution of boundary integral equations, and its Fréchet derivative adjoint states expression is deduced from the reciprocity theorem. For one-dimensional metallic surfaces, our method succeeds in inverting from synthetic data very rough surfaces with the resolutions beyond the Abbe-Rayleigh's criterion. The performance of this technique and inversion conditions clearly differ from one polarization to the other : in the TM case, interactions at longer distance than in the TE case improve yet the resolution. This work includes also an experimental validation of our inverse model on grooves in indium phosphure substrate at 633 nm.

Ce travail vise à approfondir la connaissance du comportement mécanique de polymères structurés sur de multiples échelles. Il repose sur l'acquisition simultanée d'informations couplées relatives à la mécanique, à la thermique et à la microstructure. Nous avons à cet effet étudié le comportement du Polyéthylène Haute Densité (PEHD) sous sollicitations de traction à déformation contrôlée à l'aide de trois techniques optiques utilisées in situ simultanément : (i) Le système Vidéotraction® pour l'accès aux contraintes et aux déformations vraies. (ii) La thermographie infrarouge pour l'accès aux sources thermiques. (iii) Une technique de rétrodiffusion de lumière (ISLT) pour l'accès à une forme d'endommagement. Le corps du travail porte sur la reconstruction des sources thermiques 2D à partir de cartographies des champs de température mesurés au cours d'essais. Le problème inverse de reconstruction des champs de sources ther-miques est résolu par deux méthodes différentes : minimisation sous contrainte avec formulation adjointe et décomposition modale. Des techniques de régularisation du problème sont décrites en détail pour chaque méthode. D'autre part, la technique ISLT nous permet de suivre l'endommagement : l'apparition progressive de "mi-crocavités" isotropes puis un développement d'une forte anisotropie pendant la phase de durcissement hyperélastique. Nous montrons que ces informations sur la microstructure viennent corroborer les phénomènes observés sur la source thermique. Nous replaçons alors l'ensemble des informations obtenues dans le contexte d'une modélisation thermodynamique des lois de comportement. Une étude de sensibilité aux paramètres du modèle est développée, ce qui conduit à une réduction du modèle adaptée à leur bonne estimation / This research aims at looking further into current knowledge of the mechanical behavior of semi-crystalline polymers at multiple scales. It relies on the simultaneous acquisition of coupled information relating to me-chanical, thermal and microstructural behaviors. For this purpose, we have studied the response of High Density Polyethylene (HDPE) under tensile tests with controlled strain rates, using three optical techniques : (i) Vidéotraction®, a videoextensemeter to access the true stress and strain. (ii) Infrared thermography to access the thermal sources evolution. (iii) Backscattering of Incoherent Light (ISLT technique) to access information on damage mechanisms. The core of the research concerns the reconstruction of the 2D thermal sources starting from temperature fields maps measured during the test. This ill-posed problem has been solved with two different methods : A constrained optimization based on the adjoint formulation and a minimization based on a reduced spectral model. The regularization tools used in the problem are detailed for each method. The ISLT technique allows to monitor damaging processes due to crazing : the creation of isotropic micro-cavities followed by the development of a strong anisotropy during the hardening phase. We show that these information at microstructural level confirm the evolutions of the heat source thermal. All these experimental information are finally considered within a thermodynamical framework used for deri-ving constitutive laws. A sensitivity analysis applied to the parameters of our law led to a reduced model that is more appropriate for their proper estimation

L’imagerie optique par colorant potentiométrique est une méthode d’enregistrement de l’activité corticale prometteuse, mais dont le potentiel réel est limité par la présence d’artefacts et d’interférences dans les acquisitions. À partir de modèles existant dans la littérature, nous proposons un modèle génératif du signal basé sur un mélange additif de composantes, chacune contrainte dans une union d’espaces linéaires déterminés par son origine biophysique. Motivés par le problème de séparation de composantes qui en découle, qui est un problème inverse linéaire sous-déterminé, nous développons : (1) des régularisations convexes structurées spatialement, favorisant en particulier des solutions parcimonieuses ; (2) un nouvel algorithme proximal de premier ordre pour minimiser efficacement la fonctionnelle qui en résulte ; (3) des méthodes statistiques de sélection de paramètre basées sur l’estimateur non biaisé du risque de Stein. Nous étudions ces outils dans un cadre général, et discutons leur utilité pour de nombreux domaines des mathématiques appliqués, en particulier pour les problèmes inverses ou de régression en grande dimension. Nous développons par la suite un logiciel de séparation de composantes en présence de bruit, dans un environnement intégré adapté à l’imagerie optique par colorant potentiométrique. Finalement, nous évaluons ce logiciel sur différentes données, synthétiques et réelles, montrant des résultats encourageants quant à la possibilité d’observer des dynamiques corticales complexes. / Voltage-sensitive dye optical imaging is a promising recording modality for the cortical activity, but its practical potential is limited by many artefacts and interferences in the acquisitions. Inspired by existing models in the literature, we propose a generative model of the signal, based on an additive mixtures of components, each one being constrained within an union of linear spaces, determined by its biophysical origin. Motivated by the resulting component separation problem, which is an underdetermined linear inverse problem, we develop: (1) convex, spatially structured regularizations, enforcing in particular sparsity on the solutions; (2) a new rst-order proximal algorithm for minimizing e›ciently the resulting functional; (3) statistical methods for automatic parameters selection, based on Stein’s unbiased risk estimate.We study thosemethods in a general framework, and discuss their potential applications in variouselds of applied mathematics, in particular for large scale inverse problems or regressions. We develop subsequently a soŸware for noisy component separation, in an integrated environment adapted to voltage-sensitive dye optical imaging. Finally, we evaluate this soŸware on dišerent data set, including synthetic and real data, showing encouraging perspectives for the observation of complex cortical dynamics.

Les travaux présentés dans cette thèse se focalisent sur l'étude de l’effet de dégradations des composants de puissance, plus particulièrement au niveau de l’environnement proche des puces (métallisations, connexions, brasures puces/DCB), sur le comportement électrique et thermique des puces ainsi que de leur assemblage. Pour ce faire nous avons cherché à étudier la répartition des courants et des températures à la surface de la puce à l’aide d’un modèle électrothermique 2D distribué. Nous avons aussi évalué l’effet de la dégradation des brasures dans le volume de l’assemblage, à l’aide cette fois d’un modèle thermique relié à la constitution de l’assemblage. La première partie de cette thèse consiste à mettre en place un modèle électrothermique distribué de puce MOSFET, qui tient compte à la fois du caractère distribué de la dissipation de la puissance et du couplage électrothermique en régime transitoire. Ce modèle électrothermique s’appuie sur un modèle électrique aux variables d’états et un modèle thermique par éléments finis couplé au modèle électrique. Les modèles électriques et thermiques ont été développés respectivement sous Matlab et sous CAST3M, et le couplage des deux modèles a été fait sous Simulink. Dans une deuxième partie, pour la validation des résultats des températures et pour l’analyse de l’effet du vieillissement et des dégradations (sur la distribution et la dynamique de température de la surface supérieure de la puce), une méthodologie de mesure rapide de température et un banc expérimental pour thermographie infrarouge ont été mis en place. Les difficultés rencontrées lors des mesures thermiques IR sous variation rapide de la température nous ont poussé à envisager d’autres méthodes d’analyse thermique. Enfin, nous avons cherché à évaluer la réponse impulsionnelle du composant testé en estimant, par des simulations thermiques, la fonction de transfert dans le domaine fréquentiel à l’aide du logiciel COMSOL Multiphysics. Nous avons également étudié la pertinence de modèles RC équivalents (réseau RC de Cauer). Ces modèles ont ensuite été utilisés pour rendre compte de différents modes de dégradation notamment cette fois au niveau des couches de brasures entre puce et DCB et entre DCB et semelle. Mots clef : Modules de puissance à semi-conducteur, Vieillissement, Métallisation, Modélisation électrothermique, Court-circuit, Distribution de courant et de température, Problème inverse, Caméra IR, Réseaux de Cauer. / The work presented in this thesis focus on the study of the effect of degradation of power components, especially at the near environment of chips (metallization, connections, solder chips / DCB), on the electrical and thermal behavior of the chips and their assembly. As a consequence, we studied the distribution of currents and temperatures on the chip surface with a 2D electrothermal distributed model. We also evaluated the effect of solder degradation in the volume of the assembly. Firstly, we developed an electrothermal distributed model of the MOSFET chip, which takes into account both the distributed power dissipation and the electrothermal coupling transient. This electrothermal model is based on an electrical model of state variables and thermal finite element model coupled to the electric model. Electrical and thermal models were developed respectively in Matlab and CAST3M whereas the two models coupling was done in Simulink . In the second part, to validate the results of temperatures and to analyze the effect of ageing and degradation on the distribution and dynamics of temperature of the upper surface of the chip, methodology rapid temperature measurement and an experimental bench for infrared thermography were established. The difficulties encountered in IR thermal measurements with rapid temperature change led us to consider other thermal analysis methods. Eventually, we assessed the impulse response of the tested component by estimating with thermal simulations, the transfer function in the frequency domain using the COMSOL Multiphysics software. Moreover we evaluated the relevance of RC equivalent models (RC Cauer network). These models were then used to account for different modes of degradation this time especially on the solder layer between the chip and DCB and between the DCB and sole. Keywords: Power Modules semiconductor, Ageing, Metallization, electrothermal modeling, Short Circuit, Power and temperature distribution, inverse problem, IR Camera, Cauer networks.

Ce travail concerne l'imagerie micro-onde en vue d'application à l'imagerie biomédicale. Cette technique d'imagerie a pour objectif de retrouver la distribution des propriétés diélectriques internes (permittivité diélectrique et conductivité) d'un objet inconnu illuminé par une onde interrogatrice connue à partir des mesures du champ électrique dit diffracté résultant de leur interaction. Un tel problème constitue un problème dit inverse par opposition au problème direct associé qui consiste à calculer le champ diffracté, l'onde interrogatrice et l'objet étant alors connus.La résolution du problème inverse nécessite la construction préalable du modèle direct associé. Celui-ci est ici basé sur une représentation intégrale de domaine des champs électriques donnant naissance à deux équations intégrales couplées dont les contreparties discrètes sont obtenues à l'aide de la méthode des moments. En ce qui concerne le problème inverse, hormis le fait que les équations physiques qui interviennent dans sa modélisation directe le rendent non-linéaire, il est également mathématiquement mal posé au sens de Hadamard, ce qui signifie que les conditions d'existence, d'unicité et de stabilité de la solution ne sont pas simultanément garanties. La résolution d'un tel problème nécessite sa régularisation préalable qui consiste généralement en l'introduction d'information a priori sur la solution recherchée. Cette résolution est effectuée, ici, dans un cadre probabiliste bayésien où l'on introduit une connaissance a priori adaptée à l'objet sous test et qui consiste à considérer ce dernier comme étant composé d'un nombre fini de matériaux homogènes distribués dans des régions compactes. Cet information est introduite par le biais d'un modèle de « Gauss-Markov-Potts ». De plus, le calcul bayésien nous donne la distribution a posteriori de toutes les inconnues connaissant l'a priori et l'objet. On s'attache ensuite à déterminer les estimateurs a posteriori via des méthodes d'approximation variationnelles et à reconstruire ainsi l'image de l'objet recherché. Les principales contributions de ce travail sont d'ordre méthodologique et algorithmique. Elles sont illustrées par une application de l'imagerie micro-onde à la détection du cancer du sein. Cette dernière constitue en soi un point très important et original de la thèse. En effet, la détection du cancer su sein en imagerie micro-onde est une alternative très intéressante à la mammographie par rayons X, mais n'en est encore qu'à un stade exploratoire. / This work concerns the problem of microwave tomography for application to biomedical imaging. The aim is to retreive both permittivity and conductivity of an unknown object from measurements of the scattered field that results from its interaction with a known interrogating wave. Such a problem is said to be inverse opposed to the associated forward problem that consists in calculating the scattered field while the interrogating wave and the object are known. The resolution of the inverse problem requires the prior construction of the associated forward model. This latter is based on an integral representation of the electric field resulting in two coupled integral equations whose discrete counterparts are obtained by means of the method of moments.Regarding the inverse problem, in addition to the fact that the physical equations involved in the forward modeling make it nonlinear, it is also mathematically ill-posed in the sense of Hadamard, which means that the conditions of existence, uniqueness and stability of the solution are not simultaneously guaranteed. Hence, solving this problem requires its prior regularization which usually involves the introduction of a priori information on the sought solution. This resolution is done here in a Bayesian probabilistic framework where we introduced a priori knowledge appropriate to the sought object by considering it to be composed of a finite number of homogeneous materials distributed in compact and homogeneous regions. This information is introduced through a "Gauss-Markov-Potts" model. In addition, the Bayesian computation gives the posterior distribution of all the unknowns, knowing the a priori and the object. We proceed then to identify the posterior estimators via variational approximation methods and thereby to reconstruct the image of the desired object.The main contributions of this work are methodological and algorithmic. They are illustrated by an application of microwave imaging to breast cancer detection. The latter is in itself a very important and original aspect of the thesis. Indeed, the detection of breast cancer using microwave imaging is a very interesting alternative to X-ray mammography, but it is still at an exploratory stage.

La machine à collecteur mécanique reste, dans de nombreuses applications, une alternative crédible par rapport aux autres types de moteurs électriques. Dans le domaine automobile, ce sont les coûts de fabrication et surtout la fiabilité dû à l'absence d'électronique de puissance qui rendent cette solution si intéressante. Cela est également son point faible : l'absence de dispositif de régulation fragilise l'électronique de bord en cas d'appel de courant trop important au démarrage. L'usure du système balai-collecteur est également un point clé pour un démarreur d'autant plus que, pour assurer un fonctionnement « Stop-Start » en plus du fonctionnement traditionnel, le nouveau cahier des charges impose jusqu'à 450 000 cycles de démarrage (contre 30 000 auparavant). Ceci justifie la recherche actuelle sur ces dispositifs. En effet, une partie importante de cette thèse traite de la modélisation et de la mesure de la commutation dans le démarreur pour mieux évaluer les phénomènes qui y interagissent. L'arc électrique en fin de commutation joue un rôle important dans l'usure du système balais-collecteur. Une meilleure connaissance des facteurs influents permet, par la suite, de concevoir des appareils plus robustes. Dans une deuxième partie de la thèse, on s'intéresse au dimensionnement optimal d'un démarreur. Après avoir passé en revue les différentes méthodes d'optimisation pour traiter un tel problème et avoir sélectionné les plus adaptées, on traite un cas pratique de dimensionnement optimal. Cette problématique est fondamentale pour Valeo au vue des quantités d'appareils fabriqués et des futurs marchés offerts par le développement du « micro-hybride » qui se généralise / Nowadays, brushed DC motors are still used in many applications. In automotive industry, manufacturing cost and reliability are the two most important key indicators. That's why DC motors are preferred to others solutions requiring power electronics supply. But, there are also some drawbacks ; when the inrush current occurs (when the power-contact is switched on) the DC voltage of the battery decrease and can damage electronic components of the car. The wear of the brush-segments system is also a key issue for the life time of the starter. The new specifications for a « Stop-Start » starter use lead to a significant increase of the number of re-start cycles to be done by the device. We consider that a car achieve more than 450000 stop-and-start cycles during its lifetime whereas this number is only 30000 for the standard vehicles. These specifications justify the research work lead by Valeo on such device. An important part of the thesis concerns the modeling and the measurement of the commutation phenomena. The electric arc occurring at the brush trailing edges (when the physical contact ends) causes the wear of the brush-commutator system. A better knowledge of influent factors can lead to more robust device development. In a second part, we study the optimal design of a car starter regarding specifications. After a review of different optimization methods for dealing with such a problem, we select the most relevant ones and solve different practical problem around the starter. This issue is very important regarding the quantities of device manufactured and the market of micro-hybrid which is growing up very fast

L'Imagerie par Résonance Magnétique (IRM) présente deux caractéristiques principales. La première, sa capacité à manipuler le contraste, constitue son principal avantage par rapport aux autres modalités d'imagerie. Cela permet d'obtenir des informations complémentaires pour une meilleure détectabilité et une meilleure précision dans le diagnostic. Cela est particulièrement appréciable pour les pathologies du myocarde. La seconde caractéristique de l'IRM est également l'un de ces principaux inconvénients : le processus d'acquisition est relativement lent. De ce fait, les mouvements du patient constituent un obstacle important puisqu'ils perturbent ce processus d'acquisition, ce qui se traduit par des artéfacts dans l'image reconstruite. L'imagerie cardiaque et abdominale sont donc particulièrement sensibles à cette problématique du mouvement. L'objectif de cette thèse est donc de proposer une méthode de correction de mouvement intégrable dans un contexte multi-contraste. Nous avons étudié dans un premier temps la question de la correction de mouvement seule. Pour cela, nous nous sommes plus particulièrement intéressés à la méthode GRICS déjà développée au laboratoire IADI. Cette méthode permet la reconstruction conjointe d'une image sans artéfact et d'un modèle de mouvement non rigide permettant de corriger les déplacements qui surviennent pendant l'acquisition. Le premier apport majeur de cette thèse a consisté à améliorer la méthode GRICS, notamment pour l'adapter à l'imagerie volumique 3D. Il s'agit d'une nouvelle méthode de régularisation adaptative particulièrement adaptée au problème inverse posé dans GRICS. Le second apport majeur de cette thèse a consisté à gérer la correction de mouvement GRICS de manière conjointe sur des acquisitions présentant des contrastes différents. Il s'agit de concevoir l'examen IRM comme un tout et d'exploiter au mieux les informations partagées entre les différents contrastes. Toutes ces méthodes ont été appliquées et validées par des simulations, des tests sur fantôme, sur volontaires sains et sur des patients dans la cadre d'études cliniques. L'application cardiaque a été particulièrement visée. Les méthodes développées ont permis d'améliorer le processus d'acquisition et de reconstruction dans le contexte clinique réel / Magnetic Resonance Imaging (MRI) has two main features. The first one, its ability to manipulate contrast, is a major advantage compared to the other imaging modalities. It allows to access complementary information for a better detectability and a diagnostic more accurate. This is especially useful for myocardium pathologies. The second feature of MRI is also one of its main drawbacks: the acquisition process is slow. Therefore, patient motion is a significant obstacle because it disturbs the acquisition process, which leads to artifacts in the reconstructed image. Cardiac and thoracic imaging are particularly sensitive to this motion issue. The aim of this thesis is to develop a new motion correction method that can be integrated in a multi-contrast workflow. In a first phase, we studied apart the motion correction problem. To do so, we focused more particularly on the GRICS method which was already developed in the IADI laboratory. This method allows the joint reconstruction of an image free from artifact and a non-rigid motion model that describes the displacements occurring during the acquisition. The first major contribution of this thesis is an improvement of the GRICS method consisting mainly in adapting it to the 3D imaging. This was achieved with a new adaptive regularization method that perfectly suits the inverse problem posed in GRICS. The second major contribution of this thesis consists in the simultaneous management of the motion correction on multiple acquisitions with different contrasts. To do so, the MRI examination is considered as a whole. Thus we make the most of information shared between the different contrasts. All these methods have been applied and validated by simulations, tests on phantom, on healthy volunteers and on patients as part of clinical studies. We aimed more particularly at cardiac MR. Finally the developed methods improve the acquisition and reconstruction workflow in the framework of a real clinical routine

Les simulations numériques des écoulements sanguins cardiovasculaires peuvent combler d’importantes lacunes dans les capacités actuelles de traitement clinique. En effet, elles offrent des moyens non invasifs pour quantifier l’hémodynamique dans le cœur et les principaux vaisseaux sanguins chez les patients atteints de maladies cardiovasculaires. Ainsi, elles permettent de recouvrer les caractéristiques des écoulements sanguins qui ne peuvent pas être obtenues directement à partir de l’imagerie médicale. Dans ce sens, des simulations personnalisées utilisant des informations propres aux patients aideraient à une prévision individualisée des risques. Nous pourrions en effet, disposer des informations clés sur la progression éventuelle d’une maladie ou détecter de possibles anomalies physiologiques. Les modèles numériques peuvent fournir également des moyens pour concevoir et tester de nouveaux dispositifs médicaux et peuvent être utilisés comme outils prédictifs pour la planification de traitement chirurgical personnalisé. Ils aideront ainsi à la prise de décision clinique. Cependant, une difficulté dans cette approche est que, pour être fiables, les simulations prédictives spécifiques aux patients nécessitent une assimilation efficace de leurs données médicales. Ceci nécessite la solution d’un problème hémodynamique inverse, où les paramètres du modèle sont incertains et sont estimés à l’aide des techniques d’assimilation de données.Dans cette thèse, le problème inverse pour l’estimation des paramètres est résolu par une méthode d’assimilation de données basée sur un filtre de Kalman d’ensemble (EnKF). Connaissant les incertitudes sur les mesures, un tel filtre permet la quantification des incertitudes liées aux paramètres estimés. Un algorithme d’estimation de paramètres, basé sur un filtre de Kalman d’ensemble, est proposé dans cette thèse pour des calculs hémodynamiques spécifiques à un patient, dans un réseau artériel schématique et à partir de mesures cliniques incertaines. La méthodologie est validée à travers plusieurs scenarii in silico utilisant des données synthétiques. La performance de l’algorithme d’estimation de paramètres est également évaluée sur des données expérimentales pour plusieurs réseaux artériels et dans un cas provenant d’un banc d’essai in vitro et des données cliniques réelles d’un volontaire (cas spécifique du patient). Le but principal de cette thèse est l’analyse hémodynamique spécifique du patient dans le polygone de Willis, appelé aussi cercle artériel du cerveau. Les propriétés hémodynamiques communes, comme celles de la paroi artérielle (module de Young, épaisseur de la paroi et coefficient viscoélastique), et les paramètres des conditions aux limites (coefficients de réflexion et paramètres du modèle de Windkessel) sont estimés. Il est également démontré qu’un modèle appelé compartiment d’ordre réduit (ou modèle dimension zéro) permet une estimation simple et fiable des caractéristiques du flux sanguin dans le polygone de Willis. De plus, il est ressorti que les simulations avec les paramètres estimés capturent les formes attendues pour les ondes de pression et de débit aux emplacements prescrits par le clinicien. / Cardiovascular blood flow simulations can fill several critical gaps in current clinical capabilities. They offer non-invasive ways to quantify hemodynamics in the heart and major blood vessels for patients with cardiovascular diseases, that cannot be directly obtained from medical imaging. Patient-specific simulations (incorporating data unique to the individual) enable individualised risk prediction, provide key insights into disease progression and/or abnormal physiologic detection. They also provide means to systematically design and test new medical devices, and are used as predictive tools to surgical and personalize treatment planning and, thus aid in clinical decision-making. Patient-specific predictive simulations require effective assimilation of medical data for reliable simulated predictions. This is usually achieved by the solution of an inverse hemodynamic problem, where uncertain model parameters are estimated using the techniques for merging data and numerical models known as data assimilation methods.In this thesis, the inverse problem is solved through a data assimilation method using an ensemble Kalman filter (EnKF) for parameter estimation. By using an ensemble Kalman filter, the solution also comes with a quantification of the uncertainties for the estimated parameters. An ensemble Kalman filter-based parameter estimation algorithm is proposed for patient-specific hemodynamic computations in a schematic arterial network from uncertain clinical measurements. Several in silico scenarii (using synthetic data) are considered to investigate the efficiency of the parameter estimation algorithm using EnKF. The usefulness of the parameter estimation algorithm is also assessed using experimental data from an in vitro test rig and actual real clinical data from a volunteer (patient-specific case). The proposed algorithm is evaluated on arterial networks which include single arteries, cases of bifurcation, a simple human arterial network and a complex arterial network including the circle of Willis.The ultimate aim is to perform patient-specific hemodynamic analysis in the network of the circle of Willis. Common hemodynamic properties (parameters), like arterial wall properties (Young’s modulus, wall thickness, and viscoelastic coefficient) and terminal boundary parameters (reflection coefficient and Windkessel model parameters) are estimated as the solution to an inverse problem using time series pressure values and blood flow rate as measurements. It is also demonstrated that a proper reduced order zero-dimensional compartment model can lead to a simple and reliable estimation of blood flow features in the circle of Willis. The simulations with the estimated parameters capture target pressure or flow rate waveforms at given specific locations.

Cette thèse a pour objectif de fournir aux utilisateurs de SMPS (Scanning Mobility Particle Sizer) une méthodologie pour calculer les incertitudes associées à l’estimation de la granulométrie en nombre des aérosols. Le résultat de mesure est le comptage des particules de l’aérosol en fonction du temps. Estimer la granulométrie en nombre de l’aérosol à partir des mesures CNC revient à considérer un problème inverse sous incertitudes.Une revue des modèles existants est présentée dans le premier chapitre. Le modèle physique retenu fait consensus dans le domaine d’application.Dans le deuxième chapitre, un critère pour l’estimation de la granulométrie en nombre qui couple les techniques de régularisation et de la décomposition sur une base d’ondelettes est décrit.La nouveauté des travaux présentés réside dans l’estimation de granulométries présentant à la fois des variations lentes et des variations rapides. L’approche multi-échelle que nous proposons pour la définition du nouveau critère de régularisation permet d’ajuster les poids de la régularisation sur chaque échelle du signal. La méthode est alors comparée avec la régularisation classique. Les résultats montrent que les estimations proposées par la nouvelle méthode sont meilleures (au sens du MSE) que les estimations classiques.Le dernier chapitre de cette thèse traite de la propagation de l’incertitude à travers le modèle d’inversiondes données. C’est une première dans le domaine d’application car aucune incertitude n’est associée actuellement au résultat de mesure. Contrairement à l’approche classique qui utilise un modèle fixe pour l’inversion en faisant porter l’incertitude sur les entrées, nous proposons d’utiliser un modèle d’inversion aléatoire (tirage Monte-Carlo) afin d’intégrer les erreurs de modèle. Une estimation moyenne de la granulométrie en nombre de l’aérosol et une incertitude associée sous forme d’une région de confiance à 95 % est finalement présentée sur quelques mesures réelles. / This thesis aims to provide SMPS (Scanning Mobility Particle Sizer) users with a methodology to compute the uncertainties associated with the estimation of aerosol size distributions. Recovering aerosol size distribution yields to consider an inverse problem under uncertainty.The first chapter of his thesis presents a review of physical models and it shows that competitive theories exist to model the physic.A new criterion that couples regularization techniques and decomposition on a wavelet basis is described in chapter 2 to perform the estimation of the size distribution.Main improvement of this work is brought when size distributions to be estimated show both broad and sharp profiles. The multi-scale approach helps to adjust the weights of the regularization on each scale of the signal. The method is then tested against common regularization and shows better estimates (in terms of the mean square error).Last chapter of this thesis deals with the propagation of the uncertainty through the data inversion process.Results from SMPS measurements are not given with any uncertainty at this time so providing end-users with an uncertainty is already a real improvement. Common approach is to consider a fixed physical model and to model the inputs (particle count) as random. We choose to consider both the physical model as well as the inputs as random to account for the model error. The result is expressed as a mean estimate of the size distribution with a 95% coverage region. The all methodology is finally tested on real measurements.