On s'interesse dans ce travail a la stabilite des ondes de <br />choc pour des systemes hyperboliques de lois de conservation <br />multidimensionnels. Ce probleme a ete traite par Andrew <br />Majda sous une hypothese, dite de stabilite uniforme, qui <br />intervient de facon cruciale dans son analyse. Cette hypothese <br />est cependant mise en defaut dans certains exemples, par exemple <br />dans l'etude des transitions de phase liquide-vapeur. Nous <br />examinons ici la stabilite des interfaces qui ne verifient <br />pas l'hypothese de stabilite uniforme, et montrons comment les <br />resultats de Majda s'etendent a de telles discontinuites. <br /><br /><br />On commence par montrer la stabilite lineaire des chocs plans <br />faiblement stables, a l'aide d'un symetriseur de Kreiss <br />degenere qui tient compte des modes neutralement instables. <br />Cette premiere etape etablit un compte precis des pertes de <br />derivees intervenant dans les estimations d'energie. Dans un <br />second temps, nous montrons que ces estimations d'energie demeurent <br />valables lorsque l'on etudie la stabilite des interfaces (non <br />planes) proches d'un choc plan. L'utilisation du calcul <br />paradifferentiel nous permet de traiter des perturbations <br />peu regulieres du choc plan initial. Sous une hypothese de <br />petitesse sur le comportement global des courbes bicaracteristiques, nous montrons une estimation d'energie semblable a celle etablie pour le probleme linearise a coefficients constants. Ce resultat devrait permettre de montrer la stabilite non lineaire des ondes de choc faiblement stables.

Cette thèse traite de la modélisation des infrastructures autoroutières et de leur gestion par des méthodes de régulation telles que le contrôle d'accès. L'approche retenue est macroscopique et conduit à des modèles distribués sous forme d'équations aux dérivées partielles non linéaires. Nous apportons plusieurs éclairages sur l'analyse et la résolution de ces modèles (condition d'entropie pour les rampes d'accès, discrétisation simpliée) et proposons une interprétation hybride des inhomogénéités (conditions aux limites, rampes d'accès et de sorties, variations brutales des paramètres) adaptée aux problèmes de contrôle. Deux nouvelles méthodologies calculatoires sont ensuite introduites pour concevoir des contrôleurs dynamiques s'appliquant à la gestion du trafic. La première est formulée comme un problème de commande optimale en boucle ouverte et nécessite l'adaptation de la méthode adjointe traditionnelle en raison de l'irrégularité des solutions. La seconde repose sur une discrétisation sous la forme d'un système affine commuté et une synthèse boucle fermée utilisant la dissipativité et les inégalités matricielles linéaires.

.

Les lois d'endommagement incrémentales présentent de nombreux avantages en terme de modélisation thermomécanique. S'affranchissant de la notion de cycle en fatigue, elles s'appliquent naturellement aux chargements complexes, anisothermes. Des extensions à l'anisotropie du modèle d'endommagement isotrope de Lemaitre sont proposées et identifiées sur le Haynesl88, superalliage à base cobalt utilisé pour la réalisation de chambres de combustion de turbomachines, permettant de déterminer les conditions d'amorçage de fissure par la méthode des éléments finis. Le cas des hautes températures est traité via un couplage viscoplasticité/endommagement dans le cadre de la thermodynamique des matériaux solides. Le seuil d'endommagement en énergie stockée et l'extension du modèle aux conditions unilatérales de refermeture des micro défauts sont présentés. Différents schémas numériques de résolution sont proposés, dans le cas de calculs couplés, ou de post traitements de calculs 3D viscoplastiques sans endommagement. Des applications anisothermes sont simulées: fluage à température variable, chargements biaxiaux, non proportionnels, aléatoires... Les apports de l'anisotropie de l'endommagement et des conditions unilatérales de refermeture des microdéfauts sont étudiés. Enfin, la robustesse du modèle est testée sur un essai original de fatigue thermique structural sur éprouvette multiperforée, reproduisant les sollicitations subies par une pièce réelle. La corrélation entre résultats numériques et expérimentaux en terme de comportement et de durée de vie permet la validation de la méthodologie d'identification et d'utilisation du modèle incrémental, en vue de son industrialisation.

.

Durant un tremblement de terre, la rupture grandit et se propage sur la faille. Lorsqu'elle s'arrête, le séisme atteint sa taille finale. Comprendre ce qui, dans la nature, détermine la capacité à se propager ou à s'arrêter de la rupture est fondamental en sismologie, puisqu'il existe des séismes de toutes les tailles, mais que seuls les plus grands sont dévastateurs. Dans cette thèse, nous étudions l'impact de différentes façons d'arrêter ou de perturber la propagation de la rupture, à travers des études numériques dynamiques. <br /><br />Tout d'abord, nous avons inclus une limite à l'élasticité du milieu entourant la faille, afin de simuler la propagation de la rupture dans un milieu fracturé, qui dissipe une partie de l'énergie libérée. Nous avons, pour la première fois, inclus et étudié l'impact de cette dissipation dans un modèle de rupture 3D. La rupture, dans ces conditions, est beaucoup plus sensible aux barrières, et s'arrête plus facilement. La cinématique de la rupture est remarquablement modifiée (vitesse de propagation plus lente, vitesse de dislocation maximale limitée). Les mouvements engendrés en surface sont atténués. La plasticité est, de fait, un phénomène crucial à prendre en compte dans la modélisation de la rupture sismique.<br /><br />Dans une deuxième partie, nous avons étudié l'impact d'une hétérogénéité spatiale de résistance à la rupture sur la faille. L'introduction de l'hétérogénéité permet d'obtenir des profils de glissement dont la forme se rapproche des formes observées dans les cas naturels. La propagation et l'arrêt de la rupture perdent leur caractère prédictible lorsque ce sont de petites barrières qui arrêtent la rupture progressivement. On peut obtenir une grande variété de tailles d'événements pour une même statistique de taille des barrières. L'obtention d'une loi puissance de type Gutemberg-Richter sur toute la gamme des tailles d'événements est conditionnée par l'augmentation progressive de l'énergie de fracturation avec la taille de la rupture, d'une façon similaire à ce qui est obtenu en considérant un comportement plastique du milieu. <br /><br />Enfin, l'étude des relations glissement final - taille de l'aspérité rompue, dans des modèles lisses de type aspérité/barrière, a montré que la dynamique contrôlait une partie de la loi d'échelle du glissement maximum, et que la segmentation des failles modifie sensiblement la loi d'échelle.

Nous étudions les propriétés de modèles de séries chronologiques vectorielles à composantes binaires en vue de la modélisation du phénomène climatologique de la succession des jours selon leur caractère sec ou humide dans un réseau de stations de mesures. Nous menons une étude expérimentale par simulations des modèles considérés et envisageons leur utilisation dans des essais de modélisation concernant certains sous-réseaux francais de stations météorologiques

Cette thèse expérimentale porte sur l'évaporation d'une gouttelette mouillante : il s'agit d'un problème de ligne de contact mobile dans lequel deux processus de relaxation aux effets antagonistes (étalement et évaporation) s'opèrent de façon concomitante. <br /> Le travail qui est rapporté fait suite à la thèse de C. Poulard soutenue en 2005 au sein de notre groupe. Il s'est dégagé de ce premier travail de grandes tendances : notamment, la phase de rétraction est caractérisée par l'existence de lois de puissance pour le rayon et l'angle de contact. Ce caractère particulier a donné lieu à un modèle élaboré en collaboration avec M. Ben Amar. Le domaine de validité de ce modèle peut toutefois être amélioré. Le manuscrit prend ainsi la forme d'une discussion illustrée de faits expérimentaux, où les hypothèses sont reprises, les processus invoqués, testés et au cours de laquelle l'accent est mis sur les points importants en vue d'une meilleure description de la dynamique de ces objets volatils.<br /> Soulignons qu'en raison de la forme particulière du taux d'évaporation, la ligne de contact recule alors que l'angle de contact n'est pas nul. Or, cet angle est particulièrement significatif dans les dynamiques. Un accent particulier est porté sur la forme de l'interface liquide au niveau de la ligne de contact, gouvernée par un réel problème d'hydrodynamique. Notamment, des tests précis sur le couplage entre l'évaporation et les écoulements hydrodynamiques à l'extension maximale sont menés et les résultats de cette étude sont appliqués avec succès à la remise à l'échelle des dynamiques.<br /> Nous abordons dans une seconde partie, la question peu typique de l'évaporation d'une goutte de mélange. Dans un tel système, l'évaporation provoque un défaut du composé le plus volatil au niveau de la ligne de contact. Nous verrons dans cette partie que les hétérogénéités ainsi créées sont à l'origine de gradients de tension de surface qui contrôlent les dynamiques.

Nous étudions des méthodes d'approximation particulaire de solutions d'équations aux dérivées partielles décrivant l'état macroscopique de certains systèmes physiques. Elles consistent en l'introduction d'un grand nombre N de particules fictives évoluant selon des équations différentielles couplées, ordinaires ou stochastiques, dans un sens plus simple à résoudre que l'équation macroscopique; l'état de ce système de particules est décrit par une mesure de probabilité, dite mesure empirique. La validité de la méthode est donnée par la convergence, quand N tend vers l'infini, de cette mesure empirique vers la solution macroscopique originale, appelée limite de champ moyen. Nous cherchons principalement à en donner des estimations explicites, quantifiant ainsi la précision de l'approximation.<br /><br />Dans ce cadre nous étudions l'approximation des équations de transport de Vlasov et d'Euler par des systèmes de particules déterministes en interaction. Le problème de la convergence de la méthode se ramène à un problème de stabilité de solutions que nous traitons par des propriétés de type contraction pour des distances (de Wasserstein) liées à la théorie du transport optimal de mesures. Nous établissons aussi une propriété analogue de contraction pour des lois de conservation scalaires. <br /><br />Nous étudions également l'approximation d'équations de diffusion de McKean-Vlasov par des systèmes de particules stochastiques. Nous en donnons l'erreur de manière quantitative à l'aide de techniques de couplage, d'estimations de propagation du chaos et d'inégalités de concentration ou de déviation.<br /><br />De façon plus systématique nous nous intéressons à de telles inégalités de concentration pour des mesures de probabilité et à leurs relations avec des inégalités de transport (liant distances de Wasserstein et entropie) et de Sobolev logarithmiques. En particulier nous établissons de telles inégalités pour certaines classes de lois de variables dépendantes.

Le but de cette thèse est de trouver une méthodologie numérique pour simuler le comportement et la fragmentation des matériaux quasi-fragiles, tels que les roches et le béton, sous sollicitations dynamiques. Une telle méthodologie aidera à une meilleure compréhension des mécanismes de fragmentation des roches à l'explosif et, par conséquent, contribuera dans l'optimisation des performances du tir à l'explosif. Dans cette étude, on considère que le processus de fragmentation est une extension naturelle de celui de la rupture. Afin de décrire ce dernier processus, un modèle de comportement phénoménologique adapté aux matériaux considérés est développé et implémenté dans un code de calcul par éléments finis. Quant à la fragmentation dynamique, elle est traitée par une analyse en post-traitement basée sur l'historique de l'état thermodynamique du matériau. Pour reproduire le comportement macroscopique des matériaux quasi-fragiles incluant l'anisotropie induite par le chargement, la méthode des variables internes basée sur la thermodynamique des milieux continus est employée. Une variable interne de type scalaire est introduite pour modéliser, au niveau macroscopique, l'adoucissement du matériau suite à des chargements de compression. Sous des chargements de traction, un tenseur symétrique du second ordre est utilisé pour décrire l'endommagement anisotrope induit. Sous des chargements complexes, ces deux modèles sont couplés et l'effet de fermeture-réouverture des fissures est également traité. L'intégration du modèle développé dans le code éléments finis VIPLEF3D a conduit à l'élaboration d'une méthode de relaxation caractérisée par une actualisation explicite des variables d'état. Concernant la fragmentation dynamique, à partir d'essais de fragmentation au laboratoire, une formulation générale permettant de prédire la distribution des fragments est fournie. Dans cette formulation une taille moyenne de fragments est liée à une grandeur mécanique, supposée être l'origine du processus de fragmentation, par une fonction intrinsèque qui peut être identifiée en utilisant des essais adéquats de fragmentation. Cette grandeur mécanique est donnée par la résolution du problème aux limites où le modèle rhéologique est utilisé. Enfin, l'approche complète est alors appliquée à la modélisation des essais de rupture et de fragmentation en chambre, d'échantillons cylindriques de roches.