Les turbines à gaz doivent satisfaire des normes d'émission polluantes toujours en baisse. La formation de polluants est directement liée à la qualité du mélange d'air et de carburant en amont du front de flamme. Ainsi, leur réduction implique une meilleure prédiction de la formation du spray et de son interaction avec l'écoulement gazeux. La Simulation aux Grandes Échelles (SGE) semble un outil numérique approprié pour étudier ces mécanismes. Le but de cette thèse est de développer des modèles phénoménologiques décrivant la phase liquide notamment le film et son atomisation en bout de lèvre d'injecteur, dans un contexte SGE. Ces modèles sont validés ou calibrés sur l'expérience académique réalisée par l'Institut für Thermische Strömungsmaschinen (ITS) de l'université technologique de Karlsruhe (KIT), et appliqués dans une configuration réelle de moteur d'hélicoptère. Dans un premier temps, le film liquide mince est décrit par une approche Lagrangienne. Les particules de film représentent un volume élémentaire de liquide adhérent à la paroi. L'équation du mouvement est donnée par l'intégration des équations de Saint-Venant sur l'épaisseur du film. La dynamique du film est donnée par le gradient de pression longitudinal, le cisaillement interfacial du gaz et la gravité. Dans un second temps, l'atomisation du film est caractérisée par la distribution de taille de gouttes du spray généré. Celle-ci est décrite par une distribution de Rosin-Rammler dont les coefficients sont paramétrés par la vitesse du gaz, la tension superficielle du liquide et l'épaisseur de la lèvre de l'injecteur. Les constantes de ce modèle, baptisé PAMELA, sont calibrées sur l'expérience du KIT-ITS. La simulation de l'expérience KIT-ITS permet de valider le modèle de film, de vérifier la robustesse du modèle PAMELA, et de comparer qualitativement l'angle du spray généré. L'application de ces modèles dans une configuration réelle partiellement instrumentée permet de valider PAMELA sans modification de ses constantes, et d'étudier leur impact sur la structure de flamme, comparé aux méthodes traditionnelles d'injection liquide.

Le présent travail porte sur l'étude expérimentale, analytique et numérique du comportement mécanosorptif et viscoélastique couplé du bois. Le comportement mécanosorptif intègre notamment le phénomène hygroverrou qui se caractérise par un blocage des déformations en phase de séchage sous contrainte. Afin d'aborder cette étude, la structure, la microstructure ainsi que le comportement hydromécanique du bois sont rappelés. Parallèlement, la théorie de la viscoélasticité linéaire ainsi que les méthodes analytiques propres à la résolution de tels problèmes sont présentées. Ensuite, des essais mécanosorptifs préliminaires de caractérisation des propriétés mécaniques sont effectués sur du Sapin Blanc (Abies alba Mill.) du Massif central grâce à un dispositif expérimental approprié. La description du dispositif expérimental, du protocole utilisé ainsi que les dimensions des échantillons testés sont présentées. Des essais mécanosorptifs en relaxation et en fluage sous chargement et humidité relative variables sont également effectués. En particulier, deux essais en relaxation sous déformation imposée constante et humidité relative variable et un essai de fluage sous contrainte et humidité relative variables de façon cyclique, sont exposés. Sur la base des données issues des essais préliminaires et viscoélastiques, la déformation hygroverrou est isolée via la partition de la déformation selon des parts élastique, visqueuse et hydrique. Par la suite, un modèle analytique est proposé pour la modélisation de la déformation hygroverrou. Un modèle pseudo élastique, basé sur les lois de Hooke (en humidification) et de Hooke restreint (en séchage), est formulé en rigidité et en complaisance. Dans ce cas, plusieurs modèles rhéologiques viscoélastiques mécanosorptifs sont proposés. Ces modèles sont classés en fonction du couplage de leurs effets et du principe de partition formulés en contrainte ou en déformation. Des lois de comportement basées sur ces modèles et formulées sous forme d'équations intégrales et incrémentales adaptées à la mise en oeuvre numérique sont finalement développées. Enfin, des exemples illustratifs montrent l'équivalence des deux expressions pseudo élastiques formulées en rigidité et en complaisance ; parallèlement, une validation des modèles viscoélastiques mécanosorptifs proposés, en relaxation et en fluage, par comparaison aux résultats expérimentaux est effectuée. Ces résultats montrent l'efficacité des modèles établis ainsi que la pertinence des essais expérimentaux réalisés pour la mise en évidence et la prise en compte des comportements mécanosorptif et viscoélastique couplés.

Le CEA travaille sur des matériaux poreux - alvéolaires, composites, céramiques, etc. - et cherche à optimiser leurs propriétés pour des utilisations spécifiques. Ces matériaux, souvent composés de plusieurs constituants, ont en général une structure complexe avec une taille de pores de quelques dizaines de microns. Ils sont mis en oeuvre dans des systèmes de grande échelle, supérieure à leurs propres échelles caractéristiques, dans lesquels on les considère comme équivalents à des milieux homogènes, sans prendre en compte sa microstructure locale, pour simuler leur comportement dans leur environnement d'utilisation.Nous nous intéressons donc à la caractérisation des propriétés thermiques effectives de matériaux à microstructure hétérogène en cherchant à déterminer par méthode inverse en fonction de la température la diffusivité thermique qu'ils auraient s'ils étaient homogènes.L'identification de la diffusivité de matériaux poreux et/ou semi-transparents est rendue difficile par le couplage conducto-radiatif fort qui peut se développer rapidement dans ces milieux avec une augmentation de la température. Nous avons donc modélisé le transfert de chaleur couplé conducto-radiatif en fonction de la température au sein de matériaux poreux multiconstituants à partir de leur microstructure numérisée en voxels. Notre démarche consiste à nous appuyer sur la microstructure 3D obtenue par tomographie. Ces microstructures servent de support numérique à cette modélisation qui permet d'une part de simuler tout type d'expériences thermiques numériques - en particulier la méthode flash dont les résultats nous permettent de déduire la diffusivité thermique -, et d'autre part de reproduire le comportement thermique de ces échantillons dans leur condition d'utilisation.

La simulation des écoulements cavitants est confrontée à des difficultés de modélisation et de résolution numérique provenant des caractéristiques particulières de ces écoulements : changement de phase, gradient de masse volumique important, variation du nombre de Mach, turbulence diphasique, instationnarités. Dans cette thèse, nous nous sommes appliqués à dériver proprement le modèle de mélange homogène 1-fluide couplé à une modélisation RANS de la turbulence. A partir des termes contenus dans ces équations et de la nature des écoulements cavitants étudiés, plusieurs modèles de turbulence basés sur la notion de viscosité turbulente ont été testés : modèles faiblement non-linéaires (corrections SST et de réalisabilité), ajout des termes de turbulence compressible, application de la correction de Reboud, modèles hybrides RANS/LES (DES, SAS). Ces modèles ont été incorporés dans un code compressible qui fait appel à une résolution implicite en pas de temps dual des équations de conservation avec une technique de pré-conditionnement bas-Mach pour traiter les zones incompressibles. Les simulations 2D et 3D ont porté sur deux géométries de type Venturi caractérisées par la présence d'une poche de cavitation instationnaire due à l'existence d'un jet rentrant liquide/vapeur le long de la paroi. Elles montrent que l'ensemble des modèles sont capables de capturer le jet rentrant. En revanche, la dynamique de la poche varie entre les modèles et le manque de données expérimentales ne permet pas de discriminer les modèles entre eux. Il apparaît à la vue des résultats que les approches avec la correction de Reboud ou les modèles SAS améliorent la simulation des écoulements.

On étudie quelques problèmes singuliers d'homogénéisation en élasticité linéarisée ou électricité non linéaire où en plus de la période de distributions des hétérogénéités apparaissent d'autres paramètres comme la très forte (ou faible) rigidité (ou permittivité) d'une des phases et la taille relative de celle-ci. Sont visés des milieux fibrés ou stratifiés et des maçonneries planes et minces.

Le rampement est un mode de locomotion fondamental de nombreuses cellules eucaryotes, engagé dans des mécanismes aussi importants que embryogenèse, la réponse immunitaire et la cicatrisation. Son dérèglement provoque de graves maladies, en particulier des cancers. La compréhension mécanique de ce mode de locomotion présente également un grand intérêt pour la confection de robots opérant à l'échelle cellulaire. Le schéma classique du rampement cellulaire met en jeu la polymérisation du réseau d'actine dans la partie frontale de la cellule couplée avec l'activation des points d'adhésion focaux liant la cellule à son substrat, alors que la partie postérieure de la cellule se détache du substrat sous l'effet de la contraction engendrée par les molécules de myosine. De manière simplifiée, on peut voir la partie motrice d'une cellule eucaryote comme un gel actif dont les fonctions sont contrôlées par des processus chimiques et mécaniques. En particulier, les mouvements coordonnés de ce gel engendrant le rampement impliquent une auto organisation spatiale et temporelle du cytosquelette et demandent un apport continu d'énergie. Si les bases biochimiques de la motilité cellulaire sont connues, la compréhension qualitative des interactions mécaniques entre les différents acteurs rentrant en jeu dans le rampement n'est que très limitée, et ce malgré les récents efforts visant à construire des modèles complets et exhaustifs du phénomène. Cette thèse présente l'analyse d'un modèle simple et unidimensionnel expliquant le rampement cellulaire. La première partie de la thèse est dédiée à l'analyse inverse du problème d'optimisation en vitesse et en efficacité mécanique du rampement. Notre analyse montre que les distributions optimales de contraintes contractiles et de friction avec le substrat sont en bonnes accord avec les distributions observées. Dans une seconde partie, nous proposons un mécanisme de motilité cellulaire spontanée centré sur la contraction et ignorant la polymérisation et la dépolymérisation de l'actine. A l'origine de la polarisation, l'anti-diffusion auto amplifiée des moteurs pilotant la contraction déstabilise la configuration initialement symétrique de la cellule. L'apparition de cette instabilité morphologique est pilotée par le ratio entre la diffusion et la contractilité des moteurs générant un flot convergent qui, lui-même, transporte les moteurs. Par l'étude unidimensionnelle du phénomène, nous montrons que le flot ainsi produit peut générer un mouvement de translation de la cellule qui reproduit des observations concernant la motilité spontanée des fragments de keratocytes. La troisième partie de la thèse concerne la motilité cellulaire basée sur les propriétés de polymérisation et de dépolymérisation active de l'actine qui permettent non seulement l'autopropulsion de la cellule mais aussi le mécanisme de poussée (d'obstacles) et de tirée (de noyau cellulaire par exemple) de charges données . Nous utilisons un modèle minimaliste pour montrer que la relation force-vitesse dans le cas de la poussée est essentiellement réminiscente du mécanisme de protrusion piloté par la polymérisation alors que la relation force-vitesse du tirage d'une charge ne repose sur le mécanisme de protrusion que pour des charges faibles, le mécanisme de contraction prenant le relais pour des charges plus grandes.

La Fabrication Directe par Projection Laser (FDPL) est un procédé novateur qui permet la réalisation de pièces de géométrie complexe à partir d'un jet de poudre. Les pièces sont obtenues en superposant différentes couches de matière à l'aide d'une buse coaxiale avec un faisceau laser. Cette technique de fabrication, très rapide, reste cependant peu répandue, en particulier du fait de l'irrégularité des surfaces. L'objet de ce travail est de développer des modèles numériques capables de prédire le paramètre d'ondulations caractérisant l'état de surface à partir uniquement des paramètres opératoires et des propriétés thermophysiques du substrat et de la poudre, dans le but de mieux comprendre les mécanismes responsables de l'état de surface dégradé. Une première étude a porté sur la fusion locale d'un barreau vertical sous l'action d'un laser afin de valider la mise en données du modèle numérique. Ce modèle thermohydraulique inclut les effets de tension superficielle et utilise la méthode ALE (Arbitrary Lagrangian-Eulerian) pour suivre la déformation de la surface libre. Ce modèle prédictif a été transposé au procédé FDPL en intégrant l'apport de matière dû au jet de poudre. Les caractéristiques de ce jet de poudre sont obtenues grâce à un modèle 3D prédisant la trajectoire des particules et leur échauffement sous l'action du faisceau laser. Ces données servent à définir les conditions aux limites de modèles 2D thermohydrauliques du bain fondu. Nous avons ainsi retrouvé, pour un substrat mince, la corrélation établie expérimentalement entre l'amplitude des ondulations et le taux de dilution. Il a été montré que, dans le cas de l'alliage de titane Ti-6Al-4V, un meilleur état de surface est obtenu pour une forte puissance laser, une vitesse de défilement élevée et un faible débit massique. La simulation 2D de dépôts multicouches a permis d'étudier l'influence de la stratégie de balayage et du temps de pause entre couches. Une étude de sensibilité a également été menée afin d'analyser le rôle de paramètres tels que le coefficient thermocapillaire ou la viscosité. L'étude a par la suite été étendue à un cas 3D pour étudier l'effet de la distribution énergétique du faisceau laser ainsi que les propriétés thermophysiques. La prédiction des différents modèles numériques développés à l'aide du code de calcul COMSOL Multiphysics® est validée grâce à des données expérimentales.

Dans l'industrie d'extraction pétrolière, l'efficacité des séparateurs eau-huile pour la production offshore est cruciale. L'objet de ce travail est de mettre en place les outils numériques nécessaires à la modélisation du fonctionnement de ces systèmes. Les phénomènes physiques entrant en jeu sont principalement : la présence d'interfaces entre des fluides non miscibles, la viscosité de ces fluides, et les effets de tension superficielle. Les modèles physiques et numériques correspondants ont été implémentés dans le cadre de la méthode numérique SPH (Smoothed Particle Hydrodynamics) développée au L.M.F.. Cette méthode numérique appartient à la classe des méthodes particulaires (sans maillage), suivant une approche d'écoulement compressible et avec une résolution explicite. Pour modéliser au mieux les écoulements bifluides, la formulation historique de la SPH a été enrichie par deux approches différentes, développées simultanément. Chacune d'entre elles a été validée séparément. La physique supplémentaire a été rajoutée par des modèles communs qui ont été validés sur différents cas tests tels que l'écoulement de Poiseuille, les instabilités de Rayleigh-Taylor, des cas d'envahissement ou l'évolution de bulles dans un liquide. Ce dernier cas a permis la comparaison aux outils de conception utilisés dans le procédé d'ingénierie de SAIPEM S.A., par l'intermédiaire d'une validation sur la loi de Stokes. Finalement, les capacités de la méthode sont illustrées sur la séparation eau-huile dans un séparateur de géométrie simplifiée.

" Modélisation numérique multi-échelles des contacts lubrifiés - Aux frontières de la mécanique (des fluides, des solides), de la thermique et de la physique de la matière ". Derrière ce titre se cache un but double : apporter plus de physique aux modèles actuels de lubrification (rhéologie du lubrifiant, effets thermiques), et prendre en compte les effets locaux, difficilement modélisables dans le cadre de la mécanique des milieux continus. Les contraintes industrielles actuelles (concurrentielles, financières, environnementales) tendent en effet à réduire considérablement l'épaisseur du film lubrifiant, atteignant l'ordre du nanomètre. Une approche purement mécanique ne tient plus car elle omettrait des phénomènes importants : il est donc nécessaire d'envisager le contact lubrifié sous un nouvel angle. Cette démarche originale propose donc une approche à la fois discrète et continue (mêlant des modèles Eléments Finis et Dynamique Moléculaire) des contacts lubrifiés pour mieux rendre compte de la physique aux petites échelles : le lubrifiant confiné se stratifie, sa rhéologie n'est plus décrite par des lois volumiques, du glissement apparaît à l'interface fluide/surface solide, etc. Mes recherches ont notamment montré que la prise en compte du glissement à l'échelle du nanomètre (caractérisé par Dynamique Moléculaire) dans les modèles à l'échelle de contacts lubrifiés (de l'ordre du millimètre) conduisait à la fois à la formation de géométries particulières de l'interface et à une évolution du frottement corrélée aux observations expérimentales.

Ce travail de thèse s'insère dans les thèmes développés dans l'équipe "MécaWet" du laboratoire PMMH. Le fil conducteur de ces activités est l'étude de structures élancées (tiges, plaques, coques) en lien avec les autres thèmes développés au laboratoire. En particulier, nous nous intéresserons à la déformation de tubes flexibles par des forces capillaires, à l'écoulement de petits grains ou au contraire de gros objets dans des tubes souples et au flambement de ces tubes sous l'effet d'une dépression. L'approche expérimentale est fidèle à l'esprit du laboratoire. Nous développerons ainsi des expériences macroscopiques modèles afin d'extraire les ingrédients physiques d'un problème. Des expé- riences à taille humaine permettent en effet de contrôler facilement les paramètres expérimentaux sur plusieurs ordres de grandeur.