Spelling suggestions: "subject:"méthode dde newton"" "subject:"méthode dde mewton""
11 |
Simulation numérique de la solidification avec réduction de modèle PGD appliquée à la fonderie / Numerical simulation of solidification with reduced model order PGD applied foundryDespret, Pierre 08 October 2015 (has links)
La thèse CIFRE s'est déroulée dans un contexte de métallurgie industrielle et de simulation numérique. La modélisation de la solidification, via l'équation de la chaleur et avec des méthodes de réduction de modèle, était un objectif majeur. L'entreprise Montupet, spécialisée dans la fonderie d'aluminium, est le porteur du projet et financeur de la thèse. L'université de Technologie de Compiègne (UTC) a réalisé l'accompagnement académique. La méthode PGD “Proper General Decomposition”, basée sur une séparation de variables, est l'objet de nombreuses recherches. Nous avons proposé, concernant des propriétés matériaux non-linéaires, une discrétisation spatio-temporelle des matrices matériaux. Avec une formulation en température, sans chaleur latente, les gains sont élevés. L'introduction de la chaleur latente réduit fortement les gains. Nous formulons l'hypothèse que la difficulté de convergence de la méthode PGD dans le cas de la solidification repose sur une formulation en température inadaptée. Nous décidons d'opter une formulation en enthalpie. Il s'avère que cette formulation offre des perspectives encourageantes, mais nécessite encore beaucoup de développements. En parallèle de ces développements, un séjour de 5 mois aux États-Unis a été réalisé afin d'obtenir une meilleure caractérisation de la fraction solide. La recherche s'est portée sur l'évolution de la fraction solide en fonction de la vitesse de refroidissement. Sous réserve de mesures complémentaires, les essais ont mis en évidence une modification de la courbe de fraction solide en fonction de la vitesse de refroidissement, notamment un agrandissement de l'intervalle de solidification. / The PhD Thesis was carried out in a metallurgy and numerical simulation environment. The main topic was to model solidification, thought heat equation formulation and reduced order model PGD resolution. Montupet, specialized in aluminium alloys foundry hold and financed the project, the Université de Technologie de Compiègne did the acadernic supervising. The PGD method "Proper General Decomposition" is a hot topic based on variable separation. We proposed, regarding the non-linear materials, a space-time discretization of material matrix. With a temperature formulation, without latent heat, gains are high. With latent heat, gains fall drastically. We proposed the hypothesis that temperature could be an inadapted formulation. We decided to use the enthalpy formulation. This formulation offers good perspectives but needs more developments. During the thesis, five months were spent in the USA to get a better caracterisation of the solid fraction, particularly its variation in function of the cooling rate. Under reservation, the samples show a modification of solid fraction curves and particularly a change of solidification interval in function of cooling rate.
|
12 |
Modélisation, observation et commande d’une classe d’équations aux dérivées partielles : application aux matériaux semi-transparents / Modeling, analysis and control for a class of partial differential equations : application to thermoforming of semi-transparent materialsGhattassi, Mohamed 29 September 2015 (has links)
Le travail présenté dans ce mémoire nous a permis d’étudier d’un point de vue théorique et numérique le transfert de chaleur couplé par rayonnement et conduction à travers un milieu semi-transparent, gris et non diffusant dans une géométrie multidimensionnelle 2D. Ces deux modes de transfert de chaleur sont décrits par un couplage non linéaire de l’équation de la chaleur non linéaire (CT) et de l’équation du transfert radiatif (ETR). Nous avons présenté des résultats d’existence, d’unicité locale de la solution pour le système couplé avec des conditions aux limites de type Dirichlet homogènes en utilisant le théorème du point fixe de Banach. Par ailleurs, les travaux réalisés nous ont permis de mettre au point un code de calcul qui permet de simuler la température. Nous avons utilisé la quadrature S_N pour la discrétisation angulaire de l’ETR. La discrétisationde l'ETR dans la variable spatiale est effectuée par la méthode de Galerkin discontinue (DG) et en éléments finis pour l'équation de la chaleur non linéaire. Nous avons démontré la convergence du schémanumérique couplé en utilisant la méthode du point fixe discret. Le modèle discret, sous la forme d’équations différentielles ordinairesnon linéaires obtenu après une approximation nous a permis de fairel’analyse et la synthèse d’estimateurs d’état et de lois de commandepour la stabilisation. Grâce à la structure particulière du modèle età l’aide du DMVT. Nous avons proposé un observateur d’ordre réduit.D’autre part nous avons réussi à construire une matrice de gain quiassure la stabilité de l’observateur proposé. Une extension au filtrage $\mathcal{H}_{\infty}$ est également proposée. Une nouvelleinégalité matricielle (LMI) est donnée dans le cas d’une commandebasée observateur. Nous avons étendu à l’approche d’ordre réduit dans le cas de la commande basée observateur et nous avons montré la stabilité sous l’action de la rétroaction. De même une extension au filtrage $\mathcal{H}_{\infty}$ est également proposée. Tous les résultats sont validés par des simulations numériques. / This thesis investigates the theoretical and numerical analysis of coupled radiative conductive heat transfer in a semi-transparent, gray and non-scattering 2D medium. This two heat transfer modes are described by the radiative transfer equation (RTE) and the nonlinear heat equation (NHE). We proved the existence and uniqueness of the solution of coupled systems with homogeneous Dirichlet boundary conditions using the fixed-point theorem. Moreover, we developed a useful algorithm to simulate the temperature in the medium. We used the quadrature $S_{N}$ for the angular discretization of the RTE. The spatial discretization of RTE was made by the discontinuous Galerkin method (DG) and the finite element method for the non-linear heat equation. We have shown the convergence and the stability of the coupled numerical scheme using the discrete fixed point. The discrète model obtained after an approximation allowed us to do the analysis and synthesis of state estimators and feedback control design for stabilization of the system. Thanks to the special structure of the model and using the Differential Mean Value Theorem (DMVT), we proposed a reduced order observer and we construct a gain matrix, which ensures the exponential stability of the proposed observer and guarantees the boundedness of the estimate vector. An extension to $\mathcal{H}_{\infty}$ filtering is also provided. We have extended the reduced order approach in the case of the observer-based controller and we proved the exponential stability under the control feedback law. Similarly, an extension to $\mathcal{H}_{\infty}$ filtering is also provided. The obtained results were validated through several numerical simulations.
|
13 |
Méthodes d'intégration produit pour les équations de Fredholm de deuxième espèce : cas linéaire et non linéaire / Product integration methods for Fredholm integral equations of the second kind : linear case and nonlinear caseKaboul, Hanane 20 June 2016 (has links)
La méthode d'intégration produit a été proposée pour résoudre des équations linéaires de Fredholm de deuxième espèce singulières dont la solution exacte est régulière, au moins continue. Dans ce travail on adapte cette méthode à des équations dont la solution est juste intégrable. On étudie également son extension au cas non linéaire posé dans l'espace des fonctions intégrables. Ensuite, on propose une autre manière de mettre en oeuvre la méthode d'intégration produit : on commence par linéariser l'équation par une méthode de type Newton puis on discrétise les itérations de Newton par la méthode d'intégration produit / The product integration method has been proposed for solving singular linear Fredholm equations of the second kind whose exact solution is smooth, at least continuous. In this work, we adapt this method to the case where the solution is only integrable. We also study the nonlinear case in the space of integrable functions. Then, we propose a new version of the method in the nonlinear framework : we first linearize the eqaution by a Newton type method and then discretize the Newton iterations by the product integration method
|
14 |
Contribution à la Résolution Numérique de Problèmes Inverses de Diffraction Élasto-acoustique / Contribution to the Numerical Reconstruction in Inverse Elasto-Acoustic ScatteringAzpiroz, Izar 28 February 2018 (has links)
La caractérisation d’objets enfouis à partir de mesures d’ondes diffractées est un problème présent dans de nombreuses applications comme l’exploration géophysique, le contrôle non-destructif, l’imagerie médicale, etc. Elle peut être obtenue numériquement par la résolution d’un problème inverse. Néanmoins, c’est un problème non linéaire et mal posé, ce qui rend la tâche difficile. Une reconstruction précise nécessite un choix judicieux de plusieurs paramètres très différents, dépendant des données de la méthode numérique d’optimisation choisie.La contribution principale de cette thèse est une étude de la reconstruction complète d’obstacles élastiques immergés à partir de mesures du champ lointain diffracté. Les paramètres à reconstruire sont la frontière, les coefficients de Lamé, la densité et la position de l’obstacle. On établit tout d’abord des résultats d’existence et d’unicité pour un problème aux limites généralisé englobant le problème direct d’élasto-acoustique. On analyse la sensibilité du champ diffracté par rapport aux différents paramètres du solide, ce qui nous conduit à caractériser les dérivées partielles de Fréchet comme des solutions du problème direct avec des seconds membres modifiés. Les dérivées sont calculées numériquement grâce à la méthode de Galerkine discontinue avec pénalité intérieure et le code est validé par des comparaisons avec des solutions analytiques. Ensuite, deux méthodologies sont introduites pour résoudre le problème inverse. Toutes deux reposent sur une méthode itérative de type Newton généralisée et la première consiste à retrouver les paramètres de nature différente indépendamment, alors que la seconde reconstruit tous les paramètre en même temps. À cause du comportement différent des paramètres, on réalise des tests de sensibilité pour évaluer l’influence de ces paramètres sur les mesures. On conclut que les paramètres matériels ont une influence plus faible sur les mesures que les paramètres de forme et, ainsi, qu’une stratégie efficace pour retrouver des paramètres de nature distincte doit prendre en compte ces différents niveaux de sensibilité. On a effectué de nombreuses expériences à différents niveaux de bruit, avec des données partielles ou complètes pour retrouver certains paramètres, par exemple les coefficients de Lamé et les paramètres de forme, la densité, les paramètres de forme et la localisation. Cet ensemble de tests contribue à la mise en place d’une stratégie pour la reconstruction complète des conditions plus proches de la réalité. Dans la dernière partie de la thèse, on étend ces résultats à des matériaux plus complexes, en particulier élastiques anisotropes. / The characterization of hidden objects from scattered wave measurements arises in many applications such as geophysical exploration, non destructive testing, medical imaging, etc. It can be achieved numerically by solving an Inverse Problem. However, this is a nonlinear and ill-posed problem, thus a difficult task. A successful reconstruction requires careful selection of very different parameters depending on the data and the chosen optimization numerical method.The main contribution of this thesis is an investigation of the full reconstruction of immersed elastic scatterers from far-field pattern measurements. The sought-after parameters are the boundary, the Lamé coefficients, the density and the location of the obstacle. First, existence and uniqueness results of a generalized Boundary Value Problem including the direct elasto-acoustic problem are established. The sensitivity of the scattered field with respect to the different parametersdescribing the solid is analyzed and we end up with the characterization of the corresponding partial Fréchet derivatives as solutions to the direct problem with modified right-hand sides. These Fréchet derivatives are computed numerically thanks to the Interior Penalty Discontinuous Galerkin method and the code is validated thanks to comparison with analytical solutions. Then, two solution methodologies are introduced for solving the inverse problem. Both are based on an iterative regularized Newton-type methodology and the first one consists in retrieving the parameters of different nature independently, while the second one reconstructs all parameters together. Due to the different behavior of the parameters, sensitivity tests are performed to assess the impact of the parameters on the measurements. We conclude that material parameters have a weaker influence on the measurements than shape parameters, and therefore, a successful strategy to retrieve parameters of distinct nature should take into account these different levels of sensitivity. Various experiments at different noise levels and with full or limited aperture data are carried out to retrieve some of the physical properties, e.g. Lamé coefficients with shape parameters, density with shape parameters a, density, shape and location. This set of tests contributes to a final strategy for the full reconstruction and in more realistic conditions. In the final part of the thesis, we extend the results to more complex material parameters, in particular anisotropic elastic.
|
15 |
Sur des systèmes dynamiques dissipatifs de type gradient. Applications en Optimisation.Bolte, Jérôme 06 January 2003 (has links) (PDF)
L'étude et l'introduction de nouveaux systèmes dynamiques<br /> de type gradient sont l'objet central de cette thèse. Le<br /> caractère dissipatif de telles dynamiques est au coeur de<br /> nombreux domaines en mathématiques : optimisation,<br /> mécanique, équations d'évolutions en dimension infinie.<br /><br />Dans une première partie, les champs de gradients (ou de sous-différentiels<br /> de fonction convexe) sont contrôlés à l'aide d'opérateurs-barrières. <br />La motivation essentielle est d'obtenir<br /> des méthodes intérieures de descente en vue d'optimiser<br /> une fonction sous des contraintes convexes. Le cadre<br /> d'étude proposé permet d'unifier dans un même formalisme de nombreuses<br /> méthodes continues : gradient projeté, plus grande pente riemannienne,<br /> méthode continue de Newton... Parmi les conséquences de <br />la généralisation proposée, on peut, par exemple, évoquer des <br /> résultats abstraits de viabilité et de convergence globale. Toujours <br />dans cette <br />perspective, les fonctions de Legendre jouent un rôle crucial~:<br /> elles permettent d'une part de donner lieu à des structures<br /> riemanniennes possédant de nombreuses propriétés - parmi lesquelles une<br /> propriété d'intégration caractéristique remarquable -, et d'autre part, <br /> elles fournissent en dimension infinie un cadre intéressant<br /> pour l'étude de certaines équations d'évolution de type<br /> parabolique.<br /><br />La deuxième partie est consacrée à l'étude de systèmes<br /> dynamiques du second ordre en temps avec une dissipation géométrique<br /> de type hessien. Outre leur intérêt en optimisation<br /> et leurs liens avec les méthodes de type Newton, ces systèmes<br /> sont d'une grande souplesse et permettent d'approcher certains <br />phénomènes non-lisses en mécanique unilatérale. En guise d'application,<br /> il est en effet prouvé que les systèmes considérés permettent <br />d'obtenir à la limite des dynamiques <br />satisfaisant des lois de chocs inélastiques. Les<br /> perspectives de cette étude ouvrent en particulier la voie à une approche <br />alternative de certains systèmes d'inégalités variationnelles de type <br />hyperbolique.<br /><br /><br />L'une des préoccupations majeures de cette thèse est la question<br /> de la convergence des orbites des systèmes étudiés. Dans le <br /> cadre de la minimisation convexe, quasi-convexe, ou analytique, de nombreux<br /> résultats sont proposés : convergence globale, , <br />vitesse de convergence, contrôle asymptotique, attractivité des <br /> minima sous contraintes en dimension infinie.
|
16 |
Méthode de Newton revisitée pour les équations généralisées / Newton-type methods for solving inclusionsNguyen, Van Vu 30 September 2016 (has links)
Le but de cette thèse est d'étudier la méthode de Newton pour résoudre numériquement les inclusions variationnelles, appelées aussi dans la littérature les équations généralisées. Ces problèmes engendrent en général des opérateurs multivoques. La première partie est dédiée à l'extension des approches de Kantorovich et la théorie (alpha, gamma) de Smale (connues pour les équations non-linéaires classiques) au cas des inclusions variationnelles dans les espaces de Banach. Ceci a été rendu possible grâce aux développements récents des outils de l'analyse variationnelle et non-lisse tels que la régularité métrique. La seconde partie est consacrée à l'étude de méthodes numériques de type-Newton pour les inclusions variationnelles en utilisant la différentiabilité généralisée d'applications multivoques où nous proposons de linéariser à la fois les parties univoques (lisses) et multivoques (non-lisses). Nous avons montré que, sous des hypothèses sur les données du problème ainsi que le choix du point de départ, la suite générée par la méthode de Newton converge au moins linéairement vers une solution du problème de départ. La convergence superlinéaire peut-être obtenue en imposant plus de conditions sur l'approximation multivaluée. La dernière partie de cette thèse est consacrée à l'étude des équations généralisées dans les variétés Riemaniennes à valeurs dans des espaces euclidiens. Grâce à la relation entre la structure géométrique des variétés et les applications de rétractions, nous montrons que le schéma de Newton converge localement superlinéairement vers une solution du problème. La convergence quadratique (locale et semi-locale) peut-être obtenue avec des hypothèses de régularités sur les données du problème. / This thesis is devoted to present some results in the scope of Newton-type methods applied for inclusion involving set-valued mappings. In the first part, we follow the Kantorovich's and/or Smale's approaches to study the convergence of Josephy-Newton method for generalized equation (GE) in Banach spaces. Such results can be viewed as an extension of the classical Kantorovich's theorem as well as Smale's (alpha, gamma)-theory which were stated for nonlinear equations. The second part develops an algorithm using set-valued differentiation in order to solve GE. We proved that, under some suitable conditions imposed on the input data and the choice of the starting point, the algorithm produces a sequence converging at least linearly to a solution of considering GE. Moreover, by imposing some stronger assumptions related to the approximation of set-valued part, the proposed method converges locally superlinearly. The last part deals with inclusions involving maps defined on Riemannian manifolds whose values belong to an Euclidean space. Using the relationship between the geometric structure of manifolds and the retraction maps, we show that, our scheme converges locally superlinearly to a solution of the initial problem. With some more regularity assumptions on the data involved in the problem, the quadratic convergence (local and semi-local) can be ensured.
|
17 |
Les méthodes numériques de transport réactif / Numerical methods for reactive transportSabit, Souhila 27 May 2014 (has links)
La modélisation du transport réactif du contaminant en milieu poreux est un problème complexe cumulant les difficultés de la modélisation du transport avec celles de la modélisation de la chimie et surtout du couplage entre les deux. Cette modélisation conduit à un système d'équations aux dérivées partielles et algébriques dont les inconnues sont les quantités d'espèces chimiques. Une approche possible, déjà utilisée par ailleurs, est de choisir la méthode globale DAE : l'utilisation d'une méthode de lignes, correspondant à la discrétisation en espace seulement, conduit à un système différentiel algébrique (DAE) qui doit être résolu par un solveur adapté. Dans notre cas, on utilise le solveur IDA de Sundials qui s'appuie sur une méthode implicite, à ordre et pas variables, et qui requiert à chaque pas de temps la résolution d'un grand système non linéaire associé à une matrice jacobienne. Cette méthode est implémentée dans un logiciel qui s'appelle GRT3D (Transport Réactif Global en 3D). Le présent travail présente une amélioration de la méthode GDAE, du point de vue de la performance, de la stabilité et de la robustesse. Nous avons ainsi enrichi les possibilités de GRT3D, par la prise en compte complète des équations de précipitation-dissolution permettant l'apparition ou la disparition d'une espèce précipitée. En complément de l'étude de la méthode GDAE, nous présentons aussi une méthode séquentielle non itérative (SNIA), qui est une méthode basée sur le schéma d'Euler explicite : à chaque pas de temps, on résout explicitement l'équation de transport et on utilise ces calculs comme données pour le système chimique, résolu dans chaque maille de façon indépendante. Nous présentons aussi une comparaison entre cette méthode et l'approche GDAE. Des résultats numériques pour deux cas tests, celui proposé par l'ANDRA (cas-test 2D) d'une part, celui proposé par le groupe MoMas (Benchmark "easy case") d'autre part, sont enfin présentés, commentés et analysés. / Modeling reactive transport of contaminants in porous media is a complex problem combining the difficulties of modeling the trasport with those of modeling the chemistry and especially the coupling between the two .This model leads to a system of partial differential equations and algebraic equations whose unknowns are the quantities of chemical species. One approach , already used elsewhere , is choosing the global DAE method : using the method of lines, discretization in space only, leads to a differential algebraic system (DAE ) to be solved by a suitable solver . In our case , the solver IDA Sundials relies on an implicit method, order is used but not variables, and requires at each time solving a large nonlinear system associated with a Jacobian matrix . This method is implemented in a software called GRT3D (Global Reactive Transport in 3D). This paper presents an improved GDAE method , from the standpoint of performance, the stability and robustness. We have enriched the possibilities of GRT3D , by taking full account of the equations of dissolution – precipitation for the appearance or disappearance of precipitated species. In addition to the study of the GDAE method, we also present a non-iterative sequential method ( SNIA ) which is a method based on the explicit Euler scheme : at each time step, we explicitly solve the transport equation and we use these calculations as data for the chemical system which is resolved in each cell independently. We also present a comparison between this method and GDAE approach . Numerical results for two test cases , one proposed by ANDRA ( 2D test case ) on one hand and one proposed by the group MOMAS ( Benchmark "easy case" ) on the other hand, are finally presented , discussed and analyzed.
|
18 |
GPU-enhanced power flow analysis / Calcul de Flux de Puissance amélioré grâce aux Processeurs GraphiquesMarin, Manuel 11 December 2015 (has links)
Cette thèse propose un large éventail d'approches afin d'améliorer différents aspects de l'analyse des flux de puissance avec comme fils conducteur l'utilisation du processeurs graphiques (GPU). Si les GPU ont rapidement prouvés leurs efficacités sur des applications régulières pour lesquelles le parallélisme de données était facilement exploitable, il en est tout autrement pour les applications dites irrégulières. Ceci est précisément le cas de la plupart des algorithmes d'analyse de flux de puissance. Pour ce travail, nous nous inscrivons dans cette problématique d'optimisation de l'analyse de flux de puissance à l'aide de coprocesseur de type GPU. L'intérêt est double. Il étend le domaine d'application des GPU à une nouvelle classe de problème et/ou d'algorithme en proposant des solutions originales. Il permet aussi à l'analyse des flux de puissance de rester pertinent dans un contexte de changements continus dans les systèmes énergétiques, et ainsi d'en faciliter leur évolution. Nos principales contributions liées à la programmation sur GPU sont: (i) l'analyse des différentes méthodes de parcours d'arbre pour apporter une réponse au problème de la régularité par rapport à l'équilibrage de charge ; (ii) l'analyse de l'impact du format de représentation sur la performance des implémentations d'arithmétique floue. Nos contributions à l'analyse des flux de puissance sont les suivantes: (ii) une nouvelle méthode pour l'évaluation de l'incertitude dans l'analyse des flux de puissance ; (ii) une nouvelle méthode de point fixe pour l'analyse des flux de puissance, problème que l'on qualifie d'intrinsèquement parallèle. / This thesis addresses the utilization of Graphics Processing Units (GPUs) for improving the Power Flow (PF) analysis of modern power systems. Currently, GPUs are challenged by applications exhibiting an irregular computational pattern, as is the case of most known methods for PF analysis. At the same time, the PF analysis needs to be improved in order to cope with new requirements of efficiency and accuracy coming from the Smart Grid concept. The relevance of GPU-enhanced PF analysis is twofold. On one hand, it expands the application domain of GPU to a new class of problems. On the other hand, it consistently increases the computational capacity available for power system operation and design. The present work attempts to achieve that in two complementary ways: (i) by developing novel GPU programming strategies for available PF algorithms, and (ii) by proposing novel PF analysis methods that can exploit the numerous features present in GPU architectures. Specific contributions on GPU computing include: (i) a comparison of two programming paradigms, namely regularity and load-balancing, for implementing the so-called treefix operations; (ii) a study of the impact of the representation format over performance and accuracy, for fuzzy interval algebraic operations; and (iii) the utilization of architecture-specific design, as a novel strategy to improve performance scalability of applications. Contributions on PF analysis include: (i) the design and evaluation of a novel method for the uncertainty assessment, based on the fuzzy interval approach; and (ii) the development of an intrinsically parallel method for PF analysis, which is not affected by the Amdahl's law.
|
19 |
Méthodes itératives pour la résolution d'équations matricielles / Iterative methods fol solving matrix equationsSadek, El Mostafa 23 May 2015 (has links)
Nous nous intéressons dans cette thèse, à l’étude des méthodes itératives pour la résolutiond’équations matricielles de grande taille : Lyapunov, Sylvester, Riccati et Riccatinon symétrique.L’objectif est de chercher des méthodes itératives plus efficaces et plus rapides pour résoudreles équations matricielles de grande taille. Nous proposons des méthodes itérativesde type projection sur des sous espaces de Krylov par blocs Km(A, V ) = Image{V,AV, . . . ,Am−1V }, ou des sous espaces de Krylov étendus par blocs Kem(A, V ) = Image{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V } . Ces méthodes sont généralement plus efficaces et rapides pour les problèmes de grande dimension. Nous avons traité d'abord la résolution numérique des équations matricielles linéaires : Lyapunov, Sylvester, Stein. Nous avons proposé une nouvelle méthode itérative basée sur la minimisation de résidu MR et la projection sur des sous espaces de Krylov étendus par blocs Kem(A, V ). L'algorithme d'Arnoldi étendu par blocs permet de donner un problème de minimisation projeté de petite taille. Le problème de minimisation de taille réduit est résolu par différentes méthodes directes ou itératives. Nous avons présenté ainsi la méthode de minimisation de résidu basée sur l'approche global à la place de l'approche bloc. Nous projetons sur des sous espaces de Krylov étendus Global Kem(A, V ) = sev{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V }. Nous nous sommes intéressés en deuxième lieu à des équations matricielles non linéaires, et tout particulièrement l'équation matricielle de Riccati dans le cas continu et dans le cas non symétrique appliquée dans les problèmes de transport. Nous avons utilisé la méthode de Newtown et l'algorithme MINRES pour résoudre le problème de minimisation projeté. Enfin, nous avons proposé deux nouvelles méthodes itératives pour résoudre les équations de Riccati non symétriques de grande taille : la première basée sur l'algorithme d'Arnoldi étendu par bloc et la condition d'orthogonalité de Galerkin, la deuxième est de type Newton-Krylov, basée sur la méthode de Newton et la résolution d'une équation de Sylvester de grande taille par une méthode de type Krylov par blocs. Pour toutes ces méthodes, les approximations sont données sous la forme factorisée, ce qui nous permet d'économiser la place mémoire en programmation. Nous avons donné des exemples numériques qui montrent bien l'efficacité des méthodes proposées dans le cas de grandes tailles. / In this thesis, we focus in the studying of some iterative methods for solving large matrix equations such as Lyapunov, Sylvester, Riccati and nonsymmetric algebraic Riccati equation. We look for the most efficient and faster iterative methods for solving large matrix equations. We propose iterative methods such as projection on block Krylov subspaces Km(A, V ) = Range{V,AV, . . . ,Am−1V }, or block extended Krylov subspaces Kem(A, V ) = Range{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V }. These methods are generally most efficient and faster for large problems. We first treat the numerical solution of the following linear matrix equations : Lyapunov, Sylvester and Stein matrix equations. We have proposed a new iterative method based on Minimal Residual MR and projection on block extended Krylov subspaces Kem(A, V ). The extended block Arnoldi algorithm gives a projected minimization problem of small size. The reduced size of the minimization problem is solved by direct or iterative methods. We also introduced the Minimal Residual method based on the global approach instead of the block approach. We projected on the global extended Krylov subspace Kem(A, V ) = Span{V,A−1V,AV,A−2V,A2V, · · · ,Am−1V,A−m+1V }. Secondly, we focus on nonlinear matrix equations, especially the matrix Riccati equation in the continuous case and the nonsymmetric case applied in transportation problems. We used the Newton method and MINRES algorithm to solve the projected minimization problem. Finally, we proposed two new iterative methods for solving large nonsymmetric Riccati equation : the first based on the algorithm of extended block Arnoldi and Galerkin condition, the second type is Newton-Krylov, based on Newton’s method and the resolution of the large matrix Sylvester equation by using block Krylov method. For all these methods, approximations are given in low rank form, wich allow us to save memory space. We have given numerical examples that show the effectiveness of the methods proposed in the case of large sizes.
|
20 |
Contribution à l'analyse mathématique et à la résolution numérique d'un problème inverse de scattering élasto-acoustiqueEstecahandy, Elodie 19 September 2013 (has links) (PDF)
La détermination de la forme d'un obstacle élastique immergé dans un milieu fluide à partir de mesures du champ d'onde diffracté est un problème d'un vif intérêt dans de nombreux domaines tels que le sonar, l'exploration géophysique et l'imagerie médicale. A cause de son caractère non-linéaire et mal posé, ce problème inverse de l'obstacle (IOP) est très difficile à résoudre, particulièrement d'un point de vue numérique. De plus, son étude requiert la compréhension de la théorie du problème de diffraction direct (DP) associé, et la maîtrise des méthodes de résolution correspondantes. Le travail accompli ici se rapporte à l'analyse mathématique et numérique du DP élasto-acoustique et de l'IOP. En particulier, nous avons développé un code de simulation numérique performant pour la propagation des ondes associée à ce type de milieux, basé sur une méthode de type DG qui emploie des éléments finis d'ordre supérieur et des éléments courbes à l'interface afin de mieux représenter l'interaction fluide-structure, et nous l'appliquons à la reconstruction d'objets par la mise en oeuvre d'une méthode de Newton régularisée.
|
Page generated in 0.0613 seconds