LOGIN TO YOUR ACCOUNT

Username
Password
Remember Me
Or use your Academic/Social account:

CREATE AN ACCOUNT

Or use your Academic/Social account:

Congratulations!

You have just completed your registration at OpenAire.

Before you can login to the site, you will need to activate your account. An e-mail will be sent to you with the proper instructions.

Important!

Please note that this site is currently undergoing Beta testing.
Any new content you create is not guaranteed to be present to the final version of the site upon release.

Thank you for your patience,
OpenAire Dev Team.

Close This Message

CREATE AN ACCOUNT

Name:
Username:
Password:
Verify Password:
E-mail:
Verify E-mail:
*All Fields Are Required.
Please Verify You Are Human:
fbtwitterlinkedinvimeoflicker grey 14rssslideshare1
Segui, Lionel (2000)
Languages: French
Types: Article
Subjects: Problème inverse, Equation de Helmholtz, Indice de réfraction, Unicité, Milieu PML de Bérenger, Méthode de Gauss-Newton, Régularisation, Semi-norme BV, 510, Inverse problem, Helmholtz equation, Refractive index, Uniqueness, Bérenger’s Perfectly Matched Layer, Gauss-Newton method, Regularization, BV seminorm
La présence de fortes hétérogénéités à l'intérieur d'un obstacle diélectrique ne permet pas la détermination de son indice de réfraction par des méthodes optiques. C'est pourquoi il est nécessaire de se placer dans le domaine de la diffraction. L'indice de réfraction est alors lié à des mesures du champ diffracté par l'équation de Helmholtz. Du point de vue mathématique, ce problème inverse est à la fois non-linéaire et mal posé, que les données soient de type champ proche ou champ lointain. L'objet de cette thèse est de proposer une méthode de résolution approchée pour ce problème inverse, à partir de données du champ proche sur une surface entourant l'obstacle. Contrairement aux problèmes de contrôle dont le but est de minimiser un critère sans pour autant se soucier de l'unicité du minimum, on se préoccupe ici de l'identification de ce minimum, et l'unicité de la solution est de ce fait essentielle. L'indice est ainsi identifiable de manière unique dans L∞(R³ ; C) à partir de la mesure du champ proche pour toute direction incidente de S². Développée dans un premier temps pour le cas de la dimension deux, la méthode doit rester applicable à la dimension trois sans que les temps de calcul ne deviennent pour autant prohibitifs. C'est pourquoi on choisit d'évaluer les sensibilités des mesures à des variations élémentaires de l'indice par une méthode d'approximation de l'EDP plutôt que par la méthode intégrale volumique (équation de Lippmann-Schwinger) habituellement utilisées dans ce type de problème. La méthode choisie inverse des matrices creuses en lieu et place des matrices pleines de la méthode intégrale. L'approximation des conditions aux limites est effectuée de manière très précise grâce à l'introduction, autour du domaine de calcul, d'un milieu fictif absorbant PML de Bérenger. L'algorithme d'inversion choisi est celui de Gauss-Newton, pour lequel on recherche une forme de stabilité en testant successivement deux régularisations de Tychonov (en norme L² et en semi-norme H¹), puis une pénalisation par une régularisée de la semi-norme BV. Cette dernière s'avère, au vu des différents cas tests étudiés, la mieux adaptée pour résoudre ce problème. Due to high heterogeneities that can exist inside a dieletric obstacle, generally optical methods cannot be used in order to retrieve the refractive index of a body. This implies to use scattering methods, where the refractive index is linked to measurements of the scattered field via the Helmholtz equation. From a mathematical point of view, this problem is both nonlinear and ill-posed, whether we consider near field or far field data. The purpose of this thesis is to describe a numerical method to identify the refractive index from near field data on a surface surrounding the obstacle. Whereas control problems only require the minimization of a criteria with no necessary information on the uniqueness of this minimum, here we need to identify this minimum, and therefore the uniqueness of the solution is essential. We study the uniqueness of the solution in L∞(R³ ; C) from near field data for all incident directions in S². Although developped for the two dimensional case, the numerical method was built in order to be fast enough to have applications in the 3D case. We choose to calculate the measurements sensitivities with respect to elementary variations of the index by using a variational formulation of the PDE rather than by the volumic integral method (the Lippmann-Schwinger equation) usually used to solve this kind of problem. The variational method requires the inversion of sparse matrices instead of full matrices in the integral method. The radiation condition is accurately approximated by a curvilinear coordinates PML formulation around the computation domain. The inversion is performed through the Gauss-Newton algorithm, and the stability is assured by means of penalization with two different types of Tychonov regularization (L² norm and H¹ seminorm), and by a BV regularized seminorm. This last penalization appears to be the best fitted to solve our problem.
  • No references.
  • No related research data.
  • No similar publications.

Share - Bookmark

Cite this article