645阅读
0回复

[求助]这个程序的边界吸收不好，能否指点一下，如何设置PML边界

I do

UID ：91162

注册：2012-03-29
登录：2013-03-15
发帖：19
等级：仿真新人

0楼发表于: 2012-04-30 17:29:21

关键词：设置

for (n=1;n<=NSTEPS;n++)
        {
            T++;

            for (k=1;k<KE; k++)
            {
                dx[k]=dx[k]+0.5*(hy[k-1]-hy[k]);
            }

            pulse=exp(-.5*(pow((t0-T)/spread,2.0)));    /* put a Gassian pulse at the low end */
            dx[5]=dx[5]+pulse;
            printf("%5.1f  %6.2f %6.2f \n",T,pulse,dx[5]);


            for (k=0;k<KE-1; k++)
            {
                ex[k]=ga[k]*(dx[k]-ix[k]);
            }

            for (k=0;k<KE-1; k++)
            {
                ix[k]=ix[k]+gb[k]*ex[k];
            }


/* Boundary conditions */
            ex[0]=ex_low_m2;
            ex_low_m2=ex_low_m1;
            ex_low_m1=ex[1];

            ex[KE-1]=ex_high_m2;
            ex_high_m2=ex_high_m1;
            ex_high_m1=ex[KE-2];         /* calculate the hy field */
            for (k ..

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

发帖回复

[求助]这个 程序的边界吸收不好，能否指点一下，如何设置PML边界

[求助]这个程序的边界吸收不好，能否指点一下，如何设置PML边界