[求助]关于大地上方雷电通道的电磁场问题

可以将雷电通道理解为电流通道，我想用fdtd计算雷电通道产生的感应电磁场，但现在有疑问： >z,Y%A  
1、如何将源，也就是电流激励加入fdtd； Cq-99@&;  
2、如何考虑空气和大地相交的边界问题。 x/0x&la  
s"8z q;)  

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

你想求的应该是这个雷电通道旁边的感应电场，对大地那块应该不感兴趣吧。 w *o _s  
如果要考虑反射的话，你得找一下资料，看一下土壤的介电系数，电阻率啥的是多少，当作一个导体好了。 ;!m_RQPFF  
这个雷电的电流到大地的时候，是不是可以理解为全部进入大地了？也就是不需要考虑反射了。 mn\GLR.  
>EgMtZ88.<  
你的模型是什么样的？雷电的电流等效成什么样的电流源？

图片:g1.jpg

c(!{_+q"  
B,ZLX/c9  
对雷电电流应该是一个瞬时电流，包含频率非常高，所以你的时间步长取多少呢？你要研究的区域范围大概是多少？想得到什么样的结果? n
)>nfnh  
Qx[ nR/  
个人感觉应该加上一些自己的东西，这样比较好一些，才容易出文章。 %ZZ
W p%uf  
比如： Jgf=yri  
加入雷电电流源时，计算电压，如果电压差超过一个阈值，就认为空气被击穿，更新为模型的常数从空气到导电介质，这样应该更能反映真实的情况，而且可以得到数值计算所不能得到的结果。 irzWk3@:  
@jH8x!5u:  
楼主认为如何？

楼主看一下前面的帖子有个说用高频电磁波探测井的硕士论文，看有没有对大地（土壤）等一些参数的设置。

非常谢谢gwzhao，我的电流就是一个瞬态电流，双指数波，但频率不是特别高。我打算计算电流通道1000m，半径100m范围的电磁场，采用二维fdtd，具体在我上传的论文中进行了解释，希望有时间看一看论文，给予指导，万分感谢。 7-hSso.'  
我的问题是 }2eP
~3  
1、当我知道土壤的电导率和介电常数，那么在利用fdtd处理时候，应当如何处理大地和空气的边界呢？只是改变电导率和介电常数就可以吗？如果是的话，从哪点开始迭代呢，是从选的大地的底部，还是大地表面向大地深处和空气高处两个方向迭代？ J 4EG  
2、以上问题应该选择柱坐标，那么在柱坐标下的吸收边界应该如何考虑？我上传一篇论文，其中135页讨论了fdtd的处理，136页的“2) At the r direction”下面的关于吸收边界问题的公式我不明白，也就是柱坐标的吸收边界问题，请高手指点。

还有就是斑竹有没有qq，我想和你请教一下。我的是23986928

文章中说电流的波长大约300m，雷电的频率这么低么，我不是很理解，不过不管怎么样，如果波长是300米的话，你要模拟1000×100应该没问题。 8Ilg[Drj*  
如果你不怕多点内存的话，你把介质系数用一个二维数组表示，这样就可以不对介质边界另外处理了。或者你根据FDTD的书，在边界的时候就是调用的系数有一点点差别，不复杂的。你找一本FDTD的书看一下，里面都有公式推导的，一时间光打字，不发图我也说不明白。 }-:s9Lt  
`)[bu  
FDTD是数值模拟，原则上是和迭代顺序没有关系的，无论从哪个点开始都一样，如果开始迭代点不一样，模拟结果也不一样，这明显不符合实际情况，FDTD也就失去了意义。 S/?KC^JP  
`4& GumG  
对二维的，迭代都是这样的。 R30{/KK  
for i = 0 to imax 4A_}:nU  
  for j = 0 to jmax !`yg bI.  
  ...//TODO h-Q3q:  
  end 4 |E`  
end

哈，刚刚发现上面一个网友发的FDTD硕士论文，里面就用的柱坐标，你可以看一下。

谢谢斑竹，能否说一下介质二维数组是什么？也就是行代表什么，列代表什么？ oW]~\vp^0  
我的理解是 z Qtg]@S  
对于介电常数而言，epsilon[M][N],行：横坐标的介电常数，列：纵坐标的介电常数，即： G 
@ib  
epsilon[j]代表某点的介电常数；对于电导率而言，sigma[j]代表某离散点的电导率。是这个意思吗？

没有仿真过雷电 W10=SM}  
雷电会使空气强烈离子话，可能是离子仿真吧，我认为不能单纯用一束电流束，而是粒子束