S
somada141
Guest
Hi everyone,
Jaz sem trenutno bori s svojo tezo projekt, ki ima veze s 3D FDTD koda na odprtih geometrij.
Ker implicitne rabim ABCs simulirati prostora in standardne raztopine za FDTD je
2. Odredba Mur Absorpcijska robni pogoji.
Vendar tudi po obravnavi Tavlove
je knjigo o FDTD (vključno poglavja o ABCs) še vedno sem na izgubo, kako jih uporabljati.
Če želite biti bolj specifični sem programirana moj FDTD kodo brez ABCs in deluje kot čar.Vendar pa je trenutek, ko sem poskusil uveljaviti mejne pogoje, vse gre Nakriv.Ali ima kdo izkušnje s temi oznakami, tako da mi lahko pomaga ali priročnik mi nekaj koristnih (in po možnosti razumljiv) Mur kod?Obstajajo točke na Mur formule, kjer še nisem grapsed idejo.
and where do we put the Mur code in the standard FDTD formulae?Npr, kako smo naj program 3 čas potrebne ukrepe
in, kjer bomo dal Mur kodo v standardnih FDTD formule?Če nekaj lahko pomaga prosim odgovorite tukaj ali pa se obrnite po elektronski pošti na somada141 (at) yahoo.gr
Vsaka pomoč bi cenil,
Adam
PS Tu je čisti FDTD del moje kode (Python):
Koda:
# Hx
za i v xrange (1, Hx.shape [0] -1):
za j pri xrange (1, Hx.shape [1] -1):
za k, v xrange (1, Hx.shape [2] -1):
Hx [i, j, k] = Da * Hx [i, j, k] Db * (EY [i, j, k]-EY [i, j, k-1] Ez [i, j, k ] Ez [i, j-1, k])
# Hy
za j pri xrange (1, Hy.shape [1] -1):
za i v xrange (1, Hy.shape [0] -1):
za k, v xrange (1, Hy.shape [2] -1):
Hy [i, j, k] = Da * hy [i, j, k] Db * (Ez [i, j, k] Ez [i-1, j, k]-Ex [i, j, k ] Ex [i, j, k-1])# Hz
za k, v xrange (1, Hz.shape [2] -1):
za i v xrange (1, Hz.shape [0] -1):
za j pri xrange (1, Hz.shape [1] -1):
Hz [i, j, k] = Da * Hz [i, j, k] Db * (Ex [i, j, k]-Ex [i, j-1, k]-EY [i, j, k ] EY [i-1, j, k])# Ex
za i v xrange (Ex.shape [0]):
za j pri xrange (Ex.shape [1]):
za k, v xrange (Ex.shape [2]):ii = i-1
če (ii ==- 1):
ii = 0Ex [i, j, k] = Ca [ii, j, k] * Ex [i, j, k] Cb [ii, j, k] * (Hz [i, j 1, k]-Hz [ i, j, k]-hy [i, j, k 1] hy [i, j, k])# EY
za i v xrange (Ey.shape [0]):
za j pri xrange (Ey.shape [1]):
za k, v xrange (Ey.shape [2]):jj = j-1
če (jj ==- 1):
jj = 0EY [i, j, k] = Ca [i, jj, k] * EY [i, j, k] Cb [i, jj, k] * (Hx [i, j, k 1]-Hx [ i, j, k]-Hz [i 1, j, k] Hz [i, j, k])# Ez
za i v xrange (Ez.shape [0]):
za j pri xrange (Ez.shape [1]):
za k, v xrange (Ez.shape [2]):kk = k-1
if (kk ==- 1):
kk = 0Ez [i, j, k] = Ca [i, j, kk] * Ez [i, j, k] Cb [i, j, kk] * (hy [i 1, j, k]-hy [ i, j, k]-Hx [i, j 1, k] Hx [i, j, k])
Jaz sem trenutno bori s svojo tezo projekt, ki ima veze s 3D FDTD koda na odprtih geometrij.
Ker implicitne rabim ABCs simulirati prostora in standardne raztopine za FDTD je
2. Odredba Mur Absorpcijska robni pogoji.
Vendar tudi po obravnavi Tavlove
je knjigo o FDTD (vključno poglavja o ABCs) še vedno sem na izgubo, kako jih uporabljati.
Če želite biti bolj specifični sem programirana moj FDTD kodo brez ABCs in deluje kot čar.Vendar pa je trenutek, ko sem poskusil uveljaviti mejne pogoje, vse gre Nakriv.Ali ima kdo izkušnje s temi oznakami, tako da mi lahko pomaga ali priročnik mi nekaj koristnih (in po možnosti razumljiv) Mur kod?Obstajajo točke na Mur formule, kjer še nisem grapsed idejo.
and where do we put the Mur code in the standard FDTD formulae?Npr, kako smo naj program 3 čas potrebne ukrepe
in, kjer bomo dal Mur kodo v standardnih FDTD formule?Če nekaj lahko pomaga prosim odgovorite tukaj ali pa se obrnite po elektronski pošti na somada141 (at) yahoo.gr
Vsaka pomoč bi cenil,
Adam
PS Tu je čisti FDTD del moje kode (Python):
Koda:
# Hx
za i v xrange (1, Hx.shape [0] -1):
za j pri xrange (1, Hx.shape [1] -1):
za k, v xrange (1, Hx.shape [2] -1):
Hx [i, j, k] = Da * Hx [i, j, k] Db * (EY [i, j, k]-EY [i, j, k-1] Ez [i, j, k ] Ez [i, j-1, k])
# Hy
za j pri xrange (1, Hy.shape [1] -1):
za i v xrange (1, Hy.shape [0] -1):
za k, v xrange (1, Hy.shape [2] -1):
Hy [i, j, k] = Da * hy [i, j, k] Db * (Ez [i, j, k] Ez [i-1, j, k]-Ex [i, j, k ] Ex [i, j, k-1])# Hz
za k, v xrange (1, Hz.shape [2] -1):
za i v xrange (1, Hz.shape [0] -1):
za j pri xrange (1, Hz.shape [1] -1):
Hz [i, j, k] = Da * Hz [i, j, k] Db * (Ex [i, j, k]-Ex [i, j-1, k]-EY [i, j, k ] EY [i-1, j, k])# Ex
za i v xrange (Ex.shape [0]):
za j pri xrange (Ex.shape [1]):
za k, v xrange (Ex.shape [2]):ii = i-1
če (ii ==- 1):
ii = 0Ex [i, j, k] = Ca [ii, j, k] * Ex [i, j, k] Cb [ii, j, k] * (Hz [i, j 1, k]-Hz [ i, j, k]-hy [i, j, k 1] hy [i, j, k])# EY
za i v xrange (Ey.shape [0]):
za j pri xrange (Ey.shape [1]):
za k, v xrange (Ey.shape [2]):jj = j-1
če (jj ==- 1):
jj = 0EY [i, j, k] = Ca [i, jj, k] * EY [i, j, k] Cb [i, jj, k] * (Hx [i, j, k 1]-Hx [ i, j, k]-Hz [i 1, j, k] Hz [i, j, k])# Ez
za i v xrange (Ez.shape [0]):
za j pri xrange (Ez.shape [1]):
za k, v xrange (Ez.shape [2]):kk = k-1
if (kk ==- 1):
kk = 0Ez [i, j, k] = Ca [i, j, kk] * Ez [i, j, k] Cb [i, j, kk] * (hy [i 1, j, k]-hy [ i, j, k]-Hx [i, j 1, k] Hx [i, j, k])