[求助]FDTD 连接边界

离线wh0802

UID ：19837

注册：2008-10-22
登录：2024-05-04
发帖：54
等级：仿真一级

0楼发表于: 2009-07-17 21:08:06

请教一下编程时连接边界处的场值修改是怎样的，哪里可以查到理 ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册

离线baidu

111

UID ：6612

注册：2007-12-13
登录：2024-11-16
发帖：49
等级：仿真新人

1楼发表于: 2009-07-17 23:12:23

bbs.rfeda.cn/read-htm-tid-12293.html

请牛人指教一下，三维程序中均匀平面波入射时连接边界这样加入有问题吗？

里面的计算过程根据葛老师编的书上公式！（前4个式子计算E0，后面两个坐标转换，不知道这样对不对）请教牛人指导，磁场也按同样方法加入！
ie为x方向格子数+1，je，ke分别为y、z方向
left_x(i,1,k)存储左侧面ex值，其他的类似！

for i=1:ie
        for k=1:ke
            temp1(1)=i*sin(theta)*cos(psi)+jets*cos(theta)*sin(psi)+k*cos(psi);
            temp1(2)=floor(temp1(1));                                %左表面  电场
            temp1(3)=temp1(1)-temp1(2);
            temp1(4)=temp1(2)+1;
            left_x(i,1,k)=(1-temp1(3))*ein(temp1(2))+temp1(3)*ein(temp1(4));
            left_x(i,1,k)=left_x(i,1,k)*(-sin(psi)*sin(alpha)+cos(theta)*cos(psi)*cos(alpha));
            left_z(i,1,k)=-left_x(i,1,k)*sin(theta)*cos(alpha);

            temp1(1)=i*sin(theta)*cos(psi)+jete*cos(theta)*sin(psi)+k*cos(psi);
            temp1(2)=floor(temp1(1));                                %右边面  电场
            temp1(3)=temp1(1)-temp1(2);
            temp1(4)=temp1(2)+1;
            right_x(i,1,k)=(1-temp1(3))*ein(temp1(2))+temp1(3)*ein(temp1(4));
            right_x(i,1,k)=right_x(i,1,k)*(-sin(psi)*sin(alpha)+cos(theta)*cos(psi)*cos(alpha));
            right_z(i,1,k)=-right_x(i,1,k)*sin(theta)*cos(alpha);
        end
    end

    for i=1:ie
        for j=1:je
            temp1(1)=i*sin(theta)*cos(psi)+j*cos(theta)*sin(psi)+kets*cos(psi);
            temp1(2)=floor(temp1(1));                                %下边面  电场
            temp1(3)=temp1(1)-temp1(2);
            temp1(4)=temp1(2)+1;
            bottom_x(i,j,1)=(1-temp1(3))*ein(temp1(2))+temp1(3)*ein(temp1(4));
            bottom_x(i,j,1)=bottom_x(i,j,1)*(-sin(psi)*sin(alpha)+cos(theta)*cos(psi)*cos(alpha));
            bottom_y(i,j,1)=bottom_x(i,j,1)*(cos(psi)*sin(alpha)+cos(theta)*sin(psi)*cos(alpha));

            temp1(1)=i*sin(theta)*cos(psi)+j*cos(theta)*sin(psi)+kete*cos(psi);
            temp1(2)=floor(temp1(1));                                %上边面  电场
            temp1(3)=temp1(1)-temp1(2);
            temp1(4)=temp1(2)+1;
            top_x(i,j,1)=(1-temp1(3))*ein(temp1(2))+temp1(3)*ein(temp1(4));
            top_x(i,j,1)=top_x(i,j,1)*(-sin(psi)*sin(alpha)+cos(theta)*cos(psi)*cos(alpha));
            top_y(i,j,1)=top_x(i,j,1)*(cos(psi)*sin(alpha)+cos(theta)*sin(psi)*cos(alpha));
        end
    end

    for j=1:je
        for k=1:ke
            temp1(1)=iets*sin(theta)*cos(psi)+j*cos(theta)*sin(psi)+k*cos(psi);
            temp1(2)=floor(temp1(1));                                %后边面  电场
            temp1(3)=temp1(1)-temp1(2);
            temp1(4)=temp1(2)+1;
            back_y(1,j,k)=(1-temp1(3))*ein(temp1(2))+temp1(3)*ein(temp1(4));
            back_y(1,j,k)=back_y(1,j,k)*(cos(psi)*sin(alpha)+cos(theta)*sin(psi)*cos(alpha));
            back_z(1,j,k)=-back_y(1,j,k)*sin(theta)*cos(alpha);

            temp1(1)=iete*sin(theta)*cos(psi)+j*cos(theta)*sin(psi)+k*cos(psi);
            temp1(2)=floor(temp1(1));                                %前边面  电场
            temp1(3)=temp1(1)-temp1(2);
            temp1(4)=temp1(2)+1;
            front_y(1,j,k)=(1-temp1(3))*ein(temp1(2))+temp1(3)*ein(temp1(4));
            front_y(1,j,k)=front_y(1,j,k)*(cos(psi)*sin(alpha)+cos(theta)*sin(psi)*cos(alpha));
            front_z(1,j,k)=-front_y(1,j,k)*sin(theta)*cos(alpha);
        end
    end

共1条评分

gwzhao

rf币 +3

积极参与

2009-07-29

1111

离线cem-uestc

UID ：9061