• Refine Query
  • Source
  • Publication year
  • to
  • Language
  • 8
  • 8
  • 1
  • Tagged with
  • 17
  • 17
  • 9
  • 7
  • 5
  • 5
  • 4
  • 4
  • 4
  • 4
  • 4
  • 4
  • 4
  • 3
  • 3
  • About
  • The Global ETD Search service is a free service for researchers to find electronic theses and dissertations. This service is provided by the Networked Digital Library of Theses and Dissertations.
    Our metadata is collected from universities around the world. If you manage a university/consortium/country archive and want to be added, details can be found on the NDLTD website.
1

Vers un critère d'erreur énergétique pour le maillage adaptatif en régime magnétodynamique

Ladas, Dimitrios 07 May 2008 (has links) (PDF)
Dans le domaine de la simulation numérique en électromagnétisme, l'évaluation d'erreur en régime magnétodynamique n'a pas trouvée de réponse pleinement satisfaisante jusqu'à maintenant. Après avoir introduit la méthode des éléments finis à partir d'une approche thermodynamique de l'électromagnétisme, nous proposons un nouveau critère d'erreur énergétique, qui permet d'évaluer l'erreur dans les parties conductrices sièges de courants de Foucault. Les caractéristiques de ce critère en fonction de la formulation et de l'ordre des éléments sont détaillées. Nous présentons ensuite une procédure de maillage adaptatif reposant sur ce critère appliquée à un cas en régime magnéto-harmonique 2D. Une extension aux régimes transitoires est introduite. La convergence des grandeurs énergétiques globales valide la pertinence du critère d'erreur proposé.
2

Multi-scale simulation of automotive catalytic converters / Simulation multi-échelle de l'écoulement dans les systèmes de post-traitement des gaz d'échappement automobile

Ozhan, Cansu 21 November 2014 (has links)
L'utilisation croissante de véhicules au cours des dernières décennies a causé une augmentation dans les émissions de gaz d'échappement nocifs provoquant des problèmes de santé et d'environnement. Cette problématique a conduit les gouvernements à mettre en place des limites d'émissions plus strictes. Afin de respecter ces limites, il est nécessaire de développer des systèmes de post-traitement plus performants. Parmi plusieurs solutions possibles, l'analyse et l'optimisation de l'impact de l'écoulement sur les réactions chimiques est une approche importante pour résoudre ce problème complexe. Au point de vue expérimental, il est très difficile de mesurer le champ de vitesse détaillé et la distribution de la température et de la concentration. Les simulations numériques peuvent fournir des informations supplémentaires pour comprendre l'interaction entre la distribution de l'écoulement et l'efficacité des réactions. La simulation numérique de ces systèmes est très coûteuse due aux phénomènes physiques et chimiques complexes ayant lieu simultanément dans tout le système. Afin de diminuer ce coût, on peut développer des modèles physiques et chimiques simplifiés ainsi que des techniques numériques spécifiques pour simuler le système tout en réduisant le temps de calcul. Dans cette thèse, nous développons une approche générale pour modéliser et simuler le système de post-traitement des gaz d'échappement automobile contenant les aspects physiques et chimiques. L'approche présente combine des modèles simplifiés et des techniques numériques de multi-résolutions afin de capturer correctement les caractéristiques de l'écoulement dans le système. Alors que la méthode de raffinement de maillage adaptatif (AMR) est optimisée afin de minimiser l'effort de calcul dans les régions divergente et convergente, un modèle de sous-maille est développé pour décrire l'écoulement dans les micro-canaux du substrat catalytique placé entre la région divergente et la région convergente. La performance du modèle est validée par rapport aux résultats expérimentaux obtenus par Benjamin et al. (2002). Le couplage effectif de méthode AMR et le modèle de sous-maille permet de capturer les caractéristiques de l'écoulement dans le système avec un temps de calcul réduit de manière significative. L'impact de l'écoulement pulsé et de l'écoulement en expansion sur le taux de réaction dans un canal de monolithe est étudié numériquement. La méthode AMR est montrée pour capturer les couches mince de diffusion près de la paroi à l'interface solide-liquide. Sur la base des résultats numériques, nous proposons un modèle simplifié de transport capturant les effets de l'écoulement, la diffusion et la réaction catalytiques à la paroi. Le modèle simplifié de transport peut être directement appliqué dans la forme de modèle de sous-maille pour la description complète de tous les processus physiques et chimiques ayant lieu à l'intérieur du système. Grâce aux approches de modélisation numériques physiques et chimiques développées, il est possible de simuler le système en trois dimensions avec un temps de calcul raisonnable capturant encore la physique principale du problème. / The increasing use of passenger vehicles over the past decades has caused an increase in harmful exhaust gas emissions which give rise to environmental and health problems. This problematic has led governments to establish very stringent emission limits. The emission restrictions require more performing after-treatment systems. Among many other solutions, analysing and optimising the flow impact on the conversion efficiency is an important step towards the solution of the complex engineering problem. Detailed velocity, temperature and concentration distributions are very difficult to measure experimentally. Numerical simulations can provide additional information to understand the interaction of flow distribution and conversion efficiency. The simulation of these systems is computationally very expensive due to complex physical and chemical phenomena occurring simultaneously throughout the system. To overcome this cost, one can resort to some simplified physical and chemical models together with specific numerical techniques to simulate the system with a reduced computational time. In this dissertation, we develop a general approach to model and simulate the automotive catalytic converter system including all the physical and chemical processes. The present approach combines sub-grid models and numerical multi-resolution techniques in order to correctly capture the flow characteristics inside an automotive catalytic converter. While Adaptive Mesh Refinement (AMR) techniques are optimized in order to minimise the computational effort in the divergent and convergent regions, a sub-grid model is developed to describe the flow inside the catalytic substrate placed between the convergent and divergent regions. The performance of the sub-grid model is validated against the experimental results obtained by Benjamin et al. (2002). The effective coupling of AMR techniques and the sub-grid model allows to capture the flow features with significantly reduced computational time. The impact of pulsating and expansion flow on the conversion efficiency within a single monolith channel is investigated numerically. AMR techniques are shown to capture the small boundary layers near the wall at the solid-fluid interface. Based on the numerical results, we propose a simplified transport model that captures the effects of flow, diffusion and catalytic wall reactions. The simplified transport model can be directly applied as a sub-grid model for the complete description of all the physical and chemical processes taking place inside the system. The developed physical, chemical and numerical modelling approaches make the three-dimensional simulations possible with a reasonable computational time still capturing the main physics of the problem.
3

Simulation of continuous damage and fracture in metal-forming processes with 3D mesh adaptive methodology / Simulation numérique d’endommagement continu et de fissure dans les procédés de mise en forme de métal avec 3D maillage méthodologie adaptatif

Yang, Fangtao 10 November 2017 (has links)
Ces travaux s'inscrivent dans le cadre des recherches menées dans le cadre d'une collaboration entre le laboratoire Roberval de l'Université de Technologie de Compiègne et l'équipe dans le cadre du projet ANR-14-CE07-0035 LASMIS de l'Institut Charles Delaunay de l'Université de Technologie de Troyes. Nous présentons dans ces travaux une h-méthodologie adaptative tridimensionnelle des éléments finis afin de représenter l'initiation et la propagation des fissures dans des matériaux ductiles. Un modèle élasto-plastique couplé à l'endommagement isotrope proposé par l'équipe du LASMIS/UTT est utilisé. Les applications visées à terme concernent principalement la mise en forme des métaux. Dans ce contexte, une formulation Lagrangienne actualisée est employée et des remaillages fréquents s'avèrent essentiels afin d'une part d'éviter les fortes distorsions d'éléments dues aux grandes déformations plastiques et d'autre part de suivre les modifications de la topologie résultant de la création de fissures. La taille du nouveau maillage doit permettre à moindre coût représenter avec précision l'évolution des gradients des quantités physiques représentatives des phénomènes étudiées (plasticité, endommagement...). Nous proposons des indicateurs empiriques de taille d'éléments basés sur la déformation plastique ainsi que sur l'endommagement. Une courbe définie par morceau représente l'évolution de la taille d'élément suivant la sévérité de la plasticité et le cas échéant de l'endommagement. Les fissures sont représentées par une méthode de destruction d'éléments qui permet une description aisée de la géométrie de ces dernières et une gestion simplifiée de la fissuration sans nul besoin de critères additionnels. En revanche, pour permettre une description réaliste des fissures, ces dernières doivent être représentées par l'érosion des éléments de plus petite taille. Un solveur ABAQUS/Explicit® est utilisé avec des éléments tétraédriques quadratiques (C3D10M) évitant notamment les problèmes de verrouillage numérique survenant lors de l'analyse de structures en matériau compressible ou quasi-incompressible. Le contrôle de la plus petite taille du maillage est important dans un contexte explicite. De surcroît, pour les phénomènes adoucissant, la solution dépend de la taille de maille considérée alors comme un paramètre intrinsèque. Une étude nous a permis de constater que lorsque le maillage est suffisamment raffiné, les effets de la dépendance au maillage se réduisaient. Dans la littérature, les coûts de maillage ou de remaillage fréquents sont souvent considérés comme prohibitifs et de nombreux auteurs s'appuient sur cet argument pour introduire, avec succès certes, des méthodes alternatives qui limitent le coût des opérations de remaillage sans toutefois les éliminer (XFEM par exemple). Nos travaux montrent que le coût d'un remaillage local est négligeable par rapport au calcul. Compte tenu de la complexité de la géométrie et de la nécessité de raffiner le maillage, la seule alternative à ce jour est d'utiliser un mailleur en tétraèdres. La stratégie de remaillage local en tétraèdre s'appuie sur une méthode de bisection suivie si nécessaire d'une optimisation locale du maillage proposé par A. Rassineux en 2003. Le remaillage, même local, doit s'accompagner de procédures de transfert de champ des variables nodales et aux points d'intégration. Les variables nodales sont, comme le fait la plupart des auteurs, transférées en utilisant les fonctions de forme éléments finis. Le transfert de champ en 3D aux points de Gauss et les nombreux problèmes sous-jacents ont été relativement peu abordés dans la littérature. / This work is part of the research carried out in the framework of a collaboration between the Roberval laboratory of the Compiègne University of Technology and the team within the framework of the project ANR-14-CE07-0035 LASMIS of the Charles Delaunay Institute of Technology University of Troyes. In this work, we present a three-dimensional adaptive Pi-methodology of finite elements to represent the initiation and propagation of cracks in ductile materials. An elastoplastic model coupled with the isotropic damage proposed by the LASMIS / UTT team is used. The targeted applications will mainly concern the metal forming. In this context, an updated Lagrangian formulation is used and frequent remeshing is essential in order to avoid the strong distortion of elements due to large plastic deformations and to follow the modifications of the topology resulting in the creation of cracks. The size of the new mesh must allow at a lower cost to accurately represent the evolution of the gradients of the physical quantities representative of the studied phenomena (plasticity, damage ...). We propose empirical indicators of size of elements based on the plastic deformation as well as on the damage. A piecewise defined curve represents the evolution of the element size according to the severity of the plasticity and, if appropriate, the damage. The cracks are represented by a method of destruction of elements which allows an easy description of the geometry and a simplified treatment of the cracking without any need for additional criteria. On the other hand, to allow a realistic description of the cracks, the latter must be represented by erosion smaller elements. An ABAQUS / Explicit@ solver is used with quadratic tetrahedral elements (C3DIOM), avoiding in particular the problems of numerical locking occurring during the analysis of structures in compressible or quasi-incompressible material. The control of the smaller mesh size is important in an explicit context. In addition, for softening phenomena, the solution depends on the mesh size considered as an intrinsic parameter. A study has shown that when the mesh is sufficiently refined, the effects of mesh dependence are reduced. In the literature, the costs of frequent meshing or remeshing are often considered prohibitive and many authors rely on this argument to introduce, with success, alternative methods that limit the cost of remeshing operations without eliminating them ( XFEM for example). Our work shows that the cost of local remeshing is negligible compared to the calculation. Given the complexity of the geometry and the need to refine the mesh, the only alternative to date is to use a mesh in tetrahedra. The strategy of local remeshing tetrahedron is based on a bisection method followed if necessary by a local optimization of the grid proposed by A. Rassineux in 2003. The remeshing, even local, must be accompanied by field transfer procedures on both nodal variables and integration points. Node variables are, as most authors do, transferred using finite element shape functions. The 3D field transfer at Gauss points and the many underlying problems have been relatively untouched in the literature. The main difficulties to be solved in order to ensure the "quality" of the transfer concern the limitation of numerical diffusion, the lack of information near borders, the respect of boundary conditions, the equilibrium, the calculation costs, the filtering of the information points, crucial problems in 3D where the number of Gauss points used is several hundred. We propose a so-called "hybrid" method which consists, initially, in extrapolating the data at the Gauss points, in the nodes by diffuse interpolation and then in using the finite element form functions to obtain the value at the point considered.
4

Équilibrage dynamique de charge sur supercalculateur exaflopique appliqué à la dynamique moléculaire / Dynamic load balancing on exaflop supercomputer applied to molecular dynamics

Prat, Raphaël 09 October 2019 (has links)
Dans le contexte de la dynamique moléculaire classique appliquée à la physique de la matière condensée, les chercheurs du CEA étudient des phénomènes physiques à une échelle atomique. Pour cela, il est primordial d'optimiser continuellement les codes de dynamique moléculaire sur les dernières architectures de supercalculateurs massivement parallèles pour permettre aux physiciens d'exploiter la puissance de calcul pour reproduire numériquement des phénomènes physiques toujours plus complexes. Cependant, les codes de simulations doivent être adaptés afin d'équilibrer la répartition de la charge de calcul entre les cœurs d'un supercalculateur.Pour ce faire, dans cette thèse nous proposons d'incorporer la méthode de raffinement de maillage adaptatif dans le code de dynamique moléculaire ExaSTAMP. L'objectif est principalement d'optimiser la boucle de calcul effectuant le calcul des interactions entre particules grâce à des structures de données multi-threading et vectorisables. La structure permet également de réduire l'empreinte mémoire de la simulation. La conception de l’AMR est guidée par le besoin d'équilibrage de charge et d'adaptabilité soulevé par des ensembles de particules se déplaçant très rapidement au cours du temps.Les résultats de cette thèse montrent que l'utilisation d'une structure AMR dans ExaSTAMP permet d'améliorer les performances de celui-ci. L'AMR permet notamment de multiplier par 1.31 la vitesse d'exécution de la simulation d'un choc violent entraînant un micro-jet d'étain de 1 milliard 249 millions d'atomes sur 256 KNLs. De plus, l'AMR permet de réaliser des simulations qui jusqu'à présent n'étaient pas concevables comme l'impact d'une nano-goutte d'étain sur une surface solide avec plus 500 millions d'atomes. / In the context of classical molecular dynamics applied to condensed matter physics, CEA researchers are studying complex phenomena at the atomic scale. To do this, it is essential to continuously optimize the molecular dynamics codes of recent massively parallel supercomputers to enable physicists to exploit their capacity to numerically reproduce more and more complex physical phenomena. Nevertheless, simulation codes must be adapted to balance the load between the cores of supercomputers.To do this, in this thesis we propose to incorporate the Adaptive Mesh Refinement method into the ExaSTAMP molecular dynamics code. The main objective is to optimize the computation loop performing the calculation of particle interactions using multi-threaded and vectorizable data structures. The structure also reduces the memory footprint of the simulation. The design of the AMR is guided by the need for load balancing and adaptability raised by sets of particles moving dynamically over time.The results of this thesis show that using an AMR structure in ExaSTAMP improves its performance. In particular, the AMR makes it possible to execute 1.31 times faster than before the simulation of a violent shock causing a tin microjet of 1 billion 249 million atoms on 256 KNLs. In addition, simulations that were not conceivable so far can be carried out thanks to AMR, such as the impact of a tin nanodroplet on a solid surface with more than 500 million atoms.
5

Analyse a posteriori d'algorithmes itératifs pour des problèmes non linéaires. / A posteriori analyses of iterative algorithm for nonlinear problems.

Dakroub, Jad 07 October 2014 (has links)
La résolution numérique de n’importe quelle discrétisation d’équations aux dérivées partielles non linéaires requiert le plus souvent un algorithme itératif. En général, la discrétisation des équations aux dérivées partielles donne lieu à des systèmes de grandes dimensions. Comme la résolution des grands systèmes est très coûteuse en terme de temps de calcul, une question importante se pose: afin d’obtenir une solution approchée de bonne qualité, quand est-ce qu’il faut arrêter l’itération afin d’éviter les itérations inutiles ? L’objectif de cette thèse est alors d’appliquer, à différentes équations, une méthode qui nous permet de diminuer le nombre d’itérations de la résolution des systèmes en gardant toujours une bonne précision de la méthode numérique. En d’autres termes, notre but est d’appliquer une nouvelle méthode qui fournira un gain remarquable en terme de temps de calcul. Tout d’abord, nous appliquons cette méthode pour un problème non linéaire modèle. Nous effectuons l’analyse a priori et a posteriori de la discrétisation par éléments finis de ce problème et nous proposons par la suite deux algorithmes de résolution itérative correspondants. Nous calculons les estimations d’erreur a posteriori de nos algorithmes itératifs proposés et nous présentons ensuite quelques résultats d’expérience numériques afin de comparer ces deux algorithmes. Nous appliquerons de même cette approche pour les équations de Navier-Stokes. Nous proposons un schéma itératif et nous étudions la convergence et l’analyse a priori et a posteriori correspondantes. Finalement, nous présentons des simulations numériques montrant l’efficacité de notre méthode. / The numerical resolution of any discretization of nonlinear PDEs most often requires an iterative algorithm. In general, the discretization of partial differential equations leads to large systems. As the resolution of large systems is very costly in terms of computation time, an important question arises. To obtain an approximate solution of good quality, when is it necessary to stop the iteration in order to avoid unnecessary iterations? A posteriori error indicators have been studied in recent years owing to their remarkable capacity to enhance both speed and accuracy in computing. This thesis deals with a posteriori error estimation for the finite element discretization of nonlinear problems. Our purpose is to apply a new method that allows us to reduce the number of iterations of the resolution system while keeping a good accuracy of the numerical method. In other words, our goal is to apply a new method that provides a remarkable gain in computation time. For a given nonlinear equation we propose a finite element discretization relying on the Galerkin method. We solve the discrete problem using two iterative methods involving some kind of linearization. For each of them, there are actually two sources of error, namely discretization and linearization. Balancing these two errors can be very important, since it avoids performing an excessive number of iterations. Our results lead to the construction of computable upper indicators for the full error. Similarly, we apply this approach to the Navier-Stokes equations. Several numerical tests are provided to evaluate the efficiency of our indicators.
6

Adaptation de maillage anisotrope en trois dimensions applications aux simulations instationnaires en mécanique des fluides /

Alauzet, Frédéric Mohammadi, Bijan January 2003 (has links)
Thèse doctorat : Mathématiques appliquées : Montpellier 2 : 2003. / Thèse soutenue à l'Université des Sciences et Techniques du Languedoc. Bibliogr. p. 181-188.
7

Propagation d'une onde de choc en présence d'une barrière de protection

Eveillard, Sébastien 12 September 2013 (has links) (PDF)
Les travaux de thèse présentés dans ce mémoire s'inscrivent dans le cadre du projet ANR BARPPRO. Ce programme de recherche vise à étudier l'influence d'une barrière de protection face à une explosion en régime de détonation. L'objectif est d'établir des méthodes de calcul rapides de classement des zones d'effets pour aider les industriels au dimensionnement des barrières de protection. L'une à partir d'abaques, valable pour des configurations en géométrie 2D, sur des plages spécifiées de paramètres importants retenus, avec une précision de +/- 5%. L'autre à partir d'une méthode d'estimation rapide basée notamment sur les chemins déployés, valable en géométrie 2D et en géométrie 3D, mais dont la précision estimée est de +/- 30%. Afin d'y parvenir, l'étude s'appuie sur trois volets : expérimental, simulation numérique et analytique. La partie expérimentale étudie plusieurs géométries de barrière de protection à petites échelles pour la détonation d'une charge gazeuse (propane-oxygène à la stoechiométrie). Les configurations expérimentées servent à la validation de l'outil de simulation numérique constitué du solveur HERA et de la plateforme de calcul TERA 100. Des abaques d'aide au dimensionnement ont pu être réalisés à partir de résultats fournis par l'outil de simulation (3125 configurations de barrière de protection, TNT). L'étude des différents phénomènes physiques présents a également permis de mettre en place une méthode d'estimation rapide basée sur des relations géométriques, analytiques et empiriques. L'analyse de ces résultats a permis d'établir quelques recommandations dans le dimensionnement d'une barrière de protection. Les abaques et le programme d'estimation rapide permettent à un ingénieur de dimensionner rapidement une barrière de protection en fonction de la configuration du terrain et de la position de la zone à protéger en aval du merlon.
8

Direct Numerical Simulation of bubbles with Adaptive Mesh Refinement with Distributed Algorithms / Simulation numérique directe de bulles sur maillage adaptatif avec algorithmes distribués

Talpaert, Arthur 24 February 2017 (has links)
Ce travail de thèse présente l'implémentation de la simulation d'écoulements diphasiques dans des conditions de réacteurs nucléaires à caloporteur eau, à l'échelle de bulles individuelles. Pour ce faire, nous étudions plusieurs modèles d'écoulements thermohydrauliques et nous focalisons sur une technique de capture d'interface mince entre phases liquide et vapeur. Nous passons ainsi en revue quelques techniques possibles de maillage adaptatif (AMR) et nous fournissons des outils algorithmiques et informatiques adaptés à l'AMR par patchs dont l'objectif localement la précision dans des régions d'intérêt. Plus précisément, nous introduisons un algorithme de génération de patchs conçu dans l'optique du calcul parallèle équilibré. Cette approche nous permet de capturer finement des changements situés à l'interface, comme nous le montrons pour des cas tests d'advection ainsi que pour des modèles avec couplage hyperbolique-elliptique. Les calculs que nous présentons incluent également la simulation du système de Navier-Stokes incompressible qui modélise la déformation de l'interface entre deux fluides non-miscibles. / This PhD work presents the implementation of the simulation of two-phase flows in conditions of water-cooled nuclear reactors, at the scale of individual bubbles. To achieve that, we study several models for Thermal-Hydraulic flows and we focus on a technique for the capture of the thin interface between liquid and vapour phases. We thus review some possible techniques for Adaptive Mesh Refinement (AMR) and provide algorithmic and computational tools adapted to patch-based AMR, which aim is to locally improve the precision in regions of interest. More precisely, we introduce a patch-covering algorithm designed with balanced parallel computing in mind. This approach lets us finely capture changes located at the interface, as we show for advection test cases as well as for models with hyperbolic-elliptic coupling. The computations we present also include the simulation of the incompressible Navier-Stokes system, which models the shape changes of the interface between two non-miscible fluids.
9

Schémas numériques instationnaires pour des écoulements multiphasiques multiconstituants dans des bassins sédimentaires

Nadau, Lionel 22 September 2003 (has links) (PDF)
Un bassin sédimentaire est un milieu poreux de grande dimension (plusieurs dizaines de kilomètres de long et de large pour une profondeur d'environ cinq kilomètres) qui évolue au cours du temps par les effets de compaction et de sédimentation. Au cours de cette évolution, des hydrocarbures vont se former et s'écouler dans le bassin. On établit alors un modèle permettant de simuler cette évolution de bassin ainsi que la création, la migration et le piégeage des hydrocarbures dans des roches appelées roches magasins. Ces phénomènes se déroulant sur des centaines de millions d'années, on s'est attaché à étudier principalement une discrétisation temporelle de ces équations. On a ainsi mis en avant un raffinement local du pas de temps dont le principe est de recalculer la solution sur une zone jugée "mauvaise". A l'extérieure de cette zone, la solution est admissible. La difficulté vient de la détermination de la zone qui doit - être suffisamment "grande" pour avoir une bonne qualité de la solution, mais suffisamment "petite" pour obtenir un gain calcul. Les estimateurs a posteriori permettent de contourner cette difficulté. On a donc entrepris une étude théorique de ces estimateurs a posteriori dans le cas des équations linéaires elliptique et parabolique. Des simulations numériques montrent l'efficacité de ces estimateurs dans des cas académiques.
10

Méthodes d'assemblage rapide et de résolution itérative pour un solveur adaptatif en équations intégrales de frontières destiné à l'électromagnétisme

Haghi-Ashtiani, Bidjan 07 May 1998 (has links) (PDF)
Dans ce travail nous avons ouvert des voies nouvelles pour la méthode des équations intégrales de frontière appliquée à la résolution des problèmes tridimensionnels de l'électromagnétisme en basse fréquence. Elles permettent de réaliser un maillage adaptatif associé à des résolutions approchées, locales, du système intégral. Dans le premier chapitre nous discutons brièvement les différentes méthodes numériques appliquées à l'électromagnétisme. Le deuxième chapitre est consacré à la méthode des équations intégrales de frontière. Nous y justifions nos choix particuliers (type et ordre des éléments discrétisant les frontières) et décrivons les algorithmes réalisés, par exemple pour la prise en compte des symétries ou périodicités des systèmes étudiés. Dans le troisième chapitre, nos méthodes de calcul des intégrales élémentaires, combinant l'analytique et le numérique, sont décrites et validées. Dans le chapitre quatre, nous présentons un estimateur d'erreur basé sur l'écart entre la valeur interpolée du potentiel et sa valeur calculée par la méthode intégrale elle-même. Pour finir, nous utilisons dans le chapitre cinq cet estimateur d'erreur pour réaliser automatiquement un. maillage adapté à chaque problème traité, en autorisant aussi un maillage « non-conforme ». Une méthode de résolution locale permet un gain important en temps de calcul pendant cette phase de maillage adaptatif. Le maillage final est séparé en domaines, ce qui permet d'améliorer le résultat par itérations entre ces domaines, et sans recourir à une résolution globale, très coûteuse pour des problèmes de grande taille.

Page generated in 0.0829 seconds