Li Dan, Zhang Yijun, Lü Weitao, et al. A 3D self-consistent propagation model of the lightning leader. J Appl Meteor Sci, 2015, 26(2): 203-210. DOI:  10.11898/1001-7313.20150208.
Citation: Li Dan, Zhang Yijun, Lü Weitao, et al. A 3D self-consistent propagation model of the lightning leader. J Appl Meteor Sci, 2015, 26(2): 203-210. DOI:  10.11898/1001-7313.20150208.

A 3D Self-consistent Propagation Model of the Lightning Leader

DOI: 10.11898/1001-7313.20150208
  • Received Date: 2014-04-03
  • Rev Recd Date: 2014-12-02
  • Publish Date: 2015-03-31
  • The research on simulation of lightning discharge and propagation is one of key issues nowadays in the knowledge of lightning physics. To study propagation characteristics of the lightning leader and the interaction between the lightning leader and structures, a 3D self-consistent propagation model of lightning leader is proposed considering the limitation of the 2D random model that cannot describe physical features of lightning discharge three-dimension well. This model proposed can properly simulate the downward stepped leader in 3D space based on the initial conditions, such as the background electric field and the geometry of ground installations. Also, the upward continuous progression of positive connecting leaders from its inception to the final jump as the downward negative stepped leader approaches the ground can be simulated. Meanwhile, this model can self-consistently calculate physical parameters of the leader while it develops in three dimensional space based on the transient electric field under the influence of the downward stepped leader and the geometry of ground structures, such as tip location every step, leader velocity, current density, charge per unit length and the length every step moved. Compared with the optical and electrical data of the natural and artificially triggered lightning detected in Guangzhou experiment base, as well as research results available in literature, simulation results of the self-consistent propagation model show that the initiation of the positive upward leader mainly depends on the electric field intensity and the geometry of the structure in the simulated domain. As the upward leader initiates and propagates towards the downward stepped leader, its velocity increases from 104 m·s-1 to 105 m·s-1 order of magnitude gradually along with time during the first 0-600 μs, and then its velocity obviously increases to 106 m·s-1 order of magnitude when the upward positive leader connects to the downward stepped leader, during which the average velocity reaches about 5×105 m·s-1. Furthermore, the current intensity in the channel also increases with the upward positive leader moving forward, and the trend has a good coherence with the variation of the relative brightness of the leader, which is consistent with the optical data and existing research. Through simulation results of lightning flashes striking on high structure, it is also derived that the charge per unit length in the propagation channel of upward positive leader initiated from high structure reaches about 50.0-108.0 μC·m-1, with an average value of 64.3 μC·m-1. In addition, the average length of upward positive leader is approximately 417 m. The research of 3D self-consistent propagation model of lightning leader can provide reference for further study of characteristics of lightning flashes striking on structures.
  • Fig. 1  A 3D self-consistent propagation model of lightning leader

    Fig. 2  The spatial discretization model

    (a) gridding of the simulation field, (b)7-point finite difference scheme

    Fig. 3  The numerical simulation flowchart of calculating the upward leader inception condition

    Fig. 4  The geometry algorithm of CSM

    Fig. 5  Simulation cases of lightning flashes striking on tall structures

    (H is the height of the structure, R is the horizontal distance between the upward leader and the center of the roof, α is the linear increation coefficient of the electric field intensity)

    Fig. 6  Variations of the upward leader speed under different conditions

    Fig. 7  Variations of currents of the upward leader under different conditions

    Table  1  Lengths and charge per-unit length of the upward leader under different conditions

    序号 先导发展总长度/m 平均线电荷密度/(μC·m-1)
    个例1 412 50.0
    个例2 431 64.8
    个例3 373 59.6
    个例4 421 67.8
    个例5 407 35.6
    个例6 458 108.0
    DownLoad: Download CSV
  • [1]
    Carrara G, Thione L.Switching surge strength of large air gaps:A physical approach.IEEE Transactions on Power Apparatus and Systems, 1976, 95(2):512-524. doi:  10.1109/T-PAS.1976.32131
    [2]
    Rizk F.A model for switching impulse leader inception and breakdown of long air-gaps.IEEE Transactions on Power Delivery, 1989, 4(1):596-603. doi:  10.1109/61.19251
    [3]
    Petrov N I, Waters R T.Determination of the striking distance of lightning to earthed structures.Proceedings of the Royal Society A:Mathematical, Physical and Engineering Sciences, 1995, 450(1940):589-601. doi:  10.1098/rspa.1995.0102
    [4]
    Akyuz M, Cooray V.The Franklin lightning conductor:conditions necessary for the initiation of a connecting leader.Journal of Electrostatics, 2001, 51-52:319-325. doi:  10.1016/S0304-3886(01)00113-9
    [5]
    Eriksson A J.The Lightning Ground Flash-An Engineering Study.Faculty of Engineering, University of Natal, Pretoria, 1979.
    [6]
    Dellera L, Garbagnati E.Lightning strike simulation by means of the Leader Progression Model, Part I:Description of the model and evaluation of free-standing structures.IEEE Transactions on Power Delivery, 1990, PWRD-5:2009-2023. doi:  10.1088/0022-3727/39/16/028/meta
    [7]
    Dellera L, Garbagnati E.Lightning strike simulation by means of the Leader Progression Model:Ⅱ.Exposure and shielding failure evaluation of overhead lines with assessment of application graphs.IEEE Transactions on Power Delivery, 1990, PWRD-5:2023-2029. https://www.researchgate.net/publication/3271538_Lightning_stroke_simulation_by_means_of_the_leader_progression_model_II_Exposure_and_shielding_failure_evaluation_of_overhead_lines_with_assessment_of_application_graphs
    [8]
    Rizk F.Modeling of transmission lines:Exposure to direct lightning strokes.IEEE Transactions on Power Delivery, 1990, PWRD-5:1983-1989. https://www.researchgate.net/publication/3271533_Modeling_of_Transmission_Line_Exposure_to_Direct_Lightning_Strokes
    [9]
    Becerra M, Cooray V.Dynamic Modeling of the Lightning Upward Connecting Leader Inception.28th ICLP, 2006:543-548. http://industry.wanfangdata.com.cn/dl/Detail/NSTLQK?id=NSTLQK_NSTL_QKJJ0213906559
    [10]
    Becerra M, Cooray V.A simplified physical model to determine the lightning upward connecting leader inception.IEEE Transactions on Power Delivery, 2006, 21(2):897-908. doi:  10.1109/TPWRD.2005.859290
    [11]
    Becerra M, Cooray V.On the interaction of lightning upward connecting positive leaders with humans.IEEE Transactions on Electromagnetic Compatibility, 2009, 51(4):1001-1008. doi:  10.1109/TEMC.2009.2033265
    [12]
    Becerra M, Cooray V.A self-consistent upward leader propagation model.Journal of Physics D:Applied Physics, 2006, 39:3708-3715. doi:  10.1088/0022-3727/39/16/028
    [13]
    谢施君, 何俊佳, 陈维江, 等.避雷针迎面先导发展物理过程仿真研究.中国电机工程学报, 2012, 32(10):32-40. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGDC201210005.htm
    [14]
    任晓毓, 张义军, 吕伟涛, 等.闪电先导随机模式的建立与应用.应用气象学报, 2011, 22(2):194-202. doi:  10.11898/1001-7313.20110208
    [15]
    任晓毓, 张义军, 吕伟涛, 等.雷击建筑物的先导连接过程模拟.应用气象学报, 2010, 21(4):450-457. doi:  10.11898/1001-7313.20100408
    [16]
    Singer H, Steinbigler H, Weiss P.A charge simulation method for the calculation of high voltage fields.IEEE Transactions on Power Apparatus and Systems, 1974, PAS-93:1660-1668. doi:  10.1109/TPAS.1974.293898
    [17]
    Lu W, Chen L, Ma Y, et al.Lightning attachment process involving connection of the downward negative leader to the lateral surface of the upward connecting leader.Geophys Res Lett, 2013, 40:1-5. doi:  10.1029/2012GL054022
    [18]
    Lu W, Zhang Y, Chen L, et al.Attachment Processes of Two Natural Downward Lightning Flashes Striking on High Structures.30th International Conference on Lightning Protection, 2010. https://www.researchgate.net/publication/224142653_Observation_and_preliminary_analysis_on_the_attachment_process_of_lightning_flashes_striking_on_high_structures
    [19]
    Lu W, Chen L, Zhang Y, et al.Characteristics of unconnected upward leaders initiated from tall structures observed in Guangzhou.J Geophys Res, 2012, 117(D19):doi: 10.1029/2012JD018035.
    [20]
    McEachron K B.Lightning to the Empire State building.Journal of The Franklin Institution, 1939, 227:149-217. doi:  10.1016/S0016-0032(39)90397-2
    [21]
    Berger K.Novel observations of lightning discharges:Results of research on Mount San Salvatore.Journal of The Franklin Institution, 1967, 283:478-525. doi:  10.1016/0016-0032(67)90598-4
    [22]
    Kito Y, Horii K, Higashiyama Y, et al.Optical aspects of winter lightning discharges triggered by the rocket-wire technique in Hukuriku district of Japan.J Geophys Res, 1985, 90:6147-6157. doi:  10.1029/JD090iD04p06147
    [23]
    Wada A, Asakawa A, Shindo T, et al.Leader and Return Stroke Speed of Upward-initiated Lightning.Proceedings of International Conference on Atmospheric Electricity (ICAE), 2003, C3-20.
    [24]
    张义军, 杨少杰, 吕伟涛, 等.2006—2011年广州人工触发闪电观测试验和应用.应用气象学报, 2012, 23(5):513-522. doi:  10.11898/1001-7313.20120501
    [25]
    李俊, 张义军, 吕伟涛, 等.一次多回击自然闪电的高速摄像观测.应用气象学报, 2008, 19(4):401-411. doi:  10.11898/1001-7313.20080403
    [26]
    李俊, 吕伟涛, 张义军, 等.一次多分叉多接地的空中触发闪电过程.应用气象学报, 2010, 21(1):95-100. doi:  10.11898/1001-7313.20100113
    [27]
    张义军, 周秀骥.雷电研究的回顾和进展.应用气象学报, 2006, 17(6):829-834. doi:  10.11898/1001-7313.20060619
    [28]
    Lalande P.Study of the Lightning Stroke Conditions on a Grounded Structure.Paris:University of Paris, 1996.
    [29]
    郑栋, 张义军, 吕伟涛, 等.先导-回击模型与人工触发闪电特征参数计算.中国电机工程学报, 2006, 26(23):151-157. doi:  10.3321/j.issn:0258-8013.2006.23.027
  • 加载中
  • -->

Catalog

    Figures(7)  / Tables(1)

    Article views (3103) PDF downloads(1240) Cited by()
    • Received : 2014-04-03
    • Accepted : 2014-12-02
    • Published : 2015-03-31

    /

    DownLoad:  Full-Size Img  PowerPoint