1 |
Convergence du schéma Marker-and-Cell pour les équations de Navier-Stokes incompressible / Convergence of the mac scheme for the incompressible navier-stokes equationsMallem, Khadidja 14 December 2015 (has links)
Le schéma Marker-And-Cell (MAC) est un schéma de discrétisation des équations aux dérivées partielles sur maillages cartésiens, très connu en mécanique des fluides. Nous nous intéressons ici à son analyse mathématique dans le cadre des écoulements incompressibles sur des maillages cartésiens non-uniformes en dimension 2 ou 3. Dans un premier temps nous discrétisons les équations de Navier-Stokes pour un écoulement incompressible stationnaire; nous établissons des estimations a priori sur les suites de vitesses et pressions approchées qui permettent d’une part d'établir l’existence d’une solution au schéma, et d’obtenir la compacité de ces suites lorsque le pas d’espace tend vers 0. Nous montrons alors la convergence de ces suites (à une sous-suite près) vers une solution faible du problème continu, ce qui nécessite une analyse fine du terme de convection non linéaire. Nous nous intéressons ensuite aux équations de Navier-Stokes en régime instationnaire avec une discrétisation en temps implicite. Nous démontrons que le schéma préserve les propriétés de stabilité du problème continu et obtenons ainsi l’existence d’une solution au schéma. Puis, grâce à des techniques de compacité et en passant à la limite dans le schéma, nous démontrons qu’une suite de vitesses approchées converge. Si l’on se restreint au problème de Stokes, et en supposant de plus que la condition initiale de la vitesse est dans H 1 , nous obtenons une estimation sur la pression qui permet de montrer la convergence forte des pressions approchées. Enfin nous étendons l’analyse aux écoulements incompressibles à masse volumique variable. On montre la convergence du schéma. / The Marker-And-Cell (MAC) scheme is a discretization scheme for partial derivative equations on Cartesian meshes, which is very well known in fluid mechanics. Here we are concerned with its mathematical analysis in the case of incompressible flows on two or three dimensional non-uniform Cartesian grids. We first discretize the steady-state incompressible Navier-Stokes equations. We show somea priori estimates that allow to show the existence of a solution to the scheme and some compactness and consistency results. By a passage to the limit on the scheme, we show that the approximate solutions obtained with the MAC scheme converge (up to a subsequence) to a weak solution of the Navier-Stokes equations, thanks to a careful analysis of the nonlinear convection term. Then, we analyze the convergence of the unsteady-case Navier-Stokes equations. The algorithm is implicit in time. We first show that the scheme preserves the stability properties of the continuous problem, which yields, the existence of a solution. Then, invoking compactness arguments and passing to the limit in the scheme, we prove that any sequence of solutions (obtained with a sequence of discretizations the space and time step of which tend to zero) converges up to the extraction of a subsequence to a weak solution of the continuous problem. If we restrict ourselves to the Stokes equations and assume that the initial velocity belongs to H 1, then we obtain estimates on the pressure and prove the convergence of the sequences of approximate pressures. Finally, we extend the analysis of the scheme to incompressible variable density flows. we show the convergence of the scheme.
|
2 |
Méthodes de correction de pression pour les équations de Navier-Stokes compressibles / Pressure correction schemes for compressible flowsKheriji, Walid 28 November 2011 (has links)
Cette thèse porte sur le développement de schémas semi-implicites à pas fractionnaires pour les équations de Navier-Stokes compressibles ; ces schémas entrent dans la classe des méthodes de correction de pression.La discrétisation spatiale choisie est de type "à mailles décalées :éléments finis mixtes non conformes (éléments finis de Crouzeix-Raviart ou Rannacher-Turek) ou schéma MAC classique.Une discrétisation en volumes finis décentrée amont du bilan de masse garantit la positivité de la masse volumique.La positivité de l'énergie interne est obtenue en discrétisant le bilan d'énergie interne continu, par une méthode de volumes finis décentrée amont, enfin, et en couplant ce bilan d'énergie interne discret à l'étape de correction de pression.On effectue une discrétisation particulière en volumes finis sur un maillage dual du terme de convection de vitesse dans le bilan de quantité de mouvement et une étape de renormalisation de la pression; ceci permet de garantir le contrôle au cours du temps de l'intégrale de l'énergie totale sur le domaine.L'ensemble de ces estimations a priori implique en outre, par un argument de degré topologique, l'existence d'une solution discrète. L'application de ce schéma aux équations d'Euler pose une difficulté supplémentaire.En effet, l'obtention de vitesses de choc correctes nécessite que le schéma soit consistant avec l'équation de bilan d'énergie totale, propriété que nous obtenons comme suit. Tout d'abord, nous établissons un bilan discret (local) d'énergie cinétique.Ce dernier comporte des termes sources, que nous compensons ensuite dans le bilan d'énergie interne. Les équations d'énergie cinétique et interne sont associées au maillage dual et primal respectivement, et ne peuvent donc être additionnées pour obtenir un bilan d'énergie totale ; cette dernière équation est toutefois retrouvée, sous sa forme continue, à convergence : si nous supposons qu'une suite de solutions discrètes converge lorsque le pas de temps et d'espace tendent vers 0,, nous montrons en effet, en 1D au moins, que la limite en satisfait une forme faible.Ces résultats théoriques sont confortés par des tests numériques.Des résultats similaires sont obtenus pour les équations de Navier-Stokes barotropes. / This thesis is concerned with the development of semi-implicit fractional step schemes, for the compressible Navier-Stokes equations; these schemes are part of the class of the pressure correction methods.The chosen spatial discretization is staggered: non conforming mixed finite elements (Crouzeix-Raviart or Rannacher-Turek) or the classic MAC scheme. An upwind finite volume discretization of the mass balanced guarantees the positivity of the density. The positivity of the internal energy is obtained by discretising the internal energy balance by an upwind finite volume scheme and by coupling the discrete internal energy balance with the pressure correction step.A special finite volume discretization on dual cells is performed for the convection term in the momentum balance equation, along with a renormalization of the pressure; this allows to guarantee the control in time of integral of the total energy over the domain.All these a priori estimates implies lead to the existence of a discrete solution by a topological degree argument.The application of this scheme the equations of Euler yields an additional difficulty.Indeed, obtaining correct shock speeds requires that the scheme be consistent with the total energy balance,, property which we obtain as follows.First of all, a local discrete kinetic energy balance is established; it contains source terms which are compensated by adding some source terms in the internal energy balance. The kinetic and internal energy equations are associated with the dual and primal meshes respectively, and thus cannot be added to obtain a balance total energy balance; its continuous counterpart is however recovered at the limit: if we suppose that a sequence of discrete solutions converges when the space and time steps tend to 0, we indeed show, in 1D at least, that the limit satisfies a weak form of the equation. These theoretical results are comforted by numerical tests.Similar results are obtained for the barotropic Navier--Stokes equations
|
Page generated in 0.044 seconds