1/2/2023 0 Comments Fdtd algorithm![]() ![]() In addition, how to conduct direction transformation of the fields on the z-axis is also a research difficult in the 3D cylindrical coordinate system. How to reduce this kind of error in the iterative FDTD calculation has become a difficult problem. In fact, due to the existence of the term in the 3D cylindrical coordinate system, the decomposition of the fields along the direction will lead the splitting error. If we adopt the conventional FDTD method to discretize the cylindrical structure with the Cartesian grid, a significant staircasing error appears. Moreover, the geometry of interest may consist of fine structures. However, in many applications, we have to deal with 3D cylindrical structures such as in optical fiber communication, integrated optics, and defense industry. Obviously, the good feasibility of the efficient Laguerre scheme had been conformed in the rectangular coordinate system. To further reduce the splitting error, the modified perturbation term was introduced in meanwhile, the Gauss–Seidel iterative procedure was applied, which made the new Laguerre-based FDTD scheme more efficient and accurate. This provided excellent computational accuracy and can be efficiently parallel-processed on a computing cluster. To overcome this problem, a factorization-splitting scheme was used to resolve the huge sparse matrix into six tridiagonal matrixes, and then the chasing method was used to solve them. The main drawback of the conventional Laguerre-based FDTD algorithm is that it requires solving a very large sparse matrix. Thus, Laguerre-based FDTD formulation has the advantage of less numerical error when a larger time step is used. By transforming the time-domain problem to the Laguerre domain using the temporal Galerkin’s testing procedure, the transient solution is independent of time discretization. In recent years, the unconditionally stable Laguerre-based finite-difference time-domain (FDTD) algorithm has been applied to simulate transient electromagnetic problems in the Cartesian coordinate. Meanwhile, the proposed algorithm is extremely useful for the problems with fine structures in the 3D cylindrical coordinate system. The computation results show that the proposed algorithm can be better than the ADI-FDTD algorithm in terms of efficiency and accuracy. To verify the performance of the proposed algorithm, two scattering numerical examples are given. Such a treatment scheme can reduce the splitting error to a low level and obtain a good convergence in other words, it can improve the efficiency and accuracy than other algorithms. Different from the previously developed iterative procedure used in the efficient FDTD algorithm, a new perturbation term combined with the Gauss–Seidel iterative procedure is introduced to form the new Laguerre-based FDTD algorithm in the 3D cylindrical coordinate system. Here an efficient Laguerre-based finite-difference time-domain iterative algorithm is proposed. ![]()
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |