La radiothérapie interne sélective est une thérapie émergente très peu invasive du carcinome hépatocellulaire, quatrième cause de décès par cancer dans le monde. Des millions de microsphères chargées en Yttrium 90 sont injectées dans l'artère hépatique par un cathéter. Actuellement, leur distribution lors d'une injection est estimée par l'injection préliminaire d'un radiomarqueur, ce qui peut se révéler trop approximatif. Un traitement personnalisé permettrait une concentration des radiations à la tumeur tout en épargnant le tissu sain environnant. Dans ce travail je me suis intéressée au développement d'un modèle numérique, pour une simulation spécifique à chaque patient des trajectoires des microsphères, dans le but d'optimiser le traitement. Le protocole clinique d'imagerie a été exploité et optimisé pour l'extraction de données spécifiques patients telles que le foie, les tumeurs, l'artère hépatique et le flux sanguin. Les tissus et l'artère hépatique (jusqu'à un diamètre de 0.05 mm) sains et malins ont été simulés. Cela nous permet de simuler la distribution des microsphères dans le tissue hépatique, validée grâce à la scintigraphie post-traitement. Il est supposé ici que les microsphères se distribuent de façon proportionnelle au flux sanguin, lequel est modélisé par la loi de Poiseuille. Des simulations plus approfondies en mécanique de fluides numérique du flux sanguin ont ensuite été réalisées dans l'artère hépatique du patient. Pour cela nous avons utilisé et comparé les méthodes des Volumes Finis (Ansys Fluent) et de Lattice Boltzmann (programme développé dans le laboratoire). Le transport des microsphères a été simulé dans l'artère hépatique du patient avec la méthode des volumes finis, et dans une géométrie simplifiée avec la méthode de Lattice Boltzmann. Une séquence IRM de contraste de phase a aussi été optimisée pour l'extraction de la vitesse du sang dans l'artère hépatique, dans le but de valider le modèle numérique. / Selective internal radiation therapy using Yttrium-90 loaded glass microspheres injected in the hepatic artery is an emerging, minimally invasive therapy of hepatocellular carcinoma, which is the fourth cause of mortality in the world. Currently, microspheres distribution can be only approximately predicted by the injection of a radiotracer, whose behaviour may be different. A personalised intervention can lead to high concentration dose in the tumour, while sparing the surrounding parenchyma. This work is concerned with the development of a patient-specific numerical model for the simulation of microspheres trajectories and treatment optimisation. Clinical imaging protocol is utilised and optimised in order to extract patient’s specific data such as liver, tumours, hepatic artery and blood flow. Normal and malignant hepatic arterial vasculature and tissues are simulated down to a vessels diameter of 0.05 mm. A preliminary simulation of microspheres distribution in liver tissue is proposed and validated against post-treatment scintigraphy. Microspheres are here supposed to distribute proportionally to blood flow, which is computed based on Poiseuille’s law. More precise computational fluid dynamics (CFD) simulations of blood flow in the patient’s segmented arteries are performed. The Finite Volume Method (Ansys Fluent) and the Lattice Boltzmann Method (in-house developed software) are used to this purpose and their efficacy is compared. Microspheres transport is simulated in the patient’s hepatic artery using the FVM, and in a representative geometry using the LBM method. A phase contrast MRI sequence has been optimised in order to extract blood velocity from the hepatic artery and validate CFD simulations.

Dans ces travaux de recherche nous avons cherché à développer et optimiser un procédé d'oxydation avancée dédié au traitement de l'eau. Pour cela, nous avons développé un procédé combinant la génération de radicaux hydroxyles par décharge plasma à pression atmosphérique et un catalyseur fixe. L'originalité de ce travail a été d'utiliser un système d'écoulement permettant l'obtention d'un film liquide (épaisseur<1mm) afin de permettre une interaction entre le plasma, le liquide et la surface catalytique. Les résultats ont montré que l'utilisation d'une décharge de type DBD était préférable au Glidarc, thermique, énergétique et produisant trop de NOx. Les valeurs d'efficacité respectives étaient de 0,67g.kWh-1 pour la DBD et 0,23g.kWh-1 pour le Glidarc. Les résultats ont également montré que le matériau de la contre-électrode/plaque d'écoulement pouvait influencer l'efficacité du procédé. Ainsi dans le cas du laiton, la présence de réactions de corrosion produisant des NO2- réduit l'efficacité de moitié comparée à une plaque d'acier. Le développement d'un modèle CFD 1D à partir des données expérimentales du procédé a été effectué. Ce modèle a permis d'estimer que l'efficacité énergétique de production des HO est de 4,4.10-9 molHO.J-1 pour la décharge DBD. Il a également permis de montrer que le facteur majoritairement limitant était la diffusion moléculaire des espèces dans le liquide. En effet les réactions chimiques ne dure pas plus de 1µs après l'impact du streamer et n'ont lieu que dans les premiers 5% d'épaisseur de liquide. Enfin, un dépôt catalytique de DLC dopé à l'argent et réalisé par rf-PECVD a permis d'améliorer de près de 10% l'efficacité maximale du procédé. / The aim of this research lies in the development and optimization of an advanced oxidation process dedicated to wastewater treatment. This process combines an atmospheric plasma discharge with an immobilized catalyst (Ag-DLC) in order to generate hydroxyl radicals. The novelty of this work relies in the use of a falling thin film system (thickness < 1mm) to promote the interaction between plasma discharge, the solution and the catalytic surface. The results demonstrate the interest of employing a Dielectric Barrier Discharge configuration instead of a Glidarc system. Indeed, the Glidarc system leads to a too high production level of NOx. The efficiencies of these two processes were respectively 0,67 g.kWh-1 and 0,23 g.kWh-1 for the DBD and the Glidarc system. Moreover, the influence of the counter-electrode / flowing plate material on the process efficiency was pointed out. Compared to stainless steel, the presence of corrosion reactions on the brass surface produced NO2- and thus, reduced the process efficiency. The production of HO??radicals at the plasma / water interface was estimated to be 4,4.10-9 molHO.J-1 by the use of 1D CFD modeling. Side reactions such as HO? recombinaison appear as the major limiting factors. Futhermore the model demonstrates that the time needed to complete all the chemical reactions was less than 1 µs and these reactions occurred only in the first 5% of the top liquid film. Finally with an Ag-DLC based catalytic coating elaborated by RF-PECVD, the efficiency was increased to 10% compared to the best efficiency observed with the non-catalytic system.

Dans ces travaux de recherche nous avons cherché à développer et optimiser un procédé d'oxydation avancée dédié au traitement de l'eau. Pour cela, nous avons développé un procédé combinant la génération de radicaux hydroxyles par décharge plasma à pression atmosphérique et un catalyseur fixe. L'originalité de ce travail a été d'utiliser un système d'écoulement permettant l'obtention d'un film liquide (épaisseur<1mm) afin de permettre une interaction entre le plasma, le liquide et la surface catalytique. Les résultats ont montré que l'utilisation d'une décharge de type DBD était préférable au Glidarc, thermique, énergétique et produisant trop de NOx. Les valeurs d'efficacité respectives étaient de 0,67g.kWh-1 pour la DBD et 0,23g.kWh-1 pour le Glidarc. Les résultats ont également montré que le matériau de la contre-électrode/plaque d'écoulement pouvait influencer l'efficacité du procédé. Ainsi dans le cas du laiton, la présence de réactions de corrosion produisant des NO2- réduit l'efficacité de moitié comparée à une plaque d'acier. Le développement d'un modèle CFD 1D à partir des données expérimentales du procédé a été effectué. Ce modèle a permis d'estimer que l'efficacité énergétique de production des HO est de 4,4.10-9 molHO.J-1 pour la décharge DBD. Il a également permis de montrer que le facteur majoritairement limitant était la diffusion moléculaire des espèces dans le liquide. En effet les réactions chimiques ne dure pas plus de 1µs après l'impact du streamer et n'ont lieu que dans les premiers 5% d'épaisseur de liquide. Enfin, un dépôt catalytique de DLC dopé à l'argent et réalisé par rf-PECVD a permis d'améliorer de près de 10% l'efficacité maximale du procédé.

Le travail réalisé dans le cadre de cette thèse porte sur l'effet que les bâtiments exercent sur l'atmosphère urbaine et notamment sur les échanges énergétiques qui s'opèrent entre les deux systèmes. Afin de modéliser plus finement les effets thermiques du bâtiment sur les écoulements atmosphériques lors de simulations réalisées par le logiciel de CFD Code_Saturne, nous procédons au couplage de cet outil avec le modèle de bâtiment BuildSysPro. Cette bibliothèque fonctionne sous Dymola et peut calculer des matrices descriptives du bâtiment utilisables ensuite en dehors du logiciel. Ce sont donc ces matrices qui sont utilisées pour le couplage par l'intermédiaire d'un code assurant l'échange de données entre les calculs de thermique du bâtiment et ceux de CFD. Après une revue des phénomènes physiques en lien avec l'atmosphère urbaine et des modèles existants, nous nous intéressons aux interactions entre l'atmosphère et le milieu urbain, notamment les bâtiments. Ceux-ci peuvent avoir un impact sur les écoulements aussi bien dynamique, en tant qu'obstacles, que thermique, via leurs températures de parois. Parallèlement à la mise en place du couplage entre les deux logiciels, nous étudions les données de la campagne de mesures EM2PAU que nous utilisons pour notre validation. EM2PAU, réalisée en 2011 à Nantes, représente une rue canyon idéalisée par deux rangées de conteneurs. La campagne a pour spécificité de prendre simultanément les mesures de températures d'air et de parois ainsi que les vitesses du vent de référence et des écoulements dans le canyon par un anémomètre sonique placé à 10 m d'altitude et six autres positionnés en six emplacements dans le canyon. Nous cherchons donc à mettre en évidence les effets dynamiques et thermiques des bâtiments sur les écoulements à partir des résultats de cette campagne, pour ensuite les simuler. Puis la modélisation numérique des écoulements sur le domaine de EM2PAU est réalisée. L'objectif de ce travail est de mettre en évidence l'influence des effets thermiques des parois sur les flux atmosphériques. Nous comparons des simulations avec différentes méthodes pour donner les valeurs des températures de surface des conteneurs. La première méthode consiste à imposer ces températures d'après les mesures ; ainsi la température de chaque paroi sera fixée à la température de surface mesurée lors de l'instrumentation de EM2PAU. Quant à la deuxième méthode, on impose la température de l'air extérieur mesurée à l'instant simulé à toutes les parois des conteneurs, afin de créer un cas où l'on n'observe que peu ou pas d'échanges de chaleur. Enfin la troisième méthode est la simulation couplée de Code_Saturne et BuildSysPro. Les résultats des différentes simulations sont alors comparés afin de distinguer les effets thermiques des parois des bâtiments sur les écoulements d'air. Nous observons que les effets dynamiques sont primordiaux et peuvent engendrer des vitesses verticales de l'écoulement dans le canyon de l'ordre plusieurs mètres par seconde, tandis que des écarts de températures de surface de l'ordre de 15°C peuvent modifier les vitesses verticales du vent de moins de 0, 5 mètres par seconde. Si ces effets thermiques sont difficiles à isoler sur des mesures en raison des autres phénomènes susceptibles d'influencer les écoulements atmosphériques, les études numériques peuvent toutefois mieux quantifier ces différences / This thesis work is about the effect of buildings on the urban atmosphere and more precisely the energetic exchanges that take place between these two systems. In order to model more finely the thermal effects of buildings on the atmospheric flows in simulations run under the CFD software Code_Saturne, we proceed to couple this tool with the building model BuildSysPro. This library is run under Dymola and can generate matrices describing the building thermal properties that can be used outside this software. In order to carry out the coupling, we use these matrices in a code that allows the building thermal calculations and the CFD to exchange their results. After a review about the physical phenomena and the existing models, we explain the interactions between the atmosphere and the urban elements, especially buildings. The latter can impact the air flows dynamically, as they act as obstacles, and thermally, through their surface temperatures. At first, we analyse the data obtained from the measurement campaign EM2PAU that we use in order to validate the coupled model. EM2PAU was carried out in Nantes in 2011 and represents a canyon street with two rows of four containers. Its distinctive feature lies in the simultaneous measurements of the air and wall temperatures as well as the wind speeds with anemometers located on a 10 m-high mast for the reference wind and on six locations in the canyon. This aims for studying the thermal influence of buildings on the air flows. Then the numerical simulations of the air flows in EM2PAU is carried out with different methods that allow us to calculate or impose the surface temperature we use, for each of the container walls. The first method consists in imposing their temperatures from the measurements. For each wall, we set the temperature to the surface temperature that was measured during the EM2PAU campaign. The second method involves imposing the outdoor air temperature that was measured at a given time to all the surfaces, reducing every heat exchange to almost zero. The third method at last is the coupled simulation of Code_Saturne and BuildSysPro where BuildSysPro calculates the wall temperature from the Code_Saturne data. . The results of these different ways of modelling the wall temperatures are then compared in order to show the thermal effects of building wall heating on the air flows. We notice that the dynamic effects are dominant and can generate vertical wind speed that can pass several meters per second. On the other hand, differences of surface temperatures higher than 15°C can influence the vertical wind speed for less than 0.5 meters per second. These thermal effects are not easily highlighted with measured data because of the other phenomena that can impact the air flows. However they can be quantified with numerical studies

Cette thèse porte sur l'influence d'un jet plan pariétal sur une nappe de retour (backlayering). Cette étude peut être divisée en deux parties : une étude expérimentale et une étude numérique. Les expériences ont permis d'identifier le comportement de la nappe de retour en présence du jet. L'étude numérique débute par une validation du modèle numérique à partir des résultats de l'étude expérimentale. Une étude paramétrique est réalisée afin de déterminer l'influence du jet sur les caractéristiques et la stratification de la nappe de retour. La présence du jet plan pariétal induit une diminution de la longueur de la nappe de retour et une variation de l'épaisseur de la nappe de retour du fait de l'introduction du jet :pour une vitesse débitante fixée, l'augmentation du rapport entre la vitesse d'éjection, du jet et la vitesse débitante induit une augmentation de l'épaisseur de la nappe de retour ; pour un rapport entre la vitesse d'éjection du jet et la vitesse débitante constant, l'augmentation de la vitesse débitante conduit à un affinement de la nappe de retour. La stratification de la nappe de retour est conservée en présence du jet plan pariétal / The aim of this study is to determine the influence of a wall attached plane jet on the backlayering. This study is composed of two parts. The experimental part consists in identifying the behaviour of the backlayering in presence of the jet. The numerical part started with a validation of the numerical model using experimental study results. Then, a parametric study enables to determine the influence of the jet on the backlayering. The presence of the parietal plan jet induces a decrease of the backlayering length. Two main conclusions are achieved for the backlayering thickness change due to the jet : for a constant air flow velocity, the increase of the ratio between jet velocity and air flow velocity induces an increase of the backlayering thickness ; for a constant ratio between jet velocity and air flow velocity, the increase of the air flow velocity leads to an increase of the backlayering thickness. Finally, the backlayering stratification is maintained on the presence of the jet

Les systèmes fixes par aspersion d'eau sont souvent cités comme un moyen pour améliorer le niveau de sécurité dans les tunnels routiers. Cependant, leurs effets sur un incendie en tunnel et les interactions possibles avec le système de ventilation sont encore mal connus. Le travail présenté ici vise à étudier l'action de ces systèmes par une approche phénoménologique de l'évaporation et l'utilisation d'un modèle CFD tridimensionnel. L'emploi d'un code CFD incluant des modèles diphasiques repose sur la simulation de deux situations expérimentales : une expérience académique de spray évaporant et des essais d'incendie avec aspersion dans un tunnel à échelle réduite. Les résultats obtenus sont relativement satisfaisants, en utilisant une modélisation Euler-Lagrange pour l'écoulement diphasique. L'influence du système d'aspersion sur la propagation de la couche de fumée est plus particulièrement étudiée. Cependant, ce type de modélisation nécessite des données d'entrée d'un niveau difficilement accessible sans avoir recours à une expérience de caractérisation des buses. Une approche plus phénoménologique, à intégrer dans un modèle unidimensionnel, est envisagée comme alternative.

Ce travail vise à modéliser la propagation du vent par des techniques de simulation de CFD dans des tissus urbains contrastés à l’échelle du quartier (450m × 450m) afin d'évaluer les impacts de la morphologie urbaine sur le potentiel énergétique éolien. Face aux problèmes environnementaux locaux et globaux, le développement des énergies renouvelables devient de plus en plus urgent à l’échelle mondiale. Notre objectif est d’adapter la morphologie urbaine au développement du captage de l’énergie éolienne. Pour ce faire, deux questions fondamentales sont posées : quels sont les impacts de la morphologie urbaine sur l'écoulement du vent dans la ville ? Comment évaluer le potentiel énergétique éolien de la ville ? Nous partons d’une synthèse bibliographie sur les trois domaines : l'énergie éolienne urbaine, la morphologie urbaine et la simulation aéraulique. Les avantages et les inconvénients des éoliennes urbaines sont discutés. Des indicateurs de morphologie urbaine sont ainsi proposés. Les méthodes de simulation CFD (Mécanique des fluides numérique) sur l'écoulement du vent autour de bâtiments sont présentées. Des nouvelles technologies et des méthodes d'intégration architecturale de petites éoliennes sont ainsi présentées. Pour valider notre modèle, nous faisons une étude approfondie des paramètres du code ANSYS FLUENT avec la référence des données d’un benchmark issu d'expérimentations en soufflerie. Plusieurs méthodes de vérification avant et après les simulations sont proposées. Après validation de notre modèle, des simulations de formes simples sont faites. Avec un bâtiment, deux bâtiments et trois bâtiments, les impacts de différents paramètres simples de configuration sont analysés. En raison de leur efficacité et de leur faisabilité, l'évaluation énergétique est principalement faite au-dessus des toits de bâtiments. Les effets de la concentration du vent au-dessus des toits sont ainsi évalués. Nous simulons ensuite des modèles morphologiques plus complexes, par l’évaluation du potentiel éolien dans six quartiers de formes urbaines typiques dans le monde et favorables pour exploiter l'énergie éolienne. Treize indicateurs morphologiques ayant potentiellement une influence sur l'écoulement du vent sont proposés. Des analyses de cohérence entre les indicateurs morphologiques et les indicateurs d'évaluation de l'énergie éolienne au-dessus des toits sont faites. Ensuite, des variantes de formes urbaines typiques sont étudiées. L’influence de certaines modifications de ces prototypes sur le potentiel éolien est étudiée et les prototypes les plus favorables pour l'exploitation de l'énergie éolienne sont proposés. L'énergie éolienne dans un quartier existant à Pekin est ensuite évaluée. Des conditions locales (formes physiques, climat, aspects socio-économiques) sont considérées pour évaluer la faisabilité du développement de l’énergie éolienne. L'emplacement et le choix des techniques d’éoliennes urbaines adaptées sont discutés. Enfin, une méthode mathématique basée sur des paramètres morphologiques simples est présentée pour évaluer le potentiel éolien au-dessus des toits. L'indicateur de l'impact d'obstacle des bâtiments est défini. La cohérence entre les résultats de la méthode mathématique et ceux de la simulation CFD est enfin analysée. / This work concentrates on wind flow simulation by CFD techniques in different urban fabrics at the neighborhood scale (450m × 450m), in order to assess the impacts of urban morphology on wind energy potential. Facing local and global environmental problems, the development of renewable energy is becoming increasingly urgent and important. Our goal is to adapt the urban morphology in the development of urban wind energy. To do so, two fundamental questions are aroused: What are the impacts of urban morphology on the wind flow in the city? How to evaluate the urban wind power? We start with a summary of bibliography on three areas: urban wind energy, urban morphology and CFD (Computational fluid dynamics) simulation. The advantages and disadvantages of different urban wind turbines are discussed. Some new technologies and methods of architectural integration of small wind turbine are presented. Indicators of urban morphology are proposed. The methods of CFD simulation of wind flow around buildings are presented. To validate our model of air flow simulation, we start with a study of the numerical simulation parameters with the reference of a benchmark of wind tunnel experiments. Several methods of verification before and after the CFD simulation are presented. After validation of our model, we begin to simulate some simple forms. With one, two and three buildings, the impacts of different simple configurations parameters are analyzed. Considering the effectiveness and feasibility, wind energy assessment is mainly taken over the roofs of buildings. The combined effects of altitude and wind concentration above the roof, are evaluated. We then simulate more complex morphological models, by assessing the wind potential in six neighborhoods that are urban forms both worldwide typical and favorable for wind energy development. Morphological indicators that potentially affect the wind flow are proposed. Analysis of coherence are made between these indicators and the values of the coefficient of increase of wind speed and of the wind power indicator above the roof. Then, some prototypes of basic urban forms from real neighborhoods are studied. The influence of some modifications of these prototypes on the wind energy potential is analyzed, and the most favorable conditions for the exploitation of wind power prototypes are recommended. Then, wind energy in a real neighborhood in Beijing is evaluated. Local conditions (physical forms, climate, social-economic aspects) are considered to assess the feasibility of developing wind energy. The location and the choice of techniques suitable urban wind turbines are discussed. Finally, a mathematical method based on some simple morphological parameters is presented, in order to assess wind potential above the rooftops. The obstacle impact indicator of the buildings is defined. The coherence between the results of the simplified mathematical method and those of CFD simulation are analyzed.

Le contrôle précis des modèles géométriques joue un rôle important dans de nombreux domaines. Pour l’optimisation de forme en CFD, le choix des paramètres de contrôle et la technique de déformation de forme est critique. Nous proposons un modeleur paramétrique avec une nouvelle méthode de déformation d’objets, ayant pour objectif d’être intégré dans une boucle d’optimisation automatique de forme avec un solveur CFD. Notre méthodologie est basée sur une double paramétrisation des objets : géométrique et architecturale. L’approche géométrique consiste à décrire les formes par un squelette, composé d’une famille de courbes B-Splines, appelées courbes génératrice et courbes de section. Le squelette est paramétré avec une approche architecturale. Au lieu d’utiliser les points de contrôle de la représentation classique par courbes B-Splines, la géométrie est contrôlée par ces paramètres architecturaux. Cela permet de réduire considérablement le nombre de degrés de liberté utilisés dans le problème d’optimisation de forme, et permet de maintenir une description haut niveau des objets. Notre technique intègre un contrôle de forme et un contrôle de régularité, permettant d’assurer la génération de nouvelles formes valides et réalistes. Les déformations de la géométrie sont réalisées en posant un problème inverse : déterminer une géométrie correspondant à un jeu de paramètres cibles. Enfin, une technique de reconstruction de surface est proposée. Nous illustrons le modeleur paramétrique développé et intégré dans une boucle d’optimisation automatique de forme sur trois cas : un profil d’aile d’avion, un foil AC45 d’un voilier de course et un bulbe de chalutier de pêche. / The precise control of geometric models plays an important role in many domains. For shape optimisation in CFD, the choice of control parameters and the way to deform a shape are critical. In this thesis, we propose a new approach to shape deformation for parametric modellers with the purpose of being integrated into an automatic shape optimisation loop with a CFD solver. Our methodology is based on a twofold parameterisation: geometrical and architectural. The geometrical approach consist of a skeleton-based representation of object. The skeleton is made of a family of B-Spline curves, called generating curve and section curves. The skeleton is parametrised with an architectural approach: meaningful design parameters are chosen on the studied object. Thus, instead of using the control points of a classical B-spline representation, we control the geometry in terms of architectural parameters. This reduce the number of degrees of freedom and maintain a high level description of shapes. We ensure to generate valid shapes with a strong shape consistency control based on architectural considerations. Deformations of the geometry are performed by solving optimisation problems on the skeleton. Finally, a surface reconstruction method is proposed to evaluate the shape’s performances with CFD solvers. We illustrate the parametric modeller capabilities on three problems, performed with an automatic shape optimisation loop: the wind section of an plane (airfoil), the foil of an AC45 racing sail boat and the bulbous bow of a fishing trawler.

L’utilisation de techniques numériques lors de la conception d’un produit s’est largement généralisée au cours des 30 dernières années. Pourtant, la lenteur des calculs et la spécialisation des modèles numériques restent problématiques.Nous avons donc choisi de développer des modèles réduits du comportement de balayage sur un moteur Diesel 2-temps à lumières. Ces modèles sont analytiques, génériques, rapides d’utilisation et permettent d’éliminer les problématiques de traitement numérique. Ils sont aussi des instruments performants dans la recherche de solutions de conception. Une modélisation CFD 2D a tout d’abord été développée pour servir de bases de données, avec la définition des paramètres primordiaux à suivre pour quantifier un balayage optimal.Le travail de recherche a dévoilé une méthodologie nouvelle fondée sur un méta-modèle du comportement dit « neuro-séparé » comprenant un modèle neuronal d’état, un modèle neuronal pseudo-dynamique et un modèle à variables séparées. Ensuite, un processus d'aide à la décision exploitant les modèles précédents a été mis en place au travers d’un processus d’optimisation évolutionnaire (fondé sur les algorithmes génétiques) puis de la simulation comportementale rapide des solutions optimales de conception par un krigeage.La démarche de conception multipoints de vue, multi-critères et multi-physiques appliquée au moteur intègre aussi une dimension cognitive : l’exploration évolutionnaire des espaces de solutions a été menée de façon libre et forcée. Afin de valider notre approche, nous avons mis en place des critères de qualification appliqués à chacun de nos modèles, permettant de quantifier les écarts visà-vis de la base initiale CFD qui a fondé nos modèles réduits.Notre démarche a mené à la création d’un outil d’aide à la modélisation et à la décision exploitant les modules Python et Matlab développés. / The use of numerical methods to design a product became more and more commonover the past 30 years. However, numerical models are still specialized and they do not run fastwhich make their use problematic. So some reduced models of scavenging have been developed. These models are analytical andgeneric; they run quickly and avoid the numerical treatment problems. They are also some efficienttools in the search of design solutions.The work carried out has led to a new methodology based on a behavioral meta-model called“neuro-separated” including a neuronal model of state, a pseudo-dynamic neuronal model and amodel with separated variables. Then, a process of decision aids exploiting the models previouslydeveloped in evolutionary algorithms (genetic algorithms) and the fast behavioral simulation of theoptimal design solutions thanks to the kriging approach.This design approach is multi-viewpoints, multi-criteria and multi-physics. It also includes acognitive dimension: both free and controlled evolutionary explorations of solution spaces have beendone. To validate the method, some qualification criteria have been evaluated for each model. Theyallow to understand and to assume the gap between the reduced models and the initial CFD base(where the model are coming from). Our approach has led to the development of a tool of model and decision aids using Python and Matlab software programs.

L’évaluation du débit en réseaux hydrauliques à surface libre est une problématique actuelle sur le plan scientifique, à forts enjeux technologiques, économiques et écologiques. Dans cette thèse, de nouvelles méthodologies d’instrumentation, basées sur une synergie entre mesures non intrusives de hauteur d’eau et modélisation numérique ont été développées. Celles-ci s’appliquent d’une part à des collecteurs dont le fonctionnement hydraulique est complexe et, d’autre part, à des ouvrages non-standard (Venturi, déversoirs d’orage). Ce travail de thèse multidisciplinaire vise une meilleure compréhension de l’écoulement pour en déduire des relations Q=f(hi) plus robustes, spécifiques à chaque site et associées à leurs incertitudes; mais également l’identification de possibles modifications du site de mesure afin d’améliorer l’estimation du débit. Au final, l’applicabilité des méthodologies développées a été éprouvée au travers de plusieurs études sur sites réels. / The evaluation of the flow rate in free surface water systems is a current scientific problem, related to high technological, economical and ecological issues. In this study, new methods of instrumentation based on a synergy between non-intrusive water level measurements and numerical modeling have been developed. These methods are applied first to sewer pipes with complex hydraulic conditions then to non-standard hydraulic structures (Venturi flumes, Combined Sewer Overflows). This multidisciplinary work aims at a better understanding of the flow to identify more robust site-specific Q=f(hi) relationships related to their uncertainties. It also aims at the identification of possible modification of the measurement site in order to improve the flow rate evaluation. Finally, the applicability of the developed methodologies has been tested through several real site studies.