3950阅读
23回复

[求助]新手请教Traflove书中附带的2D程序

方恨少

UID ：17098

注册：2008-08-24
登录：2019-01-09
发帖：1374
等级：荣誉管理员

20楼发表于: 2008-12-22 16:04:19

在家不太方便，先解释一点吧。
%     FRONT region
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化
cbexbcf(1:iefbc,1)=0.0;
%为什么ca,cb在PML中的系数这么初始化
因为为了不反射，所以介质是不能发生突变的，突变就意味着有很大的反射。
你可以先看一下真空中的ca,cb是多少。然后再理解为什么这么初始化。
(1:iefbc,1)表示i 从1 到iefbc，j =1，那应该就是模型中最左边的那个竖边了。
for j=2:jebc
  y1=(jebc-j+1.5)*dx; %1。5和0。5是怎么来的啊，这y1,y2又什么参数
  y2=(jebc-j+0.5)*dx;
  sigmay=bcfactor*(y1^(orderbc+1)-y2^(orderbc+1));

1.5和0.5以及y1,y2都在sigmay中用到的，而sigmay是PML的介电系数，所以这可能是个经验公式啥的，因为毕竟在推导PML时，只是说只要满足那个公式就可以了，但电导率具体怎么设置，在PML层中的分布没说，而是只是大概说了一下，一般2次方或三次方就差不多了。
  ca1=exp(-sigmay*dt/(epsz*eps(1)));
  cb1=(1.0-ca1)/(sigmay*dx);
  caexbcf(1:iefbc,j)=ca1;
  cbexbcf(1:iefbc,j)=cb1;
end

共1条评分

cem-uestc

技术分 +2

积极参与讨论+技术分论坛感谢您的参与

2008-12-22

逆流而上

离线wq_463

UID ：20925

注册：2008-11-06
登录：2021-04-22
发帖：227
等级：仿真三级

21楼发表于: 2008-12-22 19:31:08

yRe:新手请教Traflove书中附带的2D程序

引用第20楼gwzhao于2008-12-22 16:04发表的 : 87}&`
在家不太方便，先解释一点吧。 tt%MoQ)
% FRONT region A*./,KT
caexbcf(1:iefbc,1)=1.0; %为什么ca,cb在PML中的系数这么初始化 ;fw1
cbexbcf(1:iefbc,1)=0.0; {X2`&<i6
%为什么ca,cb在PML中的系数这么初始化 _|72r}j
.......

首先谢谢师兄的解答
这个所谓的的经验公式是不是调试程序过程中,自己总结出来的,那么他的可用性有多大呢,就是说别人能不能直接引用呢
还有既然这段程序已经算过,为什么下面这段程序又算了一边ca1,cb1,而且这是算FRONT-PML的为什么把LEFT和RIGHT的一起也算了,而且后面的程序也有重复
sigmay = bcfactor*(0.5*dx)^(orderbc+1);
ca1=exp(-sigmay*dt/(epsz*eps(1)));
cb1=(1-ca1)/(sigmay*dx);
caex(1:ie,1)=ca1;
cbex(1:ie,1)=cb1;
caexbcl(1:iebc,1)=ca1;
cbexbcl(1:iebc,1)=cb1;
caexbcr(1:iebc,1)=ca1;
cbexbcr(1:iebc,1)=cb1;