Lightning space electromagnetic field has always been a hot topic in the electromagnetic field. Mature algorithms such as the Finite-Difference Time-Domain (FDTD), Transmission-Line Modeling (TLM) and frequency domain Method of Moments (MoM) have a wide range of applications in computing lightning problems. Since the Euler equation in Computational Fluid Dynamics (CFD) and the Maxwell equation in electromagnetics have the same conservation form and the Discontinuous Galerkin (DG) method has been widely tried in the flow field problem, a DG method based on CFD is proposed to discretize the time domain Maxwell equation. The grid partition parallel technology is used to accelerate the calculation, and the DG-based spherical RCS example is used for testing. The numerical results show the feasibility of the DG algorithm in solving the electromagnetic field problem. Calculation results of the electric field distribution in a near-field lightning channel are compared with those using the analytical solution and the simulation solution of a certain algorithm. The data basically agrees, showing that the method proposed is applicable for calculation of lightning electromagnetic field.
WANG Zihang
,
LYU Hongqiang
. Research on lightning electromagnetic field based on CFD high-precision algorithm[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2021
, 42(S1)
: 726366
-726366
.
DOI: 10.7527/S1000-6893.2021.26366
[1] QIAN L, ZHAO Z B, CUI X. Calculation and research of electromagnetic fields in space around lightning path with finite difference time domain method[C]//2007 International Symposium on Electromagnetic Compatibility. Piscataway: IEEE Press, 2007: 424-427.
[2] HU W, CUMMER S A. An FDTD model for low and high altitude lightning-generated EM fields[J]. IEEE Transactions on Antennas and Propagation, 2006, 54(5): 1513-1522.
[3] NODA T. A numerical simulation of transient electromagnetic fields for obtaining the step response of a transmission tower using the FDTD method[J]. IEEE Transactions on Power Delivery, 2008, 23(2): 1262-1263.
[4] 任合明, 周璧华. 地闪电磁脉冲对近地电缆外导体的耦合研究[J]. 电波科学学报, 2006, 21(5): 750-755. REN H M, ZHOU B H. Study on coupling of LEMP to cables very close to ground[J]. Chinese Journal of Radio Science, 2006, 21(5): 750-755(in Chinese).
[5] BABA Y, RAKOV V A. Electric fields at the top of tall building associated with nearby lightning return strokes[C]//2007 18th International Zurich Symposium on Electromagnetic Compatibility. Piscataway: IEEE Press, 2007: 179-182.
[6] RACHIDI F, JANISCHEWSKYJ W, HUSSEIN A M, et al. Current and electromagnetic field associated with lightning-return strokes to tall towers[J]. IEEE Transactions on Electromagnetic Compatibility, 2001, 43(3): 356-367.
[7] BERENGER J P. Long range propagation of lightning pulses using the FDTD method[J]. IEEE Transactions on Electromagnetic Compatibility, 2005, 47(4): 1008-1011.
[8] SARTO M S. Innovative absorbing-boundary conditions for the efficient FDTD analysis of lightning-interaction problems[J]. IEEE Transactions on Electromagnetic Compatibility, 2001, 43(3): 368-381.
[9] 代健,苏东林,赵小莹.基于FDTD的雷电脉冲对飞机介质舱体内干扰作用的研究[J].电子与信息学报,2009, 31(9): 2093-2098. DAI J, SU D L, ZHAO X Y. A research of lightning pulse interference with the medium cabin in airplane based on FDTD[J]. Journal of Electronics & Information Technology, 2009, 31(9): 2093-2098(in Chinese).
[10] YEE K E. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media[J]. IEEE Transactions on Antennas and Propagation, 1966, 14(3): 302-307.
[11] CANGELLARIS A C, WRIGHT D B. Analysis of the numerical error caused by the stair-stepped approximation of a conducting boundary in FDTD simulations of electromagnetic phenomena[J]. IEEE Transactions on Antennas and Propagation, 1991, 39(10): 1518-1525.
[12] COHEN G, MONK P. Mur-Nedelec finite element schemes for Maxwell’s equations[J]. Computer Methods in Applied Mechanics and Engineering, 1999, 169(3-4): 197-217.
[13] PIPERNO S, REMAKI M, FEZOUI L. A nondiffusive finite volume scheme for the three-dimensional Maxwell’s equations on unstructured meshes[J]. SIAM Journal on Numerical Analysis, 2002, 39(6): 2089-2108.
[14] 高煜堃, 陈红全. 基于非结构网格格点FVTD算法的电磁散射模拟[J]. 南京航空航天大学学报, 2013, 45(3): 415-423. GAO Y K, CHEN H Q. Electromagnetic scattering simulation based on cell-vertex unstructured-grid FVTD algorithm[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2013, 45(3): 415-423(in Chinese).
[15] 郭勇颜. 电磁散射数值模拟方法及其应用研究[D].绵阳: 中国空气动力研究与发展中心, 2007. GUO Y Y. The numerical simulation method of electromagnetic scattering and its application research[D].Mianyang: China Aerodynamics Research and Development Center, 2007(in Chinese).
[16] 陈刚. 结构/直角网格下时域有限体积法在计算电磁中的应用研究[D].西安: 西北工业大学, 2006. CHEN G. The application of finite volume method in time domain in computational electromagnetics under structure/right-angle mesh[D]. Xi'an: Northwestern Polytechnical University, 2006(in Chinese).
[17] COCKBURN B, KARNIADAKIS G E, SHU C W. Discontinuous Galerkin methods[M].Berlin, Heidelberg: Springer Berlin Heidelberg, 2000.
[18] COCKBURN B, SHU C W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws II: General framework[J]. Mathematics of Computation, 1989, 52(186): 411-435.
[19] COCKBURN B, LIN S Y, SHU C W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: One-dimensional systems[J]. Journal of Computational Physics, 1989, 84(1): 90-113.
[20] COCKBURN B, HOU S, SHU C W. The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. IV: The multidimensional case[J]. Mathematics of Computation, 1990, 54(190): 545-581.
[21] LU H Q, XU Y D, GAO Y K, et al. A CFD-based high-order Discontinuous Galerkin solver for three dimensional electromagnetic scattering problems[J]. Advances in Engineering Software, 2015, 83: 1-10.
CJA