L’objectif de ce travail est de proposer un schéma numérique pertinent pour les équations de Saint-Venant avec termes source de topographie et de friction de Manning. Le premier chapitre est dédié à l’étude du système de Saint- Venant muni des termes source. Dans un premier temps, les propriétés algébriques de ce système sont obtenues. Dans un second temps, nous nous intéressons à ses états stationnaires, qui sont étudiés pour les termes source individuels de topographie et de friction. Le deuxième chapitre permet de rappeler des notions concernant la méthode des volumes finis. Nous évoquons des schémas aux volumes finis pour des systèmes de lois de conservation unidimensionnels et bidimensionnels, et nous en proposons une extension permettant d’assurer un ordre élevé de précision. Le troisième chapitre concerne la dérivation d’un schéma numérique pour les équations de Saint-Venant avec topographie et friction. Ce schéma permet : de préserver tous les états stationnaires ; de préserver la positivité de la hauteur d’eau ; d’approcher les transitions entre zones mouillées et zones sèches, et ce même en présence de friction. Des cas-tests mettant en lumière les propriétés du schéma sont présentés. Le quatrième chapitre permet d’étendre le schéma proposé précédemment, pour prendre en compte des géométries bidimensionnelles et pour assurer un ordre élevé de précision. Des cas-tests numériques sont aussi présentés, y compris des simulations de phénomènes réels. / This manuscript is devoted to a relevant numerical approximation of the shallow-water equations with the source terms of topography and Manning friction. The first chapter concerns the study of the shallow-water equations, equipped with the aforementioned source terms. Algebraic properties of this system are first obtained. Then, we focus on its steady state solutions for the individual source terms of topography and friction. The second chapter introduces the finite volume method, which is used throughout the manuscript. One-dimensional and two-dimensional systems of conservation laws are studied, and a high-order strategy is presented. The third chapter deals with the numerical approximation of the shallow-water equations with topography and friction. We derive a scheme that: • preserves all the steady states; • preserves the non-negativity of the water height; • is able to deal with transitions between wet and dry areas. Relevant numerical experiments are presented to exhibit these properties. The fourth chapter is dedicated to extensions of the scheme developed in the third chapter. Namely, the scheme is extended to two space dimensions, and we suggest a highorder extension. Numerical experiments are once again provided, including real-world simulations.

La dissipation d'énergie turbulente est un processus clef dans le milieu interstellaire (MIS) froid, non seulement pour comprendre les voies de formation des étoiles, mais aussi en tant que source d'énergie supra-thermique et donc capable d'ouvrir de nouveaux chemins chimiques pour expliquer les abondances d'espèces soumises à des barrières endoénergétiques, telles que $\mathrm{CH}^+$ et $\mathrm{SH}^+$ qui sont observées dans le MIS. Dans ce contexte, l'intermittence spatio-temporelle du taux de dissipation d'énergie joue un rôle crucial car elle conduit à une injection d'énergie supra-thermique qui peut localement être très intense par rapport au taux moyen. Les caractéristiques détaillées de la distribution spatiale et les propriétés géométriques des lieux de dissipation intense peut fournir des indications précieuses pour les modèles chimiques.Nous étudions ici ces structures à l'aide de simulations numériques directes, avec un soin sans précédent pour résoudre les processus dissipatifs numériquement. Le nombre de Mach dans le MIS diffus peut prendre des valeurs aussi bien grandes que petites et nous encadrons les effets de la compression entre deux catégories de turbulence en déclin. Dans le cas extrême incompressible, nous réalisons des simulations pseudo-spectrales de magnétohydrodynamique visqueuse et résistive, avec un éclairage particulier sur la diffusion ambipôlaire due à la vitesse relative qui existe entre les ions et les neutres. Du côté de la compressibilité extrême, nous considérons des simulations isothermes basées sur des schémas sur grille (type Godunov) qui incluent dissipation visqueuse et résistive: nous nous intéressons ici particulièrement à la dissipation numérique.Nos simulations incompressibles montrent que la diffusion ambipôlaire conduit le champ magnétique à petite échelle dans une configuration libre de force de Lorentz. En conséquence, l'échelle caractéristique du chauffage par friction ion-neutre se déplace à plus grande échelle, aux échelles inertielles de la turbulence, bien plus grandes que la taille estimée par le raisonnement dimensionnel classique. Les structures dissipatives dans nos simulations sont des feuilles cohérentes spatiallement, chacune portant une nature dissipative bien distincte (visqueuse, ohmique ou bien ambipôlaire). Nous révélons les lois statistiques qui gouvernent leurs caractéristiques et nous calculons les exposants des fonctions de structure qui quantifient l'intermittence.Nous montrons que les simulations compressibles sont sujettes à une forte dissipation numérique: dans nos expériences, à peu près la moitié de la dissipation totale est attribuée aux termes de dissipation physique, le reste est produit par le schéma numérique. Nous avons mis au point une méthode pour estimer localement l'énergie perdue dans le schéma et nous l'utilisons pour examiner la structure en feuilles du champ de dissipation comme dans les simulations incompressibles. Par exemple, bien que nous confirmons que la dissipation physique visqueuse est dominée par les feuilles de cisaillement plutôt que par les chocs, nous ne pouvons exclure que la dissipation numérique ne renverse cet équilibre si celle-ci était dominée par les chocs.Pour finir, nous examinons l'efficacité de diagnostiques observationnels variés pour tracer les structures de forte dissipation. En particulier, nous trouvons que les incréments de centroïdes de vitesse ou des paramètres de Stokes corrèle très bien sur le plan du ciel avec certaines structures de forte dissipation. Nous calculons aussi les exposants de l'intermittence mesurés pour ces mêmes traceurs et nous trouvons qu'ils s'étendent sur une large plage de valeurs. Enfin, nous mélangeons les phases dans l'espace de Fourier associé à la boîte de simulation périodique, et nous démontrons le rôle crucial de la cohérence de phase pour produire la structure filamentaire observée dans les cartes d'incréments commecelles obtenues récemment par la collaboration Planck / Turbulent energy dissipation is a key process in the cold interstellar medium (ISM), not only on the road to star formation but also as a source of suprathermal energy able to open new chemical routes, otherwise inactive at the low gas temperature. Such routes are required, though, to explain the high abundance of species such as $\mathrm{CH}^+$ and $\mathrm{SH}^+$ observed in the ISM. In this context, the space-time intermittency of energy dissipation is particularly relevant because it drives injection of suprathermal energy in the ISM locally far above the average level. The detailed characteristics of the spatial distribution and the geometrical properties of the energy dissipation rate can provide valuable inputs to chemical models. We study them here with the aid of direct numerical simulations with unprecedented dedication to resolve the dissipation processes numerically. As the sonic Mach number in the diffuse interstellar medium can take values in a wide range, we bracket the possible physics by considering two categories of decaying turbulence models. On the incompressible extreme, we perform pseudo-spectral simulations of viscous and resistive magnetohydrodynamics, with a particular emphasis on ambipolar diffusion due the ion-neutral drift. On the compressible extreme we consider grid-based (Godunov) simulations of isothermal resisitive and viscous magnetohydrodynamics where our focus is on numerical dissipation. Our incompressible simulations show that ambipolar diffusion leads to force-free magnetic fields at small scales. As a result, the typical scale of ion-neutral friction heating is displaced to large scales in the inertial range, much greater than dimensional analysis would predict. The structures of high dissipation are spatially coherent sheets, each with a single nature of dissipation (viscous, ohmic or ambipolar). We reveal their statistical scaling laws and compute their intermittency exponents. We show that compressible simulations are subject to a lot of numerical dissipation: in our set up, less than half of the total dissipation is accounted for by the physical terms, the rest is produced by the numerical scheme. We design a method to recover locally the energy lost in the scheme and we use it to examine the sheet-like structure of the dissipation field as in our incompressible simulations. We show that numerical dissipation prevents us to assess the nature of dissipative structures. For instance, although we confirm previous results that physical dissipation in shearing sheets rapidly dominates over shocks, the balance could be reversed if numerical dissipation were shock dominated. Finally, we examine the efficiency of various observational tracers to characterize the structures of high dissipation. In particular, we find that increments of molecular line centroid velocity or of polarization Stokes parameters correlate very well on the plane of the sky with specific structures of high dissipation. We also compute the intermittency exponents measured for these tracers and find they span a broad range of possible values. At last, we mix the Fourier phases to demonstrate the crucial role of coherence in producing the filamentary structure of observable maps of increments such as recently produced by the Planck collaboration.

La Physique des Hautes Densités d’Énergies (HEDP) est caractérisée par desécoulements multi-matériaux fortement compressibles. Le domaine contenant l’écoulementsubit de grandes variations de taille et est le siège d’ondes de chocs et dedétente intenses. La représentation Lagrangienne est bien adaptée à la descriptionde ce type d’écoulements. Elle permet en effet une très bonne description deschocs ainsi qu’un suivit naturel des interfaces multi-matériaux et des surfaces libres.En particulier, les schémas Volumes Finis centrés Lagrangiens GLACE (GodunovtypeLAgrangian scheme Conservative for total Energy) et EUCCLHYD (ExplicitUnstructured Cell-Centered Lagrangian HYDrodynamics) ont prouvé leur efficacitépour la modélisation des équations de la dynamique des gaz ainsi que de l’élastoplasticité.Le travail de cette thèse s’inscrit dans la continuité des travaux de Maireet Nkonga [JCP, 2009] pour la modélisation de l’hydrodynamique et des travauxde Kluth et Després [JCP, 2010] pour l’hyperelasticité. Plus précisément, cettethèse propose le développement de méthodes robustes et précises pour l’extension3D du schéma EUCCLHYD avec une extension d’ordre deux basée sur les méthodesMUSCL (Monotonic Upstream-centered Scheme for Conservation Laws) et GRP(Generalized Riemann Problem). Une attention particulière est portée sur la préservationdes symétries et la monotonie des solutions. La robustesse et la précision duschéma seront validées sur de nombreux cas tests Lagrangiens dont l’extension 3Dest particulièrement difficile. / High Energy Density Physics (HEDP) flows are multi-material flows characterizedby strong shock waves and large changes in the domain shape due to rarefactionwaves. Numerical schemes based on the Lagrangian formalism are good candidatesto model this kind of flows since the computational grid follows the fluid motion.This provides accurate results around the shocks as well as a natural tracking ofmulti-material interfaces and free-surfaces. In particular, cell-centered Finite VolumeLagrangian schemes such as GLACE (Godunov-type LAgrangian scheme Conservativefor total Energy) and EUCCLHYD (Explicit Unstructured Cell-CenteredLagrangian HYDrodynamics) provide good results on both the modeling of gas dynamicsand elastic-plastic equations. The work produced during this PhD thesisis in continuity with the work of Maire and Nkonga [JCP, 2009] for the hydrodynamicpart and the work of Kluth and Després [JCP, 2010] for the hyperelasticitypart. More precisely, the aim of this thesis is to develop robust and accurate methodsfor the 3D extension of the EUCCLHYD scheme with a second-order extensionbased on MUSCL (Monotonic Upstream-centered Scheme for Conservation Laws)and GRP (Generalized Riemann Problem) procedures. A particular care is taken onthe preservation of symmetries and the monotonicity of the solutions. The schemerobustness and accuracy are assessed on numerous Lagrangian test cases for whichthe 3D extensions are very challenging.