Spelling suggestions: "subject:"fokkerplanck""
1 |
Modélisation et simulation numérique multi-échelle du transport cinétique électroniqueDuclous, Roland 24 November 2009 (has links)
Ce manuscrit est dédié au transport relativiste cinétique sous influence de champs magnétiques, identifié comme obstacle pour la modélisation et la simulation intégrée, dans le cadre de la Fusion par Confinement Inertiel (FCI). Une réalisation importante concerne le développement d'un code déterministe de référence, 2Dx-3Dv, de type Maxwell-Fokker-Planck-Landau, permettant la prise en compte de fonctions de distribution à large degré d'anisotropie. Ce travail se situe à l'interface de l'analyse numérique, des mathématiques appliquées, et de la physique des plasmas. Un deuxième résultat marquant concerne la dérivation d'un modèle collisionel multi-échelle, pour le transport d'électrons relativistes dans la matière dense. Des processus importants sont mis en évidence pour la FCI, et une analogie est menée vis-à-vis des processus de transport collisionels connus en radiothérapie. Enfin, un modèle mésoscopique aux moments angulaires, avec fermeture entropique, a été dérivé et utilisé pour le dépôt de dose pour la radiothérapie. Des schémas numériques précis, d'ordre élevé, et robustes, ont été développé dans ce cadre. / This manuscript is dedicated to the relativistic kinetic transport, under the influence of the magnetic field, identified as a barrier for the modeling and integrated simulations, in the frame of the Inertial Confinement Fusion (ICF). An important achievement concerns the development of a deterministic, reference code, 2Dx-3Dv, of Maxwell-Fokker-Planck-Landau type, that permits the treatment of distribution functions with large anisotropy degree. This work is at the interface between the numerical analysis, applied mathematics, and plasma physics. Another milestone result concerns the derivation of a multi-scale, collisional model, for the transport of relativistic electrons in dense matter. A set of processes is demonstrated to be of importance for ICF, and an analogy is conducted with respect to well-known collisional transport processes in radiotherapy. Finally, a mesoscopic angular moment model, with entropic closure, is derived and employed for radiotherapy dose computation. High order precise and robust numerical schemes are then developed in this framework.
|
2 |
Modélisation et Simulation Numérique multi-échelle du transport cinétique électroniqueDuclous, Roland 24 November 2009 (has links) (PDF)
Ce manuscrit est dédié au transport relativiste cinétique sous influence de champs magnétiques, identifié comme obstacle pour la modélisation et la simulation intégrée, dans le cadre de la Fusion par Confinement Inertiel (FCI). Une réalisation importante concerne le développement d'un code déterministe de référence, 2Dx-3Dv, de type Maxwell-Fokker-Planck-Landau, permettant la prise en compte de fonctions de distribution à large degré d'anisotropie. Ce travail se situe à l'interface de l'analyse numérique, des mathématiques appliquées, et de la physique des plasmas. Un deuxième résultat marquant concerne la dérivation d'un modèle collisionel multi-échelle, pour le transport d'électrons relativistes dans la matière dense. Des processus importants sont mis en évidence pour la FCI, et une analogie est menée vis-à-vis des processus de transport collisionels connus en radiothérapie. Enfin, un modèle mésoscopique aux moments angulaires, avec fermeture entropique, a été dérivé et utilisé pour le dépôt de dose pour la radiothérapie. Des schémas numériques précis, d'ordre élevé, et robustes, ont été développé dans ce cadre.
|
3 |
Contribution à la modélisation et à la simulation numérique multi-échelle du transport cinétique électronique dans un plasma chaudMallet, Jessy 01 October 2012 (has links)
En physique des plasmas, le transport des électrons peut être décrit d'un point de vue cinétique ou d'un point de vue hydrodynamique.En théorie cinétique, une équation de Fokker-Planck couplée aux équations de Maxwell est utilisée habituellement pour décrire l'évolution des électrons dans un plasma collisionnel. Plus précisément la solution de l'équation cinétique est une fonction de distribution non négative f spécifiant la densité des particules en fonction de la vitesse des particules, le temps et la position dans l'espace. Afin d'approcher la solution de ce problème cinétique, de nombreuses méthodes de calcul ont été développées. Ici, une méthode déterministe est proposée dans une géométrie plane. Cette méthode est basée sur différents schémas numériques d'ordre élevé . Chaque schéma déterministe utilisé présente de nombreuses propriétés fondamentales telles que la conservation du flux de particules, la préservation de la positivité de la fonction de distribution et la conservation de l'énergie. Cependant, le coût de calcul cinétique pour cette méthode précise est trop élevé pour être utilisé dans la pratique, en particulier dans un espace multidimensionnel.Afin de réduire ce temps de calcul, le plasma peut être décrit par un modèle hydrodynamique. Toutefois, pour les nouvelles cibles à haute énergie, les effets cinétiques sont trop importants pour les négliger et remplacer le calcul cinétique par des modèles habituels d'Euler macroscopiques. C'est pourquoi une approche alternative est proposée en considérant une description intermédiaire entre le modèle fluide et le modèle cinétique. Pour décrire le transport des électrons, le nouveau modèle réduit cinétique M1 est basé sur une approche aux moments pour le système Maxwell-Fokker-Planck. Ce modèle aux moments utilise des intégrations de la fonction de distribution des électrons sur la direction de propagation et ne retient que l'énergie des particules comme variable cinétique. La variable de vitesse est écrite en coordonnées sphériques et le modèle est défini en considérant le système de moments par rapport à la variable angulaire. La fermeture du système de moments est obtenue sous l'hypothèse que la fonction de distribution est une fonction d'entropie minimale. Ce modèle satisfait les propriétés fondamentales telles que la conservation de la positivité de la fonction de distribution, les lois de conservation pour les opérateurs de collision et la dissipation d'entropie. En outre une discrétisation entropique avec la variable de vitesse est proposée sur le modèle semi-discret. De plus, le modèle M1 peut être généralisé au modèle MN en considérant N moments donnés. Le modèle aux N-moments obtenu préserve également les propriétés fondamentales telles que les lois de conservation et la dissipation de l'entropie. Le schéma semi-discret associé préserve les propriétés de conservation et de décroissance de l'entropie. / In plasma physics, the transport of electrons can be described from a kinetic point of view or from an hydrodynamical point of view.Classically in kinetic theory, a Fokker-Planck equation coupled with Maxwell equations is used to describe the evolution of electrons in a collisional plasma. More precisely the solution of the kinetic equations is a non-negative distribution function f specifying the density of particles as a function of velocity of particles, the time and the position in space. In order to approximate the solution of such problems, many computational methods have been developed. Here, a deterministic method is proposed in a planar geometry. This method is based on different high order numerical schemes. Each deterministic scheme used presents many fundamental properties such as conservation of flux particles, preservation of positivity of the distribution function and conservation of energy. However the kinetic computation of this accurate method is too expensive to be used in practical computation especially in multi-dimensional space.To reduce the computational time, the plasma can be described by an hydrodynamic model. However for the new high energy target drivers, the kinetic effects are too important to neglect them and replace kinetic calculus by usual macroscopic Euler models.That is why an alternative approach is proposed by considering an intermediate description between the fluid and the kinetic level. To describe the transport of electrons, the new reduced kinetic model M1 proposed is based on a moment approach for Maxwell-Fokker-Planck equations. This moment model uses integration of the electron distribution function on the propagating direction and retains only the energy of particles as kinetic variable. The velocity variable is written in spherical coordinates and the model is written by considering the system of moments with respect to the angular variable. The closure of the moments system is obtained under the assumption that the distribution function is a minimum entropy function. This model is proved to satisfy fundamental properties such as the non-negativity of the distribution function, conservation laws for collision operators and entropy dissipation. Moreover an entropic discretization in the velocity variable is proposed on the semi-discrete model. Moreover the M1 model can be generalized to the MN model by considering N given moments. The N-moments model obtained also preserves fundamental properties such as conservation laws and entropy dissipation. The associated semi-discrete scheme is shown to preserve the conservation properties and entropy decay.
|
4 |
Analyse mathématique et numérique de modèles gyrocinétiques / Mathematical and numerical analysis of gyro-kinetic modelsCaldini-Queiros, Céline 15 November 2013 (has links)
Cette thèse porte sur les équations gyro-cinétiques et traite un développement rigoureux deslimites de l'équation de Vlasov avec différents opérateurs de collision dans un champ magnétiquefort, ainsi que du développement de méthodes numériques.On commence par une étude de l'opérateur de moyenne. L'opérateur de moyenne a été développé parM. Bostan dans le cadre général d'une équation pour laquelle une partie du transport estfortement pénalisée. Puis, on applique ces résultats généraux aux deux régimes limites que nousétudions : le régime du rayon de Larmor fini et le régime centre-guide.On s'intéresse au calcul précis et explicite de la moyenne de l'opérateur de Fokker-Planck-Landau. On se place pour cela dans le cas du régime du rayon de Larmor fini. Avant de réaliserles calculs sur l'opérateur de Fokker-Planck-Landau, qui contient des convolutions et des termesde diffusion, il semble raisonnable de calculer la moyenne de l'opérateur de relaxation deBoltzmann, dont l'expression est plus simple.On se place ensuite dans le cas du régime centre-guide et on présente un schéma numérique basésur une décomposition micro-macro de la fonction de distribution des particules qui provientd'un travail en collaboration avec N. Crouseilles et M. Lemou. On obtient un schéma uniformémentconsistant avec le modèle continu, pour tout ordre du champ magnétique. Des simulationsnumériques, basées sur cette approche, ont été réalisées à l'aide d'un code de calcul 2D quel'on a développé durant cette thèse.On présente ensuite un projet réalisé dans le cadre du Cemracs 2012, consacré à la modélisationdes écoulements sanguins dans le réseau veineux cérébral. / The main subject of this thesis is the gyro-kinetic equation. We present a rigourousdeveloppement of the Vlasov equation limits with different collision operator in a strongmagnetic field and numerical methods.We start with a study of the gyro-average operator. The average operator has been introduced byM. Bostan in the case of an equation where part of the transport is highly penalised. Then weapply our results at the two approximation we study : the finite Larmor radius approximation andthe guiding-center approximation.We first focus on the precise and explicit computation of the Fokker-Planck-Landau operatoraverage in the finite Larmor radius approximation. The Fokker-Planck-Landau operator containsconvolution and diffusion terms, it is then reasonable to first compute the average of theBoltzmann relaxation operator.We then focus on the guiding-center approximation and present a numerical scheme based on amicro-macro decomposition of the particles distribution fonction which comes from a joint workwith N. Crouseilles and M. Lemou. We obtain a scheme which is uniformly consistant with thecontinuous model for any order of the magnetic field. Numerical simulation based on thisapproach are presented.The last chapter of this thesis presents a project which was realised during the Cemracs 2012concerning the modelisation of blood flow in cerebral veins.
|
5 |
Étude et discrétisation de modèles cinétiques et de modèles fluides à bas nombre de MachDellacherie, Stéphane 02 February 2011 (has links) (PDF)
Ce mémoire résume les travaux que nous avons réalisés de 1995 à 2010. Ces travaux ont eu pour thème l'étude et la discrétisation, d'une part, de modèles cinétiques de type Fokker-Planck ou de type Boltzmann semi-classiques et, d'autre part, de modèles fluides de type Euler ou de type Navier-Stokes à bas nombre de Mach. L'équation de Fokker-Planck étudiée modélise les collisions entre ions et électrons dans un plasma chaud, et concerne ici la fusion par confinement inertiel. Les équations de Boltzmann semi-classiques étudiées sont de deux types. Le premier type modélise la réaction de fusion thermonucléaire entre un ion deuterium et un ion tritium donnant une particule alpha et un neutron, et concerne également ici la fusion par confinement inertiel. Le deuxième type - connu sous le nom d'équations de Wang-Chang & Uhlenbeck - modélise ici les transitions d'énergie quantique dans les couches électroniques d'atomes d'uranium et de fer provoquées par les collisions entre ces mêmes atomes au sein du procédé SILVA de Séparation Isotopique par Laser sur Vapeur Atomique. Nous avons étudié les propriétés de base de ces deux types d'équations de Boltzmann semi-classiques, et, dans le cas des équations de Wang-Chang & Uhlenbeck, nous avons proposé un algorithme de couplage cinétique-fluide. L'étude de cet algorithme nous a incité à étudier la notion de relaxation dans un mélange binaire de gaz et de fluides non-miscibles, et à souligner les points communs de cette approche avec la théorie cinétique standard. L'étude de modèles moyennés pour des mélanges de fluides non-miscibles nous a amené à proposer et à discrétiser un modèle sans ondes acoustiques modélisant la déformation d'une interface entre deux fluides non-miscibles provoquée par de forts gradients thermiques à bas nombre de Mach. Puis, afin d'améliorer la précision des calculs tout en en maîtrisant le coût, nous avons également étudié la possibilité de résoudre sur un maillage dynamique de type AMR un modèle simplifié de déformation d'interface. Ces études à bas nombre de Mach nous ont également incités à analyser sur maillage cartésien le mauvais comportement à bas nombre de Mach des schémas de type Godunov appliqués au système d'Euler compressible. Enfin, nous avons justifié l'algorithme LBM dans le cas de l'équation de la chaleur.
|
Page generated in 0.0503 seconds