Voir la source de Méthode des points sources distribués

{{Voir homonymes|DPSM}}
{{ébauche|analyse}}

En [[analyse numérique]], la '''méthode des points sources distribués ''', plus connue sous son acronyme anglais '''DPSM''' (''Distributed Point Source Method''), est une méthode permettant de transformer un problème continu (par exemple une [[équation différentielle]]) en un problème discret, de type [[méthode des éléments finis de frontière]]. Cette méthode a été mise au point par les chercheurs [[Dominique Placko]] (1984) et Tribikram Kundu.

Le principe de la méthode consiste à calculer numériquement des valeurs (telles que pression, potentiel, vitesse...) dans un domaine donné à partir des équations différentielles définissant le système, des conditions imposées par l'utilisateur aux frontières de ce domaine et du choix du recul des sources. Cette méthode diffère des méthodes des éléments finis de frontières dans la mesure où le recul des sources supprime les singularités des fonctions de Green sur les frontières, inhérentes à ce genre de méthode.

Cette méthode est couramment utilisée dans les disciplines telles que l'électrostatique, l'électromagnétisme, les ultrasons et le [[contrôle non destructif]]. 

== La méthode ==

=== Discrétisation ===

On part de la formulation faible du problème et on utilise une [[méthode de Galerkine]] sur un maillage du domaine d'étude. La solution appartient à un espace fonctionnel satisfaisant des propriétés de régularité bien définies.

D'un point de vue plus formel, on écrit la formulation faible sous la forme : 
:Trouver <math>u \in V</math> telle que <math>\forall v \in V, a(u,v) = L(v)</math>

où {{math|''a''}} est une forme bilinéaire, et {{math|''L''}} une forme linéaire.

L'ensemble {{math|''V''}} étant généralement de dimension infinie, on construit un espace <math>V_h \subset V </math> avec <math>{\rm dim} V_h < + \infty </math>, et on réécrit le problème de la façon suivante : 

:Trouver <math>u_h \in V_h </math> telle que <math>\forall v_h \in V_h, a(u_h,v_h) = L(v_h)</math>

Dans la pratique, on se restreint aux fonctions <math>V</math> définissant les frontières du domaine d'étude. Ces fonctions sont appelées [[Condition_aux_limites|conditions aux limites]] du problème.

L'espace <math>V_h</math> considéré est l'ensemble des [[Fonction_de_Green|fonctions de Green]] solutions ponctuelles pour chaque point source issue du maillage. 


=== Forme matricielle du problème ===

Du fait que l'espace d'approximation utilisé <math>V_h</math> est de dimension finie <math>n < +\infty</math>, on peut décomposer la solution du problème sur une [[Base (algèbre linéaire)|base]] de fonctions de Green <math>\left\lbrace g_i \right\rbrace_{1 \leqslant i \leqslant n}</math> de <math>V_h</math> :
:<math>u_h = \sum_{j=1}^n {\lambda}_j g_j</math>

Ainsi, en écrivant le problème en choisissant les fonctions condition aux limites <math>v_h = cl_i</math>, il vient :
:<math>m(u_h,e_i) = \sum_{j=1}^n {\lambda}_j a(g_j,cl_i) = L(cl_i) \ , \ \forall i \in [ \! [ 1 , n ] \! ]</math>

On obtient ainsi un [[système d'équations linéaires]] de la forme <math>M{\lambda} = l</math>, en notant
:<math>M = \left( a(g_j,cl_i) \right)_{1 \leqslant i,j \leqslant n}</math>, <math>{\lambda} = \left( {\lambda}_i \right)_{1 \leqslant i \leqslant n}</math>, <math>l = \left( L(cl_i) \right)_{1 \leqslant i \leqslant n}</math>

Cette forme matricielle est plus lisible sous la forme de [[Matrice_par_blocs|matrices par blocs]]. Les sources et les points tests sont regroupés par frontières du domaine. 

<math>m(u_h,e_i) = \sum_{j=1}^n {\lambda}_j a(g_j,cl_i) = \sum_{Fr=1}^p {\sum_{j=1}^{n_p} {\lambda}_j a(g_j,cl_i)} </math>


== Illustration de la méthode ==

=== Le modèle étudié ===
[[File:Points DPSM.png|thumb|Modélisation d'un transducteur à ultrasons dans un milieu contenant 2 bulles]]
Exemple de calcul de la [[pression acoustique]] générée par un [[transducteur]] à [[Ultrason|ultrasons]] dans un milieu contenant 2 bulles. 


=== Les éléments de la méthode ===

Les "points test" sont des points géométriques posés sur le maillage des frontières du domaine d'étude. 

Les "conditions aux limites" sont des fonctions définies sur les points tests. Elles sont scalaires ou vectorielles .

Les "points sources"  sont des points géométriques positionnés à l'aplomb des points tests.

Le "milieu" est l'espace délimité par les frontières. Le milieu est supposé homogène et il est défini par ses caractéristiques (comme la pression, la permittivité...)

Les fonctions sources de cette modélisation définies entre une source S et un point test P sont 

<math>a_{SP} = \frac {\mu}{4 \pi}\frac {\mathrm e^{-ikR_{SP}}}{R_{SP}}   </math> où <math> R_{SP}</math> est la distance entre le point source S et le point test P.

Ces fonctions sont les coefficients des matrices. Ils correspondent à la contribution des sources sur les points tests.

=== Génération de la matrice ===
[[File:DPSM matrice.png|center|936x176px|Représentation des différents blocs de la matrice de couplage pour un transducteur acoustique et 2 bulles]]

La matrice est générée bloc par bloc. Chaque bloc correspond à un couple de frontières <nowiki><Sources/Points test></nowiki>. Chaque bloc est une matrice de dimension de NxP coefficients, N étant le nombre de sources et P le nombre de points tests. 

Les blocs diagonaux sont appelés matrices d'autocouplage, ils correspondent à la contribution des sources d'une frontière sur les points tests de la même frontière. Les autres blocs sont les matrices d'intercouplage, correspondant à la contribution des sources d'une frontière sur les points tests d'une autre frontière. 

La ligne de bloc s'écrit : <math>M.\lambda = Cl </math>


Dans les cas où les frontières sont internes, c'est-à-dire qu'elles séparent 2 milieux, les conditions aux limites utilisateurs sont remplacées par des conditions de continuité comme la continuité de la pression modulo la pression de chaque milieu.

La ligne de bloc s'écrit : <math>M^+.\lambda^+ -M^-.\lambda^- = 0 </math>  où  <math>M^+ </math> correspond à la contribution des sources d'un côté de la frontière interne orientée et <math>M^- </math> au côté opposé.

=== les résultats ===
[[File:DPSM resultat US.png|320x211px|center|Pression acoustique générée par un transducteur vibrant à {{unité|1|MHz}} et {{unité|2|MHz}} dans de l'eau contenant 2 bulles d'air]]

Pression acoustique générée par un transducteur vibrant à {{unité|1|MHz}} et {{unité|2|MHz}} dans de l'eau contenant 2 bulles d'air

== Voir aussi ==
*[[Méthode des éléments finis]]
*[[Méthode de Galerkine]]
*[[méthode des éléments finis de frontière]]

== Référence ==
* {{en}} Advanced Ultrasonic Methods for Material and Structure Inspection, Auteurs:R.AHMAD, T. KUNDU éditeur Hermes Science Publishing, Royaume-Uni, {{date-|janvier 2007}}, {{Chap.}} 5 {{p.|159-187}}. 

* {{en}} [http://www.jpier.org/PIER/pier148/18.14032405.pdf A Differential DPSM Based Modeling Applied to Eddy Current Imaging Problems (Thierry Bore)]

* {{Ouvrage
| langue=en
| prénom1=Dominique
| nom1=Placko
| lien auteur1=Dominique_Placko
| lien auteur2=Tribikram Kundu
| titre=Dpsm for Modeling Engineering Problems
| éditeur=
| année=2007
| mois=juillet
| isbn=978-0-471-73314-0
}}
* {{Ouvrage
| langue=en
| prénom1=Tribikram
| nom1=Kundu
| titre=Ultrasonic and Electromagnetic Nde for Structure and Material Characterization
| sous-titre=Engineering and Biomedical Applications
| éditeur=[[CRC Press]]
| année=2012
| mois=Aout
| pages totales=890
| isbn=978-1-4398-3663-7
| lire en ligne=https://books.google.com/books?id=yhP2FJgn25wC&printsec=frontcover
}}

{{Portail|Analyse}}

{{DEFAULTSORT:Methode des Points sources distribués}}
[[Catégorie:Méthode mathématique|Galerkine]]
[[Catégorie:Méthode des éléments finis]]