Tesi sul tema "Instabilités de type Weibel"

Les instabilités de type Weibel naissent si la distribution des vitesses du plasma présente une anisotropie. Elles entraînent la génération d’un champ magnétique dû à la formation de filaments de courant ainsi qu’une activité électrostatique importante. Ces phénomènes de base apparaissent dans de nombreuses situations, naturelles (vent solaire, jets relativistes) ou expérimentales (interaction laser-plasma) : les plasmas dans lesquels ils naissent peuvent être relativistes ou non, magnétisés ou non, collisionnels ou non, ce qui pose la question du choix du modèle à utiliser pour les décrire. La théorie cinétique est le cadre le plus complexe dans lequel nous travaillerons. De par sa complexité, il est intéressant de développer des modèles réduits. Un premier travail mené au cours de cette thèse est l’utilisation d’un modèle fluide incluant la dynamique du tenseur de pression pour modéliser la phase linéaire des instabilités de type Weibel. On discute le rôle essentiel joué par les composantes hors diagonale du tenseur dans la génération du champ magnétique, puis la capacité du modèle à reproduire quantitativement ou qualitativement les résultats cinétiques en introduisant la notion de limite hydrodynamique. La seconde partie de la thèse est ciblée sur le développement du code semi-lagrangien relativiste VLEM utilisant une méthode de décomposition de domaine : on présente les principales méthodes mathématiques utilisées dans le code, puis on aborde la problématique de la conservation de la charge à laquelle on apporte une réponse reposant sur une adaptation de la méthode d’Esirkepov. Le code est enfin validé grâce à plusieurs simulations d’instabilités de type Weibel<br>Weibel-type instabilities occurs when the velocity distribution function of the charged particles displays a pronounced anisotropy. A long-lasting magnetic field is generated due to the formation of current filaments, and it is accompanied by an important electrostatic activity. These ``basic’’ phenomena have been greatly investigated because of their involvement in many physical problems, natural (solar wind, relativistic jets) or experimental (laser-plasma interaction) : they occurs in plasmas which can be collisional or not, magnetised or not, relativistic or not. One needs to choose a suitable model for their description. The kinetic theory is the most complete and somewhat complex theoretical framework which we will consider. Due to its complexity, it may be interesting to develop reduced models. The first work realised during this thesis is the utilisation of a non-relativistic fluid description, including the dynamics of the pressure tensor, in order to model the linear Weibel-type instabilities. We put in evidence the effect of the non-diagonal components of the tensor on the magnetic field generation. We discuss the ability of the model to reproduce quantitatively or qualitatively the kinetic results by introducing the hydrodynamics limit. The second part of this thesis work is dedicated to the development of the relativistic semi-lagrangian code VLEM, using a domain decomposition scheme : we present the main mathematical tools used in the code, then we deal with the problem of the charge conservation and propose a solution for VLEM, based on an adaptation of the Esirkepov method. Finally, we validate the code through simulations of Weibel-type

L'origine de champs magnétiques observés dans les plasmas de laboratoire et d'astrophysique est l'un des problèmes récurrents en physique des plasmas. À cet égard, les instabilités de type Weibel sont considérées d'une grande importance. Ces instabilités ont pour origine une anisotropie de température (instabilité de Weibel) et des moments des électrons (instabilité de filamentation de courant). L'objectif principal de cette thèse est l'étude théorique et numérique de ces instabilités dans un plasma non collisionnel en régime relativiste. Le premier aspect de ce travail est l'étude du régime non-linéaire de ces instabilités et du rôle des effets cinétiques et relativistes sur la structure des champs électromagnétiques auto-cohérents. Dans ce cadre, un problème essentiel pour les applications et la théorie, concerne l'identification et l'analyse des structures cohérentes développées spontanément dans le régime non-linéaire sur des échelles cinétiques. Un deuxième aspect du travail est le développement de techniques analytiques et numériques pour l'étude des plasmas non collisionnels. Le modèle mathématique de référence, à la base des études des plasmas chauds, est le modèle Vlasov-Maxwell, où l'équation de Vlasov (théorie des champs moyens) est couplée aux équations de Maxwell de façon auto-cohérente. Un modèle unidimensionnel, le modèle multi-faisceaux, a également été introduit durant cette thèse. Basé sur une technique de réduction en dimension, il est à la fois un modèle analytique "simple" présentant l'avantage de pouvoir résoudre une équation de Vlasov 1D pour chaque faisceau de particules, et un modèle numérique moins coûteux qu'un modèle complet<br>The origin of magnetic fields observed in laboratory and astrophysical plasmas is one ofthe most challenging problems in plasma physics. In this respect, the Weibel type instabilities are considered of key importance. These instabilities are caused by a temperature anisotropy (Weibel instability) and electron momentum (current filamentation instability). The main objective of this thesis is the theoretical and numerical study of these instabilities in a collisionless plasma in the relativistic regime. The first aspect of this work is to study the nonlinear regime of these instabilities and the role of kinetic and relativistic effects on the structure of self-consistent electromagnetic fields. In this context, a key problem for the theory and applications, is the identification and analysis of coherent structures developed spontaneously in the nonlinear regime of kinetic scales. A second aspect of the work is the development of analytical and numerical techniques for the study of collisionless plasmas. A mathematical model of reference is the Vlasov-Maxwell model, where the Vlasov equation (mean field theory) is coupled to the Maxwell equations in a self-consistent way. A one-dimensional model, the multi-stream model, is also introduced. Based on a dimensional reduction technique, it is both an analytical model "simple" having the advantage of being able to solve a 1D Vlasov equation for each particle beam, and a numerical model less expensive than a complete model

Les plasmas sont des systèmes gazeux d'ions et d'électrons qui interagissent avec les champs électromagnétiques et affichent des propriétés collectives. Parmi ceux-ci, il y a la notion de "connexion" de lignes magnétiques. Ceci exprime le fait que, dans des régimes dans lesquels les particules chargées s'enroulent suffisamment vite le long des lignes d'induction magnétique, ces dernières sont liées au mouvement du plasma massif et acquièrent une identité topologique qui leur interdit de se rompre, se croiser et se reconnecter. Cette identité topologique peut cependant être localement violée grâce à un certain nombre d'effets cinétiques, comme les collisions entre les particules, lorsque les courants dans le plasma sont suffisamment intenses : on parle de “reconnexion magnétique”. La reconnexion magnétique est un ingrédient important de l'auto-organisation du plasma et a une importance pour les plasmas spatiaux et de laboratoire, car elle est à la base de phénomènes naturels comme les éruptions solaires et les aurores polaires, ou de processus disruptifs dans les expériences de fusion thermonucléaire. Un problème de longue date dans l'étude des plasmas de laboratoire et astrophysiques est de comprendre les mécanismes d'accélération des électrons et des ions, lorsqu’un champ magnétique se reconnecte et libère de l'énergie. Dans ce travail, nous avons étudié les effets cinétiques sur les instabilités de reconnexion se développant spontanément dans les nappes de courant statique (modes de déchirement) et en combinaison avec une classe d'instabilités cinétiques (instabilités de Weibel) qui sont pertinentes à la fois pour les jets de plasma astrophysiques et pour les expériences d'interaction laser-plasma. Nous avons effectué cette étude en utilisant des modèles fluides réduits et cinétiques, et nous avons étudié la concurrence entre les modes de type déchirement et les instabilités de type Weibel au moyen de simulations cinétiques complètes avec codes semi-lagrangiennes Vlasov-Maxwell et de type “Particle-In-Cell“<br>Plasmas are gaseous systems of ions and electrons which interact via electromagnetic fields and display collective properties. Among these, is the notion of the magnetic line "connection". This expresses the fact that, in regimes in which charged particles spiral sufficiently fast along lines of magnetic induction, the latter is linked to the bulk plasma motion and acquire a topological identity which forbids them to break, intersect and reconnect. This topological identity, however, can be locally violated thanks to a number of kinetic effects, such as particle collisions, when the currents in the plasma are sufficiently intense: one speaks of "magnetic reconnection". Magnetic reconnection is an important ingredient of the plasma self-organization and has significance for both space and laboratory plasmas since it is at the basis of natural phenomena like solar flares and polar lights, or of disruptive processes in thermonuclear fusion experiments. A long-standing problem in the study of laboratory and astrophysical plasmas is to understand the mechanisms of acceleration of electrons and ions, as a magnetic field reconnect and release energy. In this work, we studied kinetic effects on reconnection instabilities developing spontaneously in static current sheets (tearing modes) and in combination with a class of kinetic instabilities (Weibel instabilities) that are relevant both to astrophysical plasma jets and to laser-plasma interaction experiments. We performed this study using reduced-fluid and kinetic models and we investigated the competition between tearing-type modes and Weibel-type instabilities by means of both semi-lagrangian full kinetic Vlasov-Maxwell simulations and particles in cell simulations

Nous avons modélisé numériquement la stabilité d'un système de deux couches d'un même fluide pur supercritique soumis à une différence de température initiale à l'interface. La grande compressibilité et la faible diffusivité thermique des fluides critiques entraînent une instabilité gravitationnelle de type Rayleigh-Taylor de la couche de diffusion. Cette instabilité est similaire à celle que l'on retrouve dans le cas de deux fluides miscibles à condition de remplacer le coefficient de diffusion moléculaire par celui de la diffusion thermique. Nos résultats numériques semblent être consistants, d'un point de vue de la relation de dispersion dans le cas linéaire, avec ceux de Duff et al. [Duf62] considérés dans le cas de deux fluides miscibles. Nous avons aussi montré que, lorsque l'épaisseur de la couche inférieure devient plus petite que l'épaisseur de la couche de diffusion thermique basée sur le taux de croissance maximal alors le système devenait stable. Un diagramme de stabilité a été établi en fonction de trois paramètres: l'épaisseur de la couche inférieure, la différence de densité entre les deux couches et la distance au point critique. Lorsque l'on s'approche du point critique, la stratification devient de plus importante (à cause de la forte compressibilité) et tend à stabiliser la configuration. Par ailleurs, le filtrage acoustique initialement utilisé afin de réduire les coûts de calcul, ne s'est pas avéré nécessaire puisque la méthode des volumes finis est naturellement filtrante (formulation intégrale ou variationnelle des équations).

Les chocs non-collisionnels jouent un rôle majeur dans de nombreux événements astrophysiques à haute densité d'énergie (sursauts gamma, restes de supernovæ, vents de pulsar...), et seraient responsables de la génération de particules supra-thermiques et de radiations. Les simulations ont démontré qu'en l’absence de champs magnétiques externes, des instabilités électromagnétiques peuvent prendre place lors de la collision de plasmas à haute vitesse. Les instabilités du type Weibel sont en effet capables de faire croître, dans ces milieux, une turbulence électromagnétique potentiellement en mesure de défléchir et d'accélérer des particules par des processus du type Fermi. En plus d'une compréhension théorique toujours croissante, la génération expérimentale de tels chocs est maintenant étudiée à l'aide de lasers de puissance. Les fluctuations thermiques électromagnétiques constituent les germes des instabilités se développant dans un plasma. Nous nous sommes attelés à leur description dans le cas d’un plasma relativiste régi par une fonction de distribution de type Maxwell-Jüttner. Des formules exactes de la densité spectrale ont pu être obtenues pour différentes orientations du vecteur propre. Ces résultats ont pu être confrontés aux prédictions d’un code de simulation particle-in-cell (PIC). Un très bon accord a été démontré.Ces résultats ont été exploités lors d'une collaboration internationale dont le but était d'estimer le temps de saturation de l'instabilité cinétique de Weibel, générant des fluctuations magnétiques. Les estimations obtenues ont pu être validées par des simulations PIC sur trois ordres de grandeur d'énergie de dérive.Nous avons ensuite mené une étude théorique et numérique des collisions de plasma d'électrons-ions en régime non-collisionnel ayant lieu lors d'événements astrophysiques tels que les restes de supernovæ. Par-delà un intérêt académique pour la compréhension des processus de transfert/transport d’énergie au sein des plasmas, la récente génération de tels plasmas en laboratoire ouvre des perspectives inédites en astrophysique des hautes énergies. La zone de recouvrement de ces faisceaux de particules est sujette à des instabilités cinétiques du type Weibel, générant des champs magnétiques intenses.Nous avons modélisé l'évolution non-linéaire d'un système soumis à l'instabilité de Weibel, et obtenu des formules analytiques de l'évolution des paramètres plasmas (températures et vitesse de dérive) et des champs magnétiques. Le modèle prédit ainsi l’évolution du système jusqu’à un stade proche de l’isotropisation complète des populations de particules et donc jusqu'à la formation d’un choc non-collisionnel. Ce modèle, en accord avec des simulations du type « particle-in-cell », pu aussi être comparé à des résultats expérimentaux récents. L'étude de la propagation des chocs non-collisionnels, m'a permis de généraliser le précédent modèle au cas de la turbulence magnétique ayant lieu en amont du front de choc.Nous nous sommes consacrés enfin aux chocs non-collisionnels créés dans un plasma dense (opaque) irradié par un laser intense. L’interaction laser-plasma qui en résulte donne lieu à un important courant d'électrons relativistes qui sont à l’origine d’instabilités cinétiques (de filamentation notamment) susceptibles d'évoluer en choc non-collisionnel. Une observation originale, contrastant avec les premières publications sur le sujet est que pour les paramètres considérés (un laser d’éclairement ~1021 Wcm-2, interagissant avec une cible solide), le choc résulte de la turbulence magnétique produite par l’instabilité électronique, plutôt que par l’instabilité ionique (dont la croissance est plus tardive). En d’autres termes, compte tenu de l’énergie très élevée des électrons accélérés par le laser, la turbulence qu'ils génèrent s’avère assez forte pour rapidement défléchir les ions<br>Collisionless shocks play a major role in powerful astrophysical objects (e.g., gamma-ray bursts, supernova remnants, pulsar winds, etc.), where they are thought to be responsible for non-thermal particle acceleration and radiation. Numerical simulations have shown that, in the absence of an external magnetic field, these self-organizing structures originate from electromagnetic instabilities triggered by high-velocity colliding flows. These Weibel-like instabilities are indeed capable of producing the magnetic turbulence required for both efficient scattering and Fermi-type acceleration. Along with rapid advances in their theoretical understanding, intense effort is now underway to generate collisionless shocks in the laboratory using energetic lasers. In a first part we study the (w,k)-resolved electromagnetic thermal spectrum sustained by a drifting relativistic plasma. In particular, we obtain analytical formulae for the fluctuation spectra, the latter serving as seeds for growing magnetic modes in counterstreaming plasmas. Distinguishing between subluminal and supraluminal thermal fluctuations, we derived analytical formulae of their respective spectral contributions. Comparisons with particle-in-cell (PIC) simulations are made, showing close agreement in the subluminal regime along with some discrepancy in the supraluminal regime. Our formulae are then used to estimate the saturation time of the Weibel instability of relativistic pair plasmas. Our predictions are shown to match 2-D particle-in-cell (PIC) simulations over a three-decade range in flow energyWe then develop a predictive kinetic model of the nonlinear phase of the Weibel instability induced by two counter-streaming, symmetric and non-relativistic ion beams. This self consistent, fully analytical model allows us to follow the evolution of the beams' properties up to a stage close to complete isotropization and thus to shock formation. Its predictions are supported by 2D and 3D particle-in-cell (PIC) simulations of the ion Weibel instability in uniform geometries, as well as shock-relevant non-uniform configurations. Moreover, they are found in correct agreement with a recent laser-driven plasma collision experiment. Along with this comparison, we pinpoint the important role of electron screening on the ion-Weibel dynamics, which may affect the results of simulations with artificially high electron mass. We subsequently address the shock propagation resulting from the magnetic Weibel turbulence generated in the upstream region. Generalizing the previous symmetric-beam model to the upstream region of the shock, the role of the magnetic turbulence in the shock-front has been analytically and self-consistently characterized. Comparison with simulations validates the model. The interaction of high-energy, ultra-high intensity lasers with dense plasmas is known to produce copious amounts of suprathermal particles. Their acceleration and subsequent transport trigger a variety of Weibel-like electromagnetic instabilities, acting as additional sources of slowing down and scattering. Their understanding is important for the many applications based upon the energy deposition and/or field generation of laser-driven particles. We investigate the ability of relativistic-intensity laser pulses to induce Weibel instability-mediated shocks in overdense plasma targets, as first proposed by Fiuza in 2012. By means of both linear theory and 2D PIC simulations, we demonstrated that in contrast to the standard astrophysical scenario previously addressed, the early-time magnetic fluctuations (Weibel instability) generated by the suprathermal electrons (and not ions) are strong enough to isotropize the target ions and, therefore, induce a collisionless electromagnetic shock

Dans cette thèse de doctorat nous présentons une étude théorique de deux types d'instabilitésferroélectriques: celles apparaissant dans des géométries confinés et celles induites par le magnétismedans dans composés massifs de structure perovskite. Dans une première partie nous abordons leproblème des instabilités ferroélectriques apparaissant dans des nanotubes et des nanocoquillesoù nous développons un modèle théorique phénoménologique approprié à ces structures. Nousétudions comment l'émergence de la polarisation est affectée par (i) l'épaisseur des nanostructures,(ii) par la réponse diélectrique des matériaux environant la couche ferroélectrique et (iii) les conditionsaux interfaces. Nous observons un effet de taille finie topologique qui peut promouvoirune compétition inhabituelle entre deux types de distribution de la polarization, irrotationel eten vortex, dans la limite des très petites épaisseurs. Dans une deuxième partie nous utilisons descalculs ab-initio à base de la théorie de la fonctionnelle de la densité pour étudier les instabilitésferroélectriques des perovskites manganites à base de terres rares (RMnO3). A partir de ces calculsnous prédisons qu'il est possible d'induire une transition de phase sous pression dans EuMnO3 lefaisant transiter d'un ordre antiferromagnétique de type A isolant vers un ordre ferromagnétiquemétallique sous pression. Ce type de transition n'avait jamais été reporté précédemment dans lesmatériaux RMnO3. Nous étendons ensuite cette analyse à l'étude des effets de strain épitaxial dansles films minces de TbMnO3 et EuMnO3. Nos résultats montrent que le diagramme de phase souscontrainte d'épitaxie est bien plus riche que celui sous pression hydrostatique. Nous trouvons queles types antiferromagnétiques E-AFM et E*-AFM sont stabilisés dans le cas de TbMnO3, où letype E*-AFM est une phase métallique polaire. Dans le cas de EuMnO3, nous trouvons une phaseantiferromagnétique de type E qui n'a pas été observée sous pression hydrostatique<br>In this thesis, we present a theoretical study of two types of ferroic instabilities: the ferroelectric instability in novel confined geometries and magnetic instabilities controlled by the distortion of the underlying crystal lattice. On the one hand, we consider in detail the ferroelectric instability, specifically, in the nanotubes and the spherical nanoshells and develop a phenomenological theory for describing such an instability. We determine how the emergence of polarization is affected bythe thickness of the nanoparticle, the dielectric properties of the surrounding media and the interfacial boundary conditions. We finnd an intriguing topological finite-size effect that can promote an unexpected competition between two different types of distribution of polarization - irrotational and vortex-like - in the ultra-thin limit. One the other hand, we employ a different formalism to investigate the structural, electronic and magnetic properties of the rare-earth manganites. Specifically,we conduct a theoretical investigation from first-principles calculations. First, we predict a pressure-induced A-AFM insulator to FM metal transition on EuMnO3 under hydrostatic pressure, that is unprecedented in the multiferroic rare-earth manganites RMnO3. This investigation is extended to the study to the epitaxial strain effects on both EuMnO3 and TbMnO3 thin films. We show that epitaxial strain generates a much richer phase diagram compared to hydrostatic pressure. We predict novel magnetically-induced insulator { metal and polar { non-polar transitions. More specifically, we find that both the multiferroic E-AFM order and the polar metallic E*-AFM state are stabilized in TbMnO3 by means of epitaxial strain. In the contrast, we find a novel epitaxial-strain-induced multiferroic E-AFM state in EuMnO3 that cannot be obtained by means of just hydrostatic pressure

Etude de la transition paraélectrique-ferroélectrique des composés KH₂PO₄ au moyen d'un interféromètre à balayage qui permet l'observation du spectre vibrationnel complet. Les fréquences et amortissements des modes transverses et longitudinaux optiques sont obtenus en simulant les spectres expérimentaux au moyen de la forme factorisée de la fonction diélectrique. Par comparaison des résultats, on identifie les modes de vibration. On montre que la transition de phase est liée à la déstabilisation de l'un des modes internes du tétraèdre. La transition est displacive, mais induite par un phénomène de désordre des protons.

Les propulseurs plasmiques de type Hall possèdent des caractéristiques fonctionnelles qui répondent particulièrement bien aux besoins en propulsion des nouvelles générations de satellites. Le fonctionnement de ces propulseurs est basé sur la création et l'accélération d'un plasma de xénon dans une décharge en champs croisés E x B. Dans le cadre d'un groupement de recherche constitué du CNRS, du CNES, de la SNECMA et de l'ONERA, un effort particulier est porté pour améliorer la compréhension des mécanismes physiques qui gère le fonctionnement de ces propulseurs. L'étude expérimentale de la réponse du moteur à des extinctions de courant obtenues avec un interrupteur ultra-rapide a pris une place particulière dans cette thèse. Cette approche a permis d'apporter un nouvel éclairage sur les processus collisionnels de création du plasma et a rendu possible une caractérisation fine du comportement dynamique de la décharge. Le transport des électrons dans la décharge à travers le champ magnétique reste encore un problème ouvert. Des mesures à l'aide d'antennes large-bande ont permis d'observer la corrélation entre le développement d'instabilités haute fréquence (6 MHz) et une augmentation macroscopique du courant. Ces mesures mettent en évidence l'existence d'un transport turbulent des électrons dans la décharge.

Les satellites artificiels baignent dans un environnement radiatif hostile qui conditionne en partie leur fiabilité et leur durée de vie en opération : les ceintures de Van Allen. Afin de les protéger, il est nécessaire de caractériser la dynamique des électrons énergétiques piégés dans les ceintures radiatives. Elle est déterminée essentiellement par les interactions entre les électrons énergétiques et les ondes électromagnétiques existantes. <br /><br />Le travail de cette thèse a consisté à concevoir un schéma numérique original pour la résolution du système d'équations modélisant ces interactions : les équations de Vlasov-Maxwell relativistes. Notre choix s'est orienté vers des méthodes d'intégration directe. Nous proposons trois nouvelles méthodes spectrales pour discrétiser en impulsion les équations : une méthode de Galerkin et deux méthodes de type collocation. Ces approches sont basées sur des fonctions de Hermite qui ont la particularité de dépendre d'un facteur d'échelle permettant d'obtenir une bonne résolution en vitesse. <br /><br />Nous présentons dans ce manuscript les calculs conduisant à la discrétisation et à la résolution du problème de Vlasov-Poisson monodimensionnel ainsi que les résultats numériques obtenus. Puis nous étudions les extensions possibles des méthodes au problème complet relativiste. Afin de réduire les temps de calcul, une parallélisation et une optimisation des algorithmes ont été mises en \oe uvre. Enfin, les calculs de validation du code 1Dx-3Dv, à partir d'instabilités de types Weibel et whistlers, à une ou deux espèces d'électrons, sont détaillés.

Cette thèse s’intéresse à la dynamique des fronts d’ablation accélérés se développant dans les expériences de Fusion par Confinement Inertiel. La FCI est un procédé destiné à obtenir l’implosion d’une capsule de deutérium ou de tritium en la comprimant par irradiation de sa couche extérieure, l’ablateur. Le front d’ablation parcourant cette coquille superficielle est déformé par des perturbations hydrodynamiques pouvant compromettre l’allumage d’une réaction de fusion. Ces travaux proposent donc une étude numérique de l’apparition d’instabilités de Rayleigh-Taylor dans le processus d’ablation, étudié en deux dimensions d’espace. La résolution numérique s’appuie sur une modélisation asymptotique du phénomène dans la limite d’un grand rapport de température entre le plasma chaud vaporisé par le laser et l’ablateur encore froid. Cette modélisation est effectuée en deux étapes. Tout d’abord, la partie thermo-diffusive est approchée par un système de type Hele-Shaw. Nous considérons ensuite la partie hydrodynamique comme une perturbation du système précédent. Asymptotiquement, il s’agit de résoudre l’advection d’une nappe de vorticité localisée sur le front. La discrétisation du système limite passe par une formulation eulérienne couplée avec un suivi de front par des marqueurs. La partie thermique est résolue par la méthode de la Frontière Elargie récemment développée et la partie hydrodynamique par une méthode de Volumes Finis<br>This thesis deals with the dynamics of accelerated ablative front spreading in Inertial Confinement Fusion experiments. ICF is designed for the implosion of a mdeuterium or tritium spherical target. The outer shell, the ablator, is irradiated providing a high level pressure inside the target. During this first stage, the ablation front propagating inward is perturbed by hydrodynamics instabilities, which can prevent the fusion reaction in the decelerated stage. We propose here a study on Rayleigh-Taylor instabilities during ablation process, in the two dimensional case. In order to obtain a numerical solution, we perform an asymptotic analysis in the limit of a high temperature ratio, between the remaining cold ablator and the hot ablated plasma. This study is divided in two steps. First, the thermo-diffusive part of the set of equations is approximated by a Hele-Shaw model, which is then perturbed by the hydrodynamics part. Using a vortex method, we have to solve the advection of a vortical sheet moving with the ablation front. We compute the numerical solution on an eulerian mesh coupled with a front tracking method. The thermal part is calculated by implementing the Fat Boundary Method, recently developped. The hydrodynamic part is obtained from a Finit Volume scheme

Notre etude a pour objectif de mieux caracteriser l'ecoulement dans l'entrefer des moteurs electriques rapides de type ouvert. Ceci afin d'ameliorer le refroidissement de ces moteurs en optimisant les echanges thermiques par convection forcee. La connaissance de donnees physiques (conditions aux limites et parametres a priori mal connus), obtenues experimentalement est indispensable au developpement et la mise en uvre de calculs permettant la modelisation de l'ecoulement dans cette geometrie particuliere (entrefer encoche). Nos essais ont consiste en des visualisations et des mesures de vitesses. La methode de mesure utilise l'anemometrie laser tandis que les visualisations s'effectuent par tomographie laser ou par eclairage classique au moyen de spot halogene puissant. Deux series d'essais ont ete effectuees: l'une correspondant a un ecoulement de rotation sans debit axial, permet de faire apparaitre l'influence de la geometrie sur la structure de l'ecoulement. En effet, on s'est attache a indiquer les modifications resultant de la presence des encoches par rapport au probleme analogue de couette-taylor entre cylindres lisses. L'autre avec rotation et debit axial a pour but d'obtenir en premieres approches les principales caracteristiques de l'ecoulement dans des conditions qui se rapprochent de celles d'un moteur reel de type ouvert. Les evolutions des grandeurs instationnaires (vitesses moyennes, fluctuations moyennes) ont ete determinees en diverses sections puis interpretees. Une estimation du debit d'encoche a ete effectuee afin de permettre d'evaluer le role de l'encoche dans le transfert de masse dont dependent les echanges convectifs. Une confrontation de nos resultats avec ceux de la modelisation a permis d'estimer et d'orienter le code de calcul mis au point pour etudier l'ecoulement et les transferts thermiques dans cette generation de moteurs

La thèse est consacrée aux instabilités de type Rayleigh-Taylor. Elle est inspirée par les progrès récents dans la physique de la compression des plasmas, surtout avec des dispositifs comme les Z-pinches. L'instabilité Rayleigh-Taylor (IRT) joue un rôle important dans l'évolution des plasmas magnétisés dans les expériences, aussi bien que dans les étoiles et dans les fluides classiques. Dans les phénomènes plasmas associés à la fusion thermonucléaire, l'IRT est souvent le phénomène qui limite les possibilités de compression. On essaie d'étudier des particularités de ces instabilités afin d'en profiter et d'éliminer leurs effets nuisibles. Dans ce travail on étudie aussi bien le cas général de l'instabilité Rayleigh-Taylor " classique " (dans les fluides incompressibles) que les cas plus spécifiques des instabilités de type Rayleigh-Taylor dans les plasmas magnétisés, dans les implosions des coquilles à fils etc.<br /> On a étudié l'influence de la diffusion Hall du champ magnétique sur le taux de croissance de l'instabilité. On a obtenu des solutions auto-similaires pour l'élargissement du profil initial et pour l'onde de pénétration du champ magnétique. Puis, on a étudié l'évolution postérieure du champ magnétique dans les commutateurs à ouverture de plasmas (COP) et on a montré la possibilité d'existence d'une onde de raréfaction forte, dans les cas avec et sans collisions. Cette onde peut expliquer le phénomène d'ouverture de COP.<br />Le phénomène de suppression de l'IRT par des oscillations forcées de la frontière entre deux fluides permet de proposer quelques idées pour les expériences de fusion inertielle. On considère le cas général de l'instabilité, c'est-à-dire deux fluides incompressibles visqueux dans un champ de gravitation. On obtient l'expression précise analytique du taux de croissance et on analyse l'influence des paramètres de " pompage " du système sur l'instabilité. Les résultats peuvent être appliqués à une grande diversité de situations, à partir de l'hydrodynamique classique, jusqu'aux plasmas astrophysiques. <br />Le schéma des coquilles à fils est récemment devenu une méthode très populaire d'obtention de rayonnement haute puissance ou d'implosion de haute qualité dans le domaine des Z-pinches (plasmas de striction magnétique). Il se trouve que les résultats expérimentaux sont bien meilleurs dans le cas d'une implosion de fils fins multiples situés sur un cylindre que dans le cas du schéma ordinaire du " liner ". On a examiné un problème modélisant la stabilisation de IRT dans ce cas due aux modulations régulières de la surface de séparation entre le plasma et le champ magnétique. Ce type de modulation peut apparaître grâce aux explosions initiales et à l'évolution postérieure des fils solides. On a montré que le couplage des modes de l'instabilité en présence du champ magnétique peut effectivement réduire le taux de croissance.

Dans certains matériaux contenant des atomes de solutés diffusant vers les dislocations, on observe une déformation macroscopique hétérogène lors de sollicitations de traction pour certaines vitesses de chargement et températures d'essais, cette hétérogénéité macroscopique de la déformation est appelée « phénomène Portevin - Le Chatelier » ; elle se traduit par la présence de décrochements ou de paliers sur les courbes de traction. Chaque décrochement ou palier correspond à la propagation de bandes de déformation localisée le long de l'éprouvette de traction. L'interprétation microscopique de l'effet PLC est souvent associée à une interaction dynamique entre les nuages d'impuretés ou de solutés et les dislocations mobiles. Ces instabilités qui se développent lors de procédés de formage conduisent à une baisse de ductilité du matériau ainsi qu'à un mauvais état de surface du produit fini. Le présent travail consiste à étudier l'influence de la température et de la vitesse de chargement sur les instabilités plastiques de type Portevin-Le Chatelier dans un alliage d'aluminium de la série 5000 dans le cas d'essais de traction uniaxiale effectués à l'aide de machines dure et souple respectivement à vitesses de déformation et de contrainte imposées. Le domaine d'existence de l'effet PLC dans un diagramme log [varepsilon]. - 1/T, puis log [sigma]. - 1/T est ferme aux hautes températures et vitesses de chargement élevées. Ce résultat, prévu par des approches théoriques, est confirmé par les résultats expérimentaux de la présente étude. Les mesures expérimentales de la déformation critique [varepsilon]c à laquelle apparait l'effet PLC montrent que celle-ci présente un comportement normal aux vitesses de chargement élevées. Cependant, un comportement inverse est observe aux faibles vitesses de chargement. Les chutes de contrainte [delta][sigma] et le temps de rechargement tR augment en fonction de la déformation et tendent en général vers une valeur maximale de saturation. D'autre part, la moyenne des valeurs de [delta][sigma] à la saturation présente une décroissance lorsque la vitesse de déformation augmente avec deux valeurs de pentes assez proches de celles prévues par la théorie de la diffusion. Les résultats obtenus à partir d'essais réalisés sur machine molle nous ont permis de bien étudier les différentes caractéristiques des bandes de déformation localisée et de discuter leur évolution au cours de la déformation en fonction de la température et de la vitesse de chargement. Les variations des déformations [delta][varepsilon]b contenues dans les bandes de déformation présentent une similitude avec celles de la déformation critique [varepsilon]c : on observe pour [delta][varepsilon]b la même transition entre comportements inverse et normal. Les vitesses de propagation des bandes Vb, qui présentent les mêmes variations avec la déformation que les largeurs de bandes Wb, diminuent lorsque la vitesse de contrainte [sigma]. Augmente. Certains résultats expérimentaux ont été comparés à ceux prévus par les modèles théoriques de Hähner, Kubin-Estrin et Bréchet-Estrin<br>In certain material containing solutes diffusing towards dislocations, one observes a heterogeneous macroscopic deformation during tensile tests for certain loading rates and temperatures, this macroscopic heterogeneity of the deformation is called “Portevin – Le Chatelier” effect ; it results in the presence of serrations or stairs on the traction diagrams. Each serration of stair corresponds to the propagation of deformation bands located along the tensile specimen. The microscopic interpretation of PLC effect is often associated with a dynamic interaction between the clouds of impurities or solute atoms and mobile dislocations. These instabilities which develop during forming process lead to a fall of ductility of material as well as in a bad surface quality of the finished product. This work consists in studying the influence of the temperature and loading rate on plastic instabilities of type Portevin-Le Chatelier in an aluminium alloy of 5000 series in the case of uniaxial tensile tests carried out using hard and soft tensile machines respectively at imposed strain and stress rates. The domain of occurrence of PLC effect in log [varepsilon]. - 1/T, then log [sigma]. - 1/T diagram is closed at high temperatures and loading rates. This result, envisaged by theoretical approaches, is confirmed by the experimental results of this stuudy. This result, envisaged by theoretical approaches, is confirmed by the experimental results of this study. Experimental measurements of the critical strain [varepsilon]c for the oneset of repeated yielding show that [varepsilon]c presents a normal behaviour for high loading rates. However, an inverse behavior is observed at low loading rates. The stress drop [delta][sigma] and reloading times tR increase with strain and generally tend to a maximum value of saturation. Further, the average value of saturation presents a decrease when strain rate increases with two values of slopes rather close to those envisaged by the theory of the diffusion. The results obtained starting from tests carried out on soft machine enabled us to investigate different characteristics of localised deformation bands and to discuss their evolution during the deformation as function of stres rate and variopus test temperatures. The variations of band strain [delta][varepsilon]b present a similarity with those of the critical strain [varepsilon]c : one observes for [delta][varepsilon]b the same transition between inverse and normal behaviours. Band propagation velocities Vb, which present the same variations with strain that the bandwidths Wb, decrease when the stress rate [sigma]. Increases. Certain experimental results were compared with those envisaged by the theoretical models of Hähner, Kubin-Estrin and Brechet-Estrin

L'extrusion est un procédé de mise en forme très répandu dans des industries variées telles que la plasturgie, la métallurgie, les céramiques, ou encore l'industrie agroalimentaire. Les matériaux doivent cependant être soumis à des niveaux de pression et de contrainte élevés pour être extrudés à des débits suffisants.<br />La présente thèse a pour objectif de mieux connaître l'influence des paramètres de pression et contrainte sur les propriétés des polymères et la stabilité de leurs écoulements, et ce en relation avec leur structure mésoscopique.<br />Les effets de la pression ont été étudiés sur quatre polyéthylènes (PE) de structures différentes. D'abord, les conditions expérimentales précises ont été définies pour isoler les effets de la pression des effets de la température et pour différencier les écoulements stables ou instables. Nos résultats montrent que les effets de la pression sont les mêmes en cisaillement et en élongation, sauf pour l'un des PE où ils sont 30% supérieurs en cisaillement. En outre, c'est un critère de contrainte critique qui caractérise l'apparitions des instabilités viscoélastiques quelle que soit la pression moyenne.<br />Dans la deuxième partie de cette thèse, trois copolymères de la famille des SEBS montrant une séparation de phase aux échelles nanoscopiques ont été considérés. L'étude des défauts macroscopiques d'extrusion, et de la propagation de fissures surfaciques, de ces copolymères à blocs a permis d'identifier l'origine du défaut dit de « refente d'extrudat » ainsi qu'un nouveau régime d'extrusion : le « pelage continu ». Ces défauts ont été mis en relation avec la structure mésoscopique par des essais de diffusion de rayons-X aux petits angles.

Aujourd'hui traverser les océans peut s'effectuer sans difficulté en utilisant toute sorte de véhicule marin. Cependant, il n'en est pas de même pour l'exploration des fonds marins. Face à ce milieu hostile et dangereux, mais potentiellement riche tant sur le plan biologique que sur le plan d'éventuelles ressources exploitables, l'homme a besoin d'assistance dans sa découverte des profondeurs, l'intervention de robots sous-marins est une alternative pour écarter les dangers qui peuvent subir les plongeurs. Dans ce contexte, le présent manuscrit porte sur la modélisation et la commande d'un ROV (Remotely Operated Vehicle) destiné pour l'observation des sites archéologies. Après avoir identifié les différentes variables caractérisantes de la géométrie fixe, nous détaillons, dans un premier temps, la modélisation cinématique et dynamique du véhicule. Il est important de noter qu'un jeux de paramètres d'inertie, masse ajoutée et coefficients de traînées est identifié moyennant les caractéristiques géométriques du robot. Un modèle hydrodynamique non-linéaire et complet du sous-marin a ainsi pu être développé. Pour réussir une opération d'observation, en utilisant les caméras embarquées, la deuxième partie de la thèse traite le problème de stabilisation du modèle du ROV à l'équilibre. Nous proposons une commande instationnaire explicite dépendante à la fois de l'état et du temps, suivie d'une étude de robustesse de la commande par rapport aux perturbations extérieures vérifiant certains degrés d'homogénéité. L'autonomie d'une opération d'observation nécessite aussi de contrôler les déplacements du Rov tout au long d'une trajectoire de référence. Nous avons traité dans la troisième partie de ce manuscrit le problème de stabilité et de stabilisation d'un système d'erreur entre la position réelle du véhicule et la position d'un modèle de référence. Afin d'exploités les résultats théoriques de la thèse, en bénéficiant de l'expertise de l'équipe IRA2 dans le domaine de la Réalité Virtuelle (RV), en dépit des simulations classiques réalisées sous Matlab, nous proposons de construire un environnement sous-marin (ou piscine) qui intègre la CAD du Rov: simulation interactive pour la navigation. Ainsi, on a posé les problématiques liées aux capteurs virtuels et la construction des observateurs, interfaçage des boucles de commande (à travers simulink) et la plateforme virtuelle qui fait appel à Virtools<br>Today cross oceans can be done easily. However, it is not the same case for the seabed exploration. As this hostile and dangerous environment can be biologically rich and has exploitable resources, the man needs help in his discovery of the depths. Therefore, the intervention of underwater robots was a solution. In this context, the present manuscript deals with modeling and control of a ROV (Remotely Operated Vehicle). After identifying the different variables characterizing the fixed geometry, we study, at first, kinematics and dynamic modeling of the ROV. It is important to note that sets of inertia parameters of added mass and streaks coefficients is identified by means of the geometrical characteristics of the robot. A full nonlinear dynamic model of the submarine has been established. The second part of the thesis deals with the stabilization problem of the ROV’s model. We offer an explicit unsteady dependent control of both the state and time. A robust study of the control relative to external interference checking certain degree of homogeneity has been established. The autonomy of Rov also requires control movement along a reference path. We treated in the third part the Rov’s stability problem to ensure the tracking of a reference trajectory. These results are operated on a virtual platform, and implemented on the dedicated Virtools software for this application. To lighten the structure in terms of sensors and because of the high prices of various sensors, it is necessary to design a system called auxiliary observer who charge rebuild unmeasurable states using available information. A nonlinear observer has been proposed to the estimation of linear and non-measurable angular velocity, which will be considered as virtual sensors. These sensors will be implemented on the platform that will be used to animate the ROV in its virtual world

Le transport de grains cisaillés par un fluide en écoulement est fréquemment rencontré dans la nature et dans l'industrie. Sur un sol non érodable, si le débit de grains transportés par charriage est assez faible et l'écoulement du fluide est dans un sens unique, le lit n'est pas continu, mais constitué de dunes isolées en forme de croissant: les barkhanes. On présente ici une étude expérimentale de la dynamique de dunes du type barkhanes sous un écoulement turbulent d'eau, dans une conduite de section rectangulaire. D'une part, on s'intéresse à la morphologie et à la migration des barkhanes, d'autre part nous mesurons la perturbation causée par la dune et identifions le frottement qu'il exerce sur sa surface<br>The transport of solid particles entrained by a fluid flow as bed-load is frequently found in nature and in industrial environments. If bed-load takes place over a non-erodible ground, and if the particles flow rate is small enough and the fluid flow one-directional, an initial thin continuous layer of particles becomes discontinuous and composed of isolated dunes with a crescentic shape: the barchans. We present here an experimental study concerning barchans dynamics when submitted to a rectangular closed-conduit turbulent fluid flow. On the one hand, we investigate the migration of barchan dunes. On the other hand, we investigate the fluid flow perturbation caused by the dune, in particular the shear stress acting on its surface

On s'est intéressé dans ce travail à la résolution numérique des équations d'Euler linéarisées autour d'un écoulement stationnaire, subsonique et assez régulier. Dans le but d'obtenir des matrices symétriques dans ces équations, et in fine, une équation d'équilibre énergétique, nous considérons la linéarisation d'une forme symétrique des équations d'Euler tridimensionnelles. Nous proposons un schéma non-diffusif de type Galerkin discontinu en domaine temporel (GDDT) s'appuyant sur une formulation centrée-élément avec des ux numériques totalement centrés et un schéma en temps explicite de type saute-mouton, ce qui permet d'obtenir une approximation sans dissipation et fournit une estimation précise des variations de l'énergie aéroacoustique. En effet, dans le cas général de la linéarisation autour d'un écoulement non-uniforme, il existe une équation d'équilibre énergétique au niveau continu que nous vérifions au niveau discret. Nous montrons qu'il existe un terme source discret permettant de conserver exactement l'énergie ce qui permet de prouver la stabilité de notre schéma. Ainsi, notre schéma non-diffusif de type GDDT fournit un outil précis pour contrôler des phénomènes tel que les instabilités de Kelvin-Helmholtz. Nous illustrons la capacité de notre méthode aussi bien sur plusieurs cas tests académiques que sur différentes configurations complexes grâce à une implémentation parallèle.

Cette thèse s'inscrit dans le cadre de l'astrophysique de laboratoire. Nous abordons divers aspects de la physique des chocs non-collisionels en présence de flots de plasma relativistes dans des configurations d'intérêt pour les communautés astrophysique et de l’interaction laser-plasma (ILP). Notre approche repose sur la modélisation analytique et la simulation cinétique haute-performance, outil central pour décrire les processus d'ILP et la physique non linéaire à l'origine des chocs étudiés. Le code Particle-in-Cell SMILEI a été largement utilisé et développé au cours ce travail. Trois configurations physiques sont étudiées. L’instabilité Weibel en présence de faisceaux d'électrons contre-propagatifs alignés avec un champ magnétique externe est décrite. Les phases linéaires et non linéaires sont expliquées à l’aide de modèles théoriques confirmés par des simulations. La génération de chocs non-collisionels lors de l’interaction de deux plasmas relativistes de paires est étudiée en présence d’un champ magnétique perpendiculaire. L’accent est mis sur la comparaison des prédictions théoriques sur les grandeurs macroscopiques avec les simulations, ainsi que sur la définition du temps de formation du choc, l’ensemble de ces grandeurs étant d’une grande importance pour de futures expériences. Enfin, nous proposons un schéma permettant de recréer, en laboratoire, l’instabilité Weibel ionique par l'utilisation d'un laser intense. Les flots de plasmas produits ici sont plus rapides et denses que dans les expériences actuelles, conduisant à un taux de croissance et des champs magnétiques plus élevés. Ces résultats sont également important pour l’ILP à très haute intensité<br>The work presented in this thesis belongs to the general framework of Laboratory Astrophysics. We address various aspects of the physics of collisionless shocks developing in the presence of relativistic plasma flows, in configurations of interest for the astrophysical and the laser-plasma interaction (LPI) communities. The approach used throughout this thesis relied on both analytical modeling and high-performance kinetic simulations, a central tool to describe LPI processes as well as the non-linear physics behind shock formation. The PIC code SMILEI has been widely used and developed during this work. Three physical configurations are studied. First we consider the Weibel instability driven by two counter-streaming electron beams aligned with an external magnetic field. The linear and non-linear phases are explained using theoretical models confirmed by simulations.Then the generation of non-collisional shocks during the interaction of two relativistic plasma pairs is studied in the presence of a perpendicular magnetic field. We focus on the comparison of theoretical predictions for macroscopic variables with the simulation results, as well as on the definition and measurement of the shock formation time, all of which are of great importance for future experiments.Finally, we proposed a scheme to produce, in the laboratory, the ion-Weibel-instability with the use of an ultra-high-intensity laser. The produced flows are faster and denser than in current experiments, leading to a larger growth rate and stronger magnetic fields. These results are important for the LPI at very high intensity

Le couplage entre la diffusion microscopique et la circulation méridienne dans l'intérieur radiatif des étoiles en rotation est encore mal connu. Dans cette thèse, nous nous intéressons plus particulièrement aux effets sur la circulation méridienne du gradient de poids moléculaire moyen dû à la diffusion microscopique de l'hélium. Pendant longtemps ces effets ont été négligés mais on sait aujourd'hui que leurs conséquences sont importantes et peuvent être directement comparées aux observations d'abondances chimiques et d'oscillations stellaires. Au cours de cette étude, nous avons mis l'accent sur différents points : 1) Compréhension et modélisation du couplage diffusion-circulation : simulation 2D du couplage pour une meilleure compréhension des processus physiques mis en jeu, modélisation du couplage à une dimension pour les calculs d'évolution stellaire. 2) Introduction dans les calculs d'évolution stellaire de la modélisation 1D et application à différents types d'étoiles (étoiles de pop. I et II). .<br>The coupling between microscopic diffusion and meridional circulation is not yet well-known. In this PhD, we are interested in the effects of molecular weight gradients induced by microscopic diffusion on the meridional circulation. These effects had been neglected for a long time but we now know that they can be important and compared to the observations of chemical abundances and stellar oscillations. During the PhD we concentrated on several points : 1) Understanding and modeling the diffusion-circulation coupling using a 2D simulation. 2) Introduction of a deduced 1D modeling in the Toulouse-Geneva stellar evolution code and applications to different types of stars (pop. .

Avec les nouveaux enjeux en matière de pollution et de gain de puissance des turbomachines modernes, il est aujourd’hui nécessaire de maîtriser de façon optimale la combustion du carburant. Les nouvelles performances recherchées imposent donc aux motoristes un cahier des charges strict sur les systèmes d'injection. C’est dans le but d’optimiser le développement des systèmes d’injection à effet de pression de type tourbillonnaire double-débit utilisés dans les chambres de moteurs aéronautiques qu’un banc expérimental reproduisant, à grande échelle, les différents circuits constitutifs d’un tel injecteur a été mis en place à l'ONERA Centre de Toulouse. Une campagne d’essais a été effectuée afin de fournir une banque de données permettant une meilleure compréhension de la structure de l'écoulement interne d’un injecteur tourbillonnaire et de son influence sur l'atomisation du jet en sortie d’orifice. Les maquettes développées sont réalisées entièrement en matériau transparent, et sont constituées de différents modules interchangeables afin d’explorer plusieurs géométries. Des campagnes de mesure par fluorescence laser, par LDA et par PIV ont foumi une cartographie de l'écoulement interne révélant une architecture complexe fortement instationnaire. Schématiquement, le liquide à pulvériser entre dans une chambre dite à swirl via plusieurs fentes placées tangentiellement. L'application d’un important mouvement tourbillonnaire induit une structure de l'écoulement de type vortex de Rankine et se caractérise aussi par la formation d’une colonne d’air tout le long de l’axe. Par ailleurs, les instabilités de ce noyau d’air central ont pu être caractérisées : le caractère instationnaire des jets d'alimentation a une importance essentielle sur les origines des instabilités se propageant dans l'injecteur. En parallèle, ces données ont été exploitées avec l'objectif de valider un modèle 2D axisymétrique et un modèle 3D réalisés avec les logiciels Fluent 5. 7 et 6 suivant une méthode diphasique VOF d’ordre 2 avec reconstruction d'interface.

Cette thèse a pour objectif de détecter les structures turbulentes aux grandes échelles présentes dans une couche de mélange de type Rayleigh-Taylor incompressible à faible nombre d'Atwood. Diverses grandeurs statistiques conditionnées par la présence de ces structures ont été obtenues, et il est désormais possible de les comparer avec les résultats des modèles statistiques turbulents dits bi-structure, tel le modèle 2SFK développé au CEA. Afin de réaliser les simulations numériques directes du mélange turbulent, un code numérique tridimensionnel incompressible à densité variable a été développé. Ce code a été parallélisé dans les trois directions. Plusieurs méthodes de détection de structure ont été conçues et testées. Bien que toutes ces méthodes présentent différents intérêts, seule la plus efficace vis-à-vis de nos critères de détection a été gardée pour faire des simulations à forte résolution (plus d'un milliard de mailles, 1024^3). Un filtrage temporel de la vitesse verticale est utilisé dans cette méthode de détection afin de : 1) corriger les distorsions dues aux points d'arrêt et zones de recirculation dans l'écoulement, 2) minimiser l'effet de la turbulence aux petites échelles et mieux mettre en évidence les grandes échelles, 3) introduire un effet mémoire permettant de prolonger la bimodalité du champ de détection depuis les zones laminaires extérieures jusqu'au centre de la zone de mélange turbulent. Plusieurs simulations numériques directes 1024^3 ont été effectuées. Les résultats viennent conforter ceux obtenus avec le modèle bi-structure 2SFK et justifient une étude plus poussée des grandeurs statistiques en vue de sa validation.

Les travaux de recherche présentés dans cette HDR sont répartis en trois thématiques complémentaires. La première concerne la modélisation des instabilités locales dans les matériaux métalliques (localisation, striction ...). Il s'agit de développer des outils théoriques et numériques de prédiction de ces phénomènes et leur validation au travers de courbes limites de formage pour des aciers ferritiques ou dual-phase. Deux approches complémentaires ont été adoptées : la première est micromécanique se basant sur la plasticité cristalline couplée à des techniques de transition d'échelles, tandis que la seconde est phénoménologique utilisant des modèles de comportement avancés que nous avons couplé à de l'endommagement. Le critère de bifurcation de Rice, qui est relié à la perte d'ellipticité des équations gouvernant le problème aux limites, a été particulièrement analysé au travers de ces deux approches. L'influence sur l'apparition d'une localisation de plusieurs paramètres matériau ou liés à la microstructure ainsi que l'effet de certains mécanismes déstabilisants ont été analysés. Il est mis en évidence, notamment, le rôle tout à fait essentiel de la formation de points de vertex aux points courants de la surface de plasticité dans la détection d'une localisation de la déformation plastique en bande de cisaillement. Le deuxième axe de recherche est relatif aux instabilités de type structure (flambage, plissement ...). Ces phénomènes se manifestent le plus souvent en présence de structures minces ou relativement élancées. Dans ce cadre, nous proposons une approche basée sur l'étude de stabilité des évolutions quasi statiques. Pour des matériaux sensibles à la vitesse de déformation, l'absence de bifurcation ainsi que d'états d'équilibre nous amènent naturellement à poser le problème comme celui de la stabilité de trajectoires quasi statiques. Pour les matériaux élasto-plastiques, cette approche est justifiée par le fait que bien souvent nous sommes en présence d'une évolution quasi statique, pour un trajet de chargement donné, même si chaque point de cette évolution représente un état d'équilibre. Le critère de stabilité unifié que nous proposons est donné par la positivité de la seconde variation de l'énergie totale et est valable pour des solides visco-élastiques, visco-plastiques ou élasto-plastiques. Plus récemment, nous avons étendu ce critère à des modèles à gradients. Enfin, la comparaison de ce critère avec les résultats existants relatifs au flambage plastique nous a permis de rediscuter le critère de non-bifurcation de Hill en relation avec le choix du modèle de plasticité approprié. Il est mis en évidence, de nouveau, le rôle important joué par le développement de points de vertex sur la surface de plasticité lors de la prédiction du flambage avec la théorie du module tangent élasto-plastique. Le troisième axe de recherche concerne le développement d'éléments finis de type coques volumiques. L'idée de ce nouveau concept est de combiner les avantages des éléments coques et 3D pour formuler des éléments particulièrement adaptés à la simulation des structures minces. Ces nouveaux éléments de coques volumiques ont l'avantage de prendre en compte la flexion suivant une direction choisie tout en conservant la formulation classique des éléments volumiques. Ils permettent ainsi de modéliser des structures comportant des parties minces et des zones épaisses, sans les difficultés habituelles de raccord de maillages contenant des éléments coques et 3D. Leur utilisation en emboutissage est très prometteuse, en particulier pour des problèmes où les effets dans l'épaisseur sont d'importance majeure. Étant sous-intégrés, pour améliorer leurs performances, les modes de hourglass générés sont alors efficacement stabilisés par des techniques récentes. Les modes de verrouillages (membrane et cisaillement transverse) sont éliminés par des techniques de projection pouvant se mettre sous le formalisme " Assumed Strain Method ".

Les interactions hydrodynamiques sont souvent négligées dans le transport de particules à cause de leur complexité. Aux fortes concentrations, ces interactions se manifestent par une différence entre les cinématiques de ces particules et celle du fluide non perturbé contrairement aux hypothèses communément admises. Aujourd'hui, seuls de lourds calculs numériques (DNS) peuvent résoudre ces problèmes. Nous développons ici une méthode alternative, basée sur la matrice de résistance généralisée, appliquée au transport d'une particule cylindrique confinée. Pour la sédimentation, elle chute plus vite en position dissymétrique. En déterminant la cinématique de la particule en écoulement de Poiseuille nous montrons la non validité de l'hypothèse évoquée. Ensuite nous étudions les phénomènes de transferts sur le même cylindre. Selon la condition aux limites, on peut soit réduire, soit augmenter le flux sur le cylindre.

碩士<br>國立虎尾科技大學<br>光電與材料科技研究所<br>99<br>Even though the field of EM wave amplification for microwaves has been rather mature and enjoying many available technologies, amplifying a lightwave already left the laser resonant cavity can be a tough problem with very little choice. Long before erbium-doped optic fiber amplifier (EDFA) has become a reality and mass produced, other lightwave amplification schemes have been considered, even realized, and tested. Among them, the most noticed are the Stimulated Raman Scattering (SRS) and the Stimulated Brillouin Scattering (SBS). Comparing with the EDFA, a major difference in them both is that in their lightwave amplification mechanisms, no atomic (molecular) energy levels are involved, and only pure plasma physics processes are relevant. However, up to this day, it has been concluded that both the SBS and SRS amplification effects are so much weaker than the gain that stimulated emission provides in a doped-fiber amplifier that Raman and Brillouin amplifiers tend to involve very long distances and very high pump powers. On the other hand, since the EDFA amplification approach is entirely limited by quantum transition physics among energy levels, even though it may at best provide 980 nm and 1490 nm spectrum intervals for communication purposes, overall, it fails to give a real wide-band, flat-top working spectrum. Furthermore, apparently EDFA can only amplify lightwaves guided by optic fibers, not those flying in the open space and of arbitrary wavelengths and amplitudes. Nevertheless, mankind in fact has been longing for a lightwave-enhancing technology that is essentially unrestricted by available atomic energy levels, and at the same time can be applied on lights propagating in free space or within optic fibers. In man-made plasma sources, and in fusion experimental machines (such as the Tokamak), we often witness the working of a plasma instability called Weibel EM instability. However, in these cases, Weibel instability is an undesirable path for system energy loss. In natural environments, Weibel instability also plays a major role in causing gamma ray bursts observed by satellites. Here, we intend to direct the application of Weibel instability to a new direction, viz., amplifying essentially lightwaves of arbitrary wavelengths and amplitudes in either open or fiber-confined space, in a fashion more like the aforementioned SRS and SBS. The approach adopted is using controlled “plasma”, in the form of vertically (with respect to the incident light) oscillating electrons, to trigger the Weibel instability to further cause exponential growth of the incident lightwave amplitude. This current research mainly aims to test the feasibility of such lightwave Weibel amplification theory as a preliminary step toward the ultimate goal of air-borne lightwave amplification.