# 9月13日 于成都黄龙溪 1 #!/usr/bin/perl # Author : Leon Email: yangli0534@gmail.com # fdtd simulation , plotting with gnuplot, writting in perl # perl and gnuplot software packages should be installed before running this program #use Time::HiRes qw(sleep);…

我用perl和python写了相同功能的一段程序,计算一维fdtd,用gnuplot动态显示,可是python的数据没有显示出来,看横纵坐标的变化数据是正确收到了的,如最后的图片,求大神指点,谢谢. #!/usr/bin/perl # Author : Leon Email: yangli0534@gmail.com # fdtd simulation , plotting with gnuplot, writting in perl # perl and gnuplot software pa…

源文件 #!/usr/bin/env python from os import popen class gnuplot_leon: # Author : Leon Email: yangli0534@gmail.com # a gnuplot api of python def __init__(self): self.gnuplot = popen('gnuplot','w') self.write = self.gnuplot.write self.flush = self.gnuplot…

源代码 #!/usr/bin/env python from math import exp from gnuplot_leon import * imp0 = 377.0 class fdtd_leon: # Author : Leon Email: yangli0534@gmail.com # fdtd simulation #initialization def __init__(self,size=400,time=0,MaxTime=1000,delay = 30, width = 1…

背景: 基于公式1.42(Ez分量).1.43(Hy分量)的1D FDTD实现. 计算电场和磁场分量,该分量由z方向的电流片Jz产生,Jz位于两个理想导体极板中间,两个极板平行且向y和z方向无限延伸.           平行极板相距1m,差分网格Δx=1mm. 电流面密度导致分界面(电流薄层)磁场分量的不连续,在两侧产生Hy的波,每个强度为5×10^(-4)A/m.因为电磁波在自由空间中传播,本征阻抗为η0.           显示了从左右极板反射前后的传播过程.本例中是PEC边界,正切电场…

post post Table of Contents 1. Post-processing 1.1. Reverse flow 1.1.1. reasons 1.1.2. solutions 1.2. variable definitions– ch34, field function definitions, user guide, fluent 1.3. Residuals 1.3.1. access to the residual values for each cell in my f…

Table of Contents 1. contacts 2. Paper digest 2.1. LES vs. RANS 2.2. Dynamics of Transient Fuel Injection, Mixing, and Auto-Ignition 2.3. rothamer research group, Experimental Studies of Transient Jets 2.4. the transient start of supersonic jets 2.5.…

Octave是一个旨在提供与Matlab语法兼容的开放源代码科学计算及数值分析的工具,是Matlab商业软件的一个强有力的竞争产品. 参考:[ML:Octave Installation] General Installation files for all platforms are available at the GNU Octave Repository on SourceForge. The Gnu Octave Wiki has installation instructions f…

待续 %1D FDTD simulation with a simple absorbing boundary condition % and a TFSF boundary between hy[] and ez[]. SIZE = ; ez = zeros(,SIZE); hy = zeros(,SIZE); imp0=377.0; maxTime = ; : : maxTime hy(SIZE)=hy(SIZE-);% : : SIZE - hy(mm) = hy(mm) + (ez(mm…

链接 There are several resources. But it depends on what you actually want to learn...Let me explain: First of all, There are many Numerical Modeling Techniques in Electromagnetic, and people are always confused about which one to choose. But if you ar…