A method and system for fast calculation of discharge photon emission spectrum and total light intensity in a gas medium

Abstract

The application relates to a fast calculation method and system of discharge photon radiation spectrum and total light intensity in a gas medium, which comprises the following steps: step 1, establishing a discharge process numerical model of a to-be-detected discharge object, and solving discharge parameter distribution at each spatial position in a discharge area; step 2, calculating an electron temperature field at each spatial position based on the discharge parameter distribution; step 3, calculating a local radiation spectrum at each spatial position based on the electron temperature and combining spectral parameters of light-emitting particles at the spatial position, and obtaining an overall radiation spectrum and total light intensity; and step 4, repeating steps 2-3 to obtain full-time-space evolution results of the overall radiation spectrum and the total light intensity in the discharge process, and completing the calculation process. Compared with the prior art, the application has the advantages of efficient calculation of full-time-space spectrum evolution in a complex discharge process.

Description

Technical Field

[0001] This invention relates to the field of gas discharge diagnosis and radiation characteristics, and in particular to a method and system for rapidly calculating the discharge photon radiation spectrum and total light intensity in a gaseous medium. Background Technology

[0002] In the fields of gas discharge diagnosis and radiation characteristic analysis, accurately calculating the gas radiation spectrum generated during the discharge process is of great significance for the study of discharge mechanisms, analysis of active particles, interpretation of experimental results, and optical diagnosis. Gas discharge is accompanied by complex physical processes such as electron collision excitation, dissociation, ionization, recombination, and radiative transitions. Different particles at different energy level transitions will form emission spectra with characteristic information. Obtaining the local spectra and overall radiation results at various times and spatial locations during the discharge process will help to more comprehensively reveal the spatiotemporal evolution laws of complex discharge processes.

[0003] In existing technologies, one type of method focuses on the numerical simulation of the discharge process. This typically involves solving the Poisson equation, the particle continuity equation, and related auxiliary equations to obtain the spatiotemporal distribution of parameters such as electric field, particle concentration, and electron temperature within the discharge region. This type of method can describe the evolutionary characteristics of the discharge process itself relatively well, but it usually cannot directly output the corresponding local radiation spectrum results. Another type of method focuses on the calculation of gas emission spectra. This typically involves calculating the intensity and broadening / superimposing the emission spectral lines of the target particles under given temperature, particle population, or other preset parameters to obtain theoretical spectral results. While this type of method can achieve spectral solutions, its input parameters are usually derived from given assumptions or simplified conditions, lacking direct coupling with the spatiotemporal evolution results of the actual discharge process. Therefore, it is difficult to reflect the spectral evolution characteristics across the entire spatiotemporal range during complex discharge processes.

[0004] Furthermore, when it is necessary to further calculate the local spectra at various spatial locations within the entire discharge field at each time point during a complex discharge process, and obtain the overall spectrum or total light intensity based on these local spectra, existing technologies typically face the following problems: First, discharge simulation and spectral calculation are separated, lacking a unified and direct coupling bridge (e.g., patent application CN104360253A); second, traditional spectral calculations are often based on single-point parameters or overall assumptions, making them difficult to apply to point-by-point local spectral reconstruction across the entire field; third, when dealing with the full spatiotemporal spectral evolution of complex discharge processes, the calculation process is cumbersome and computationally intensive, hindering efficient analysis and engineering applications. Therefore, it is necessary to propose a new method for calculating discharge radiation spectra to achieve a more efficient and complete characterization of the radiation characteristics of complex discharge processes. Summary of the Invention

[0005] The purpose of this invention is to provide a rapid calculation method and system for the emission spectrum and total intensity of discharge photons in a gaseous medium that can truly reflect the spatiotemporal evolution under complex discharge conditions.

[0006] The objective of this invention can be achieved through the following technical solutions: A rapid method for calculating the emission spectrum and total intensity of discharge photons in a gaseous medium includes the following steps: Step 1: Establish a numerical model of the discharge process of the object under test, and solve for the distribution of discharge parameters at each spatial location within the discharge region; Step 2: Calculate the electron temperature at each spatial location based on the discharge parameter distribution; Step 3: Based on the electron temperature, calculate the local radiation spectrum at each spatial location by combining the spectroscopic parameters of the luminescent particles at the spatial location, and obtain the overall radiation spectrum and total light intensity. Step 4: Repeat steps 2-3 to obtain the spatiotemporal evolution results of the overall radiation spectrum and total light intensity during the discharge process, thus completing the calculation process.

[0007] Furthermore, the discharge object to be tested is a gas discharge structure under different insulation defect conditions, or a discharge system formed under different electrode forms, different insulating media, different gas pressure conditions, and different voltage excitation methods.

[0008] Furthermore, the numerical model of the discharge process includes the particle continuity equation and the Poisson equation, wherein the particle continuity equation is expressed as: , , , In the formula, the subscript , , They are electrons, cations, and anions, respectively. These represent the number densities of electrons, cations, and anions, respectively. These represent the mobility of electrons, cations, and anions, respectively. These are the diffusion coefficients for electrons, cations, and anions, respectively. For the electric field in space, For gradient operators, For photoionization source terms, The ionization coefficient, The adhesion coefficient, The recombination coefficient of electrons and cations. The recombination coefficient of anions and cations; The Poisson equation is expressed as: , , In the formula, This is the space electromotive force. For elementary charge, is the vacuum permittivity.

[0009] Furthermore, the calculation expression for the photoionization source term is: , In the formula, Coordinates are Photoionization source term at the location, The second ionization coefficient of photoionization. for and The solid angle between them , They are respectively in Electron density and electron migration rate at that location Let be the excitation coefficient of the excited state. For the first The absorption coefficient of each excited state.

[0010] Furthermore, the step of obtaining the electron temperature includes: The reduced electric field is extracted from the discharge parameter distribution; Based on the pre-defined correspondence between the reduced electric field and the average electron energy, the average electron energy of the reduced electric field at each spatial location is obtained; Based on the average electron energy, the relationship between the average electron energy and the electron temperature is established using the Maxwell distribution, expressed as: , In the formula, The average electron energy, Boltzmann's constant, For electron temperature; Based on the relationship between the average electron energy and the electron temperature, the electron temperature at each spatial location was calculated. ,in For any given moment, For any spatial location.

[0011] Furthermore, the step of calculating the local radiation spectrum at each spatial location includes: Read any time point and any spatial location electron temperature at the location and obtain spatial location Information on the particle concentration of the luminescent particles; According to the electron temperature Based on particle concentration information and combined with the energy level parameters of the luminescent particles, the population distribution of the corresponding excited-state particles is calculated. The population distribution of the excited-state particles is described using particle number density. For the ... There are 1 excited energy level, and its particle number density is expressed as: , In the formula, For the first Particle number density at each excitation energy level The total number density of luminescent particles. For the first The statistical weights of each excitation energy level To be at electron temperature The partition function under, For the first The energy of an excited energy level Boltzmann's constant; Based on the population distribution of the excited-state particles and the transition parameters between the corresponding excited energy levels, the local emission intensity of each target spectral line at each spatial location is calculated and expressed as: , in: , In the formula, To be caused by the excitation energy level To excitation energy level The local emission intensity of the target spectral line produced by the transition. To be caused by the excitation energy level To excitation energy level The spontaneous emission transition probability corresponding to the transition. Let be Planck's constant. To be caused by the excitation energy level To excitation energy level The radiation frequency corresponding to the transition. The speed of light in a vacuum. To be caused by the excitation energy level To excitation energy level The center wavelength corresponding to the transition; For each target spectral line at each spatial location, the spectral lines are superimposed based on the corresponding spectral line center position and the local emission intensity to obtain the local radiation spectrum within the target wavelength range at each spatial location, as shown below: , In the formula, For the first A spatial location is at wavelength The local spectral intensity at that location, i.e., the local radiation spectrum. For the first At the spatial location of the first The local emission intensity that triggers energy level transitions. This is the line shape function corresponding to the target spectral line. To superimpose all target spectral lines.

[0012] Furthermore, the steps of obtaining the overall radiation spectrum and total light intensity include: For any time Suppose the discharge region is discretized as Each spatial calculation point is used as a spatial location; The first Each spatial computation point at time... The local radiation spectrum below is denoted as The local radiation spectra corresponding to all spatial calculation points at the same time are summed to obtain the time. The overall radiation spectrum of the entire discharge region is expressed as: , In the formula, For a moment The entire discharge region at wavelength The overall radiation spectrum at that location; Based on the overall radiation spectrum, the target band Integrating the light intensity within the target wavelength band yields the total light intensity, expressed as: , In the formula, For a moment Total light intensity within the target wavelength band These are the starting and ending wavelengths of the target band; Based on the overall radiation spectrum, across the entire wavelength range Integrating, we obtain the total light intensity across the entire wavelength range, expressed as: , In the formula, For a moment Total light intensity across the entire wavelength range. , These represent the lower and upper limits for the entire wavelength range.

[0013] Furthermore, the step of obtaining the spatiotemporal evolution results of the overall radiation spectrum and total light intensity during the discharge process includes: The total computation time of the discharge process is discretized into... Each time step is denoted as _____. ; For each moment Perform steps 2-3 to obtain the overall radiation spectrum and total light intensity; Based on the overall radiation spectrum at each moment, a sequence of the temporal variation of the overall radiation spectrum is formed. ,in, For a moment The entire discharge region at wavelength The overall radiation spectrum at that location; Based on the total light intensity at each moment, a sequence of total light intensity changing over time is generated, including the sequence of total light intensity changing over time within the target wavelength band. and the sequence of total light intensity over time across the entire wavelength range ,in, For a moment Total light intensity within the target wavelength band For a moment The total light intensity across the entire wavelength range.

[0014] Furthermore, the electron temperature is replaced by one of the following: reduced electric field, electric field strength, electron density, ion concentration, and particle concentration.

[0015] The present invention also provides a fast calculation system for the emission spectrum and total intensity of discharge photons in a gaseous medium, including a memory, a processor, and a program stored in the memory. When the processor executes the program, it implements the fast calculation method for the emission spectrum and total intensity of discharge photons in a gaseous medium as described above.

[0016] Compared with the prior art, the present invention has the following beneficial effects: (1) This invention obtains the distribution of discharge parameters at various spatial locations within the discharge area by establishing a numerical model of the discharge process and extracts the electron temperature. The electron temperature is used as a coupling bridge between the numerical simulation of the discharge process and the radiation spectrum calculation, so that the radiation spectrum calculation no longer depends on a single given assumed parameter, but can directly and realistically reflect the spatiotemporal evolution characteristics of the discharge process. By combining the simulation of the discharge process with the calculation of the local radiation spectrum, the calculation of the local radiation spectrum, the overall radiation spectrum and the total light intensity in the entire spatiotemporal range during the complex discharge process can be realized. This solves the technical problem that traditional spectral calculation methods usually rely on given assumed parameters and are difficult to reflect the spatiotemporal evolution of the real discharge process.

[0017] (2) In view of the problem that traditional spectral calculation methods have large computational load and low efficiency in full-field point-by-point spectral reconstruction under complex discharge conditions, this invention relies on the numerical model of the discharge process to directly output the discharge parameter distribution, and to construct the field and solve the spectrum. Without the need for additional complex experimental calibration, it can realize the efficient calculation of the full-time and space-time spectral evolution in complex discharge processes, and improve the completeness and engineering applicability of discharge radiation characteristic analysis.

[0018] (3) The local radiation spectrum calculation process of this invention uses electron temperature as the core input parameter and directly maps the local energy state during the discharge process to the corresponding local luminescence characteristics. Compared with the traditional method of calculating a single theoretical spectrum only under the overall assumption, this invention can realize the time-by-time and space-by-space location local spectrum solution in complex discharge processes, thereby more accurately characterizing the spatiotemporal evolution law of radiation characteristics during the discharge process.

[0019] (4) In this invention, for any time step and any spatial location, the corresponding local radiation spectrum can be calculated based on the electron temperature and particle concentration information at that location, combined with the energy level parameters and transition parameters of the target luminescent particles. This allows the invention to move beyond theoretical spectral solutions under single-time, single-point, or overall average conditions, and instead construct a local spectral field that corresponds one-to-one with the spatiotemporal distribution of the discharge process. Based on this technical approach, this invention can achieve point-by-point calculation of local radiation spectra at each time and spatial location during complex discharge processes, thereby solving the problem of difficulty in describing the full spatiotemporal spectral evolution in existing technologies. By characterizing the spatial non-uniformity and temporal dynamics of radiation characteristics during complex discharge processes, the completeness of radiation analysis results is improved.

[0020] (5) By uniformly calculating and integrating the local radiation spectra at various spatial locations, this invention can further obtain the overall radiation spectrum and total light intensity results, solving the problem in the prior art that local parameters are difficult to directly convert into overall observation results. Furthermore, it is applicable to the analysis of radiation characteristics under different insulation defects, different dielectric environments, different discharge structures, and different optical detection conditions. Therefore, this invention can not only describe local luminescence behavior but also output overall radiation results that better meet the actual optical detection requirements, thus forming a complete mapping relationship between internal discharge parameters and external observables.

[0021] (6) This invention repeatedly performs the processes of constructing the electron temperature field, calculating the local spectrum, and solving the overall radiation results at multiple time steps, thereby forming an evolution sequence of the overall spectrum and the evolution sequence of the total light intensity over time. This allows for further analysis of the radiation enhancement, attenuation, peak occurrence time, and changes in different spectral line components during the discharge development process. Compared to technical solutions that can only obtain static spectral results at a single moment, this invention provides continuous dynamic information for analyzing the discharge development process by obtaining the time evolution results of the radiation spectrum and light intensity during the discharge process, thus revealing the dynamic luminescence mechanism in complex discharge processes more comprehensively. Attached Figure Description

[0022] Figure 1 This is a schematic diagram of the method flow of the present invention. Detailed Implementation

[0023] The present invention will now be described in detail with reference to the accompanying drawings and specific embodiments. These embodiments are based on the technical solution of the present invention and provide detailed implementation methods and specific operating procedures. However, the scope of protection of the present invention is not limited to the following embodiments.

[0024] Example 1 This embodiment provides a rapid calculation method for the photon radiation spectrum and total light intensity of discharge in a gaseous medium. This addresses the problem in existing technologies where numerical simulation of the discharge process and calculation of gas emission spectra are separated, lacking a direct coupling mechanism, and making it difficult to characterize the spectral evolution characteristics across the entire spatiotemporal range during complex discharge processes. This method uses the electron temperature during the discharge process as a bridge to couple the discharge simulation results with local spectral calculations, enabling the calculation of local radiation spectra at various times and spatial locations within the discharge field, and further obtaining the overall spectrum or total light intensity results. Specifically, as... Figure 1 As shown, the method includes the following steps: S1. Establish a numerical model of the discharge process and determine the distribution of discharge parameters at each time and spatial location within the discharge region; use the particle continuity equation and Poisson equation to simulate the entire discharge process.

[0025] In this embodiment, the discharge object under test can be a gas discharge structure under different insulation defect conditions, or a discharge system formed under different electrode forms, different insulating media, different gas pressure conditions, and different voltage excitation methods. First, a numerical model of the discharge process of the discharge object under test is used, and the relevant control equations of the discharge process (including the Poisson equation, the particle continuity equation, and auxiliary equations describing the nonlocal discharge process) are solved to obtain the distribution of discharge parameters at each time and spatial location within the discharge region. The discharge parameters include at least electron temperature, and preferably also include reduced electric field, electric field strength, electron density, ion concentration, particle concentration, or other parameters reflecting the local discharge state. Specifically, this step includes: Particle continuity equations: The entire model can be viewed as a fluid model, in which there are processes of particle generation and loss. These processes are usually described by continuity equations for various types of particles (electrons, cations, and anions) and Poisson's equation for the electric field. The continuity equations for electrons, positive ions, and negative ions are written as follows: (1) (2) (3) in, For the electric field in space, For photoionization source terms, subscript , and These represent electrons, positive ions, and negative ions, respectively. Among the coefficients containing subscripts... Number density, For mobility, The ionization coefficient, The adhesion coefficient, Represents the recombination coefficient of electrons and positive ions. This represents the recombination coefficient between negative and positive ions. The gradient operator is given by equations (1) to (3). Equations (1) to (3) describe the movement of electrons, positive ions and negative ions in space under the consideration of various reactions such as photoionization, migration, diffusion and recombination, attachment and ionization.

[0026] Regarding the photoionization source term, Dhali et al.'s research indicates that there is very little research and information on photoionization and light emission, with only Morrow et al. providing a calculation formula for photoionization in one-dimensional discharge: (4) in, Coordinates are Photoionization source term at the location, The secondary ionization coefficient represents the photoionization. for and The solid angle between them , They represent in Electron density and electron migration rate at that location Indicates the first The absorption coefficient of each excited state, Let be the excitation coefficient for this excited state. Upon verification, it was found that in the formula... and It is difficult to confirm these results, and this formula can only serve as an ideal model and cannot be used for numerical calculations. Based on previous research, under the premise of simulating photoionization with a fixed initial value, background ionization can be used to replace the photoionization process. That is, a suitable constant is added to the continuity equation, assuming that a fixed photoionization source always exists, and a specified number of electrons and positive ions are generated by photoionization every second. This constant is taken as... .

[0027] Poisson's equation: Simultaneously, the Poisson equation is introduced to describe the relationship between electric potential (electric field in space) and space charge density: (5) (6) in, It is the vacuum permittivity. It is space electromotive force. It is the elementary charge.

[0028] The above particle continuity equation and Poisson equation constitute a numerical model of the discharge process. By substituting equation (6) into equation (5) and solving equations (1) to (5), the process of particle density change with time at each point in space can be analyzed.

[0029] S2. Extract electron temperature from the discharge parameter distribution and construct an electron temperature field corresponding to the spatiotemporal evolution of the discharge process.

[0030] After completing the numerical simulation of the discharge process, it is necessary to further construct the electron temperature field corresponding to the spatiotemporal distribution of the discharge region, as an input parameter for subsequent local radiation spectrum calculations. The electron temperature can be directly derived from the solution results of the discharge model, or it can be calculated by reducing the electric field, average electron energy, or other intermediate physical quantities. Specifically, this step includes: In this step, the reduced electric field distribution at each time point and spatial location is first extracted from the discharge simulation results (discharge parameter distribution). The reduced electric field is usually expressed as the ratio of electric field strength to gas number density, and is used to characterize the strength of the electric field effect on electrons under unit particle number density conditions. Since the reduced electric field can reflect the energy acquisition ability of electrons in a partial discharge environment, it can be used as an important intermediate parameter to characterize the energy state of electrons.

[0031] Furthermore, based on the pre-defined correspondence between the reduced electric field and the average electron energy, the average electron energy at each time and spatial location is calculated. The reduced electric field-average electron energy relationship can be obtained through calculations based on electron transport characteristics, or it can be derived from a pre-established lookup table, numerical fitting relationship, or the solution results of the Boltzmann equation. In a preferred embodiment, a one-to-one correspondence between the reduced electric field and the average electron energy is pre-established according to the type of discharge medium and operating conditions; after obtaining the reduced electric field at a certain spatial location, the corresponding average electron energy is obtained through table lookup, interpolation, or function mapping.

[0032] After obtaining the average electron energy, the relationship between the average electron energy and the electron temperature is established based on Maxwell's distribution. For an electronic system that satisfies the Maxwell's distribution assumptions, the average electron energy and the electron temperature satisfy the following relationship: (7) in: Indicates the average energy of an electron; Represents the Boltzmann constant; Indicates electron temperature.

[0033] As shown in the above equation, given the average electron energy at each spatial location, the corresponding electron temperature can be calculated. Therefore, for any given moment... and any spatial location The average electron energy can be determined first from the reduced electric field value at that location, and then the electron temperature can be calculated to obtain the corresponding electron temperature value. .

[0034] By repeating the above calculation process for all spatial locations and all time steps within the discharge region, the electron temperature field corresponding to the spatiotemporal evolution of the discharge process can be constructed. The electron temperature field is essentially a parametric field that varies with time and space, and can be expressed as: (8) For a two-dimensional discharge model, the electron temperature field is expressed as a function of two-dimensional spatial coordinates and time; for a three-dimensional discharge model, the electron temperature field can be further expressed as spatial coordinates. and time The function.

[0035] Through the above processing, an electron temperature field data set corresponding to the spatiotemporal distribution of the discharge process can be established.

[0036] In this embodiment of the invention, the electron temperature field serves as a crucial bridge connecting the simulation of the discharge process with local spectral calculations. On one hand, the electron temperature field originates from the spatiotemporal simulation results of the discharge process, reflecting the non-uniform distribution of electron energy states at different locations and times during the discharge. Compared to methods that rely solely on given uniform parameters for spectral calculations, it more realistically characterizes the discharge radiation properties. On the other hand, the electron temperature can be directly used to solve for subsequent luminescent particle population distribution and local radiation spectra. Therefore, by constructing the electron temperature field through the above steps, a direct mapping from the discharge process parameter field to the local spectral field can be achieved, providing a foundation for calculating the full spatiotemporal evolution of radiation spectra during complex discharge processes.

[0037] S3. Using the electron temperature at each spatial location as the input parameter for local spectral calculation, and combining the energy level parameters, transition parameters, and particle concentration information of the target luminescent particles, calculate the local radiation spectrum at each spatial location.

[0038] After constructing the electron temperature field corresponding to each time and spatial location within the discharge region, it is further necessary to calculate the corresponding local radiation spectrum based on the electron temperature at each spatial location, thereby establishing a mapping relationship between the discharge parameter field and the local spectral field. In this embodiment, the luminescent particles can be atoms, ions, molecules, or combinations thereof; the spectral calculation object can be a single characteristic spectral line (target spectral line), multiple characteristic spectral lines, a spectral band within a specific wavelength range, or the radiation spectrum within the entire target wavelength range. Specifically: In this step, for any moment within the discharge region and any spatial location First, read the electron temperature at that location. Simultaneously, the particle concentration information of the luminescent particles (as target luminescent particles) at this spatial location is obtained. The particle concentration can be directly derived from the discharge simulation results, or it can be calculated from relevant parameters in the discharge simulation results. In a preferred embodiment, the particle concentration is preferably characterized using positive ion concentration, as a basic parameter for subsequent excited-state particle population calculations.

[0039] Furthermore, based on the electron temperature and particle concentration, and combined with the energy level parameters of the target luminescent particles, the population distribution of the corresponding excited-state particles is calculated. In a preferred embodiment, the population relationship of excited-state particles is described according to the Boltzmann distribution, that is, for the ... The particle number density of an excited energy level can be expressed as: (9) in: Indicates the first Particle number density at each excitation energy level; This represents the total number density of the target luminescent particles; Indicates the first Statistical weights of each excitation energy level; Indicates the electron temperature The partition function below; Indicates the first The energy of an excited energy level; Represents the Boltzmann constant; This indicates the electron temperature at that spatial location.

[0040] Furthermore, the Boltzmann distribution described above can be replaced by the collision radiation approximation, the empirical population model, or other methods that can be used to describe the distribution of excited-state particles, in order to describe the population distribution of excited-state particles.

[0041] It can be seen that the population distribution of excited-state particles changes at different spatial locations and times due to different electron temperatures, thus characterizing the spatiotemporal nonuniformity of local luminescence capability during complex discharge processes.

[0042] After obtaining the population distribution of excited-state particles, the local emission intensity of each target spectral line is calculated by further combining the transition parameters between the corresponding energy levels. For the emission from the upper energy level... Downward energy level The transition of , and the local emission intensity of the corresponding target spectral line can be expressed as: (10) in: Indicates energy level Towards energy level The local emission intensity of the target spectral line produced by the transition; This represents the spontaneous emission transition probability corresponding to the transition; Denotes Planck's constant; This indicates the radiation frequency corresponding to the transition.

[0043] If expressed in wavelength form, the radiation frequency and wavelength satisfy the following relationship: (11) in: This represents the speed of light in a vacuum. This indicates the center wavelength corresponding to the transition.

[0044] Therefore, the local emission intensity expression for the target spectral line mentioned above can also be written as: For all the transition spectral lines of the target luminescent particle to be calculated, repeat the above process to obtain the discrete radiation intensity (local emission intensity) of each target spectral line at that spatial location: (12) In a preferred embodiment, the target luminescent particle can be an atom, ion, molecule, or a combination thereof; the transition parameters include at least the spectral center wavelength, energy level parameters, statistical weights, and transition probabilities. These parameters can be derived from standard spectral databases, experimental measurement results, or preset parameter files. In one specific embodiment, the luminescent particle is preferably the dominant luminescent ion in the target discharge medium, and its transition parameters can be derived from the NIST atomic spectral database.

[0045] To construct a local emission spectrum, for each target spectral line at any spatial location, the spectral lines can be superimposed based on their corresponding center positions and local emission intensities. This yields the local spectral distribution, or local emission spectrum, of that spatial location within the target wavelength range, expressed as: (13) in: Indicates the first A spatial location is at wavelength The local radiation spectrum at that location; Indicates the first At the spatial location of the first The local emission intensity of the transition; This represents the line shape function of the corresponding target spectral line; This indicates that all target spectral lines are superimposed.

[0046] Through the above steps, the corresponding local radiation spectrum can be obtained for each spatial location and each time step within the discharge region. This forms a local spectral field corresponding to the electron temperature field, providing a foundation for subsequent integral solutions of the overall radiation spectrum and total light intensity.

[0047] In this embodiment of the invention, the local radiation spectrum calculation process uses electron temperature as the core input parameter, directly mapping the local energy state during the discharge process to the corresponding local luminescence characteristics. Compared with the traditional method of calculating a single theoretical spectrum only under overall assumptions, this embodiment of the invention can achieve time-wise and spatial-positional local spectral solutions during complex discharge processes, thereby more accurately characterizing the spatiotemporal evolution of radiation characteristics during the discharge process.

[0048] In some embodiments, in addition to electron temperature, the local spectral calculation results can be further corrected by combining electron density, electric field strength, ion concentration, particle collision parameters or experimental correction coefficients; however, regardless of the specific correction method used, as long as the local radiation spectrum at each spatial location is calculated based on the spatiotemporal distribution parameters of the discharge process, it falls within the protection scope of this invention.

[0049] S4. Accumulate or integrate the local radiation spectra at each spatial location to obtain the overall radiation spectrum and total light intensity at the corresponding time.

[0050] After obtaining the local radiation spectra corresponding to each spatial location within the discharge region in step S3, the local radiation spectra at all spatial locations at the same time are further summarized to obtain the overall radiation spectrum of the entire discharge region at that time. Based on this, the total light intensity or the light intensity result within a specified wavelength band is calculated. Specifically, this step includes: For any time step Suppose the discharge region is discretized as The nth spatial calculation point, then the nth The local radiation spectrum corresponding to each spatial calculation point (each spatial calculation point is considered as a spatial location) at this moment is denoted as . By directly summing the local radiation spectra of all spatial points at the same moment, the overall radiation spectrum of the entire discharge region at that moment can be obtained, and its expression is: (14) in: Indicates time step The entire discharge region at wavelength The overall radiation spectrum intensity at that location; Indicates time step Next A spatial point at wavelength Local radiation spectral intensity at the location; This represents the total number of spatial points involved in the calculation within the discharge area.

[0051] Furthermore, the above point-by-point direct summation can be replaced by weighted integration based on grid area or volume, integration based on observation path, or integration based on specified area, in order to calculate the overall radiation spectrum.

[0052] After obtaining the overall radiation spectrum, different forms of overall radiation results can be calculated according to the analysis needs.

[0053] When it is necessary to obtain the total light intensity within a specified wavelength range, the overall radiation spectrum can be integrated over the target wavelength band. If the target wavelength range is... Then the total light intensity within the target wavelength band is: (15) in: Indicates time step The total light intensity within a specified wavelength range; and These represent the starting and ending wavelengths of the target band, respectively.

[0054] When the total luminous intensity needs to be obtained over the entire target spectral range, the overall radiation spectrum can be integrated over the entire calculation wavelength range. Let the target wavelength range be... The total light intensity can then be expressed as: (16) in: Indicates time step The total light intensity across the entire target wavelength range; and These represent the lower and upper limits of the target wavelength range, respectively.

[0055] Through the above calculations, this step can output the overall radiation spectrum of the entire discharge region at the corresponding time, the total light intensity of the specified band, the integral intensity of the characteristic spectral lines, and the overall total light intensity, thereby realizing the conversion from local radiation spectrum to overall observation results.

[0056] S5. Repeat steps S2 to S4 for multiple time steps to obtain the spatiotemporal evolution results of the radiation spectrum and total light intensity during the discharge process.

[0057] In this step, the local radiation spectra at various spatial locations within the discharge region at the same time are accumulated or integrated to obtain the overall radiation spectrum of the entire discharge region at the corresponding time. For specific application scenarios, the overall radiation spectrum can also be integrated at specified wavelength positions or within specified wavelength ranges to obtain characteristic spectral line intensities, band light intensities, or total light intensities. Repeating the above calculations for different time steps can further obtain the evolution results of the overall radiation spectrum and total light intensity over time during the discharge process, realizing the characterization of the spatiotemporal spectral evolution during complex discharge processes. Specifically: After obtaining the overall radiation spectrum and total light intensity results at each time step in step S4, the calculation results at different time steps are further organized and analyzed in a time sequence to construct the time evolution process of radiation spectrum and light intensity during the discharge process.

[0058] In this step, the total calculation time of the discharge process is discretized as... Each time step is denoted as: (17) For each time step ( Following steps S2 to S4, the electron temperature field, local radiation spectrum, overall radiation spectrum, and total light intensity at each given moment are obtained sequentially. This allows for the formation of a sequence showing the overall radiation spectrum changing over time. (18) in: Indicates time step The entire discharge region at wavelength The overall radiation spectrum intensity at that location; This represents the time evolution sequence of the overall radiation spectrum throughout the entire discharge process.

[0059] Furthermore, a sequence of total light intensity changes over time can be formed: (19) in: Indicates time step The total light intensity across the entire target wavelength range; This represents the time evolution sequence of the total light intensity.

[0060] When it is necessary to analyze the evolution of a specific spectral line or a particular band, its intensity sequence over time can be constructed. For example, the integral intensity of a characteristic spectral line can be expressed as: (20) in, For time step Lower characteristic spectral line intensity.

[0061] The total light intensity within a specified wavelength band can be expressed as: (twenty one) Therefore, the variation patterns of characteristic spectral line intensity, specific band light intensity, and overall total light intensity during the entire discharge process can be obtained.

[0062] After the above time series is formed, it can be further visualized or dynamically analyzed. For example, the overall radiation spectrum at different time steps can be plotted as multiple time-series spectral curves to characterize the enhancement, attenuation, or migration of different wavelength components during the discharge development process; the total light intensity, characteristic spectral line intensity, or band light intensity can also be plotted as time-varying curves to characterize the onset, development, peak, and attenuation characteristics of the overall luminescence intensity during the discharge process.

[0063] This step allows for the organization of the radiation results calculated independently at each moment into a complete temporal evolution process, thereby enabling dynamic characterization of radiation characteristics during complex discharge processes. Compared to methods that only obtain spectral results at a single moment, this step reveals the continuous law of radiation spectrum changes over time during the discharge process, making it more suitable for analyzing the discharge development process, changes in the luminescence stage, and discharge characteristics under different insulation defect conditions.

[0064] In this embodiment of the invention, step S5 extends the radiation results from a single moment to dynamic radiation results throughout the entire process. This enables the entire technical approach to not only output instantaneous local and global spectra, but also to further obtain full-time-domain spectral evolution information during complex discharge processes, thereby providing a complete data foundation for subsequent optical detection, discharge diagnosis, and state identification.

[0065] Example 2 This embodiment provides a fast calculation system for the emission spectrum and total intensity of discharge photons in a gaseous medium, including a memory, a processor, and a program stored in the memory. When the processor executes the program, it implements the fast calculation method for the emission spectrum and total intensity of discharge photons in a gaseous medium as described in Embodiment 1 above.

[0066] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0067] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. The solutions in the embodiments of the present invention can be implemented using various computer languages, such as the object-oriented programming language Java and the interpreted scripting language JavaScript.

[0068] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0069] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0070] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0071] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the invention.

[0072] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.

Claims

1. A rapid method for calculating the discharge photon radiation spectrum and total light intensity in a gaseous medium, characterized in that, Includes the following steps: Step 1: Establish a numerical model of the discharge process of the object under test, and solve for the distribution of discharge parameters at each spatial location within the discharge region; Step 2: Based on the discharge parameter distribution, calculate the electron temperature at each spatial location. The steps for obtaining the electron temperature include: The reduced electric field is extracted from the discharge parameter distribution; Based on the pre-defined correspondence between the reduced electric field and the average electron energy, the average electron energy of the reduced electric field at each spatial location is obtained; Based on the average electron energy, the relationship between the average electron energy and the electron temperature is established using the Maxwell distribution, expressed as: ， In the formula, The average electron energy, Boltzmann's constant, For electron temperature; Based on the relationship between the average electron energy and the electron temperature, the electron temperature at each spatial location was calculated. ,in For any given moment, For any spatial location; Step 3: Based on the electron temperature and the spectroscopic parameters of the luminescent particles at each spatial location, calculate the local radiation spectrum at each spatial location, and obtain the overall radiation spectrum and total light intensity. The step of calculating the local radiation spectrum at each spatial location includes: Read any time point and any spatial location electron temperature at the location and obtain spatial location Information on the particle concentration of the luminescent particles; According to the electron temperature Based on particle concentration information and combined with the energy level parameters of the luminescent particles, the population distribution of the corresponding excited-state particles is calculated. The population distribution of the excited-state particles is described using particle number density. For the ... There are several excited energy levels, and their particle number density is expressed as: ， In the formula, For the first Particle number density at each excitation energy level The total number density of luminescent particles. For the first The statistical weights of each excitation energy level To be at electron temperature The partition function under, For the first The energy of an excited energy level Boltzmann's constant; Based on the population distribution of the excited-state particles and the transition parameters between the corresponding excited energy levels, the local emission intensity of each target spectral line at each spatial location is calculated and expressed as: ， in: ， In the formula, To be caused by the excitation energy level To excitation energy level The local emission intensity of the target spectral line produced by the transition. To be caused by the excitation energy level To excitation energy level The spontaneous emission transition probability corresponding to the transition. is Planck's constant. To be caused by the excitation energy level To excitation energy level The radiation frequency corresponding to the transition. The speed of light in a vacuum. To be caused by the excitation energy level To excitation energy level The center wavelength corresponding to the transition; For each target spectral line at each spatial location, the spectral lines are superimposed based on the corresponding spectral line center position and the local emission intensity to obtain the local radiation spectrum within the target wavelength range at each spatial location, as shown below: ， In the formula, For the first A spatial location is at wavelength The local spectral intensity at that location, i.e., the local radiation spectrum. For the first At the first spatial location The local emission intensity that triggers energy level transitions. This is the line shape function corresponding to the target spectral line. To superimpose all target spectral lines; Step 4: Repeat steps 2-3 to obtain the spatiotemporal evolution results of the overall radiation spectrum and total light intensity during the discharge process, thus completing the calculation process.

2. The method for rapidly calculating the discharge photon radiation spectrum and total light intensity in a gaseous medium according to claim 1, characterized in that, The discharge object to be tested is a gas discharge structure under different insulation defect conditions, or a discharge system formed under different electrode forms, different insulating media, different gas pressure conditions, and different voltage excitation methods.

3. The method for rapid calculation of the discharge photon radiation spectrum and total light intensity in a gaseous medium according to claim 1, characterized in that, The numerical model of the discharge process includes the particle continuity equation and the Poisson equation. The particle continuity equation is expressed as follows: ， ， ， In the formula, the subscript , , They are electrons, cations, and anions, respectively. These represent the number densities of electrons, cations, and anions, respectively. These represent the mobility of electrons, cations, and anions, respectively. These are the diffusion coefficients for electrons, cations, and anions, respectively. For the electric field in space, For gradient operators, For photoionization source terms, The ionization coefficient, The adhesion coefficient, The recombination coefficient of electrons and cations. The recombination coefficient of anions and cations; The Poisson equation is expressed as: ， ， In the formula, This is the space electromotive force. For elementary charge, is the vacuum permittivity.

4. The method for rapid calculation of the discharge photon radiation spectrum and total light intensity in a gaseous medium according to claim 3, characterized in that, The calculation expression for the photoionization source term is: ， In the formula, Coordinates are Photoionization source term at the location, The second ionization coefficient of photoionization, for and The solid angle between them , They are respectively in Electron density and electron migration rate at that location Let be the excitation coefficient of the excited state. For the first The absorption coefficient of each excited state.

5. The method for rapid calculation of the discharge photon radiation spectrum and total light intensity in a gaseous medium according to claim 1, characterized in that, The steps for obtaining the overall radiation spectrum and total light intensity include: For any time Suppose the discharge region is discretized as Each spatial calculation point is used as a spatial location; The first Each spatial computation point at time... The local radiation spectrum below is denoted as The local radiation spectra corresponding to all spatial calculation points at the same time are summed to obtain the time. The overall radiation spectrum of the entire discharge region is expressed as: ， In the formula, For a moment The entire discharge region at wavelength The overall radiation spectrum at that location; Based on the overall radiation spectrum, the target band Integrating the light intensity within the target wavelength band yields the total light intensity, expressed as: ， In the formula, For a moment Total light intensity within the target wavelength band These are the starting and ending wavelengths of the target band; Based on the overall radiation spectrum, across the entire wavelength range Integrating the light, we obtain the total light intensity across the entire wavelength range, expressed as: ， In the formula, For a moment Total light intensity across the entire wavelength range. , These represent the lower and upper limits for the entire wavelength range.

6. The method for rapid calculation of the discharge photon radiation spectrum and total light intensity in a gaseous medium according to claim 1, characterized in that, The steps for obtaining the spatiotemporal evolution results of the overall radiation spectrum and total light intensity during the discharge process include: The total computation time of the discharge process is discretized into... Each time step is denoted as _____. ; For each moment Perform steps 2-3 to obtain the overall radiation spectrum and total light intensity; Based on the overall radiation spectrum at each moment, a sequence of the temporal variation of the overall radiation spectrum is formed. ,in, For a moment The entire discharge region at wavelength The overall radiation spectrum at that location; Based on the total light intensity at each moment, a sequence of total light intensity changing over time is generated, including the sequence of total light intensity changing over time within the target wavelength band. and the sequence of total light intensity over time across the entire wavelength range ,in, For a moment Total light intensity within the target wavelength band For a moment The total light intensity across the entire wavelength range.

7. The method for rapid calculation of the discharge photon radiation spectrum and total light intensity in a gaseous medium according to claim 1, characterized in that, The electron temperature can be replaced by one of the following: reduced electric field, electric field strength, electron density, ion concentration, and particle concentration.

8. A system for rapidly calculating the emission spectrum and total intensity of discharge photons in a gaseous medium, comprising a memory, a processor, and a program stored in the memory, characterized in that, When the processor executes the program, it implements a method for rapidly calculating the discharge photon radiation spectrum and total light intensity in a gaseous medium as described in any one of claims 1-7.

Images