Re: cylinder in FDTD
ravi if my program are rigth plz send me a visible men head model i asked you this in this topic
the program in Matlab you can to conver it in C++ if you use it
the speace for simulation are borned by : ia,ib;ja,jb;ka,kb (sulivan notation) dont forget PML
ic,jc,kc are the center cylinder
h=height_cylinder;
radius= radius_cylinder;
program :
clc
t=1000;
epsi0=8.8419e-12;
N=124;
ia=9;
ja=9;
ka=9;
ic=N/2;
jc=N/2;
kc=N/2;
ib=N-ia+1;
jb=N-ja+1;
kb=N-ka+1;
ddx=0.0025;
dt=ddx/6e8;
epsilon=44;
sigma=0.9;
for j=1:1:N
for k=1:1:N
for i=1:1:N
ga(i,j,k)=1;
gb(i,j,k)=0;
end
end
end
h=40;
radius= 10;
for i=ia:1:ib
for j=ja:1:jb
for k=ka:1:kb
xdist=(ic-i);
ydist=(jc-j);
kdist=abs(kc-k);
dist=sqrt(xdist^2+ydist^2);
if (dist<=radius)&&(kdist<=h/2)
ga(i,j,k)=1/(epsilon+(sigma*dt/epsi0));
gb(i,j,k)=sigma*dt/epsi0;
end
end
end
end
i test it its work
it is not necessery simulate coaxial. you simulate juste your antenna and you fixe the E filed or D filed directly in the cell wher the coaxial connected that is wat the people do
this document will help you
bat remenber to send me the visible man head please
email :
tassoki@yahoo.fr
or post it in this topic plz
you use matlab or c++ or other
if you have any difficulte to convert the program on c++ i will help you
excuse my english i dont speak it