最新时域有限差分方法林志立PPT课件.ppt

《最新时域有限差分方法林志立PPT课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、时域有限差分方法林志立电磁学基本方程麦克斯韦方程组(MaxwellEquations)JamesClerkMaxwell(1831–1879)(法拉第感应定律)(安培环路定律)(高斯定律－电场)(高斯定律－磁场)物质本构关系(ConstitutiveRelations)（欧姆定律）（电极化）（磁化）一切电磁场和电磁波问题均可由以上方程，以及各类具体的边界条件所决定！麦克斯韦方程组中的运算符散度（Divergence)旋度（Curl)连续函数的偏微分运算FDTD空间偏微分的近似以Hz为例：类似地，可实现各电磁场分量的空间偏微分计算。FDTD时间偏微分的近似以Hz为例：t=(

2、n+1/2)Δtt=nΔtt=(n-1/2)Δt时间上的推移(1)电场在时间上取整数倍的Δt;t=n*Δt;(2)磁场在时间上取(整数＋1/2)倍的Δt;t=(n+1/2)*Δt;麦克斯韦方程的离散化近似以Hz为例：上式即为Hz的更新方程，由前一时刻的磁场和前半时刻的临近空间格点的电场即可求出最新时刻的磁场。麦克斯韦方程的离散化近似采取类似的步骤，可以推导出其它场量的更新表达式：例如，对于Ez：FDTD的离散参数的选择★元胞尺寸：边长小于最短波长的1/10,以减小数值色散。数值色散方程：理想色散方程：要求：例如，取例如，取FDTD的离散参数的稳定性条件★时间步长：Cour

3、ant稳定性条件对于非色散介质，时间步长不能大于以下表达式：（vonNeumannmethod）为了保持稳定性，该方程的所有解的模必须小于1。Ｚ域数值色散方程：介质电磁参量的设定不同的元胞的电磁参量应设置为所在空间所代表的介质的介电常数和磁导率。长方体，物质I长方体，物质II球体，物质III空气场量与介质参数要对应色散介质的FDTD模拟以Lorentz介质为例：更新方程：近似，求系数Z-变换时域：频域：色散介质的FDTD模拟模拟Lorentz色散介质的不同方法：MSEapproach:FDTD编程流程主循环初始化输出结果编程举例1：一维FDTD问题x基本旋度方程：X向电导

4、率X向磁电导率X向电流X向磁流编程举例1：一维FDTD问题Matlab程序代码：%Defineinitialconstantseps_0=8.854187817e-12;%permittivityoffreespacemu_0=4*pi*1e-7;%permeabilityoffreespacec=1/sqrt(mu_0*eps_0);%speedoflight%Defineproblemgeometryandparametersdomain_size=1;%1Dproblemspacelengthinmetersdx=1e-3;%cellsizeinmetersdt=3

5、e-12;%durationoftimestepinsecondsnumber_of_time_steps=2000;%numberofiterationsnx=round(domain_size/dx);%numberofcellsin1Dproblemspacesource_position=0.5;%positionofthecurrentsourceJz1.定义物理常量2.定义问题的参量和结构尺寸%InitializefieldandmaterialarraysCeze=zeros(nx+1,1);Cezhy=zeros(nx+1,1);Cezj=zeros(nx

6、+1,1);Ez=zeros(nx+1,1);Jz=zeros(nx+1,1);eps_r_z=ones(nx+1,1);%freespacesigma_e_z=zeros(nx+1,1);%freespaceChyh=zeros(nx,1);Chyez=zeros(nx,1);Chym=zeros(nx,1);Hy=zeros(nx,1);My=zeros(nx,1);mu_r_y=ones(nx,1);%freespacesigma_m_y=zeros(nx,1);%freespace编程举例1：一维FDTD问题（续）3.初始化场量和介质参量阵列电场电流部分磁场磁流部

7、分%CalculateFDTDupdatingcoefficientsCeze=(2*eps_r_z*eps_0-dt*sigma_e_z)..../(2*eps_r_z*eps_0+dt*sigma_e_z);Cezhy=(2*dt/dx)..../(2*eps_r_z*eps_0+dt*sigma_e_z);Cezj=(-2*dt)..../(2*eps_r_z*eps_0+dt*sigma_e_z);Chyh=(2*mu_r_y*mu_0-dt*sigma_m_y)..../(2*mu_r_y*mu_0+dt*sigma