[求助]求FDTD三维柱坐标半径等于0处奇异点的解决方案

大家好： w%xCTeK[  
       最近遇到棘手的问题，希望各位高手不吝赐教！在三维柱坐标下，当半径r=0时是奇异点，此处的网格跟其它的不同，不知道该如何处理？ !wh&>3~  
      在r=0附近的网格很小，会影响时间步长，不知道有没有 .. 'fY9a(Xt.  
G/{ ~_&t
 

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

额，难度贼大

我班门弄斧一下吧： w"zE_9I\  
_D;@v?n6!O  
（1） 首先要把场分量在网格中的具体位置标出， 一种放置是磁场全部被安置在网格表面的中心点，磁场分量方向垂直于其表面。电场分量被放置在网格线上。柱坐标中心处周围就是尖顶网格。电场分量Ez，E_theta （ 或者叫E_phi）， 及磁场Hr无限趋进在Z坐标轴上。 *@S@x{{s  
（2） 由此放置， 你会发现你是不需要计算中心处的E_theta 和Hr， （ 窍决就在构建的FDTD方程中） eSy(~Y
 
（3） 但是需要计算轴线上的Ez， 这个比较简单， 应用积分形式的Maxwell's eqaution中安培定律，时变电场Ez在圆环上的面积分等于环绕的磁场H_theta的环线积分。 W>-B [5O&[  
Ciihsm  

好像是这么回事，一说就有些印象了。