Pulsed arc plasma modeling method and device based on high-speed image

Abstract

The invention provides a pulsed arc plasma modeling method and device based on a high-speed image. The method comprises the steps that pulsed arc image data, voltage and current are collected; calculating a brightness center position and a pulse arc calculation matrix; gridding the pulsed arc and obtaining the plasma density three-dimensional distribution of the pulsed arc; calculating the electron migration speed of each horizontal layer; calculating the three-dimensional distribution condition of the conductivity, and obtaining the three-dimensional distribution condition of the electron collision frequency and the electron temperature; and calculating total energy of pulsed arc particles and total macroscopic energy of the pulsed arc, and when the difference between the total energy of the particles and the total macroscopic energy meets a convergence condition, completing modeling. According to the scheme, the modeling method is simple, model parameters are more real and effective, time-varying characteristics are achieved, and application is wide.

Description

Technical Field

[0001] The present invention relates to the fields of plasma modeling and simulation, image data processing and pulse arc data modeling, and in particular to a method for constructing a three-dimensional transient plasma model of a pulse arc based on high-speed images and related equipment. Background Art

[0002] Arc discharge pulses are widely used, for example, in the fields of plasma treatment, material processing, welding, etc. In plasma treatment, arc discharge pulses can be used for surface treatment, cleaning, modification, etc. In material processing, arc discharge pulses can be used for cutting, drilling, marking, etc. In welding, arc discharge pulses can be used for welding, cutting, etc. These applications require arc discharge pulses to generate high-energy plasma to achieve specific processing effects. Modeling and analyzing pulsed arcs is an important means to explore the mechanism of pulsed arc characteristics and is a theoretical support for achieving precise control of pulsed arcs.

[0003] In the past, when modeling pulse arcs, if the finite element analysis method is used based on magnetohydrodynamic theory, although three-dimensional data can be obtained, the calculation model is complex, huge computing resources are consumed, and a large number of parameters cannot be accurately obtained in the process. In addition, most models are steady-state models or quasi-steady-state models, and it is difficult to obtain continuous changes within the pulse cycle, which is inconsistent with the rapid change characteristics of the pulse arc current. If the Meyer model or the Cauchy model is used, although the calculation method is simple, the calculated model is a one-dimensional mean model, and microscopic parameters such as plasma density and plasma temperature cannot be obtained. In addition, the above two methods are separated from the actual occurrence and development process of the pulse arc, which limits the accuracy and effectiveness of the model.

[0004] Therefore, how to establish an accurate and effective three-dimensional transient plasma model of pulsed arc is a problem that needs to be solved urgently and has great practical application significance. Summary of the invention

[0005] In view of the deficiencies in the prior art, the present invention proposes a pulse arc plasma modeling method and device based on high-speed images. Specifically, the present invention discloses the following technical solutions: In one aspect, the present invention provides a pulse arc plasma modeling method based on high-speed images, the method comprising: Step 1, initializing a high-speed camera for acquiring high-speed images, and determining a pulse arc frequency; Step 2, collecting multiple pulse arc image data within the pulse period, and collecting the voltage and current at the time of image data collection within the pulse period; Step 3: determine the pulse arc coordinate origin based on the pulse arc image data, and determine the coordinate system and the effective brightness range matrix with the pulse arc coordinate origin as the origin; perform mean processing on the effective brightness range matrix to obtain the pulse arc calculation matrix; Step 4, dividing the pulse arc into a three-dimensional structure consisting of multiple concentric rings in the horizontal direction and multiple horizontal layers in the vertical direction, so as to perform three-dimensional gridding on the pulse arc; based on the pulse arc calculation matrix and the initial value of the plasma density, obtaining the three-dimensional distribution of the pulse arc plasma density; Step 5, calculating the electron migration velocity of each horizontal layer based on the three-dimensional distribution of pulse arc plasma density; Step 6: Combine the results of step 5 to calculate the three-dimensional distribution of conductivity, and obtain the three-dimensional distribution of electron collision frequency and electron temperature; Step 7: Combine the results of step 6 to calculate the total energy of the pulse arc particles from the perspective of microscopic particles. W total , the total energy of the particle W total Includes total ion energy and total electron energy; Step 8: Combine the results of step 6 to calculate the total macro energy of the pulse arc W 等离子体 , the macroscopic total energy W 等离子体 Calculation based on energy balance and time cumulative effect; Step 9: Determine the total energy of the particle W total Total macroscopic energy W 等离子体 Whether the percentage of the difference meets the convergence condition, if it does, then end, if not, adjust the initial value of the plasma density and return to step 4 until the convergence condition is met.

[0006] Preferably, in step 3, the pulse arc coordinate origin is determined by: First calculate the brightness center position ( C x , C y ): ; ; in, H arc ( k , w ) represents the first k Column, No. w the brightness of the row; X max is the maximum horizontal pixel coordinate of the pulse arc image,Y max is the maximum value of the vertical pixel coordinate of the pulse arc image; by( C x , Y 阴 ) is the pulse arc coordinate origin, where, Y 阴 is the longitudinal coordinate of the upper surface of the cathode of the pulse arc emission source.

[0007] Preferably, in step 3, H arc ( k , w ) is calculated as follows: ;

[0008] Right now H arc ( k , w ) represents the pulse arc image k Column, No. w The brightness of the row. R ( k , w )、 G ( k , w )、 B ( k , w ) represents the RGB three-color value matrix of each pulse arc image.

[0009] Preferably, in step 3, the pulse arc calculation matrix is ​​determined as follows: First, in the coordinate system with the pulse arc coordinate origin as the origin, determine the effective brightness range matrix of the pulse arc J ( b , j ), b ∈[- H 边 ,1]∪[1, H 边 ], j ∈[1, H 纵 ],in, H 纵 Indicates the number of pixels from the upper surface of the cathode of the pulse arc emission source to the lower surface of the anode of the pulse arc emission source in the longitudinal direction; H 边 Indicates that the pulse arc coordinate origin is the center ， Brightness decay to alpha W max The maximum value of the horizontal coordinate of the pixel;W max express H arc ( k , w ) in the image; α represents the attenuation ratio; J ( b , j ) indicates the b Column, No. j The brightness value of the row; b、j Take an integer to represent the element subscript; Secondly, the effective brightness range matrix is ​​processed by averaging to obtain the pulse arc calculation matrix: K ( i , j ) = [ J (- i , j )+ J ( i , j )]´0.5; in, , i ∈[1, H 边 ], j ∈[1, H 纵 ] .

[0010] Preferably, in step 4, during the three-dimensional gridding process, the number of horizontal layers is H 纵 , the number of concentric rings is H 边 ; The electron density matrix N( i , j ) represents the three-dimensional distribution of pulse arc plasma density, N( i , j ) is calculated as: , i ∈[1, H 边 ], j ∈[1, H 纵 ]; in, n 0 is the initial value of plasma density, T abel for H 边 × H 边 The Abel coefficient matrix, K(i , j ) is the pulse arc calculation matrix.

[0011] Preferably, in step 5, the electron migration speed distribution matrix v 迁 Represents the electron migration speed of each horizontal layer, among which the vertical j Electron migration speed in horizontal layers v 迁 ( j ) is calculated as: ;

[0012] Among them, N ( i , j ) represents the electron density matrix, which represents the i Layer concentric rings, j The electron density of the horizontal layer; S ( i ) is the i The area of ​​the concentric rings, e is the electron charge, I arc It is the arc current corresponding to the moment when the pulse arc image exposure is generated.

[0013] Preferably, in step 6, the three-dimensional distribution of conductivity is reflected by the conductivity at various locations inside the pulse arc, and the conductivity at various locations inside the pulse arc is calculated as follows: ;

[0014] Among them, δ ( i , j ) indicates the j The horizontal layer i The conductivity of the concentric rings of the layer, Δ l is the actual height corresponding to the length of each pixel in the z-axis direction, N ( i , j ) represents the electron density matrix, v 迁 ( j ) is the longitudinal j The electron migration velocity distribution matrix of the horizontal layer, e is the electron charge, Δ U is the voltage drop of each vertical horizontal layer. l The value is equal to the vertical distance between the lower surface of the anode of the pulse arc emission source and the upper surface of the cathode of the pulse arc emission source and H 纵 The ratio, H 纵Indicates the number of pixels from the upper surface of the cathode of the pulse arc emission source to the lower surface of the anode of the pulse arc emission source in the longitudinal direction.

[0015] Preferably, Δ U The calculation is as follows: ; in, U arc The arc voltage at the time corresponding to the pulse arc image exposure generation, U 阴极鞘层 is the cathode sheath voltage of the pulse arc emission source, U 阳极鞘层 is the anode sheath voltage of the pulse arc emission source. U 阴极鞘层 and U 阳极鞘层 The value of is related to the electrode material and has little to do with the arc current and frequency, and can be obtained by looking up the table.

[0016] Preferably, in step 6, the three-dimensional distribution of the electron collision frequency and the electron temperature is obtained by calculating the electron collision frequency and the electron temperature based on the following conductivity characteristics: ;

[0017] in, f ei ( i , j ) is the electron collision frequency, lnΛ is the Coulomb logarithm, m e is the electron mass, Z i is the relative charge number of plasma, T e ( i , j ) is the electron temperature.

[0018] Preferably, in step 7, the total energy of the pulse arc particles is W total The calculation method is: ;

[0019] in, i ∈[1, H 边 ], j ∈[1, H 纵 ], N ( i , j ) represents the electron density matrix, T e ( i ,j ) is the electron temperature, Z i is the relative charge number of plasma, V arc ( i , j ) represents the first j The horizontal layer i The volume of the concentric rings can be combined Figure 5 The pulse arc structure diagram and the actual size represented by each pixel are calculated.

[0020] Preferably, in step 8, the pulse arc macroscopic total energy W 等离子体 The calculation method is: ; ; ; ;

[0021] in, W 等离子体 ( t ) is the total energy of the pulse arc at the current time t, W 等离子体 ( t -1) is the total energy of the pulse arc at the previous moment, Δ W 等离子体 ( t ) is the change in energy within this time step, W input ( t ) is the input energy, W loss ( t ) is the energy loss, W f ( t ) is the energy lost by the particle, W s ( t ) is the radiation loss energy, Δ t is the time step value, U arc ( t ) is the arc voltage at the current time t, U 阴极鞘层 is the cathode sheath voltage of the pulse arc emission source, U 阳极鞘层 is the anode sheath voltage of the pulse arc emission source.

[0022] Preferably, the particle loses energy W f (t ) is solved as follows: the statistical range of ion loss energy includes the layer corresponding to the upper surface, the layer corresponding to the lower surface and the layer corresponding to the side of the cylinder. W f ( t ) is calculated as follows: ; Among them, Shell represents the range composed of the layer corresponding to the upper surface of the pulse arc, the layer corresponding to the lower surface and the layer corresponding to the side of the cylinder. is the electron flux density of the outermost layer of the pulse arc, is the ion flux density of the outermost layer of the pulse arc, T e ( i , j ) is the electron temperature, Indicates the surface area corresponding to the outermost position of the pulse arc.

[0023] Preferably, W s ( t ) according to the energy loss W loss ( t ) is calculated as 15% of the

[0024] On the other hand, the present invention also provides a pulse arc plasma modeling device based on high-speed images, the device comprising: a timing control module, a pulse arc oscillation circuit, a high-speed camera, a host computer and a voltage and current measuring device; The pulse arc oscillation circuit generates a pulse arc based on a timing control module, the high-speed camera is used to collect multiple pulse arc image data within a pulse period and upload them to a host computer, and the voltage and current measuring device is used to collect the voltage and current at the time of image data collection within the pulse period and upload them to the host computer; The host computer at least includes a processor and a memory, and the processor calls computer instructions in the memory to execute the pulse arc plasma modeling method based on high-speed images as described above.

[0025] Preferably, the pulse arc oscillation generating circuit is composed of a pulse energy circuit, an arc striking circuit, an absorption protection circuit and a pulse arc emission source.

[0026] Preferably, the pulse energy circuit consists of a resistor R 0 ,inductance L 1 and capacitor C 0 and bidirectional triode thyristor VT 0 The bidirectional triode thyristor VT 0 Responsible for the conduction control of the control circuit; bidirectional triode thyristor VT 0 The T1 pole connection inductance L 1 one end of the inductor; L 1 The other end of the capacitor is connected C 0 The positive electrode of the capacitor C 0 Negative connection resistance R 0 one end of the resistor R 0 The other end of the bidirectional three-stage thyristor is grounded. VT 0 The T2 pole is connected to the positive pole of the pulse arc emission source.

[0027] Preferably, the arc striking circuit is composed of an electrolytic capacitor C 1 , Current limiting resistor R 1 , thyristor VT 1 And power diode VD 1 Composition: Thyristor VT 1 The positive terminal is connected to the current limiting resistor R 1 One end of the current limiting resistor R 1 The other end of the electrolytic capacitor is connected C 1 The positive electrode and power diode VD 1 The negative electrode of the electrolytic capacitor C 1 The cathode of the power diode VD 1 The positive pole of the thyristor is grounded. 1 The negative electrode is connected to the positive electrode of the pulse arc emission source.

[0028] Preferably, the absorption protection circuit includes a capacitor C 2 ,resistance R 2 and varistor R v ;resistance R 2 With capacitor C 2 Series, varistor R v With resistance R 2 and capacitor C 2The series circuit is connected in parallel, that is, the varistor R v Connect one end of the resistor R 2 Keep away from the capacitor C 2 One end of the varistor R v The other end of the capacitor is grounded; C 2 Keep away from the resistor R 2 One end is grounded.

[0029] Compared with the prior art, this solution has at least the following advantages: 1. The modeling method is simple, and the plasma density, temperature and electron migration velocity distribution results of the three-dimensional pulse arc can be obtained without using the finite element analysis method; 2. The modeling process of this model is based on the real data of arc images, pulse arc voltage and current, so the model parameters are more real and effective; 3. This model is time-varying and is a dynamic model within the entire cycle of the pulse arc. The obtained microscopic parameters also change dynamically. Theoretically, as long as the shooting speed is fast enough and the voltage and current sampling accuracy is high enough, even if the pulse arc frequency is very high, a non-steady-state three-dimensional model can be established. BRIEF DESCRIPTION OF THE DRAWINGS

[0030] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0031] Figure 1 A schematic diagram of a model building method flow chart of an embodiment of the present invention; Figure 2 It is a schematic diagram of the structure of a pulse arc generation measurement system according to an embodiment of the present invention; Figure 3 A timing diagram of a pulse arc generation measurement system according to an embodiment of the present invention; Figure 4 The principle of generating high-speed images of pulse arcs in a certain horizontal direction according to an embodiment of the present invention; Figure 5 A schematic diagram of a multi-layer structure of a pulse arc according to an embodiment of the present invention; Figure 6 This is a schematic diagram of a 3kA pulse arc image according to an embodiment of the present invention; Figure 7It is a schematic diagram of the three-dimensional distribution result of the electron density corresponding to the 3kA pulse arc of an embodiment of the present invention; Figure 8 It is a schematic diagram of the three-dimensional calculation results of the conductivity corresponding to the 3kA pulse arc according to an embodiment of the present invention; Fig. 9 This is a schematic diagram of the electron migration velocity along the z-axis direction corresponding to a 3kA pulse arc according to an embodiment of the present invention; Fig.10 It is a schematic diagram of the three-dimensional distribution result of electron temperature corresponding to the 3kA pulse arc of an embodiment of the present invention. DETAILED DESCRIPTION

[0032] The embodiments of the present invention are described in detail below in conjunction with the accompanying drawings. It should be clear that the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0033] Those skilled in the art should know that the following specific embodiments or specific implementations are a series of optimized settings listed by the present invention to further explain the specific content of the invention, and these settings can be combined or used in association with each other, unless the present invention clearly states that some or a specific embodiment or implementation cannot be associated or used together with other embodiments or implementations. At the same time, the following specific embodiments or implementations are only used as the most optimized settings, and are not to be understood as limiting the scope of protection of the present invention.

[0034] In a specific embodiment, the present invention establishes a pulse arc microscopic parameter model with time-varying characteristics and three-dimensional distribution based on the arc image and voltage and current parameters within the pulse arc cycle. The specific process of the modeling method is as follows: Figure 1 shown.

[0035] In this embodiment, the generation and parameter measurement of the pulse arc rely on a pulse arc generation measurement system, which consists of a pulse arc oscillation circuit, a high-speed camera, an oscilloscope, a host computer and a voltage and current measurement device. The circuit topology is as follows: Figure 2 The timing control module is responsible for controlling the pulse arc oscillation circuit, camera, and oscilloscope action timing. The specific timing control is as follows: Figure 3 The pulse arc is generated by the pulse arc oscillation circuit. The voltage and current signals of the pulse arc are collected by the oscilloscope and uploaded to the host computer system. The high-speed camera takes the high-speed image of the pulse arc. The host computer system integrates the pulse arc image and the voltage and current signals to model the pulse arc.

[0036] In a more preferred embodiment, in combination Figure 2 As shown, the pulse arc oscillation generating circuit is composed of a pulse energy circuit, an arc starting circuit, an absorption protection circuit and a pulse arc emission source.

[0037] For the pulse energy circuit, the resistor R 0 ,inductance L 1 and capacitor C 0 and bidirectional triode thyristor VT 0 , where the bidirectional triode thyristor VT 0 Responsible for the conduction control of the control circuit. VT 0 The T1 pole connection inductance L 1 one end of the inductor; L 1 The other end of the capacitor is connected C 0 The positive electrode of the capacitor C 0 Negative connection resistance R 0 one end of the resistor R 0 The other end of the bidirectional three-stage thyristor is grounded. VT 0 The T2 pole is connected to the positive pole of the pulse arc emission source.

[0038] Combined with the circuit structure of the pulse energy circuit, it can be seen that the resistance R 0 ,inductance L 1 and capacitor C 0 The voltage relationship satisfies equations (1) and (2): (1) (2) Among them, U c Indicates capacitance C 0 Voltage, U R Represents resistance R 0 Voltage, U L Inductance L 1 Voltage. When the resistance of the circuit R 0 ,inductance L 1 and capacitorC 0 When the relationship in equation (3) is satisfied, the main circuit generates a non-oscillating discharge process, and the oscillating current i(t) can be described as equation (4): U 0 is the pre-charge voltage of the main circuit capacitor, ω is the angular frequency of the second-order oscillator circuit, δ is the time constant, and the calculation formulas are: (3) (4) (5) (6) The arc-starting circuit is used to generate a low-current DC arc so that the medium in the electrode gap remains conductive before the pulse current is applied. In this embodiment, a specific circuit design is used as an example. The arc-starting circuit consists of a 50mF / 220V electrolytic capacitor. C 1 , current limiting resistor with a resistance of 2.5Ω R 1 , thyristor VT with capacity of 1800V / 200A 1 And power diode VD 1 Composition. Thyristor VT 1 The positive terminal is connected to the current limiting resistor R 1 One end of the current limiting resistor R 1 The other end of the electrolytic capacitor is connected C 1 The positive pole and power diode VD 1 Negative electrode of electrolytic capacitor C 1 The cathode of the power diode VD 1 The positive pole of the thyristor is grounded. 1 The negative electrode is connected to the positive electrode of the pulse arc emission source. 1 When conducting, the electrolytic capacitor C 1 Through the current limiting resistor R 1 Discharge, so that current flows through the positive and negative electrodes of the pulse spectrum emission source. Power diode VD 1 It protects the electrolytic capacitor. The time constant of the discharge process is 125ms and the maximum current is 80A.

[0039] For the absorption protection circuit, including capacitor C 2 ,resistance R 2 and varistorR v .resistance R 2 With capacitor C 2 Series, varistor R v With resistance R 2 and capacitor C 2 The series circuit is connected in parallel, that is, the varistor R v Connect one end of the resistor R 2 Keep away from the capacitor C 2 One end of the varistor R v The other end of the capacitor is grounded; C 2 Keep away from the resistor R 2 One end of the capacitor is grounded. C 2 and resistor R 2 Used to adjust the recovery voltage rise rate d after the pulse current passes zero U / d t , Varistor R v It is used to limit the amplitude of the recovery voltage. The absorption protection circuit can prevent the positive and negative electrodes of the pulse spectrum source from being broken down to avoid the emission of unplanned spectrum, and it also has a withstand voltage protection effect on each device in the line.

[0040] The pulse arc emission source is mainly composed of an electromagnetic action mechanism and positive and negative electrodes. The electromagnetic action mechanism is operated by a timing control module according to a preset timing logic (refer to Figure 3 As shown), to control the arc between the positive and negative electrodes.

[0041] In this embodiment, combined with the above-mentioned pulse arc generation measurement system, the overall process of the pulse arc plasma modeling method is as follows: Figure 1 As shown, the details are as follows: STEP 1: Determine the frequency of the pulse arc and adjust the inductance of the pulse energy loop L 1 and capacitor C 0 value; determine the current level of the pulse arc, and set the capacitor C 0 Charge to preset value; arc-starting circuit capacitance C 1Charge to 200V; set the shooting speed of the high-speed camera, adjust the camera lens to face the center of the gap between the positive and negative electrodes of the pulse arc emission source, and adjust the camera aperture so that the R, G, and B values ​​of the captured arc image are all less than the maximum digit value (for example, if the camera is 8-bit, the maximum value is 255; if the camera is 16-bit, the maximum value is 65535). If the camera aperture value is adjusted to the minimum and still exceeds the maximum digit, add a grayscale film in front of the camera lens, and the grayscale film has a consistent percentage of radiation intensity attenuation in the visible light range.

[0042] STEP 2: Press Figure 3 The timing logic shown in the figure is used to perform pulse arc experiment. The host computer collects the pulse period. p A pulse arc image ".RAW" format file, p The value is determined by the pulse frequency and the camera shooting speed; the oscilloscope acquires the voltage and current values ​​within the pulse arc cycle. It is understandable that the specific file format of the pulse arc image can adopt other common types except the ".RAW" format, and the file format here is not understood as limiting the scope of protection of the present invention.

[0043] Combination Figure 3 As shown, the timing in this embodiment is set as follows: thyristor VT 1 The trigger is used as a timing start signal. After a time t 1 , the electromagnetic action mechanism in the pulse arc emission source is triggered, and the electrode gap is turned on; at time t 2 Node, high-speed camera is triggered, camera CCD is started, and then to t 3 Node, the bidirectional triode thyristor VT 0 Trigger, pulse arc is put into operation; until time node t 4 When the bidirectional triode thyristor VT 0 Turn off and the pulse arc is extinguished.

[0044] We take 3kA pulse arc as an example, and the pulse arc image obtained is as follows Figure 6 shown.

[0045] STEP 3: By reading the ".RAW" format file, the RGB three-color value matrix of each pulse arc image can be obtained, specifically the matrix R ( k , w )、 G ( k , w )、 B ( k , w ),in, k ∈[1, Xmax ], w ∈[1, Y max ], X max is the maximum horizontal pixel coordinate of the arc image, Y max is the maximum value of the vertical pixel coordinate of the arc image. The origin of the coordinate can generally be the lower left corner (or upper right corner) of the image. Taking a 320 (total number of columns) ´240 (total number of rows) pixel matrix as an example, the origin is the pixel corresponding to the 1st column and the 1st row, and the origin coordinate is (1,1). X max is 320, Y max is 240, then R ( k , w ) also represents the k Column, No. w The pixel corresponding to the row R value, G ( k , w )and B ( k , w ) has a similar meaning. The grayscale of each pulse arc image is calculated. In this embodiment, the calculation standard refers to the photometric calculation standard of "Photometric / Digital ITU BT.709", and the pixel position is ( k , w ) of the arc brightness pixel matrix H arc ( k , w ) is calculated as follows: (7) Right now H arc ( k , w ) represents the pulse arc image k Column, No. w The brightness of the line.

[0046] Since the starting position of the pulse arc is random, it is necessary to determine the pulse arc brightness center in order to select the valid data range. H arc ( k , w ), calculate the pixel position of the brightness center position as ( C x , C y ), the calculation formula is as follows: (8) (9) Determine the effective brightness range of the pulse arc: The basic characteristic of the pulse arc image is that its brightness gradually decays from the center of the arc to both sides. We first determine the pulse arc coordinate origin, and then perform subsequent pulse arc image processing in a new coordinate system with the pulse arc coordinate origin as the origin. Then the effective brightness range can be determined as: ( C x , Y 阴 ) is the pulse arc coordinate origin O ; The range in the longitudinal direction is the cathode (i.e. Figure 2 The negative electrode in the Figure 2 The number of pixels on the surface below the positive electrode H 纵 , Y 阴 express H arc ( k , w ) corresponds to the vertical pixel coordinate of the arc area closest to the cathode upper surface; the horizontal range is (- H 边 , H 边 ), H 边 Origin O Centered ， The arc brightness decays to W max The horizontal coordinate of the pixel where the attenuation ratio α=10% is located. W max Represents the arc brightness pixel matrix H arc The maximum brightness value, that is, the maximum brightness value in the pulse arc image. It can be understood here that the attenuation ratio α can be selected and changed, for example, adjusted to 8%, 15%, etc. This embodiment takes 10% as an example.

[0047] Then, the pulse arc effective brightness range matrix can be obtained J ( b , j ), b ∈[- H 边 ,1]∪[1, H 边 ], j ∈[1, H 纵 ], b ,j is the pixel number, which is an integer and not equal to 0. From the above derivation process, it can be seen that the matrix J The value of each element in represents the pixel brightness at the corresponding position. The effective brightness range matrix J Has been converted to the pulse arc coordinate origin O In the coordinate system with as the origin.

[0048] In order to reduce the impact of asymmetric distribution of pulse arc, the arc brightness of the left and right halves is averaged, that is, the origin is O The pulse arc calculation matrix covering half of the arc pulse image is obtained by taking the average brightness of two pixels symmetrical in the horizontal direction as the center. K : K ( i , j ) = [ J (- i , j )+ J ( i , j )]´0.5, , i ∈[1, H 边 ], j ∈[1, H 纵 ] .

[0049] like C x、 C y、 Y 阴、 H 边 If it is not an integer, it will be rounded up. i and j Refers to the matrix element number, which is the corresponding pixel number in this embodiment. It can only be an integer, not a decimal. Taking this as an example, the pulse arc calculation matrix K The horizontal range is from 1st to H 边 Pixels, vertically from 1st to H 纵 Pixels, therefore, in this embodiment, if C x、 C y、 Y 阴、 H 边 If the number is not an integer, it will be rounded off.

[0050] STEP 4: The pulse arc can be considered as a columnar body with symmetrical distribution of characteristics along the center of the arc. Its contribution to the brightness of each pixel of the arc image is proportional to the cumulative value of the plasma density corresponding to the pixel along a certain direction. In this embodiment, the columnar body of the arc is gridded to facilitate the subsequent calculation of the three-dimensional distribution, that is, the three-dimensional grid division is performed based on the pulse arc calculation matrix in STEP3, refer to Figure 4 , Figure 5 As shown, in this embodiment, the column is divided into concentric rings in the horizontal direction, and grid division is performed according to equal-height cylinders in the vertical direction, that is, the pulse arc column can be regarded as a three-dimensional structure with m horizontal layers in the vertical direction (z-axis direction) and n concentric rings in the horizontal direction (xOy plane). In this embodiment, the grid division is divided according to the pixel as the unit, and m is H 纵 , n is H 边 , that is, after gridding, the number of concentric rings is H 边 , the number of horizontal layers is H 纵 .

[0051] In the following description, the concentric rings are numbered from the inside to the outside, starting with i Indicates the number of pulse arc cylinders from the inside to the outside. i Layer concentric rings (reference Figure 4 ), i The value range is [1, H 边 ], j Indicates the pulse arc cylinder from bottom to top j Horizontal layer, j The value range is [1, H 纵 It can be understood that the horizontal layer starting position of the pulse arc cylinder can be the uppermost or lowermost end of the pulse arc. In this embodiment, the pulse arc calculation matrix K The position corresponding to the vertical pixel coordinate of 1 is the starting position of the horizontal layer, that is, the pulse arc coordinate origin O The layer is the starting layer; in addition, the innermost 1st layer ring can be understood as a complete circle, or as a ring excluding the point on the coordinate y-axis. These understandings do not affect the calculation and use of this solution.

[0052] In this embodiment, the pulse arc calculation matrix K ( i,j ) performs an inverse Abel transform to obtain the matrix N( i,j ), N( i,j ) is the electron density matrix, which is calculated as shown in formula (10), whereT abel for H 边 × H 边 The Abel coefficient matrix can be obtained by looking up the table.

[0053] N ( i , j )= n 0 K ( i , j ) T abel (10) In the above formula, n 0 is the initial value of the plasma density, then the matrix is the electron density matrix of the pulse arc, where n = H 边 ,m= H 纵 ,in n The value changes from small to large, corresponding to the horizontal change from the center of the pulse arc to the surrounding area; m The value increases from small to large, corresponding to the vertical change from the cathode upper surface to the anode lower surface. i, j ) represents the first i Concentric rings, j The electron density of the horizontal layer.

[0054] Continuing with the 3kA pulse arc in this embodiment as an example, the corresponding three-dimensional distribution result of the electron density is as follows: Figure 7 shown.

[0055] STEP 5: After 3D meshing, the pulse arc column can be regarded as a 3D structure with m layers in the longitudinal direction (z-axis direction) and n layers in the horizontal direction (xOy plane), such as Figure 5 shown.

[0056] Voltage drop of each vertical layer Δ U are equal, and the calculation formula is as follows: (11) in, U arc The arc voltage at the time corresponding to the pulse arc image exposure generation, U 阴极鞘层 is the cathode sheath voltage of the pulse arc emission source, U 阳极鞘层is the anode sheath voltage of the pulse arc emission source, U 阴极鞘层 and U 阳极鞘层 The value of is related to the electrode material and has little to do with the arc current and frequency, and can be obtained by looking up the table.

[0057] According to the law of conservation of current, the current in each longitudinal layer is also equal. I arc The arc current at the corresponding moment of the pulse arc image exposure is generated. The pulse arc voltage and current data are collected by the oscilloscope and transmitted to the host computer, which interpolates the pulse arc image exposure time to obtain U arc , I arc .

[0058] The total current value passing through is the same and can be obtained by measurement. It can be understood that the total current value is the current value of the cathode of the pulse arc emission source, and the current value of each xOy horizontal layer is equal. In this embodiment, the current density calculation formula is as follows: (12) in, Q(i,j) represents the current density, N(i,j) is the electron density matrix (i.e., plasma density), e is the electron charge, v 迁 ( j ) is the electron migration speed.

[0059] For a certain vertical number j For the horizontal layer, the current calculation formula is as follows: (13) in, S ( i ) is on the horizontal plane xOy, i The area of ​​the concentric rings, H 边 Indicated by ( C x , Y 阴 ) as the center ， The effective brightness range decays to α W max The maximum value of the horizontal coordinate of the pixel, W max Indicates the effective brightness range J ( b , j ). According to formula (13), the vertical number can be calculated as jThe electron migration speed of the horizontal layer v 迁 ( j ),Right now: .

[0060] Furthermore, we can find m Electron migration speed of the layer v 迁 ( j ), j ∈[1, m ].

[0061] Still taking the 3kA pulse arc in this embodiment as an example, the calculation result of the electron migration velocity along the z-axis direction is as follows: Fig. 9 shown.

[0062] STEP 6: For a certain vertical number j For the horizontal layer, i The conductivity calculation formula of the concentric rings is as follows: (14) Among them, Δ l The actual height corresponding to the length of each pixel in the z-axis direction, its value is equal to L gap / m ,in L gap It represents the vertical distance between the lower surface of the positive electrode (anode) and the upper surface of the negative electrode (cathode).

[0063] According to the definition of conductivity in formula (15), the electron collision frequency can be further calculated as f ei ( i , j ) and the electron temperature T e ( i , j ) distribution: (15) in f ei ( i , j ) is the electron collision frequency, unit Hz; lnΛ is the Coulomb logarithm, which can be obtained by looking up the table; m e is the mass of the electron, which is 9.11×10 -31 kg; Z i is the relative charge number of the plasma, which can be obtained by looking up the table; T e ( i ,j ) is the electron temperature, unit is eV.

[0064] Still taking the 3kA pulse arc in this embodiment as an example, the corresponding three-dimensional distribution of conductivity is calculated as follows: Figure 8 The three-dimensional distribution of electron temperature is shown in Fig.10 shown.

[0065] STEP 7: Based on the above results, calculate the total energy from the perspective of microscopic particles inside the pulse arc W total ,by W 离子 is the total ion energy, W 电子 represents the total energy of the electron, N ( i , j ) represents the electron density matrix, T e ( i , j ) is the electron temperature matrix, Z i is the relative charge number of plasma, V arc ( i , j ) represents the first j The horizontal layer i The volume of the concentric rings can be combined Figure 5 The total energy is calculated from the pulse arc structure diagram and the actual size represented by each pixel. W total As shown in formula (16): (16) in, i ∈[1, H 边 ], j ∈[1, H 纵 ], i , j Take the integer.

[0066] STEP 8: Calculate the total energy of the pulse arc based on energy balance and time accumulation effect. Suppose the total energy of the pulse arc at the current time t is W 等离子体 ( t ), then it is equal to the total energy at the previous moment W 等离子体 ( t -1) plus the change in energy during this time step Δ W 等离子体 ( t ), the change ΔW 等离子体 ( t ) is the input energy W input ( t ) and lost energy W loss ( t ), where the input energy is calculated as (19), and the loss energy W loss ( t ) includes the energy lost by the particles W f ( t ) and radiation loss energy W s ( t ), as shown in formula (20). t is the time step value, that is, the shooting interval between two images of the high-speed camera.

[0067] (17) (18) (19) (20) The particle loses energy W f ( t ) is caused by the outermost particles (including electrons and ions) of the pulse arc moving outside the boundary. Therefore, its statistical range includes the layer corresponding to the upper surface ( H 纵 The corresponding layer of the lower surface (the first horizontal layer) and the corresponding layer of the cylinder side (the first horizontal layer of the entire cylinder) H 边 The shell range is the outermost layer of the pulse arc, that is, the values ​​of i and j are within the shell range (which can be expressed as ( i, j )∈Shell) is calculated by the following formula (21): (twenty one) Among them, Shell represents the range composed of the layer corresponding to the upper surface of the pulse arc, the layer corresponding to the lower surface and the layer corresponding to the side of the cylinder. is the electron flux density in the outermost layer (i.e., Shell range) of the pulse arc, and the calculation formula is as follows (22), where For the position ( i , j ) at which the electron escape velocity is is the ion flux density of the outermost layer (i.e., Shell range) of the pulse arc, and the calculation formula is as follows (23), where For the position ( i , j ) at which the ion escape velocity is m ion is the ion mass, T e ( i , j ) is the electron temperature, Indicates the surface area corresponding to the outermost position of the pulse arc.

[0068] (twenty two) (twenty three) in, N ( i , j ) means at position ( i , j ), W s ( t ) represents the energy dissipated by radiation, conduction and other forms. This part of energy is the total dissipated energy (i.e., lost energy). W loss ). It can be understood that the particle loss energy given in this step is W f ( t ) is a preferred implementation method, the particle loss energy W f ( t ) can also be calculated by other estimation methods or table lookup methods, and other conventional calculation methods can also achieve the calculation effect of this solution to a certain extent. Therefore, they should all be regarded as falling within the protection scope of the present invention.

[0069] STEP 9: Judgement W total and W 等离子体 ( t ) Whether the percentage difference meets the convergence condition b , in this example b Take 5%. If it is satisfied, go to step STEP10. The optimal calculation method for convergence judgment is as follows: ; If the convergence condition is not met, adjust the initial value of plasma density n 0 , and start the calculation from STEP 4 again until the convergence condition is met.

[0070] STEP 10: Proceed to the next step of simulation modeling, or end the simulation.

[0071] The modeling method of the scheme provided in this embodiment is simple, and the plasma density, temperature and electron migration velocity distribution results of the three-dimensional pulse arc can be obtained without using the finite element analysis method; during the modeling process of this scheme, the model parameters are more real and effective based on the real data of the captured arc image and the pulse arc voltage and current; the model constructed by this scheme is time-varying and is a dynamic model within the entire cycle of the pulse arc. The obtained microscopic parameters also change dynamically. In theory, as long as the shooting speed is fast enough and the voltage and current sampling accuracy is high enough, even if the pulse arc frequency is very high, a non-steady-state three-dimensional model can be established.

[0072] It should be noted that any process or method description in the flowchart or otherwise described herein may be understood to represent a module, fragment or portion of a code comprising one or more executable instructions for implementing steps of a specific logical function or process, and the scope of the preferred implementation of the present scheme includes alternative implementations, in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the implementation of the present scheme belongs.

[0073] The logic and / or steps represented in the flowchart or otherwise described herein may be embodied in any readable storage medium for use by an instruction execution system, apparatus, or device (such as a computer-based system, a system including a processor, or other system that can fetch instructions from and execute instructions on an instruction execution system, apparatus, or device), or in combination with such instruction execution systems, apparatuses, or devices.

[0074] The above is only a specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed by the present invention should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims

1. A pulse arc plasma modeling method based on high-speed images, characterized in that: The method comprises: Step 1, initializing a high-speed camera for acquiring high-speed images, and determining a pulse arc frequency; Step 2, collecting multiple pulse arc image data within the pulse period, and collecting the voltage and current at the time of image data collection within the pulse period; Step 3: determine the pulse arc coordinate origin based on the pulse arc image data, and determine the coordinate system and the effective brightness range matrix with the pulse arc coordinate origin as the origin; perform mean processing on the effective brightness range matrix to obtain the pulse arc calculation matrix; Step 4, dividing the pulse arc into a three-dimensional structure consisting of multiple concentric rings in the horizontal direction and multiple horizontal layers in the vertical direction, so as to perform three-dimensional gridding on the pulse arc; based on the pulse arc calculation matrix and the initial value of the plasma density, obtaining the three-dimensional distribution of the pulse arc plasma density; Step 5, calculating the electron migration velocity of each horizontal layer based on the three-dimensional distribution of pulse arc plasma density; Step 6: Combine the results of step 5 to calculate the three-dimensional distribution of conductivity, and obtain the three-dimensional distribution of electron collision frequency and electron temperature; Step 7: Combine the results of step 6 to calculate the total energy of the pulse arc particles from the perspective of microscopic particles. W total , the total energy of the particle W total Includes total ion energy and total electron energy; Step 8: Combine the results of step 6 to calculate the total macro energy of the pulse arc W 等离子体 , the macroscopic total energy W 等离子体 Calculation based on energy balance and time cumulative effect; Step 9: Determine the total energy of the particle W total Total macroscopic energy W 等离子体 Whether the percentage of the difference meets the convergence condition, if it does, then end, if not, adjust the initial value of the plasma density and return to step 4 until the convergence condition is met.

2. The method according to claim 1, characterized in that In step 3, the pulse arc coordinate origin is determined as follows: First calculate the brightness center position ( C x , C y ): ； ； in, H arc ( k , w ) represents the first k Column, No. w the brightness of the row; X max is the maximum horizontal pixel coordinate of the pulse arc image, Y max is the maximum value of the vertical pixel coordinate of the pulse arc image; by( C x , Y 阴 ) is the pulse arc coordinate origin, where, Y 阴 is the longitudinal coordinate of the upper surface of the cathode of the pulse arc emission source.

3. The method according to claim 2, characterized in that In step 3, the pulse arc calculation matrix is ​​determined as follows: First, in the coordinate system with the pulse arc coordinate origin as the origin, determine the effective brightness range matrix of the pulse arc J ( b , j ), b ∈[- H 边 ,1]∪[1, H 边 ], j ∈[1, H 纵 ],in, H 纵 Indicates the number of pixels from the upper surface of the cathode of the pulse arc emission source to the lower surface of the anode of the pulse arc emission source in the longitudinal direction; H 边 Indicates that the pulse arc coordinate origin is the center ， Brightness decay to alpha W max The maximum value of the horizontal coordinate of the pixel; W max express H arc ( k , w ) in the image; α represents the attenuation ratio; J ( b , j ) represents the brightness value of the b-th column and the j-th row; b、j Take the integer; Secondly, the effective brightness range matrix is ​​processed by averaging to obtain the pulse arc calculation matrix: K （ i , j ）= [ J （- i , j ）+ J （ i , j ）]×0.5； in, , i ∈[1, H 边 ], j ∈[1, H 纵 ] .

4. The method according to claim 3, characterized in that In step 4, during the three-dimensional meshing process, the number of horizontal layers is H 纵 , the number of concentric rings is H 边 ; The electron density matrix N( i , j ) represents the three-dimensional distribution of pulse arc plasma density, N( i , j ) is calculated as: ， i ∈[1， H 边 ]， j ∈[1， H 纵 ]； in, n 0 is the initial value of plasma density, T abel for H 边 × H 边 The Abel coefficient matrix, K( i , j ) is the pulse arc calculation matrix.

5. The method according to claim 4, characterized in that In step 5, the electron migration speed distribution matrix v 迁 Represents the electron migration speed of each horizontal layer, among which the vertical j Electron migration speed in horizontal layers v 迁 ( j ) is calculated as: ； Among them, N ( i , j ) is the electron density matrix, S ( i ) is the i The area of ​​the concentric rings, e is the electron charge, I arc It is the arc current corresponding to the moment when the pulse arc image exposure is generated.

6. The method according to claim 4, characterized in that In step 6, the three-dimensional distribution of conductivity is represented by the conductivity at each location inside the pulse arc, and the conductivity at each location inside the pulse arc is calculated as follows: ； Among them, δ ( i , j ) indicates the j The horizontal layer i The conductivity of the concentric rings of the layer, Δ l is the actual height corresponding to the length of each pixel in the z-axis direction, N ( i , j ) represents the electron density matrix, v 迁 ( j ) is the longitudinal j The electron migration speed of the horizontal layer, e is the electron charge, Δ U is the voltage drop of each vertical horizontal layer.

7. The method according to claim 6, characterized in that In step 6, the three-dimensional distribution of the electron collision frequency and the electron temperature is obtained by calculating the electron collision frequency and the electron temperature based on the following conductivity characteristics: ； in, f ei ( i , j ) is the electron collision frequency matrix, lnΛ is the Coulomb logarithm, m e is the electron mass, Z i is the relative charge number of plasma, T e ( i , j ) is the electron temperature matrix.

8. The method according to claim 4, characterized in that In step 7, the total energy of the pulse arc particles W total The calculation method is: ; Among them, N ( i , j ) represents the electron density matrix, T e ( i , j ) is the electron temperature matrix, Z i is the relative charge number of plasma, V arc ( i , j ) indicates the j The horizontal layer i The volume of the concentric rings.

9. The method according to claim 1, characterized in that: In step 8, the pulse arc macroscopic total energy The calculation method is: ; ; ; ; in, W 等离子体 ( t ) is the total energy of the pulse arc at the current time t, W 等离子体 ( t -1) is the total energy of the pulse arc at the previous moment, Δ W 等离子体 ( t ) is the change in energy within this time step, W input ( t ) is the input energy, W loss ( t ) is the energy loss, W f ( t ) is the energy lost by the particle, W s ( t ) is the radiation loss energy, Δ t is the time step value, U arc ( t ) is the arc voltage at the current time t, U 阴极鞘层 is the cathode sheath voltage of the pulse arc emission source, U 阳极鞘层 is the anode sheath voltage of the pulse arc emission source.

10. Pulse arc plasma modeling equipment based on high-speed images, characterized in that: The equipment includes: a timing control module, a pulse arc oscillation circuit, a high-speed camera, a host computer and a voltage and current measuring device; The pulse arc oscillation circuit generates a pulse arc based on a timing control module, the high-speed camera is used to collect multiple pulse arc image data within a pulse period and upload them to a host computer, and the voltage and current measuring device is used to collect the voltage and current at the time of image data collection within the pulse period and upload them to the host computer; The host computer includes at least a processor and a memory, and the processor calls computer instructions in the memory to execute the pulse arc plasma modeling method based on high-speed images as described in any one of claims 1-9.

Images