1d FDTD nestabilus su metalo plona plėvele

[zh168]

Sveiki, aš rašau 1D FDTD imituoti banga skleisti į aliuminio plona plėvele. Deja, kai bangos nukentėjo metalo paviršiaus, E-lauko palaipsniui didėja ir tampa nestabilus. Ar kas nors žino, kas programą ir bet koks sprendimas, kad? Čia žemiau yra parametrai, kad aš naudoju: bangos ilgis: 1um n_Al: 1,35 (matuojama) k_Al: 9,58 (matuojama) tinklo dydis: wavelength/50 time_step: grid dydis / (2 * c) Thanks!

 

[rrumpf]

Nors aš negaliu būti tikras, ką jūsų problema, phas greitis gali viršyti c tam tikromis aplinkybėmis. Tai reikalauja, kad jums peržiūrėtą stabilumo būsena apskaičiuoti stabilų laiko žingsnis. Žr EQ. (4,76) ir aplinkinių diskusiją apie mano disertaciją 103 psl. http://etd.fcla.edu/CF/CFE0001159/Rumpf_Raymond_C_200608_PhD.pdf Hope this helps! Patarimas

 

[iyami]

Tiesą sakant, manau, jūsų problema yra tiesiog, kad jūsų dielektrinė skvarba yra neigiama (raidė epsilon = n ^ 2 - k ^ 2) nurodytą dažnį. FDTD tik stabilus santykinė dielektrinė skvarba> = 1,0. Manau, kad geriausias būdas kovoti su savo problema yra naudoti standartiniai modeliai (Lorenco, Drude, Debye) dispersinis medžiagos (vieno poliaus yra pakankamai) ir atitiktų parametrus ir tada jį įgyvendinti naudojant rekursinis Konvoliucija arba ADE.