US 12,353,804 B1
Electromagnetic field simulation method based on face-centered cubic and subgrid technique
Ming Fang, Hefei (CN); Zhenshan Guan, Hefei (CN); Jian Feng, Hefei (CN); Zhixiang Huang, Hefei (CN); Ke Xu, Hefei (CN); Xuesong Deng, Hefei (CN); Yang Shi, Hefei (CN); Donghua Zhou, Hefei (CN); Yihao Li, Hefei (CN); and Wei Yang, Hefei (CN)
Assigned to Anhui University, Hefei (CN)
Filed by Anhui University, Hefei (CN)
Filed on Dec. 31, 2024, as Appl. No. 19/006,480.
Claims priority of application No. 202410019009.0 (CN), filed on Jan. 5, 2024.
Int. Cl. G06F 30/20 (2020.01)
CPC G06F 30/20 (2020.01) 7 Claims
OG exemplary drawing
 
1. An efficient three-dimensional electromagnetic field simulation method for simulating fine grides based on a face-centered cubic (FCC) and a subgrid, comprising:
based on the faced-centered cubic and the subgrid, setting periodic boundary conditions and a metal plate, and constructing an electromagnetic field simulation model;
setting a source point in a subgrid region through the electromagnetic field simulation model, selecting a detecting point, and simulating an electromagnetic field;
wherein setting the periodic boundary conditions comprises:
obtaining an updating equation of the face-centered cubic based on sampling point distribution of the face-centered cubic grid combined with the sampling point distribution of finite-difference time-domain (FDTD) and the updating equation;
obtaining a structure and a constant of a simulation region;
generating the periodic boundary conditions according to the structure and the constant of the simulation region and the updating equation of the faced-centered cubic;
for a transverse magnetic (TM) mode:
Ezn+1(i,j)=CA(m)·Ezn+1(i,j)−

OG Complex Work Unit Math
wherein CA(m) is a third coefficient, CB(m) is a fourth coefficient, Ez is an electric field value in a z direction, Hy is a value of a magnetic field in a y direction, Hx is a value of a magnetic field in an x direction, n is a space step, n−1 represents a space step at a previous moment, n+1 represents a space step at a next moment, n−½ represents a space step of a first half moment, n+½ represents a space step of a second half moment, (m) is a dielectric constant at a position m, σ(m) is a conductivity at the position m, and Δt is a moment step;
Hx and Hy in the TM mode are:

OG Complex Work Unit Math

OG Complex Work Unit Math
wherein CQ(m) is a second coefficient, i is a coordinate in an x direction, j is a coordinate in a y direction, Δx is a grid size in the x direction, Δy is a grid size in the y direction, and CP(m) is a first coefficient
wherein simulating the electromagnetic field comprises:
S1, placing the subgrid region based on the simulation region;
S2, initializing an electric field component and a magnetic field component of coarse and fine grids,
S3, obtaining an electric field component of the coarse grid at a next moment based on the FCC and an initialized the electric field component of the coarse grid;
S4, transferring the electric field component of the coarse grid at the next moment to the fine grid by a linear interpolation, and obtaining an electric field component and a magnetic field component of the fine grid at the next moment based on a size of the coarse grid of the FCC,
S5, weighting the magnetic field component of the fine grid at the next moment to obtain a magnetic field component of the coarse grid at a next moment of a junction to reduce moment consumption and memory occupation of FCC, and
S6, based on the FCC obtaining the magnetic field component of the coarse grid at the next moment, and ending if an iteration moment is greater than a preset iteration moment; otherwise, returning to the S3.