2D FDTD avec onde plane d'incidence (code en matlab)

M

MusicMan

Guest
joint à ce message est un 2D FDTD (mode TE) avec une source d'onde plane

i emprunté le 2d FDTD par susan hagness et changé à la source d'onde plane ...

s'il vous plaît
n'hésitez pas à faire des observations
s'il ya une erreur ;-)
Désolé, mais vous avez besoin de login pour afficher cette pièce jointe

 
Je tiens à poser au sujet de la mise à la source de l'incident ey_inc (3).Pourquoi devrait au moins de 3 point?et pourquoi quand je tente de passer à un autre point de plus de 3, les résultats toujours mess.
Si
quelqu'un sait,
s'il vous plaît expliquez-moi.
Merci

 
Nous vous remercions de votre code.Je pense que ce code peut supprimer pml région parce TFSF pouvez cacher la scatered domaine.J'ai une question à propos de iinterpolate 1D onde incidente à l'interface 2D à TFSF.Merci.

<img src="http://www.edaboard.com/images/smiles/icon_biggrin.gif" alt="Very Happy" border="0" />
 
khalmusa a écrit:

Je tiens à poser au sujet de la mise à la source de l'incident ey_inc (3).
Pourquoi devrait au moins de 3 point?
et pourquoi quand je tente de passer à un autre point de plus de 3, les résultats toujours mess.

Si quelqu'un sait, s'il vous plaît expliquez-moi.

Merci
 
je vous remercie de votre code, je tiens à vous demander quand je changer le fond epsilon n (aucune valeur), raison pour laquelle la propagation
n'est pas la même que dans le cas de l'espace libre?

 
epsilon spécifiques avez-vous changé?peut-être il serait mieux que vous pouvez télécharger le code et de mettre en lumière ce que vous avez changé une partie spécifique = 0

 
salut

y at-il un travail de traduction scilab?

i ont scilab, matlab pas.J'ai essayé de traduire le fichier, il a été un succès, mais je ne peux pas exécuter le code traduit, d'une erreur ...

 
Je voudrais poser une question.
i faire exactement comme le livre writed par sullivan.Je me suis un peu de réflexion quand la vague a frappé la frontière dans 1D.but l'auteur a dit il
n'y a pas quelque chose de réflexion.J'espère que
quelqu'un peut me dire why.thanks a lot!

 
Je voudrais poser une question.
i faire exactement comme le livre writed par sullivan.Je me suis un peu de réflexion quand la vague a frappé la frontière dans 1D.but l'auteur a dit il
n'y a pas quelque chose de réflexion.J'espère que
quelqu'un peut me dire why.thanks a lot!

 
si vous utilisez la PML frontière.s'il vous plaît vérifier l'interface de la région de calcul et de couches PML.
si vous utilisez la simple condition à la limite, mai être pas la magie
n'a pas d'équivalent dans le cas de diélectriques fond.
d'espace libre, la magie
n'est pas égale.

 
la vague de transmettre dans l'espace libre. merci pour le même!

 
Je suis seulement à apprendre FDTD et
c'est trés utile pour moi.

 
pourriez-vous me dire comment peut-on définir la direction de l'onde plane oblique incident

 
Je
n'ai pas regardé ou MATLAB code ci-joint, mais je mai être en mesure de répondre à certaines des questions posées dans ce débat.Je
m'excuse si l'un de mes réponses ne
s'appliquent pas à l'attaché de code.

1.Pas de reflets dans 1D FDTD?
Pour le cas particulier de la FDTD 1D il existe une condition à la limite exacte à partir de laquelle vous ne devriez pas obtenir de réflexions.Cela est possible parce que l'énergie ne peut être voyagent dans une direction et toutes les fréquences ont la même vitesse.Grosso modo, cette frontière fonctionne en stockant les valeurs de champ ultrapériphériques pour deux pas de temps et de les utiliser comme terrain de la frontière valeurs de champ à l'heure actuelle pas.

2.Oblique incidence dans FDTD?
Il
s'agit d'une question délicate.Si vous utilisez des conditions aux limites périodiques et qui souhaitent utiliser une source d'impulsion à l'incidence oblique, il
n'ya pas de bonne façon de faire.Il existe une technique appelée la "méthode de mise à jour d'angle", mais cela est limité dans les angles, il peut mettre en œuvre.Sinon, vous êtes bloqué avec seulement une incidence normale.Si vous pouvez faire à votre problème d'espace couvre plus d'une période dans sa section, il y aura d'autres angles discrète à votre disposition.Vous pouvez calculer ce que les angles en utilisant l'équation de la grille.
Si vous
n'avez pas de limites périodiques (PML par exemple) ou si vous utilisez une source sinusoïdale, vous pouvez facilement intégrer l'incidence oblique.Je vous suggère d'utiliser la technique pour ce total-field/scattered-field, mais d'autres approches existent, tels que la dispersion
en champ seule méthode.L'utilisation de TF / SF, vous pouvez régler toutes vos PML régions à la dispersion
en champ.Vous mai aussi un peu plus de cellules disséminées dans le champ si vous souhaitez plus facilement calculer le reflète l'énergie de ce que vous modélisation.

Méthodes alternatives
================
Vous voulez également mai pour examiner les différences finies domaine fréquence (FDFD) method.C'est excellent pour les problèmes 2D ou modélisation des appareils qui sont très résonance où FDTD est moins efficace.C'est la méthode la plus facile je sais de mettre en œuvre.J'enseigne un cours sur la méthode et ont fait de mes notes de cours et des exemples de codes disponibles gratuitement à partir de www.kraetonics.com.Il est d'autant MATLAB code et il peut modèle modes TE et TM, à tout angle d'incidence que vous souhaitez.
L'exemple des codes de modèle un polariseur, d'un cristal photonique, et d'un guide-mode de résonance du filtre.L'idée est simple à mettre en place des exemples afin que vous puissiez facilement les modifier au modèle que vous désirez.

Hope this helps!
-Tip

 
Citation:

Vous voulez également mai pour examiner les différences finies domaine fréquence (FDFD) method.
C'est excellent pour les problèmes 2D ou modélisation des appareils qui sont très résonance où FDTD est moins efficace.
C'est la méthode la plus facile je sais de mettre en œuvre.
J'enseigne un cours sur la méthode et ont fait de mes notes de cours et des exemples de codes disponibles gratuitement à partir de www.kraetonics.com.
Il est d'autant MATLAB code et il peut modèle modes TE et TM, à tout angle d'incidence que vous souhaitez.
 
Par exemple, supposons que vous avez un système qui ne peut stocker trois chiffres de précision et les deux numéros que vous avez souhaité ajouter sont les suivantes:

a = 2,31 e 0
b = 5,62 e -2

La vraie réponse devrait être

a b = 2,31 0,00562 = 2,31562

mais, depuis votre ordinateur stocke seulement trois chiffres, la réponse de votre ordinateur vous est

a b = 2,31

Cette opération a complètement négligé l'essentiel de la variable b.

Si vous mettez en oeuvre FDTD, vos équations sont itérées des milliers ou des millions de fois pour cette erreur
s'accumule au cours de la simulation.

Je pense que 32-bit floating point numbers, 7 chiffres de précision 64-bit nombre a quelque chose comme 15.

En électromagnétisme, les champs E et H sont liés par une impédance de la matière par le rapport E / H = Z. Pour l'espace libre Z = 377 ohms si E et H sont trois ordres de grandeur différente.Chaque fois que vous mettez en place une équation qui ajoute ou soustrait, vous sont en train de perdre de précieux chiffres de précision.Le problème
n'est pas aussi mauvais que cet argument semble indiquer
qu'il existe d'autres facteurs d'échelle dans l'équation que rendre la situation meilleure.Le point, cependant, est que chaque fois que vous vous présenter quelque chose de calculer l'erreur.Pendant des milliers ou des millions de calculs, cette erreur peut devenir assez importante.

L'erreur est encore pire si vous
n'avez pas de normaliser autant que possible.FDTD est une technique particulièrement robuste et vous pouvez sortir avec apparemment beaucoup de choses.Dans la modélisation numérique, il est toujours de bonne pratique de normaliser autant que possible.

-Tip

 
Citation:a = 2,31 e 0

b = 5,62 e -2La vraie réponse devrait êtrea b = 2,31 0,00562 = 2,31562mais, depuis votre ordinateur stocke seulement trois chiffres, la réponse de votre ordinateur vous esta b = 2,31Cette opération a complètement négligé l'essentiel de la variable b.

 
You must log in or register to reply here.

Welcome to EDABoard.com

Sponsor

Back
Top