La qualité des services écosystémiques qu'offrent les lacs est liée à la structure et au fonctionnement de leur écosystème. Protéger leur masse d'eau est devenu un objectif global qui requiert une meilleure compréhension de leur fonctionnement, un suivi et une évaluation de leur qualité. Expliquer les hétérogénéités spatio-temporelles des variables physico-chimiques et du phytoplancton est un problème récurrent rencontré en écologie et en hydrobiologie. Comprendre la dynamique de ces hétérogénéités est aussi un prérequis essentiel pour évaluer, protéger et restaurer objectivement les écosystèmes lacustres. En ce qui concerne la surveillance, les hétérogénéités spatio-temporelles introduisent des incertitudes sur la représentativité des mesures par rapport à l'entièreté de la masse d'eau qui est donc discutable et doit être vérifiée. En Europe, la directive cadre sur l'eau (DCE) initiée en 2000 définit un cadre pour la gestion et la protection des eaux. La classification des masses d'eau en fonction de leur état écologique est un point important dans l'implémentation de cette directive. Pour les lacs et les retenues, l'évaluation de cet état écologique est basée sur des paramètres biologiques, physico-chimiques et hydromorphologiques. Les indicateurs liés au phytoplancton et aux paramètres physico-chimiques sont calculés à partir de quatre prélèvements réalisés pendant la période d'activité biologique pour une année sur un plan de gestion de six ans. Dans ce contexte, la modélisation tridimensionnelle (3D) et la prise en compte des forçages qui conduisent aux hétérogénéités spatio-temporelles est une condition préalable nécessaire tant en limnologie appliquée que théorique. Cette thèse aborde la complexité du fonctionnement des lacs et la capacité des modèles 3D à reproduire leur fonctionnement. L'apport de la modélisation 3D est présenté i) pour la compréhension du fonctionnement de lacs de différentes tailles, ii) couplée aux observations satellitaires, pour l'étude de l'influence des forçages par le vent et de l'hydrodynamique sur l'abondance et la distribution spatiale de phytoplancton, iii) dans le cadre de la DCE, pour l'évaluation des incertitudes d'une évaluation de l'état écologique d'un plan d'eau. Pour cela, deux modèles 3D ont été créés et analysés, un pour le lac de Créteil (42~ha) et un autre pour le Léman (580~km$^2$). Celui du lac de Créteil a été validé à partir de données à hautes fréquences acquises en trois points du lac. Il reproduit correctement son hydrodynamique complexe, sa structure thermique, l'alternance entre les périodes de stratification et les épisodes de mélange, ainsi que les ondes internes. Le modèle du Léman a été validé en utilisant des données mensuelles et bimensuelles en deux stations de prélèvement du lac. Il reproduit aussi correctement son hydrodynamique et la variabilité saisonnière de paramètres physico-chimiques et biologiques. Les résultats des simulations mettent en avant les mécanismes physiques et hydrodynamiques responsables de l'apparition de sites où la biomasse de phytoplancton observée est plus élevée. Dans le cadre de la DCE, ces résultats montrent aussi une variabilité spatiale importante des sous-états écologiques basés sur les différents paramètres qui dépendent du choix des dates des campagnes de mesure et du point d'échantillonnage. Ces résultats ont aussi été utilisés pour estimer la représentativité d'une station de prélèvement. Les résultats de cette thèse i) confirment que le fonctionnement des plans d'eau de toute taille est complexe et que les processus physiques génèrent des hétérogénéités spatio-temporelles, ii) suggèrent que le vent et l'hydrodynamique influencent significativement l'abondance et la distribution spatiale du phytoplancton et que iii) ces hétérogénéités peuvent biaiser notre estimation du statut écologique des plans d'eau dans le cadre de la DCE / The quality of ecosystem services provided by lakes is related to the ecosystem structure and functioning. Protecting water bodies is therefore a global goal that requires a better understanding of their function, a monitoring and a water quality assessment. Explaining spatio-temporal heterogeneities of physico-chemical parameters and phytoplankton has been a recurrent ecological and hydrobiological issue. Understanding the dynamics of these heterogeneities is an essential prerequisite for objectively assessing, protecting and restoring freshwater ecosystems. Moreover, three-dimensional (3D) and taking into account and drivers of these heterogeneities are essential prerequisites for theoretical and applied limnology. Concerning the monitoring, spatio-temporal heterogeneities are responsible of uncertainties on the representativeness of the data versus the whole lake which might be questionable and needs to be verified. In Europe, the Water Framework Directive (WFD) initiated in 2000 defines a framework for managing and protecting water bodies in Europe. The classification of water bodies into ecological status is a key issue for the implementation of that framework. For lakes and reservoirs, the assessment of this status is based on biological, physico-chemical and hydro-morphological indicators. Physico-chemical and phytoplankton indicators are calculated based on four observations at an unique sampling station over the growing season, this evaluation being assessed one year for a six-year management plan. In this context, this thesis focuses on the complexity of lakes functioning and the capability of three-dimensional (3D) models to reproduce their functioning. The contribution of 3D models is presented i) for understanding the functioning of lakes of different sizes, ii) coupled to satellite observations, for studying of the influence of wind forcing and hydrodynamics on phytoplankton abundance and spatial heterogeneities, iii) in the context of the WFD, for assessing uncertainties in the lake ecological status assessment. To do that, two 3D models have been created and analyzed, one for Lake Créteil (42 ha) and another for Lake Geneva (580 km2). Lake Créteil 3D model was validated by using high frequency data recorded at three stations. It reproduces well the complex hydrodynamic functioning of the lake, its thermal structure, the alternation between thermal stratification episodes and mixing events, and internal waves. Lake Geneva 3D model was validated by using monthly and bimonthly data at two stations. It reproduces also properly the hydrodynamic functioning of the lake and the seasonal variability of biological and physico-chemical parameters. Simulation results highlight physical and hydrodynamic mechanisms responsible for the occurrence of seasonal hot-spots in phytoplankton abundance. In the context of the WFD, simulation results show also a strong spatial variability of lake ecological status depending on the timing of the four sampling dates as well as the location of the sampling station. These results were also used to assess to representativeness of sampling stations. The results of this thesis suggest that i) the functioning of lakes of different sizes is complex and physical processes generates spatio-temporal heterogeneities, ii) wind and hydrodynamics influence the abundance and the spatial distribution of phytoplankton et iii) spatio-temporal heterogeneities can bias our evaluation of lake ecological status in the WFD

Afin de cibler une meilleure rentabilité des gisements d'hydrocarbures, la simulation numérique présente un intérêt grandissant dans le secteur pétrolier. Dans ce contexte, deux principaux types d'applications sont distingués : l'ingénierie du réservoir, où les modèles géologiques sont définis comme statiques et l'exploration des gisements par la simulation de la formation des hydrocarbures. Cette dernière application nécessite des modèles de bassins dynamiques. Pour ce type de simulation, la modélisation mathématique et numérique a connu une avancée importante durant les dernières années. En revanche, la construction de maillages indispensables pour les simulations reste une tâche lourde dans le cas de géométries complexes et dynamiques. La difficulté se traduit par la particularité du domaine à mailler, qui est dictée par la mécanique des milieux granulaires, c'est-à-dire des milieux quasi-incompressibles et hétérogènes. De plus, les données sismiques sont des surfaces non-implicites modélisant des blocs volumiques. Dans ce cadre, nous nous intéressons à la génération de maillages lagrangiens hexa-dominants des bassins à géométrie complexe. Le maillage souhaité doit contenir une couche de mailles principalement hexaédriques entre deux horizons et respecter les surfaces failles traversant ces horizons. Une méthodologie originale basée sur la construction d’un espace déplié est proposée. Le maillage souhaité est alors traduit par une grille 3D contrainte dans cet espace. Plusieurs techniques d’optimisation de maillages sont aussi proposées / To target more profitable oil and gas deposits, the numerical simulation is of growing interest in the oil sector. In this context, two main types of applications are distinguished: reservoir engineering where geological models are defined as static and exploration for the simulation of hydrocarbons formation. The latter application requires dynamic basins models. For this type of simulation, mathematical and numerical modeling has been an important advance in recent years. However, the construction of meshes needed for the simulations is a difficult task in the case of complex and dynamic geometries. The difficulty is reflected by the characteristic of the domain to mesh, which is defined from the mechanics of granular media which are almost incompressible and heterogeneous environments. In addition, the seismic data represent non-implicit surfaces modeling volume blocks. In this context, we focus on Lagrangian hex-dominant mesh generation of basins with complex geometry. The desired mesh must contain a layer of almost hexahedral meshes between two horizons and conform to fault surfaces through these horizons. A novel methodology based on the construction of an unfolded space is introduced. The desired mesh is then seen as a constrained 3D grid in this novel space. Several mesh optimization techniques are also proposed

La mise en registre 3D (opération parfois appelée alignement) est un processus de transformation d’ensembles de données 3D dans un même système de coordonnées afin d’en aligner les éléments communs. Deux ensembles de données alignés ensemble peuvent être les scans partiels des deux vues différentes d’un même objet. Ils peuvent aussi être deux modèles complets, générés à des moments différents, d’un même objet ou de deux objets distincts. En fonction des ensembles de données à traiter, les méthodes d’alignement sont classées en mise en registre rigide ou non-rigide. Dans le cas de la mise en registre rigide, les données sont généralement acquises à partir d’objets rigides. Le processus de mise en registre peut être accompli en trouvant une seule transformation rigide globale (rotation, translation) pour aligner l’ensemble de données source avec l’ensemble de données cible. Toutefois, dans le cas non-rigide, où les données sont acquises à partir d’objets déformables, le processus de mise en registre est plus difficile parce qu’il est important de trouver à la fois une transformation globale et des déformations locales. Dans cette thèse, trois méthodes sont proposées pour résoudre le problème de mise en registre non-rigide entre deux ensembles de données (représentées par des maillages triangulaires) acquises à partir d’objets déformables. La première méthode permet de mettre en registre deux surfaces se chevauchant partiellement. La méthode surmonte les limitations des méthodes antérieures pour trouver une grande déformation globale entre deux surfaces. Cependant, cette méthode est limitée aux petites déformations locales sur la surface afin de valider le descripteur utilisé. La seconde méthode est s’appuie sur le cadre de la première et est appliquée à des données pour lesquelles la déformation entre les deux surfaces est composée à la fois d’une grande déformation globale et de petites déformations locales. La troisième méthode, qui se base sur les deux autres méthodes, est proposée pour la mise en registre d’ensembles de données qui sont plus complexes. Bien que la qualité que elle fournit n’est pas aussi bonne que la seconde méthode, son temps de calcul est accéléré d’environ quatre fois parce que le nombre de paramètres optimisés est réduit de moitié. L’efficacité des trois méthodes repose sur des stratégies via lesquelles les correspondances sont déterminées correctement et le modèle de déformation est exploité judicieusement. Ces méthodes sont mises en oeuvre et comparées avec d’autres méthodes sur diverses données afin d’évaluer leur robustesse pour résoudre le problème de mise en registre non-rigide. Les méthodes proposées sont des solutions prometteuses qui peuvent être appliquées dans des applications telles que la mise en registre non-rigide de vues multiples, la reconstruction 3D dynamique, l’animation 3D ou la recherche de modèles 3D dans des banques de données. / Three-dimensional registration (sometimes referred to as alignment or matching) is the process of transforming many 3D data sets into the same coordinate system so as to align overlapping components of these data sets. Two data sets aligned together can be two partial scans from two different views of the same object. They can also be two complete models of an object generated at different times or even from two distinct objects. Depending on the generated data sets, the registration methods are classified into rigid registration or non-rigid registration. In the case of rigid registration, the data is usually acquired from rigid objects. The registration process can be accomplished by finding a single global rigid transformation (rotation, translation) to align the source data set with the target data set. However, in the non-rigid case, in which data is acquired from deformable objects, the registration process is more challenging since it is important to solve for both the global transformation and local deformations. In this thesis, three methods are proposed to solve the non-rigid registration problem between two data sets (presented in triangle meshes) acquired from deformable objects. The first method registers two partially overlapping surfaces. This method overcomes some limitations of previous methods to solve large global deformations between two surfaces. However, the method is restricted to small local deformations on the surface in order to validate the descriptor used. The second method is developed from the framework of the first method and is applied to data for which the deformation between the two surfaces consists of both large global deformation and small local deformations. The third method, which exploits both the first and second method, is proposed to solve more challenging data sets. Although the quality of alignment that is achieved is not as good as the second method, its computation time is accelerated approximately four times since the number of optimized parameters is reduced by half. The efficiency of the three methods is the result of the strategies in which correspondences are correctly determined and the deformation model is adequately exploited. These proposed methods are implemented and compared with other methods on various types of data to evaluate their robustness in handling the non-rigid registration problem. The proposed methods are also promising solutions that can be applied in applications such as non-rigid registration of multiple views, 3D dynamic reconstruction, 3D animation or 3D model retrieval.

Les alliages d'aluminium utilisés dans les structures aéronautiques (fuselage, voilure) sont soumis à des chargements cycliques, faisant de la fatigue l'un des facteurs dimensionnant. Dans cette thèse, les mécanismes d'amorçage de ces fissures de fatigue – au niveau des particules intermétalliques – et de micropropagation sont étudiés expérimentalement et numériquement sur les alliages 2050-T8 et 7050-T7451. Les analyses des premiers chapitres portent sur la description des particules intermétalliques qui sont les plus susceptibles de donner lieu à une amorce de fissure dans la matrice d'aluminium. Les effets de la nature des particules et de leur taille sont quantifiés. La proximité entre les particules intermétalliques et les pores y est décrite. Cette analyse expérimentale fait intervenir des observations de surface en microscopie électronique à balayage (MEB), ainsi que des caractérisations tridimensionnelles (3D) réalisées à l'aide de la tomographie par rayonnement synchrotron. Dans les chapitres suivants, les analyses traitent de la prévision des chemins de fissuration à l'échelle d'un grain. Elles s'appuient sur des observations expérimentales de surface et des simulations de plasticité cristalline 3D pour comprendre les chemins de fissuration. Un modèle d'endommagement tenant compte de la cristallographie est alors proposé pour simuler la propagation de fissure par éléments finis. / Aluminium alloys used for aerospace structures (wing, fuselage) are subjected to cyclic loading. Fatigue properties of such alloys are therefore taken into account for the design of such parts. In this thesis, initiation mechanisms of fatigue cracks – near intermetallic particles – and micropropagation are studied experimentally and numerically on alloys 2050-T8 and 7050-T7451. In the first chapters, the analysis focuses on intermetallic particles which are most prone to initiate a fatigue crack in the aluminium matrix. The effects of the nature of particles as well as their size are quantified. The proximity between intermetallic particles and pores is described. This experimental analysis use surface observations obtained with a scanning electron microscope (SEM), and three-dimensional (3D) characterizations using synchrotron tomography. In the last chapters, analysis are oriented towards the prediction of crack paths at the grain size. They rely on surface experimental observations and 3D crystal plasticity modelling in order to understand crack paths. A damage model taking into account crystallography is proposed to simulate crack propagation using the finite element method.

L'objectif de cette thèse est, l'étude fine des effets d'amplification suite à la géologie locale, notamment dans des bassins sédimentaires, ainsi que l'évaluation des avantages et Iimites des méthodes d'estimations. La partie expérimentale est basée sur des données acquises au site test EURO-SEISTEST, le premier site-test en Europe occidentale proche de Thessalonique (Grèce). Après une évaluation des effets d'amplification au site même et une étude de leur stabilité en fonction de plusieurs paramètres, les estimations obtenues par cinq méthodes différentes (la méthode des rapports spectraux classique, la méthode des fonctions récepteurs, une inversion généralisée, une méthode de la Coda et la méthode de Nakamura) sont comparées entre elles. Les avantages et limites de chaque méthode sont évaluées. Les résultats sont dans la suite confrontés et comparés avec des modélisations uni - et bidimensionnelles. La variabilité spatiale à travers la vallée examinée est étudiée par la cohérence spatiale. La grande couverture spatiale par des stations installées au site a permis de séparer deux zones d'un comportement différent de la cohérence entre deux signaux dépendant de la distance qui les sépare. Ceci n'a jamais été remarqué dans des études précédents à cause d 'une couverture spatiale trop faible. Finalement, l'étude paramétrique des effets d'amplifications dans un milieu hétérogène en trois dimensions montre que l'effet géométrique est important et peut provoquer des amplifications significatives. La réponse d'un milieu plus complexe n'est difficilement à prévoir, dépend de l 'azimut de l'onde incidente et montre des amplifications importantes en présence d'un forte contraste d'impedence présent.

L'objectif principal de cette thèse est de développer des outils pour la reconstruction de l'environnement urbain à partir d'images. Les entrées typiques de notre travail est un ensemble d'images de façades, des empreintes au sol de bâtiments, et des modèles 3D reconstruits à partir d'images aériennes. Les principales étapes comprennent le calibrage des images,le recalage avec le modèle 3D, la récupération des informations de profondeur ainsi que la sémantique des façades.Pour atteindre cet objectif, nous utilisons des techniques du domaine de vision par ordinateur, reconnaissance de formes et de l'informatique graphique. Les contributions de notre approche sont présentés en deux parties.Dans la première partie, nous nous sommes concentrés sur des techniques de reconstruction multi-vues dans le but de récupérer automatiquement les informations de profondeur de façades à partir un ensemble des photographies non calibrées. Tout d'abord, nous utilisons la technique structure et mouvement pour calibrer automatiquement l'ensemble des photographies. Ensuite, nous proposons des techniques pour le recalage de la reconstruction avec un modèle 3D. Enfin, nous proposons des techniques de reconstruction 3d dense (stéréo multi-vues et voxel coloring) pour produire un maillage 3D texturé d'une scène d'un ensemble d'images calibrées.La deuxième partie est consacrée à la reconstruction à partir d'une seule vue et son objectif est de récupérer la structure sémantique d'une façade d'une image ortho-rectifiée. La nouveauté de cette approche est l'utilisation d'une grammaire stochastique décrivant un style architectural comme modèle pour la reconstruction de façades. nous combinons un ensemble de détecteurs image avec une méthode d'optimisation globale stochastique en utilisant l'algorithme Metropolis-Hastings.

Le projet de recherche doctorale porte sur la cartographie, modélisation tridimensionnelle, stratigraphie et lithostratigraphie des sédiments de surface des bassins versants des rivières Saint-François et Chaudière. Ce projet s'inscrivait dans une série de projets initiés par le Ministère du Développement Durable, des Parcs et de l'Environnement du Québec dont l'objectif ultime était de mieux connaître et protéger les eaux souterraines des bassins versants des rivières Chaudière et Saint-François, de la frontière internationale à sa confluence avec le fleuve Saint-Laurent, dans la région de la capitale nationale. Ainsi, la présente recherche constitue un projet d'étude sur la dynamique quaternaire du Sud-Est du Québec qui permet de comprendre la nature et la distribution des sédiments et des séquences quaternaires. Dans ce contexte, un des intrants majeurs dans la détermination des propriétés géologiques régionales est la compréhension de la distribution et de l'architecture tridimensionnelle des sédiments de surface. La compilation et la cartographie des formations superficielles de ces deux bassins versants ont été des objectifs incontournables et requis pour l'analyse tridimensionnelle. Dans les bassins considérés, des dizaines de coupes naturelles ont été analysées et compilées ainsi qu'un grand nombre d'exploitations anthropiques dont les informations s'ajoutent à des milliers de forages de sources diverses. Leur intégration à la stratigraphie existante a permis l'élaboration des cartes d'épaisseurs de dépôts meubles et de la construction d'un modèle tridimensionnel des unités quaternaires pour l'ensemble des deux bassins versants. Le développement du cadre chronologique vient compléter la thèse. Un outil de corrélation stratigraphique basé sur l'âge des dépôts est proposé dans cette étude. Ce programme a pour but de corréler et de distinguer les différentes unités glaciolacustres des Appalaches. Cette nouvelle approche est fondée sur la valeur des rapports Ln/Tn, qui sont l'expression de la luminescence naturelle sur la luminescence induite après irradiation en laboratoire. Malgré une remise à zéro partielle des sédiments, les âges sur les varves de la fin du Sangamonien (85 ± 2 ka), du Wisconsinien moyen (56 ± 6 ka) et du Wisconsinien supérieur (< 24 ± 2 ka) ont permis de situer les séquences glaciolacustres à l'intérieur des cycles connus (stades, interstades et interglaciaires) qui sont eux-mêmes basés sur des axiomes (stratigraphie isotopique océanique et stratigraphie continentale). Un autre objectif principal de cette étude visait à définir les principaux événements de la paléogéographie de cette région. La complexité de la séquence stratigraphique reflète des conditions uniques d'englaciation et de déglaciation dans un secteur caractérisé par des retraits glaciaires à contre-pente, permettant le développement de lacs proglaciaires et la préservation d'un volume de sédiments quaternaires antérieurs à la dernière glaciation. Ces sédiments sont essentiellement associés à l'épisode glaciolacustre de la Formation de Gayhurst. Un travail de synthèse sur l'extension et la caractérisation des sédiments associés à ce lac proglaciaire (Wisconsinien moyen) est proposé afin d'évaluer son extension dans les deux bassins versants. L'âge de 45 ± 5 ka obtenus sur un faciès deltaïque, dans le cadre du programme de datation, indique que cette séquence s'est déposée lors du Wisconsinien Moyen. ______________________________________________________________________________ MOTS-CLÉS DE L’AUTEUR : Cartographie, stratigraphique quaternaire des Appalaches, Basses-Terres du Saint-Laurent, Modélisation géologique 3D, Inlandsis laurentidien, séquences et faciès glaciolacustres, géochronologie.

La modélisation à l’aide de l’approche RANS (Reynolds Averaged Navier Stokes) a été menée en trois dimensions, en considérant les régimes permanent et non permanent, dans le but de simuler les écoulements au niveau des jonctions, bifurcations et déversoirs d’orage. A travers ces trois exemples d’étude, plusieurs stratégies de modélisation (différentes combinaison lois de paroi/modèles de turbulence, différents algorithmes de couplage pression/vitesse, différents schémas de discrétisation spatiale, différentes conditions aux limites, différents types et tailles de maille, etc.) ont été testées et évaluées à l’aide de plusieurs indicateurs de performance (de type RMS –Root Mean Square) en s’appuyant sur les données de vitesses (vitesses moyennes in situ et champs de vitesses obtenus par PIV en laboratoire), hauteurs d’eau, débits (répartition de débits en bifurcation en laboratoire ou débits déversés in situ). Les résultats obtenus sont transposables aux autres cas de singularités et ouvrages spéciaux rencontrés en réseau d’assainissement et montrent que : i) les résultats des simulations 3D sont sensibles à la rugosité, aux conditions limites de hauteur et de vitesse ; ii) dans les trois cas d’étude, les schémas de discrétisation du second ordre et l’algorithme de couplage pression/vitesse PISO sont appropriés ; ii) la loi de paroi scalable couplée aux modèles de turbulence de type k-ε pour le cas des jonctions (avec un débit latéral inférieur ou égal au débit principal) ou des déversoirs semblent convenir, tandis que le modèle de turbulence RSM associé à la loi de paroi enhanced ou scalable permet de mieux représenter les écoulements à travers les bifurcations ou au niveau des jonctions lorsque le débit latéral est dominant. Sur la base de ces résultats, une méthodologie de modélisation plus générale définie en six étapes et fondée sur le guide proposée par Jakeman et al. (2006) a été mise au point. La méthodologie ainsi définie a été utilisée pour i) améliorer l’instrumentation du déversoir OTHU (Observatoire de Terrain en Hydrologie Urbaine) situé à Ecully, à partir de la simulation de sa courbe de fonctionnement et en fournissant les incertitudes sur les débits déversés obtenus ; ii) simuler l’implémentation des capteurs débitmétriques à l’aval d’une jonction. Elle a permis de concevoir et de dimensionner un nouveau dispositif de maîtrise des flux d’eau et de polluants déversés (technologie DSM – dispositif de surveillance et de maîtrise des flux déversés au milieu naturel). Ce dispositif a fait l’objet d’un dépôt de brevet international. Enfin, la mise en œuvre de cette méthodologie a été à l’origine de la découverte d’une nouvelle structure d’écoulement dans la branche latérale d’une bifurcation à 90°. L’analyse des résultats des simulations des écoulements mettant en évidence cette nouvelle structure a montré qu’il était possible de prédire l’apparition de cette dernière à partir du rapport d’aspect et du nombre de Froude. / The understanding of sewer flows behaviour is a key component for better urban drainage monitoring and management. However, these flows are conveyed across singularities (such as bends, drops, deviations, etc.) and special structures (combined sewer overflows –CSOs–, channels junction, dividing flow structures, etc.). These singularities and specific structures exhibit complex geometries, leading to open channel turbulent, three-dimensional and multiphase (pollutants and storm and sewer waters) flows. Using three-dimensional CFD (Computational Fluid Dynamics) platform allows a better understanding of mechanisms of contaminants transport through these structures and singularities, leading to a better sewer monitoring. In this thesis, 3D-RANS (Reynolds Averaged Navier Stokes) modelling approach under steady-state conditions is used in order to study flows within CSOs, junctions and bifurcations. Through these three structures, several modelling strategies (wall law/turbulence model combination, velocity/pressure coupling algorithm, spatial discretisation schemes, boundary conditions, computational mesh –shape and size of cells–, etc.) are tested and evaluated thanks to performance indicators (such as RMS –Root Mean Square– indicators) based on velocities (in situ mean velocities and PIV velocity fields obtained in laboratory), water depths and discharge (discharge repartition for bifurcation in laboratory or in situ overflow discharge, for CSOs). Results deriving from these tests are transposable to other singularities or special structures encountered in sewer network and suggest that: i) simulated CFD results are sensitive to the roughness coefficient; ii) for the three studied structures, second-order discretisation schemes and SIMPLE or PISO velocity/pressure coupling algorithm are appropriate; iii) scalable wall function associated to the group of k-ε turbulence model for junction flows (with a lateral inflow lower or equal to the main inflow) or for CSOs is appropriate, whereas RSM turbulence model associated to enhanced wall function allows a better representation of bifurcation flows or junctions flows when the lateral inflow is greater than the main inflow. Based on these results and on Jakeman et al. (2006) guideline, a six steps-methodology focused on the using of RANS approach modelling has been proposed. This six steps-methodology is used in order to i) enhance the monitoring of an OTHU (Observatoire de Terrain en Hydrologie Urbaine) CSO located at Ecully accounting for uncertainties on overflow discharge values ;ii) simulate the performance of flowmeters downstream of a junction, defining the best location for these sensors. This methodology is used to design the new overflow discharge measurement device. This device is an international patent. Finally, the application of the methodology led to point out a new flow structure, occurring in the downstream lateral branch of a 90° bifurcation.

Les méthodes osseuses de diagnose sexuelle donnent d’excellents résultats, avec des taux de bonnes prédictions du sexe supérieurs à 90%. Néanmoins, les éléments osseux ne sont pas toujours disponibles et, lorsqu’ils le sont, selon les causes de la mort et/ou les conditions de conservation du corps, ils peuvent avoir été lourdement dégradés. Pour parvenir à déterminer le sexe d’un individu dans ces situations là, des techniques dentaires de diagnose sexuelle sont nécessaires.Contrairement à la détermination du sexe, la problématique de la détermination de l’âge a fait l’objet de nombreuses études et d’autant de propositions de techniques différentes, signe d’une insuffisance de l’ensemble de ces méthodes, que la littérature n’a pas manqué de mettre en exergue.Le premier objectif de ce travail est de proposer une technique dentaire de détermination du sexe simple et précise, de façon à pouvoir répondre à cette problématique anthropologique de façon fiable en l’absence d’autres éléments du squelette.Le second est de développer une méthode d’estimation de l’âge de mise en œuvre facile, rapide et ne nécessitant ni un plateau technique lourd, ni la détérioration du matériel, pour une application aussi bien anthropologique que médico-légale. L’échantillon d’apprentissage était composé de 210 scanners d’individus présentant chacun 4 canines saines. L’échantillon d’application était composé de 55 scanners d’individus, présentant chacun au moins 1 canine saine. Les volumes pulpaire et total de chaque dent ont été modélisés et calculés à l’aide du logiciel Mimics&#63720;. Une régression logistique binaire a été utilisée pour déterminer sept modèles de prédiction du sexe, en fonction du nombre et du type de canines exploitables. La comparaison des aires sous la courbe ROC a mis en évidence une plus grande performance du modèle utilisant les volumes des 4 canines. La méthode des moindres carrés pondérés a été utilisée pour déterminer les équations d’estimation de l’âge pour les sept mêmes modèles. Une plus grande performance du modèle utilisant les volumes des 2 canines maxillaires a été mise en évidence. Toutes ces régressions ont été testées sur l’échantillon d’application pour procéder à leur validation externe.Il s’avère que dans un contexte scientifique historique où les canines mandibulaires constituaient l’outil dentaire de référence dans la diagnose sexuelle, la modélisation tridimensionnelle souligne l’intérêt potentiel de leurs homologues maxillaires. Par ailleurs, il semblerait que l’évaluation de la réduction pulpaire physiologique, imputable à l’apposition naturelle de dentine secondaire sur les canines uniquement, soit un critère performant dans la détermination de l’âge.A défaut d’être parfaits, les modèles proposés à l’issue cette étude permettent donc toutefois de disposer d’éléments fiables, qui ont toute leur place dans l’indispensable faisceau de preuves, nécessaire à la détermination de l’âge et du sexe dans le cadre d’une identification estimative. / The bone methods of sexual diagnosis give excellent results, with rates of good predictions over 90%. Nonetheless, bones are not always available and, when they are, according to the causes of death and/or storage conditions of the body, they may have been severely degraded. To be able to determine the sex of an individual in these situations, dental techniques of sexual diagnosis are needed.In contrast to sex determination, the problem of age estimation has been the topic of many studies and so many different techniques, evidence of failure of all these methods, which the literature did not fail to highlight.The first aim of this work was to propose a simple and accurate dental technique for determining sex, in order to resolve this anthropological issue, reliably in the absence of other elements of the skeleton.The second was to develop a method for estimating the age of easy and fast implémentation, neither requiring a heavy technical support nor damage to material, for an application as well in forensics as in anthropology.The training sample consisted of 210 subjects’ CT-scans, each with four canines healthy. The validation sample was composed of 55 CT-scans, each with at least one healthy canine. Pulp and total volumes of each tooth were modeled and calculated using the software Mimics &#63720;. A binary logistic regression was used to determine seven prediction models of sex, according to the number and type of canine available in practice. The comparison of the areas under the ROC curve showed a greater performance of the model using the volumes of the four canines. The weighted least squares method was used to determine the equations for estimating the age for the same seven models. Greater performance of the model using the volumes of maxillary canines has been demonstrated. All these regressions were tested on the validation sample to perform their external validation.It turns out that in a scientific context where historical mandibular canines were the reference tool in the sex diagnosis, three-dimensional modeling emphasizes the potential value of maxillary counterparts. Moreover, it appears that the evaluation of the physiological pulp reduction, due to the natural apposition of secondary dentine on canines only, is a performance criterion in determining age.The models proposed in this study can therefore have reliable evidence, which all have all their place in the necessary body of evidence for determining age and sex, in estimating identification.

Cette thèse s’inscrit dans le contexte du développement de Bâtiments Basse Consommation. La conception de telles constructions les rend sensibles aux sollicitations internes. Aussi, les outils de thermique du bâtiment existants ne sont pas adaptés pour simuler assez fidèlement ce type de bâtiments, si bien qu’un modèle tridimensionnel et dynamique a été développé ici. Celui-ci présente plusieurs particularités : il s’appuie sur une discrétisation spatiale optimisée des parois, la tache solaire y est localisée et l’intégration des dynamiques des conditions environnementales est assurée par un solveur numérique à pas de temps adaptatif et un seul nœud d’air est considéré. La validation du modèle s’est suivant une confrontation avec des mesures en conditions réelles réalisées dans une cellule de BESTlab d’EDF R&D. Un suivi visuel de la tache solaire a permis de confirmer sa bonne localisation par notre modèle. Des mesures de température en surface complétées par des cartographies thermographiques ont été comparées aux champs de températures simulés, montrant une bonne concordance. Les comparaisons de températures d’air mesurées et simulées ont montré des résidus ne dépassant pas 1,5 ˚C, pour des erreurs moyennes de 0,5 ˚C. La pertinence des deux principales innovations du modèle a été ensuite démontrée : l’utilisation d’entrées échantillonnées à la minute associées à un solveur à pas de temps adaptatif permet de minimiser les erreurs de simulation : en mi-saison, les résidus maximaux sont respectivement de 1 ˚C et 2 ˚C pour des entrées à la minute et à l’heure. En hiver, les températures d’air simulées tendent à plus osciller autour de la consigne quand le pas d’échantillonnage des entrées s’allonge. Deux modèles unidimensionnels, représentatifs de modèles courants, M1D,sol diluant le rayonnement solaire sur le sol seul et M1D,parois le distribuant de façon homogène sur les parois au prorata de la taille de la tache solaire censée les frappées, ne dégradent que légèrement la précision des calculs de température d’air. Cependant, ces modèles 1D ne permettent pas de calcul des champs de températures sur les parois si bien qu’ils présentent des erreurs locales dépassant 20 ˚C aux endroits touchés par la tache solaire. Enfin en hiver, le modèle 3D permet de prédire des consommations de chauffage surestimées de 6,5 % quand M 1D,parois les surestime de 11 % et M1D,sol de 22 %. Les améliorations apportées par notre modèle ont été confirmées pour d’autres types de cellules. D’ailleurs des écarts plus importants entre M1D,sol et le modèle 3D ont été observés pour une cellule dont parois et sol ont des compositions très différentes, alors que l’orientation a aussi un impact. Ce travail confirme la nécessité de représenter plus finement les phénomènes physiques pour des locaux fortement isolés. Des améliorations sont à intégrer, comme la description de l’anisothermie de l’air. / Low energy building constructions become sensitive to internal gains : any internal heating source has an impact on the envelope. Therefore, it is important to evaluate the performance of current transient thermal models when adapted to low energy buildings. This work describes a numerical model to simulate a single room, using a refined spatial three-dimensional description of heat conduction in the envelope but a single air node is considered. The model has been developed for environmental conditions that vary over short time-steps and has integrated the projection of solar radiation through a window onto interior walls : the sun patch. The validation of the model has been done through a detailed comparison between model and measurements. The in situ experiment has been carried out in one of the BESTlab cells (EDF R&D). The sun patch has been followed by a camera to validate its calculated position and surface. Temperature measurements by thermocouples and by thermal cameras have been compared to the models outputs. Differences between air and surface temperatures measured and simulated were never above 1.5 ˚C and mean errors reached 0.5 ˚C. The two innovations of the model have then be proven. Using minute wise weather data and inputs associated to an adaptative solver, enabled to pull down simulation errors : in May maximal differences rised from 1 ˚C to 2 ˚C for respectivelly one minute and hourly wise inputs. More important errors are seen in summer whereas in winter, air temperatures simulated tend to more fluctuate around the set up temperature when the sampling step gets longer. Two one dimensional models, close to traditional taken simulation tools, were used. Model M 1D,sol supposed the incoming radiation to reach only the floor. A 1D model with sun patch movement, called 1D,parois , was also used. These two models evaluated the air temperature with an acceptable error. However, their surface temperatures were still subject to important errors. Thus, for temperature surfaces evaluation, both 1D model presented differences up to 20 ˚C for surfaces touched by the sun patch. In winter, the 3D model can predict heating energy consumptions overestimated by 6.5 % when M 1D,parois overestimated them by 11 % and M1D,sol by 22 %. The improvements brought by our model have been proven also for other cells with different thermal masses. For these cells, differences between M1D,sol and the 3D model could reach 4.5 ˚C. Differences seemed to be more important for low thermal mass cells, and the orientation of the building had a strong impact. This work has confirmed the necessity of representing more accuratelly the descriptions of the enveloppe for strongly insulated rooms. To improve the model, the anisothermal hypotheses of the air should be considered.