FDTD Mur: Kjente stabilitetsproblemer (Vær hjelpe folk!)

[somada141]

Hei alle sammen,

Siden mitt forrige innlegg var sett av 24 personer, og ingen gidder å hjelpe deg allikevel
vil jeg prøve å avklare mitt problem.

Jeg gjør en 3D FDTD kode med Mur ABCs.

Når jeg bare bruke rent FDTD formler alt virker fint, men når jeg aktiverer mur formler på grenser og etter en stund av normal atferd en unormal økning skjer

Her er et øyeblikksbilde:<img src="http://images.elektroda.net/25_1226519582_thumb.gif" border="0" alt=""/>

<img src="http://images.elektroda.net/55_1226520637_thumb.gif" border="0" alt=""/>

<img src="http://images.elektroda.net/71_1226520674_thumb.gif" border="0" alt=""/>

Dette er en 3D gylden waveguide med to åpning på y-aksen.Jeg har bare håndheves Mur på y = 0 for å gjøre forskjellen mellom Mur og fravær klar.

Rammen på y = 0 er 50 celler unna wave guide mening 3wavelengths bort på 8GHz som er en stor frekvens for denne waveguide.

Er det noen som hvem vet hva som kan ha gått galt?Har noen sett noe lignende.

Når det gjelder standard FDTD stabilitetsproblemer mening D (celle størrelse) og
Dt. (tid trinn) er det både dekkes av standard ligning av Allen tavlove.D er 15 ganger mindre enn bølgelengde og
Dt. er produsert av D / (SQRT (3) * C).Vennligst noen hjelpe.
Hvis noen plager som kan hjelpe meg så jeg kan hjelpe noen andre med denne metoden og dens forskjellige mangler og problemer.Takk på forhånd,
Adam

somada141 (at) yahoo.gr

 

[Winn]

Faktisk Mur ABC er en av de enkleste ABC, bør du programmet på en slik måte at:

/ / Abc på begynnelsen av simulering domene
ex [0] = ex_low_m2;
ex_low_m2 = ex_low_m1;
ex_low_m1 = ex [1];

/ / Abc på slutten av simulering domene
ex [nxsize-1] = ex_low_m2;
ex_low_m2 = ex_low_m1;
ex_low_m1 = ex [nxsize-2];

Håper dette hjelper deg!

Hilsen

 

[somada141]

Egentlig fant jeg feilen ble jeg gjør.Det var litt av en dum feil.

Jeg bare legger inn dette i tilfelle noen finner det nyttig:

Mens oppdatere den innerste feltet for E med FDTD slik at du oppdaterer den ytre feltet komponenter med mur du beregne indre fileds bare på skiltene der du håndheve Mur.

For eksempel for Ex:

Vi oppdaterer med FDTD i:

for x i (0, avslutter)
for Y i (1, end-1)
for z i (1, end-1)
FDTD

det jeg gjorde var å regne ut x 1, end-1 ... Derav ustabilitet forårsaket.

Håper dette hjelper noen der ute

 

[medaziz]

hei

Er det mer diskusjon om dette emnet

 

[somada141]

Nei, fra begynnelsen var det praktisk talt ingen interesse i dette emnet fra andre brukere.Men jeg har løst problemet mitt så hvis du har noen lignende problemer, eller hvis du har andre spørsmål vedrørende FDTD & Mur spør ...

 

[medaziz]

hei

Jeg arbeider også med FDTDsince 2003, men mine interset er focalised å studere biomedicals antenne.Sincerly Jeg har vanligvis noen FDTD problemer.

Hjertelig

 

[somada141]

hva slags problemer du møter med FDTD?

 

[somada141]

Denne saken ER RESOLVED.CHECK ut dette:http://www.edaboard.com/viewtopic.php?t=335023&highlight =