[求助]遇到问题  求助

  本人新手，刚接触等离子体模拟的知识，看了不少文献，大部分都是用FDTD的方法在解电磁场，，，搜索了下论坛 大家的FDTD程序也是如火纯青，，，，我最近看了个外文文献，上面用的PIC的方法模拟激光等离子体相互作用，，，Au+D,作用出质子（无碰撞），，，没用FDTD的方法，我也照着他的方法用Fortran编写了程序，最近调试中发现不少问题，，，我作了如下尝试， A[K:/tB  
    1.不加任何粒子，（模拟激光在真空的传播），空间中的电磁场分布正确，断定场的解法正确， _.u~)Q`6  
    2.只加入2e5的Au离子，出来的离子速度在0.05左右（激光最大场强2e13），基本认定 程序粒子处理部分正确 \?aOExG I  
    3 只加入1.9e6的电子，出现问题，，，，在最后用GUASS消去法解线性方程时，算不了多久就会出现发散的信号（Fail),尝试了很多时间步长，都不行，，，关于稳定性的问题也有很多分析的文章，我也调试了多次，可结果都是如此，，，现在还在继续调试中，也不知道是什么样的问题，，，，实在是不知道再该怎么测试了，，，， ?SBh^/zf  
  有懂的同仁们 教下在下，，，不甚感激    ！！！！ .. rD_Ss.\^g  
=#7s+d-  

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