Spelling suggestions: "subject:"équations""
261 |
Spécialisation du pseudo-groupe de Malgrange et irréductibilité / Specialisation of the Malgrange pseudogroup and irreductibilityDavy, Damien 13 December 2016 (has links)
Le pseudo-groupe de Malgrange d'un champ de vecteurs défini sur une variété est la sous-pro-variété de l'espace des jets de biholomorphismes locaux de cette variété obtenue en prenant la clôture de Zariski des flots du champ de vecteurs. Une équation différentielle ordinaire d'ordre 2 définit un champ de vecteurs sur une variété de dimension 3. Le pseudogroupe de Malgrange de ce dernier est de type différentiel d'ordre inférieur ou égal à 2. Une équation différentielle ordinaire d'ordre 2 est dite irréductible si ses solutions générales ne peuvent pas être exprimées à l'aide de solutions d'équations algébriques, différentielles linéaires ou différentielles d'ordre 1. Si le type différentiel du pseudo-groupe de Malgrange d'une équation d'ordre 2 est exactement 2 alors cette dernière est irréductible. Nous donnons plusieurs définitions du pseudo-groupe de Malgrange d'un champ de vecteurs équivalentes à la définition originale donnée par Bernard Malgrange. La définition du premier paragraphe nous permet d'appliquer un théorème de semi-continuité de la dimension des clôtures de Zariski des feuilles d'un feuilletage holomorphe de Philippe Bonnet. Nous obtenons le résultat suivant concernant les équations différentielles ordinaires dépendant de paramètres. Si le type différentiel du pseudo-groupe de Malgrange de l'équation spécialisée en une valeur des paramètres est à exactement 2 alors il en sera de même pour les pseudo-groupes de Malgrange de l'équation spécialisée en des valeurs générales des paramètres. Une première application de ce résultat est de redémontrer l'irréductibilité des équations de Painlevé pour des valeurs générales des paramètres. Une seconde application est de déterminer complètement les pseudo-groupes de Malgrange de ces équations pour des valeurs générales des paramètres. Les définitions du pseudo-groupe de Malgrange et les résultats de spécialisations s'adaptent aux équations aux q-différences. En appliquant ces résultats aux équations de Painlevé discrètes, nous obtenons le pseudo-groupe de Malgrange de ces dernières pour des valeurs générales des paramètres. / The Malgrange pseudogroup of a vector field on a variety is the sub-pro-variety of the jet space of local biholomorphisms of this variety obtained by taking the Zariski closure of the flow of the vector field. A second-order ordinary differential equation defines a vector field on a variety of dimension 3. The differential type of the Malgrange pseudogroup of this one is at most 2. A second-order ordinary differential equation is said to be irreductible if its general solutions can not be expressed using solutions of algebraic equations, linear differential equations or differential equations of order 1. If the differential type of the Malgrange pseudogroup of a second-order differential equation is exactly 2 then the latter is irreductible. We give several definitions of the Malgrange pseudogroup of a vector field which are equivalent to the original definition given by Bernard Malgrange. The definition of the first paragraph leads us to apply a semi-continuity theorem of the dimension of the Zariski closure of the leaves of a holomorphic foliation given by Philippe Bonnet. We obtain the following result about the ordinary differential equations which depend on parameters. If the differential type of the Malgrange pseudogroup of the equation specialized in one value of parameters is exactly two then it will be the same for the Malgrange pseudogroup of the equation specialized in a general value of parameters. A first application of this result is an other proof of the irreductibility of the Painlevé equations for general value of parameters. A second application is to fully determined the Malgrange pseudogroups of this equations for general value of parameters. The definitions of the Malgrange pseudogroup of a vector field and the specialisation results can be adapted the q-difference equations. By applying this results to the discret Painlevé equations, we fully determined the Malgrange pseudogroup of the latters for general value of parameters.
|
262 |
Interférométrie atomique avec un condensat de Bose-Eintein : effet des interactions internes / Atom interferometry with a Bose-Einstein condensate : effect of internal interactionsJannin, Raphaël 08 October 2015 (has links)
Le travail réalisé dans le cadre de cette thèse s'articule en deux volets. Le premier porte sur l'étude de l'effet des interactions entre atomes au sein d'un interféromètre atomique, dont la source est un condensat de Bose-Eintein. Nous présentons un modèle analytiquepermettant d'obtenir des expressions simples pour le déphasage induit par celles-ci. Ce modèle est comparé à des simulations numériques résolvant les équations de Gross-Pitaevskii couplées, et présente un excellent accord. Le second concerne la conception et la construction d'un nouveau dispositif expérimental visant à obtenir un condensat de Bose-Einteindans le but de réaliser des mesures de haute précision par interférométrie atomique. / The work performed during this thesis comprises two orientations. The first one is the study of the effect of interactions between atoms in an atom interferometer which source of atoms is a Bose-Einstein condensate. We present an analytical model allowing to obtain simple expressions for the phase shift induced by them. This model is compared to numerical simulations solving the coupled Gross-Pitaevskii equations and presents a good agreement. The second one is the design and construction of a new experimental set-up for the production of a Bose-Einstein condensate to perform high precision measurements with the use of atom interferometry.
|
263 |
Réflexion quantique sur le potentiel de Casimir-Polder / Quantum reflection from the Casimir-Polder potentialDufour, Gabriel 20 November 2015 (has links)
Les collisions entre atomes ultrafroids et surfaces matérielles sont caractérisées par la réflexion de l'onde de matière atomique sur le potentiel attractif de Casimir-Polder. Cette réflexion quantique est déterminante pour des expériences telles que GBAR, qui mesurera l'accélération d'un atome d'antihydrogène froid chutant vers une plaque de détection. Dans cette thèse, le potentiel de Casimir-Polder est calculé à partir des propriétés de diffusion électromagnétique de l'atome et de la surface. Il s'avère dépendre de la réponse diélectrique, de l'épaisseur et de la densité du milieu. Nous montrons que la réflexion sur ce potentiel est associée à une rupture de l'approximation semiclassique et qu'elle augmente pour des atomes lents et des potentiels faibles. Les transformations de Liouville relient des équations de Schrödinger avec des potentiels différents mais les mêmes amplitudes de diffusion. L'équivalence entre la réflexion quantique sur un puits de potentiel et l'effet tunnel à travers une barrière offre de nouvelles perspectives sur le problème. Nous discutons aussi des effets de la gravité sur le paquet d'onde atomique et de ses conséquences pour les expériences avec des atomes en chute libre. Associée à la réflexion quantique sur un miroir horizontal, la gravité permet de maintenir des particules dans des états à longue durée de vie aux applications prometteuses pour la métrologie. En particulier, nous proposons un système pour améliorer la précision de GBAR en réduisant la dispersion en vitesse des atomes d'antihydrogène. / Collisions between ultracold atoms and material surfaces are characterized by the reflection of the atomic matter wave from the attractive Casimir-Polder potential. This quantum reflection is particularly relevant to experiments such as GBAR, which will determine the gravitational acceleration of a cold antihydrogen atom by timing its fall onto a detection plate. In this thesis, the Casimir-Polder potential is computed from the electromagnetic scattering properties of the atom and surface and it is found to depend notably on the dielectric response, thickness and density of the medium. We show that reflection on this potential is associated with a breakdown of the semiclassical approximation and that it is enhanced for slow atoms and weak potentials. Liouville transformations relate Schrödinger equations with different potential landscapes but identical scattering properties. We gain new insights on the problem of quantum reflection on a potential well by mapping it onto an equivalent problem of tunneling through a wall. We also discuss the effect of gravity on the atomic wavepacket and its implications for free fall experiments with atoms. When combined with quantum reflection from a horizontal mirror, gravity can be used to trap particles in long lived states with promising applications for metrology. In particular, we suggest a scheme to improve the precision of the GBAR experiment by reducing the velocity dispersion of the falling atoms.
|
264 |
Limite de champ moyen pour des modèles discrets et équation de Schrödinger non linéaire discrète / Mean field limit for discrete models and nonlinear discrete Schrödinger equationPawilowski, Boris 11 December 2015 (has links)
Dans une série de travaux Zied Ammari et Francis Nier ont développé des méthodes pour étudier la dynamique de champ moyen bosonique pour des états quantiques généraux pouvant présenter des corrélations. Ils ont obtenu des formules pour décrire la dynamique des corrélations, ou plus généralement des matrices densité réduites d'ordre arbitraire. Cette thématique a été largement développée ces dernières années. Norbert Mauser en a été un des contributeurs, ainsi que sur la notion de mesure de Wigner qui est la clé de l'analyse développée par Z. Ammari et F. Nier. En général, il est admis que l'asymptotique de champ moyen est une bonne approximation du problème à N particules quand N dépasse la dizaine. Cela concerne l'asymptotique de la matrice densité réduite à une particule qui ne décrit pas la dynamique des corrélations. Un objectif est de tester la validité de la dynamique de champ moyen pour les matrices densité réduites à 2-particules. Pour des tests numériques, les modèles discrets qui n'ont pas été vraiment traités en détail dans les travaux précédents de Z. Ammari et F. Nier semblent bien adaptés. La thèse comprendra donc plusieurs étapes: adapter les résultats précédents de Z. Ammari et F. Nier à des modèles discrets , développer des méthodes numériques pour des systèmes simples mais pertinents, permettant de valider l'approximation de champ moyen et les formules pour la dynamique des corrélations. Au niveau numérique, on utilise des schémas numériques symplectiques, développés spécifiquement ces dernières années pour la discrétisation des équations hamiltoniennes. Une dernière étape concerne la combinaison des deux asymptotiques, champ moyen et approximation des modèles continus par les modèles discrets. / In a serie of works Z. Ammari and F. Nier developed methods to study the dynamics of bosonic mean field for general quantum states which can present correlations. They obtained formulas to describe the dynamics of the correlations, or more generally reduced density matrices with an arbitrary order. This topic was widely developed these last years. N.J. Mauser was one of contributors, as well as on the notion of Wigner measure which is the key of the analysis developed by Z. Ammari and F. Nier. Generally, the mean field asymptotic is admitted is a good approximation of the N-body problem when N exceed about ten. It concerns the asymptotics of the reduced density matrices for one particle which does not describe the dynamics of the correlations. An objective is to test the validity of the mean field dynamics for reduced density matrices for 2 particles. For numerical tests, the discrete models which were not really handled in detail in the previous works of Z. Ammari and F. Nier seem adapted well. The thesis will thus include several steps: adapt the previous results from Z. Ammari and F. Nier to discrete models , develop numerical methods, for simple but relevant systems, allowing to validate the approximation of mean field and the formulas for the dynamics of the correlations. About numerics, symplectic numerical scheme are used, developed specifically these last years for the discretization of the hamiltonian equations. A last possible step concerns the combination of both asymptotics, that is mean field and approximation of the continuous models by the discrete models.
|
265 |
Existence et stabilité de solutions fortes en théorie cinétique des gaz / Existence and stability of strong solutions in kinetic theoryTristani, Isabelle 22 June 2015 (has links)
Cette thèse est centrée sur l’étude d’équations issues de la théorie cinétique des gaz. Dans tous les problèmes qui y sont explorés, une analyse des problèmes linéaires ou linéarisés associés est réalisée d’un point de vue spectral et du point de vue des semi-groupes. A cela s’ajoute une analyse de la stabilité non linéaire lorsque le modèle est non linéaire. Plus précisément, dans une première partie, nous nous intéressons aux équations de Fokker-Planck fractionnaire et Boltzmann sans cut-off homogène en espace et nous prouvons un retour vers l’équilibre des solutions de ces équations avec un taux exponentiel dans des espaces de type L1 à poids polynomial. Concernant l’équation de Landau inhomogène en espace, nous développons une théorie de Cauchy de solutions perturbatives dans des espaces de type L2 avec différents poids (polynomiaux ou exponentiels) et nous prouvons également la stabilité exponentielle de ces solutions.Nous démontrons ensuite pour l’équation de Boltzmann inélastique inhomogène avec terme diffusif le même type de résultat dans des espaces L1 à poids polynomial dans un régime de faible inélasticité. Pour finir, nous étudions dans un cadre général et uniforme des modèles qui convergent vers l’équation de Fokker-Planck du point de vue de l’analyse spectrale et des semi-groupes. / The topic of this thesis is the study of models coming from kinetic theory. In all the problems that are addressed, the associated linear or linearized problem is analyzed from a spectral point of view and from the point of view of semigroups. Tothat, we add the study of the nonlinear stability when the equation is nonlinear. More precisely, to begin with, we treat the problem of trend to equilibrium for the fractional Fokker-Planck and Boltzmann without cut-off equations, proving an exponential decay to equilibrium in spaces of type L1 with polynomial weights. Concerning the inhomogeneous Landau equation, we develop a Cauchy theory of perturbative solutions in spaces of type L2 with various weights such as polynomial and exponential weights and we also prove the exponential stability of these solutions. Then, we prove similar results for the inhomogeneous inelastic diffusively driven Boltzmann equation in a small inelasticity regime in L1 spaces with polynomial weights. Finally, we study in the same and uniform framework from the spectral analysis point of view with a semigroup approach several Fokker-Planck equations which converge towards the classical one.
|
266 |
Modèles mathématiques de la dynamique des populations en environnement déterministe et stochastique / Mathematical population dynamics models in deterministic and stochastic environmentsNguyen, Trong Hieu 13 October 2014 (has links)
Dans ce travail de thèse, nous étudions des modèles mathématiques de la dynamique des populations en environnements déterministe et stochastique. Pour les environnements déterministes, nous considérons trois modèles. Le premier est un modèle intra-guilde prenant en compte des effets d'un environnement spatial hétérogène avec une migration rapide des individus entre les différents sites. Le deuxième est un modèle de pêche dans une zone constituée d’une aire marine protégée où la pêche est interdite et d’une zone où la population de poissons est pêchée. Enfin le troisième est un modèle prédateur-proie considérant une proie et deux prédateurs avec des réponses fonctionnelles de Beddington-DeAngelis. Pour les environnements stochastiques, nous étudions un modèle épidémique SIRS et un modèle prédateur-proie en prenant en compte un bruit télégraphique. Nous étudions le comportement dynamique de ces modèles et nous recherchons les conditions de maintien ou de disparition des espèces modélisées. / In this thesis, we consider mathematical population dynamics models in deterministic and stochastic environments. For deterministic environments, we study three models: an intraguild model with the effects of spatial heterogeneous environment and fast migration of individuals; a fishery model with Marine Protected Area where fishing is prohibited and an area where the fish population is harvested; a predator-prey model which has one prey and two predators with Beddington-DeAngelis functional responses. For stochastic environments, we study SIRS epidemic model and predator-prey models under telegraph noise. We try to present the dynamical behavior of these models and show out the existence or vanishing of species in the models.
|
267 |
Équation de réaction-diffusion en milieux hétérogènes : persistence, propagation et effet de la géométrie / Reaction diffusion equation in heterogeneous media : persistance, propagation and effect of the geometryBouhours, Juliette 08 July 2014 (has links)
Dans cette thèse nous nous intéressons aux équations de réaction-diffusion et à leurs applications en sciences biologiques et médicales. Plus particulièrement on étudie l'existence ou la non-existence de phénomènes de propagation en milieux hétérogènes à travers l'existence d'ondes progressives ou plus généralement l'existence de fronts de transition généralisés. On obtient des résultats d'existence de phénomènes de propagation dans trois environnements différents. Dans un premier temps on étudie une équation de réaction-diffusion de type bistable dans un domaine extérieur. Cette équation modélise l'évolution de la densité d'une population soumise à un effet Allee fort dont le déplacement suit un processus de diffusion dans un environnement contenant un obstacle. On montre que lorsque l'obstacle satisfait certaines conditions de régularité et se rapproche d'un domaine étoilé ou directionnellement convexe alors la population envahit tout l'espace. On se questionne aussi sur les conditions optimales de régularité qui garantissent une invasion complète de la population. Dans un deuxième travail, nous considérons une équation de réaction-diffusion avec vitesse forcée, modélisant l'évolution de la densité d'une population quelconque qui se diffuse dans l'espace, soumise à un changement climatique défavorable. On montre que selon la vitesse du changement climatique la population s'adapte ou s'éteint. On montre aussi que la densité de population converge en temps long vers une onde progressive et donc se propage (si elle survit) selon un profile constant et à vitesse constante. Dans un second temps on étudie une équation de réaction-diffusion de type bistable dans des domaines cylindriques variés. Ces équations modélisent l'évolution d'une onde de dépolarisation dans le cerveau humain. On montre que l'onde est bloquée lorsque le domaine passe d'un cylindre très étroit à un cylindre de diamètre d'ordre 1 et on donne des conditions géométriques plus générales qui garantissent une propagation complète de l'onde dans le domaine. On étudie aussi ce problème d'un point de vue numérique et on montre que pour les cylindres courbés la courbure peut provoquer un blocage de l'onde pour certaines conditions aux bords. / In this thesis we are interested in reaction diffusion equations and their applications in biology and medical sciences. In particular we study the existence or non-existence of propagation phenomena in non homogeneous media through the existence of traveling waves or more generally the existence of transition fronts.First we study a bistable reaction diffusion equation in exterior domain modelling the evolution of the density of a population facing an obstacle. We prove that when the obstacle satisfies some regularity properties and is close to a star shaped or directionally convex domain then the population invades the entire domain. We also investigate the optimal regularity conditions that allow a complete invasion of the population. In a second work, we look at a reaction diffusion equation with forced speed, modelling the evolution of the density of a population facing an unfavourable climate change. We prove that depending on the speed of the climate change the population keeps track with the climate change or goes extinct. We also prove that the population, when it survives, propagates with a constant profile at a constant speed at large time. Lastly we consider a bistable reaction diffusion equation in various cylindrical domains, modelling the evolution of a depolarisation wave in the brain. We prove that this wave is blocked when the domain goes from a thin channel to a cylinder, whose diameter is of order 1 and we give general conditions on the geometry of the domain that allow propagation. We also study this problem numerically and prove that for curved cylinders the curvature can block the wave for particular boundary conditions.
|
268 |
Exposants de Lyapunov d’opérateurs de Schrödinger ergodiques / Lyapunov exponents of ergodic Schrödinger operatorsMetzger, Florian 08 June 2017 (has links)
L'objectif de cette thèse est de traiter de deux aspects différents de la théorie de l'exposant de Lyapunov de cocycles de Schrödinger définis par une dynamique ergodique. Dans la première partie, on s'intéresse aux estimées de grandes déviations de type Bourgain & Goldstein pour des cocycles quasi-périodiques, puis pour ceux définis par le doublement de l'angle. Après avoir montré que seule une estimée par dessus sur une bande complexe est nécessaire pour avoir la minoration, on redémontre cette inégalité pour une dynamique quasi-périodique en utilisant des techniques de mouvement brownien en lien avec des fonctions sous-harmoniques. Ensuite on adapte la méthode au cas du doublement de l'angle pour lequel on prouve des estimées de grandes déviations sur les branches inverses de cette dynamique. Dans la deuxième partie sont étudiés des cocycles de Schrödinger dont la dynamique est une somme de dynamiques quasi-périodique et aléatoire. On démontre que, dans le régime perturbatif, les développements asymptotiques de l'exposant de Lyapunov attaché à ces cocycles sont similaires à ceux déjà démontrés dans le cas aléatoire par Figotin & Pastur ou Sadel & Schulz-Baldes. L'analyse se fait en fonction du caractère diophantien ou résonant de l'énergie par rapport à la fréquence diophantienne de la partie quasi-périodique du potentiel. / In this thesis we are interested in the Lyapunov exponent of ergodic Schrödinger cocycles. These cocycles occur in the analysis of solutions to the Schrödinger equation where the potential is defined with ergodic dynamics. We study two distinct aspects related to the the Lyapunov exponent for different kinds of dynamics. First we focus on a large deviation theorem for quasi-periodic cocycles and then for potentials defined by the doubling map. We prove that estimates of Bourgain & Goldstein type are granted if an upper estimate involved in the theorem is true on a strip of the complex plane. Then we establish a new technique to prove this upper bound in the quasi-periodic setting, based on subharmonic arguments suggested by Avila, Jitomirskaya & Sadel. We adapt afterwards the method to the doubling map and prove a large deviation theorem for the inverse branches of this dynamics. In the second part, we establish an asymptotic development similar to the results of Figotin & Pastur and Sadel & Schulz-Baldes for the Lyapunov exponent of Schrödinger cocycles at small coupling when the potential is a mixture of quasi-periodic and random. The analysis distinguishes the cases when the energy is either diophantine or resonant with respect to the frequency of the quasi-periodic part of the potential.
|
269 |
An axial polynomial expansion and acceleration of the characteristics method for the solution of the Neutron Transport Equation / Méthode accélérée aux caractéristiques pour la solution de l'équation du transport des neutrons, avec une approximation polynomiale axialeGraziano, Laurent 16 October 2018 (has links)
L'objectif de ce travail de thèse est le développement d'une approximation polynomiale axiale dans un solveur basé sur la Méthode des Caractéristiques. Le contexte, est celui de la solution stationnaire de l'équation de transport des neutrons pour des systèmes critiques, et l'implémentation pratique a été réalisée dans le solveur "Two/three Dimensional Transport" (TDT), faisant partie du projet APOLLO3®. Un solveur MOC pour des géométries en trois dimensions a été implémenté dans ce code pendant un projet de thèse antécédent, se basant sur une approximation constante par morceaux du flux et des sources des neutrons. Les développements présentés dans la suite représentent la continuation naturelle de ce travail. Les solveurs MOC en trois dimensions sont capables de produire des résultats précis pour des géométries complexes. Bien que précis, le coût computationnel associé à ce type de solveur est très important. Une représentation polynomiale en direction axiale du flux angulaire des neutrons a été utilisée pour réduire ce coût computationnel.Le travail réalisé pendant cette thèse peut être considéré comme divisé en trois parties: transport, accélération et autres. La première partie est constituée par l'implémentation de l'approximation polynomiale choisie dans les équations de transmission et de bilan typiques de la Méthode des Caractéristiques. Cette partie a aussi été caractérisée par le calcul d'une série de coefficients numériques qui se sont révélés nécessaires afin d'obtenir un algorithme stable. Pendant la deuxième partie, on a modifié et implémenté la solution des équations de la méthode d'accélération DPN. Cette méthode était déjà utilisée pour l'accélération et des itérations internes et externes dans TDT pour les solveurs deux et trois dimensionnels avec l'approximation des flux plat, quand ce travail a commencé. L'introduction d'une approximation polynomiale a demandé plusieurs développements numériques regardant la méthode d'accélération. Dans la dernière partie de ce travail on a recherché des solutions pour un mélange de différents problèmes liés aux premières deux parties. En premier lieux, on a eu à faire avec des instabilités numériques associées à une discrétisation spatiale ou angulaire pas suffisamment précise, soit pour la partie transport que pour la partie d'accélération. Ensuite, on a essayé d'utiliser différentes méthodes pour réduire l'empreinte mémoire des coefficients d'accélération. L'approche qu'on a finalement choisie se base sur une régression non-linéaire au sens des moindres carrés de la dépendance en fonction des sections efficaces typique de ces coefficients. L'approche standard consiste dans le stockage d'une série de coefficients pour chaque groupe d'énergie. La méthode de régression permet de remplacer cette information avec une série de coefficients calculés pendant la régression qui sont utilisés pour reconstruire les matrices d'accélération au cours des itérations. Cette procédure ajoute un certain coût computationnel à la méthode, mais nous pensons que la réduction de la mémoire rende ce surcoût acceptable.En conclusion, le travail réalisé a été concentré sur l'application d'une simple approximation polynomiale avec l'objectif de réduire le coût computationnel et l'empreinte mémoire associées à un solveur basée sur la Méthodes des Caractéristiques qui est utilisé pour calculer le flux neutroniques pour des géométries à trois dimensions extrudées. Même si cela ne constitue pas une amélioration radicale des performances, l'approximation d'ordre supérieur qu'on a introduit permet une réduction en termes de mémoire et de temps de calcul d'un facteur compris entre 2 et 5, selon le cas. Nous pensons que ces résultats constitueront une base fertile pour des futures améliorations. / The purpose of this PhD is the implementation of an axial polynomial approximation in a three-dimensional Method Of Characteristics (MOC) based solver. The context of the work is the solution of the steady state Neutron Transport Equation for critical systems, and the practical implementation has been realized in the Two/three Dimensional Transport (TDT) solver, as a part of the APOLLO3® project. A three-dimensional MOC solver for 3D extruded geometries has been implemented in this code during a previous PhD project, relying on a piecewise constant approximation for the neutrons fluxes and sources. The developments presented in the following represent the natural continuation of this work. Three-dimensional neutron transport MOC solvers are able to produce accurate results for complex geometries. However accurate, the computational cost associated to this kind of solvers is very important. An axial polynomial representation of the neutron angular fluxes has been used to lighten this computational burden.The work realized during this PhD can be considered divided in three major parts: transport, acceleration and others. The first part is constituted by the implementation of the chosen polynomial approximation in the transmission and balance equations typical of the Method Of Characteristics. This part was also characterized by the computation of a set of numerical coefficients which revealed to be necessary in order to obtain a stable algorithm. During the second part, we modified and implemented the solution of the equations of the DPN synthetic acceleration. This method was already used for the acceleration of both inners and outers iteration in TDT for the two and three dimensional solvers at the beginning of this work. The introduction of a polynomial approximation required several equations manipulations and associated numerical developments. In the last part of this work we have looked for the solutions of a mixture of different issues associated to the first two parts. Firstly, we had to deal with some numerical instabilities associated to a poor numerical spatial or angular discretization, both for the transport and for the acceleration methods. Secondly, we tried different methods to reduce the memory footprint of the acceleration coefficients. The approach that we have eventually chosen relies on a non-linear least square fitting of the cross sections dependence of such coefficients. The default approach consists in storing one set of coefficients per each energy group. The fit method allows replacing this information with a set of coefficients computed during the regression procedure that are used to re-construct the acceleration matrices on-the-fly. This procedure of course adds some computational cost to the method, but we believe that the reduction in terms of memory makes it worth it.In conclusion, the work realized has focused on applying a simple polynomial approximation in order to reduce the computational cost and memory footprint associated to a Method Of Characteristics solver used to compute the neutron fluxes in three dimensional extruded geometries. Even if this does not a constitute a radical improvement, the high order approximation that we have introduced allows a reduction in terms of memory and computational times of a factor between 2 and 5, depending on the case. We think that these results will constitute a fertile base for further improvements.
|
270 |
Stabilisation de quelques équations d’évolution du second ordrepar des lois de rétroaction / Stabilization of second order evolution equations with dynamical feedbacksAbbas, Zainab 02 October 2014 (has links)
Dans cette thèse, nous étudions la stabilisation de certaines équations d’évolution par des lois de rétroaction. Dans le premier chapitre nous étudions l’équation des ondes dans R avec conditions aux limites dynamiques appliquées sur une partie du bord et une condition de Dirichlet sur la partie restante. Nous fournissons des conditions suffisantes qui garantissent une stabilité polynomiale en utilisant une méthode qui combine une inégalité d’observabilité pour le problème non amorti associé avec des résultats de régularité du problème non amorti. L’optimalité de la décroissance est montrée dans certains cas à l’aide des résultats spectraux précis de l’opérateur associé. Dans le deuxième chapitre nous considérons le système sur un domaine de Rd, d ≥ 2. On trouve des conditions suffisantes qui permettent la stabilité forte. Ensuite, nous discutons de la stabilité non uniforme ainsi que de la stabilité polynomiale. L’approche en domaine fréquentiel nous permet d’établir une décroissance polynomiale sur des domaines pour lesquels l’équation des ondes avec l’amortissement standard est exponentiellement ou polynomialement stable. Dans le troisième chapitre nous considérons un cadre général d’équations d’évolution avec une dissipation dynamique. Sous une hypothèse de régularité, nous montrons que les propriétés d’observabilité pour le problème non amorti impliquent des estimations de décroissance pour le problème amorti. / In this thesis, we study the stabilization of some evolution equations by feedback laws. In the first chapter we study the wave equation in R with dynamical boundary control applied on a part of the boundary and a Dirichlet boundary condition on the remaining part. We furnish sufficient conditions that guarantee a polynomial stability proved using a method that combines an observability inequality for the associated undamped problem with regularity results of the solution of the undamped problem. In addition, the optimality of the decay is shown in some cases with the help of precise spectral results of the operator associated with the damped problem. Then in the second chapter we consider the system on a domain of Rd, d ≥ 2. In this case, the domain of the associated operator is not compactly embedded into the energy space. Nevertheless, we find sufficient conditions that give the strong stability. Then, we discuss the non uniform stability as well as the polynomial stability by two methods. The frequency domain approach allows us to establish a polynomial decay on some domains for which the wave equation with the standard damping is exponentially or polynomially stable. Finally, in the third chapter we consider a general framework of second order evolution equations with dynamical feedbacks. Under a regularity assumption we show that observability properties for the undamped problem imply decay estimates for the damped problem. We finally illustrate our general results by a variety of examples.
|
Page generated in 0.1022 seconds