Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport
Les méthodes particulaires sont des méthodes numériques adaptées à la résolution d'équations de conservation. Leur principe consiste à introduire des particules ``numériques'' conservant localement l'inconnue sur un petit volume, puis à les transporter le long de leur trajectoire...
Main Author: | |
---|---|
Other Authors: | |
Language: | fr |
Published: |
2011
|
Subjects: | |
Online Access: | http://www.theses.fr/2011GRENM023/document |
id |
ndltd-theses.fr-2011GRENM023 |
---|---|
record_format |
oai_dc |
collection |
NDLTD |
language |
fr |
sources |
NDLTD |
topic |
Méthodes particulaires Analyse numérique Equation de transport Remaillage Schémas de différences finies TVD Particle methods Numerical analysis Transport equation Remeshing Finite difference schemes TVD 510 |
spellingShingle |
Méthodes particulaires Analyse numérique Equation de transport Remaillage Schémas de différences finies TVD Particle methods Numerical analysis Transport equation Remeshing Finite difference schemes TVD 510 Magni, Adrien Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
description |
Les méthodes particulaires sont des méthodes numériques adaptées à la résolution d'équations de conservation. Leur principe consiste à introduire des particules ``numériques'' conservant localement l'inconnue sur un petit volume, puis à les transporter le long de leur trajectoire. Lorsqu'un terme source est présent dans les équations, l'évolution de la solution le long des caractéristiques est prise en compte par une intéraction entre les particules. Ces méthodes possèdent de bonnes propriétés de conservation et ne sont pas soumises aux conditions habituelles de CFL qui peuvent être contraignantes pour les méthodes Eulériennes. Cependant, une contrainte de recouvrement entre les particules doit être satisfaite pour vérifier des propriétés de convergence de la méthode. Pour satisfaire cette condition de recouvrement, un remaillage périodique des particules est souvent utilisé. Elle consiste à recréer régulièrement de nouvelles particules uniformément réparties, à partir de celles ayant été advectées à l'itération précédente. Quand cette étape de remaillage est effectuée à chaque pas de temps, l'analyse numérique de ces méthodes particulaires remaillées nécessite d'être reconsidérée, ce qui représente l'objectif de ces travaux de thèse. Pour mener à bien cette analyse, nous nous basons sur une analogie entre méthodes particulaires avec remaillage et schémas de grille. Nous montrons que pour des grands pas de temps les schémas numériques obtenus souffrent d'une perte de précision. Nous proposons des méthodes de correction, assurant la consistance des schémas en tout point de grille, le pas de temps étant contraint par une condition sur le gradient du champ de vitesse. Cette méthode est construite en dimension un. Des techniques de limitation sont aussi introduites de manière à remailler les particules sans créer d'oscillations en présence de fortes variations de la solution. Enfin, ces méthodes sont généralisées aux dimensions plus grandes que un en s'inspirant du principe de splitting d'opérateurs. Les applications numériques présentées dans cette thèse concernent la résolution de l'équation de transport sous forme conservative en dimension un à trois, dans des régimes linéaires ou non-linéaires. === Particle methods are numerical methods designed to solve advection dominated conservation equations. Their principle is to introduce ``numerical'' particles that concentrate the unknown locally on a small volume, and to transport them along their trajectories. These methods have good conservation properties and are not subject to the usual CFL conditions that can be binding for the Eulerian methods. However, an overlap condition must be satisfied between the particles to ensure convergence properties of the method. To satisfy this condition, a periodic remeshing of the particles is often used. New particles uniformly distributed are created on a regular mesh. When this remeshing step is performed at every time step, numerical analysis of particle methods needs to be revisited. This is the purpose of this thesis. To carry out this analysis, we rely on an analogy between remeshed particle methods and grid schemes. We show that for large time step the numerical schemes have a loss of accuracy. We propose correction methods wich ensure consistency at any grid point, provided the time step satisfies a condition based on the gradient of the velocity field. Limitation techniques are also introduced to remesh particles without creating any oscillations in the presence of strong variations of the solution. Finally, these methods are generalized to dimensions greater than one. Numerical example on various transport equations are given to illustrate the benefit of the proposed algorithms. |
author2 |
Grenoble |
author_facet |
Grenoble Magni, Adrien |
author |
Magni, Adrien |
author_sort |
Magni, Adrien |
title |
Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
title_short |
Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
title_full |
Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
title_fullStr |
Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
title_full_unstemmed |
Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
title_sort |
méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport |
publishDate |
2011 |
url |
http://www.theses.fr/2011GRENM023/document |
work_keys_str_mv |
AT magniadrien methodesparticulairesavecremaillageanalysenumeriquenouveauxschemasetapplicationspourlasimulationdequationsdetransport AT magniadrien particlemethodswithremeshingnumericalanalysisnewschemesandapplicationsforthesimulationoftransportequations |
_version_ |
1719018579135823872 |
spelling |
ndltd-theses.fr-2011GRENM0232019-04-16T04:03:22Z Méthodes particulaires avec remaillage : analyse numérique nouveaux schémas et applications pour la simulation d'équations de transport Particle methods with remeshing : numerical analysis, new schemes and applications for the simulation of transport equations Méthodes particulaires Analyse numérique Equation de transport Remaillage Schémas de différences finies TVD Particle methods Numerical analysis Transport equation Remeshing Finite difference schemes TVD 510 Les méthodes particulaires sont des méthodes numériques adaptées à la résolution d'équations de conservation. Leur principe consiste à introduire des particules ``numériques'' conservant localement l'inconnue sur un petit volume, puis à les transporter le long de leur trajectoire. Lorsqu'un terme source est présent dans les équations, l'évolution de la solution le long des caractéristiques est prise en compte par une intéraction entre les particules. Ces méthodes possèdent de bonnes propriétés de conservation et ne sont pas soumises aux conditions habituelles de CFL qui peuvent être contraignantes pour les méthodes Eulériennes. Cependant, une contrainte de recouvrement entre les particules doit être satisfaite pour vérifier des propriétés de convergence de la méthode. Pour satisfaire cette condition de recouvrement, un remaillage périodique des particules est souvent utilisé. Elle consiste à recréer régulièrement de nouvelles particules uniformément réparties, à partir de celles ayant été advectées à l'itération précédente. Quand cette étape de remaillage est effectuée à chaque pas de temps, l'analyse numérique de ces méthodes particulaires remaillées nécessite d'être reconsidérée, ce qui représente l'objectif de ces travaux de thèse. Pour mener à bien cette analyse, nous nous basons sur une analogie entre méthodes particulaires avec remaillage et schémas de grille. Nous montrons que pour des grands pas de temps les schémas numériques obtenus souffrent d'une perte de précision. Nous proposons des méthodes de correction, assurant la consistance des schémas en tout point de grille, le pas de temps étant contraint par une condition sur le gradient du champ de vitesse. Cette méthode est construite en dimension un. Des techniques de limitation sont aussi introduites de manière à remailler les particules sans créer d'oscillations en présence de fortes variations de la solution. Enfin, ces méthodes sont généralisées aux dimensions plus grandes que un en s'inspirant du principe de splitting d'opérateurs. Les applications numériques présentées dans cette thèse concernent la résolution de l'équation de transport sous forme conservative en dimension un à trois, dans des régimes linéaires ou non-linéaires. Particle methods are numerical methods designed to solve advection dominated conservation equations. Their principle is to introduce ``numerical'' particles that concentrate the unknown locally on a small volume, and to transport them along their trajectories. These methods have good conservation properties and are not subject to the usual CFL conditions that can be binding for the Eulerian methods. However, an overlap condition must be satisfied between the particles to ensure convergence properties of the method. To satisfy this condition, a periodic remeshing of the particles is often used. New particles uniformly distributed are created on a regular mesh. When this remeshing step is performed at every time step, numerical analysis of particle methods needs to be revisited. This is the purpose of this thesis. To carry out this analysis, we rely on an analogy between remeshed particle methods and grid schemes. We show that for large time step the numerical schemes have a loss of accuracy. We propose correction methods wich ensure consistency at any grid point, provided the time step satisfies a condition based on the gradient of the velocity field. Limitation techniques are also introduced to remesh particles without creating any oscillations in the presence of strong variations of the solution. Finally, these methods are generalized to dimensions greater than one. Numerical example on various transport equations are given to illustrate the benefit of the proposed algorithms. Electronic Thesis or Dissertation Text fr http://www.theses.fr/2011GRENM023/document Magni, Adrien 2011-07-12 Grenoble Cottet, Georges-Henri |