Le principal objectif de cette thèse est de caractériser l'écoulement dans un anévrisme abdominal aortique (AAA) sous différentes conditions physiologiques et à différents stades de son développement. Cette étude est consacrée aux AAA ax- isymétriques, modélisés comme une dilatation de profil gaussien et de section circu- laire. Ainsi, les résultats s'appliquent surtout aux étapes précoces du développement d'un AAA. Le modèle d'AAA est caractérisé par une hauteur maximale H et une largeur W , l'unité de mesure étant le diamètre d'entrée de l'artère. Pour com- mencer, la dynamique est étudiée pour les écoulements stationnaires. La stabilité globale de ces écoulements de base est analysée en calculant les valeurs propres et les fonctions propres pour des perturbations de faible amplitude. Pour comprendre les mécanismes d'instabilité, le transfer d'énergie entre l'écoulement de base et les perturbations est calculé. L'écoulement pour des AAA peu profonds (ou de grande longueur) se déstabilise par un mécanisme de 'lift-up' et les perturbations ampli- fiées sont stationnaires. Des anévrismes plus localisés (ou plus profonds) deviennent instables pour des nombres de Reynolds plus élevés, sans doute par instabilité el- liptique ; dans cette situation, les perturbations sont des modes oscillants. Dans le cas des écoulements pulsés, deux types de profil de débit physiologique ont été considérés dans cette étude, correspondant à une situation de repos ou d'exercice physique. Ces écoulements restent collés aux parois pendant la phase de systole et un écoulement décollé est généralement observé pendant la décélération après le maximum de systole. Dans cette phase, un vortex se forme à l'extrémité aval. Ce vortex s'agrandit au cours du temps et impacte l'extrémité aval de l'AAA, ce qui conduit à de forts gradients de contrainte pariétale, qui ne sont pas observés dans les cas sains. Il a été observé que les conditions d'écoulement varient significative- ment avec les nombre de Womersley (Wo) et de Reynolds (Re); l'écoulement reste attaché aux parois plus longtemps pour des nombres de Womersley croissants. Le principal effet d'une augmentation de Re est un renforcement du vortex primaire qui se forme après le maximum de systole. Les décollements de l'écoulement, l'impact de vortex au bord aval de l'AAA ou encore de faibles contraintes pariétales oscil- lantes (des caractéristiques importantes dans les cas d'anévrismes pathologiques) sont observés même pour des anévrismes de faible profondeur. Pour des anévrismes plus développés, des vortex multiples sont observés tout au long du cycle dans la cavité de l'AAA. Une analyse de stabilité de ces écoulements de base pulsés a mon- tré que le maximum des perturbations se développe vers l'extérmité aval des AAA. Cependant, les perturbations ne sont pas complètement confinées dans la cavité de l'AAA et se développent aussi au-delà en aval. On en déduit qu'une fois qu'un AAA s'est développé, les perturbations affectent aussi les artères saines en aval de l'AAA. Enfin, en considérant deux profils équivalents d'AAA, de formes sinusoï- dale et gaussienne, la sensibilité des résultats aux détails de la géométrie a pu être établie.

Les oiseaux et poissons se déplacent dans leur environnement fluide en interagissant avec l'air/eau qui les entoure. Pour des régimes inertiels, les mécanismes de propulsion se basent sur un transfert de quantité de mouvement au fluide; les battements d'ailes ou de nageoires générant un jet dans le sillage de l'animal qui le propulse vers l'avant. Pour les oiseaux comme pour les poissons, les structures utilisées possèdent une certaine flexibilité, et sont donc susceptibles de plier de façon importante. La littérature montre que ces déformations passives peuvent améliorer les performances de propulsion lorsqu'elles sont exploitées de façon constructive. Le détail des mécanismes en jeu reste cependant mal compris. L'objectif de cette thèse est d'étudier, à travers deux modèles biomimétiques, la façon dont une structure battante déformable génère des forces de propulsion. Le premier modèle est une version mécanique simplifiée d'insecte dotée d'ailes flexibles, tandis que le deuxième est un nageur dont le corps élastique reproduit le mouvement d'ondulation d'une anguille. Nous montrons que la façon dont ces systèmes se déforment passivement est déterminante pour leurs performances, et que leur réponse élastique peut être décrite par des modèles théoriques simplifiés d'oscillateurs forcés. Ces modélisations mettent par ailleurs en avant le rôle crucial joué par le frottement fluide quadratique qui s'oppose aux mouvements de battements de la structure. Ce résultat introduit l'idée, un peu contre-intuitive, qu'il peut s'avérer avantageux de dissiper une part de son énergie dans le fluide pour améliorer ses performances.

L'objectif de cette thèse est de développer un modèle numérique instationnaire permettant de simuler l'évolution du climat distribué sous serre pour des conditions de nuit en présence de chauffage et pour des conditions diurnes où l'influence du rayonnement solaire est prépondérante. La thèse s'appuie sur des approches expérimentales et numériques qui ont pour ambition d'être appliquées à deux situations, hivernales et pré-estivales correspondant aux problématiques auxquelles les producteurs doivent faire face : gestion de la température et de l'humidité avec une préoccupation de réduction des consommations énergétiques. Des campagnes de mesures ont été menées au printemps 2010 et à l'hiver 2011 pour, d'une part, améliorer la compréhension des interactions entre les paramètres climatiques, et d'autre part, fournir des données d'entrée et de validation aux modèles numériques. Les campagnes de mesures de nuit (hiver 2011) en présence de chauffage ont mis en évidence les besoins plus importants en chauffage et les températures de toiture plus basses lors des nuits claires. Les campagnes de mesure menées au printemps 2010 ont permis mettre en évidence l'influence prépondérante du rayonnement solaire sur la température et la transpiration du végétal ainsi que le phénomène de stockage-restitution de chaleur dans le sol. Enfin, la forte sensibilité du taux d'humidité ainsi que la moindre dépendance de la température vis-à-vis de la " qualité " de la ventilation a été démontrée. Les données expérimentales sont exploitées afin de mettre en œuvre les modèles de transpiration adaptés aux conditions de culture sous serre. Les bilans hydriques réalisés sur le végétal divisé en deux couches fournissent des informations précieuses sur la distribution de la transpiration et des contributions de chaque couche. Les simulations numériques sont réalisées en 2D instationnaire, en application sur une situation nocturne (nuit claire et nuit nuageuse) et une situation diurne (journée parfaitement ensoleillée). Elles intègrent deux éléments novateurs : la résolution de l'ETR dans le couvert végétal et l'utilisation d'une routine permettant de simuler la condensation. L'évolution des paramètres climatiques simulés est validée sur la base des données expérimentales. Il apparait cependant que les écarts de température et d'humidité entre l'air dans et au dessus du végétal sont sous estimés. Enfin, les modèles numériques sont exploités pour réaliser des bilans hydriques et énergétiques.

Les moteurs à plusieurs degrés de liberté angulaire sont des systèmes utilisés dans de nombreuses applications avec par exemple la réalisation d'articulation de robots, de roues omnidirectionnelles ou la stabilisation de composants optiques. La méthode classiquement utilisée pour réaliser ce type de plateformes limite cependant la capacité de ces systèmes à être miniaturisés. Une solution facilitant la réalisation de systèmes compacts à plusieurs degrés de liberté angulaire est l'utilisation d'un unique rotor sphérique. Nous proposons dans ce travail, un actionneur ultrasonique piézoélectrique adapté à la fabrication de moteurs sphériques compacts. L'innovation de l'actionneur développé tient dans l'adaptation de sa surface de contact au rotor, ce qui permet d'obtenir une surface étendue de contact et un guidage total du rotor sphérique. L'état de l'art des différents systèmes d'actionnement sphérique à un seul rotor est proposé ainsi que l'ensemble des considérations théoriques utiles à la modélisation et la réalisation de moteurs ultrasoniques. Un modèle complet pour la conception d'actionneurs ultrasoniques est proposé, il décrit l'étude de la géométrie et de ses modes de vibration, l'analyse du contact par l'application d'un modèle de transfert de force et la conversion électromécanique par méthodes semi-analytiques. Un prototype expérimental est proposé ainsi qu'une caractérisation de ses performances. Les données obtenues sont comparées au modèle théorique et permettent de valider les méthodes de conception proposées. L'effet de l'amplitude de stimulation, de la précontrainte et de l'étendue du contact y sont notamment analysés. Bien qu'appliqué au cas particulier d'un actionneur multimode incurvé, la démarche employée est générale et peut-être adaptée à la réalisation d'autres systèmes ultrasoniques.

Un sillage est la zone de perturbations turbulentes que laisse derrière lui tout corps en mouvement dans un fluide. Parmi les sillages de corps profilés auxquels on s'intéresse ici, il convient de différencier les sillages de traînée, des sillages de corps propulsés. Dans ce dernier cas, la connaissance des caractéristiques de l'écoulement présente un intérêt certain pour la détection des véhicules. L'objectif de cette étude est d'apporter les informations expérimentales les plus complètes possibles sur l'évolution du champ de vitesse. Pour cela nous avons mis en place une expérience en soufflerie aérodynamique. Le nombre de Reynolds calculé à partir de la vitesse de l'écoulement amont et du diamètre de la maquette est de 5,8.104. Un système d'acquisition automatisé des trois composantes du vecteur vitesse instantanée qui utilise une sonde anémométrique à films chauds triples a été développé spécialement. Une exploration axiale comparative entre les sillages pur et autopropulsé du même corps est effectuée pour un vaste domaine qui s'étend du sillage proche au sillage lointain (de 0,2 à 50 diamètres de maquette en aval du corps) ce qui constitue un ensemble très complet de données. La connaissance de toutes les caractéristiques temporelles du champ de vitesse permet de disposer : - des trois composantes de la vitesse moyenne - des six composantes du tenseur de Reynolds - des moments d'ordre 3 de la vitesse - des densités spectrales de puissance selon les directions axiale, radiale et tangentielle - des trois composantes de la vitesse moyenne et des six composantes du tenseur de Reynolds, analysés en phase avec la rotation de l'hélice. Un résultat original de cette étude est la mise en évidence expérimentale de l'existence d'une affinité pour un sillage lointain de corps autopropulsé par hélice. Cependant, les lois d'évolution fournies par l'analyse des équations du mouvement et reposant sur la distinction entre deux régimes extrêmes d'écoulement ne paraît pas justifiée expérimentalement. Nous avons déterminé tous les termes de l'équation d'énergie cinétique turbulente accessibles à la mesure. On a ainsi pu montrer l'affinité des échanges turbulents en champ lointain dans les deux configurations. Pour le sillage de corps autopropulsé par hélice, le fait que la production d'énergie cinétique turbulente soit négligeable dans la zone affine est un résultat expérimental important. En reprenant l'analyse dimensionnelle du problème et en ajoutant à l'approximation dite de la " couche limite " la condition d'autopropulsion, on montre que ce terme de production est effectivement négligeable devant les autres. Dès la zone intermédiaire, la forme de la densité spectrale de puissance relative se conserve selon une direction radiale, en sillage pur ou autopropulsé. Cela traduit la conservation en proportion de la répartition de l'énergie entre les différentes structures turbulentes. Dans le cas du sillage autopropulsé, l'analyse des densités spectrales de puissance met en évidence un signal " large bande " auquel viennent s'ajouter des raies spectrales. La zone de sillage proche, jusqu'à deux diamètres en aval de la maquette, est très fortement marquée par ces phénomènes d'excitation induits par le passage périodique des pales. Ce constat nous a conduits à effectuer un traitement statistique de la vitesse en phase avec la rotation de l'hélice. On a obtenu une visualisation des structures organisées de l'écoulement présentes pour une région qui s'étend jusqu'à deux diamètres du corps en aval de celui-ci.

Lors du procédé de cokéfaction, en raison de la faible largeur de la chambre de carbonisation des fours modernes, l'expansion horizontale de la pâte à coke génère une poussée sur les parois de chauffage. L'objectif de cette thèse, qui s'inscrit dans le cadre du projet européen " Swelling Pressure in a Coke Oven, Transmission on Oven walls ", est de mieux comprendre le phénomène de poussée des charbons lors de la cokéfaction et de développer un modèle permettant d'anticiper ce phénomène. Pour cela, un modèle phénoménologique prenant en compte les phénomènes physico-chimiques en présence a été développé. Une mise en équation originale est proposée pour la cinétique de condensation des goudrons et le gonflement des grains de charbon lors de la pyrolyse. Le modèle proposé est le premier reproduisant simultanément la poussée sur les piédroits et la pression des gaz produits lors de la cokéfaction. Les résultats de simulation du cas particulier de la cokéfaction du charbon Blue Creek dans le four pilote du Centre de Pyrolyse de Marienau et les mesures de pression, de température et de poussée réalisées lors des essais présentent des écarts mais demeurent en bon accord. Ces écarts sont essentiellement dus à la méconnaissance des propriétés du charbon et de son comportement mécanique. L'hypothèse d'un comportement élastique linéaire entraîne une surestimation de la poussée. L'étude de sensibilité amène, entre autres, à la conclusion que la poussée ne dépend pas directement de la pression des gaz et que le gonflement des grains de charbon joue un rôle déterminant.

A numerical model to predict train induced vibrations is presented. The dynamic computation considers mutual interactions in vehicle/track coupled system by means of a finite and discrete elements method. The vehicle is modeled by 7 bi-dimensional rigid elements representing: the body, the two boogies and the four wheels. The railway is discretized as finite Timoshenko beam elements. Axial deformation is assumed insignificant. The substructure is made-up of: rail-pads, sleepers, ballast, and background. Rail-pads are modeled as spring/damper couples without mass and sleepers are modeled as rigid elements. The rail-sleeper contact is assumed extended to a connection-area, rather than a single point assumption. To model this area many spring/damper couples are disposed along the length of sleepers. The ballast is modeled as blocks of mass made-up of rigid elements, connected to sleeper by spring/damper couples. To allow the transmission of vibrations in longitudinal direction too, spring/damper couples connect ballast elements horizontally.The dynamic interaction between the wheel-sets and the rail is accomplished by using the non-linear Hertzian model with hysteresis damping. The rail defects and the case of out-of-round wheels are considered too.A modal analysis of supporting structure is done to validate the substructure model comparing it to experimental data.Comparisons between numerical results of our model, experimental data and numerical results of others literature models are done on contact-force, rail accelerations and sleepers accelerations to validate the coupled vehicle/track system.Moreover a modal analysis of the coupled vehicle/track system is done to analyze the relationship between resonance frequencies, train velocities and ballast displacements.A sensitivity analysis is done to evaluate the variables more affecting the maintenance costs. The parameters more conditioning the ballast maintenance costs are the ballast modulus and the train mass.The effects of train velocity on the ballast displacements are analyzed in relationship with substructure properties. A new formulation to evaluate the railway toll connected to ballast wear is introduced.A new interpretation of the critical velocity in the range 100-300 km/h is proposed.

Ce mémoire de thèse a pour objectif de caractériser et de modéliser le comportement mécanique du béton à l'échelle mésoscopique. Le cadre plus général de cette étude est le développement d'un modèle mésoscopique du béton; ce modèle consiste à représenter le béton comme un milieu hétérogène en prenant en compte la différence qui existe entre les granulats et la pâte de ciment en respectant la courbe granulométrique qui joue un rôle primordial dans le comportement du béton. Les paramètres du modèle sont ceux qui décrivent les comportements mécanique et thermique de la pâte de ciment et des granulats. Nous souhaitons ainsi nous intéresser à la compréhension du comportement du béton, celui-ci étant considéré comme une structure. Un programme de tirage aléatoire de la structure granulaire valable en 2D et en 3D a été développé. Ce programme est interfacé avec le code de calcul Cast3M qui permet d'effectuer les simulations numériques. Une méthode de représentation numérique des inclusions du béton a été également développée et validée par projection de la géométrie sur les fonctions de forme, éliminant ainsi les problèmes de maillage qui rendaient la représentation de l'ensemble du squelette granulaire quasi-impossible lorsque on souhaite représenter des petits et des gros granulats, particulièrement en 3D. Dans un premier temps, le modèle est utilisé en 2D et 3D pour optimiser le modèle géométrique de structure interne du béton, la stratégie de maillage et la plus petite taille des granulats à modéliser. Les résultats du modèle 2D et 3D sont analysés et comparés dans le cas de sollicitations de traction et de compression uniaxiale. Le modèle utilisé est un modèle d'endommagement isotrope unilatéral développé par Fichant. Le modèle permet de simuler à la fois le comportement macroscopique mais aussi local avec l'étude de la distribution de la fissuration et de l'ouverture des fissures. Le modèle montre des résultats intéressants sur la transition entre l'endommagement diffus et localisé et est capable de reproduire les effets de dilatance en compression. Finalement, le modèle mésoscopique a été appliqué dans les trois cas de simulations : le calcul de la perméabilité du béton fissuré; la simulation de l'hydratation du béton au jeune âge et enfin les calculs de structures avec la mise en évidence de l'effet d'échelle sur des poutres entaillées.

La simulation numérique des systèmes multicorps en présence d'interactions complexes, dont le contact frottant, pose de nombreux défis, tant en terme de modélisation que de temps de calcul. Dans ce manuscrit de thèse, nous étudions deux familles de décomposition de domaine adaptées au formalisme de la dynamique non régulière des contacts (NSCD). Cette méthode d'intégration implicite en temps de l'évolution d'une collection de corps en interaction a pour caractéristique de prendre en compte le caractère discret et non régulier d'un tel milieu. Les techniques de décomposition de domaine classiques ne peuvent de ce fait être directement transposées. Deux méthodes de décomposition de domaine, proches des formalismes des méthodes de Schwarz et de complément de Schur sont présentées. Ces méthodes se révèlent être de puissants outils pour la parallélisation en mémoire distribuée des simulations granulaires 2D et 3D sur un centre de calcul haute performance. Le comportement de structure des milieux granulaires denses est de plus exploité afin de propager rapidement l'information sur l'ensemble des sous domaines via un schéma semi-implicite d'intégration en temps.

L'objectif de cette étude est d'analyser l'écoulement traversant le récepteur solaire haute température de la centrale PEGASE. Cet écoulement est turbulent et soumis à un très fort gradient de température. Nous utilisons la simulation des grandes échelles thermiques appliquée aux équations bas-Mach pour étudier l'écoulement turbulent dans une géométrie de canal plan bi-périodique avec températures imposées aux parois. Pour cette étude, nous avons réalisé des simulations pour deux nombres de Reynolds turbulents différents (Re_tau m = 180 et Re _taum = 395) et pour quatre rapports de température (T2=T1 = 1 ; 1; 04 ; 2 et 5). Après avoir validé notre modèle, nous avons réalisé une étude sur la modélisation sous-maille thermique qui conclut sur la nécessité d'utiliser un modèle sous-maille thermique dynamique pour un écoulement ayant une forte intensité turbulente soumis à un très fort gradient de température. Nous analysons l'influence du gradient de température sur les différents profils de vitesse et de température ainsi que sur les spectres d'énergie. Augmenter le gradient de température crée une dissymétrie de tous les profils. Une relaminarisation du côté chaud est visible. On remarque une nouvelle répartition des fluctuations à travers le canal ainsi que de l'énergie en fonction de la taille des échelles. Une création de fluctuations de vitesse et de température a également été mise évidence du côté froid et du côté chaud du canal. Ces effets sont dus, à la fois à une variation de la viscosité cinématique (effet direct de la température) et à une interaction entre les champs turbulents dynamique et thermique.