A method and system for judging dynamic super-uniformity

Abstract

The present invention discloses a method and system for determining dynamic hyperuniformity. The method comprises: S1, using high-speed video to capture the trajectory of tracer particles in a reactor, and obtaining the two-dimensional image coordinates of the tracer particle trajectory based on a visual tracking algorithm; S2, extracting two sets of two-dimensional image coordinates of the same tracer particle trajectory obtained by two high-speed cameras, performing three-dimensional reconstruction on the two sets of two-dimensional image coordinates, and then performing coordinate normalization processing; S3, using the tracer particle trajectory and the normalized three-dimensional coordinates to determine the mixing uniformity of the mixing system using a dynamic hyperuniformity method. The present invention helps improve the uniformity of the mixing process of the mixing system in industrial production, enhance the mixing effect and mass transfer efficiency of the reactor, optimize the production process, reduce production costs, and improve product quality.

Description

Technical Field

[0001] The present invention relates to the technical field of metallurgy and chemical engineering, and in particular to a method and system for determining dynamic hyperuniformity. Background Art

[0002] In industrial production, the mixing process of a mixing system is a common operation. The determination of mixing uniformity is a key issue. The uniformity of the mixing process directly affects the quality of the product and production efficiency. Mixing uniformity plays a vital role in overall performance.

[0003] Mixing is a common process in the chemical and metallurgical industries. Current methods, mostly based on empirical rules or traditional statistical methods, have limitations and cannot fully and accurately reflect mixing uniformity. Traditional methods often overlook the complexity and dynamics of mixing systems, making it difficult to fully and accurately assess mixing uniformity. Therefore, a new approach is needed to address this problem.

[0004] In a mixing system, stirring is the process of mixing various material components under the influence of external forces, achieving a uniform distribution of particles within any given volume. Stirring is a crucial step in ensuring the quality of mixed materials and improving their performance. Mixing efficiency measures the uniformity of the dispersion and blending of different components. Mixing time is a key parameter for evaluating the uniformity of mixing in a mixing system. Currently, technologies for monitoring solid-liquid mixing effects are categorized as invasive and non-invasive. Invasive testing includes measuring electrical and thermal conductivity using conductivity and thermal conductivity probes. Non-invasive testing includes image analysis and processing, mathematical quantification, ultrasonic attenuation, and electrical capacitance tomography. Mathematical quantification is a common method for quantitatively evaluating engineering problems, and many researchers have used mathematical methods to study the efficiency of multiphase mixing. The Q method (potential quality measurement method) is widely used to measure the uniformity of a set of points and can tolerate occasional close or even overlapping points. However, this research is limited to the uniformity of a single set of points, and its spatial distribution focuses on the relationships between particles, lacking consideration of the overall spatial distribution. To reduce the degree of overall misjudgment, the torque balance method divides the entire system into blocks and takes the average. This division can reduce the degree of misjudgment, but it cannot fundamentally eliminate the effects of local and overall inhomogeneity. Xu et al. proposed L2-star discrepancy (CD) and wrap-around L2-star discrepancy (WD) to evaluate the mixing behavior of 3D regions through numerical simulation, but no experimental comparison has been conducted. During mixing operations, the evaluation of mixing effects is crucial. While these methods can be used to characterize the degree of fluid mixing, their use is somewhat restrictive with respect to the working fluid. Summary of the Invention

[0005] In view of the above problems, the present invention provides a method for determining dynamic hyperuniformity, the method comprising:

[0006] S1, using high-speed camera to extract the two-dimensional image coordinates of the tracer particle trajectory in the reactor;

[0007] S2, extracting two sets of two-dimensional image coordinates of the same tracer particle trajectory acquired by two high-speed cameras, performing three-dimensional reconstruction on the two sets of two-dimensional image coordinates and then performing coordinate normalization processing;

[0008] S3. Based on the motion trajectory of the tracer particles and the normalized three-dimensional image coordinates, the dynamic hyperuniform method is used to determine the uniformity of the mixing system.

[0009] Optionally, in S1, extracting two-dimensional image coordinates of the tracer particle trajectory in the reactor using high-speed photography specifically includes:

[0010] The tracer particle trajectories in the reactor are captured using high-speed cameras, where the high-speed cameras are positioned perpendicular to each other and on the same horizontal plane;

[0011] The visual tracking algorithm is used to locate the coordinates of the tracer particle motion trajectory on the captured image to obtain the two-dimensional image coordinates of the tracer particle motion trajectory.

[0012] Optionally, in S2, the coordinate normalization processing after the three-dimensional reconstruction of the two sets of two-dimensional image coordinates specifically includes:

[0013] Assume that the two sets of extracted two-dimensional image coordinates are (x1, y1) and (x2, y2);

[0014] The two sets of 2D image coordinates are reconstructed into 3D coordinates, and the 3D coordinates are denoted as (X, Y, Z), where X = x1, Y = x2, and Z = (y1 + y2) / 2;

[0015] The obtained three-dimensional coordinates are normalized.

[0016] Optionally, in S3, judging the mixing uniformity of the mixing system using a dynamic hyperuniform method based on the tracer particle motion trajectory and the normalized three-dimensional coordinates specifically includes:

[0017] Based on the motion trajectory of the tracer particles and the normalized three-dimensional image coordinates, the dynamic hyperuniform state time I when the dynamic hyperuniform factor λ approaches the critical value 1 is calculated. The dynamic hyperuniform state time is used to characterize the mixing uniformity of the mixing system.

[0018] Optionally, the calculation method of the dynamic hyperuniform factor λ includes:

[0019] Divide the three-dimensional image into equal-sized regions using the three-dimensional image coordinates;

[0020] The number and variance of the tracer balls in each region were calculated using MATLAB, and the dynamic hyperuniform factor λ was obtained based on the number and variance of the tracer balls.

[0021] Alternatively, when λ ≥ 1, the mixture is in a hyperuniform state; λ ≤ 0, the mixture is in a disordered state; and 0 < λ < 1, the mixture is in a dynamic hyperuniform state. Specifically, when 0 < λ < 1, the magnitude of λ can qualitatively describe the strength of the dynamic hyperuniform state: the closer λ is to 0, the stronger the disordered state; the closer λ is to 1, the stronger the dynamic hyperuniform state.

[0022] The present invention also discloses a dynamic hyperuniformity judgment system, which includes: a coordinate extraction module, a coordinate preprocessing module and a uniformity judgment module;

[0023] The coordinate extraction module is used to extract the two-dimensional image coordinates of the tracer particle trajectory in the reactor using high-speed photography;

[0024] The coordinate preprocessing module is used to extract two sets of two-dimensional image coordinates of the same tracer particle trajectory obtained by two high-speed cameras, and perform coordinate normalization processing after three-dimensional reconstruction of the two sets of two-dimensional image coordinates;

[0025] The uniformity judgment module is used to judge the mixing uniformity of the mixing system using a dynamic hyperuniformity method based on the motion trajectory of the tracer particles and the three-dimensional coordinates after normalization.

[0026] Optionally, the coordinate extraction module further includes a tracer particle track image acquisition submodule and a track coordinate extraction submodule;

[0027] The tracer particle track image acquisition submodule is used to capture the tracer particle track in the reactor using high-speed cameras, with the positions of the high-speed cameras being perpendicular to each other and on the same horizontal plane;

[0028] The trajectory coordinate extraction submodule is used to use a visual tracking algorithm to locate the coordinates of the tracer particle motion trajectory on the captured image to obtain the two-dimensional image coordinates of the tracer particle motion trajectory.

[0029] Optionally, the coordinate preprocessing module further includes: a dimension conversion submodule and a coordinate normalization submodule;

[0030] The dimension conversion submodule is used to set the two sets of extracted two-dimensional image coordinates as (x1, y1) and (x2, y2) respectively; reconstruct the two sets of two-dimensional image coordinates into three dimensions, and the three-dimensional coordinates are denoted as (X, Y, Z), where X = x1, Y = x2, and Z = (y1 + y2) / 2;

[0031] The coordinate normalization submodule is used to perform normalization processing on the obtained three-dimensional coordinates.

[0032] Optionally, the uniformity judgment module calculates the dynamic hyperuniform state time I when the dynamic hyperuniform factor λ approaches the critical value 1 based on the tracer particle motion trajectory and the normalized three-dimensional image coordinates, and uses the dynamic hyperuniform state time to characterize the mixing uniformity of the mixing system.

[0033] Compared with the prior art, the present invention has the following beneficial effects:

[0034] (1) The present invention is supported by the dual-camera positioning method and the dynamic super-uniformity theory, and has high feasibility and simple operation. (2) The present invention uses the dynamic super-uniformity method to determine whether the tracer ball trajectory reaches the dynamic super-uniform state, and then evaluates the mixing uniformity of the tracer ball trajectory. This method fully considers the spatiotemporal distribution and complexity of the mixing system, and evaluates the uniformity of the mixing system by comprehensively analyzing the dynamic motion state and spatial distribution of each phase component in the mixing system. (3) The dynamic super-uniformity method proposed in the present invention evaluates the uniformity of the mixing system, and this method can quantitatively evaluate the mixing characteristics of the three-dimensional stirring system. (4) The dynamic super-uniformity method proposed in the present invention has great research significance in the field of engineering applications, and provides a basis for studying the mixing uniformity of the mixing system, enhancing the stirring effect, and optimizing the three-dimensional stirring equipment. It helps to improve the uniformity of the mixing process of the mixing system in industrial production, improve the mixing effect and mass transfer efficiency of the reactor, optimize the production process, reduce production costs, and improve product quality. (5) The dynamic superuniformity of the present invention combines the ever-changing characteristics of particles in a uniform fluid by analogy with the density fluctuations of a crystal structure, and based on the ergodic theory of mixing, strictly considers the spatial distribution of fluid particles, which can effectively judge the mixing uniformity of fluid mixing. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the technical solution of the present invention, the following briefly introduces the drawings required for use in the embodiments. 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 any creative work.

[0036] Figure 1 is a flow chart of the method of the present invention;

[0037] Figure 2 This is a schematic diagram of the dynamic super-uniform method of the present invention to characterize the uniformity of the mixing system;

[0038] Figure 3 This is a schematic diagram of the principle of the dynamic super-uniform method of the present invention. DETAILED DESCRIPTION

[0039] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0040] Example 1

[0041] The present invention provides a method for determining dynamic superuniformity, such as Figure 1-Figure 2 As shown, the method includes:

[0042] S1. Use high-speed camera to extract the two-dimensional image coordinates of the tracer particle trajectory in the reactor.

[0043] The present invention takes the mixing uniformity of the mixing system in the industrial production process as the research object; then, since the high-speed camera can quickly and effectively capture instantaneous images in the mixing system, the dual cameras capture one frame per second and the shooting area is selected as the fully developed flow area, and the visual tracking algorithm records the coordinate position of the instantaneous image captured by the high-speed camera, the two-dimensional coordinates of the tracer particles can be recorded in real time, and then the two-dimensional coordinates are recorded and saved in a computer, which is convenient for the subsequent reconstruction and calculation of the three-dimensional coordinates.

[0044] S2. Extract two sets of 2D image coordinates of the same tracer particle trajectory captured by two high-speed cameras. These two sets of 2D image coordinates are then reconstructed into three dimensions and then normalized. The specific contents of normalization are as follows: The acquired 3D coordinate data requires normalization preprocessing to enable comparison and analysis of data with different characteristics on the same scale, ensuring that the data range is between [0, 1] and eliminating undesirable effects caused by singular sample data. Dual cameras capture and record the tracer particle trajectory in real time, forming two sets of 2D coordinates (x1, y1) and (x2, y2) for the tracer particle trajectory. These two sets of 2D coordinates are then reconstructed into three dimensions (X, Y, Z). In the 3D reconstructed coordinates, X = x1 and Y = x2. Due to experimental error, this results in an average relative deviation of y1 and y2 of 0.6%. Therefore, the Z axis is selected as (y1 + y2) / 2 to ensure the accuracy of the experimental results. Finally, the reconstructed three-dimensional coordinates are normalized to eliminate the undesirable effects caused by the singular sample data.

[0045] S3. Based on the tracer particle trajectory and the normalized three-dimensional image coordinates, the dynamic hyperuniform method is used to determine the uniformity of the mixing system. In this embodiment, the three-dimensional scatter distribution of the tracer ball trajectory is used to determine the uniformity of the mixing system.

[0046] Based on the spatial distribution of the tracer particles and the normalized three-dimensional coordinates, the time it takes for the tracer ball trajectory to reach the dynamic super-uniform state is obtained, and the speed of the dynamic super-uniform state is compared to judge the uniformity of the mixing achieved by the tracer ball trajectory. Figure 3 In this embodiment, the spatial distribution is a three-dimensional stirring area, and the distribution of the tracer balls in the stirring area.

[0047] The dynamic hyperuniform state time I when the dynamic hyperuniform factor λ approaches the critical value 1 is calculated, and the dynamic hyperuniform state time is used to characterize the mixing uniformity of the mixing system.

[0048] Dynamic hyperuniformity is similar to the local density fluctuation patterns associated with ordinary lattices in arbitrary spatial dimensions. On large scales, it exhibits hyperuniformity similar to that of crystals, but on small scales, it can continuously aggregate and dissipate, exhibiting large fluctuations. Density fluctuations exhibit hyperuniformity similar to that of crystals on large scales, but on small scales, it can continuously aggregate and dissipate, exhibiting large fluctuations. This characteristic can be correlated with the trajectory motion of particles in a fluid (a solution in a stirred tank). The uniformity of the trajectory points is calculated using the dynamic hyperuniformity formula described above, and the uniformity is determined by λ.

[0049] Based on the spatial distribution of tracer particles and the normalized three-dimensional coordinates, the dynamic hyperuniform method is used to obtain the time it takes for the mixing system to reach the dynamic hyperuniform state. The time of the dynamic hyperuniform state can characterize the mixing uniformity of the system mixing, and the mixing effect of the mixing system can be judged by comparing the speed of the dynamic hyperuniform state.

[0050] The time it takes for the dynamic hyperuniformity method to reach a uniform state can be used to characterize the mixing time required for each set of operating conditions to reach a uniform state. Therefore, using a dynamic hyperuniformity method for evaluating and quantifying three-dimensional mixing uniformity can effectively assess the mixing efficiency of each operating condition.

[0051] Reference Figure 3 As shown in the figure, the principle of dynamic hyperuniform theory combines the ever-changing characteristics of particles in a uniform fluid by analogy with the density fluctuations of a crystal structure. At the same time, the idea of ergodicity in mathematics is used, that is, when a particle appears at different locations in a fluid domain over a period of time, then the particle has traversed the fluid domain and can be considered to have reached a mixed state. During the stirring process, the three-dimensional coordinates of the tracer ball are reconstructed by a dual-camera positioning method, and the three-dimensional mixed image is divided into eight equal-sized regions N. n , (n=1,2……,8), the region shape is a cuboid, use Matlab to calculate the number of tracer balls in each region, and get 8 groups of tracer balls L N , (N=1,2......,8), find the variance of β1 in 8 regions, and use β 2Indicated by variance β 2 The correlation β 2 ~N d-λ To reflect the mixing uniformity, we calculated λ, where λ represents the dynamic hyperuniform state factor, N = 8, and d represents the dimension (in three dimensions, d = 3). Analyzing the λ value, λ ≥ 1 indicates a hyperuniform mixing state; λ ≤ 0 indicates a disordered mixing state; and 0 < λ < 1 indicates a dynamically hyperuniform mixing state. Specifically, when 0 < λ < 1, the magnitude of λ can qualitatively describe the strength of the dynamic hyperuniform state: the closer λ is to 0, the stronger the disorder; the closer λ is to 1, the stronger the dynamic hyperuniform state.

[0052] Example 2

[0053] The present invention also provides a dynamic hyperuniformity determination system, comprising:

[0054] The coordinate extraction module is used to extract the two-dimensional image coordinates of the tracer particle trajectory in the reactor using high-speed photography.

[0055] The present invention takes the mixing uniformity of the mixing system in the industrial production process as the research object; then, since the high-speed camera can quickly and effectively capture instantaneous images in the mixing system, the dual cameras capture one frame per second and the shooting area is selected as the fully developed flow area, and the visual tracking algorithm records the coordinate position of the instantaneous image captured by the high-speed camera, the two-dimensional coordinates of the tracer particles can be recorded in real time, and then the two-dimensional coordinates are recorded and saved in a computer, which is convenient for the subsequent reconstruction and calculation of the three-dimensional coordinates.

[0056] The coordinate preprocessing module is used to extract two sets of two-dimensional image coordinates of the same tracer particle trajectory obtained by two high-speed cameras, and perform coordinate normalization processing after three-dimensional reconstruction of the two sets of two-dimensional image coordinates;

[0057] Dual cameras capture and record the tracer particle trajectories in real time, generating two sets of two-dimensional coordinates (x1, y1) and (x2, y2). These two sets of tracer particle trajectories are then reconstructed into three dimensions (X, Y, Z). In these reconstructed coordinates, X = x1 and Y = x2. Due to experimental error, this results in an average relative deviation of y1 and y2 of 0.6%. Therefore, the Z axis is set to (y1 + y2) / 2 to ensure the accuracy of the experimental results. Finally, the reconstructed three-dimensional coordinates are normalized to eliminate undesirable effects caused by singular sample data.

[0058] The uniformity judgment module is used to judge the mixing uniformity of the mixing system using the dynamic hyperuniform method based on the motion trajectory of the tracer particles and the normalized three-dimensional coordinates.

[0059] Based on the spatial distribution of the tracer particles and the normalized three-dimensional coordinates, the time it takes for the tracer ball trajectory to reach the dynamic super-uniform state is obtained, and the speed of the dynamic super-uniform state is compared to judge the uniformity of the mixing achieved by the tracer ball trajectory. Figure 3 shown.

[0060] The dynamic hyperuniform state time I when the dynamic hyperuniform factor λ approaches the critical value 1 is calculated, and the dynamic hyperuniform state time is used to characterize the mixing uniformity of the mixing system.

[0061] Based on the spatial distribution of tracer particles and the normalized three-dimensional coordinates, the dynamic hyperuniform method is used to obtain the time it takes for the mixing system to reach the dynamic hyperuniform state. The time of the dynamic hyperuniform state can characterize the mixing uniformity of the system mixing, and the mixing effect of the mixing system can be judged by comparing the speed of the dynamic hyperuniform state.

[0062] The time it takes for the dynamic hyperuniformity method to reach a uniform state can be used to characterize the mixing time required for each set of operating conditions to reach a uniform state. Therefore, using a dynamic hyperuniformity method for evaluating and quantifying three-dimensional mixing uniformity can effectively assess the mixing efficiency of each operating condition.

[0063] Reference Figure 3 As shown in the figure, the principle of dynamic hyperuniformity theory combines the ever-changing characteristics of particles in a uniform fluid by analogy with the density fluctuations of a crystal structure. At the same time, the idea of ergodicity in mathematics is used, that is, when a particle appears at different locations in a fluid domain over a period of time, then the particle has traversed the fluid domain and can be considered to have reached a mixed state. The method for obtaining the dynamic hyperuniformity factor λ is as follows: During the stirring process, the three-dimensional coordinates of the tracer ball are reconstructed using a dual-camera positioning method, and the three-dimensional mixed image is divided into eight equal-sized regions N. n , (n=1, 2, ....., 8), the region shape is a cuboid, and Matlab is used to calculate the number of tracer balls in each region, L N , (N=1,2......,8), find the variance of β1 in 8 regions, and use β 2 Indicated by the variance β 2 The correlation β 2 ~N d-λ To reflect the mixing uniformity, we can obtain λ, where λ represents the dynamic hyperuniform state factor, N = 8, and d represents the dimension. In three dimensions, d = 3.

[0064] Analyzing the λ value, λ ≥ 1 indicates a hyperuniform state; λ ≤ 0 indicates a disordered state; and 0 < λ < 1 indicates a dynamic hyperuniform state. Specifically, when 0 < λ < 1, the magnitude of λ can qualitatively describe the strength of the dynamic hyperuniform state: the closer λ is to 0, the stronger the disordered state; the closer λ is to 1, the stronger the dynamic hyperuniform state.

[0065] Example 3

[0066] In order to verify the technical feasibility and effect of the method proposed in this invention, a calculation example is performed. Three different experiments are selected for verification. In the mixed image, 8 equal-sized rectangular areas (N=8) are taken, and the number of tracer balls in each area is calculated using Matlab. If the number of tracer balls surrounded by the 8 areas is N1=30, N2=32, N3=35, N4=31, N5=33, N6=37, N7=34, and N8=36 respectively, the variance of β1 in the 8 areas is β 2 =5.25, d=3, N=8, and β 2 , d, and N are substituted into β 2 ~N d-λ We get λ = 2.20, which means that the mixing process belongs to the super-uniform state. If the number of tracer balls enclosed by the 8 regions is N1 = 39, N2 = 24, N3 = 35, N4 = 26, N5 = 22, N6 = 30, N7 = 38, and N8 = 25, the variance of β1 in the 8 regions is β 2 =35.78, d=3, N=8, and β 2 , d, and N are substituted into β 2 ~N d-λ We get λ = -0.77, which means that the mixing process belongs to disordered state. If the number of tracer balls surrounded by the 8 regions is N1 = 25, N2 = 30, N3 = 27, N4 = 32, N5 = 28, N6 = 33, N7 = 26, N8 = 31, the variance of β1 in the 8 regions is β 2 =7.45, d=3, N=8, and β 2 , d, and N are substituted into β 2 ~N d-λ It is obtained that λ = 0.97, that is, the mixing process belongs to a dynamic super-uniform state.

[0067] The present invention investigates the effects of impeller bottom height and motor speed on mixing performance through hydraulic experiments. The uniformity of the mixing system is determined by reconstructing the three-dimensional coordinates of the motion trajectory of tracer particles in a stirred reactor. High-speed dual cameras are used to capture images of fully developed flow patterns. A visual tracking algorithm is used to record the coordinate positions of the instantaneous images captured by the high-speed cameras. The two-dimensional image coordinates of the tracer particles are recorded in real time and then stored in a computer for later calculation of the three-dimensional coordinate reconstruction. The experiment employed six operating conditions, L1-L6, each lasting approximately five minutes. The dual cameras captured images at one frame per second, generating two sets of 300 point sets. Three-dimensional reconstruction was then performed to obtain the three-dimensional coordinates, which were then normalized to eliminate undesirable effects caused by singular sample data. A dynamic hyperuniformity method was used to calculate the time it takes for a random, irregular point set to reach a dynamic hyperuniform state. The speed of the dynamic hyperuniform state was compared to determine the mixing performance of the mixing system. The time of dynamic super-uniform state obtained by dynamic super-uniform method can represent the mixing time when each group of working conditions reaches the mixing uniform state.

[0068] The dynamic hyperuniformity method was used to determine the temporal evolution and trend of the dynamic hyperuniform state for each operating condition, as well as the time I for each operating condition to reach the dynamic hyperuniform state. A smaller I indicates a more uniform mixing system. The closer the calculated dynamic hyperuniform factor λ approaches the critical value of 1, the more the mixing approaches the dynamic hyperuniform state. Operating condition L5 has the shortest time I to reach the dynamic hyperuniform state, I = 53s, while operating condition L6 has the longest time I = 75s. According to the dynamic hyperuniformity method, operating condition L5 has the highest mixing efficiency. Correspondingly, operating condition L5 also has the shortest mixing time, measured by conductivity, at 62s, demonstrating the effectiveness and feasibility of the dynamic hyperuniformity method.

[0069] The embodiments described above are merely descriptions of preferred embodiments of the present invention and are not intended to limit the scope of the present invention. Without departing from the spirit of the present invention, various modifications and improvements made to the technical solutions of the present invention by persons skilled in the art should fall within the scope of protection defined by the claims of the present invention.

Claims

1. A method for determining dynamic hyperuniformity, characterized in that: The method comprises: S1, using high-speed camera to extract the two-dimensional image coordinates of the tracer particle trajectory in the reactor; S2, extracting two sets of two-dimensional image coordinates of the same tracer particle trajectory acquired by two high-speed cameras, performing three-dimensional reconstruction on the two sets of two-dimensional image coordinates and then performing coordinate normalization processing; S3. Based on the tracer particle motion trajectory and the normalized three-dimensional image coordinates, the dynamic hyperuniform method is used to determine the uniformity of the mixing system, specifically including: Based on the motion trajectory of the tracer particles and the normalized three-dimensional image coordinates, the dynamic hyperuniform state time I when the dynamic hyperuniform factor λ approaches the critical value 1 is calculated. The dynamic hyperuniform state time is used to characterize the mixing uniformity of the mixing system. The calculation method of the dynamic hyperuniform factor λ includes: Divide the three-dimensional image into equal-sized regions using the three-dimensional image coordinates; MATLAB is used to calculate the number and variance of the tracer balls in each region, and a dynamic hyperuniform factor λ is obtained based on the number and variance of the tracer balls; when Indicates that the mixture is in a super-uniform state; Indicates that the mixture is in a disordered state; Indicates that the mixture is in a dynamic super-uniform state; when When λ is close to 0, the size of λ can qualitatively describe the strength of the dynamic superuniform state. The closer λ is to 0, the stronger the disordered state is; the closer λ is to 1, the stronger the dynamic superuniform state is.

2. The method for determining dynamic hyperuniformity according to claim 1, wherein: In S1, the content of extracting the two-dimensional image coordinates of the tracer particle trajectory in the reactor using high-speed photography specifically includes: The tracer particle trajectories in the reactor are captured using high-speed cameras, where the high-speed cameras are positioned perpendicular to each other and on the same horizontal plane; The visual tracking algorithm is used to locate the coordinates of the tracer particle motion trajectory on the captured image to obtain the two-dimensional image coordinates of the tracer particle motion trajectory.

3. The method for determining dynamic hyperuniformity according to claim 1, wherein: In S2, the coordinate normalization processing after the three-dimensional reconstruction of the two sets of two-dimensional image coordinates specifically includes: Assume that the two sets of extracted two-dimensional image coordinates are (x1, y1) and (x2, y2); The two sets of two-dimensional image coordinates are reconstructed into three dimensions, and the three-dimensional coordinates are marked as (X, Y, Z), where X=x1, Y=x2, and Z=(y1+y2) / 2 The obtained three-dimensional coordinates are normalized.

4. A dynamic hyperuniformity determination system, the system being used to implement the determination method according to any one of claims 1 to 3, characterized in that: The system includes: a coordinate extraction module, a coordinate preprocessing module and a uniformity judgment module; The coordinate extraction module is used to extract the two-dimensional image coordinates of the tracer particle trajectory in the reactor using high-speed photography; The coordinate preprocessing module is used to extract two sets of two-dimensional image coordinates of the same tracer particle trajectory obtained by two high-speed cameras, and perform coordinate normalization processing after three-dimensional reconstruction of the two sets of two-dimensional image coordinates; The uniformity judgment module is used to judge the mixing uniformity of the mixing system using a dynamic hyperuniformity method based on the motion trajectory of the tracer particles and the three-dimensional coordinates after normalization.

5. The dynamic super-uniformity determination system according to claim 4, characterized in that: The coordinate extraction module also includes a tracer particle track image acquisition submodule and a track coordinate extraction submodule; The tracer particle track image acquisition submodule is used to capture the tracer particle track in the reactor using high-speed cameras, with the positions of the high-speed cameras being perpendicular to each other and on the same horizontal plane; The trajectory coordinate extraction submodule is used to use a visual tracking algorithm to locate the coordinates of the tracer particle motion trajectory on the captured image to obtain the two-dimensional image coordinates of the tracer particle motion trajectory.

6. The dynamic super-uniformity determination system according to claim 4, characterized in that: The coordinate preprocessing module also includes: a dimension conversion submodule and a coordinate normalization submodule; The dimension conversion submodule is used to set the two sets of extracted two-dimensional image coordinates as (x1, y1) and (x2, y2) respectively; reconstruct the two sets of two-dimensional image coordinates into three dimensions, and the three-dimensional coordinates are denoted as (X, Y, Z), where X=x1, Y=x2, and Z=(y1+y2) / 2; The coordinate normalization submodule is used to perform normalization processing on the obtained three-dimensional coordinates.

7. The dynamic super-uniformity determination system according to claim 4, characterized in that: The uniformity judgment module calculates the dynamic hyperuniform state time I when the dynamic hyperuniform factor λ approaches the critical value 1 based on the tracer particle motion trajectory and the normalized three-dimensional image coordinates, and uses the dynamic hyperuniform state time to characterize the mixing uniformity of the mixing system.

Images