Notre étude vise à modéliser l’écoulement pulsé sanguin dans le réseau vasculaire humain. Celui ci est constitué d’un très grand nombre de vaisseaux disposés dans un vaste réseau ayant différentes propriétés mécaniques. Le modèle simplifié unidimensionnel (1D) permet une étude numérique dans le réseau sanguin et plusieurs applications seront proposées.Le modèle 1D est établi grâce aux hypothèses de grande longueur d’onde de l’onde de pouls comparée aux rayons des vaisseaux et de profil de vitesse de révolution, en moyennant transversalement les équations de Navier-Stokes et de conservation de la masse. Un modèle viscoélastique de Kelvin-Voigt est adopté pour l’équation constitutive du tube. Cela conduit à un système hyperbolique-parabolique non linéaire, qui est ensuite résolu avec quatre schémas numériques, à savoir: MacCormack, Taylor-Galerkin, schéma monotone décentré pour les équations de loi de conservation (MUSCL) et Galerkin discontinu local. Les schémas sont mis en oeuvre dans un premier temps dans MATLAB et les solutions numériques sont vérifiées favorablement à des solutions semi-analytiques et des observations cliniques. Des comparaisons entre les schémas sont faites pour quatre aspects importants: la précision, la capacité de capturer desphénomènes de type choc, la vitesse de calcul et la complexité de la mise en oeuvre, enfin les conditions appropriées pour l’application de chaque système sont discutées. Après cela, un code objet général en C++ est développé et testé sur plusieurs réseaux: un cercle d’artères, un réseau systémique humain de 55 artères et un rein de souris avec plus d’un millier les segments. La répartition en fonction du temps de la pression dans les réseaux est visualisée et les modes de propagation des ondes sont bien capturés. Une bonne accélération est atteinte par parallélisation du code.Le code développé est ensuite appliqué dans trois études. En premier lieu, les coefficients de frottement du fluide et la viscosité de la paroi sontdéterminés avec des dispositifs expérimentaux bien définis constitués de tuyaux élastiques in vitro. Ces deux facteurs amortissant les ondes de pouls, ils sont difficiles à évaluer séparément. Nous les estimons par ajustement du modèle viscoélastique 1D avec les ondes de pression mesurées expérimentalement. Les valeurs ajustées des paramètres viscoélastiques sont conformes aux valeurs estimées avec d’autres méthodes. Les deux effets visqueux sont du même ordre de grandeur. In vivo, des séries chronologiques de la pression du diamètre en différents points d’un réseau artériel de mouton, sont analysées et les paramètres de viscoélasticité sont estimés. Le réseau du mouton est ensuite simulé, on montre que la viscoélasticité amortit de manière significative les hautes fréquences. En troisième lieu, la variation de la circulation induite par des anastomoses axillo- et fémoro-fémorales avec une sténose iliaquesévère est simulée. L’influence de la voie de contournement est étudié. / The vascular network consists of a very large number of segments with various properties and thus the pulsatile blood flow inside is very complicated. With the time-domain-based nonlinear 1D model, this thesis studies the blood flow in networks, focusing on the numerical computing and several applications.With assumptions of long wave and axisymmetric velocity profile, the 1D governing equations of mass and momentum are derived by integrating the continuity and Navier-Stokes equations along the radius.A Kelvin-Voigt viscoelastic model is adopted for the constitutive equation of the tube.This leads to a nonlinear hyperbolic-parabolic system, which is then solved with four numerical schemes, namely: MacCormack, Taylor-Galerkin, Monotonic Upwind Scheme for Conservation Law (MUSCL) and local discontinuous Galerkin.The schemes are implemented in MATLAB and the numerical solutions are checked favorably against analytical, semi-analytical solutions and clinical observations.Among the numerical schemes, comparisons are made in four important aspects: accuracy, ability to capture shock-like phenomena, computational speed and implementation complexity. The suitable conditions for the application of each scheme are discussed.After this, a general purpose C++ code is developed and tested on several networks:a circle of arteries, a human systemic network with 55 arteries and a mouse kidney with more than one thousand segments. The time dependent distribution of pressure in the networks is visualized and the propagation patterns of the waves are well captured.Good speedup is achieved by parallelizations of the code.The developed code is applied in three studies.First, the coefficients of fluid friction and wall viscosity are determined with aides of a well defined experimental setup.Because both the two factors damp the pulse waves, they are difficult to evaluate separately.We estimate them in pairs by fitting the 1D viscoelastic model against pressure waves measured on the experimental setup. The fitted values of viscoelastic parameters are consistent with values estimated with other methods.The effect of wall viscosity on the pulse wave has been shown in the same order of that of fluid viscosity.Second, with time series of pressure and diameter measured in several locations of the sheep arterial network,the viscoelasticity parameters are estimated.With those values, the pulse waves in the sheep network are simulated and the effect of viscoelasticity is investigated.Numerical solutions show that the viscoelasticity damps significantly the high frequency components of the pulse waves.Third, we simulate the change of blood flow induced by the axillofemoral and femoral-femoral anastomoses with a severe iliac stenosis.The influence of the bypassing path is studied.

Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark « Ecoulement diphasique » proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé / The two-phase flow in porous media is a complex phenomenon and which relate to many industrial problems. EDF works on the feasibility and the safety of a storage in deep geologic layer of nuclear waste. In this domain the simulation of the two-phase flow in porous media is particularly important in at least three domains : first of all during the phase of ventilation of the galleries of the storage which could desaturate the rock and so modify its properties, but also during the phase of resaturation of the materials and finally during the arrival of the water on the metal parts contained in the storage which will then involve phenomena of corrosion and a hydrogen release. In this context, EDF wishes to obtain robust numerical methods without restrictive condition on the mesh. This work is dedicated at first to the development of the finite volume scheme SUSHI (Scheme Using Stabilization and Hybrid Interfaces) in the code of mechanics of EDF, Code_Aster in order to simulate the two-phase flow in porous media. This scheme was developed in 2D and in 3D. At the same time a new formulation which allows to simulate in a uniform way the flows in saturated and unsaturated porous media for miscible and immiscible problems is proposed. Various studies simulating difficulties related to the problems of the storage of nuclear waste in deep geological layers were study. We can quote the study of a bi-material which advances the capillary rebalancing of a material by an other one possessing properties and initial very heterogeneous conditions in saturation. We will also quote the study of the injection of hydrogen in an porous media initially saturated in pure water which is proposed by the benchmark "two-phase Flow " proposed by the GNR MOMAS. This study had for objective to bring to light the good treatment of the appearance of a phase in a saturated porous media and thus the relevance of our new formulation to study with a way unified a problem of saturated flow and a problem of unsaturated flow

Dans cette thèse, nous nous intéressons à l’analyse mathématique théorique et numérique des équations deNavier-Stokes compressibles en régime barotrope. La plupart des travaux présentés ici combinent desméthodes d’analyse des équations aux dérivées partielles et des méthodes d’analyse numérique afin de clarifierla notion de solution faible ainsi que les mécanismes de convergence de méthodes numériques approximant cessolutions faibles. En effet les équations de Navier-Stokes compressibles sont fortement non linéaires et leuranalyse mathématique repose nécessairement sur la structure de ces équations. Plus précisément, nousprouvons dans la partie théorique l’existence de solutions faibles pour un modèle d’écoulement compressibled’entropie variable où l’entropie du système est transportée. Nous utilisons les méthodes classiques permettantde prouver l’existence de solutions faibles aux équations de Navier-Stokes compressibles en regime barotrope.Nous étudions aussi dans cette partie la réduction de dimension 3D/2D dans les équations de Navier-Stokescompressibles en utilisant la méthode d’énergie relative. Dans la partie numérique nous nous intéressons auxestimations d’erreur inconditionnelles pour des schémas numériques approximant les solutions faibles deséquations de Navier-Stokes compressibles. Ces estimations d’erreur sont obtenues à l’aide d’une versiondiscrète de l’énergie relative satisfaite par les solutions discrètes de ces schémas. Ces estimations d’erreur sontobtenues pour un schéma numérique académique de type volumes finis/éléments finis ainsi que pour le schémanumérique Marker-and-Cell. Nous prouvons aussi que le schéma Marker-and-Cell est inconditionnellement etuniformément asymptotiquement stable en régime bas Mach. Ces résultats constituent les premiers résultatsd’estimations d’erreur inconditionnelles pour des schémas numériques pour les équations de Navier-Stokescompressibles en régime barorope. / In this thesis, we deal with mathematical and numerical analysis of compressible Navier-Stokes equations inbarotropic regime. Most of these works presented here combine mathematical analysis of partial differentialequations and numerical methods with aim to shred more light on the construction of weak solutions on oneside and on the convergence mechanisms of numerical methods approximating these weak solutions on theother side. Indeed, the compressible Navier-Stokes equations are strongly nonlinear and their mathematicalanalysis necessarily relies on the structure of equations. More precisely, we prove in the theorical part theexistence of weak solutions for a model a flow of compressible viscous fluid with variable entropy where theentropy is transported. We use the classical techniques to prove the existence of weak solutions for thecompressible Navier-Stokes equations in barotropic regime. We also investigate the 3D/2D dimensionreduction in the compressible Navier-Stokes equations using the relative energy method. In the numerical wedeal with unconditionally error estimates for numerical schemes approximating weak solutions of thecompressible Navier-Stokes equations. These error estimates are obtained by using the discrete version of therelative energy method. These error estimates are obtained for a academic finite volume/finite element schemeand for the Marker-and-Cell scheme. We also prove that the Marker-and-cell scheme is unconditionally anduniformly asymptotically stable at the Low Mach number regime. These are the first results onunconditionally error estimates for numerical schemes approximating the compressible Navier-Stokesequations in barotropic regime.

Ce travail de thèse s'inscrit dans une problématique visant à étudier numériquement le bruit d'origine aérodynamique généré par les écoulements turbulents autour des véhicules en utilisant la méthode Boltzmann sur Réseau (LBM). Les objectifs de cette thèse sont l'étude des capacités aéroacoustiques de la LBM ainsi que l'élaboration d'un code de calcul tridimensionnel et parallèle.<br />Dans un premier temps, les élements historiques et théoriques de la LBM sont présentés ainsi que le développement permettant de passer de l'équation de Boltzmann aux équations de Navier-Stokes. La construction des modèles à vitesses discrètes est également décrite. Deux modèles basés sur des opérateurs de collision différents sont présentés : le modèle LBM-BGK et le modèle LBM-MRT. Pour l'étude des capacités aéroacoustiques de la LBM, une analyse de von Neumann est réalisée pour les modèles LBM-BGK et LBM-MRT ainsi que pour l'équation de Boltzmann à vitesse discrète (DVBE). Une comparaison avec les schémas Navier-Stokes d'ordre élevé est alors menée. Pour remédier aux instabilités numériques de la méthode Boltzmann sur Réseau intervenant lors de la propagation dans des directions particulières à M>0.1, des filtres sélectifs sont utilisés et leur effet sur la dissipation est étudié.<br />Dans un second temps, le code de calcul L-BEAM est présenté. La structure générale et les différentes techniques de calculs sont décrites. Un algorithme de transition de résolution est développé. La modélisation de la turbulence est abordée et le modèle de Meyers-Sagaut est implémenté dans le code. Enfin, des cas tests numériques sont utilisés pour valider le code et la simulation d'un écoulement turbulent complexe est réalisée.

Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark " Ecoulement diphasique " proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé

Les travaux de cette thèse, réalisés au sein de l’équipe CLEF/DMFN de l’ONERA (Office National d’ Etudes et de Recherches Aérospatiales) en partenariat avec le laboratoire DynFluid et le CIRT (Consortium Industrie-Recherche en Turbomachines) s’inscrivent dans une demarche d’amélioration des outils de simulations pour les turbomachines. Compte tenu de ce contexte, l’objectif de cette étude est de développer de nouvelles méthodes pour le traitement des raccords non coincidents dans les turbomachines qui soit à la fois d’ordre elevé et conservatifs. Les développements proposés sont validés et composés de configurations de difficulté croissante. / A high-order and conservative method is developed for the numerical treatment of interface conditions in patched grids, based on the use of a ctitious grid methodology. The proposed approach is compared with a non-conservative interpolation of the state variables from the neighbouring domain for selected internal fow problems.

Dans cette thèse, on étudie la stabilité orbitale d’états stationnaires de modèles mathématiques de type "Hamiltonian mean-field", dits modèles HMF. Cette étude est d’abord menée d’un point de vue théorique en utilisant des méthodes variationnelles. Puis, elle est menée d’un point de vue numérique en commençant par l’élaboration de schémas conservant exactement des états stationnaires. Le Chapitre 2 présente une étude théorique de la stabilité orbitale des états stationnaires du modèle HMF Poisson. Plus précisément, on prouve la stabilité orbitale d’une grande classe d’états stationnaires solutions du système HMF avec potentiel de Poisson. Ces états stationnaires sont des minimiseurs d’un problème à une, deux ou une infinité de contraintes d’une certaine fonctionnelle. La preuve s’appuie sur une approche variationnelle. Cependant le caractère borné du domaine empêche l’utilisation des techniques usuelles basées sur des invariances d’échelles. On introduit alors de nouvelles méthodes, spécifiques à ce problème, mais demeurant dans l’esprit des outils de réarrangements introduits pour le système de Vlasov-Poisson. En particulier, ces méthodes permettent de considérer un nombre arbitraire de contraintes et aboutissent à un résultat de stabilité pour une grande classe d’états stationnaires. Dans le Chapitre 3, on construit des schémas numériques conservant exactement des états stationnaires donnés. Ces schémas modélisent mieux la propriété de stabilité orbitale que les schémas classiques. Puis, on propose un schéma plus général en construisant un schéma qui conserve tous les états stationnaires des modèles HMF. Pour finir, à l’aide de ces schémas, est menée une étude numérique de la stabilité des états stationnaires du système de HMF Poisson qui vient compléter l’étude théorique du Chapitre 2. / In this thesis, we study the nonlinear orbital stability of steady states of "Hamiltonian mean-field" models, called HMF models. First, this study is being done theoretically by using variational methods. It is then carried out numerically by building numerical schemes wich exactly preserve steady states. Chapter 2 presents a theoretical study of the orbital stability of steady states which are solutions to the HMF Poisson system. More specifically, the orbital stability of a large class of steady states which are solutions to the HMF system with Poisson potential is proved. These steady states are obtained as minimizers of an energy functional under one, two or infinitely many constraints. The proof relies on a variational approach. However the boundedness of the space domain prevents us from using usal technics based on scale invariance. Therefore, we introduce new methods which, although specific to our context, remain somehow in the same spirit of rearrangements tools introduced for the Vlasov-Poisson system. In particular, these methods allow for the incorporation of an arbitrary number of constraints, and yield a stability result for a large class of steady states. In Chapter 3, numerical schemes exactly preserving given steady states are built. These schemes model the orbital stability property better than the classic ones. Then, a more general scheme is introduced by building a scheme wich preserves all steady states of HMF models. Lastly, by means of these schemes, we conduct a numerical study of stability of steady states solutions to HMF Poisson system. This completes the theoretical study in Chapter 2.

Cette thèse comporte trois volets indépendants en cryptographie, en théorie de Hodge p-adique et en analyse numérique.La première partie consiste en l'étude d'algorithmes performants de résolution du logarithme discret. La résolution du logarithme discret consiste à déterminer les exposants d'une famille fixée de générateurs dans la décomposition des éléments du groupe. Dans le cas des groupes multiplicatifs d'un corps fini, la complexité des calculs dépendent de la taille - dite de petite, moyenne ou grande caractéristique- de la caractéristique du corps dans lesquels on effectue les calculs.Nous présentons différents algorithmes dans chacune des caractéristiques (petite, moyenne ou grande) en précisant quel est l'algorithme le plus performant dans chacun des cas.La seconde partie s'inscrit dans le contexte du programme de Langlands p-adique. Nous présentons une généralisation de l'un des outils centraux de la théorie, les modules de Breuil-Kisin, en plusieurs variables La troisième partie est un travail effectué en collaboration avec Victor Vilaça Da Rocha, Roberta Tittarelli, Richard Sambilason Rafefimanana, Victor Michel-Dansac et Benjamin Couéraud. Il a été initié lors de la treizième SEME, Semaine d'Etudes Maths Entreprises organisée par l'Agence pour les Mathématiques en Interaction avec l'Entreprise et la Société (AMIES).L'Institut Français du Pétrole et des Energies Nouvelles nous a soumis un problème de résolution numérique d'un système d'équations modélisant la désorption d'un gaz de schiste en une dimension.Nous proposons plusieurs schémas du premier ordre recourant à un traitement implicite de l'équation de relaxation. Enfin nous présentons un schéma numérique d'ordre deux en temps. / In this thesis, we present three independent works in cryptography, p-adic Hodge theory and Numerical analysis.First we present several algorithms to solve the discrete logarithm in several characteristic finite fields. We are particularly interested with the determination of classes of polynomial functions with small coefficients.The second part of the thesis deals with one of the major object of p-adic Hodge theory. We present a multi-variable version of Breuil-Kisin modules where the Lubin-Tate tower replaces the classical cyclotomic tower. He third proposes two numerical schemes for the modelisation of desorption of shale gaz.

La Simulation aux Grandes Echelles (SGE) est de plus en plus utilisée dans les processus de développement et la conception des réacteurs aéronautiques industriels. L'une des raisons pour ce besoin résulte dans la capacité de la SGE à fournir des informations instantanées d'un écoulement turbulent augmentant la quantité des prédictions de la composition des gaz d'échappement. Ce manuscrit de thèse aborde deux sujets récurrents de la SGE. D'une part, les schémas numériques pour la SGE nécessitent certaines propriétés, notamment une précision élevée avec une diffusivité faible pour ne pas nuire aux modèles de turbulence. Afin de répondre à ce pré requis, une famille de schémas d'intégration temporelle d'ordre élevée est proposée, permettant de modifier la diffusion numérique du schéma. D'autre part, la SGE étant intrinsèquement instationnaire, elle est très consommatrice en temps CPU. De plus, une géométrie complexe prend beaucoup de temps de simulation même avec les super calculateurs d'aujourd'hui. Dans le cas particulier d'intérêt et souvent rencontré dans les applications industrielles, l'approche bas-Mach est constitue une alternative intéressante permettant de réduire le coût et le temps de retour d'une simulation LES. L'impact et la comparaison des formalismes compressible et incompressible sont toutefois rarement quantifiés, ce qui est proposé dans ce travail pour une configuration représentative d'un brûleur swirlé industriel mesuré au CORIA / Large-Eddy Simulation (LES) becomes a more and more demanded tool to improve the design of aero-engines. The main reason for this request stems from the constraints imposed on the next generation low-emission engines at the industrial development level and the ability for LES to provide information on the instantaneous turbulent flow field which greatly contributes to improving the prediction of mixing and combustion thereby offering an improved prediction of the exhaust emission. The work presented in this thesis discusses two recurring issues of LES. For one, numerical schemes for LES require certain properties, i.e. low-diffusion schemes of high order of accuracy so as not to interfere with the turbulence models. To meet this purpose in the context of fully unstructured solvers, a new family of high-order time-integration schemes is proposed. With this class of schemes, the diffusion implied by the numerical scheme become adjustable and built-in. Second, since fully unsteady by nature, LES is very consuming in terms of CPU time. Even with today's supercomputers complex problems require long simulation times. Due to the low flow velocities often occurring in industrial applications, the use of a low-Mach number solver seems suitable and can lead to large reductions in CPU time if comparable to fully compressible solvers. The impact of the incompressibility assumption and the different nature of the numerical algorithms are rarely discussed. To partly answer the question, detailed comparisons are proposed for an experimental swirled configuration representative of a real burner that is simulated by LES using a fully explicit compressible solver and an incompressible solution developed at CORIA

Mes travaux de thèse portent sur le développement de schémas numériques d'ordre élevé en temps et en espace pour la simulation de propagation des ondes. Nous avons proposé de discrétiser dans un premier temps l'équation des ondes par rapport au temps, en utilisant une technique de type équation modifiée. Puis nous avons utilisé une méthode d'éléments finis de type Galerkine discontinue pour la discrétisation en espace. En modifiant l'ordre de la discrétisation, nous avons construit des schémas tout aussi précis que ceux déjà existants pour un coût de mise en oeuvre très intéressant. Après avoir validé numériquement la nouvelle méthode, nous nous sommes intéressés à sa stabilité ainsi qu'à son adaptivité en temps et en espace. Pour arriver à cela, nous avons dû faire une étude précise de la stabilité de la méthode de Galerkine discontinue et nous avons proposé des améliorations à cette technique entraînant des gains de temps significatifs. / My work consists in developing some high order numerical schemes in time and space for the modeling of the wave propagation. We have proposed to first discretize the wave equation with respect to the time using the so called Modified Equation Technique. Then, we have used a Discontinuous Galerkine Finite Element method for the space discretization. Switching the classical discretization process, we have constructed schemes as accurate as the classical ones with a numerical cost very interesting. After the numerical validation of this method, we have focused on its stability and on its adaptibility in time and space. To reach these objectives, we have performed a stability analysis of the Discontinuous Galerkin method and we have proposed some improvements to this technique which imply very important gain in terms of computationnal time.