Reproduction de : Thèse de doctorat : Mathématiques : Faculté des sciences de Paris : 1868. / Titre provenant de la page de titre du document numérisé.

Ce mémoire contient quelques résultats sur l'intégration numérique. Ils sont liés à la célèbre formule de quadrature de K. F. Gauss. Une généralisation très intéressante de la formule de Gauss a été obtenue par P. Turán. Elle est contenue dans son article publié en 1948, seulement quelques années après la seconde guerre mondiale. Étant données les circonstances défavorables dans lesquelles il se trouvait à l'époque, l'auteur (Turán) a laissé beaucoup de détails à remplir par le lecteur. Par ailleurs, l'article de Turán a inspiré une multitude de recherches; sa formule a été étendue de di érentes manières et plusieurs articles ont été publiés sur ce sujet. Toutefois, il n'existe aucun livre ni article qui contiennent un compte-rendu détaillé des résultats de base, relatifs à la formule de Turán. Je voudrais donc que mon mémoire comporte su samment de détails qui puissent éclairer le lecteur tout en présentant un exposé de ce qui a été fait sur ce sujet. Voici comment nous avons organisé le contenu de ce mémoire. 1-a. La formule de Gauss originale pour les polynômes - L'énoncé ainsi qu'une preuve. 1-b. Le point de vue de Turán - Compte-rendu détaillé des résultats de son article. 2-a. Une formule pour les polynômes trigonométriques analogue à celle de Gauss. 2-b. Une formule pour les polynômes trigonométriques analogue à celle de Turán. 3-a. Deux formules pour les fonctions entières de type exponentiel, analogues à celle de Gauss pour les polynômes. 3-b. Une formule pour les fonctions entières de type exponentiel, analogue à celle de Turán. 4-a. Annexe A - Notions de base sur les polynômes de Legendre. 4-b. Annexe B - Interpolation polynomiale. 4-c. Annexe C - Notions de base sur les fonctions entières de type exponentiel. 4-d. Annexe D - L'article de P. Turán. / This mémoire contains some results about numerical integration. They are related to the famous quadrature formula of K. F. Gauss. A very interesting generalization of the formula of Gauss was obtained by P.Turán. It is contained in a paper that was published in 1948, only a few years after the second world war. Due to adverse circunstances he was in at the time, the author (Turán) left many details for the reader to fill in. Otherwise, the article of Turán inspired a multitude of research, and his formula has been extended in many ways and several papers have been written on this subject. However, there is no single book or paper where one can nd a clear and comprehensive account of the basic results pertaining to Turán's formula. Thus, I would like my Master's mémoire to contain enough details that can enlighten the reader and present an exposition of much that has been done on this subject. Here is how we have arranged the contents of the mémoire. 1-a. The original formula of Gauss for polynomials - statement along with a proof. 1-b. Turán's point of view - detailed account of the results contained in his paper. 2-a. A formula for trigonometric polynomials analogous to that of Gauss. 2-b. A formula for trigonometric polynomials analogous to that of Turán. 3-a. Two formulae for entire functions of exponential type, analogous to the one of Gauss for polynomials. 3-b. A formula for entire functions of exponential type, analogous to that of Turán. 4-a. Annexe A - Basic facts about Legendre polynomials. 4-b. Annexe B - Polynomial interpolation. 4-c. Annexe C - Basic facts about entire functions of exponential type. 4-d. Annexe D - Paper of P. Turán.

Ce mémoire contient quelques résultats sur l'intégration numérique. Ils sont liés à la célèbre formule de quadrature de K. F. Gauss. Une généralisation très intéressante de la formule de Gauss a été obtenue par P. Turán. Elle est contenue dans son article publié en 1948, seulement quelques années après la seconde guerre mondiale. Étant données les circonstances défavorables dans lesquelles il se trouvait à l'époque, l'auteur (Turán) a laissé beaucoup de détails à remplir par le lecteur. Par ailleurs, l'article de Turán a inspiré une multitude de recherches; sa formule a été étendue de di érentes manières et plusieurs articles ont été publiés sur ce sujet. Toutefois, il n'existe aucun livre ni article qui contiennent un compte-rendu détaillé des résultats de base, relatifs à la formule de Turán. Je voudrais donc que mon mémoire comporte su samment de détails qui puissent éclairer le lecteur tout en présentant un exposé de ce qui a été fait sur ce sujet. Voici comment nous avons organisé le contenu de ce mémoire. 1-a. La formule de Gauss originale pour les polynômes - L'énoncé ainsi qu'une preuve. 1-b. Le point de vue de Turán - Compte-rendu détaillé des résultats de son article. 2-a. Une formule pour les polynômes trigonométriques analogue à celle de Gauss. 2-b. Une formule pour les polynômes trigonométriques analogue à celle de Turán. 3-a. Deux formules pour les fonctions entières de type exponentiel, analogues à celle de Gauss pour les polynômes. 3-b. Une formule pour les fonctions entières de type exponentiel, analogue à celle de Turán. 4-a. Annexe A - Notions de base sur les polynômes de Legendre. 4-b. Annexe B - Interpolation polynomiale. 4-c. Annexe C - Notions de base sur les fonctions entières de type exponentiel. 4-d. Annexe D - L'article de P. Turán. / This mémoire contains some results about numerical integration. They are related to the famous quadrature formula of K. F. Gauss. A very interesting generalization of the formula of Gauss was obtained by P.Turán. It is contained in a paper that was published in 1948, only a few years after the second world war. Due to adverse circunstances he was in at the time, the author (Turán) left many details for the reader to fill in. Otherwise, the article of Turán inspired a multitude of research, and his formula has been extended in many ways and several papers have been written on this subject. However, there is no single book or paper where one can nd a clear and comprehensive account of the basic results pertaining to Turán's formula. Thus, I would like my Master's mémoire to contain enough details that can enlighten the reader and present an exposition of much that has been done on this subject. Here is how we have arranged the contents of the mémoire. 1-a. The original formula of Gauss for polynomials - statement along with a proof. 1-b. Turán's point of view - detailed account of the results contained in his paper. 2-a. A formula for trigonometric polynomials analogous to that of Gauss. 2-b. A formula for trigonometric polynomials analogous to that of Turán. 3-a. Two formulae for entire functions of exponential type, analogous to the one of Gauss for polynomials. 3-b. A formula for entire functions of exponential type, analogous to that of Turán. 4-a. Annexe A - Basic facts about Legendre polynomials. 4-b. Annexe B - Polynomial interpolation. 4-c. Annexe C - Basic facts about entire functions of exponential type. 4-d. Annexe D - Paper of P. Turán.

L’enrichissement des éléments finis standard est un outil performant pour améliorer la qualité d’approximation. L’idée principale de cette approche est d’ajouter aux fonctions de base un ensemble de fonctions censées améliorer la qualité des solutions approchées. Le choix de ces dernières est crucial et est en grande partie basé sur la connaissance a priori de quelques informations telles que les caractéristiques de la solution, de la géométrie du problème à résoudre, etc. L’efficacité de cette approche pour résoudre une équation aux dérivées partielles dans un maillage fixe, sans avoir recours au raffinement, a été prouvée dans de nombreuses applications dans la littérature. La clé de son succès repose principalement sur le bon choix des fonctions de base et plus particulièrement celui des fonctions d’enrichissement. Une question importante se pose alors : quelles conditions faut-il imposer sur les fonctions d’enrichissement afin qu’elles génèrent des éléments finis bien définis ?Dans cette thèse sont abordés différents aspects d’une approche générale d’enrichissement d’éléments finis. Notre première contribution porte principalement sur l’enrichissement de l’élément fini du type Q_1. Par contre, notre seconde contribution, certainement la plus importante, met l’accent sur une approche plus générale pour enrichir n’importe quel élément fini qu’il soit P_k, Q_k ou autres, conformes ou non conformes. Cette approche a conduit à l’obtention des versions enrichies de l’élément de Han, l’élément de Rannacher-Turek et l’élément de Wilson, qui font maintenant partie des codes d’éléments finis les plus couramment utilisés en milieu industriel. Pour établir ces extensions, nous avons eu recours à l’élaboration de nouvelles formules de quadrature multidimensionnelles appropriées généralisant les formules classiques bien connues en dimension 1, dites du “point milieu,” des “trapèzes” et de leurs versions perturbées, ainsi que la formule de Simpson. Elles peuvent être vues comme des extensions naturelles de ces formules en dimension supérieure. Ces dernières, en plus de leurs tests numériques implémentés sous MATLAB, version R2016a, ont fait l’objet de notre troisième contribution. Nous mettons particulièrement l’accent sur la détermination explicite des meilleures constantes possibles apparaissant dans les estimations d’erreur pour ces formules d’intégration. Enfin, dans la quatrième contribution nous testons notre approche pour résoudre numériquement le problème d’élasticité linéaire à l’aide d’un maillage rectangulaire. Nous effectuons l’analyse numérique aussi bien l’analyse de l’erreur d’approximation et résultats de convergence que l’analyse de l’erreur de consistance. Nous montrons également comment cette dernière peut être établie à n’importe quel ordre, généralisant ainsi certains travaux menés dans le domaine. Nous réalisons la mise en œuvre de la méthode et donnons quelques résultats numériques établis à l’aide de la bibliothèque libre d’éléments finis GetFEM++, version 5.0. Le but principal de cette partie sert aussi bien à la validation de nos résultats théoriques, qu’à montrer comment notre approche permet d’élargir la gamme de choix des fonctions d’enrichissement. En outre, elle permet de montrer comment cette large gamme de choix peut aider à avoir des solutions optimales et également à améliorer la validité et la qualité de l’espace d’approximation enrichie. / The enrichment of standard finite elements is a powerful tool to improve the quality of approximation. The main idea of this approach is to incorporate some additional functions on the set of basis functions. These latter are requested to improve the accuracy of the approximate solution. Their best choice is crucial and is based on the knowledge of some a priori information, such as the characteristics of the solution, the geometry of the problem to be solved, etc. The efficiency of such an approach for finding numerical solutions of partial differential equations using a fixed mesh, without recourse to refinement, was proved in numerous applications in the literature. However, the key to its success lies mainly on the best choice of the basis functions, and more particularly those of enrichment functions.An important question then arises: How to suitably choose them, in such a way that they generate a well-defined finite element ?In this thesis, we present a general approach that enables an enrichment of the finite element approximation. This was the subject of our first contribution, which was devoted to the enrichment of the classical Q_1 element, as a first step. As a second step, in our second contribution, we have developed a more general framework for enriching any finite element either P_k, Q_k or others, conforming or nonconforming. As an illustration of how to use this framework to build new enriched finite elements, we have introduced the extensions of some well-known nonconforming finite elements, notably, Han element, Rannacher-Turek element and Wilson element, which are now part of the main code of finite element methods. To establish these extensions, we have introduced a new family of multivariate versions of the classical trapezoidal, midpoint and Simpson rules. These latter, in addition to their numerical tests under MATLAB, version R2016a, have been the subject of our third contribution. They may be viewed as an extension of the well-known trapezoidal, midpoint and Simpson’s one-dimensional rules to higher dimensions. We particularly pay attention to the explicit expressions of the best possible constants appearing in the error estimates for these cubatute formulas. Finally, in the fourth contribution we apply our approach to numerically solving the linear elasticity problem based on a rectangular mesh. We carry out the numerical analysis of the approximation error and also for the consistency error, and show how the latter can be established to any order. This constitutes a generalization of some work already done in the field. In addition to our theoretical results, we have also made some numerical tests, which were achieved by using the GetFEM++ library, version 5.0. The aim of this contribution was not only to confirm our theoretical predictions, but also to show how the new developed framework allows us to expand the range of choices of enrichment functions. Furthermore, we have shown how this wide choices range can help us to improve some approximation properties and to get the optimal solutions for the particular problem of elasticity.

Le sujet de ce travail de thèse est l’application de la méthode de caractéristiques longues (MOC) pour résoudre l’équation du transport des neutrons pour des géométries à trois dimensions extrudées. Les avantages du MOC sont sa précision et son adaptabilité, le point faible était la quantité de ressources de calcul requises. Ce problème est même plus important pour des géométries à trois dimensions ou le nombre d’inconnues du problème est de l’ordre de la centaine de millions pour des calculs d’assemblage.La première partie de la recherche a été dédiée au développement des techniques optimisées pour le traçage et la reconstruction à-la-volé des trajectoires. Ces méthodes profitent des régularités des géométries extrudées et ont permis une forte réduction de l’empreinte mémoire et une réduction des temps de calcul. La convergence du schéma itératif a été accélérée par un opérateur de transport dégradé (DPN) qui est utilisé pour initialiser les inconnues de l’algorithme itératif and pour la solution du problème synthétique au cours des itérations MOC. Les algorithmes pour la construction et la solution des opérateurs MOC et DPN ont été accélérés en utilisant des méthodes de parallélisation à mémoire partagée qui sont le plus adaptés pour des machines de bureau et pour des clusters de calcul. Une partie importante de cette recherche a été dédiée à l’implémentation des méthodes d’équilibrage la charge pour améliorer l’efficacité du parallélisme. La convergence des formules de quadrature pour des cas 3D extrudé a aussi été explorée. Certaines formules profitent de couts négligeables du traitement des directions azimutales et de la direction verticale pour accélérer l’algorithme. La validation de l’algorithme du MOC a été faite par des comparaisons avec une solution de référence calculée par un solveur Monte Carlo avec traitement continu de l’énergie. Pour cette comparaison on propose un couplage entre le MOC et la méthode des Sous-Groupes pour prendre en compte les effets des résonances des sections efficaces. Le calcul complet d’un assemblage de réacteur rapide avec interface fertile/fissile nécessite 2 heures d’exécution avec des erreurs de quelque pcm par rapport à la solution de référence.On propose aussi une approximation d’ordre supérieur du MOC basée sur une expansion axiale polynomiale du flux dans chaque maille. Cette méthode permet une réduction du nombre de mailles (et d’inconnues) tout en gardant la même précision.Toutes les méthodes développées dans ce travail de thèse ont été implémentées dans la version APOLLO3 du solveur de transport TDT. / The topic of our research is the application of the Method of Long Characteristics (MOC) to solve the Neutron Transport Equation in three-dimensional axial geometries. The strength of the MOC is in its precision and versatility. As a drawback, it requires a large amount of computational resources. This problem is even more severe in three-dimensional geometries, for which unknowns reach the order of tens of billions for assembly-level calculations.The first part of the research has dealt with the development of optimized tracking and reconstruction techniques which take advantage of the regularities of three-dimensional axial geometries. These methods have allowed a strong reduction of the memory requirements and a reduction of the execution time of the MOC calculation.The convergence of the iterative scheme has been accelerated with a lower-order transport operator (DPN) which is used for the initialization of the solution and for solving the synthetic problem during MOC iterations.The algorithms for the construction and solution of the MOC and DPN operators have been accelerated by using shared-memory parallel paradigms which are more suitable for standard desktop working stations. An important part of this research has been devoted to the implementation of scheduling techniques to improve the parallel efficiency.The convergence of the angular quadrature formula for three-dimensional cases is also studied. Some of these formulas take advantage of the reduced computational costs of the treatment of planar directions and the vertical direction to speed up the algorithm.The verification of the MOC solver has been done by comparing results with continuous-in-energy Monte Carlo calculations. For this purpose a coupling of the 3D MOC solver with the Subgroup method is proposed to take into account the effects of cross sections resonances. The full calculation of a FBR assembly requires about 2 hours of execution time with differences of few PCM with respect to the reference results.We also propose a higher order scheme of the MOC solver based on an axial polynomial expansion of the unknown within each mesh. This method allows the reduction of the meshes (and unknowns) by keeping the same precision.All the methods developed in this thesis have been implemented in the APOLLO3 version of the neutron transport solver TDT.