[RFEDA原创]发一个电磁波在梯度等离子体中传播的程序回馈一下论坛

可以看到的。感觉问题还是很多啊。 w%iwxo  
1）如果是一维的，比如说电磁场变量为Ex和Hy，一般的推导过程是：i表示位置，n表示n时刻。 c(2?
./\|  
Ex(i,n + 1) = Ex(i,n) + f(Hy(i,n),Hy(i+1,n)); _Z9d.-  
Hy(i,n+1) = Hy(i,n) + f(Ex(i,n),Ex(i+1,n)); *>mjUT}cP  
你相当于直接写成: hi/d%lNZ  
Ex(i,n+1) = g(Ex(i,n),Ex(i,n-1),Ex(i+1,n),Ex(i+1,n-1))的差分格式。 >JsVIfAF
 
我们做模拟的时候，用电场和磁场变量原因很多，一个就是可以随时的观察电场和磁场的变化，相位差，幅度什么的。如果程序里面只有一个un的话，感觉是求解一个方程的稳定解，而不是数值模拟了，对吧。 a@|H6:|  
2）在电子运动方程中 e*E = me * d(u_e) / dt ， 这里面就没有考虑到电子的碰撞，要加上的话，应该是： '\op$t/  
me*d(u_e)/dt = e*E - me* ve * u_e;记不太清楚了，ve是电子碰撞频率，从公式中也可以看到这一项是阻止电子运动的，对应着能量损耗。 +75"Q:I  
3）这个方程默认前提是在等离子体密度均匀的前提下得到的，如果等离子体密度不均匀，这个是不成立的。比如说电子密度不均匀，那电子受力分析就要加上电子密度梯度的影响，这样方程中应该含有对电子密度的梯度了。当然如果根据分析，这个密度梯度很小，对电子运动的影响可以忽略不计，那也可以。

谢谢提点。用一个波动方程来表示确实有很多问题，因为波在等离子的反射是因为粒子碰撞的影响，而且等离子介电常数的实部和虚部影响着波的透射和反射，这些都没有体现到。迟点我再做个加入磁场变量的仿真和这个波动方程的结果做个比较。

你结果应该能体现反射的一些结果吧，但是衰减应该是没有的，你可以把频谱的能量加起来看看，是不是没变？

确实如此啊，虽然波峰的峰值衰减了，但主要是有色散引起的，传播区域内的能量基本没变化啊

是啊，这是由你的方程决定的，要是能量总和有变化的话，那就奇怪了。

我也在研究等离子体方向的课题，与各位交流下。如有相关程序，可否共享，谢谢。QQ：465109005