PythonでFDTD(3)で導出した電界と磁界に関する式をプログラムし易いように書き直すと、 となる. 時間ステップは、Courantの安定条件より算出している. 実際のコード # Python 3.13.7 # NumPy 2.4.0 # Matplotlib 3.10.8 # scipy 1.16.3 # cupy-cuda13x 13.6.0 # 2026.2.16 # ©ぐらすあー ...
Maxwell's curl equations (1D, linear, non-magnetic): dE_x/dt = (1/eps) * dH_y/dz dH_y/dt = (1/mu0) * dE_x/dz Yee grid (1D): E[k] is located at z = k*dz H[k] is located at z = (k + 0.5)*dz E and H are ...
一次元FDTDの最終的なコードなんだけど、cupyを使うと何故か吸収境界が正常に動作しないようで、反射が起きてしまう.とりあえずnumpyを使ったバージョンを掲載しておく. # Python 3.13.7 # NumPy 2.4.2 # Matplotlib 3.10.8 # scipy 1.17.02 # cupy-cuda13x 14.0.0 # 2026.2.17 # ©ぐらす ...
%% code for 1D FDTD (gaussian pulse hitting a lossy dielectric medium-pml boundary condition-using electric flux density) Ex_low_m1 = 0; %leftmost boundary condition Ex_low_m2 = 0; %leftmost boundary ...
Abstract: Complex or asymmetrical geometries can be difficult to model using ordinary, 1D, circuit modelling techniques. Simulating these geometries using 2D or 3D FDTD (Finite Difference Time Domain) ...
一部の結果でアクセス不可の可能性があるため、非表示になっています。
アクセス不可の結果を表示する