Modelado y estudio del líder descendente negarivo en una descarga eléctrica atmosférica y su influencia en la localización del punto de impacto del rayo
Authors:   , Hermoso Costa Juan Ramón 
Published: July 2023

In lightning protection techniques a basic aspect is the identification of the impact point possible, which is determined by the evolution of the initial phase of the lightning, known as the stepped leader, the ground geometry and structures on it. In general, the available theoretical models are consistent with ideal situations and do not include situations such as non-vertical trajectories, non-flat terrain or complex structures.

This work reviews and compares the most outstanding models in the literature on this subject and proposes a model for the study of the stepped leader based on the Finite Element Method (FEM). With the intention of bringing the process to the real situation of progress of a downward negative leader, who can follow a trajectory 3D, several modeling aspects are contemplated, as: adjusting the charge density in the channel, the length and direction of jump, which can be modified depending on the height that is the leader, and the propagation delays caused by the distance from the observation point. On the other hand, the model allows the implementation of 3D geometries and the incorporation of complex structures and non-flat terrain, providing results in situations when it is not possible to have theoretical approaches.

Simulations have been performed with different proposals for charge distribution in the channel, in situations of leader vertical, inclined and 3D trajectory. The model has been evaluated, comparing their results with theoretical ideal cases and reproducing some real situations from field measurements. The application of the model in simulating real cases for downward negative leader provides parameter values that are not available through the recorded measurements and are useful for studying this type of discharge. In this work we have focused on some of them, and we have obtained interesting results as: estimating the speed of the leader according of height (from the measures of ¿E (t)), the definition of new relationships between charge density in the channel, ¿Emax (leader close to ground), the return peak current and the observation point distance, as well as we can verify the conditions of inception for a discharge into a storm cloud.

Finally, we use the model to evaluate the striking distance, a basic concept on protection techniques. It is inferred a new equation for its calculus and the results are compared with proposals available in the literature, showing good agreement with the recommendation of IEC 62305-1.

The results encourage to continue the study and they open the possibility to extend the method to other stages, not only on the downward leader, but also on other phenomena related to the formation and evolution of storm situations.
En las técnicas de protección contra el rayo un aspecto básico es la identificación del posible punto de impacto, que viene determinado por la evolución de la fase inicial de la descarga atmosférica, conocida como líder escalonado (stepped leader), y por la geometría del terreno y estructuras sobre él. En general, los modelos teóricos disponibles se ajustan a situaciones ideales y no contemplan situaciones como trayectorias no verticales, terreno no plano o estructuras complejas.

En este trabajo se revisan y comparan los modelos más destacados en la bibliografía sobre el tema y se propone un modelo para el estudio del stepped leader basado en el Método de los Elementos Finitos (FEM). Con la intención de aproximar el proceso a la situación real de avance de un líder negativo descendente, que puede seguir una trayectoria 3D, en el modelado se contemplan aspectos como: el ajuste de la densidad de carga en el canal, la longitud y dirección del salto, que pueden ser modificadas en función de la altura a la que se encuentre el líder, y los retardos de propagación provocados por la distancia del punto de observación. Por otra parte, el modelo permite la implementación de geometrías en 3D y la incorporación de estructuras complejas y terreno no plano, facilitando obtener resultados en situaciones en las que no es posible disponer de aproximaciones teóricas.

Se han realizado simulaciones con diferentes propuestas de distribución de carga en el canal y para situaciones de líder vertical, inclinado y con trayectoria 3D. El modelo ha sido evaluado, contrastando sus resultados con casos teóricos ideales y reproduciendo situaciones reales a partir de diversas medidas de campo. La aplicación del modelo en la simulación de casos reales para líder negativo descendente nos proporciona valores de parámetros que no están disponibles a través de las medidas registradas y que son de utilidad para el estudio de las descargas negativas descendentes. En este trabajo nos hemos centrado en algunos de ellos y se han obtenido resultados de interés como: una estimación de la velocidad de avance del líder en función de la altura (a partir de medidas disponibles de ΔE(t)), la definición de nuevas relaciones entre densidad de carga en el canal, ΔEmax (líder próximo a tierra), la corriente de pico de retorno y la distancia del punto de observación, además de verificar las condiciones de inicio de descargas en una nube de tormenta.

Por último, se utiliza el modelo para la evaluación de la distancia de cebado (striking distance), concepto básico en técnicas de protección. Se infiere una nueva ecuación para su cálculo y se comparan los resultados con las propuestas disponibles en la literatura, mostrando una buena concordancia con la recomendación de la norma IEC 62305-1.

Los resultados obtenidos animan a continuar con el estudio y abren la posibilidad de ampliar la aplicación del método a otras etapas, no sólo en la descarga descendente, sino también de otros fenómenos relacionados con la formación y evolución de las situaciones de tormenta.
En les tècniques de protecció contra el llamp, un aspecte bàsic és la identificació del possible punt d’impacte, que ve determinat per l’evolució de la fase inicial de la descàrrega atmosfèrica, coneguda com a líder esglaonat (stepped leader), i per la geometria del terreny i les estructures que hi ha a sobre d’ell. En general, els models teòrics dsiponibles s’ajusten a situacions ideals i no contemplen situacions com trajectòries no verticals, terreny no pla o estructures complexes.

En aquest treball es revisen i comparen els models més destacats en la bibliografia sobre el tema i es proposa un model per a l’estudi del stepped leader basat en el Mètode dels Elements Finits (FEM). Amb la intenció d’apropar el procés a la situació real d’avançament d’un líder negatiu descendent, el qual pot seguir una trajetòria 3D, en el modelat es contemplemn aspectes com: l’ajustament de la densitat de càrrega al canal, la longitud i la direcció del salt, que poden ser modificades en funció a l’altura on es trobi el líder i als retards de propagació provocats per la distància del punt d’observació. D’una altra banda, el model permet la implementació de geometries en 3D i la incorporació d’estructures complexes i terreny no pla, facilitant així l’obtenció de resultats en situacions on no es possible disposar d’aproximacions teòriques.

S’han realitzat simulacions amb diferents propostes de distribució de càrrega al canal i per a situacions de líder vertical, inclinat i amb trajectòria 3D. El model ha estat avaluat, contrastant els seus resultats amb casos teòrics ideals i reproduint situacions reals a partir de diferents mesures de camp. L’aplicació del model en la simulació de casos reals per a líder negatiu descendent ens proporciona valors de paràmetres que no estan disponibles a travès de les mesures registrades i que són d’utilitat per a l’estudi de les descàrregues negatives descendents. En aquest treball ens n’hem centrat en alguns i s’han obtingut resultats d’interès com: una estimació de la velocitat d’avançament del líder en funció d’altura (només a partir de mesures disponibles de ΔE(t)), la definició de les noves relacions entre densitat de càrrega al canal, ΔEmax (líder proper a terra), el corrent del pic de retorn i la distància del punt d’observació, a més de verificar les condicions d’inici de descàrregues en un núvol de tempesta.

Per últim, s’utilitza el model per a l’avaluació de la distància de d'encebament (striking distance), concepte bàsic en tècniques de protecció. S’infereix una nova equació per al seu càlcul i es comparen els resultats amb les propostes disponibles en la literatura, mostrant una bona concordància amb la recomanació de la norma IEC 62305-1.
Els resultats obtinguts animen a continuar amb l’estudi i obren la posibilitat d’ampliar l’aplicació del mètode a d’altres etapes, no només en la descàrrega descendent, sinó també d’altres fenòmens relacionats amb la formació i evolució de les situacions de tempesta.