微波EDA网,见证研发工程师的成长!
首页 > 研发问答 > 微波和射频技术 > 电磁仿真讨论 > 3D FDTD model with a second order PML in Matlab

3D FDTD model with a second order PML in Matlab

时间:03-31 整理:3721RD 点击:
I'm have developed a 3D FDTD model with a second order ABC in matlab.

I wish to replace the ABC with a PML.

Does anyone have some sample matlab code for this?

Thanks,
Martin

Anyone?

Has anyone a complete set of Tafloves and hagnesses code. Im getting the book next week, but I'd like to get a start with the code...


Thanks

sorry, what is PML?

May I find any refrence for PML on the edaboard

Hi,

here is the MatLab code from A. Tafloves book (2ed edition);
it are 3 examples for 1d, 2d and 3d FDTD.

best: ricy

I've downloaded this code but in the 3D example there is not the PML condition. In FDTD3D file there is the simulation of a rectangular cavity resonator with PEC boundaries:
ex(i,j,k)=0 on the j=1, j=jb, k=1, and k=kb planes
ey(i,j,k)=0 on the i=1, i=ib, k=1, and k=kb planes
ez(i,j,k)=0 on the i=1, i=ib, j=1, and j=jb planes

There's an example with 3D PML condition?

Please help me!

Pippo!

can you send me this program
and i will send you the programes in PML
email: ahmedsk_zf81@yahoo.com

who can provide me with 3d pml formula in brief

Now i'm trying to write 3D PML code but I have a ploblem in the PMl region (between main and PMl) And I'm confuse in coefficient in pml
Does anyone have some sample 3DPML matlab code or brief of method to write it?
Thank you

how do you incorporate a plane wave incident field in FDTD?

You can use the 2D implementation in the 3D UPML implementation.
What I mean is that you could do a kind of copy-paste from the 2D PML to the 3D UPML and get a 3D PML working code.
Hope this help
regards

The code is very usefull for my project

With regards
vasi

hello, i am working in a project of a master of science and I need FDTD 3D to simulate propagation in an indoor environment. I have no longer time, i have problems when i integrate PML please could some one help me?

I will be really very gratefull.

Dear Musicman,

A plane wave can be implemented in the FDTD code through either of the two ways:
1. Scattered Field approach given in Kunze and Luebber's text book (it has a working fortran code also-FDTDA.for)
2. The Total Field Scattered Field (TFSF) approach given in Taflove's Book. It is a bit complicated.

PML (Perfectly Matched Layer) is an absorbing material boundary condition that first introduced by J.P. Berenger.

J.P. Berenger, ?A perfectly matched layer for the absorption of electromagnetic waves,? Jour. Comp. Phys., vol.114, Aug. 1994, pp.185-200

can anyone give me an example with 2.5D(which is BOR-body of revolution) FDTD? thank you

Added after 5 minutes:

btw thanks for the reply to my first question about TF/SF field. i was able to do it in 2D but i am having problems implementing it in 2.5D(or BOR) FDTD(with pml ABC). if anyone knows pls help me. thanks a lot =)

One of the major problems in the standard FDTD is that the requirement for artificial mesh truncation condition. These aritficial termination conditions are known as absorbing boundary condition(ABCs) as they theoretically absorb incident and scattered fields. The accuracy of the ABCdictates the accruracy of the FDTD method.
The need for accurate ABCs has resulted in various types of ABCs, Berenger's perfectly matched layer(PML) is the most widely accepted.
In the perfectly matched layer truncation technique, an artificial layer of absorbing material is placed around the outer boundary of the computational domain. The goal is to ensure that a plane wave that is incident from FDTD free space to the PML region at an arbitrary angle is complete transmission of the incident plane wave at the interface between free space and the PML region. Thus the FDTD and the PML region are said to be perfectly matched.

The 3D FDTD codes I've seen around here so far do not include (Berenger's) PML conditions.
Does anyone have the right code (or at least the steps for implementing the PML conditions) for the 3D case?

Thanks

You could take my implemented PML's (Split-Field, UPML, CPML, NPML) for 2D TMz from https://www.edaboard.com/viewtopic.p...highlight=winn and understand the working principle, after that it should be no problem to extend the scheme to 3D.

I kindly request u to forward 3D FDTD model with a second order ABC in matlab code to my mail. My ID is ravidurbha@gmail.com.

Thanx in Advance,
Ravi Durbha

Copyright © 2017-2020 微波EDA网 版权所有

网站地图

Top