💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。
⛳️座右铭：行百里者，半于九十。
📋📋📋本文目录如下：🎁🎁🎁
目录
⛳️赠与读者
👨💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑。哲学是科学之母，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。
或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎
💥1 概述
"二维FDTD（有限差分时域）解决完全电导体边界条件问题"
该程序的目标是解决包含沿xy平面极化磁场组分Hy和Hx以及沿z方向极化电场Ez的TM波的Maxwell方程。在一个空间内，所有自由空间的物理参数都不被标准化为1，而是给定实数和已知值。更新是使用标准的更新方程进行的，这些方程是从Maxwell的旋度方程的差分形式获得的，其中具有非常低的电导率和磁导率（4x10^(-4)单位）的电场和磁场被纳入其中。场点由Yee算法描述的网格中定义。H场被定义在每一半个坐标步长的空间点。更确切地说，Hx部分被定义在每一半个y坐标和完整x坐标处，Hy部分被定义在每一半个x坐标和完整y坐标处，E场即Ez部分在每个完整的x和y坐标点被定义。此外，在这里，空间步长被取为1微米，而不是在之前程序中假设的无量纲域中的1单位。同样，时间步进是使用Leapfrog时间步进进行的。在此，H场（即Hx和Hy）在每一半个时间步被更新，E场（即Ez）在每一个完整时间步被更新。这通过两种交替的向量更新进行，跨越空间网格的部分，在那个特定时间点，波从源点开始到达，避免在不必要的情况下更新网格中的所有点。这些空间更新在主要的时间更新循环内，跨越整个时间网格。此外，在循环开始之前，这里用作更新方程的乘法因子的矩阵被初始化，以避免在每次循环迭代中重复计算相同的内容，这是对优化的细微尝试。这里的边界条件是完全电导体边界，即无论外部场的影响如何，边界网格点的电场值均为零。
在空间范围中心定义了一个电场源，这是一个硬源，因此其值不会由于外部场的干扰而发生变化，也就是说，该源是一个完全导电体。可以使用正弦、高斯和脉冲变量改变源的形式。源有四种标准形式 - 单位时间步、脉冲、高斯和正弦形式。在每个时间步骤上显示Ez场的颜色缩放图覆盖整。
二维有限差分时域（Finite-Difference Time-Domain, FDTD）方法是现代电磁场数值模拟中的一个重要工具，尤其适用于解决高频电磁波在复杂结构中的传播与散射问题。在考虑完全电导体（Perfect Electric Conductor, PEC）边界条件时，这种模拟技术的应用尤为关键，因为PEC边界条件在射频、微波、以及天线设计等领域极为常见，它要求电场在边界上必须垂直于边界，以模拟实际中金属等完美导体表面的无损反射特性。
PEC边界条件在FDTD中的实现
在二维FDTD模拟中，实现PEC边界通常涉及以下几点：
-
更新电场边界规则：在PEC边界上，电场分量必须满足垂直于边界的条件。这意味着在每次时间步进时，位于边界上的电场分量要么被强制设为零（对于法向分量），要么其更新规则需要特殊处理，以保证电场沿边界法向的连续性。
-
使用辅助场：一种常用的技术是引入所谓的“辅助场”，在边界附近增加一层额外的网格点，通过这些辅助场的更新来间接满足PEC边界条件，同时确保计算区域内的场值不受影响。
-
镜像法：另一种经典方法是镜像法，即在边界外侧虚拟地放置一组电荷或电流源，它们的配置与内部实际源关于边界镜像对称，由此产生的场正好满足PEC边界的要求。
-
精度与稳定性：实现PEC边界时，需要特别关注算法的稳定性与精度，确保在模拟过程中不会因边界条件的处理不当引入不必要的误差或导致计算不稳定。
研究挑战与前沿
-
提高计算效率：随着计算域的增大或时间分辨率的提高，PEC边界的处理成为限制FDTD计算效率的一个重要因素。研发更高效的算法，减少计算量，同时保持高精度，是一个持续的挑战。
-
多物理场耦合：在某些应用场景中，电磁场与热传导、流体力学等其他物理过程耦合，如何在考虑这些耦合作用的同时，准确实现PEC边界，是一个复杂且前沿的研究方向。
-
复杂几何形状的处理：真实世界中的电导体结构往往具有复杂的几何形状，如何在FDTD框架内高效且精确地处理非规则PEC边界，是一个亟待解决的问题。
-
多尺度问题：在纳米尺度到微米尺度的器件设计中，精确模拟电磁场与PEC表面的相互作用，要求FDTD方法能够有效处理尺度不连续性带来的挑战。
总之，二维FDTD方法在解决完全电导体边界条件问题上的研究，不仅是电磁学和向量场理论的基础应用，也是推动现代科技发展中多个关键技术进步的关键因素。随着算法的不断优化和计算能力的增强，这一领域的研究正持续向着更高精度、更广应用领域迈进。
📚2 运行结果
部分代码：
% Initialization of permittivity and permeability matrices
epsilon=epsilon0*ones(xdim,ydim);
mu=mu0*ones(xdim,ydim);
% Initializing electric and magnetic conductivity matrices
sigma=4e-4*ones(xdim,ydim);
sigma_star=4e-4*ones(xdim,ydim);
%Choice of nature of source
gaussian=0;
sine=0;
% The user can give a frequency of his choice for sinusoidal (if sine=1 above) waves in Hz
frequency=1.5e+13;
impulse=0;
%Choose any one as 1 and rest as 0. Default (when all are 0): Unit time step
%Multiplication factor matrices for H matrix update to avoid being calculated many times
%in the time update loop so as to increase computation speed
A=((mu-0.5*deltat*sigma_star)./(mu+0.5*deltat*sigma_star));
B=(deltat/delta)./(mu+0.5*deltat*sigma_star);
%Multiplication factor matrices for E matrix update to avoid being calculated many times
%in the time update loop so as to increase computation speed
C=((epsilon-0.5*deltat*sigma)./(epsilon+0.5*deltat*sigma));
D=(deltat/delta)./(epsilon+0.5*deltat*sigma);
% Update loop begins
for n=1:1:time_tot
%if source is impulse or unit-time step
if gaussian==0 && sine==0 && n==1
Ez(xsource,ysource)=1;
end
% Setting time dependent boundaries to update only relevant parts of the
% vector where the wave has reached to avoid unnecessary updates.
🎉3 参考文献
文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。
[1]毛云龙.基于表面阻抗的时域有限差分方法边界条件研究[D].哈尔滨工程大学[2024-07-03].DOI:CNKI:CDMD:1.1018.291752.
[2]李响,杨明武,王伟吉,等.时域有限差分电磁建模系统的研究与设计[J].Journal of Terahertz Science and Electronic Information Technology, 3(4)[2024-07-03].
[3]贾洪祥.时域有限差分法电磁场计算采用吸收边界条件的一种二维角点处理[C]//第五届全国青年计算物理学术交流会.0[2024-07-03].DOI:ConferenceArticle/5aa0c3a9c095d722207faa1b.
🌈4 Matlab代码实现
资料获取，更多粉丝福利，MATLAB|Simulink|Python资源获取
762
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
