[求助]一维FDTD的无量纲化问题？

参数设置： !RJuH;8  
oOU?6nq  
时空步长：dx = 0.000199862, dt = 6e-013 0)?.rthk4S  
)l|/lj  
Omega0 = 9.4247778e+10 Rpi@^~aPE  
^Iz(V2  
V\ 7O)g
 
;ZSJ-r  
Y@+e)p{  
我现在想将我的程序无量纲化， 目前的进展是，我先用阻抗Z_0将EH归一化，已经不成问题。 YXdd=F  
w[A$bqz  
现在的问题是用什么合适的常量去将时间步长dt，空间步长dx无量纲化？ `h:$3a:5  
J'%  
我的做法是用点源Sin或SinGauss型中的Omega0将dt无量纲化，即Dt = dt*Omega0, )G a%Eg9  
_Kw<4$0<p  
与此同时用波矢K0将dx无量纲化，即Dx = dx * k0 = dx*2*pi / lambda = dx * Omega0 / speed_c. UZ`GS$D@  
+-VkRr#  
但是在计算时EH出现无穷，可能因为此时计算H时的系数K_H _b = Dt*K0*Z_0 / mu_0 = 1.69529e+0070是个很大的值。 2[#7YWs  
wTL&m+xr  
ZE!dg^-L  
)Ycjx~   
Wd R~  
请问到底这样做对吗？怎样做才合适呢？ =xa:>Vh#  
qNH= W?T8.  
9qHbV 9,M  
[KT'aGK$  
&6yh4-(7  
\}:&Hl+  
]4$t'wI.  
~Lyy7B9  
YCWt%a*I'  
t81}jD  
cr&sI=i  
exnFy-  
K_fQFuj+  

ncGg@$E  
[xTu29X.  

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

前几天讨论过无量纲化的问题  就是用高斯单位制  不知道dx dy为什么也要无量纲化？

我这里的无量纲就是将所有的有关物理量的单位全置换为没有量纲的。 因为我还要添加等离子体所以必须要无量纲化。 :epjJ1mW  
AXQG  
高斯单位制尚且使用不熟。要大幅改动程序实属不易，现在我就想找个量消去相关物理量的单位即可。 %+1;iuDL  
`H^?jX>7  
cP]5Qz  
Y!+q3`-%T  
t$*V*gK{  
hPM:=@N$  
`akbzHOM  
1z#0CX}Y/H  
===================================================================================

使用CFL number. S ＝ c *dt/dx. 对应的更新H的 系数为S/Z_O. 对应更新E的系数为 Z_O/S。

用高斯制看看。 >};6>)0  
公式推导好了再说，不要这边改一点，那边改一点，最后也不知道系数哪里有问题了。