1026阅读
3回复

[讨论]FDTD三维程序主解释

离线roy_1

FDTD好玩

UID ：80037

注册：2011-07-10
登录：2011-07-22
发帖：15
等级：仿真新人

0楼发表于: 2011-07-11 14:52:14

关于FDTD三维的程序，请大神给点注解，小弟感激不尽。
#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#define ie 26
#define je 26
#define ke 26
main()
{
float gax[ie][je][ke],gay[ie][je][ke],gaz[ie][je][ke];
float dx[ie][je][ke],dy[ie][je][ke],dz[ie][je][ke];
float ex[ie][je][ke],ey[ie][je][ke],ez[ie][je][ke];
float hx[ie][je][ke],hy[ie][je][ke],hz[ie][je][ke];
int l,n,i,j,k,ic,jc,kc,nsteps,npml;
float ddx,dt,T,epsz,pi,epsilon,sigma,eaf;
// float xn,xxn,xnum,xd,curl_e;
float t0,spread,pulse;
FILE *fp,*fopen();
ic = ie/2;
jc = je/2;
kc = ke/2;
ddx = 0.01;
dt = ddx/6e8;
epsz = 8.8e-12;
pi = 3.14159;
// printf("jc=%d \n",jc);
for(k=0;k<ke;k++)
{
  for(j=0;j<je;j++)
  {
   for(i=0;i<ie;i++)
   {
    ex[j][k] = 0.0; ZPFTNwf
    ey[j][k] = 0.0; k!rz8S"
    ez[j][k] = 0.0; l0o_C#"<S
    dx[j][k] = 0.0; H*<dte<
    dy[j][k] = 0.0; e;\c=J,eE
    dz[j][k] = 0.0; mjc:0hH
    hx[j][k] = 0.0; AE~}^(G`
    hy[j][k] = 0.0; +#9 (T
    hz[j][k] = 0.0; ?Cl"jcQ*
    gax[j][k] = 1.0; Unk+@$E&
    gay[j][k] = 1.0; k82LCV+6
    gaz[j][k] = 1.0; |bUmkw
   } ;f*xOdi*k
  } u>@G:kt8
} ^DJU99
for(k=7;k<21;k++) uQ)]g
{ eHIC'b.
  gaz[ic][jc][k] = 0.0; $^GnY7$!>
} Qk?Jy<Ra
gaz[ic][jc][kc] = 1.0; &WS%sE{p_
t0 = 20.0; =,9'O/br
spread = 6.0; ~)(Dm+vZ
T = 0; *Ux"3IXO
nsteps = 1; 1.CYs<
while (nsteps>0) l0)6[yXK
{ ,\"gN5[$(
  printf("nsteps-->"); yFv3>\
  scanf("%d",&nsteps); yx4c+(J^8
  printf("%d \n",nsteps); < ,n4|z)
  for(n=1;n<=nsteps;n++) VNfx>&`
  { A}O9e
   T = T+1; ]P5u:~U
   SJe;T
  //calculate the Dx field KO7cZME
  for(k=1;k<ke;k++) P]^] T}5
  { "=O)2}
   for(j=1;j<je;j++) 2j|Eh
   { 6jA Q
    for(i=1;i<ie;i++) ObnB6ShKi
    { m\Nc}P_"p
     dx[j][k] = dx[j][k] +0.5*(hz[j][k] - hz[j-1][k] - hy[j][k] + hy[j][k-1]); *8+YR
    } rZ7 Ihof
   } 'JVvL
  } ->UrWW^
  //calculate the Dy field 6UTdy1Qq>
  for(k=1;k<ke;k++) efm<bJB2
  { T9yW# .
   for(j=1;j<je;j++) }2"k:-g
   { UNHHzTsr?
    for(i=1;i<ie;i++) H|?`n uiD
    { MR: {Ps&,
     dy[j][k] = dy[j][k] +0.5*(hx[j][k] - hx[j][k-1] - hz[j][k] + hz[i-1][j][k]); 4T-AWk
    } ;anG F0x
   } Qmn5-yiw1d
  } cgAcAcmY
   //calculate the Dz field ^%.<(:k[L
  for(k=1;k<ke;k++) '-qc\6UY
  { ?U[AE -*
   for(j=1;j<je;j++) %kT:"j(xW
   { 0>Nq$/!
    for(i=1;i<ie;i++) j/\XeG>
    { irS62Xe
     dz[j][k] = dz[j][k] +0.5*(hy[j][k] - hy[i-1][j][k] - hx[j][k] + hx[j-1][k]); jHzb,&
    } j=LF1dG"
   } "a7d`l:
  } n9yxZu
otx7J\4
  pulse = exp(-0.5*(pow((t0-T)/spread,2.0))); .Dz /MSl
  dz[ic][jc][kc] = pulse; mB`r6'#=
//  dz[0][0][0] = pulse; YFY)Z7fK
  printf("******************************************************\n"); 4)z](e$
  printf("pulse=%6.3f\n",dz[ic][jc][kc]); wXnt3)e
  printf("******************************************************\n"); lw{|~m5`
1lM0pl6M
   //calculate the E from D field bzS [X
  for(k=1;k<ke-1;k++) *!kg@ _0K
  { VdgPb (
   for(j=1;j<je-1;j++) a6OT2B
   { R _%pR_\
    for(i=1;i<ie-1;i++) 1d6pQ9 N
    { -0{WB(P
     ex[j][k] = gax[j][k] * dx[j][k]; ?u|g2!{_
     ey[j][k] = gay[j][k] * dy[j][k]; TM;)[R@
     ez[j][k] = gaz[j][k] * dz[j][k]; .N2Yxty8>
    } E'}$'n?:
   } .[! ^L
  } .nN7*))Fj
-1:asM7
   //calculate the HX field TMrmyvv
  for(k=1;k<ke-1;k++) xUUp?]9y
  { pY@+.V`a
   for(j=1;j<je-1;j++) J^T66}r[f,
   { &R]G)f#w%*
    for(i=1;i<ie;i++) 9"e!0Q40
    { K.SHY!U}
     hx[j][k] = hx[j][k] +0.5*(ey[j][k+1] - ey[j][k] - ez[j+1][k] + ez[j][k]); ~Q0}>m,S
    } n1|%xQBU@
   } ]SI`fja/
  } eAj}/2y"
P:+:Cm<
   //calculate the HY field N C&1l]
  for(k=1;k<ke-1;k++) RPu-E9g@
  { Fk\xq`3'c
   for(j=1;j<je;j++) z&8#1'
   { Q\Eq(2p
    for(i=1;i<ie-1;i++) bhRpYP%x
    { (46)v'?
     hy[j][k] = hy[j][k] +0.5*(ez[i+1][j][k] - ez[j][k] - ex[j][k+1] + ex[j][k]); ?5D7n"jY
    } 9azPUf) C
   } e;$s{CNo
  } ORO~(%-(e
   //calculate the HZ field $A,=z
  for(k=1;k<ke;k++) Yp8$0KK
  { *MNY1+RJ
   for(j=1;j<je-1;j++) D{mu2'q
   { n9%rjS$
    for(i=1;i<ie-1;i++) -Y6JU
    { FVHL;J]nf1
     hz[j][k] = hz[j][k] +0.5*(ex[j+1][k] - ex[j][k] - ey[i+1][j][k] + ey[j][k]); ~H.;pJ{ 8
    } \[BnAgsF
   } 0;9LIL5
  } B.o&%5dG
  }//end for R?(j#bk
  printf("EZ \n"); MD,-<X)Qy
  for(k=1;k<ke;k++) Gu-Sv!4p
  { K(?7E6\vO
   printf("%2d \n",k); B4?P"|
   printf("\n"); W*0KAC`m
   .. z'm;H{xf
[3s~Z8 pP
未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线lx84

呵呵

UID ：30566

注册：2009-04-22
登录：2020-09-16
发帖：866
等级：积极交流六级

1楼发表于: 2011-07-14 14:29:15

这个写得挺清楚的不知道楼主要什么注解？是不是得具体些。

继续！

离线roy_1

FDTD好玩

UID ：80037

注册：2011-07-10
登录：2011-07-22
发帖：15
等级：仿真新人

2楼发表于: 2011-07-16 15:59:40

我想要一些关于程序的解释，如
for(k=7;k<21;k++)
{
gaz[ic][jc][k] = 0.0;
}
gaz[ic][jc][kc] = 1.0;
为什么要从7开始。还有为什么要引入D等等。
不知道你是否懂，能不能给点提示。感激不尽。

离线jiandan1986

UID ：77042

注册：2011-05-09
登录：2012-04-16
发帖：25
等级：仿真新人

3楼发表于: 2011-07-24 19:36:16

回 2楼(roy_1) 的帖子

从7开始应该是和边界条件有关系吧，D好像是公式中的一个系数还是什么的，具体忘记了，你可以参考一下其它人分享的资料，里面有解释的。

发帖回复