FDTD MUR = Instabilité heure tardive (résolu)

S

somada141

Guest
Eh bien bonjour à tous,

J'ai posté deux fois demander aux gens s'ils savaient quelque chose au sujet instabilités causées par la combinaison de MUR FDTD.

Bien que certaines personnes ont répondu et je les remercie pour leur époque, personne ne savait vraiment pourquoi y aurait-il des instabilités, mais plutôt proposé d'utiliser la PML au lieu de MUR.

Après avoir lutté pendant un mois avec Mur (et bien que je pensais que je résolu mon problème comme indiqué dans un précédent post, qui du lien est:

http://www.edaboard.com/viewtopic.php?t=333674&highlight =

Check it Out il clarifie le problème)

J'ai découvert quel est le problème.

Il ne suffit pas ABC MUR sont vraiment faibles comparant à d'autres tels que ABC Liao ou les amoureux de la couche PML elle a plus trop de problèmes.

Lorsqu'il est combiné avec FDTD Mur souffre ce qu'on appelle le temps d'instabilité ultérieure qui est exactement mon problème (voir le post ci-dessus) où, après un certain temps il ya une explosion des valeurs sur certaines des limites du domaine de calcul.

Cependant j'ai trouvé quelques articles sur le net sur ce point et surtout un grand donnant une solution ou alternative à la réparation Mur.

Le titre du journal est la suivante:Les études de la stabilité des FDTD avec la condition aux limites absorbante Mur de second ordre en 3-D des problèmes de diffusion.Bien que je n'ai pas essayé leur solution sur mon code encore (je vais bien), il semble fonctionner très bien et sans avoir à apporter des changements majeurs dans votre code existant déjà.

J'ai posté ce parce que j'espère que quelqu'un va être aidé par la présente.

Vive FDTD

PS.Le papier est un peu inaccessible mais j'ai reçu d'un ami avec un abonnement grand.S'il vous plaît me dire si je peux le télécharger ici (je n'ai pas pris la peine de lire les règles).Si vous avez besoin de contact par e-mail (somada141 (at) yahoo.gr) et je peux vous l'envoyer.PS2.
Il semble que je besoin d'un moyen de savoir quand mon code a atteint l'état d'équilibre ... Est-ce que quelqu'un sait comment faire cela?
Aide si vous pouvez<img src="http://www.edaboard.com/images/smiles/icon_biggrin.gif" alt="Very Happy" border="0" />

Ajouté après 10 minutes:Il y aura des mises à jour sur ce que dès que je tente cette méthode ....
Dernière édition par somada141 le 27 novembre 2008 23:32, édité 1 fois au total

 
Salut je peux avoir le papier? = pi besoin Parce que je suis en train de simuler une structure à long
Adresse e-mail est mon brrrrrgrrrrr (at) yahoo.com (5r de ^ _ ^)

 
Mes bonnes gens,

Vous verrez le miracle de l 'Anti - LTI comme je l'ai appelé cette chose (Anti-heure tardive instabilité).

Les gens qui ont vu mes messages au cours du mois dernier saura j'ai eu un énorme proplem avec les instabilités causées par l'utilisation à long terme du 2e ordre ABC MUR (s'il vous plaît lire les messages respectifs).

Toutefois Zhang et Wang Yusheng Wenbing m'a montré le chemin.

Le document ci-dessus a la solution pour vous.Je suggère fortement de vous causer à cela fonctionne vraiment.

Ok Let's Get est sérieuse nous?

J'ai passé la journée à mettre en œuvre et de tester le rapprochement donnée par ce document.

Ils utilisent le rapprochement vent pour descritizing la manière des équations d'ondes et un Mur utilisés dans ABC du 2e ordre.L'équation qui se produit fonctionne mieux que je pouvais espérer.

Toutefois photos parlent mieux que les mots.Pendant tout ce temps je me débattais pour obtenir le MUR de travailler J'étais surtout des essais avec une structure de guide d'ondes ouvert sur deux se termine par un hard-source d'excitation (voir le post précédent).

L'image ci-dessous est tirée de la simulation à No300 pas de temps.Aucun anti-LTI a été utilisée ici.Cependant, il ya de 2e ordre MUR sur les faces intérieures et le 1er MUR ordre sur les bords et des coins (voir Kunz)<img src="http://images.elektroda.net/98_1227820946.gif" border="0" alt="FDTD MUR = Late Time Instability (SOLVED)" title="FDTD MUR = Instabilité heure tardive (résolu)"/>Vous pouvez voir tout ce qui est allé en enfer.Jusqu'à l'étape 250 tout a fonctionné à trouver mais le coup de pied dans le LTI et le résultat est indiqué ci-dessus.Par étape 300, la solution est corrompu.

Ci-dessous est la même simulation, mais avec l'utilisation des anti-LTI.En voici j'utilise 1ère commande Mur sur les bords et les coins comme auparavant, mais l'approche proposée par le papier dans les parties intérieures des visages.Les résultats parlent d'eux-mêmes:<img src="http://images.elektroda.net/53_1227821239_thumb.gif" border="0" alt="FDTD MUR = Late Time Instability (SOLVED)" title="FDTD MUR = Instabilité heure tardive (résolu)"/> Comme vous pouvez le voir tout fonctionne très bien.En dehors de quelques réflexions numériques en raison de la faible distance entre les geomtery et les limites (pour la simulation rapide) qui devraient se produit pas d'instabilité.

Voici ci-dessus après 50 étapes plus:<img src="http://images.elektroda.net/37_1227821333_thumb.gif" border="0" alt="FDTD MUR = Late Time Instability (SOLVED)" title="FDTD MUR = Instabilité heure tardive (résolu)"/> Comme vous pouvez le voir tout fonctionne encore.J'ai simulé plusieurs centaines de marches avant et pas LTI a eu lieu.Conclusion: Avez-vous des problèmes d'instabilité avec ABC Mur?
Si c'est le cas cet article donne la solution.
Impossible d'obtenir le papier?
Demandez-moi et je vais vous l'envoyer dès que possible.Vive FDTD ...

Cheers,
AdamAjouté après 14 minutes:Depuis que je suis heureux de ce travaillé pour moi je partage les équations complet qui est proposé par le papier pour FDTD 3D avec des cellules rectangulaires (ce qui recouvre la donnée cubiques) pour les 6 limites (il peut y avoir des erreurs - si vous voyez tout contact moi plz ):<img src="http://images.elektroda.net/100_1227822247_thumb.gif" border="0" alt="FDTD MUR = Late Time Instability (SOLVED)" title="FDTD MUR = Instabilité heure tardive (résolu)"/>Désolé, mais vous avez besoin de connexion pour voir cette pièce jointe

 
il suffit d'utiliser ces que vous utilisez le classique de 2e ordre Mur.Ils fonctionnent comme un charme!

 
Jusqu'à présent, je pensais que c'était contre les règles du forum pour poster les documents eux-mêmes, mais j'ai rencontré plusieurs personnes qui l'ont fait et sont encore en vie.

C'est pourquoi il est ici.Le papier avec le LTI-Mur anti ABC
Désolé, mais vous avez besoin de connexion pour voir cette pièce jointe

 
Cher modérateur / administrateur

Si ce n'est pas la loi abbiding après tout puis retirez la pièce jointe, mais s'il vous plaît épargnez-moi

 

Welcome to EDABoard.com

Sponsor

Back
Top