Method and device for predicting dynamic behavior information of cavitation bubbles near particle clusters

Abstract

The present invention provides a method and device for predicting the dynamic behavior information of cavitation bubbles near particle clusters, comprising: determining the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when a particle cluster exists around the cavitation bubble according to the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles; determining the velocity distribution of any position point in the fluid field outside the cavitation bubble according to the velocity potential of any position point outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when a particle cluster exists around the cavitation bubble; determining the impulse distribution of any position point in the fluid field outside the cavitation bubble according to the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle, so as to realize accurate analysis of the morphological evolution of the cavitation bubble growth and collapse stages, and facilitate real-time observation and research improvement of material damage.

Description

Technical Field

[0001] The invention relates to the field of fluid mechanics, and in particular to a method and a device for predicting information on the dynamic behavior of cavitation bubbles near a particle cluster. Background Art

[0002] Cavitation refers to the process of formation, development and collapse of cavitation bubbles inside the liquid or at the liquid-solid interface when the local pressure inside the liquid drops below the saturated vapor pressure. During the whole process, cavitation bubbles are rapidly generated, expanded and collapsed, forming shock waves or high-speed jets in the liquid. Cavitation is the phenomenon that shock waves and microjets formed after the collapse of cavitation bubbles damage the walls of nearby fluid machinery. Cavitation is widely present in hydraulic machinery and will have a serious impact on the efficiency and life of related equipment such as water pumps and turbines. At the same time, the river contains relatively more sand, and sediment wear has become one of the important issues that cannot be ignored in the operation and management of hydraulic machinery. Sediment wear will not only cause erosion to the flow-through parts of hydraulic machinery and reduce the service life of hydraulic machinery, but also when cavitation occurs at the same time, more destructive abrasion phenomena will occur.

[0003] In the abrasive erosion phenomenon, particles always exist in the form of particle clusters in sandy water flow. The combined effect of cavitation bubbles and particle clusters can seriously damage the material surface of hydraulic machinery flow channel components, resulting in high material protection costs. Therefore, understanding the interaction between cavitation bubbles and particle clusters is crucial for predicting the abrasive damage mechanism.

[0004] Current research focuses on the interaction between single particles and cavitation bubbles. Through theoretical, experimental, and numerical simulation methods, the evolution of cavitation bubble morphology, jet behavior during the collapse of cavitation bubbles, and particle motion driven by cavitation bubbles are studied. The influence of particle groups on the dynamic behaviors of cavitation bubble growth, collapse, jet flow, and motion is closer to reality and more complex. Summary of the invention

[0005] In view of the problems existing in the prior art, the present invention provides a method and device for predicting the dynamic behavior information of cavitation bubbles near a particle cluster.

[0006] In a first aspect, the present invention provides a method for predicting information on the dynamic behavior of cavitation bubbles near a particle cluster, comprising:

[0007] Acquire a fluid scene graph, wherein the fluid scene graph contains particle clusters and cavitation bubbles, identify the location regions of the particle clusters and cavitation bubbles on the fluid scene graph, and construct the location regions of virtual cavitation bubbles and uniform line sinks based on the location region of each particle on the fluid scene graph;

[0008] Constructing a coordinate system in the fluid scene graph, and determining the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system;

[0009] Determine the position distribution information of the virtual cavitation bubble, the particle and the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, wherein the position distribution information includes the area, radius and spacing;

[0010] According to the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble is determined;

[0011] According to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined;

[0012] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined.

[0013] In one embodiment, determining the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the cavitation bubble includes:

[0014] Determine the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the cavitation bubble;

[0015] Accordingly, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble.

[0016] In one embodiment, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble, including:

[0017] Based on Weiss's theorem, a model is constructed for the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble to generate an additional velocity potential calculation model;

[0018] The coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

[0019] In one embodiment, the method of determining the velocity distribution of any position in the fluid field outside the cavitation bubble according to the velocity potential of any position outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble includes:

[0020] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0021] According to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and the total additional velocity potential caused by the particle cluster at any position point in the fluid field outside the cavitation bubble and each particle, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined by the following calculation formula;

[0022]

[0023] Among them, n x The unit vector representing the x-axis component of the liquid velocity at any point, n y represents the unit vector of the component of the liquid velocity at any position in the y-axis direction, Φ represents the velocity potential at any position in the fluid field outside the cavitation bubble when there are particle clusters around the cavitation bubble; u represents the velocity distribution at any position in the fluid field outside the cavitation bubble;

[0024]

[0025] in, represents the velocity potential of any point in the fluid field outside the cavitation bubble in the infinite flow field, It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle;

[0026]

[0027] in,

[0028]

[0029] Where m represents the source intensity, that is, the volume of fluid released from the point source per unit time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, and π is the pi;

[0030]

[0031]

[0032] in, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, p0 represents the ambient pressure, σ represents the tension coefficient, κ represents the polynomial index, and t represents the time.

[0033] In one embodiment, determining the impulse distribution at any position in the fluid field outside the cavitation bubble based on the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle includes:

[0034] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0035] According to the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined by the following calculation formula;

[0036]

[0037]

[0038]

[0039] Where T represents the cavitation bubble period, r0 represents the initial position coordinate of the cavitation bubble, and r is the coordinate of any position point in the fluid field outside the cavitation bubble and each particle. It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, t represents the time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, ρ represents the fluid density, and m represents the source intensity, that is, the volume of fluid released from the point source per unit time.

[0040] In a second aspect, the present invention provides a device for predicting information on the dynamic behavior of cavitation bubbles near a particle cluster, comprising:

[0041] An identification module is used to obtain a fluid scene graph, wherein the fluid scene graph contains particle clusters and cavitation bubbles, identify the location areas of the particle clusters and cavitation bubbles on the fluid scene graph, and construct the location areas of virtual cavitation bubbles and uniform line sinks based on the location area of ​​each particle on the fluid scene graph;

[0042] A construction module, configured to construct a coordinate system in the fluid scene graph, and determine the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system;

[0043] A determination module, for determining position distribution information between virtual cavitation bubbles, particles and cavitation bubbles according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, wherein the position distribution information includes area, radius and spacing;

[0044] A calculation module, for determining, based on the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there is a particle cluster around the cavitation bubble;

[0045] The first prediction module is used to determine the velocity distribution of any position point in the fluid field outside the cavitation bubble based on the velocity potential of any position point outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble;

[0046] The second prediction module is used to determine the impulse distribution of any position in the fluid field outside the cavitation bubble according to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle.

[0047] In one embodiment, the determining module is specifically used to:

[0048] Determine the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the cavitation bubble;

[0049] Accordingly, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble.

[0050] In one embodiment, the calculation module is specifically used for:

[0051] Based on Weiss's theorem, a model is constructed for the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble to generate an additional velocity potential calculation model;

[0052] The coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

[0053] In a third aspect, the present invention provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of the method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster as described in any one of the first aspects are implemented.

[0054] In a fourth aspect, the present invention provides a non-transitory computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster as described in any one of the first aspects.

[0055] The method and device for predicting the dynamic behavior information of cavitation bubbles near particle clusters provided by the present invention determine the coordinate information of the position area of ​​the particle clusters and the cavitation bubbles on the coordinate system, and calculate and determine the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle based on the coordinate information, thereby determining the velocity distribution and impulse distribution of any position point in the fluid field outside the cavitation bubble based on the additional velocity potential and the velocity potential of any position point outside the cavitation bubble in the infinite flow field, thereby achieving accurate analysis of the morphological evolution of the cavitation bubble growth and collapse stages, and facilitating real-time observation and research improvement of material damage. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0057] Figure 1 It is a schematic flow chart of the method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster provided by the present invention;

[0058] Figure 2 It is a schematic diagram of the layout of three particles and a single cavitation bubble provided by the present invention;

[0059] Figure 3It is a cloud diagram of the velocity distribution field around particles and cavitation bubbles during the cavitation bubble expansion stage provided by the present invention;

[0060] Figure 4 It is a cloud diagram of the velocity distribution field around particles and cavitation bubbles in the cavitation bubble collapse stage provided by the present invention;

[0061] Figure 5 It is the impulse distribution vector diagram near three particles provided by the present invention;

[0062] Figure 6 It is a schematic diagram of the structure of the device for predicting the dynamic behavior information of cavitation bubbles near a particle cluster provided by the present invention;

[0063] Figure 7 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0064] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are 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.

[0065] Combine the following Figures 1 to 7 The invention describes a method, device, electronic device and storage medium for predicting the dynamic behavior information of cavitation bubbles near a particle cluster.

[0066] Figure 1 A schematic diagram showing a flow chart of a method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster provided by the present invention is shown in FIG. Figure 1 , the method comprising:

[0067] 11. Obtain a fluid scene graph, which contains particle clusters and cavitation bubbles, identify the location areas of the particle clusters and cavitation bubbles on the fluid scene graph, and construct the location areas of virtual cavitation bubbles and uniform line sinks based on the location area of ​​each particle on the fluid scene graph;

[0068] 12. Construct a coordinate system in the fluid scene graph, and determine the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system;

[0069] 13. According to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, the position distribution information between the virtual cavitation bubble, the particles and the cavitation bubble is determined, and the position distribution information includes area, radius and spacing;

[0070] 14. According to the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles, determine the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble;

[0071] 15. Determine the velocity distribution of any point in the fluid field outside the cavitation bubble based on the velocity potential of any point outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble;

[0072] 16. Based on the additional velocity potential caused by each particle at any point in the fluid field outside the cavitation bubble and each particle, determine the impulse distribution at any point in the fluid field outside the cavitation bubble.

[0073] With respect to steps 11 to 16, it should be noted that, in the present invention, cavitation refers to the process of formation, development and collapse of cavitation bubbles inside the liquid or at the liquid-solid interface when the local pressure in the liquid is reduced to below the saturated vapor pressure. During the whole process, cavitation bubbles are rapidly generated, expanded and collapsed, forming shock waves or high-speed jets in the liquid. Cavitation is a phenomenon in which shock waves and microjets formed after the collapse of cavitation bubbles damage the walls of fluid machinery nearby. Cavitation is widely present in hydraulic machinery and will have a serious impact on the efficiency and life of related equipment such as water pumps and turbines. At the same time, the sand content of rivers is relatively high, and sediment wear has become one of the important issues that cannot be ignored in the operation and management of hydraulic machinery. Sediment wear will not only cause erosion to the flow-through parts of hydraulic machinery and reduce the service life of hydraulic machinery, but also when cavitation occurs at the same time, more destructive abrasion phenomena will occur.

[0074] In the abrasive erosion phenomenon, particles always exist in the form of particle clusters in sandy water flow. The combined effect of cavitation bubbles and particle clusters can seriously damage the material surface of hydraulic machinery flow channel components, resulting in high material protection costs. Therefore, understanding the interaction between cavitation bubbles and particle clusters is crucial for predicting the abrasive damage mechanism.

[0075] In the present invention, photography technology is used to obtain a fluid scene image, in which particle clusters and cavitation bubbles are contained. To this end, image recognition technology can be used to identify the location area of ​​the particle clusters and cavitation bubbles on the fluid scene image.

[0076] In the present invention, the velocity potential of the fluid near the cavitation bubble when the particle exists can be solved based on the Weiss theorem. At this time, the influence of a single particle on the flow field (i.e., the fluid scene) can be equivalent to a combination of a virtual cavitation bubble and a uniform line sink.

[0077] A particle cluster is an accumulation of multiple particles. In view of the difference between particles and cavitation bubbles (such as size, shape, etc.), particles can be equalized to cavitation bubbles. To this end, a virtual cavitation bubble is constructed based on the position area of ​​each particle on the fluid scene graph, and the position area of ​​the uniform line sink is constructed.

[0078] The following takes three particles near the cavitation bubble as an example to explain the formation of the virtual cavitation bubble and the uniform line sink. Figure 2 As shown in the figure, the center of the three particles is taken as the origin of the coordinate system, L represents the distance between the centers of each particle, and the three larger circles on the periphery represent the three radii, R p1 , R p2 and R p3 The smaller circle at the center represents a particle with a maximum radius of R. max and the cavitation bubble with the initial center position r0=(x0,y0). l1 represents the distance between particle 1 and the center of the cavitation bubble, l2 represents the distance between particle 2 and the center of the cavitation bubble, and l3 represents the distance between particle 3 and the center of the cavitation bubble. The three dotted small circles in the three larger circles represent three virtual cavitation bubbles, and the three dotted lines represent three uniform line sinks. Figure 2 The center is the line segment from the center of the virtual cavitation bubble to the center of the particle.

[0079] In the present invention, a coordinate system is constructed in the fluid scene graph (see Figure 2 ), the center of the cavitation bubble can be taken as the origin of the coordinate system. The coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble are determined on the coordinate system.

[0080] The use of the coordinates of each object on the fluid scene graph can be used to calculate the position distribution information of each object. The position distribution information may be the area, radius, and spacing of the object, but is not limited thereto.

[0081] To this end, the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble is determined based on the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble.

[0082] In the present invention, the position distribution information of virtual cavitation bubbles, particles and cavitation bubbles is mainly reflected in: the regional area of ​​cavitation bubbles, the radius of cavitation bubbles, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble.

[0083] The purpose of the present invention is to obtain the dynamic behavior information of cavitation bubbles near particle clusters. The dynamic behavior information of cavitation bubbles includes: the velocity distribution of any position point in the fluid field outside the cavitation bubble when there are particle clusters around the cavitation bubble, and the impulse distribution of any position point in the fluid field outside the cavitation bubble. In other words, the dynamic behavior information of cavitation bubbles is equivalent to the distribution of velocity and impulse at different positions in the fluid scene diagram.

[0084] In the present invention, since particles always exist in the form of particle clusters in the sand-containing water flow, the presence of particles will also affect the velocity of the surrounding fluid. Therefore, by analyzing and calculating the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles, the velocity potential of any position point in the fluid field outside the cavitation bubble is determined when there are particle clusters around the cavitation bubble.

[0085] The determination of the above “velocity potential” is further explained as follows:

[0086] According to the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0087] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0088] Obtain the velocity potential of any point in the fluid field outside the cavitation bubble in the infinite flow field;

[0089] Based on the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle and the velocity potential at any position in the fluid field outside the cavitation bubble in the infinite flow field, the velocity potential at any position in the fluid field outside the cavitation bubble when there are particle clusters around the cavitation bubble is determined.

[0090] In this regard, it should be noted that, in the present invention, since particles always exist in the form of particle clusters in the sand-laden water flow, it is necessary to determine the effect of each particle on the velocity of the surrounding fluid. Here, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined.

[0091] Then, based on the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined.

[0092] In the present invention, the infinite flow field corresponds to a scenario where there is no boundary around the cavitation bubble (eg, no particles, etc.). In this case, the velocity potential of any position point in the fluid field outside the cavitation bubble in the infinite flow field needs to be calculated.

[0093] Then, based on the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle and the velocity potential at any position in the fluid field outside the cavitation bubble in the infinite flow field, the velocity potential at any position in the fluid field outside the cavitation bubble when there is a particle cluster around the cavitation bubble is determined. For example, the total additional velocity potential and the velocity potential at any position in the fluid field outside the cavitation bubble in the infinite flow field are summed.

[0094] In the present invention, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined based on the velocity potential of any position point outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle. That is, based on the additional velocity potential and the velocity potential, the velocity distribution of any position point in the fluid field outside the cavitation bubble can be obtained according to a preset calculation method.

[0095] In the present invention, the determination process of "velocity distribution" is further explained as follows:

[0096] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0097] According to the velocity potential of any position outside the cavitation bubble in the infinite flow field, and the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble, the velocity distribution of any position in the fluid field outside the cavitation bubble is determined by the following calculation formula;

[0098]

[0099] Among them, n x The unit vector representing the x-axis component of the liquid velocity at any point, n y represents the unit vector of the component of the liquid velocity at any position in the y-axis direction, Φ represents the velocity potential at any position in the fluid field outside the cavitation bubble when there are particle clusters around the cavitation bubble; u represents the velocity distribution at any position in the fluid field outside the cavitation bubble;

[0100]

[0101] in, represents the velocity potential of any point in the fluid field outside the cavitation bubble in the infinite flow field, It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle;

[0102]

[0103] in,

[0104]

[0105] Where m represents the source intensity, that is, the volume of fluid released from the point source per unit time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, and π is the pi;

[0106]

[0107]

[0108] in, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, p0 represents the ambient pressure, σ represents the tension coefficient, k is the polynomial index, and t represents the time.

[0109] In the present invention, based on the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined. Based on the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position point in the fluid field outside the cavitation bubble and each particle is determined. Then, based on the total additional velocity potential caused by the particle cluster at any position point in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position point in the fluid field outside the cavitation bubble is determined.

[0110] In the present invention, by determining the velocity distribution and impulse distribution of any position point in the fluid field outside the cavitation bubble, accurate analysis of the morphological evolution of the cavitation bubble during the growth and collapse stages is achieved.

[0111] Take the fluid scene with three particles as an example, see Figure 3 , Figure 4 and Figure 5 , which are respectively the velocity distribution field cloud diagram of particles and the surrounding areas of cavitation bubbles in the expansion stage of cavitation bubbles, the velocity distribution field cloud diagram of particles and the surrounding areas of cavitation bubbles in the collapse stage of cavitation bubbles, and the impulse distribution vector diagram near the three particles.

[0112] The method for predicting the dynamic behavior information of cavitation bubbles near particle clusters provided by the present invention identifies the position areas of particle clusters and cavitation bubbles on a fluid scene diagram, determines the coordinate information of the position areas of particle clusters and cavitation bubbles on a coordinate system, and calculates and determines the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle based on the coordinate information, thereby determining the velocity distribution and impulse distribution of any position point in the fluid field outside the cavitation bubble based on the additional velocity potential and the velocity potential of any position point outside the cavitation bubble in the infinite flow field, thereby achieving accurate analysis of the morphological evolution of the cavitation bubble growth and collapse stages, and facilitating real-time observation and research improvement of material damage.

[0113] In a further method of the above method, the processing process of determining the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble is explained as follows:

[0114] Based on Weiss's theorem, the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble are modeled to generate an additional velocity potential calculation model.

[0115] The coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

[0116] In this regard, it should be noted that in the present invention, since the number of particles present in the area surrounding the cavitation bubble may be different, it may cause different effects on the velocity caused by any position point in the fluid field, that is, the same calculation method cannot be used to calculate the velocity potential in each scenario. Therefore, it is necessary to perform model analysis on a large amount of historical data and summarize an additional velocity potential calculation model. The model can generate different calculation methods based on different particle conditions. To this end, the coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model, and the coordinates can be calculated according to a reasonable calculation method to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

[0117] The following is an explanation of the fluid scene with three particles, as follows:

[0118] See also Figure 2 , r i1 Represents the position coordinates of the virtual cavitation bubble assumed for particle 1, r p1The center of particle 1 is r p1 to r i1 The position range of the average line sink is as follows:

[0119]

[0120]

[0121] The additional velocity potential caused by particle 1 for:

[0122]

[0123] r i2 Represents the position coordinates of the virtual cavitation bubble assumed for particle 2, r p2 The center of particle 2 is r p2 to r i2 The position range of the average line sink is as follows:

[0124]

[0125]

[0126] The additional velocity potential caused by particle 2 for:

[0127]

[0128]

[0129] r i3 Represents the position coordinates of the virtual cavitation bubble assumed by particle 3, r p3 The center of particle 3 is r p3 to r i3 The position range of the average line sink is as follows:

[0130]

[0131]

[0132] The additional velocity potential caused by particle 3 for:

[0133]

[0134] Based on the velocity potential superposition theorem, the additional velocity potential of the liquid caused by the particle can be calculated from the additional velocity potential caused by particle 1. The additional velocity potential caused by particle 2 and the additional velocity potential caused by particle 3 Superposition gives:

[0135]

[0136] In a further method of the above method, the processing process of determining the impulse distribution at any position in the fluid field outside the cavitation bubble based on the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is explained as follows:

[0137] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0138] According to the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined by the following calculation formula;

[0139]

[0140]

[0141]

[0142] Where T represents the cavitation bubble period, r0 represents the initial position coordinate of the cavitation bubble, and r is the coordinate of any position point in the fluid field outside the cavitation bubble and each particle. It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, t represents the time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, ρ represents the fluid density, and m represents the source intensity, that is, the volume of fluid released from the point source per unit time.

[0143] The present invention identifies the position areas of particle clusters and cavitation bubbles on a fluid scene diagram, determines the coordinate information of the position areas of particle clusters and cavitation bubbles on a coordinate system, and calculates and determines the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle based on the coordinate information, thereby determining the velocity distribution and impulse distribution of any position point in the fluid field outside the cavitation bubble based on the additional velocity potential and the velocity potential of any position point outside the cavitation bubble in an infinite flow field, thereby achieving accurate analysis of the morphological evolution of the cavitation bubble growth and collapse stages, and facilitating real-time observation and research improvement of material damage.

[0144] The following is a description of the device for predicting the dynamic behavior information of cavitation bubbles near particle clusters provided by the present invention. The device for predicting the dynamic behavior information of cavitation bubbles near particle clusters described below and the method for predicting the dynamic behavior information of cavitation bubbles near particle clusters described above can be referenced to each other.

[0145] Figure 6 The schematic diagram of the structure of a device for predicting the dynamic behavior information of cavitation bubbles near a particle cluster provided by the present invention is shown in FIG. Figure 6 The device includes a data identification module 61, a construction module 62, a determination module 63, a first calculation module 64, a second calculation module 65, a first prediction module 66 and a second prediction module 67, wherein:

[0146] The identification module 61 is used to obtain a fluid scene graph, wherein the fluid scene graph contains particle clusters and cavitation bubbles, identify the location areas of the particle clusters and cavitation bubbles on the fluid scene graph, and construct the location areas of virtual cavitation bubbles and uniform line sinks based on the location area of ​​each particle on the fluid scene graph;

[0147] A construction module 62, configured to construct a coordinate system in the fluid scene graph, and determine the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system;

[0148] A determination module 63 is used to determine the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, wherein the position distribution information includes area, radius and spacing;

[0149] A calculation module 64 is used to determine, based on the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there is a particle cluster around the cavitation bubble;

[0150] The first prediction module 65 is used to determine the velocity distribution of any position point in the fluid field outside the cavitation bubble based on the velocity potential of any position point outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble;

[0151] The second prediction module 66 is used to determine the impulse distribution of any position in the fluid field outside the cavitation bubble according to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle.

[0152] In a further device of the above device, the determining module is specifically used for:

[0153] According to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the cavitation bubble, the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble are determined;

[0154] Accordingly, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble.

[0155] In a further device of the above device, the calculation module is specifically used in the process of determining the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle according to the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble, to:

[0156] Based on Weiss's theorem, the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble are modeled to generate an additional velocity potential calculation model.

[0157] The coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

[0158] In a further device of the above device, the first prediction module is specifically used for:

[0159] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0160] According to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and the total additional velocity potential caused by the particle cluster at any position point in the fluid field outside the cavitation bubble and each particle, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined by the following calculation formula;

[0161]

[0162] Among them, n x The unit vector representing the x-axis component of the liquid velocity at any point, n y represents the unit vector of the component of the liquid velocity at any position in the y-axis direction, Φ represents the velocity potential at any position in the fluid field outside the cavitation bubble when there are particle clusters around the cavitation bubble; u represents the velocity distribution at any position in the fluid field outside the cavitation bubble;

[0163]

[0164] in, represents the velocity potential of any point in the fluid field outside the cavitation bubble in the infinite flow field, It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle;

[0165]

[0166] in,

[0167]

[0168] Where m represents the source intensity, that is, the volume of fluid released from the point source per unit time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, and π is the pi;

[0169]

[0170]

[0171] in, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, p0 represents the ambient pressure, σ represents the tension coefficient, κ represents the polynomial index, and t represents the time.

[0172] In a further device of the above device, the second prediction module is specifically used for:

[0173] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined;

[0174] According to the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined, including:

[0175] According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined by the following calculation formula;

[0176]

[0177]

[0178]

[0179] Where T represents the cavitation bubble period, r0 represents the initial position coordinate of the cavitation bubble, and r is the coordinate of any position point in the fluid field outside the cavitation bubble and each particle. It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, t represents the time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, ρ represents the fluid density, and m represents the source intensity, that is, the volume of fluid released from the point source per unit time.

[0180] Since the principle of the device described in the embodiment of the present invention is the same as that of the method described in the above embodiment, a more detailed explanation will not be repeated here.

[0181] The device for predicting the dynamic behavior information of cavitation bubbles near particle clusters provided by the present invention identifies the position areas of particle clusters and cavitation bubbles on a fluid scene diagram, determines the coordinate information of the position areas of particle clusters and cavitation bubbles on a coordinate system, and calculates and determines the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle based on the coordinate information, thereby determining the velocity distribution and impulse distribution of any position point in the fluid field outside the cavitation bubble based on the additional velocity potential and the velocity potential of any position point outside the cavitation bubble in an infinite flow field, thereby achieving accurate analysis of the morphological evolution of the cavitation bubble growth and collapse stages, and facilitating real-time observation and research improvement of material damage.

[0182] Figure 7 An example of a physical structure diagram of an electronic device is shown in FIG. Figure 7As shown, the electronic device may include: a processor 71, a communication interface 72, a memory 73 and a communication bus 74, wherein the processor 71, the communication interface 72 and the memory 73 communicate with each other through the communication bus 74. The processor 71 may call a computer program in the memory 73 to execute the steps of a method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster, for example, including: obtaining a fluid scene graph, the fluid scene graph contains particle clusters and cavitation bubbles, identifying the position areas of the particle clusters and cavitation bubbles on the fluid scene graph, and constructing the position areas of virtual cavitation bubbles and uniform line sinks on the fluid scene graph based on the position areas of each particle; constructing a coordinate system in the fluid scene graph, and determining the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble based on the coordinate system; determining the virtual cavitation bubble, the particle and the uniform line sink according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble; and determining the virtual cavitation bubble, the particle and the uniform line sink according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble. Position distribution information between cavitation bubbles, including area, radius and spacing; when there are particle clusters around the cavitation bubbles, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles; when there are particle clusters around the cavitation bubbles, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the velocity potential at any position outside the cavitation bubble in the infinite flow field, and when there are particle clusters around the cavitation bubbles, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined, the velocity distribution of any position in the fluid field outside the cavitation bubble is determined; when there is an additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution of any position in the fluid field outside the cavitation bubble is determined.

[0183] In addition, the logic instructions in the above-mentioned memory 73 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or the part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0184] On the other hand, the present invention also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium, and the computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the steps of a method for predicting the dynamic behavior information of cavitation bubbles near particle clusters, for example, including: obtaining a fluid scene graph, the fluid scene graph contains particle clusters and cavitation bubbles, identifying the position areas of the particle clusters and cavitation bubbles on the fluid scene graph, and constructing the position areas of virtual cavitation bubbles and uniform line sinks on the fluid scene graph based on the position areas of each particle; constructing a coordinate system in the fluid scene graph, and determining the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system; according to the coordinates of the center of the virtual cavitation bubble, the coordinates of the uniform line sink The position distribution information of virtual cavitation bubbles, particles and cavitation bubbles is determined by the range, coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, and the position distribution information includes area, radius and spacing; when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined according to the position distribution information of virtual cavitation bubbles, particles and cavitation bubbles; when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined; when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble is determined according to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined; according to the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle, the impulse distribution of any position point in the fluid field outside the cavitation bubble is determined.

[0185] On the other hand, an embodiment of the present invention further provides a processor-readable storage medium, wherein the processor-readable storage medium stores a computer program, wherein the computer program is used to enable the processor to execute the steps of a method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster, for example, including: obtaining a fluid scene graph, wherein the fluid scene graph contains particle clusters and cavitation bubbles, identifying the position areas of the particle clusters and cavitation bubbles on the fluid scene graph, and constructing the position areas of virtual cavitation bubbles and uniform line sinks on the fluid scene graph based on the position areas of each particle; constructing a coordinate system in the fluid scene graph, and determining the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system; and The coordinates of the center of the cavitation bubble are used to determine the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble, and the position distribution information includes the area, radius and spacing; when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined according to the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble; when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined; when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble is determined according to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and when there are particle clusters around the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle is determined, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined; according to the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle, the impulse distribution of any position point in the fluid field outside the cavitation bubble is determined.

[0186] The processor-readable storage medium can be any available medium or data storage device that can be accessed by the processor, including but not limited to magnetic storage (such as floppy disks, hard disks, magnetic tapes, magneto-optical disks (MO), etc.), optical storage (such as CD, DVD, BD, HVD, etc.), and semiconductor storage (such as ROM, EPROM, EEPROM, non-volatile memory (NANDFLASH), solid-state drive (SSD)), etc.

[0187] The device embodiments described above are merely illustrative, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without paying creative labor.

[0188] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0189] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster, characterized in that: include: Acquire a fluid scene graph, wherein the fluid scene graph contains particle clusters and cavitation bubbles, identify the location regions of the particle clusters and cavitation bubbles on the fluid scene graph, and construct the location regions of virtual cavitation bubbles and uniform line sinks based on the location region of each particle on the fluid scene graph; Constructing a coordinate system in the fluid scene graph, and determining the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system; Determine the position distribution information of the virtual cavitation bubble, the particle and the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, wherein the position distribution information includes the area, radius and spacing; According to the position distribution information between virtual cavitation bubbles, particles and cavitation bubbles, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble is determined; According to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined; According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined; The method of determining the impulse distribution at any position in the fluid field outside the cavitation bubble according to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle comprises: According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined; According to the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined by the following calculation formula; Where T represents the cavitation bubble period, r0 represents the initial position coordinate of the cavitation bubble, and r is the coordinate of any position point in the fluid field outside the cavitation bubble and each particle. It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, t represents the time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, ρ represents the fluid density, and m represents the source intensity, that is, the volume of fluid released from the point source per unit time.

2. The method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster according to claim 1, characterized in that: The determining of the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the cavitation bubble comprises: Determine the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the cavitation bubble; Accordingly, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble.

3. The method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster according to claim 2, characterized in that: Determining the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble, comprises: Based on Weiss's theorem, a model is constructed for the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble to generate an additional velocity potential calculation model; The coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

4. The method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster according to claim 3, characterized in that: The method of determining the velocity distribution of any position in the fluid field outside the cavitation bubble according to the velocity potential of any position outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble includes: According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined; According to the velocity potential of any position point outside the cavitation bubble in the infinite flow field, and the total additional velocity potential caused by the particle cluster at any position point in the fluid field outside the cavitation bubble and each particle, the velocity distribution of any position point in the fluid field outside the cavitation bubble is determined by the following calculation formula; Among them, n x The unit vector representing the x-axis component of the liquid velocity at any point, n y represents the unit vector of the component of the liquid velocity at any position in the y-axis direction, Φ represents the velocity potential at any position in the fluid field outside the cavitation bubble when there are particle clusters around the cavitation bubble; u represents the velocity distribution at any position in the fluid field outside the cavitation bubble; in, represents the velocity potential of any point in the fluid field outside the cavitation bubble in the infinite flow field, It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle; in, Where m represents the source intensity, that is, the volume of fluid released from the point source per unit time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, and π is the pi; in, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, p0 represents the ambient pressure, σ represents the tension coefficient, κ represents the polynomial index, and t represents the time.

5. A device for predicting the dynamic behavior information of cavitation bubbles near a particle cluster, characterized in that: include: An identification module is used to obtain a fluid scene graph, wherein the fluid scene graph contains particle clusters and cavitation bubbles, identify the location areas of the particle clusters and cavitation bubbles on the fluid scene graph, and construct the location areas of virtual cavitation bubbles and uniform line sinks based on the location area of ​​each particle on the fluid scene graph; A construction module, configured to construct a coordinate system in the fluid scene graph, and determine the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the center of the cavitation bubble based on the coordinate system; A determination module, for determining position distribution information between virtual cavitation bubbles, particles and cavitation bubbles according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle and the coordinates of the center of the cavitation bubble, wherein the position distribution information includes area, radius and spacing; A calculation module, for determining, based on the position distribution information between the virtual cavitation bubble, the particle and the cavitation bubble, the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there is a particle cluster around the cavitation bubble; The first prediction module is used to determine the velocity distribution of any position point in the fluid field outside the cavitation bubble based on the velocity potential of any position point outside the cavitation bubble in the infinite flow field and the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle when there are particle clusters around the cavitation bubble; The second prediction module is used to determine the impulse distribution of any position in the fluid field outside the cavitation bubble according to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle; The method of determining the impulse distribution at any position in the fluid field outside the cavitation bubble according to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle comprises: According to the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle, the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle is determined; According to the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, the impulse distribution at any position in the fluid field outside the cavitation bubble is determined by the following calculation formula; Where T represents the cavitation bubble period, r0 represents the initial position coordinate of the cavitation bubble, and r is the coordinate of any position point in the fluid field outside the cavitation bubble and each particle. It represents the total additional velocity potential caused by the particle cluster at any position in the fluid field outside the cavitation bubble and each particle, t represents the time, R represents the instantaneous radius of the cavitation bubble, It represents the first-order derivative of the instantaneous radius of the cavitation bubble with respect to time, that is, the radial velocity at the wall of the cavitation bubble, ρ represents the fluid density, and m represents the source intensity, that is, the volume of fluid released from the point source per unit time.

6. The device for predicting the dynamic behavior information of cavitation bubbles near particle clusters according to claim 5, characterized in that: The determination module is specifically used for: Determine the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble according to the coordinates of the center of the virtual cavitation bubble, the coordinate range of the uniform line sink, the coordinates of the center of each particle, and the coordinates of the cavitation bubble; Accordingly, the additional velocity potential caused by each particle at any position in the fluid field outside the cavitation bubble and each particle is determined based on the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble.

7. The device for predicting the dynamic behavior information of cavitation bubbles near particle clusters according to claim 6, characterized in that: The calculation module is specifically used for: Based on Weiss's theorem, a model is constructed for the area of ​​the cavitation bubble, the radius of the cavitation bubble, the distance between the centers of each particle in the particle cluster, the radius of each particle, and the distance between the center of each particle and the center of the cavitation bubble to generate an additional velocity potential calculation model; The coordinates of any position point in the fluid field outside the cavitation bubble and each particle are input into the additional velocity calculation model to obtain the additional velocity potential caused by each particle at any position point in the fluid field outside the cavitation bubble and each particle.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps of the method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster as described in any one of claims 1 to 4 are implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for predicting the dynamic behavior information of cavitation bubbles near a particle cluster according to any one of claims 1 to 4 are implemented.

Images