Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronics Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

3D FDTD model with a second order PML in Matlab

Status
Not open for further replies.

king_of_inismac

Newbie level 2
Joined
Mar 7, 2006
Messages
2
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,314
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
 

king_of_inismac

Newbie level 2
Joined
Mar 7, 2006
Messages
2
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,314
perfectly matched layer matlab code

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
 

media

Newbie level 3
Joined
Mar 10, 2006
Messages
3
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,300
pml condition three dimension

sorry, what is PML?
 

fasaghandis

Member level 5
Joined
Jul 10, 2002
Messages
81
Helped
1
Reputation
2
Reaction score
0
Trophy points
1,286
Location
tehran
Activity points
638
matlab pml

May I find any refrence for PML on the edaboard
 

ricy

Member level 1
Joined
Feb 1, 2005
Messages
37
Helped
12
Reputation
24
Reaction score
3
Trophy points
1,288
Activity points
315
3d pml implementation

Hi,

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

best: ricy
 

pippopamp

Newbie level 4
Joined
Jun 26, 2007
Messages
6
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,304
matlab plane wave 3d

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!!
 

ahmedsk

Newbie level 5
Joined
Jul 13, 2007
Messages
9
Helped
1
Reputation
2
Reaction score
0
Trophy points
1,281
Activity points
1,330
derivation of the berenger scheme for pml

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

shenhongj2004

Member level 1
Joined
Aug 20, 2006
Messages
35
Helped
2
Reputation
4
Reaction score
1
Trophy points
1,288
Activity points
1,468
pml conditions

who can provide me with 3d pml formula in brief
 

sangwoo

Junior Member level 1
Joined
May 21, 2007
Messages
16
Helped
2
Reputation
4
Reaction score
2
Trophy points
1,283
Activity points
1,388
3d plane matlab

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
 

MusicMan

Member level 2
Joined
Apr 18, 2007
Messages
42
Helped
9
Reputation
18
Reaction score
4
Trophy points
1,288
Activity points
1,527
pml book by j.p.berenger

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

grendhell

Member level 1
Joined
May 26, 2007
Messages
40
Helped
3
Reputation
6
Reaction score
1
Trophy points
1,288
Activity points
1,608
plane wave and pml - matlab

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
 

vasikara

Junior Member level 2
Joined
Sep 11, 2007
Messages
22
Helped
6
Reputation
12
Reaction score
6
Trophy points
1,283
Activity points
1,407
a. tafloves

The code is very usefull for my project

With regards
vasi
 

goldenboy1981

Newbie level 1
Joined
Jan 23, 2007
Messages
1
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,285
pml lame matlab

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.
 

jithesh

Member level 5
Joined
Oct 11, 2007
Messages
92
Helped
41
Reputation
80
Reaction score
29
Trophy points
1,298
Location
India
Activity points
1,921
example of pml boundary condition

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.
 

scheila

Newbie level 3
Joined
Sep 22, 2006
Messages
4
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,304
pml implementation matlab

media said:
sorry, what is PML?

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
 

MusicMan

Member level 2
Joined
Apr 18, 2007
Messages
42
Helped
9
Reputation
18
Reaction score
4
Trophy points
1,288
Activity points
1,527
3d pml

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 =)
 

llsgz85

Newbie level 2
Joined
Nov 28, 2007
Messages
2
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,295
pml 3dimentional

media said:
sorry, what is PML?
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.
 

bfcgg7

Newbie level 1
Joined
Feb 11, 2008
Messages
1
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,285
pml matlab code

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
 

Winn

Newbie level 5
Joined
Aug 21, 2008
Messages
10
Helped
3
Reputation
6
Reaction score
1
Trophy points
1,283
Activity points
1,353
fdtd+pml+these+tmz

You could take my implemented PML's (Split-Field, UPML, CPML, NPML) for 2D TMz from and understand the working principle, after that it should be no problem to extend the scheme to 3D.
 

fireblade

Newbie level 6
Joined
Nov 1, 2006
Messages
13
Helped
2
Reputation
4
Reaction score
2
Trophy points
1,283
Activity points
1,356
Re: pml matlab

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
 

Status
Not open for further replies.

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Top