Author: Edouard Demaldent
Publisher:
ISBN:
Category :
Languages : fr
Pages : 218
Get Book Here
Book Description
Cette thèse s'inscrit dans le domaine de la simulation numérique, et concerne l'étude des phénomènes de diffraction électromagnétique en régime harmonique. Nous nous intéressons plus particulièrement aux méthodes de représentation intégrale et aux simulations qui nécessitent l'usage d'un solveur direct. Leur domaine d'application est rapidement restreint avec les schémas d'approximation classiques, car ceux-ci requièrent un grand nombre d'inconnues pour obtenir un résultat précis. Pour remédier à ce problème, nous nous proposons d'adapter la méthode des éléments finis spectraux aux équations intégrales de l' électromagnétisme, puis au couplage intégro-différentiel. Notre approche préserve la conformité de l'espace d'approximation dans Hdiv(dans Hdiv-Hrotpour le couplage), et découple le temps d'assemblage de l'ordre d'approximation. Elle autorise ainsi une montée en ordre significative qui résulte en une réduction spectaculaire du nombre d'inconnues et des coûts de calcul, tout en assurant la précision du résultat. Une autre originalité de notre étude réside dans le développement d'éléments finis hexaédriques d'ordre anisotrope, pour traiter des obstacles métalliques recouverts d'une fine couche de matériau.
Author: Sophie Depeyre
Publisher:
ISBN:
Category :
Languages : fr
Pages : 238
Get Book Here
Book Description
NOUS NOUS SOMMES INTERESSES A LA CONSTRUCTION ET A L'ETUDE D'UNE CLASSE DE SCHEMAS D'ORDRE TROIS OU QUATRE EN TEMPS ET EN ESPACE, BASES SUR DES FORMULATIONS -SHEMAS DE TYPE VOLUMES FINIS OU ELEMENTS FINIS, POUR DES MAILLAGES BIDIMENSIONNELS EN RECTANGLES OU EN TRIANGLES. NOUS CONSIDERONS DANS UN PREMIER TEMPS DES PROBLEMES HYPERBOLIQUES LINEAIRES, COMME L'EQUATION D'ADVECTION ET LE SYSTEME DE MAXWELL. UNE ETUDE DE STABILITE ET DE PRECISION, A L'AIDE DES EQUATIONS EQUIVALENTES A ETE PRESENTEE, AFIN DE COMPARER LES SCHEMAS ET DE RETENIR LES PLUS PRECIS. EN PARTICULIER, POUR LE SYSTEME DE MAXWELL, UNE CONDITION NECESSAIRE ET SUFFISANTE DE STABILITE A ETE DEMONTREE POUR LE SCHEMA DECENTRE D'ORDRE UN, SUR UN MAILLAGE EN RECTANGLES. NOUS AVONS AUSSI PROPOSE UNE NOUVELLE FORMULATION DU SYSTEME DE MAXWELL, EN RAJOUTANT UN TERME DE VISCOSITE DANS LES EQUATIONS, AFIN QUE NOS SCHEMAS PRENNENT MIEUX EN COMPTE LES RELATIONS DE DIVERGENCE. UNE ETUDE DE STABILITE A PERMIS DE DETERMINER LE PARAMETRE DE VISCOSITE N'INTRODUISANT AUCUNE CONTRAINTE SUPPLEMENTAIRE SUR LE PAS DE TEMPS, ET NOUS AVONS MONTRE A L'AIDE DE RESULTATS NUMERIQUES, POURQUOI LA NOUVELLE FORMULATION ETAIT MEILLEURE. DANS LA DEUXIEME PARTIE, NOUS NOUS SOMMES INTERESSES A DES MODELES HYPERBOLIQUES NON LINEAIRES, COMME LES EQUATIONS D'EULER. NOUS AVONS CHERCHE A CONSTRUIRE DES LIMITEURS D'ORDRE ELEVE AFIN DE RENDRE NOS SCHEMAS POSITIFS. EN PARTICULIER, NOUS AVONS PRESENTE UN NOUVEAU LIMITEUR D'ORDRE TROIS, QUI S'EST AVERE STABLE ET ROBUSTE, POUR DES CALCULS DE TUBE A CHOC ET D'ECOULEMENTS TRANSSONIQUES STATIONNAIRES. NOUS AVONS FINALEMENT CONSIDERE UN MODELE EULERIEN D'ECOULEMENT DIPHASIQUE, HYPERBOLIQUE ET CONSERVATIF, COMPORTANT UN TERME SOURCE RAIDE. LA METHODE CLASSIQUE D'INTEGRATION EN TEMPS EST UNE METHODE DE PAS FRACTIONNAIRES ; TOUTEFOIS, ELLE COMPORTE PLUSIEURS FAIBLESSES, ET NOUS AVONS PROPOSE UNE METHODE COUPLEE, QUI S'AVERE PLUS PRECISE LORSQUE LE RAYON DES PARTICULES DEVIENT PETIT.
Author: Marcella Bonazzoli
Publisher:
ISBN:
Category :
Languages : en
Pages : 0
Get Book Here
Book Description
The time-harmonic Maxwell's equations present several difficulties when the frequency is large, such as the sign-indefiniteness of the variational formulation, the pollution effect and the problematic construction of iterative solvers. We propose a precise and efficient solution strategy that couples high order finite element (FE) discretizations with domain decomposition (DD) preconditioners. High order FE methods make it possible for a given precision to reduce significantly the number of unknowns of the linear system to be solved. DD methods are then used as preconditioners for the iterative solver: the problem defined on the global domain is decomposed into smaller problems on subdomains, which can be solved concurrently and using robust direct solvers. The design, implementation and analysis of both these methods are particularly challenging for Maxwell's equations. FEs suited for the approximation of the electric field are the curl-conforming or edge finite elements. Here, we revisit the classical degrees of freedom (dofs) defined by Nédélec to obtain a new more friendly expression in terms of the chosen high order basis functions. Moreover, we propose a general technique to restore duality between dofs and basis functions. We explicitly describe an implementation strategy, which we embedded in the open source language FreeFem++. Then we focus on the preconditioning of the linear system, starting with a numerical validation of a one-level overlapping Schwarz preconditioner, with impedance transmission conditions between subdomains. Finally, we investigate how two-level preconditioners recently analyzed for the Helmholtz equation work in the Maxwell case, both from the theoretical and numerical points of view. We apply these methods to the large scale problem arising from the modeling of a microwave imaging system, for the detection and monitoring of brain strokes. In this application accuracy and computing speed are indeed of paramount importance.
Author: Marc Duruflé
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
Dans cette thèse, nous nous intéressons à la résolution des équations de Maxwell en régime fréquentiel, afin de calculer précisément la signature radar de cibles diverses. Pour avoir une grande précision nécessaire pour des expérience de grande taille, nous utilisons des méthodes d'ordre élevé. Dans le cas scalaire, les éléments finis spectraux hexaédriques avec condensation de masse, permettent d'obtenir un produit matrice vecteur rapide et peux coûteux en stockage. Dans le cas vectoriel, les hexaèdres de la première famille ne réalisent pas la condensation de masse, mais on peut écrire un algorithme rapide de produit matrice-vecteur. Des résultats numériques 3-D montrent la performance de l'algorithme proposé. Nous traitons également le cas o `u la géométrie présente une symétrie de révolution. On est alors ramenés à une succession de problèmes 2-D indépendants. Nous proposons une méthode éléments finis d'ordre élevé couplée à des équations intégrales d'ordre élevé.
Author: Ludovic Gatard
Publisher:
ISBN:
Category :
Languages : fr
Pages : 251
Get Book Here
Book Description
Pour les hautes fréquences, nous résolvons les équations de Maxwell en régime harmonique en utilisant la formulation intégrale de Després. Darrigrand a initié une méthode de résolution, notée FMD, couplant la méthode de discrétisation microlocale de Abboud, Nédélec et Zhou, et la méthode multipôle FMM. Nous avons intégré des éléments finis d'ordre élevé à la FMD afin d'améliorer sa précision et sa convergence en maillage. Nous avons proposé une optimisation de la FMD permettant d'obtenir une méthode dont la complexité est quasi-linéaire. Des tests numériques ont montré l'efficacité de la FMD optimisée d'ordre élevé.
Author:
Publisher:
ISBN:
Category : Aeronautics
Languages : en
Pages : 466
Get Book Here
Book Description
Author: Hélène Fourment-Gaible
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
CE MEMOIRE EST CONSACRE A L'ETUDE NUMERIQUE D'UNE FORMULATION MATHEMATIQUE, PERMETTANT LE CALCUL DU CHAMP ELECTROMAGNETIQUE DIFFRACTE PAR UNE STRUCTURE DIELECTRIQUE COMPLEXE, SOUMISE A UNE ONDE INCIDENTE PLANE. LE PROBLEME CONSISTE A CONSTRUIRE LA SOLUTION DES EQUATIONS DE MAXWELL EN REGIME HARMONIQUE DANS UN DOMAINE BORNE QUI SOIT LA RESTRICTION DE LA SOLUTION DANS TOUT L'ESPACE. LA METHODE UNIMOMENT CLASSIQUE CONSTITUE UNE APPROCHE POSSIBLE, MAIS TIENT EFFECTIVEMENT COMPTE DU CARACTERE BIDIMENSIONNEL D'UN PROBLEME. POUR REMEDIER A CELA, TOUT EN GARDANT L'ASPECT MIXTE DE CETTE METHODE, ON INTRODUIT D'UNE PART DES ELEMENTS FINIS MIXTES, QUI PERMETTENT LA PRISE EN COMPTE NATURELLE DES CONDITIONS D'INTERFACE, ET D'AUTRE PART UNE METHODE DE RACCORDEMENT SUR LA FRONTIERE REGULIERE BORNANT LE DOMAINE. L'APPROXIMATION ELEMENTS FINIS CONSIDEREE CONDUIT A LA RESOLUTION DE SYSTEMES LINEAIRES COMPLEXES SYMETRIQUES, D'ORDRE ELEVE. DEUX ALGORITHMES ITERATIFS, LE BI-GRADIENT CONJUGUE ET LE GRADIENT CONJUGUE CARRE SONT COMPARES. AFIN D'ACCELERER LA CONVERGENCE DE CES METHODES, DIVERS TYPES DE PRECONDITIONNEMENT SONT ENVISAGES. PAR AILLEURS, DANS L'OPTIQUE D'UN STOCKAGE PAR ELEMENT, L'ALGORITHME DE POWELL ET WELSH EST ADAPTE: IL PERMET D'AUGMENTER LE DEGRE DE VECTORISATION. POUR MINIMISER LE COUT INFORMATIQUE, LA VALIDATION DE LA METHODE A ETE REALISEE SUR DES CAS TESTS BIDIMENSIONNELS, DE SOLUTION ANALYTIQUE CONNUE. ENFIN, LES RESULTATS DES EXPERIENCES NUMERIQUES MENEES SUR DEUX CIBLES DEFINIES AU WORKSHOP JINA'90 SONT EN ACCORD AVEC CEUX OBTENUS PAR D'AUTRES APPROCHES.
Author: Adrien Catella
Publisher:
ISBN:
Category :
Languages : fr
Pages : 181
Get Book Here
Book Description
L’objectif général de cette étude est le développement et l’évaluation des schémas en temps efficaces pour des méthodes de type Galerkin discontinu (GD) en maillages tétraédriques non structurés pour la résolution numérique des équations de Maxwell en domaine temporel. Dans la première partie de cette thèse nous rappelons les équations de Maxwell et faisons une rapide revue des principales méthodes numériques utilisées pour résoudre ce système. Dans la seconde partie de cette thèse nous présentons la méthode Galerkin discontinue basée sur des approximations centrées d’ordre générique. Dans ce chapitre nous nous intéresserons qu’aux schémas en temps explicite. Nous détaillerons dans le troisième chapitre la partie principale de ce travail de thèse, c’est-à-dire les schémas implicites en temps, plus particulièrement le schéma implicite très étudié dans la littérature de Crank-Nicolsonn et dans un second temps un schéma implicite d’ordre 4 obtenu à l’aide de la technique du défaut corrigé. Nous réalisons une étude comparative de deux solveurs (direct et intératif) pour la résolution du système linéaire au chapitre 4. Pour des questions d’espace mémoire, nous nous intéressons au chapitre 5 à appliquer le schéma implicite à un sous ensemble du domaine de calcul. Pour cela nous utilisons un schéma hybride explicite/implicite. Au chapitre6, nous présentons les résultats 3D obtenus avec cette méthode. Les problèmes considérés ont plusieurs millions d’inconnues.
Author: Sophie Depeyre
Publisher:
ISBN:
Category :
Languages : fr
Pages : 0
Get Book Here
Book Description
Author: Paul Ayoub
Publisher:
ISBN:
Category :
Languages : fr
Pages : 237
Get Book Here
Book Description
Ce travail a consisté essentiellement en l’élaboration d'un nouveau solveur, pour la résolution des équations de Maxwell tridimensionnelles dans le domaine temporel, répondant aux critères suivants : une méthode de maillage non structure, une méthode d’éléments finis linéaires conformes, un schéma explicite en temps, un contrôle numérique optimal de la contrainte sur la divergence. On sait qu'on peut découpler les équations de maxwell en deux systèmes d’équations, de type équations des ondes. A partir de là, nous avons développé trois formulations différentes, toutes basées sur une régularisation de l’équation d'origine. Après, dans la seconde et la troisième formulation la contrainte de divergence est traitée en un sens faible a l'aide d'un multiplicateur de Lagrange. Dans ces deux cas, le schéma explicite est obtenu, respectivement, par l'utilisation de la méthode de compressibilité artificielle et par pénalisation de la contrainte. La stabilité du problème discret est garantie a l'aide d'une technique de stabilisation. D'après les divers tests numériques de validation effectues, nous avons conclu que la troisième formulation révèle une meilleur précision et robustesse. par conséquent, elle a fait l'objet d'une étude théorique et numérique. la discretisation temporelle est assurée par un schéma aux différences finies. Le code du calcul mi3d a été développé en c++. De nombreux cas tests numériques ont été effectués pour les géométries (conducteur parfait) suivantes : sphère, ogive, cavité cylindrique, avion de chasse et voiture ; et ceci dans le cas d'une source harmonique en temps et d'une impulsion. Finalement, nous avons applique la troisième formulation sur la résolution du problème de diffraction d'ondes électromagnétiques par une structure fine (antenne).
