Method and device for determining micro-particle distribution information, electronic equipment and storage medium

Abstract

A method, apparatus, electronic device, and storage medium for determining particle distribution information are disclosed. The method includes: acquiring the particle size of a micro-group to be tested, wherein the micro-group to be tested is pneumatically broken into broken particles having various particle sizes; acquiring the mass fraction of broken particles of each particle size; determining the total number of broken particles based on the micro-group particle size and the mass fraction of broken particles of each particle size; acquiring the dispersion angle of each broken particle, wherein the dispersion angle is the angle between the direction of motion velocity of the micro-group to be tested and the direction of motion velocity of the broken particles; determining a first proportion of broken particles located at each dispersion angle; and determining the number of particles at each dispersion angle based on the total number of particles and the first proportion. This application can accurately simulate the distribution information of each particle formed after pneumatic breakage.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of data processing, and in particular to a method and device for determining particle distribution information, an electronic device, and a storage medium. BACKGROUND

[0002] Aerodynamic breakup refers to the deformation of a micro-cluster moving in a gas flow if there is a slip velocity between the micro-cluster and the gas flow. When the relative velocity between the micro-cluster and the gas flow reaches a critical value, the micro-cluster will break up due to the viscous force of the deformed micro-cluster being greater than its surface tension, generating smaller particles.

[0003] Exploring the distribution morphology of the broken particles can have a positive effect on natural phenomena and engineering practical problems in the process of aerodynamic dispersion. For example, in the process of fuel air explosive (FAE) blasting, predicting the concentration and distribution of fuel particles generated after the aerodynamic breakup of fuel micro-clusters can guide the secondary ignition scheme and optimize the cloud explosion damage effect.

[0004] How to accurately simulate the distribution information of each particle formed after aerodynamic breakup is a problem to be solved. SUMMARY

[0005] In view of the above, the embodiments of the present application provide a method and device for determining particle distribution information, an electronic device, and a storage medium, which can accurately simulate the distribution information of each particle formed after aerodynamic breakup.

[0006] An embodiment of the present application provides a method for determining particle distribution information, comprising: obtaining a micro-cluster particle size of a micro-cluster to be measured, the micro-cluster to be measured forming broken particles with each particle size after aerodynamic breakup; obtaining a mass fraction of the broken particles with each particle size; determining a total particle number of the broken particles based on the micro-cluster particle size and the mass fraction of the broken particles with each particle size; obtaining a dispersion angle of each broken particle, the dispersion angle being an included angle between a direction of a movement speed of the micro-cluster to be measured and a direction of a movement speed of the broken particle; determining a first number ratio of the broken particles located at each dispersion angle; and determining a dispersion angle particle number located at each dispersion angle based on the total particle number and the first number ratio.

[0007] With this technical solution, the number of each broken particle formed after the breakup of the micro-cluster particle size at each dispersion angle can be obtained based on the micro-cluster particle size, the mass fraction of the broken particles with each particle size, and the number of each broken particle at each dispersion angle, thereby accurately simulating the concentration distribution of each broken particle at each dispersion angle.

[0008] In some embodiments, after the step of determining the number of particles at the scattering angle based on the total number of particles and the first number ratio, the method further comprises: obtaining a movement speed of the to-be-tested micro-cluster; obtaining an included angle of an asymptote line, the included angle of the asymptote line being established based on a movement range of the broken particle; obtaining a maximum movement speed of the broken particle at the scattering angle based on the movement speed of the to-be-tested micro-cluster, the included angle of the asymptote line, and the scattering angle of the broken particle; and determining a speed of the broken particle at the scattering angle based on the maximum movement speed.

[0009] By using the technical scheme, the initial speed of the broken particle at the scattering angle after breaking can be predicted, so that the subsequent diffusion evolution process of the broken particle can be obtained in combination with the initial speed.

[0010] In some embodiments, the maximum movement speed of the broken particle at the scattering angle is obtained based on the movement speed of the to-be-tested micro-cluster, the included angle of the asymptote line, and the scattering angle of the broken particle, and the maximum movement speed of the broken particle at the scattering angle is obtained based on a first preset formula, the first preset formula being:

[0011]

[0012] wherein, the V0 is the movement speed of the to-be-tested micro-cluster, the is a half of the included angle of the asymptote line, the θ is the scattering angle of the broken particle, and the Vmax is the maximum movement speed. θ

[0013] In some embodiments, the speed of the broken particle at the scattering angle is determined based on the maximum movement speed, and the speed of the broken particle at the scattering angle is determined by: obtaining a plurality of speed coefficients; obtaining distribution information of the plurality of speed coefficients in the broken particles, the distribution information representing a second number ratio of the number of the broken particles with the speed coefficient to the total number of particles; assigning the speed coefficient to the broken particle at the scattering angle according to the second number ratio; and multiplying the maximum movement speed by the assigned speed coefficient to obtain the speed of the broken particle at the scattering angle.

[0014] In some embodiments, the speed of the broken particle at the scattering angle is determined based on the maximum movement speed, and the speed of the broken particle at the scattering angle is determined by: obtaining a plurality of speed coefficients; obtaining distribution information of the plurality of speed coefficients in the broken particles, the distribution information representing a second number ratio of the number of the broken particles with the speed coefficient to the total number of particles; assigning the speed coefficient to the broken particle at the scattering angle according to the second number ratio; and multiplying the maximum movement speed by the assigned speed coefficient to obtain the speed of the broken particle at the scattering angle.

[0015] By using the technical scheme, the speed coefficient is multiplied by the speed of the broken particle, so that the speeds of the broken particles diffused from the same scattering angle can be accurately obtained, and the distribution of the broken particles diffused from the same scattering angle can be obtained.

[0016] ​In some embodiments, determining the total particle number of the broken particles based on the micro-cluster particle size and the mass fraction of the broken particles of each particle size comprises: determining the number of broken particles of each particle size based on a second preset formula, wherein the second preset formula is:

[0017]

[0018] wherein, i is the identification of the particle size, d i is the particle size of the broken particles of the identification i, D is the micro-cluster particle size of the micro-cluster to be measured, y i is the mass fraction of the broken particles of the particle size d i , and n i is the number of the broken particles of the particle size d i .

[0019] The total particle number is obtained based on the number of the broken particles of each particle size.

[0020] In some embodiments, obtaining the mass fraction of the broken particles of each particle size comprises: obtaining the ratio of the mass median diameter to the Sauter mean diameter of each broken particle; determining the size distribution index of the broken particles based on the ratio of the mass median diameter to the Sauter mean diameter; determining the cumulative mass fraction of the broken particles of each particle size based on the size distribution index; and obtaining the mass fraction of the broken particles of each particle size based on the cumulative mass fraction.

[0021] In some embodiments, after determining the number of particles at each dispersion angle based on the total particle number and the first number ratio, the method further comprises: obtaining the drag force applied to the broken particles; establishing the Newton motion equation of the broken particles based on the mass of the broken particles and the drag force; and obtaining the motion trajectory of the broken particles changing with time based on the Newton motion equation.

[0022] By using the technical scheme, the diffusion evolution process of the broken particles changing with time can be predicted.

[0023] An embodiment of the present application further provides a device for determining particle distribution information, comprising:

[0024] An information obtaining module is configured to obtain the micro-cluster particle size of a micro-cluster to be measured, and obtain the mass fraction of broken particles of each particle size, wherein the broken particles are formed after the micro-cluster to be measured is broken by gas.

[0025] A total number calculating module is configured to determine the total particle number of the broken particles based on the micro-cluster particle size and the mass fraction of the broken particles of each particle size.

[0026] The proportion determining module obtains a dispersion angle of each broken particle, the dispersion angle being an included angle between a moving speed direction of the to-be-tested micro-cluster and a moving speed direction of the broken particle; and determines a first number proportion of the broken particles located at each dispersion angle;

[0027] The concentration determining module determines a number of dispersion angle particles located at each dispersion angle based on the total particle number and the first number proportion.

[0028] An embodiment of the present application further provides an electronic device, which comprises a processor and a memory, the memory being used for storing instructions, and the processor being used for calling the instructions in the memory, so that the electronic device executes the method for determining particle distribution information.

[0029] An embodiment of the present application further provides a computer readable storage medium, which stores computer instructions, and when the computer instructions are executed on an electronic device, the electronic device executes the method for determining particle distribution information. BRIEF DESCRIPTION OF DRAWINGS

[0030] Figure 1 FIG. 1 is a step flow chart of the method for determining particle distribution information provided by an embodiment of the present application;

[0031] Figure 2 FIG. 2 is a sub-step flow chart of step 102 provided by an embodiment of the present application;

[0032] Figure 3 FIG. 3 is a relationship diagram of broken form, Weber number and Onegzog number provided by an embodiment of the present application;

[0033] Figure 4 FIG. 4 is a relationship diagram of size distribution index of broken particles and target diameter ratio provided by an embodiment of the present application;

[0034] Figure 5 FIG. 5 is a relationship diagram of cumulative mass fraction of broken particles and particle size provided by an embodiment of the present application;

[0035] Figure 6 FIG. 6 is a diagram showing that mass fraction of broken particles and number of broken particles change with particle size provided by an embodiment of the present application;

[0036] Figure 7 FIG. 7 is a diagram showing a scenario of aerodynamic breaking of a to-be-tested micro-cluster provided by an embodiment of the present application;

[0037] Figure 8 FIG. 8 is a concentration distribution model diagram of broken particles at each dispersion angle provided by an embodiment of the present application;

[0038] Figure 9 is a step flow chart for determining the speed of the broken particles according to an embodiment of the present application;

[0039] Figure 10 is a schematic diagram of a vector triangle according to an embodiment of the present application;

[0040] Figure 11 is a sub-step flow chart of step 904 according to an embodiment of the present application;

[0041] Figure 12 is a schematic diagram of distribution information of the speed coefficient in each broken particle according to an embodiment of the present application;

[0042] Figure 13 is a flow chart for determining the motion trajectory of each broken particle evolving over time according to an embodiment of the present application;

[0043] Figure 14 is a structural schematic diagram of a device for determining the particle distribution information according to an embodiment of the present application;

[0044] Figure 15 is a schematic diagram of an electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0045] In order to more clearly understand the above objectives, features and advantages of the present application, the present application will be described in detail below with reference to the accompanying drawings and specific embodiments. It should be noted that the embodiments of the present application and the features in the embodiments can be combined with each other without conflict.

[0046] In the following description, a large number of specific details are set forth in order to facilitate a thorough understanding of the present application. The described embodiments are merely some of the embodiments of the present application, but not all the embodiments.

[0047] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the present application.

[0048] Further, it should be noted that herein, the terms “comprising”, “including”, or any other variant thereof are intended to cover non-exclusive inclusion, so that a process, method, article or apparatus including a series of elements not only includes those elements, but also includes other elements not explicitly listed, or further includes elements inherent in such a process, method, article or apparatus. Without more limitations, the element defined by the phrase “including a…” does not exclude the presence of additional identical elements in the process, method, article or apparatus including the element.

[0049] In the present application, “at least one” refers to one or more, and “multiple” refers to two or more than two. “And / or” describes the association between the associated objects, which means that there can be three relationships, for example, A and / or B can represent: A exists alone, A and B exist together, and B exists alone, where A and B can be singular or plural. The terms “first”, “second”, “third”, “fourth” and the like (if any) in the specification and claims of the present application and the drawings are used to distinguish similar objects, and are not used to describe a specific order or sequence.

[0050] In the embodiments of the present application, the words such as “exemplary” or “for example” are used to mean an example, illustration, or description. Any embodiment or design scheme described as “exemplary” or “for example” in the embodiments of the present application should not be interpreted as more preferred or more advantageous than other embodiments or design schemes. Rather, the words such as “exemplary” or “for example” are intended to present the relevant concept in a specific manner.

[0051] In the process of fuel air explosive (FAE) explosion, FAE will rapidly disperse the solid-liquid mixed cloud fuel in a large space under the explosion driving action of the central explosive, so as to form a multi-phase cloud area with uniform dispersion and relatively stable fuel dispersion, and then form a cloud explosion through secondary ignition, causing large-scale damage.

[0052] If the characteristic time of the stable multi-phase fuel cloud area, the characteristic configuration of the cloud area, and the spatial distribution of the fuel concentration can be simulated, it is very important for formulating a secondary ignition scheme, optimizing and adjusting the structure of the warhead, predicting and optimizing the cloud explosion damage effect.

[0053] However, the above fuel explosion process is located in a transient multiphase flow field, and the diagnostic means of the transient multiphase flow field is lacking, which makes it difficult to capture the micro-evolution process of the broken fuel particles, that is, it is difficult to predict the distribution information of the fuel particles.

[0054] Therefore, how to accurately simulate the distribution information of each particle formed after the aerodynamic breakup is a problem to be solved.

[0055] In view of this, the present application provides a particle distribution information determination method, which can be applied to the above-mentioned fuel explosion application scene, and is also suitable for describing all natural phenomena and engineering practical problems involving the process of liquid droplet aerodynamic dispersion, such as engine fuel droplet atomization, molten droplet atomization and granulation, spray drying, and fire extinguishing agent throwing.

[0056] The method for determining particle distribution information provided in this application includes: obtaining the particle size of a micro-particle to be tested, wherein the micro-particle to be tested is pneumatically broken into broken particles with various particle sizes, and obtaining the mass fraction of broken particles of each particle size; determining the total number of broken particles based on the particle size of the micro-particle and the mass fraction of broken particles of each particle size; obtaining the dispersion angle of each broken particle, wherein the dispersion angle is the angle between the direction of motion velocity of the micro-particle to be tested and the direction of motion velocity of the broken particles; determining a first proportion of the number of broken particles located at each dispersion angle; and determining the number of dispersion angle particles located at each dispersion angle based on the total number of particles and the first proportion.

[0057] The embodiments of this application can obtain the number of each broken particle at each dispersion angle after the micro-aggregate is broken based on the particle size of the micro-aggregate and the mass fraction of each particle size of the broken particles, thereby accurately simulating the concentration distribution of each broken particle at each dispersion angle.

[0058] The method for determining particle distribution information in this application can be applied to one or more electronic devices. The electronic device is a device capable of automatically performing numerical calculations and / or information processing according to pre-set or stored instructions. Its hardware includes, but is not limited to, processors, microprogrammed control units (MCUs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), digital signal processors (DSPs), embedded devices, etc. The electronic device can be a portable electronic device (such as a mobile phone or tablet), a personal computer, a server, etc.

[0059] Figure 1 This is a flowchart illustrating one embodiment of the method for determining particle distribution information according to this application. The order of the steps in the flowchart can be changed, and some steps can be omitted, depending on different requirements.

[0060] See Figure 1 As shown, the method for determining the particle distribution information may include the following steps.

[0061] Step 101: Obtain the particle size of the micro-group to be tested.

[0062] The micro-particles to be measured can be liquid micro-particles, such as fuel micro-particles, fire extinguishing material micro-particles, etc., and this application does not limit this. The particle size is the scale of the space occupied by the micro-particle.

[0063] The to-be-tested micro-cluster can form broken micro-particles with different micro-particle sizes after pneumatic breaking. For example, the fuel micro-cluster breaks to form fuel micro-particles under the action of air flow.

[0064] The following takes the to-be-tested micro-cluster as a fuel micro-cluster as an example. When the fuel air explosive (FAE) explodes, the fuel micro-cluster breaks to form fuel micro-particles:

[0065] The fuel explosion is divided into a near-field stage and a far-field stage. The near-field stage generally occurs within a time scale of O(100)-(101) after initiation. In this stage, the fuel cloud ring belt expands to a space scale of O(101) times the initial fuel shell outer diameter. The fuel shell will accelerate and expand under the driving of the central explosion load and the detonation product flow field, forming a jet structure and decomposing into fuel micro-clusters. In the far-field stage, the fuel cloud area is separated from the influence of the central detonation product and shock wave flow field. A large number of non-uniformly distributed high-speed fuel micro-clusters rapidly disperse in an approximately static flow field, and at the same time, pneumatic breaking occurs to form fuel micro-particles.

[0066] The breaking effect of fuel particles is generally described by a mass shedding model. According to the shedding model proposed by Engel, due to the convective shear effect caused by the rapid flow of gas, small droplets are generated by the shedding of the surface of the liquid droplet. The shedding rate is represented by the following formula:

[0067]

[0068] wherein, ρ is the density of air, ρ L ini is the density of fuel, μ is the air viscosity coefficient, μ L is the fuel viscosity coefficient, u is the air velocity, u L is the fuel velocity, and l is the average radius of the fuel droplet.

[0069] Step 102, obtaining the mass fraction of the broken micro-particles with different micro-particle sizes.

[0070] The mass fraction of the broken micro-particles with a certain micro-particle size is the ratio of the total mass of the broken micro-particles with the micro-particle size to the total mass of the broken micro-particles formed after the pneumatic breaking of the to-be-tested micro-cluster.

[0071] In some embodiments, referring to FIG. 1, step 102 can include: Figure 2

[0072] Step 1021, obtaining the ratio of the mass median diameter of the broken micro-particles to the Sauter mean diameter.

[0073] ​In a series of crushed particles, when the total mass of particles of various sizes smaller than a certain aerodynamic diameter accounts for 50% of the total mass of all crushed particles, this diameter is called the median mass diameter. In other words, in a series of crushed particles, half of the particles have a diameter smaller than the median mass diameter, and the other half have a diameter larger than the median mass diameter. The median mass diameter can be denoted as MMD.

[0074] The Sauter Mean Diameter (SMD) is denoted as D. 32 Its physical meaning can be characterized as the ratio of the total volume of each broken particle to the total area of ​​each broken particle.

[0075] In some embodiments, step 1021 may include: obtaining the Weber number (we) and Ohnesorge number (oh) of the micro-group to be tested; and determining the ratio of the mass median diameter of each broken particle to the Sotter mean diameter (hereinafter referred to as the target diameter ratio) based on the Weber number and the Ohnesorge number.

[0076] Among them, the Weber number is the ratio of the inertial force to the surface tension of the micro-element under test; the Onezog number is the ratio of the viscous force to the surface tension of the micro-element under test.

[0077] Specifically, determining the ratio of the mass median diameter to the Sotter mean diameter of each fragmented particle based on the Weber number and the Onezoglu number may include: determining the interval in which the Weber number and the Onezoglu number of the micro-particle to be measured are located; and obtaining the target diameter ratio that matches the interval in the correspondence between each interval and the diameter ratio.

[0078] For example, refer to Figure 3 As shown, different intervals of Weber number and Onezoglu number correspond to different fragmentation forms, and different fragmentation forms correspond to different diameter ratios. The fragmentation form of the micro-element to be measured can be obtained from the intervals where the Weber number and Onezoglu number of the micro-element to be measured are located. Then, based on the correspondence between the fragmentation form and the diameter ratio, the target diameter ratio of the micro-element to be measured can be obtained.

[0079] The crushing methods may include, but are not limited to, bag breakup, shear-stripping breakup, bag-stamen breakup, and oscillatory breakup.

[0080] For example, if the breakage mode of the micro-group to be tested is determined to be shear breakage based on the Weber number and Onezoglu number, the corresponding diameter ratio is approximately equal to 1.2.

[0081] In some other embodiments, the broken form of the micro-aggregates to be measured can be directly obtained, and then the target diameter ratio matching the broken form of the micro-aggregates to be measured is obtained.

[0082] For example, the broken form of the micro-aggregates to be measured can be stored in the electronic device as a shear broken form, so that based on the corresponding relationship between the shear broken form and the diameter ratio, the target diameter ratio matching the shear broken form is obtained as 1.2.

[0083] In step 1022, the size distribution index of the broken particles is determined based on the ratio of the mass median diameter to the Sauter mean diameter.

[0084] The relationship between the diameter ratio and the size distribution index of the broken particles is shown as follows:

[0085]

[0086] wherein MMD is the mass median diameter of each broken particle, SMD is the Sauter mean diameter of each broken particle, q is the size distribution index of the broken particles, and G represents the gamma function.

[0087] For example, according to the diameter ratio and the size distribution index of the broken particles, the relationship between the size distribution index of the broken particles and the target diameter ratio can be referred to as shown in FIG. 1. Figure 4 In the case where the broken form of the micro-aggregates to be measured is a shear broken form, the MMD / SMD of the broken particles is approximately equal to 1.2, that is, the target size ratio can be regarded as 1.2, and it can be known from FIG. 1 that the size distribution index of the broken particles is 3.0. Figure 4

[0088] In step 1023, the cumulative mass fraction of the broken particles of each particle size is determined based on the size distribution index.

[0089] The cumulative mass fraction of the broken particles of a certain particle size is the ratio of the total mass of the broken particles not smaller than the particle size to the total mass of all the broken particles.

[0090] The relationship between the particle size and the cumulative mass fraction of the broken particles can be described by the following formula:

[0091]

[0092] wherein d is the particle size, SMD is the Sauter mean diameter of each broken particle, q is the size distribution index of the broken particles, and Y d is the ratio of the total mass of the broken particles with the particle size larger than d to the total mass of all the broken particles, that is, the cumulative mass fraction corresponding to the particle size d is Y d .

[0093] ​For example, in shear crushing mode, SMD is approximately equal to 0.90d0, where d0 is the particle size of the micro-particles to be measured, thus yielding the SMD value. With an SMD of 2.0, a particle size range of 0.01 mm to 0.4 mm, and a micro-particle size D of 2 mm, the cumulative mass fraction Y of the crushed particles is... d The relationship with particle size d can be found by referring to Figure 5 As shown.

[0094] Step 1024: Based on the cumulative mass fraction, obtain the mass fraction of broken particles of each particle size.

[0095] The mass fraction of broken particles with a certain particle size is the ratio of the total mass of all broken particles with that particle size to the total mass of all broken particles.

[0096] For example, let y i For particles with a diameter of d i The mass fraction of broken particles, y i =Y i -Y j , where Y i For particles with a diameter of d i The cumulative mass fraction of broken particles, Y j For particles with a diameter of d j The cumulative mass fraction of broken particles, d j Under ideal conditions, it approaches d infinitely. i In practical applications, d can be set according to requirements. i Value. For example, the mass fraction of broken particles with a particle size of 0.1 can be calculated using the following formula:

[0097] y 0.1 =Y 0.1 -Y 0.99 .

[0098] For example, by Figure 5 The cumulative mass fraction Y shown d The relationship with the particle size d can be obtained Figure 6 The curve 601 shows the relationship between the mass fraction of broken particles and the particle size.

[0099] Step 103: Determine the total number of broken particles based on the particle size of the micro-clusters and the mass fraction of broken particles of each particle size.

[0100] In some embodiments, the number of broken particles with a certain particle size can first be obtained based on the particle size of the micro-clusters and the mass fraction of broken particles with a certain particle size. Then, the number of broken particles corresponding to each particle size can be calculated. Finally, the number of broken particles corresponding to each particle size is summed to obtain the total number of particles.

[0101] Furthermore, based on the particle size of the micro-clusters and the mass fraction of broken particles with a certain particle size, the number of broken particles with that particle size can be obtained, which may include:

[0102] The particle size is determined to be d based on the second preset formula. i The number of broken particles, the second preset formula is:

[0103]

[0104] Wherein, i is the identifier of the particle size, and d i Let i be the particle size, D be the particle size of the micro-group to be measured, and y be the particle size of the micro-group to be measured. i For particles with a diameter of d i The mass fraction of broken particles, n i For particles with a diameter of d i The number of broken particles.

[0105] The above formulas can be used to calculate the number of broken particles with different particle sizes formed by the breakage of the micro-group to be tested. Then, the total number of broken particles formed by the breakage of the micro-group to be tested can be obtained by summing them up.

[0106] For example, Figure 6 Curve 602 illustrates the relationship between particle size and particle number.

[0107] Step 104: Obtain the dispersion angle of each broken particle.

[0108] The dispersion angle is the angle between the direction of the velocity of the micro-group to be measured and the direction of the velocity of the broken particles.

[0109] After the micro-group to be tested is broken, each broken particle will peel off at a different angle, as shown in the reference. Figure 7 As shown, in the experiment, a leftward airflow with a certain velocity was applied to the micro-particle to be tested, causing the micro-particle to break up. The broken particles moved to the right relative to the airflow.

[0110] Step 105: Determine the first proportion of broken particles located at each dispersion angle.

[0111] In some embodiments, step 105 may include: obtaining the number of broken particles at each dispersion angle after the micro-group to be tested undergoes aerodynamic breakage, obtaining a concentration distribution model, and obtaining a first quantity ratio based on the concentration distribution model.

[0112] Users can obtain experimental data based on experiments and fit the data to obtain a concentration distribution model, which can then be configured in the electronic device. For example, refer to... Figure 7 As shown, the microparticles to be tested are placed in an airflow with a preset flow rate to simulate the aerodynamic breakup of the microparticles. The number of broken particles located at each dispersion angle after aerodynamic breakup is counted to obtain the concentration distribution model.

[0113] Depend on Figure 7 As shown, the broken particles exhibit a normal distribution across the dispersion angles. The concentration of broken particles is high near the axis of the sample particle, and the axis direction corresponds to the velocity direction of the sample particle; therefore, the axis can be set as the 0° line. Based on the statistically significant number of broken particles at each dispersion angle, the mean and standard deviation of the normal distribution function can be obtained, thus yielding... Figure 8 The fitted concentration distribution model is shown. Wherein, Figure 8 The horizontal axis represents the dispersion angle θ, and the vertical axis represents the first quantity proportion.

[0114] Step 106: Based on the total number of particles and the first quantity ratio, determine the number of dispersion angle particles located at each dispersion angle.

[0115] After the micro-particles under test undergo aerodynamic breakup, multiple broken particles may disperse from the same dispersion angle. The number of particles at the dispersion angle is the total number of broken particles located at that dispersion angle.

[0116] Multiplying the total number of particles by the first quantity ratio yields the number of particles at the dispersion angle.

[0117] In some embodiments, after step 106, the method may further include: determining the velocity of the broken particles at each dispersion angle based on the number of particles at each dispersion angle; and determining the distribution information of the broken particles at each dispersion angle based on the velocity of the broken particles at each dispersion angle.

[0118] Among them, you can refer to Figure 9 As shown, determining the velocity of the broken particles at each dispersion angle based on the number of particles at each dispersion angle can include:

[0119] Step 901: Obtain the motion velocity of the micro-group to be tested.

[0120] For example, during the experiment, the initial velocity of the micro-particle can be obtained by applying the airflow velocity to it.

[0121] Step 902, an included angle of the asymptote line is obtained.

[0122] The included angle of the asymptote line is established based on a motion range of each broken particle.

[0123] It can be seen by observing the morphology of the particle group formed by each broken particle that the broken particle forms an asymptote line with an opening angle of about 60 degrees, that is, the motion range of the broken particle is within the included angle range of the two asymptote lines, and the dispersion angle is also within the included angle range of the asymptote line.

[0124] Step 903, based on the motion speed of the to-be-tested micro-group, the included angle of the asymptote line, and the dispersion angle of the broken particle, a maximum motion speed of the broken particle on the dispersion angle is obtained.

[0125] Reference Figure 10 As shown in the figure, Figure 10 The dashed line in the figure represents the asymptote line, and the farthest position reached by the broken particle along the direction of the dispersion angle θ after dispersion is on the asymptote line. Therefore, a vector triangle can be established based on the motion speed of the to-be-tested micro-group, the included angle of the asymptote line, and the dispersion angle of the broken particle, wherein the dispersion angle θ and half of the included angle of the asymptote line can be used as the included angle of the triangle. As the included angle of the triangle, θ and The size of the common side occupied by the two is set to the size of the motion speed of the to-be-tested micro-group; then, the size of the side of the triangle with the included angle θ with the common side is calculated, and the maximum motion speed can be obtained based on the size of the side.

[0126] In some embodiments, step 903 can include obtaining the maximum motion speed of the broken particle on the dispersion angle based on a first preset formula, the first preset formula being:

[0127]

[0128] Wherein, the V0 is the motion speed of the to-be-tested micro-group, the is half of the included angle of the asymptote line, the θ is the dispersion angle of the broken particle, and the V θ is the maximum motion speed.

[0129] Step 904, based on the maximum motion speed, the speed of the broken particle located on the dispersion angle is determined.

[0130] In some embodiments, referring to the figure, Figure 11 Step 904 can include:

[0131] Step 9041, each speed coefficient is obtained.

[0132] Wherein, the interval of the speed coefficient is between [0, 1].

[0133] Step 9042: Obtain the distribution information of each velocity coefficient in each broken particle.

[0134] Distribution information includes the second proportion of the number of broken particles with a certain velocity coefficient to the total number of particles.

[0135] This distribution information can be obtained from experimental results, and this application embodiment does not limit it. For example, refer to Figure 12 As shown, Figure 12 This is a schematic diagram showing the distribution of the velocity coefficient k among the broken particles.

[0136] Step 9043: Assign velocity coefficients to the broken particles located at the dispersion angle according to the second quantity ratio.

[0137] For example, the second quantitative ratio corresponding to the velocity coefficient k1 is p1, and the number of dispersion angle particles n1 at a certain dispersion angle can be obtained by multiplying the number of dispersion angle particles n1 by the second quantitative ratio p1, i.e., p1*n1.

[0138] By multiplying the second quantity ratio corresponding to each velocity coefficient by the number of particles at the dispersion angle, the number of broken particles corresponding to each velocity coefficient at the dispersion angle can be obtained.

[0139] Step 9044: Multiply the maximum velocity by the allocated velocity coefficient to obtain the velocity of the broken particles located at the dispersion angle. That is, the velocity of the broken particles is k*V. θ .

[0140] In this embodiment, the velocity of each broken particle at the dispersion angle is multiplied by its velocity coefficient to represent the difference between the particle's velocity and its maximum velocity. k = 1 characterizes particles reaching the asymptote; the closer k is to 1, the further forward the particle is. Figure 11 In the histogram shown, the closer k is to 1, the larger the proportion of the second quantity. This allows us to reproduce the phenomenon observed in the experiment where the concentration of particles closer to the head is higher and they are more aggregated. Thus, we can accurately present the distribution of each broken particle at the same dispersion angle.

[0141] The above steps can be used to obtain the initial position and velocity information of the broken particles formed by the breakup of a single micro-group under test. Then, the evolution process of each broken particle over time can be calculated based on the aerodynamic model to obtain the evolution process of the particle cloud formed by each broken particle.

[0142] For example, you can refer to Figure 13 As shown, obtaining the motion trajectory of each fragmented particle formed by the micro-cluster under test over time can include:

[0143] Step 1301: Obtain the drag force applied to the broken particles.

[0144] In some embodiments, the drag force F D The drag force F

[0145]

[0146] wherein C D is a drag coefficient, p fuel is the density of the broken particle, V p is the slip velocity of the broken particle in the airflow, the slip velocity being the relative velocity between the airflow and the broken particle, A p is the frontal area of the broken particle, the frontal area of the broken particle being the projected area of the broken particle in the direction of motion.

[0147] The drag coefficient can be determined according to a drag force model, which includes, but is not limited to, a Schiller-Naumann drag force model, a Moore drag force model, a Morsi-Alexander drag force model, a Clift drag force model, etc.

[0148] The Schiller-Naumann drag force model is suitable for fluid-fluid systems, and in some embodiments, the drag coefficient C D may be determined according to the Schiller-Naumann drag force model, which is shown as follows:

[0149]

[0150] wherein Re is the Reynolds number, which is a dimensionless number used to characterize the flow of a fluid.

[0151] At step 1302, a Newton motion equation of the broken particle is established based on the mass of the broken particle and the drag force.

[0152] The mass of the broken particle can be calculated according to the following formula:

[0153]

[0154] wherein m is the mass of the broken particle, is the volume of the broken particle, p fuel is the density of the broken particle, and g is the acceleration of gravity.

[0155] The established Newton motion equation of the broken particle can be shown as follows:

[0156]

[0157] wherein m is the mass of the broken particle, V pF represents the sliding velocity of the broken particles in the airflow, t represents time, and F represents the sliding velocity of the broken particles in the airflow. D The drag force applied to the broken particles is m, where m is the mass of the broken particles and g is the gravitational acceleration.

[0158] When flying in a static flow field, a near-spherical microparticle with a density of O(103) kg / m3 will be subjected to gravity and drag force, while buoyancy, pressure gradient force, virtual mass force and lift can be ignored. Therefore, by establishing Newton's equation of motion through gravity and drag force, the force analysis of the broken particles can be made more accurate, and the amount of calculation can be reduced.

[0159] In other embodiments, in order to further improve the accuracy of the force analysis of the broken particles, and thus improve the accuracy of the motion analysis of the broken particles, it is possible to establish Newton's equations of motion based on the buoyancy, pressure gradient force, virtual mass force and lift applied to the broken particles.

[0160] Step 1303: Obtain the trajectory of the broken particles over time based on Newton's equations of motion.

[0161] Using the aforementioned Newtonian equations of motion, the current acceleration of the broken particle can be obtained. The acceleration reflects the rate of change of the particle's velocity over time, thus allowing us to determine the particle's velocity over time based on the acceleration. Based on the particle's velocity over time, we can obtain the particle's trajectory. After determining the velocity and trajectory of each broken particle over time, we can obtain the morphological evolution of the cloud-like area formed by each broken particle after the aerodynamic breakup of the micro-cluster under test.

[0162] Based on the same idea as the method for determining particle distribution information in the above embodiments, this application also provides a device for determining particle distribution information, which can be used to execute the above-described method for determining particle distribution information. For ease of explanation, the schematic diagram of the embodiment of the device for determining particle distribution information only shows the parts related to the embodiments of this application. Those skilled in the art will understand that the illustrated structure does not constitute a limitation on the device, and may include more or fewer components than illustrated, or combine certain components, or have different component arrangements.

[0163] like Figure 14 As shown, the device for determining particle distribution information includes an information acquisition module 1401, a total quantity calculation module 1402, a proportion determination module 1403, and a concentration determination module 1404. In some embodiments, the above modules can be programmable software instructions stored in memory and executable by a processor. It is understood that in other embodiments, the above modules can also be program instructions or firmware embedded in the processor.

[0164] The information obtaining module 1401 is configured to obtain the micellar particle size of a to-be-tested micelle, and obtain the mass fraction of the broken particles with the micellar particle size after the to-be-tested micelle is broken by pneumatic breaking.

[0165] The total number calculating module 1402 is configured to determine the total particle number of the broken particles based on the micellar particle size and the mass fraction of the broken particles with the micellar particle size.

[0166] The proportion determining module 1403 is configured to obtain the dispersion angle of each broken particle, the dispersion angle being the included angle between the moving speed direction of the to-be-tested micelle and the moving speed direction of the broken particle, and determine the first number proportion of the broken particles located in each dispersion angle.

[0167] The concentration determining module 1404 is configured to determine the dispersion angle particle number of the broken particles located in each dispersion angle based on the total particle number and the first number proportion.

[0168] Figure 15 The schematic diagram of an embodiment of the electronic device.

[0169] The electronic device 100 includes a memory 20, a processor 30, and a computer program 40 stored in the memory 20 and executable on the processor 30. The processor 30 implements the steps in the above-described embodiment of the method for determining particle distribution information when executing the computer program 40, for example Figure 1 Steps 101-105 shown in the figure.

[0170] For example, the computer program 40 can also be divided into one or more modules / units, which are stored in the memory 20 and executed by the processor 30. The one or more modules / units can be a series of computer program instruction segments capable of completing a specific function, which are used to describe the execution process of the computer program 40 in the electronic device 100. For example, the information obtaining module 1401, the total number calculating module 1402, the proportion determining module 1403, and the concentration determining module 1404 shown in the figure can be divided. Figure 14

[0171] Those skilled in the art can understand that the schematic diagram is only an example of the electronic device 100 and does not constitute a limitation on the electronic device 100, which can include more or fewer components than the diagram, or combine certain components, or different components, for example, the electronic device 100 can also include an input / output device, a network access device, a bus, etc.

[0172] ​The processor 30 can be a central processing unit (CPU), and can also be other general-purpose processors, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic device, discrete gate or transistor logic, discrete hardware components, etc. The general-purpose processor can be a microprocessor, a microcontroller, or the processor 30 can also be any conventional processor.

[0173] The memory 20 can be used to store a computer program 40 and / or modules / units, and the processor 30 realizes various functions of the electronic device 100 by running or executing the computer program and / or modules / units stored in the memory 20, and calling data stored in the memory 20. The memory 20 can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, at least one application program required by a function (such as a sound playing function, an image playing function, etc.), etc.; and the data storage area can store data (such as audio data) created according to the use of the electronic device 100, etc. In addition, the memory 20 can include a high-speed random access memory, and can also include a nonvolatile memory, for example, a hard disk, a memory, a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, at least one disk storage device, a flash memory device, or other nonvolatile solid-state storage device.

[0174] The modules / units integrated in the electronic device 100, if implemented in the form of software function units and sold or used as independent products, can be stored in a computer-readable storage medium. Based on such understanding, all or part of the processes in the above-mentioned embodiment methods can also be implemented by a computer program instructing related hardware to complete, and the computer program can be stored in a computer-readable storage medium. When the computer program is executed by a processor, the steps of the above-mentioned various method embodiments can be implemented. The computer program includes computer program code, which can be in the form of source code, object code, executable files, or some intermediate forms. The computer-readable medium can include any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0175] In several embodiments provided in the present application, it should be understood that the disclosed electronic device and method can be implemented in other ways. For example, the above-described electronic device embodiments are only illustrative, for example, the division of the units is only a logical function division, and another division mode can be used in actual implementation.

[0176] In addition, each functional unit in each embodiment of the present application can be integrated in the same processing unit, or each unit can be physically present alone, or two or more units can be integrated in the same unit. The integrated unit can be realized in the form of hardware or hardware plus software function module.

[0177] It is obvious for those skilled in the art that the present application is not limited to the details of the above-described exemplary embodiments, and the present application can be implemented in other specific forms without departing from the spirit or essential characteristics of the present application. Therefore, the embodiments should be regarded as exemplary and non-limiting in any respect. In addition, it is obvious that the word "comprise" does not exclude other units or steps, and the singular does not exclude the plural. The plurality of units or electronic devices stated in the electronic device claims can also be implemented by the same unit or electronic device through software or hardware. The words first, second, etc. are used to indicate names, and do not indicate any specific order.

[0178] Finally, it should be noted that the above examples are only used to illustrate the technical solutions of the present application and are not limiting. Although the present application has been described in detail with reference to the above examples, those skilled in the art should understand that the technical solutions of the present application can be modified or equivalently replaced without departing from the spirit and scope of the technical solutions of the present application.

Claims

1. A method of determining information on a particle size distribution, characterized by, The method comprises: acquiring the micellar particle size of a to-be-tested micelle, the to-be-tested micelle forming broken particles with respective particle sizes after pneumatic breaking; acquiring the mass fraction of the broken particles with respective particle sizes, comprising: acquiring the ratio of the mass median diameter to the Sauter mean diameter of each broken particle; determining the size distribution index of the broken particle based on the ratio of the mass median diameter to the Sauter mean diameter; determining the cumulative mass fraction of the broken particles with respective particle sizes based on the size distribution index; and obtaining the mass fraction of the broken particles with respective particle sizes based on the cumulative mass fraction; Based on the micro-cluster particle size and the mass fraction of the broken particles of each micro-particle size, the total number of the broken particles is determined, including: determining the number of the broken particles of each micro-particle size based on a second preset formula, wherein the second preset formula is: wherein, is the identification of the micro-particle size, is the micro-particle size identified as , is the micro-cluster particle size of the micro-cluster to be measured, is the mass fraction of the broken particles with the micro-particle size of , is the number of the broken particles with the micro-particle size of ; and the total number of the broken particles is obtained based on the number of the broken particles of each micro-particle size. acquiring the dispersion angle of each broken particle, the dispersion angle being the included angle between the moving speed direction of the to-be-tested micelle and the moving speed direction of the broken particle; determining the first number proportion of the broken particles located at respective dispersion angles; determining the number of dispersion angle particles located at the respective dispersion angles based on the total particle number and the first number proportion.

2. The method of determining particle size distribution information according to claim 1, wherein After the step of determining the number of dispersion angle particles located at the respective dispersion angles based on the total particle number and the first number proportion, the method further comprises: acquiring the moving speed of the to-be-tested micelle; acquiring the included angle of the asymptote, the included angle of the asymptote being established based on the moving range of the broken particles; obtaining the maximum moving speed of the broken particle at the dispersion angle based on the moving speed of the to-be-tested micelle, the included angle of the asymptote, and the dispersion angle of the broken particle; determining the speed of the broken particle located at the dispersion angle based on the maximum moving speed.

3. The method of claim 2, wherein the step of determining the particle size distribution information is performed by a method comprising: The step of obtaining the maximum moving speed of the broken particle at the dispersion angle based on the moving speed of the to-be-tested micelle, the included angle of the asymptote, and the dispersion angle of the broken particle comprises: obtaining the maximum moving speed of the broken particle at the dispersion angle based on a first preset formula, the first preset formula being: ， Wherein, the is the motion speed of the micro-cluster to be measured, the is half of the included angle of the asymptote, the is the dispersion angle of the broken particles, the is the maximum motion speed.

4. The method of claim 2, wherein the step of determining the particle size distribution information is performed by a method comprising: The step of determining the speed of the broken particle located at the dispersion angle based on the maximum moving speed comprises: acquiring respective speed coefficients; acquiring distribution information of the respective speed coefficients in the broken particles, the distribution information representing the second number proportion of the number of broken particles with the speed coefficient in the total particle number; allocating speed coefficients to the broken particles located at the dispersion angle according to the second number proportion; multiplying the maximum moving speed by the allocated speed coefficient to obtain the speed of the broken particle located at the dispersion angle.

5. The method of determining the information of the particle distribution according to any one of claims 1 to 4, wherein, After the step of determining the number of dispersion angle particles located at the respective dispersion angles based on the total particle number and the first number proportion, the method further comprises: acquiring the drag force applied to the broken particle, establishing the Newton motion equation of the broken particle based on the mass of the broken particle and the drag force; obtaining the motion trajectory of the broken particle changing with time based on the Newton motion equation.

6. A device for determining information on a particle size distribution, characterized by The device comprises: The information acquisition module is configured to acquire the micellar particle size of a to-be-tested micelle, and the to-be-tested micelle is broken by pneumatic breaking to form broken particles with each particle size; and to acquire a mass fraction of the broken particles with each particle size, including: acquiring a ratio of a mass median diameter to a Sauter mean diameter of each broken particle; determining a size distribution index of the broken particle based on the ratio of the mass median diameter to the Sauter mean diameter; determining a cumulative mass fraction of the broken particles with each particle size based on the size distribution index; and obtaining the mass fraction of the broken particles with each particle size based on the cumulative mass fraction; The total number calculation module determines the total number of broken particles based on the micro-cluster particle size and the mass fraction of the broken particles of each particle size, and includes: determining the number of broken particles of each particle size based on a second preset formula, wherein the second preset formula is: wherein, is the identification of the particle size, is the particle size identified as , is the micro-cluster particle size of the micro-cluster to be measured, is the mass fraction of the broken particles with the particle size of , is the number of broken particles with the particle size of ; and the total number of broken particles is obtained based on the number of broken particles of each particle size. The proportion determination module is configured to acquire a dispersion angle of each broken particle, the dispersion angle being an included angle between a moving speed direction of the to-be-tested micelle and a moving speed direction of the broken particle; and to determine a first number proportion of the broken particles located at each dispersion angle; The concentration determination module is configured to determine a dispersion angle particle number located at each dispersion angle based on the total particle number and the first number proportion. 7.An electronic device comprising a processor and a memory, wherein The memory is configured to store instructions, and the processor is configured to invoke the instructions in the memory, so that the electronic device executes the method for determining particle distribution information according to any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions, and when the computer instructions run on the electronic device, the electronic device executes the method for determining particle distribution information according to any one of claims 1 to 5.

Images