Induced eddy current calculation model construction method, calculation method and electronic equipment

Abstract

The invention discloses an induced eddy current calculation model construction method, a calculation method and electronic equipment, and is suitable for calculation of induced eddy current of a vacuum wall in a magnetic confinement nuclear fusion device.The model construction method comprises the steps that calibration discharge data of a plurality of coil calibration cannons are obtained, the calibration discharge data comprises coil current calibration data and magnetic diagnosis calibration data of magnetic diagnosis equipment; determining induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data; coil current calibration data is used as an input quantity, induced eddy current calibration data is used as an output quantity, and an induced eddy current calculation model constructed based on a transfer function is used for fitting a response coefficient between coil current and induced eddy current to serve as a model parameter of the calculation model. Based on the transmission function, the response coefficient between the coil current and the induced eddy current is directly fitted in a data processing mode from real data, and the induced eddy current of the coil is efficiently and accurately calculated.

Description

Technical Field

[0001] This application relates to the field of nuclear fusion technology, and particularly relates to a method for constructing a calculation model of induced eddy current, a calculation method, and an electronic device. Background Art

[0002] A magnetic confinement nuclear fusion device is a device that confines and controls high-temperature plasma through a magnetic field to achieve controlled nuclear fusion. Taking a tokamak device as an example, the vacuum vessel wall of the tokamak device is usually composed of a metal material. During the operation of the tokamak device, the changing currents in various coils (such as poloidal field coils, toroidal field coils, central solenoid coils, etc.) will generate induced eddy currents on the vacuum vessel wall. The existence of the induced eddy currents has various impacts on the operation of the device, such as changing the shape of the breakdown region, changing the plasma configuration, and introducing errors into the acquisition results of various diagnostic measurement devices. Therefore, quantitatively analyzing the distribution of these induced eddy currents, especially quickly obtaining the distribution of induced eddy currents during the discharge operation of the device, is of great significance for the normal operation of the tokamak device and the research of many physical problems.

[0003] However, due to the continuous and closed structure of the vacuum vessel wall, it is difficult to obtain the distribution of induced eddy currents through actual measurement. Therefore, in related technologies, computer simulation is often used to obtain the induced eddy currents on the vacuum vessel wall.

[0004] In one of the related technologies, the L-R differential equation solving method is adopted. Specifically, the vacuum vessel wall is divided into numerous small blocks. Since the eddy current satisfies the L-R differential equation, the differential equation is solved to calculate the eddy current. However, this method has a large amount of calculation and a long calculation time.

[0005] In another related technology, the equilibrium inversion method is used to solve the induced eddy current, specifically as follows: equations are constructed using the coil current, plasma current, and the relationship between the induced eddy current and the magnetic diagnostic device to obtain the unknowns: the induced eddy current and the plasma current. Since the plasma current is also an unknown, it cannot be guaranteed that the solution obtained by this method is the optimal solution, and this method has poor stability and large calculation errors.

[0006] Therefore, how to quickly and accurately calculate the induced eddy current has become a technical problem to be solved urgently. Summary of the Invention

[0007] This application provides a method for constructing a calculation model of induced eddy current, a calculation method, and an electronic device to at least solve the technical problem of how to quickly and accurately calculate the induced eddy current in related technologies.

[0008] According to the first aspect of the present application, a method for constructing a calculation model of induced eddy current is provided, which is applicable to the calculation of induced eddy current of the vacuum vessel in a magnetic confinement fusion device. The method for constructing the calculation model of induced eddy current includes: obtaining calibration discharge data of multiple coil calibration guns, where the calibration discharge data includes coil current calibration data and magnetic diagnosis calibration data of a magnetic diagnosis device; determining induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data; using the coil current calibration data as an input quantity and the induced eddy current calibration data as an output quantity, and fitting the response coefficient between the coil current and the induced eddy current by using an induced eddy current calculation model constructed based on a transfer function as the model parameter of the calculation model.

[0009] In one embodiment, using the coil current calibration data as an input quantity and the induced eddy current calibration data as an output quantity, and fitting the response coefficient between the coil current and the induced eddy current by using an induced eddy current calculation model constructed based on a transfer function as the model parameter of the calculation model includes: using a transfer function and a linear fitting function to build a calculation model to be fitted, where the input parameter of the calculation model to be fitted is the coil current and the output parameter is the induced eddy current; substituting multiple sets of coil current calibration data and the induced eddy current calibration data into the calculation model to be fitted to determine the response coefficient.

[0010] In one embodiment, obtaining calibration discharge data of multiple coil calibration guns includes: performing calibration gun discharge on each coil according to a calibration gun discharge mode in which one coil among multiple coils discharges and the other coils among multiple coils are charged; respectively obtaining discharge data of multiple calibration gun discharges of each coil.

[0011] In one embodiment, substituting multiple sets of coil current calibration data and the induced eddy current calibration data into the calculation model to be fitted to determine the response coefficient includes: splicing the coil current calibration data in the discharge data of multiple calibration guns of multiple coils to obtain the input quantity of the calculation model to be fitted; splicing the induced eddy current calibration data corresponding to the coil current calibration data in the discharge data of multiple calibration guns of multiple coils to obtain the output quantity of the calculation model to be fitted; substituting the input quantity and the output quantity into the calculation model to be fitted to fit the response coefficient.

[0012] In one embodiment, determining the induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data includes: calculating the Green's function between the coil, the magnetic diagnosis device, and the vacuum vessel; constructing a magnetic diagnosis balance equation based on the Green's function, where the calibrated magnetic diagnosis data in the magnetic diagnosis balance equation includes a first magnetic diagnosis contribution share corresponding to the coil current calibration data and a second magnetic diagnosis contribution share corresponding to the induced eddy current calibration data; substituting the coil current calibration data and the magnetic diagnosis data into the magnetic diagnosis balance equation to calculate the induced eddy current calibration data.

[0013] In one embodiment, calculating the Green's function between the coil, the magnetic diagnosis device, and the vacuum vessel includes: using the finite element method to mesh the vacuum vessel to obtain a plurality of equivalent mesh regions; calculating the Green's function between each of the equivalent mesh regions, the coil, and the magnetic diagnosis device respectively.

[0014] In one embodiment, obtaining the calibration discharge data of multiple coil calibration shots further includes: screening out the invalid data of the magnetic diagnosis calibration data in the calibration discharge data.

[0015] According to a second aspect, an embodiment of the present application provides an induced eddy current calculation method, including: respectively obtaining the coil current of each coil; inputting the coil current into the calculation model constructed by the calculation model construction method of the induced eddy current according to any one of the first aspects above to obtain a sub-induced eddy current corresponding to each coil; superimposing each sub-induced eddy current to obtain the induced eddy current on the vacuum vessel.

[0016] According to a third aspect, an embodiment of the present application provides an electronic device, including: a processor and a memory; the memory is used to store a computer program; the processor is used to execute the computer program and implement the calculation model construction method of the induced eddy current according to any one of the first aspects above and / or execute the induced eddy current calculation method in the second aspect above when executing the computer program.

[0017] According to a third aspect, an embodiment of the present application provides a computer-readable storage medium, on which a computer instruction is stored, and the computer instruction is used to cause a computer to execute the calculation model construction method of the induced eddy current according to any one of the first aspects above and / or execute the induced eddy current calculation method in the second aspect above.

[0018] The present application has at least the following beneficial effects compared with the prior art:

[0019] In this application, an induction eddy current calculation model is constructed through a transfer function. The calibrated coil current data is used as the input quantity, and the calibrated induction eddy current data is used as the output quantity to fit the response coefficient between the coil current and the induction eddy current it generates on the vacuum wall. Different from the related art where the mutual inductance model needs to be expressed as a differential equation based on physical formulas and the eddy current is calculated by means of the resistance, inductance, mutual inductance coefficient, and mutual inductance formula of the loop, this application directly fits the response coefficient between the coil current and its induction eddy current by means of data processing, without the need to solve numerous and complex physical differential equations. The eddy current induced by the coil current on the vacuum wall can be calculated only through the calibrated response coefficient and the coil current, with low calculation complexity, saving calculation time. By using data-driven and the response coefficient fitted from experimental data, the calculated eddy current is closer to the actual situation, so the credibility will be higher. At the same time, the algorithm of fitting the induction eddy current by the transfer function can not only accurately simulate the eddy current, but also simplify the calculation method and has a fast calculation speed, providing technical support for the subsequent real-time control of the plasma. Accurately fitting the induction eddy current is of great significance for the calculation of the plasma current distribution and the plasma equilibrium inversion. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The drawings described herein are used to provide a further understanding of the present application and form a part of the present application. The schematic embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation of the present application. In the drawings:

[0021] Figure 1 is a schematic flow chart of the method for constructing the calculation model of the induction eddy current provided by the present application;

[0022] Figure 2 is a schematic flow chart of the induction eddy current calculation method provided by the present application;

[0023] Figure 3 is a schematic diagram of a calculated induction eddy current provided by the present application;

[0024] Figure 4 is a schematic diagram of the electronic device provided by the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0025] In order to more clearly illustrate the overall concept of the present application, the following will be described in detail by way of examples in conjunction with the drawings of the specification.

[0026] To enable those skilled in the art to better understand the solution of this application, the technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of this application. Obviously, the described embodiments are only a part of the embodiments of this application, rather than all the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in this application without creative efforts shall fall within the scope of protection of this application.

[0027] It should be noted that the terms "first", "second", etc. in the specification and claims of this application and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily need to be used to describe a specific order or sequence. It should be understood that such data can be interchanged under appropriate circumstances so that the embodiments of this application described here can be implemented in an order other than those illustrated or described here. In addition, the terms "comprising" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products, or devices.

[0028] This application proposes a method for constructing a calculation model of induced eddy current. The calculation model of induced eddy current is used to calculate the induced eddy current induced on the vacuum vessel during the discharge process of a magnetic confinement fusion device. The calculation model constructed in this application uses data-driven, fits the response coefficient through experimental data, and the induced eddy current calculated by it is closer to the actual situation, so the credibility will be higher.

[0029] See Figure 1 As shown, the method for constructing a calculation model of induced eddy current may include the following steps S10 to S30:

[0030] S10, obtaining the calibration discharge data of multiple coil calibration guns, where the calibration discharge data includes coil current calibration data and magnetic diagnosis calibration data of magnetic diagnosis equipment. A magnetic confinement fusion device may include a tokamak, a stellarator, and a magnetic mirror. In this embodiment, a tokamak is taken as an example for illustration.

[0031] It should be noted that a tokamak may include various types of coils, such as poloidal field (PF) coils, toroidal field (TF) coils, center solenoid (CS) coils, and merging-compression (MC) coils. The magnetic diagnosis equipment may include magnetic probes and flux loops.

[0032] In this embodiment, there may be 9 coils including 3 central solenoid coils, 5 poloidal field coils and 1 MC coil, and the magnetic diagnostic device may include 64 magnetic diagnostic devices such as 49 magnetic probes and 15 flux loops.

[0033] In this embodiment, in order to fit the induced eddy current of each coil current on the vacuum vessel wall, during the calibration shot discharge, each of the above-mentioned coils is controlled to discharge alone, while the other coils are charged to resist the eddy current. Specifically, during each calibration shot discharge, each coil is subjected to calibration shot discharge according to the calibration shot discharge method of controlling one coil among the multiple coils to discharge and charging the other coils among the multiple coils until the discharge of all coils is completed, and the discharge data of multiple calibration shot discharges of each coil are respectively obtained.

[0034] In this embodiment, the discharge data may include coil current calibration data and magnetic diagnostic calibration data of the magnetic diagnostic device. Among them, the coil current calibration data can be obtained by using a current detection device to detect the coil current in real time during the coil discharge to obtain the coil current as the coil current calibration data. The magnetic diagnostic calibration data is the magnetic flux of the coil current and the induced eddy current acting on the magnetic diagnostic device together.

[0035] In this embodiment, the discharge can be carried out according to the above-mentioned calibration shot discharge requirements, and the current data of 9 coils including 3 central solenoid coils, 5 poloidal field coils and 1 MC coil are respectively obtained as the coil current calibration data, and the total magnetic flux data of 64 magnetic diagnostic devices including 49 magnetic probes and 15 flux loops are obtained as the magnetic diagnostic calibration data. When collecting the discharge data, n milliseconds can be used as the discharge time window, and the discharge window time is equally spaced and divided into m segments for sampling. Together with the two moments at the beginning and end of the time window, there are a total of m + 1 sampling moments.

[0036] Sampling is carried out based on m + 1 sampling moments, and (m + 1)×9 coil currents and (m + 1)×64 interpolation results of the total magnetic flux of the magnetic diagnostic device are respectively obtained. In this embodiment, the coil current calibration data and the magnetic diagnostic calibration data can be arranged in the form of a matrix, and can also be presented in the form of a sequence. For the convenience of calculation, in this embodiment, the arrangement in the form of a matrix is taken as an example for description.

[0037] In one embodiment, various factors such as the failure of the magnetic diagnostic device itself and the abnormality of the magnetic diagnostic device will cause the collected magnetic diagnostic data to be invalid. Therefore, it is necessary to mark some invalid magnetic diagnostic data in the discharge data of each calibration shot to avoid the interference of invalid magnetic diagnostic devices on the response coefficient of calibration and the calculation of induced eddy current.

[0038] S20. Determine the induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data. In this embodiment, the coil current calibration data and the magnetic diagnosis calibration data are known quantities collected in actuality. Therefore, based on the relationship among the coil current, the induced eddy current, and the magnetic diagnosis data, the induced eddy current can be solved through the known quantities to obtain the induced eddy current calibration data. Exemplarily, the Green's function between the coil, the magnetic diagnosis device, and the vacuum vessel can be calculated. That is, the mutual inductance coefficients between the coil, the magnetic diagnosis device, and the vacuum vessel are calculated respectively. In this embodiment, after the fusion device is built, when the positions of its coil and magnetic diagnosis device and the structure of the fusion device are determined, the mutual inductance coefficients between the coil, the magnetic diagnosis device, and the vacuum vessel are often determined and unchanged. Taking 9 coils such as 3 central solenoid coils, 5 poloidal field coils, and 1 MC coil, and the magnetic diagnosis device including 64 magnetic diagnosis devices such as 49 magnetic probes and 15 flux loops as an example for illustration. Obtain the coordinates of each coil, the installation coordinates of the magnetic diagnosis device, and the coordinates of the vacuum vessel, and calculate the first mutual inductance matrix M1 between the coil and the vacuum vessel, the second mutual inductance matrix M2 between the coil and the magnetic diagnosis device, and the third mutual inductance matrix M3 between the vacuum vessel and the magnetic diagnosis device respectively, to obtain the Green's function between the coil, the magnetic diagnosis device, and the vacuum vessel.

[0039] In one embodiment, due to the complex structure of the vacuum vessel, it is difficult to accurately calculate the eddy current in each region of the vacuum vessel by globally modeling the Tokamak device. Therefore, in this embodiment, the finite element method is used to mesh the vacuum vessel to obtain multiple equivalent mesh regions. When calculating the Green's function, the vacuum vessel is equivalent to multiple equivalent mesh regions. Taking the number of equivalent mesh regions as Q for example, the mutual inductance coefficients between the coil, the magnetic diagnosis device, and the vacuum vessel can be calculated by obtaining the coordinates of each coil, the installation coordinates of the magnetic diagnosis device, and the coordinates of each equivalent mesh region, and calculating the first mutual inductance matrix M1 (with a size of 9×Q) between the coil and multiple equivalent mesh regions, the second mutual inductance matrix M2 (with a size of 9×64) between the coil and the magnetic diagnosis device, and the third mutual inductance matrix M3 (with a size of Q×64) between the equivalent mesh regions and the magnetic diagnosis device respectively, to obtain the Green's function between the coil and the vacuum vessel, between the coil and the magnetic diagnosis device, and between the vacuum vessel and the magnetic diagnosis device.

[0040] After obtaining the Green's function, since the magnetic diagnosis data of the calibration shot is only contributed by two parts: the coil current and the eddy current induced by the coil current on the vacuum vessel, and the coil current is a known quantity, the magnetic diagnosis balance equation can be constructed through the Green's function, and the equation can be solved to obtain the induced eddy current calibration data induced by the coil current calibration data.

[0041] S30. Take the coil current calibration data as the input quantity and the induced eddy current calibration data as the output quantity, and use the induced eddy current calculation model constructed based on the transfer function to fit the response coefficient between the coil current and the induced eddy current as the model parameter of the calculation model. The transfer function is the ratio of the Laplace transform (or z-transform) of the response quantity of a linear system under zero initial conditions to the Laplace transform of the excitation quantity, and can describe the dynamic characteristics of the linear system. Based on this, in this application, an induced eddy current calculation model is constructed through the transfer function, taking the coil current calibration data as the input quantity and the induced eddy current calibration data as the output quantity, and fitting the response coefficient between the coil current and the induced eddy current induced on the vacuum wall. Different from the related art where the mutual inductance model needs to be expressed as a differential equation based on physical formulas, and the eddy current is calculated by means of the resistance, inductance, mutual inductance coefficient and mutual inductance formula of the loop; in this application, the response coefficient between the coil current and its induced eddy current is directly fitted by means of data processing, without the need to solve numerous and complex physical differential equations. The eddy current induced by the coil current on the vacuum wall can be calculated only through the calibrated response coefficient and the coil current, with low calculation complexity, saving calculation time. Using data-driven, the response coefficient fitted by experimental data makes the calculated eddy current closer to the actual situation, so the credibility will be higher. At the same time, the algorithm of fitting the induced eddy current by the transfer function can not only accurately simulate the eddy current, but also simplify the calculation method, with fast calculation speed, providing technical support for the subsequent real-time control of the plasma; accurately fitting the induced eddy current is of great significance for the calculation of the plasma current distribution and the plasma equilibrium inversion.

[0042] In one embodiment, taking the coil current calibration data as the input quantity and the induced eddy current calibration data as the output quantity, and using the induced eddy current calculation model constructed based on the transfer function to fit the response coefficient between the coil current and the induced eddy current as the model parameter of the calculation model includes: using the transfer function and the linear fitting function to build the calculation model to be fitted. In this embodiment, the calculation model to be fitted is based on the transfer function, and the linear relationship between the input quantity and the output quantity is fitted through the fitting function. In this embodiment, the calculation model to be fitted can be expressed by the following formula:

[0043]

[0044] where a0……an, b0……bn are the coefficients of the transfer function, all of which are real numbers; s is the Laplace operator.

[0045] In this embodiment, multiple groups of coil current calibration data and the induced eddy current calibration data are substituted into the to-be-fitted calculation model to determine the response coefficient. Specifically, an initial value of the response coefficient can be specified in advance. A transfer function is constructed with the response coefficient, and the coil current is used as the input of the transfer function to obtain the response signal. The linear fitting function is used to calculate the error between the response signal and the induced eddy current. If the error is large, the fitting function will give a new response coefficient according to the gradient of the error change. Then, the response signal corresponding to the coil current is calculated again using the transfer function until the error is less than the tolerable range. The response coefficient at this time is the response coefficient between the coil current and the induced eddy current. Similarly, the response coefficient between the coil current and the magnetic diagnosis calibration data can also be fitted.

[0046] In one embodiment, to avoid interference from other unknown factors on the data of a single shot, which may affect the accuracy of the data, at least two shots of magnetic diagnosis data with different current values are taken for each power supply to ensure the correctness of the calibration result. When fitting the response coefficient, the coil current calibration data in the discharge data of multiple calibration shots of multiple coils are spliced to obtain the input of the to-be-fitted calculation model; the induced eddy current calibration data corresponding to the coil current calibration data in the discharge data of multiple calibration shots of multiple coils are spliced to obtain the output of the to-be-fitted calculation model; the input and the output are substituted into the to-be-fitted calculation model to fit the response coefficient. By loading the discharge data of the calibration shots of each coil discharged separately, the coil current calibration data and the induced eddy current calibration data induced by the coil current corresponding to each coil discharged separately are spliced to obtain the excitation signal (composed of the spliced coil current calibration data of multiple calibration shots) and the response signal (composed of the spliced induced eddy current calibration data induced by the coil current of multiple calibration shots). The transfer function and the linear fitting function are used to build the fitting algorithm for calibrating the response coefficient (see the specific calibration algorithm in the above embodiment). The spliced excitation signal and response signal are used as the input and the output respectively. After analysis and processing by the fitting algorithm, the response coefficient between the coil current and the induced eddy current of the coil can be obtained.

[0047] The embodiment of the present application also provides an induced eddy current calculation method, as Figure 2 shown. The induced eddy current calculation method includes the following steps S100 to S300:

[0048] S100. Obtain the coil current of each coil respectively. In this embodiment, during the discharge process of the nuclear fusion device, the coil current of each coil is obtained in real time. Exemplarily, the coil currents of 9 coils such as 3 central solenoid coils, 5 poloidal field coils, and 1 MC coil can be obtained respectively.

[0049] S200. Input the coil current into the calculation model constructed by the method for constructing the calculation model of induced eddy current to obtain the sub-induced eddy current corresponding to each coil. This calculation model is pre-constructed, and its construction method is as follows: starting from the actually collected data, construct equations between the coil current, its induced eddy current, and each magnetic diagnostic device, and use the transfer function and fitting function to calibrate the response coefficients between the coil current and the eddy current it contributes, and the response coefficients between the coil current and each magnetic diagnostic signal, so as to obtain the above calculation model.

[0050] S300. Superimpose each sub-induced eddy current to obtain the induced eddy current on the vacuum vessel wall. Construct a calculation model by means of the calibrated response coefficient between the coil current and the response eddy current, and then calculate the eddy current induced by each coil current on the vacuum vessel wall respectively. The eddy currents induced by each coil are superimposed to be the induced eddy current of the coil current on the vacuum vessel wall, as Figure 3 the schematic diagram of the induced eddy current on the vacuum vessel wall calculated by the calculation model constructed based on the response coefficient between the coil current and the induced eddy current shown.

[0051] Furthermore, if the response coefficients between the coil current and the magnetic diagnostic data are marked, the contributions of the coil current and its eddy current to the magnetic diagnostic device can be calculated, simplifying the plasma equilibrium inversion calculation. Without changing the positions of the coils and magnetic diagnostic devices and the device structure, the induced eddy current of each coil in any discharge mode can be calculated through the calibrated response coefficients. Conversely, it is necessary to re-calibrate the gun data and fit new response coefficients.

[0052] According to another aspect of the embodiments of the present application, there is also provided an electronic device for implementing the above method for constructing the calculation model of induced eddy current and / or the method for calculating induced eddy current. This electronic device can be a server, a terminal, or a combination thereof.

[0053] Figure 4 is a structural block diagram of an electronic device according to an embodiment of the present application, as Figure 4 shown, including a processor 401, a communication interface 402, a memory 403, and a communication bus 404. Among them, the processor 401, the communication interface 402, and the memory 403 complete communication with each other through the communication bus 404. Among them,

[0054] The memory 403 is used to store computer programs;

[0055] The processor 401, when executing the computer program stored on the memory 403, realizes the following steps:

[0056] Obtain the calibrated discharge data of multiple coil calibration guns, and the calibrated discharge data includes coil current calibration data and magnetic diagnostic calibration data of magnetic diagnostic devices;

[0057] Determine the induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data;

[0058] Use the coil current calibration data as the input quantity and the induced eddy current calibration data as the output quantity, and use the induced eddy current calculation model constructed based on the transfer function to fit the response coefficient between the coil current and the induced eddy current as the model parameter of the calculation model.

[0059] In one embodiment, when the processor 401 executes the computer program stored on the memory 403, the following steps can also be implemented:

[0060] Obtain the coil current of each coil respectively;

[0061] Input the coil current into the calculation model constructed by the calculation model construction method of the induced eddy current to obtain the sub-induced eddy current corresponding to each coil;

[0062] Superimpose each sub-induced eddy current to obtain the induced eddy current on the vacuum vessel wall.

[0063] The calculation model construction method and the induced eddy current calculation method of the induced eddy current in the above embodiments can be executed in the same electronic device or in different electronic devices. In one embodiment, the above communication bus can be a PCI (Peripheral Component Interconnect) bus, an EISA (Extended Industry Standard Architecture) bus, etc. This communication bus can be divided into an address bus, a data bus, a control bus, etc. For the sake of simplicity of representation, Figure 4 only a thick line is shown in the figure, but it does not mean that there is only one bus or one type of bus.

[0064] The communication interface is used for communication between the above electronic device and other devices.

[0065] The memory can include RAM and can also include non-volatile memory, for example, at least one disk memory. In addition, the memory can also be at least one storage device located far from the aforementioned processor.

[0066] The above-mentioned processor can be a general-purpose processor, which may include but is not limited to: CPU (Central Processing Unit, central processing unit), NP (Network Processor, network processor), etc.; it can also be a DSP (Digital Signal Processing, digital signal processor), ASIC (Application Specific Integrated Circuit, application-specific integrated circuit), FPGA (Field-Programmable Gate Array, field programmable gate array) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0067] Specific examples in this embodiment may refer to the examples described in the above-mentioned embodiment, and will not be repeated here.

[0068] Those of ordinary skill in the art can understand that Figure 4 The structure shown is only schematic, and the device for implementing the method for constructing the calculation model of induced eddy current and / or the method for calculating induced eddy current can be a terminal device. Figure 4 It does not limit the structure of the above-mentioned electronic device. For example, the terminal device may further include more or fewer components (such as a network interface, a display device, etc.) than those shown in Figure 4 or have a different configuration from that shown in Figure 4 Those shown.

[0069] Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above-mentioned embodiments can be completed by instructing the relevant hardware of the terminal device through a program, and the program can be stored in a computer-readable storage medium. The storage medium may include: a flash drive, a ROM, a RAM, a magnetic disk, or an optical disc, etc.

[0070] According to another aspect of the embodiments of the present application, a storage medium is further provided. The above-mentioned storage medium can be used to execute the program code of the method for constructing the calculation model of induced eddy current and / or the method for calculating induced eddy current.

[0071] In one embodiment, the above-mentioned storage medium may be located on at least one of the multiple network devices in the network described above.

[0072] Further, the storage medium is configured to store program code for performing the following steps:

[0073] Obtain the calibration discharge data of multiple coil calibration guns, where the calibration discharge data includes coil current calibration data and magnetic diagnosis calibration data of a magnetic diagnosis device;

[0074] Determine the induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data;

[0075] Use the coil current calibration data as the input quantity and the induced eddy current calibration data as the output quantity, and use the induced eddy current calculation model constructed based on the transfer function to fit the response coefficient between the coil current and the induced eddy current as the model parameter of the calculation model.

[0076] In one embodiment, the storage medium can also be set to store program code for performing the following steps:

[0077] Obtain the coil current of each coil respectively;

[0078] Input the coil current into the calculation model constructed by the calculation model construction method of the induced eddy current to obtain the sub-induced eddy current corresponding to each coil;

[0079] Superimpose each sub-induced eddy current to obtain the induced eddy current on the vacuum vessel wall.

[0080] The program code corresponding to the calculation model construction method and the induced eddy current calculation method in the above embodiments can be stored in the same storage medium or in different storage media.

[0081] The specific examples in this embodiment can refer to the examples described in the above embodiments, and will not be elaborated herein.

[0082] Furthermore, the above storage medium can include but is not limited to: various media such as USB flash drives, ROM, RAM, mobile hard disks, magnetic disks or optical discs that can store program code.

[0083] The serial numbers of the embodiments of the present application above are only for description and do not represent the advantages and disadvantages of the embodiments.

[0084] If the integrated unit in the above embodiments is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in the above computer-readable storage medium. Based on this understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, or all or part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in the storage medium and includes several instructions for causing one or more computer devices (which can be personal computers, servers or network devices, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present application.

[0085] In the above embodiments of the present application, the descriptions of the various embodiments have their own emphases. For the parts not detailed in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0086] In several embodiments provided by this application, it should be understood that the disclosed client can be implemented in other ways. Among them, the device embodiments described above are merely illustrative. For example, the division of the units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of units or modules can be in electrical or other forms.

[0087] The units described as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place or distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution provided in this embodiment.

[0088] In addition, the functional units in each embodiment of this application can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.

[0089] Each embodiment in this specification is described in a progressive manner. The same or similar parts between each embodiment can be referred to each other. The key point of each embodiment is to illustrate the differences from other embodiments.

[0090] The above is only the embodiments of this application and is not used to limit this application. For those skilled in the art, this application can have various changes and modifications. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of this application shall be included within the scope of the claims of this application.

Claims

1. A method for constructing a computational model of induced eddy currents, which is applicable to the computation of induced eddy currents of vacuum walls in magnetic confinement nuclear fusion devices, and is characterized in that: include: Acquiring calibration discharge data of a plurality of coil calibration guns, wherein the calibration discharge data includes coil current calibration data and magnetic diagnosis calibration data of a magnetic diagnosis device; determining induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data; The coil current calibration data is used as input, the induced eddy current calibration data is used as output, and the response coefficient between the coil current and the induced eddy current is fitted by using an induced eddy current calculation model constructed based on a transfer function as a model parameter of the calculation model.

2. The method for constructing a computational model of induced eddy current according to claim 1, characterized in that: The coil current calibration data is used as input, the induced eddy current calibration data is used as output, and the response coefficient between the coil current and the induced eddy current is fitted by using the induced eddy current calculation model constructed based on the transfer function as the model parameter of the calculation model, including: Using the transfer function and the linear fitting function to build a calculation model to be fitted, wherein the input parameter of the calculation model to be fitted is the coil current, and the output parameter is the induced eddy current; Substituting multiple sets of coil current calibration data and the induced eddy current calibration data into the calculation model to be fitted to determine the response coefficient.

3. The method for constructing a computational model of induced eddy current according to claim 2, characterized in that: The obtaining of calibration discharge data of multiple coil calibration guns comprises: Calibration gun discharge is performed on each coil in a calibration gun discharge mode in which one coil among the plurality of coils is controlled to discharge and other coils among the plurality of coils are charged; The discharge data of multiple calibration gun discharges of each coil are obtained respectively.

4. The method for constructing a computational model of induced eddy current according to claim 3, characterized in that: Substituting the multiple sets of coil current calibration data and the induced eddy current calibration data into the calculation model to be fitted to determine the response coefficient comprises: splicing the coil current calibration data in the discharge data of multiple calibration guns of the plurality of coils to obtain the input quantity of the calculation model to be fitted; splicing the induced eddy current calibration data corresponding to the coil current calibration data in the discharge data of multiple calibration guns of the plurality of coils to obtain the output of the calculation model to be fitted; Substitute the input quantity and the output quantity into the calculation model to be fitted to fit the response coefficient.

5. The method for constructing a computational model of induced eddy current according to claim 1, characterized in that: The determining of the induced eddy current calibration data based on the coil current calibration data and the magnetic diagnosis calibration data comprises: calculating a Green's function between the coil, the magnetic diagnostic device and the vacuum wall; Constructing a magnetic diagnostic balance equation based on the Green's function, wherein the calibrated magnetic diagnostic data in the magnetic diagnostic balance equation includes a first magnetic diagnostic contribution share corresponding to the coil current calibration data and a second magnetic diagnostic contribution share corresponding to the induced eddy current calibration data; The coil current calibration data and the magnetic diagnosis data are substituted into the magnetic diagnosis balance equation to calculate and obtain the induced eddy current calibration data.

6. The method for constructing a computational model of induced eddy current according to claim 5, characterized in that: The calculating of the Green's function between the coil, the magnetic diagnostic device and the vacuum wall comprises: Using a finite element method to mesh the vacuum wall, obtaining a plurality of equivalent mesh regions; The Green's functions between the equivalent grid regions, the coils and the magnetic diagnostic device are calculated respectively.

7. The method for constructing a computational model of induced eddy current according to claim 1, characterized in that: The obtaining of calibration discharge data of a plurality of coil calibration guns also includes: Invalid data of the magnetic diagnosis calibration data in the calibration discharge data is filtered out.

8. A method for calculating induced eddy current, characterized in that: include: respectively obtaining the coil current of each coil; Inputting the coil current into the calculation model constructed by the method for constructing a calculation model of induced eddy current according to any one of claims 1 to 7 to obtain a sub-induced eddy current corresponding to each coil; The induced eddy current on the vacuum wall is obtained by superimposing each sub-induced eddy current.

9. An electronic device, characterized in that: include: Processor, memory; The memory is used to store computer programs; The processor is used to execute the computer program and implement the method for constructing a calculation model of the induced eddy current according to any one of claims 1 to 7 and / or execute the method for calculating the induced eddy current according to claim 8 when executing the computer program.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for constructing a calculation model of an induced eddy current according to any one of 1 to 7 and / or execute the method for calculating an induced eddy current according to claim 8.

Images