Method and device for calculating plasma current flat-top time of nuclear fusion device

Abstract

This application belongs to the field of nuclear fusion technology, specifically relating to a method and apparatus for calculating the plateau time of plasma current in a nuclear fusion device. The method includes: acquiring waveform data of plasma current magnitude and time; dividing the data points in the waveform data into multiple data groups based on the magnitude of the plasma current corresponding to each data point, obtaining the number of data points in each data group; different data groups correspond to different magnitudes of plasma current; and obtaining the plateau time of the plasma current based on the time corresponding to the second and third largest data groups. This method solves the problems of accurate algorithm implementation for obtaining the plateau time of plasma current in controlled nuclear fusion in tokamak and the scalability for subsequent algorithm modifications.

Description

Technical Field

[0001] This application belongs to the field of nuclear fusion technology, specifically relating to a method and apparatus for calculating the plasma current flat-top time of a nuclear fusion device. Background Technology

[0002] A tokamak is a toroidal container that uses magnetic confinement to achieve controlled nuclear fusion. It consists of a central toroidal vacuum chamber surrounded by coils. When energized, a powerful helical magnetic field is generated inside the tokamak, heating the plasma to extremely high temperatures to achieve nuclear fusion.

[0003] In a tokamak device, the current variation of an ohmic coil provides the plasma current for generating, establishing, and maintaining it; its physical waveform is as follows: Figure 1 As shown. Whether the plasma current is stable is one of the most important experimental conclusions of tokamak nuclear fusion experiments. Its quantitative index is the plateau time used to describe the sustained stability of the plasma current intensity, i.e. Figure 1 The area shown is within the dashed box. From a macroscopic perspective, it's clear that the waveform here approximates a square wave, indicating that the current is constant. The time domain between the start and end of the square wave is the flat-top time.

[0004] The current mainstream method for obtaining the flat-top time involves analyzing a two-dimensional array of generated waveforms. Starting from a point where the Y-axis is not zero, the actual Y-value is compared to the maximum value. If the Y-axis value exceeds n consecutive points but is less than a set difference, the first point is taken as the starting point of the flat-top time; if it exceeds n consecutive points but is less than the set difference, the last point is taken as the ending point. This calculation method has limitations when the macroscopic appearance appears flat, but there are microscopic fluctuations in the plasma current, such as… Figure 2 As shown, the accurate time of the flattening cannot be obtained. Summary of the Invention

[0005] The purpose of this application is to provide a method and apparatus for calculating the flat-top time of plasma current in a nuclear fusion device, which solves the problem in the prior art where the plasma current waveform is macroscopically flat-topped but microscopically fluctuating, making it difficult to calculate experimental parameters.

[0006] The technical solution to achieve the purpose of this application is as follows:

[0007] The first aspect of this application provides a method for calculating the plateau time of plasma current in a nuclear fusion device, the method comprising:

[0008] Acquire waveform data of plasma current magnitude and time;

[0009] Based on the magnitude of the plasma current corresponding to the data points in the waveform data, the data points are divided into multiple data groups to obtain the number of data points in each data group; the magnitude of the plasma current corresponding to different data groups is different.

[0010] The plateau time of the plasma current is obtained based on the time corresponding to the second and third largest data groups.

[0011] Optionally, the step of dividing the data points into multiple data groups based on the magnitude of the plasma current corresponding to the data points in the waveform data, and obtaining the number of data points in each data group, specifically includes:

[0012] Determine the maximum value of the plasma current in the waveform data. max and minimum value data min ;

[0013] Based on the number of data groups, data max and data min This gives the data span for each data group.

[0014] The number of data points in each data group is obtained based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data.

[0015] Optionally, based on the number of data groups, data max and data min This yields the data span for each data group, specifically including:

[0016] The data span gap for each data group is obtained according to the following formula (1);

[0017] gap=data max / bins (1)

[0018] Here, bins represents the number of data sets.

[0019] Optionally, the step of obtaining the number of data points in each data group based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data specifically includes:

[0020] Based on the data span gap and the magnitude of the plasma current corresponding to the data point. j Determine the data group to which the data point belongs, and obtain the number of data points in each data group.

[0021] Optionally, obtaining the plateau time of the plasma current based on the time corresponding to the second and third largest data groups specifically includes:

[0022] Determine the earliest and latest times for each data point in the second and third largest data sets;

[0023] The earliest time and the latest time are used as the two boundaries of the flat-top time.

[0024] A second aspect of this application provides a plasma current flat-top time calculation device for a nuclear fusion device, the device comprising:

[0025] The waveform acquisition module is used to acquire waveform data of plasma current magnitude and time.

[0026] The quantity determination module is used to divide the data points into multiple data groups according to the magnitude of the plasma current corresponding to the data points in the waveform data, and to obtain the number of data points in each data group; the magnitude of the plasma current corresponding to different data groups is different;

[0027] The time determination module is used to obtain the plateau time of the plasma current based on the time corresponding to the second and third largest data groups.

[0028] Optionally, the quantity determination module is specifically used for:

[0029] Determine the maximum value of the plasma current in the waveform data. max and minimum value data min ;

[0030] Based on the number of data groups, data max and data min This gives the data span for each data group.

[0031] The number of data points in each data group is obtained based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data.

[0032] Optionally, the quantity determination module is specifically used for:

[0033] The data span gap for each data group is obtained according to the following formula (1);

[0034] gap=data max / bins (1)

[0035] Here, bins represents the number of data sets.

[0036] Optionally, the quantity determination module is specifically used for:

[0037] Based on the data span gap and the magnitude of the plasma current corresponding to the data point. j Determine the data group to which the data point belongs, and obtain the number of data points in each data group.

[0038] Optionally, the time determination module is specifically used for:

[0039] Determine the earliest and latest times for each data point in the second and third largest data sets;

[0040] The earliest time and the latest time are used as the two boundaries of the flat-top time.

[0041] The beneficial technical effects of this application are as follows:

[0042] This application provides a method and apparatus for calculating the plateau time of plasma current in a nuclear fusion device. It separates the algorithm from the calculation program and provides an accurate algorithm and calculation tools. It solves the problems of correctly obtaining the plateau time of plasma current in controlled nuclear fusion of tokamak, the scalability of subsequent algorithm changes, and provides experimental tools for selecting algorithms for the plateau time of plasma in controlled nuclear fusion of tokamak, using algorithms for data budgeting, and verifying and modifying the accuracy of historical experimental data. Attached Figure Description

[0043] Figure 1 This is a waveform diagram of a plasma current.

[0044] Figure 2 This is a waveform diagram of another type of plasma current;

[0045] Figure 3 A flowchart illustrating a method for calculating the plasma current flat-top time of a nuclear fusion device, provided in an embodiment of this application;

[0046] Figure 4 A schematic diagram showing the distribution of the number of data points in a data group in a method for calculating the plasma current flat-top time of a nuclear fusion device provided in an embodiment of this application;

[0047] Figure 5 This is a schematic diagram of the plateau time in a method for calculating the plateau time of plasma current in a nuclear fusion device provided in an embodiment of this application. Detailed Implementation

[0048] To enable those skilled in the art to better understand this application, the technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the embodiments described below are only a part of the embodiments of this application, and not all of them. Based on the embodiments described in this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0049] See Figure 3 The figure is a flowchart illustrating a method for calculating the plasma current flat-top time of a nuclear fusion device according to an embodiment of this application.

[0050] This application provides a method for calculating the plasma current flat-top time of a nuclear fusion device, including:

[0051] Step S101: Obtain waveform data of plasma current magnitude and time.

[0052] Step S102: Based on the magnitude of the plasma current corresponding to the data points in the waveform data, divide the data points into multiple data groups to obtain the number of data points in each data group.

[0053] In the embodiments of this application, the magnitudes of plasma currents corresponding to different data groups are different.

[0054] Step S103: Obtain the plateau time of the plasma current based on the time corresponding to the data groups with the second and third largest number of data points.

[0055] The calculation of the plateau time of plasma current in a tokamak fusion device is mainly used to calculate the duration of the macroscopically stable discharge portion of the plasma current generated in controlled nuclear fusion in a tokamak. Due to the inherent physical complexity of controlled nuclear fusion, the plasma current exhibits macroscopic stability but microscopic fluctuations, necessitating an effective algorithm and reliable computational tools to address the computational challenges of calculating its plateau time.

[0056] In this embodiment, the data points are divided into multiple data groups based on the magnitude of the plasma current corresponding to the data points in the waveform data. The peak time of the plasma current is obtained based on the time corresponding to the second and third largest data groups, effectively avoiding the interference of micro-fluctuations on the peak time calculation.

[0057] In practical implementation, the computational tool is used to apply the flat-top time algorithm to actual experimental data. Because large-scale controlled nuclear fusion experiments require continuous testing, both historical and real-time data exist. The computational tool is used to perform real-time calculations on the real-time data and to estimate and compare historical data, thereby obtaining more accurate experimental conclusions. Its functionalities include:

[0058] (1) Algorithm selection: The calculation tool supports the selection of different algorithms to calculate the flat-top time.

[0059] The basic algorithm library is a collection of algorithms primarily used to perform single-value selection or basic operations on experimental data during computation, enabling their application in subsequent complex algorithms. The system supports major mathematical calculation algorithms, as shown in Table 1.

[0060] Table 1

[0061]

[0062]

[0063] (2) Data budget: This means calling the algorithm to pre-calculate the flat-top time for a specified single experiment, storing the calculation results in memory and displaying them, without directly filling in the result data.

[0064] (3) Data recalculation: After estimating the results of historical or current experiments and determining that the calculation results are valuable and real data, the algorithm can be called to calculate the actual experimental data, and then the recalculated experimental data can be stored. Data recalculation supports segmented recalculation of historical experiments (i.e., different experiments using different algorithms or different parameters).

[0065] By separating the algorithm from the computational program and providing accurate algorithms and computational tools, this paper solves the problems of algorithmic implementation for correctly obtaining the plateau time of the plasma current in controlled nuclear fusion of tokamak, scalability for subsequent algorithm changes, and provides experimental tools for algorithm selection for the plateau time of controlled nuclear fusion of tokamak plasma, data budgeting using the algorithm, and accuracy verification and modification of historical experimental data.

[0066] In some possible implementations of this application, step S102 may specifically include:

[0067] Determine the maximum value of the plasma current in the waveform data. max and minimum value data min ;

[0068] Based on the number of data groups, data max and data min This gives the data span for each data group.

[0069] The number of data points in each data group is obtained based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data.

[0070] In one example, based on the number of data groups, data maxand data min This yields the data span for each data group, specifically including:

[0071] The data span gap for each data group is obtained according to the following formula (1);

[0072] gap=data max / bins (1)

[0073] Here, bins represents the number of data sets.

[0074] In another example, determining the number of data points in each data group based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data specifically includes:

[0075] Based on the data span gap and the magnitude of the plasma current corresponding to the data point. j Determine the data group to which the data point belongs, and obtain the number of data points in each data group.

[0076] Figure 4 An example is shown illustrating the distribution of the number of data points in a dataset. From Figure 4 As can be seen, the data group with the largest number of occurrences actually corresponds to data points where the current is zero, while the data groups with the second and third largest number of occurrences correspond to data points where the peak is flat.

[0077] Therefore, in some possible implementations of the embodiments of this application, step S103 specifically includes:

[0078] Determine the earliest and latest times for each data point in the second and third largest data sets;

[0079] The earliest time and the latest time are used as the two boundaries of the flat-top time.

[0080] Figure 5 An example diagram illustrates the plateau time obtained using a method for calculating the plateau time of plasma current in a nuclear fusion device provided in an embodiment of this application.

[0081] Based on the above-described embodiment of a method for calculating the plasma current flat-top time of a nuclear fusion device, this application also provides a device for calculating the plasma current flat-top time of a nuclear fusion device.

[0082] This application provides a plasma current flat-top time calculation device for a nuclear fusion device, comprising:

[0083] The waveform acquisition module is used to acquire waveform data of plasma current magnitude and time.

[0084] The quantity determination module is used to divide the data points into multiple data groups according to the magnitude of the plasma current corresponding to the data points in the waveform data, and to obtain the number of data points in each data group; the magnitude of the plasma current corresponding to different data groups is different;

[0085] The time determination module is used to obtain the plateau time of the plasma current based on the time corresponding to the second and third largest data groups.

[0086] In some possible implementations of the embodiments of this application, the quantity determination module is specifically used for:

[0087] Determine the maximum value of the plasma current in the waveform data. max and minimum value data min ;

[0088] Based on the number of data groups, data max and data min This gives the data span for each data group.

[0089] The number of data points in each data group is obtained based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data.

[0090] In one example, the quantity determination module is specifically used for:

[0091] The data span gap for each data group is obtained according to the following formula (1);

[0092] gap=data max / bins (1)

[0093] Here, bins represents the number of data sets.

[0094] In another example, the quantity determination module is specifically used for:

[0095] Based on the data span gap and the magnitude of the plasma current corresponding to the data point. j Determine the data group to which the data point belongs, and obtain the number of data points in each data group.

[0096] In some possible implementations of the embodiments of this application, the time determination module is specifically used for:

[0097] Determine the earliest and latest times for each data point in the second and third largest data sets;

[0098] The earliest time and the latest time are used as the two boundaries of the flat-top time.

[0099] The present application has been described in detail above with reference to the accompanying drawings and embodiments. However, the present application is not limited to the above embodiments. Within the scope of knowledge possessed by those skilled in the art, various changes can be made without departing from the spirit of the present application. All content not described in detail in this application can be derived from existing technology.

Claims

1. A method for calculating the plateau time of plasma current in a nuclear fusion device, characterized in that, The method includes: Acquire waveform data of plasma current magnitude and time; Based on the magnitude of the plasma current corresponding to the data points in the waveform data, the data points are divided into multiple data groups to obtain the number of data points in each data group; the magnitude of the plasma current corresponding to different data groups is different. The plateau time of the plasma current is obtained based on the time corresponding to the second and third largest data groups. The step of dividing the data points into multiple data groups based on the magnitude of the plasma current corresponding to the data points in the waveform data, and obtaining the number of data points in each data group, specifically includes: Determine the maximum value of the plasma current in the waveform data. and minimum value ; Based on the number of data sets, and This yields the data span for each data group. The number of data points in each data group is obtained based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data. Based on the number of data sets, and This yields the data span for each data group, specifically including: The data span gap for each data group is obtained using the following formula; in, The number of data sets; The step of determining the number of data points in each data group based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data specifically includes: Based on the data span gap and the magnitude of the plasma current corresponding to the data point Determine the data group to which the data point belongs, and obtain the number of data points in each data group; The step of obtaining the plateau time of the plasma current based on the time corresponding to the second and third largest data groups in terms of the number of data points specifically includes: Determine the earliest and latest times for each data point in the second and third largest data sets; The earliest time and the latest time are used as the two boundaries of the flat-top time.

2. A plasma current flat-top time calculation device for a nuclear fusion device, characterized in that, The device includes: The waveform acquisition module is used to acquire waveform data of plasma current magnitude and time. The quantity determination module is used to divide the data points into multiple data groups according to the magnitude of the plasma current corresponding to the data points in the waveform data, and obtain the number of data points in each data group; the magnitude of the plasma current corresponding to different data groups is different; The time determination module is used to obtain the plateau time of the plasma current based on the time corresponding to the second and third largest data groups in terms of the number of data points. The quantity determination module is specifically used for: Determine the maximum value of the plasma current in the waveform data. and minimum value ; Based on the number of data sets, and This yields the data span for each data group. The number of data points in each data group is obtained based on the data span and the magnitude of the plasma current corresponding to the data points in the waveform data. The quantity determination module is specifically used for: The data span gap for each data group is obtained using the following formula; in, The number of data sets; The quantity determination module is specifically used for: Based on the data span gap and the magnitude of the plasma current corresponding to the data point Determine the data group to which the data point belongs, and obtain the number of data points in each data group; The time determination module is specifically used for: Determine the earliest and latest times for each data point in the second and third largest data sets; The earliest time and the latest time are used as the two boundaries of the flat-top time.

Images