La turbine basse-pression est un composant essentiel d'un turboréacteur car elle entraine la soufflante qui génère la plus grande partie de la poussée dans la configuration actuelle des turboréacteurs à double flux. Dans la perspective d'accroître son rendement en termes de consommation de carburant, il y a une recherche permanente dans la réduction du nombre d'aubage (c'est-à-dire la réduction de la masse) qui implique un chargement plus élevé par aube de rotor. Cet environnement est caractérisé par un écoulement dont le nombre de Reynolds est faible ainsi qu'une large diffusion le long de la partie aval de l'extrados. Par conséquent, l'écoulement le long de cette surface est potentiellement sujet à une séparation laminaire qui, suivant le statut de la bulle de recirculation, pourrait causer une diminution de la performance aérodynamique (sillages plus larges et plus profonds). La présente thèse de doctorat se concentre sur l'investigation du phénomène de la transition induite par séparation dans les écoulements de turbines basse-pression. L'accent est mis sur les prédictions numériques basées sur une approche CFD RANS utilisant le modèle innovant de transition γ-Reθt à deux équations de transport (la première équation pour l'intermittence numérique et la seconde équation pour le nombre de Reynolds dont la longueur caractéristique est l'épaisseur de quantité de mouvement au début de transition γ-Reθt). Neuf aubes différentes de rotor de turbine basse-pression constituent une base de données de référence et couvrent les plages de fonctionnement de différents nombres de Reynolds de sortie isentropique, de nombres de Mach de sortie isentropique, d'intensités de turbulence d'entrée, avec ou sans sillage provenant d'une rangée d'aubes amont et avec deux configurations de rugosité locale. Une première analyse de cette base de données met en évidence l'effet de la séparation sur le début de la transition et sur les performances. La définition d'une corrélation a été tentée et permet de lier le taux de diffusion d'un aubage au nombre de Reynolds de sortie isentropique fia la condition de " Bursting ". Une méthodologie numérique fiable et robuste a été établie afin de prédire la transition dans le cas d'un écoulement amont uniforme. Les résultats sont en bon accord avec les mesures expérimentales même si il a été nécessaire d'adapter les conditions limites dans le but de prédire une séparation laminaire numériquement pour des aubages fortement chargés et fia fort taux de diffusion uniquement. La résolution des profils de vitesse de la couche limite permet d'obtenir une évaluation détaillée des paramètres de la topologie de l'écoulement. Cela fournit une information sur l'épaisseur de quantité de mouvement qui est le paramètre principal définissant les corrélations de transition. La technique " Chimère " des maillages recouvrants est utilisée pour faciliter la modélisation des moyens de contrôle passif pour déclencher la transition. C'est une technique appropriée pour l'implémentation de géométries simples ou plus élaborées.

Afin d'étudier la technologie du stockage géologique du CO2, l'intégrité dans le temps de la formation hôte doit être assurée, notamment par l'évaluation de l'altération des propriétés mécaniques des géomatériaux en présence de CO2. Le scénario à l'origine des ces travaux de thèse illustre la possibilité d'une remontée de CO2 le long d'une faille de la roche de couverture. Ce problème géologique est complexe en raison du nombre de paramètres à considérer : température, confinement, effet d'échelle, hétérogénéités du géomatériau, composition du géofluide, réactions chimiques... Ces travaux se cantonnent à l'étude du comportement d'une discontinuité dans un géomatériau sous sollicitations chemo-mécaniques à l'échelle du laboratoire. Ils ont été menés en cotutelle entre le laboratoire SIAME de l'Université de Pau et des Pays de l'Adour (France) et le laboratoire de mécanique des roches et de géologie appliquée de l'Université de Sherbrooke (Québec, Canada). Le programme expérimental a été défini pour caractériser d'une part la refermeture d'une fissure sous effort uni-axial cyclique, et d'autre part le comportement en cisaillement d'un joint rocheux dégradé chimiquement. Les données obtenues ont été utilisés pour modéliser le comportement d'une discontinuité. L'étude du comportement d'une fissure sous effort normal a été conduite par un essai cyclique de traction-compression. Il précise l'évolution des déformations inélastiques générées lors de cycles de fermeture-ouverture d'une discontinuité dans un géomatériau de synthèse (béton). Les données obtenues confirment la présence d'une hystérésis lors des cycles de chargement. L'analyse du champ de déplacements par corrélation d'images lors de la refermeture de fissure permet d'en déduire que les déformations inélastiques seraient gouvernées en partie par les frottements générés lors du ré-emboîtement des lèvres de la discontinuité. La non-correspondance des profils rugueux, entraînant la dissipation d'énergie par frottements illustrée par l'hystérésis, serait directement liée à la présence de contraintes internes. Le comportement d'un joint rocheux altéré sous sollicitation tangentielle utilise l'essai de cisaillement direct. Les surfaces rugueuses de roche sont dégradées chimiquement par immersion dans des solutions acides. La numérisation de ces surfaces par un profilomètre laser indique de légères modifications de la rugosité par dissolution. Les essais révèlent une profonde modification du comportement des joints attaqués : le pic de contrainte disparaît, la rigidité du joint diminue avec l'angle de dilatance et la phase de contractance est accrue. Ces évolutions s'expliquent par l'endommagement chimique qui accentue la non-correspondance des profils rugueux des joints et diminue les propriétés mécaniques du matériau de part et d'autre de la discontinuité. D'un point de vue numérique, l'apport de ces travaux de thèse réside dans la modélisation du comportement d'une discontinuité par le couplage d'un modèle élasto-plastique endommageable continu avec une résolution discrète du problème de contact/frottement. Les modélisations sont réalisées avec le code de calcul aux éléments finis Cast3M. Les géométries des lèvres de discontinuités sont directement modélisées à partir des profils de rugosité issus des numérisations expérimentales. Les résultats numériques représentent correctement les phénomènes de frottements constatés expérimentalement. Enfin, une modélisation de l'essai de cisaillement des joints altérés est réalisée en couplant le modèle mécanique avec un modèle d'endommagement chimique.

David Katz a fait l’observation que le mouvement entre la peau et l’objet est aussi important pour le sens du toucher que la lumière l’est pour la vision. Un stimulus tactile déplacé sur la peau active toutes les afférences cutanées. Les signaux résultants sont très complexes, covariant avec différents facteurs dont la vitesse, mais aussi la texture, la forme et la force. Cette thèse explore la capacité des humains à estimer la vitesse et la rugosité de surfaces en mouvements. Les bases neuronales de la vitesse tactile sont aussi étudiées en effectuant des enregistrements unitaires dans le cortex somatosensoriel primaire (S1) du singe éveillé. Dans la première expérience, nous avons montré que les sujets peuvent estimer la vitesse tactile (gamme de vitesses, 30 à 105 mm/s) de surfaces déplacées sous le doigt, et ceci sans indice de durée. Mais la structure des surfaces était essentielle (difficulté à estimer la vitesse d’une surface lisse). Les caractéristiques physiques des surfaces avaient une influence sur l’intensité subjective de la vitesse. La surface plus rugueuse (8 mm d’espacement entre les points en relief) semblait se déplacer 15% plus lentement que les surfaces moins rugueuses (de 2 et 3 mm d’espacement), pour les surfaces périodiques et non périodiques (rangées de points vs disposition aléatoire). L’effet de la texture sur la vitesse peut être réduit en un continuum monotonique quand les estimés sont normalisés avec l’espacement et présentés en fonction de la fréquence temporelle (vitesse/espacement). L'absence de changement des estimés de vitesse entre les surfaces périodiques et non périodiques suggère que les estimés de rugosité devraient aussi être indépendants de la disposition des points. Dans la deuxième expérience, et tel que prévu, une équivalence perceptuelle entre les deux séries de surfaces est obtenue quand les estimés de la rugosité sont exprimés en fonction de l'espacement moyen entre les points en relief, dans le sens de l'exploration. La troisième expérience consistait à rechercher des neurones du S1 qui pourraient expliquer l’intensité subjective de la vitesse tactile. L’hypothèse est que les neurones impliqués devraient être sensibles à la vitesse tactile (40 à 105 mm/s) et à l’espacement des points (2 à 8 mm) mais être indépendants de leur disposition (périodique vs non périodique). De plus, il est attendu que la fonction neurométrique (fréquence de décharge/espacement en fonction de la fréquence temporelle) montre une augmentation monotonique. Une grande proportion des cellules était sensible à la vitesse (76/119), et 82% d’entres elles étaient aussi sensibles à la texture. La sensibilité à la vitesse a été observée dans les trois aires du S1 (3b, 1 et 2). La grande majorité de cellules sensibles à la vitesse, 94%, avait une relation monotonique entre leur décharge et la fréquence temporelle, tel qu’attendu, et ce surtout dans les aires 1 et 2. Ces neurones pourraient donc expliquer la capacité des sujets à estimer la vitesse tactile de surfaces texturées. / David Katz showed that movement between the skin and an object is as important for touch as light is to vision. Moving tactile stimuli activate all of the cutaneous afferents involved in discriminative touch. The resultant signals are complex, varying with multiple factors including speed and also texture, local shape, and force. This thesis explored the human ability to estimate the speed and roughness of moving tactile stimuli. The neuronal basis underlying tactile speed perception was investigated using single unit recordings from primary somatosensory cortex (S1) in awake monkeys. In the first psychophysical experiment, we showed that subjects (n=26) can scale tactile speed (range, 30-105 mm/s), and this, contrary to previous studies, in a situation in which the duration of each trial was constant across all speeds tested. Surface structure was, in contrast, essential since subjects had difficulty scaling the speed of a smooth surface. Moreover, the physical characteristics of the surfaces influenced tactile speed perception. The roughest surface (8 mm raised-dot spacing) seemed to move 15% slower than the smoother surfaces (2 and 3 mm spacing), and this independently of dot disposition (periodic: rectangular array of raised dots vs non periodic: random dots). The effects of surface texture on speed were reduced to a single continuum when the estimates were normalized by dot spacing and plotted as a function of temporal frequency (speed/dot spacing). The absence of any difference in speed scaling as a function of dot disposition (periodic vs non periodic) suggested that tactile roughness should also be independent of dot disposition. A second psychophysical experiment (n=15) confirmed our hypothesis, showing perceptual equivalence for the periodic and non periodic surfaces when these were matched for dot spacing in the direction of the scan. The third experiment investigated the neuronal mechanisms that underlie subjective tactile speed perception, by recording the responses of cutaneous neurones in the hand representation of S1 cortex to the displacement of textured surfaces under the finger tips of two awake rhesus monkeys. The hypothesis was that neurones implicated in tactile speed perception should be sensitive to tactile speed (similar range to that used in the human experiments) and dot spacing, but be independent of dot disposition (periodic vs non periodic). Furthermore, we predicted that the neurometric function (discharge frequency/dot spacing as a function of temporal frequency) would show a monotonic relation. A large proportion of S1 neurones were sensitive to speed (76/119); 82% of these were also sensitive to texture. Speed sensitivity was widely distributed across the three areas that comprise the cutaneous hand representation, areas 3b, 1, and 2. Of 94 neurons fully tested (periodic and nonperiodic surfaces), the large majority of speed-sensitive cells (60/64) showed a significant monotonic relation with temporal frequency for both surfaces when discharge frequency was normalized by dot spacing. The neurones with the strongest relation to temporal frequency were concentrated in caudal S1, areas 1 and 2, and likely contribute to the human ability to scale tactile speed.

Depuis de nombreuses années, le givrage a été identifié comme un danger dans le domaine de l’aéronautique.L’accrétion de givre se produit lorsque des gouttelettes d’eau surfondue se déposent sur une surface, enparticulier le bord d’attaque d’une aile ou la lèvre d’entrée d’air moteur, et gèlent après l’impact. Ceci peutensuite engendrer une dégradation des performances aérodynamiques, un dysfonctionnement des sondes ouencore un endommagement du moteur. C’est pourquoi cette problématique est étudiée avec attention. Lesessais en vol et en soufflerie étant longs et coûteux, la simulation numérique de l’accrétion de givre est devenueun outil nécessaire dans le processus de conception et de certification des avions. Cette thèse s’inscrit dans le contexte de la modélisation 3D de l’accrétion de givre, et plus particulièrement des couches limites dynamique et thermique qui se développent autour du corps givré. Les outils numériques devant être rapides et robustes, l’approche proposée dans cette thèse pour le calcul aérodynamique est une méthode couplée Euler/couche limite intégrale. Ainsi, un modèle intégral est développé pour représenterle développement de la couche limite dynamique. La partie thermique est modélisée soit par une méthodesimplifiée basée sur des approches algébriques, soit par une méthode intégrale. Cette modélisation des coucheslimites dynamique et thermique est valable sur paroi lisse ou rugueuse et permet de fournir notamment lecoefficient de frottement et le coefficient d’échange thermique qui sont nécessaires pour un calcul d’accrétion degivre. Les équations intégrales de couche limite, associées à leurs relations de fermeture, sont ensuite résoluespar une méthode Volumes-Finis sur maillage surfacique non structuré, qui est bien adaptée pour les géométriescomplexes. De plus, des traitements numériques spécifiques sont mis en œuvre pour améliorer la précision dela méthode au voisinage du point d’arrêt et pour rendre le code robuste au passage du décollement.Après la validation de la méthode de couche limite, le code est utilisé dans les chaînes de givrage 2D et 3Dde l’ONERA pour des applications d’accrétion de givre. Ceci permet de montrer l’intérêt de la méthode entermes de robustesse et de précision par rapport aux codes de couche limite habituellement utilisés dans lescodes de givrage actuels. / Icing has since long been identified as a serious issue in the aeronautical world. Ice accretion occurs whensupercooled water droplets impinge on a surface, particularly the leading edge of a wing or an engine inlet, andfreeze after the impingement. This can lead to degradation of aerodynamic performances, sensor malfunctionor engine damage. This is why this issue is being carefully studied. The lengthy and costly flight and windtunnel tests have made numerical simulation of ice accretion a necessary tool in the aircraft design andcertification process. The present work deals with the 3D numerical modeling of ice accretion, and more particularly the modeling of the dynamic and thermal boundary layers that develop around an iced body. Since numerical tools must befast and robust, the approach proposed in this thesis for aerodynamic computation is a coupled Euler/integralboundary layer method. Thus, an integral model is developed to represent the development of the dynamicboundary layer. The thermal part is modeled either by a simplified method based on algebraic approaches,or by an integral method. This modeling of the dynamic and thermal boundary layers is valid on smoothor rough walls and provides the friction coefficient and heat exchange coefficient that are necessary for thecalculation of ice accretion. The integral boundary layer equations, associated with their closure relations,are then solved by a Finite-Volume method on unstructured surface mesh, that is well suited for complexgeometries. In addition, specific numerical treatments are implemented to improve the accuracy of the methodin the vicinity of the stagnation point and to make the code robust to separated boundary layers.After validation of the boundary layer method, the code is used in ONERA’s 2D and 3D icing tools foricing applications. This demonstrates the value of the method in terms of robustness and accuracy comparedto the boundary layer codes more commonly used in current icing tools.

Pour améliorer l'efficacité des procédés thermiques de craquages et réduire les phénomènes de cokage liés à la température de paroi trop élevée, l'utilisation de tubes nervurés est une technique potentiellement car elle permet d'améliorer le mélange et d'augmenter les transferts de chaleur. Cependant, la perte de charge est significativement augmentée. En raison de la complexité de l'écoulement turbulent, du système chimique et du couplage turbulencechimie, il est difficile d'estimer a priori la perte réelle en termes de sélectivité des tubes nervurés. Les expériences représentatives de laboratoire combinant turbulence, transferts de chaleur et chimie sont très rares et trop coûteuses à l'échelle industrielle. Dans ce travail, l'approche simulation aux grandes échelles résolue à la paroi (WRLES) est utilisée pour étudier écoulement non-réactif puis réactif dans des tubes à la fois lisses et nervurés, pour quantifier leur impact sur la turbulence et sur la chimie. Le code AVBP, qui résout les équations de Navier-Stokes compressibles pour les écoulements turbulents, est utilisé avec des schémas chimique réduites du craquage de l'éthane puis du butane. L'écoulement à la paroi est analysé en détail et comparé pour les deux géométries, fournissant des informations utiles pour le développement ultérieur de modèles de parois pour ce type de rugosité. L'impact de la résolution du maillage et du schéma numérique est également discuté, pour trouver le meilleur compromis entre coût et précision de calcul pour une application industrielle. L'impact des structures d'écoulement turbulent ainsi que leurs effets sur le transfert thermique et le mélange sur les réactions chimique sont étudiés à la fois pour les tubes lisses et les tubes nervurés. Perte de pression, transfert de chaleur et conversion chimique sont finalement comparés. / To improve the efficiency of thermal-cracking processes, and to reduce the coking phenomena due to high wall temperature, the use of ribbed tubes is an interesting technique as it allows better mixing and heat transfer. However it also induces significant increase in pressure loss. The complexity of the turbulent flow, the chemical system, and the chemistry-turbulence interaction makes it difficult to estimate a priori the real loss of ribbed tubes in terms of selectivity. Experiments combining turbulence, heat transfer and chemistry are very rare in laboratories and too costly at the industrial scale. In this work, Wall-Resolved Large Eddy Simulation (WRLES) is used to study non-reacting and reacting flows in both smooth and ribbed tubes, to show the impact of the ribs on turbulence and chemistry. Simulations were performed with the code AVBP, which solves the compressible Navier-Stokes equations for turbulent flows, using reduced chemistry scheme of ethane and butane cracking for reacting cases. Special effort was devoted to the wall flow, which is analyzed in detail and compared for both geometries, providing useful information for further development of roughness-type wall models. The impact of grid resolution and numerical scheme is also discussed, to find the best trade-off between computational cost and accuracy for industrial application. Results investigate and analyze the turbulent flow structures, as well as the effect of heat transfer efficiency and mixing on the chemical process in both smooth and ribbed tubes. Pressure loss, heat transfer and chemical conversion are finally compared.

L'étude concerne les couches limites turbulentes dans le cas d'écoulements à surface libre sur fond rugueux homogène. Afin de réaliser cette étude, deux dispositifs de mesure PIV par stéréoscopie (PIV 2D-3C) ont été mis en place avec comme double objectif de fournir les lignes directrices au design d'un système stéréoscopique PIV in situ et d'étudier l'influence de la faible submersion des éléments rugueux sur la structure universelle de la couche limite sur fond rugueux, i.e. pour des submersions h/D comprises entre 0,33 et 0,66 (avec h la hauteur des rugosités et D la hauteur d'eau). Pour le premier objectif, le dispositif de mesure a permis un accès optique facilité par un point de vue incliné des caméras à la zone proche des rugosités dans des conditions naturelles d'écoulement (turbidité et éclairement naturels). Les mesures de vitesse ont été faites dans un canal hydraulique de petite dimension (12 m x 0,5 m x 0,25 m) rempli d'hémisphères positionnées en quinconce. Une étude paramétrique de l'influence de l'inclinaison des caméras ainsi que de la turbidité de l'eau sur la qualité des mesures de vitesse a été entreprise suivie par une nouvelle méthodologie basée sur l'analyse de l'intensité lumineuse dans le système. Il a été montré que l'écoulement est correctement résolu jusqu'à une turbidité d'environ 25 NTU avec un angle d'inclinaison par rapport au plan vertical de mesure de 25°. Pour le second objectif, les investigations expérimentales ont été réalisées dans une veine hydraulique de plus grande dimension (26 m x 1,10 m x 0,50 m), dont le fond rugueux est constitué par des cubes en PVC de 2 cm de côté comme dans l'étude de Florens et al. (2013). Les résultats mettent en évidence que l'étendue de la sous-couche rugueuse augmente avec la submersion pour finalement occuper toute la colonne d'eau dans le cas de la plus faible submersion (h/D=0,66). Malgré cela, une loi logarithmique est tout de même observée, et ce, quelle que soit la submersion étudiée. / This work deals with turbulent boundary layers in open-channel flows over rough homogeneous beds. The objectives of this work are, first, to provide some guidance for the design of an efficient in situ stereoscopic PIV measurements system (SPIV), and, second, to assess the effect of the relative submergence on the universal turbulent boundary layer structure for very high relative submergence ratios, i:e: 0.33 < h/D <0.66 (where h is the roughness height and D the water depth). For the first objective, a stereoscopic PIV configuration was set-up with steeply inclined camera viewpoints in order to improve the image quality and the optical access into the bed canopy under naturally occurring turbid conditions. Velocity measurements were undertaken in a 12 m x 0.5 m x 0.25 m open-channel flume filled with staggered hemispheres as surrogates for bed river peebles. A parametric study was then carried out to both analyze how the turbidity and camera angle impact the quality of PIV measurements. An innovative light intensity-based methodology was developed and applied to perform data analysis. The latter shows good PIV results up to 25 NTU with an optimal camera angle with respect to the vertical PIV measurements plane of 25. The SPIV measurements for the second objective were performed in a 26-m-long, 1.10-m-wide and 0.50-m-deep steep open channel filled with 2-cm cubes as in Florens et al. (2013). The results show that the extent of the roughness sublayer increases with the relative submergence to fill the entire water column for the highest relative submergence investigated. Despite this, the logarithmic law is still observed even for the highest relative submergence studied (h/D=0.66).

Le collage par adhérence moléculaire est un collage basé sur la mise en contact de deux surfaces sans l'utilisation de colle ou matériaux additionnels. Ce procédé de collage est utilisé dans de nombreuses applications, notamment dans les domaines de l'optique terrestre et spatiale. Bien qu'un prototype ait déjà passé avec succès l'environnement spatial - où les contraintes d'utilisations sont différentes de celles rencontrées sur Terre - la spatialisation de cette technologie nécessite une caractérisation plus fine du procédé ainsi qu'une amélioration de la tenue mécanique des interfaces adhérées afin de valider les normes de l'Agence spatiale Européenne. Pour répondre à cette problématique de spatialisation de la technologie, des essais mécaniques ainsi que des analyses chimiques ont été réalisés dans le but d'étudier l'influence de certains paramètres du procédé ) sur la tenue mécanique et l'énergie de collage. Ces essais ont également été réalisés afin de comparer les deux matériaux étudiés : le verre de silice et le Zérodur vis-à-vis de l'adhésion. A l'issue de ces essais, les paramètres à appliquer permettant de doubler la tenue mécanique des interfaces adhérées ont été déterminés. Parallèlement, une loi phénoménologique reliant l'énergie de collage aux précédents paramètres du procédé a été développée ainsi qu'un modèle macroscopique visant à décrire l'intensité d'adhésion. Ces deux modèles une fois couplés permettent de modéliser le comportement normal de l'interface en fonction des paramètres du procédé. Enfin, ces deux lois sont implémentées dans un code éléments finis afin de simuler la propagation de la fissure lors de l'essai de clivage au coin. / Direct bonding consists in joining two surfaces without the use of any adhesive or additional material. This process is used in several applications, particularly in terrestrial and spatial optics. Although a prototype passed with success spatial environment - where constraints involved are very different from those encountered on Earth - this technology requires a more detailed characterization and an improvement of the mechanical strength of bonded interfaces in order to validate the European Space Agency standards. To address this issue, mechanical tests (double shear tests, cleavage tests and wedge tests) and chemical analysis (wetting tests and XPS spectroscopy) were performed in order to study the influence of some process parameters (roughness, relative air humidity during room temperature bonding, the annealing temperature and time) on the mechanical strength and the bonding energy. These tests compared the two materials used: fused silica glass and Zerodur glass. As a result of these tests, optimal parameters doubling the mechanical strength were also obtained. In the same time, a phenomenological law relating the bonding energy to the previous parameters is developed as well as a macroscopic model to describe the adhesion intensity. Both models when coupled describe the normal behavior of the bonded interface depending on the process parameters. Then, the both laws are implemented in a finite elements model in order to simulate the crack propagation during the wedge test.

Le comportement mécanique en cisaillement sous contrainte normale constante de joints rocheux est étudié en utilisant une approche numérique par éléments discrets (DEM Discrete Element Model). Les influences respectives de la rugosité des surfaces des joints, de l'élasticité des épontes, de la rupture des aspérités de surface et du niveau de contrainte de compression sur les comportements en fermeture et cisaillement des joints rocheux sont particulièrement analysées. Pour la première fois la rugosité des joints considérée comme auto-affine est utilisée avec DEM pour étudier le frottement des joints rocheux. Cette rugosité est décrite par l’intermédiaire de trois paramètres :exposant de rugosité auto-affine, longueur de corrélation auto-affine et variance des fluctuations de hauteur. Sur la base d’un algorithme fondé sur la méthode spectrale, huit surfaces auto-affines isotropes correspondant à différentes rugosités ont été générées. Ces surfaces numériques sont utilisées comme moules permettant de générer les surfaces composées d’éléments discrets utilisées dans la suite de l’étude. La modélisation par éléments discrets s’appuie sur une calibration des propriétés élastiques effectuée à partir d’un volume élémentaire représentatif suivie de l’implémentation d’un critère elliptique de contraintes de rupture (au niveau des lois d’union entre éléments) permettant de simuler les grandes lignes du comportement quasi-fragile d’un mortier(utilisé lors d’expérimentations antérieures). Sur cette base et une fois les surfaces rugueuses implémentées dans les modèles DEM, les essais de fermeture (test de compression) des huit joints sont effectués sous deux niveaux de contrainte de compression : 14 MPa et 21 MPa. Par la suite, les joints sont cisaillés selon deux directions perpendiculaires. Pour chaque direction de cisaillement et chaque niveau de contrainte de compression, les joints sont testés en utilisant trois modèles mécaniques différents : 1) modèle rigide dans lequel, à l’exception des surfaces de joint en contact,les épontes ne peuvent pas se déformer, 2) modèle élastique dans lequel les épontes peuvent se déformer dans leur volume et 3) modèle élastique-fracture dans lequel les épontes peuvent se déformer dans leur volume et les liens entre les particules peuvent rompre selon le critère elliptique de contrainte. L'utilisation de ces trois modèles mécaniques différents permet d'étudier de façon systématique l'influence de la rugosité seule (modèle rigide), l'influence de l'élasticité et de la rugosité (modèle élastique) et enfin, l'effet combiné de la rugosité, de l'élasticité et de la rupture(modèle élastique-fracture). L’étude des résultats obtenus lors des simulations DEM est accompagnée d’une analyse énergétique permettant d’estimer l’évolution de l’énergie élastique stockée dans le système, de l’énergie de friction, du travail associé à la dilatance du joint et de l’énergie dissipée au cours de l’essai de cisaillement. / The shear behavior of rock joints under constant normal stress is studied using Discrete Element Method (DEM). The respective influences of joint surface roughness, elasticity of medium, fracture of surface asperities, and level of compression load on the closure and shear behaviors of rock joints are particularly analyzed. For the first time the roughness of the joints considered as self-affine is use dwith DEM to study the friction of rock joints, the roughness is described through three parameters:self-affine roughness exponent, self-affine correlation length and height variance. Using a numerical algorithm based on spectral method, eight isotropic self-affine surfaces corresponding to different roughness are generated. Latter, numerical surfaces are used as molds to generate the discrete elements surfaces. The discrete element modeling is premised on a preliminary calibration of the elastic properties performed on a representative elementary volume and on the implementation of the fracture properties (elliptic fracture criterion expressed in stress) describing with a reasonable accuracy the quasi-brittle fracture behavior of mortar (used in previous experimental tests). On this basis and once the roughness surfaces implemented in DEM, the simulations of the compression/closure test are performed on the eight joints and this for two compression stress levels: 14 MPa and 21 MPa. Then, the eight DEM joints are sheared along two perpendicular directions. For each shear direction and each level of compression stress, the joints are tested through three different mechanical models: 1) rigid model in which the medium cannot deform excepted at the contact surface of joints, 2) elastic model in which the medium can deform in its volume and 3) elastic-fracture model in which the medium can deform in its volume and the bondsbetween discrete elements can failed according to the elliptic fracture criterion. The use of these three mechanical models allows studying systematically the influence of the roughness alone (rigidmodel), the influence of elasticity and roughness (elastic model) and finally, the combined effect ofthe joint roughness, of the elasticity and of the fracture (elastic-fracture model). The study of the results obtained from the DEM simulations is followed by an energetic analysis allowing theestimation of the evolutions, as a function of the shear displacement, of the elastic energy stored inthe system, of the friction energy, of the work related to the joint dilatancy and of the energy dissipated by internal damping of the DEM. / En esta tesis se estudia la fricción en juntas rocosas utilizando el Método de Elementos Discretos (DEM). En particular, se estudia la influencia de la rugosidad de las superficies de la junta, la elasticidad, la fractura, y el nivel de carga de compresión sobre el comportamiento de cierre y de cizalla de las juntas rocosas. Por primera vez la rugosidad de las juntas considerada como auto-afín esutilizada para estudiar la fricción de juntas rocosas, la rugosidad se describe mediante tres parámetros: el exponente de rugosidad, la longitud de correlación auto-afín y la varianza de alturas. Mediante un algoritmo de computadora basado en métodos espectrales, ocho superficies autoafines isotrópicas con diferente rugosidad fueron creadas. Posteriormente, las ocho superficies fueron utilizadas como moldes para generar las juntas utilizando elementos discretos. Antes de realizar las simulaciones de compresión y cizallaura, se calibraron las propiedades elásticas y defractura (criterio de fractura elíptico basado en esfuerzos) de las juntas numéricas a los datos experimentales (obtenidos previamente) de unas muestras de mortero mediante la utilización de un volumen elemental representativo (REV). Una vez que las propiedades mecánicas de las juntas se obtuvieron mediante la calibración del REV, se realizaron las pruebas de cierre (prueba de compresión) de las ocho juntas DEM. Se utilizaron dos niveles de esfuerzo de compresión para laspruebas de cierre: 14 MPa y 21 MPa. Después, las ocho juntas DEM fueron cizalladas en dos direcciones mutuamente perpendiculares. Para cada dirección de cizalla y cada nivel de esfuerzo decompresión (14 y 21 MPa), las juntas fueron cizalladas usando uno de los tres modelos mecánicos siguientes: 1) un modelo rígido, en el que las juntas no se pueden deformar, excepto en su superficie,2) un modelo puramente elástico, en el que las juntas se pueden deformar en todo su volumen y 3)un modelo elástico con fractura en el que las juntas se pueden deformar en su volumen y, si elesfuerzo sobre las uniones entre partículas excede cierto nivel de esfuerzo máximo, las uniones se rompen de una manera irreversible. El uso de estos tres modelos mecánicos nos permitirá estudiar de manera sistemática: la influencia de la rugosidad (modelo rígido), la influencia de la elasticidad y rugosidad (modelo puramente elástico) y, finalmente, el efecto combinado de la rugosidad de las juntas, la elasticidad y la fractura (modelo elástico con fractura). El estudio de los resultados obtenidos de las simulaciones DEM es seguido por una análisis energético el cual permite estudiar la evolución de los diferentes tipos de energía en función del desplazamiento de cizalla: energía elástica almacenada en el sistema, energía de fricción entre elementos discretos, el trabajo relacionado conla dilatación de la junta y la energía disipada por el amortiguamiento interno del DEM.

Au cours de cette thèse ont été développées de nouvelles techniques de mesure des propriétés mécaniques par nanoindentation, adaptées à des conditions expérimentales difficiles : les très faibles enfoncements, les surfaces rugueuses, et la température. Une simulation numérique par éléments finis d’un échantillon de silice chauffé à 80°C, indenté par un diamant Berkovich dont la température initiale est de 25°C, dans de l’air à 60 °C, a mis en évidence la nécessité de chauffer l’indenteur. De plus, les essais expérimentaux à haute température effectués sur ce même échantillon ont permis de montrer que le signal de déplacement de l’indenteur est fortement perturbé par de faibles variations de température (<0,1°C), rendant le calcul des propriétés mécaniques imprécis avec la technique de mesure classiquement utilisée en nano-indentation. Une nouvelle technique, reposant sur la mesure de l'amplitude de la seconde harmonique du signal de déplacement, a été développée pendant cette thèse. Elle permet la détermination des propriétés mécaniques indépendamment de la mesure de l’enfoncement. Elle est donc adaptée pour des tests à haute température. Elle été expérimentée sur des matériaux homogènes (silice, PMMA), sur un matériau (monocristal de calcite) possédant une dureté plus élevée en surface (Indentation Size Effect), et sur des couches minces de PMMA déposées sur un substrat de silicium, à température ambiante. Les résultats ont montré que les propriétés mécaniques sont mesurées de façon plus précise aux faibles enfoncements. Les essais expérimentaux sur la calcite ont mis en évidence que l'Indentation Size Effect observé sur cet échantillon est mieux détecté avec la technique de la seconde harmonique. Elle permet aussi de calculer la pénétration de l'indenteur a posteriori. Une seconde technique de mesure des propriétés mécaniques, reposant sur le calcul de la dérivée de la hauteur de contact en fonction de la pénétration, a été développée. Elle permet une mesure des propriétés mécaniques des matériaux par nano-indentation plus précise que la méthode classique à température ambiante aux faibles pénétrations. Les mesures sont aussi améliorées sur les échantillons rugueux et pour les essais réalisés à haute température. / The aim of this thesis was the development of new measurement techniques based on nano-indentation, adapted for difficult experimental conditions: small penetration depths, rough surfaces, and high temperature. The thermal contact between a Berkovich indenter initially at 25 °C and a fused silica sample heated at 80 °C, in air at 60 °C, was numerically simulated by finite elements method. The results showed the necessity to heat the indenter in order to avoid effects due to the difference of temperature between the two solids. Furthermore, high temperature nanoindentation tests showed the displacement signal is greatly influenced by temperature variations (<0.1°C), resulting in imprecise mechanical properties calculation. A new experimental technique, based on the measurement of the amplitude of the second harmonic of displacement, was developed. With this method, the determination of the mechanical properties is independent of the indentation depth measurement. So, the second harmonic method is adapted to high temperature tests. It was tested on homogeneous materials (fused silica and PMMA), on a sample which is known to exhibit an Indentation Size Effect (calcite), and on thin PMMA layers deposited onto silicon wafers, at room temperature. With the second harmonic method, the mechanical properties are measured more precisely at small penetration depths. Experiments performed on the calcite sample showed that the Indentation Size Effect is more precisely measured with this new method. Furthermore, the indentation depth can be calculated “a posteriori” with second harmonic method. A second new measurement technique, based on the derivative of the contact depth with respect to the indentation depth, was developed. With this simple method, the mechanical properties are more precisely measured at room temperature at small indentation depths. Measurements are also improved on rough samples and at high temperature.

Dans le contexte de l’amélioration des modèles de prévision en acoustique extérieure, ces travaux de thèse se focalisent sur la modélisation des effets des irrégularités de la surface du sol sur la propagation acoustique. Pour ce faire, des méthodes numériques temporelles sont utilisées : d’une part, la méthode FDTD basée sur la résolution des équations de l’acoustique par des schémas aux différences finies, et d’autre part, la méthode des lignes de transmission (TLM). La modélisation des effets de la rugosité de surface est abordée en considérant le formalisme de l’impédance effective. Deux modèles d’impédance effective sont étudiés : le premier caractérise les effets d’une rugosité déterministe constituée de diffuseurs de géométrie constante, le second caractérise les effets moyens d’une rugosité aléatoire définie par un spectre de rugosité. Ce second modèle est validé expérimentalement par une campagne de mesures en salle semi-anéchoïque, audessus de surfaces rugueuses dont la rugosité a été définie très précisément. Les deux modèles d’impédance effective sont également validés par des simulations numériques FDTD et TLM. La possibilité d’implémenter ces conditions d’impédance effective dans les deux codes temporels est ainsi montrée, ce qui permet de modéliser les effets de la rugosité sans avoir à réaliser un maillage précis du profil des hauteurs de la surface du sol. Une campagne de mesures de l’impédance de différents terrains est réalisée afin d’étudier les effets de la variabilité spatiale et saisonnière de l’impédance sur la prévision des niveaux sonores. Les impédances mesurées lors de cette campagne sont également utilisées comme des données d’entrée réalistes pour le code TLM, afin de simuler et d’étudier les effets de la propagation acoustique au-dessus d’un sol hétérogène présentant une impédance spatialement variable. / In the context of prediction models improvement for outdoor sound propagation, this work focuses on the modelling of the effects of ground irregularities on sound propagation. Time-domain numerical methods are used: on one hand, the solving of the governing equations by finite difference schemes (FDTD method), and on the other hand, the transmission line matrix (TLM) method. Effective impedance is considered to model the effects of surface roughness. Two effective impedance models are studied: the first one takes into account the effects of a deterministic roughness formed by scatterers of constant geometry, the second one takes into account the mean effects of a random roughness defined by a roughness spectrum. This second model is validated experimentally by a measurements campaign carried out in a semi-anechoic chamber, above rough surfaces whose roughness profiles were precisely designed. The two effective impedance models are also validated by FDTD and TLM simulations. The possibility to use the effective impedances directly into the numerical methods is then shown, allowing the modelling of roughness effects without meshing the exact height profile of the ground surface. A measurements campaign of the impedance of different grounds is performed in order to assess the effects of space and seasonal variability of the impedance on the sound levels predictions. The results of this campaign are also used as realistic entry data for the TLM code, and the propagation above a heterogeneous ground showing spatially variable impedance is simulated.