Method, device, equipment and medium for analyzing suspended sediment transport distribution data

Abstract

The present invention relates to the field of data processing, and discloses a method, device, equipment and medium for analyzing the transport distribution data of suspended sediment, which are used to improve the accuracy of the transport distribution data analysis of suspended sediment. The method includes: obtaining first flow field information preset in the flow field in the first time step; performing force analysis on a particle group preset in the flow field according to the first flow field information, determining the water flow force and the interaction force between particles corresponding to the particle group, and generating acceleration data of the particle group according to the water flow force and the interaction force between particles; performing motion speed analysis on the particle group according to the acceleration data, determining the corresponding motion speed, and determining the transport distribution data of the particle group in the flow field according to the motion speed; determining the force corresponding to the flow field through the water flow force to generate the second flow field information of the flow field in the second time step; and re-determining the transport distribution data of the particle group in the flow field based on the second flow field information.

Description

Technical Field

[0001] The present invention relates to the field of data processing technology, and in particular to a method, device, equipment and medium for analyzing the transport distribution data of suspended sediment. Background Art

[0002] Suspended sediment in sandy rivers accounts for more than half of the total sediment volume, and seriously affects river water quality. Therefore, in-depth research on the movement characteristics of suspended sediment is of great significance. In the field of water conservancy, the widely used suspended sediment movement model generally generalizes the sediment model into a continuous medium with a certain concentration, and constructs a mathematical model based on the scalar transport equation. That is, the distribution of suspended sediment is assumed to be the concentration of a continuous medium, and the suspended sediment transport process is solved by coupling the convection-diffusion equation of the concentration with the convection-diffusion equation of the flow field.

[0003] However, the assumption of continuous medium concentration deviates from the essence of the discrete characteristics of sediment to a certain extent, and there are difficulties in the selection of the diffusion coefficient. The accuracy of the numerical simulation of suspended sediment movement based on the scalar transport equation is heavily dependent on the selection of the suspended sediment diffusion coefficient. However, the sediment concentration diffusion coefficient is a physical quantity that characterizes the diffusion intensity of suspended sediment in the water body. It is heavily dependent on hydrodynamic conditions and the physical properties of sediment and is often difficult to determine. Therefore, the numerical simulation of suspended sediment based on the scalar transport equation often results in inaccurate predictions. At the same time, sediment particles have natural discreteness. Ignoring the discrete characteristics of sediment in the scalar transport model based on sediment concentration will introduce errors again. Therefore, it is urgent to carry out in-depth and refined simulations to improve simulation accuracy. Summary of the invention

[0004] In view of this, an embodiment of the present invention provides a method, device, equipment and medium for analyzing the transport distribution data of suspended sediment, which solves the technical problem of low accuracy when analyzing the transport distribution data of suspended sediment.

[0005] A first aspect of the present invention provides a method for analyzing the transport distribution data of suspended sediment, the method comprising: acquiring first flow field information of a flow field preset in a first time step; performing a force analysis on a particle group preset in the flow field according to the first flow field information, determining the water flow force and the interaction force between particles corresponding to the particle group, and generating acceleration data of the particle group according to the water flow force and the interaction force between particles; performing a motion speed analysis on the particle group according to the acceleration data, determining the motion speed corresponding to the particle group, and determining the transport distribution data of the particle group in the flow field according to the motion speed; determining the force corresponding to the flow field through the water flow force, and generating second flow field information of the flow field in a second time step through the force; using the second flow field information as new first flow field information, and re-determining the transport distribution data of the particle group in the flow field based on the second flow field information.

[0006] The suspended sediment transport distribution data analysis method provided by the present invention has a discrete processing of sediment particles, which makes it more advantageous than the continuous method in calculating and analyzing the behavior of sediment particles. A small number of calculation particles replace a large number of particles with the same properties, and collisions between particles are simulated, which reduces the amount of calculation and improves the calculation efficiency and accuracy, avoids the problem of selecting the diffusion coefficient of the numerical simulation of suspended sediment movement, and improves the accuracy of the simulation.

[0007] In combination with the first aspect, in a first implementation of the first aspect, obtaining the first flow field information preset in the flow field in the first time step includes: obtaining velocity information of a particle group preset in the flow field and velocity information of the flow field; performing flow field velocity analysis on the flow field through the velocity information to generate a flow field velocity corresponding to the flow field; and using the flow field velocity and the pressure information as the first flow field information preset in the flow field in the first time step.

[0008] In this scheme, by updating the flow field velocity and pressure information of the flow field in the current time step and using the flow field velocity and pressure information as the flow field information corresponding to the flow field, the quality of flow field initialization is improved and the efficiency of numerical solution is improved.

[0009] In combination with the first implementation of the first aspect, in the second implementation of the first aspect, the force analysis of the particle group preset in the flow field is performed according to the first flow field information to determine the water flow force and the inter-particle interaction force corresponding to the particle group, and the acceleration data of the particle group is generated according to the water flow force and the inter-particle interaction force, including: performing momentum increment analysis on the particle group through the flow field velocity and the pressure information to determine the momentum increment corresponding to the particle group; performing water flow force analysis on the particle group based on the momentum increment to generate the water flow force corresponding to the particle group; gridding the particle group to generate particle group grid information, and calculating the particle phase volume fraction of the particle group through the particle group grid information to obtain the particle phase volume fraction corresponding to the particle group; performing interaction force analysis on the particle group through the particle phase volume fraction to determine the inter-particle interaction force corresponding to the particle group; calculating the total force on the particle group according to the water flow force and the inter-particle interaction force; and calculating the acceleration of the particle group according to the total force.

[0010] In this scheme, the interaction force between particles is introduced, and particles and grids are mixed for operation, which improves the efficiency and accuracy of particle motion state analysis.

[0011] In combination with the second implementation of the first aspect, in the third implementation of the first aspect, the calculation of the total force exerted on the particle group based on the water flow force and the inter-particle interaction force includes: calculating the flow field force exerted on the particle group based on the water flow force; calculating the stress within the particle group exerted on each particle in the particle group based on the inter-particle interaction force; and calculating the total force exerted on the particle group based on the flow field force and the stress within the particle group exerted on each particle.

[0012] In this scheme, the shear test results of coarse particles suspended between concentric cylinders, that is, the expressions of particle shear stress and inter-particle discrete stress in the inertial action zone, are used for calculation, and then the inter-particle interaction force corresponding to the particle group is obtained, further improving the accuracy of the motion simulation of sediment particles.

[0013] In combination with the first aspect, in a fourth implementation of the first aspect, the movement speed analysis of the particle group is performed according to the acceleration data to determine the movement speed corresponding to the particle group, and the transport distribution data of the particle group in the flow field is determined according to the movement speed, including: performing a time integration operation on the acceleration data to obtain the movement speed of the particle group; performing a time integration operation on the movement speed to obtain the transport distribution data of the particle group in the flow field.

[0014] In combination with the first aspect, in the fifth implementation of the first aspect, the force corresponding to the flow field is determined by the water flow force, and the second flow field information within the second time step is generated by the force, including: determining the force of the particle group acting on the flow field according to the water flow force; converting the force into a volume force to determine the volume force corresponding to the flow field; generating flow field data for the flow field within the second time step based on the volume force to obtain the second flow field information.

[0015] In combination with the first aspect or any one of the first to fifth embodiments of the first aspect, in the sixth embodiment of the first aspect, the transport distribution data analysis method of suspended sediment also includes: taking the second time step as a new first time step, and within a specified time step, recording the transport distribution data of the particle group in the flow field within each first time step to obtain multiple transport distribution data; performing data analysis on the multiple transport distribution data to obtain transport distribution data analysis results of suspended sediment; performing data conversion on the transport distribution data analysis results based on preset spatial position update parameters to obtain spatial transport distribution data of the suspended sediment.

[0016] In this scheme, by coupling the sediment movement velocity with the sediment diffusion coefficient, the poor analysis effect of sediment transport distribution data caused by the selection of the sediment diffusion coefficient can be avoided, and the accuracy of suspended sediment transport distribution data analysis can be further improved.

[0017] According to a second aspect, an embodiment of the present invention provides a suspended sediment transport distribution data analysis device, comprising:

[0018] An acquisition module, used for acquiring first flow field information preset in a first time step of the flow field;

[0019] A generating module, configured to perform force analysis on a particle group preset in the flow field according to the first flow field information, determine a water flow force and an interaction force between particles corresponding to the particle group, and generate acceleration data of the particle group according to the water flow force and the interaction force between particles;

[0020] An analysis module, configured to analyze the movement speed of the particle group according to the acceleration data, determine the movement speed corresponding to the particle group, and determine the transport distribution data of the particle group in the flow field according to the movement speed;

[0021] A determination module, used to determine the force corresponding to the flow field through the water flow force, and generate second flow field information of the flow field in a second time step through the force;

[0022] An updating module is used to use the second flow field information as new first flow field information, and to redetermine the transport distribution data of the particle group in the flow field based on the second flow field information.

[0023] According to the third aspect, an embodiment of the present invention provides an electronic device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the computer instructions to execute the method for analyzing the transport distribution data of suspended sediment described in the first aspect or any one of the embodiments of the first aspect, or to execute the method for analyzing the transport distribution data of suspended sediment described in the second aspect or any one of the embodiments of the second aspect.

[0024] According to the fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the suspended sediment transport distribution data analysis method described in the first aspect or any one of the embodiments of the first aspect, or to execute the suspended sediment transport distribution data analysis method described in the second aspect or any one of the embodiments of the second aspect. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0026] Figure 1 It is a flow chart of a method for analyzing suspended sediment transport distribution data in an embodiment of the present invention;

[0027] Figure 2 Schematic diagram of the comparison between the vertical sediment concentration distribution of different sections and the analytical solution in the embodiment of the present invention;

[0028] Figure 3 It is a flow chart of calculating the acceleration of a particle group according to the water flow force and the interaction force between particles in an embodiment of the present invention;

[0029] Figure 4 A flow chart of updating the transport distribution data of a particle group in a flow field according to the movement speed in an embodiment of the present invention;

[0030] Figure 5 A flowchart of calculating the second flow field information in an embodiment of the present invention;

[0031] Figure 6Schematic diagram of a suspended sediment transport distribution data analysis device according to an embodiment of the present invention;

[0032] Figure 7 A schematic diagram of the hardware structure of an electronic device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0033] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0034] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicating the orientation or positional relationship, are based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance.

[0035] In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0036] For ease of understanding, the specific process of the embodiment of the present invention is described below. Figure 1 , Figure 1 FIG. 1 is a flow chart of a method for analyzing suspended sediment transport distribution data according to an embodiment of the present invention. Figure 1 As shown, the flow chart includes the following steps:

[0037] Step S101: obtaining first flow field information preset in a first time step;

[0038] It is understandable that the execution subject of the present invention may be a suspended sediment transport distribution data analysis device, or a terminal or a server, which is not limited here. The embodiment of the present invention is described by taking a server as the execution subject as an example.

[0039] It should be noted that when obtaining the first flow field information preset in the first time step, the server mainly performs information analysis and calculation through the fluid motion equation, which refers to the Navier-Stokes equation of incompressible fluid. The server calculates the pressure and flow velocity of the flow field preset in the current time step through the Navier-Stokes equation, where the Navier-Stokes equation is as follows:

[0040]

[0041]

[0042] Among them, x i (i=1,2,3) represent the longitudinal, transverse and vertical coordinates respectively. i (i=1,2,3) represents the longitudinal, transverse and vertical velocity components respectively; t is time; ρ is the fluid density; Fb i The volume force of the sediment particles on the fluid (i=1, 2, 3) represents the longitudinal, lateral and vertical components respectively. In this embodiment, based on the influence of the particle group on the dissipation of turbulent kinetic energy in the water flow field, this influence is reflected in the model by modifying the traditional k-ε turbulence model, which can improve the accuracy of the k-ε turbulence model. At the same time, the server calculates the flow velocity and pressure of the flow field through the equation to determine the corresponding first flow field information, so as to facilitate the subsequent analysis of the motion state of the sediment particles and improve the accuracy of the motion state analysis of the sediment particles.

[0043] In one example, the first flow field information may be initial flow field information of the flow field, and the initial state of the flow field may be determined through the first flow field information.

[0044] Step S102: performing force analysis on a preset particle group in the flow field according to the first flow field information, determining the water flow force and the interaction force between particles corresponding to the particle group, and generating acceleration data of the particle group according to the water flow force and the interaction force between particles;

[0045] It should be noted that the flow field information mainly includes velocity and pressure. The force exerted by the water flow on the particle group is calculated. The water flow force mainly includes drag force, buoyancy and gravity. The interaction force between particles is considered based on the solid phase normal stress model. The force exerted by the flow field on the particle group and the stress exerted by the particles in the particle group are summed to obtain the total force exerted on the particles, and the acceleration of the particle group is deduced based on the total force. In this embodiment, the interaction force between particles is introduced, and the particles and the grid are mixed and operated, which improves the efficiency and accuracy of the particle motion state analysis.

[0046] Step S103: analyzing the movement speed of the particle group according to the acceleration data, determining the movement speed corresponding to the particle group, and determining the transport distribution data of the particle group in the flow field according to the movement speed;

[0047] Specifically, the new movement speed of the particle group is updated by time integration of the acceleration of the particle group, and the distribution of the particle group in the flow field is predicted by time integration of the movement speed of the particle group again, so as to obtain the transport distribution data of the particle group in the flow field.

[0048] Step S104: determining the force corresponding to the flow field through the water flow force, and generating second flow field information of the flow field in the second time step through the force;

[0049] It should be noted that, since the flow field receives the force of the particle group, its magnitude is the same as the force of the water flow on the particle group, but in the opposite direction. At the same time, the server adds this force to the right side of the above fluid motion equation in the form of a volume force source term, and uses the fluid motion equation containing the volume force to calculate the flow field information of the next time step, that is, the second flow field information in the above second time step. That is, the second time step can be understood as the next time step of the first time step. The second flow field information and the first flow field information respectively represent the flow field information of the flow field in different time steps. The first flow field information is the flow field information of the flow field in the first time step, and the second flow field information is the flow field information of the flow field in the second time step.

[0050] Step S105: taking the second flow field information as new first flow field information, and re-determining the transport distribution data of the particle group in the flow field based on the second flow field information.

[0051] Specifically, by iterative calculations at different time steps, the movement state of suspended sediment is simulated, thereby avoiding the selection of sediment diffusion coefficient and improving the accuracy of sediment numerical simulation. For example, in a straight river channel, the average inflow velocity is 0.1735m / s, and the sediment particle D is 0.0005m. The upstream and downstream of the straight channel adopt periodic boundary conditions. After the calculation time exceeds 50L / U (where L is the channel length and U is the average flow velocity), the water and sediment flow field is fully developed, such as Figure 2 As shown, the vertical sediment concentration distribution is statistically analyzed at different river sections, which is consistent with the analytical solution. Therefore, it is verified that this embodiment can be applied to the calculation of suspended sediment concentration, and there is no need to manually select the sediment diffusion coefficient.

[0052] By executing the above steps, the discrete processing of sediment particles makes it more advantageous than the continuous method in calculating and analyzing the behavior of sediment particles. Replacing a large number of particles with the same properties with a small number of calculation particles and simulating the collision between particles with the solid phase stress model significantly reduces the amount of calculation and improves the calculation efficiency and accuracy.

[0053] In a specific embodiment, the process of executing step S101 may specifically include the following steps:

[0054] (1) obtaining the velocity information of the particle group preset in the flow field in the first time step and the pressure information of the flow field;

[0055] (2) performing flow field velocity analysis on the flow field through velocity information to generate a flow field velocity corresponding to the flow field;

[0056] (3) The flow velocity and pressure information of the flow field is used as the first flow field information preset in the first time step.

[0057] Specifically, the server first obtains the velocity information of the particle group preset in the flow field and the pressure information of the flow field. In one example, in order to improve the accuracy of the pressure information, the pressure information can be obtained in the following manner: according to the velocity information of the particle group and the fluid density of the flow field, the initial value of the pressure information in the current time step in the state to be initialized is calculated, and the initial value of the fluid velocity under the incoming flow condition in the state to be initialized is calculated, so as to correct the initial value of the pressure information based on the initial value of the flow field velocity. The initial value of the flow field temperature and the initial value of the flow field velocity under the incoming flow condition in the state to be initialized are calculated, and the initial value of the flow field velocity vector under the incoming flow condition in the state to be initialized is calculated, and then the above formula (1) and formula (2) are used for interpolation calculation to obtain the pressure information and velocity vector on the corresponding grid in the current time step, and the flow field velocity and pressure information of the flow field in the current time step are updated at the same time, and the flow field velocity and pressure information of the flow field are used as the first flow field information preset in the first time step of the flow field, thereby improving the quality of flow field initialization and the efficiency of numerical solution.

[0058] In a specific embodiment, if Figure 3 As shown, the process of executing step S102 may specifically include the following steps:

[0059] S201: analyzing the momentum increment of the particle group through the flow field velocity and pressure information to determine the momentum increment corresponding to the particle group;

[0060] S202: analyzing the water flow force on the particle group based on the momentum increment, and generating the water flow force corresponding to the particle group;

[0061] Specifically, the server initializes the position and velocity of particles in the particle swarm, and updates the velocity and position of individual particles and particle swarms. The current environment calculates the momentum increment corresponding to the particle swarm, and determines the corresponding water flow force through the momentum increment. It should be noted that the momentum increment refers to the momentum of the particle in the final state minus the momentum of the initial state. According to the law of conservation of energy, the momentum increment can be used to determine the momentum change caused by the water flow force, and then solve the initial water flow force.

[0062] S203: gridding the particle swarm to generate particle swarm grid information, and calculating the particle phase volume fraction of the particle swarm using the particle swarm grid information to obtain the particle phase volume fraction corresponding to the particle swarm;

[0063] S204: performing interaction force analysis on the particle group according to the particle phase volume fraction, and determining the inter-particle interaction force corresponding to the particle group;

[0064] Specifically, the water flow force is calculated based on the first flow field information, wherein it should be noted that when there are multiple particles in the fluid, the force of the particles will be different from that of a single particle. The movement of any particle may be affected by other particles. The main forms of interaction between particles are contact, position exchange and collision between particles. At the same time, the presence of a large number of particles will affect the flow characteristics of the liquid phase, and the change of the latter will in turn affect the movement of the particles. Therefore, in this embodiment, the particle group is gridded. For the situation where multiple particles exist, the Lagrange acceleration equation needs to be corrected so that the interaction between particles can be considered, and the particle phase volume fraction is introduced for force analysis. Specifically, in this embodiment, the particles are corrected by the inertial effect of the fluid and the particle concentration. For the interaction between particles, the shear test results of coarse particles suspended between concentric cylinders, that is, the expression of particle shear stress and inter-particle discrete stress in the inertial action area, are used for calculation, and then the particle interaction force corresponding to the particle group is obtained, which further improves the accuracy of the motion simulation of sediment particles.

[0065] S205: Calculate the total force on the particle group according to the water flow force and the interaction force between particles;

[0066] S206: Calculate the acceleration of the particle group according to the total force.

[0067] Specifically, all particles are first checked one by one to find the maximum and minimum values ​​of the position coordinates of these particles in the three directions of X, Y, and Z. The distribution span of the particle set in the three directions is calculated by the difference between the maximum and minimum values. Then, the distribution span in each direction is divided into equal intervals with the length of the smooth core radius as the spacing. Finally, a three-dimensional space grid with the smooth core radius as the grid unit length is established. At this time, the number of sediment particles is determined. For a regular water body, the number of particles distributed in each direction of X, Y, and Z is determined, and then the number of the three directions is multiplied. The number of particle distributions can be selected arbitrarily. Increasing the number can make the effect realistic, but it will lead to increased resource consumption. Reducing the number can speed up the processing speed, but it will lead to less realistic simulation effects. The specific number can be determined based on the resources used and the expected effect. At the same time, the total force on the particle group is calculated based on the water flow force and the interaction force between particles, and the acceleration of the particle group is calculated based on the total force.

[0068] In a specific embodiment, the process of executing step S205 may specifically include the following steps:

[0069] (1) Calculate the flow field force on the particle group based on the water flow force;

[0070] Specifically, since the properties of each particle are obtained by interpolating other particles in the particle support domain, in order to facilitate the management of particles in the support domain, we establish an adjacent particle list for each particle, which contains the index values ​​of all particles located in the particle support domain, and obtain the flow field force corresponding to the particle group according to the water flow force determined in the above steps. In this embodiment, the space where the particle group is located is gridded to accelerate particle search and improve calculation efficiency.

[0071] (2) calculating the stress within the particle group on each particle in the particle group based on the interaction force between the particles;

[0072] (3) Calculate the total force on the particle group based on the flow field force and the stress on each particle within the particle group.

[0073] Specifically, the force on each particle can be decomposed into pressure, gravity and viscosity. The stress on each particle in the particle group is calculated by the above calculation formula. It should be noted that in order to complete the calculation of the force on the particle, the distribution of the pressure field, that is, the pressure at the location of the particle, needs to be obtained. Therefore, this application uses a modified form of the ideal gas equation to calculate the pressure of the particle. The specific process is as follows: first calculate the density of the particle, and at the same time calculate the pressure, gravity, pressure, viscosity at the location of the particle, and the total force on the particle.

[0074] In a specific embodiment, if Figure 4 As shown, the above step S103 specifically includes the following steps:

[0075] S301: Perform time integration operation on the acceleration data to obtain the movement speed of the particle group;

[0076] S302: Perform time integration operation on the movement speed to obtain the transport distribution data of the particle group in the flow field.

[0077] It should be noted that acceleration is the ratio of the change in velocity to the time taken for this change to occur. It is a physical quantity that describes how fast an object's velocity changes. Acceleration is a vector, and its direction is the direction of the change in the velocity of the object, which is the same as the direction of the resultant external force. Suppose that we want to calculate the area of ​​a<=x<=b under any curve f(x), and divide this interval into n slices, each with a width of (ba) / n. The left coordinate of the i-th slice is xi, and the height is f(xi). The total area is obtained by adding up all the small areas in an integral way. This total area is the change in the velocity of the particle group. At the same time, according to the initial velocity of the particle group, the velocity of the particle group is determined, and then the velocity is integrated over time to obtain the transport distribution data of the particle group in the flow field.

[0078] In a specific embodiment, if Figure 5 As shown, the above step S104 specifically includes the following steps:

[0079] S401: Determine the force exerted by the particle group on the flow field according to the force of the water flow;

[0080] S402: converting the acting force into a volume force to determine the volume force corresponding to the flow field;

[0081] S403: Generate flow field data for the flow field in the second time step based on the body force to obtain second flow field information.

[0082] It should be noted that the force exerted by the particle group on the flow field is the same in magnitude but opposite in direction to the force exerted by the water flow on the particle group. Therefore, the force exerted by the particle group on the flow field can be determined based on the force exerted by the water flow. At the same time, the presence of the particle group affects the dissipation of turbulent kinetic energy in the water flow field. This influence is reflected in the model by modifying the traditional k-ε turbulence model. The modified k-ε turbulence model is shown in the following formula:

[0083]

[0084]

[0085] The empirical constant C μ =0.09, C 1ε =1.44, C 2ε =1.92,σ k =1.0,σ ε =1.3, the last term of the control equation (3) is the effect of the particle on the turbulent kinetic energy and dissipation rate transport equation of water flow. It should be noted that, among them, f i =(u i -u pi ) / Γ p , Γ pis the time that particle motion lags behind water motion, Γ p =4ρ p d p / 3(ρC D |uu p |), The drag force on the sand particles in is the actual drag force, is the Stokes drag force, g represents the acceleration due to gravity 9.81m / s2, ν t is the eddy viscosity coefficient, k is the turbulent energy, and ε is the turbulent energy dissipation rate. The volume force corresponding to the force is determined through the above turbulence model control equation, and then the flow field information in the next time step is calculated according to the volume force and the fluid motion equation.

[0086] In a specific embodiment, after the above step S105, the following steps are also included:

[0087] (1) taking the second time step as a new first time step, and within a specified time step, recording the transport distribution data of the particle group in the flow field within each first time step to obtain a plurality of transport distribution data;

[0088] Specifically, an iterative method is used to calculate the transport distribution data in the next time step based on the transport distribution data in the current first time step, obtain the transport distribution data in the second time step, and use the transport distribution data in the second time step as the transport distribution data in the first time step of the next round. And so on, and in the process of calculation, it is determined whether the time step reaches the specified time step. If it reaches the specified time step, the transport distribution data of the particle group in the flow field in each first time step is output to obtain multiple transport distribution data. In one example, if the time step does not reach the specified time step, the iteration continues.

[0089] (2) Analyze multiple transport distribution data to obtain the transport distribution data analysis results of suspended sediment

[0090] Specifically, data analysis can be performed by fitting multiple water and sediment flows to multiple transport distribution data, thereby obtaining the displacement data analysis results of suspended sediment;

[0091] (3) Based on the preset spatial position update parameters, the transport distribution data analysis results are converted to obtain the spatial transport distribution data of suspended sediment.

[0092] Specifically, the control equation for updating the spatial position of sediment particles in the next time step is shown in formula (5).

[0093]

[0094] Among them, X p0 is the position of the sediment particle at time t, is the velocity of the sediment particles at time t+Δt. is the sediment moving speed U p and sediment diffusion coefficient U' p The sum of U p is obtained by the Lagrangian method. Since the random movement speed of diffusion satisfies the mean of zero and the standard deviation of Gaussian distribution, specifically, the server, through multiple calculation results, performs water and sediment multi-flow fitting, obtains the suspended sediment transport distribution data analysis results, and converts the transport distribution data analysis results through the above control equation, obtains the suspended sediment spatial transport distribution data, and by coupling the sediment movement speed with the sediment diffusion coefficient, it can avoid the poor analysis effect of sediment transport distribution data caused by the selection of sediment diffusion coefficient, and further improve the accuracy of suspended sediment transport distribution data analysis. At the same time, in the first few time steps of the calculation, because the water level and flow velocity are assumed at the beginning of the calculation, there is a big difference from the actual situation, which often causes the program to diverge or cause certain calculation errors. After the simulated river section is basically stable, several time steps are calculated to enable the water flow and sediment to exchange information in real time. When the flow velocity and water level have changed very little in each time step, the water and sediment multi-flow fitting is performed to obtain the suspended sediment transport distribution data analysis results.

[0095] The embodiment of the present invention also provides a suspended sediment transport distribution data analysis device, such as Figure 6 As shown, the suspended sediment transport distribution data analysis device specifically includes:

[0096] An acquisition module 501 is used to acquire first flow field information preset in a first time step of the flow field;

[0097] A generation module 502 is used to perform force analysis on a particle group preset in the flow field according to the first flow field information, determine the water flow force and the interaction force between particles corresponding to the particle group, and generate acceleration data of the particle group according to the water flow force and the interaction force between particles;

[0098] The analysis module 503 is used to analyze the movement speed of the particle group according to the acceleration data, determine the movement speed corresponding to the particle group, and determine the transport distribution data of the particle group in the flow field according to the movement speed;

[0099] A determination module 504, configured to determine a force corresponding to the flow field through a water flow force, and generate second flow field information of the flow field in a second time step through the force;

[0100] The updating module 505 is used to use the second flow field information as the new first flow field information, and re-determine the transport distribution data of the particle group in the flow field based on the second flow field information.

[0101] The further functional description of each of the above modules is the same as that of the above corresponding method embodiments and will not be repeated here.

[0102] Through the coordinated cooperation of the above-mentioned components, the fluid motion equation containing volume force is used to calculate the flow field information of the next time step, that is, the flow field information in the second time step mentioned above. Through iterative calculations of different time steps, the movement state of suspended sediment is simulated, thereby avoiding the selection of sediment diffusion coefficient and improving the accuracy of sediment numerical simulation.

[0103] The embodiment of the present invention further provides an electronic device, such as Figure 7 As shown, the electronic device may include a processor 601 and a memory 602, wherein the processor 601 and the memory 602 may be connected via a bus or other means. Figure 7 The example of connecting through bus is taken in the following.

[0104] The processor 601 can be a central processing unit (CPU). The processor 601 can also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components and other chips, or a combination of the above-mentioned various types of chips. The memory 602, as a non-transient computer-readable storage medium, can be used to store non-transient software programs, non-transient computer executable programs and modules, such as program instructions / modules corresponding to the method in the embodiment of the present invention. The processor 601 executes various functional applications and data processing of the processor by running the non-transient software programs, instructions and modules stored in the memory 602, that is, implementing the above-mentioned method.

[0105] The memory 602 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required by at least one function; the data storage area may store data created by the processor 601, etc. In addition, the memory 602 may include a high-speed random access memory, and may also include a non-transitory memory, such as at least one disk storage device, a flash memory device, or other non-transitory solid-state storage device. In some embodiments, the memory 602 may optionally include a memory remotely arranged relative to the processor 601, and these remote memories may be connected to the processor 601 via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0106] One or more modules are stored in the memory 602, and when executed by the processor 601, the above method is performed.

[0107] The embodiment of the present invention further provides a non-transitory computer storage medium, which stores computer executable instructions, and the computer executable instructions can execute the personnel counting method in any of the above method embodiments. Among them, the storage medium can be a disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory (Flash Memory), a hard disk (HDD) or a solid-state drive (SSD), etc.; the storage medium can also include a combination of the above types of memory.

[0108] Those skilled in the art will appreciate that all or part of the processes in the above-mentioned embodiment method can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, wherein the storage medium can be a magnetic disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk drive (HDD) or a solid-state drive (SSD), etc.; the storage medium can also include a combination of the above-mentioned types of memories.

[0109] The above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit the same. Although the present invention has been described in detail with reference to the above embodiments, a person skilled in the art should understand that the specific implementation modes of the present invention can still be modified or replaced by equivalents, and any modification or equivalent replacement that does not depart from the spirit and scope of the present invention should be included in the scope of the claims of the present invention.

Claims

1. A suspended sediment transport distribution data analysis method, It is characterized in that include: Acquire velocity information of a particle group preset in the flow field in a first time step and pressure information of the flow field; Performing flow field velocity analysis on the flow field using the velocity information to generate a flow field velocity corresponding to the flow field; Using the pressure information and the flow field velocity as first flow field information preset in the flow field in a first time step; Performing momentum increment analysis on the particle group according to the flow field velocity and the pressure information to determine the momentum increment corresponding to the particle group; Performing water flow force analysis on the particle group based on the momentum increment to generate a water flow force corresponding to the particle group; Performing grid division on the particle group to generate particle group grid information, and calculating the particle phase volume fraction of the particle group using the particle group grid information to obtain the particle phase volume fraction corresponding to the particle group; Correcting the inertial effect of the fluid and the particle concentration on the particle group, analyzing the interaction force of the particle group according to the particle phase volume fraction, and determining the inter-particle interaction force corresponding to the particle group; Calculating the total force exerted on the particle group according to the water flow force and the interaction force between the particles; Calculating the acceleration of the particle group according to the total force; Performing a motion speed analysis on the particle group according to the acceleration data to determine a motion speed corresponding to the particle group, and determining transport distribution data of the particle group in the flow field according to the motion speed; Based on the influence of the particle group on the water flow turbulent kinetic energy and dissipation rate transport equation, the force of the particle group on the flow field is determined according to the water flow force; Converting the acting force into a volume force to determine the volume force corresponding to the flow field; Generate flow field data for the flow field in a second time step based on the body force to obtain second flow field information; using the second flow field information as new first flow field information, and re-determining the transport distribution data of the particle group in the flow field based on the second flow field information; The method of analyzing the interaction force of the particle group by the volume fraction of the particle phase to determine the inter-particle interaction force corresponding to the particle group includes: obtaining the shear test results of coarse particles suspended between concentric cylinders, and calculating the inter-particle interaction force corresponding to the particle group by using the expressions of particle shear stress and inter-particle discrete stress in the inertial action zone.

2. The suspended sediment transport distribution data analysis method according to claim 1, It is characterized in that The calculating the total force exerted on the particle group according to the water flow force and the inter-particle interaction force comprises: Calculating the flow field force on the particle group according to the water flow force; Calculating the stress within the particle group on each particle in the particle group according to the inter-particle interaction force; The total force on the particle group is calculated according to the flow field force and the stress on each particle in the particle group.

3. The suspended sediment transport distribution data analysis method according to claim 1, It is characterized in that The step of analyzing the movement speed of the particle group according to the acceleration data, determining the movement speed corresponding to the particle group, and determining the transport distribution data of the particle group in the flow field according to the movement speed includes: Performing a time integration operation on the acceleration data to obtain the movement speed of the particle group; A time-integrated operation is performed on the movement speed to obtain the transport distribution data of the particle group in the flow field.

4. The suspended sediment transport distribution data analysis method according to any one of claims 1 to 3, It is characterized in that The suspended sediment transport distribution data analysis method further includes: The second time step is used as a new first time step, and within a specified time step, the transport distribution data of the particle group in the flow field in each first time step is recorded to obtain a plurality of transport distribution data; Performing data analysis on the plurality of transport distribution data to obtain a transport distribution data analysis result of suspended sediment; The transport distribution data analysis result is converted based on preset spatial position update parameters to obtain the spatial transport distribution data of the suspended sediment.

5. A suspended sediment transport distribution data analysis device, It is characterized in that The suspended sediment transport distribution data analysis device comprises: An acquisition module is used to acquire velocity information of a particle group preset in a flow field in a first time step and pressure information of the flow field; perform flow field velocity analysis on the flow field through the velocity information to generate a flow field velocity corresponding to the flow field; and use the pressure information and the flow field velocity as first flow field information preset in the flow field in the first time step; A generation module is used to perform momentum increment analysis on the particle group through the flow field velocity and the pressure information to determine the momentum increment corresponding to the particle group; perform water flow force analysis on the particle group based on the momentum increment to generate the water flow force corresponding to the particle group; perform grid division on the particle group to generate particle group grid information, and calculate the particle phase volume fraction of the particle group through the particle group grid information to obtain the particle phase volume fraction corresponding to the particle group; correct the inertia of the fluid and the particle concentration on the particle group, and correct the particle phase volume fraction according to the particle phase volume fraction. The particle group is subjected to interaction force analysis to determine the inter-particle interaction force corresponding to the particle group; the total force on the particle group is calculated according to the water flow force and the inter-particle interaction force; the acceleration of the particle group is calculated according to the total force; the interaction force analysis of the particle group by the particle phase volume fraction to determine the inter-particle interaction force corresponding to the particle group includes: obtaining the shear test results of coarse particles suspended between concentric cylinders, calculating the inter-particle interaction force corresponding to the particle group by the expression of particle shear stress and inter-particle discrete stress in the inertial action area; An analysis module, configured to analyze the movement speed of the particle group according to the acceleration data, determine the movement speed corresponding to the particle group, and determine the transport distribution data of the particle group in the flow field according to the movement speed; A determination module is used to determine the force of the particle group on the flow field according to the water flow force based on the influence of the particle group on the water flow turbulent kinetic energy and dissipation rate transport equation; convert the force into a volume force to determine the volume force corresponding to the flow field; generate flow field data for the flow field in the second time step based on the volume force to obtain second flow field information; An updating module is used to use the second flow field information as new first flow field information, and to redetermine the transport distribution data of the particle group in the flow field based on the second flow field information.

6. An electronic device, It is characterized in that include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the suspended sediment transport distribution data analysis method according to any one of claims 1 to 4 by executing the computer instructions.

7. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the suspended sediment transport distribution data analysis method according to any one of claims 1 to 4.

Images