Bonjour,
Voilà j'aurais plusieurs questions sur l'algorithme de Yee (FDTD appliqué à la résolution des équations de maxwell).
J'ai un document vraiment bien rédigé -M. D. Sullivan- "Electromagnetic simulation using the FDTD method". Ce livre présente entre autre des exemples écrit en C permettant de simuler toute sortes de cas (vraiment génial est très pédagogique comme livre).
J'ai récris ces exemples dans le langage octave. Et je dispose maintenant d'un code permettant de calculer le mode TE en 2D avec PML. Mais je m'interesse plutôt au cas TM. Et ce cas n'est pas vraiment discuter dans le livre.
L’algorithme de base pour le cas TE (non magnétique) consiste en :
calculer D, avec : dD/dt=rotH
Inserer les terme de source D=D+source
Calculer E, avec E=D/epsilon (seule équation qui dépends du probleme)
Calculer H avec dH/dt=-1/mu0*rotE
et reprendre... jusqu’à l'instant finale (chaque boucle correspondant à 1 step temporel)
Question :
Cette architecture de base s'applique telle dans le cas TM ?
Notamment le terme de source... Comment faire ?
Je pensais faire :
calculer D, avec : dD/dt=rotH
Calculer E, avec E=D/epsilon (seule équation qui dépends du probleme)
Calculer H avec dH/dt=-1/mu0*rotE
Inserer les terme de source H=H+source
Es-ce l'idée ?
Si vous vous y connaissez bien en FDTD on pourra aussi discuter par mp.
Mon objectif est de faire un code libre simulant,
incidence oblique quelconque, gaussien, TM... (j'ai aussi des questions).
impactant réseau diélectrique ou métallique de forme quelconque.
Merci pour l'aide.
Cordialement,
-----