Particle motion trail simulation method and electronic equipment

Abstract

The invention relates to a particle motion trail simulation method and electronic equipment, and relates to the technical field of pollutant trail simulation, and the method comprises the steps: determining a rebound position after a particle collides with a shoreline under the condition that a first motion trail of the particle intersects with the shoreline; according to the rebound position, obtaining a target position where the particle arrives after moving along the shoreline; according to the initial position and the target position of the particle, obtaining a corrected second motion track of the particle; the starting position is a starting point of the first motion track. By using the particle motion trail simulation method provided by the invention, a more accurate motion trail of the particle can be obtained.

Description

Technical Field

[0001] The present disclosure relates to the technical field of pollutant trajectory simulation, and in particular, to a method for simulating the movement trajectory of particles and an electronic device. Background Art

[0002] Currently, due to reasons such as industrial emissions and transportation accidents, sudden water environment pollution accidents occur. It is possible to simulate the diffusion and migration of pollutants in water bodies, thereby tracking the movement trajectories of pollutants in water bodies, predicting the distribution and concentration changes of pollutants in water bodies, and providing a basis for the treatment of pollutants.

[0003] In the related art, during the process of simulating the movement trajectory of particles in pollutants, there are errors in the simulated movement trajectory of the particles. Summary of the Invention

[0004] The purpose of the present disclosure is to provide a method for simulating the movement trajectory of particles and an electronic device to solve the above technical problems.

[0005] To achieve the above purpose, the present disclosure provides a method for simulating the movement trajectory of particles, including: When the first movement trajectory of the particle intersects with the shoreline, determining the rebound position of the particle after colliding with the shoreline; Based on the rebound position, obtaining the target position reached by the particle after moving along the shoreline; Based on the starting position and the target position of the particle, obtaining the corrected second movement trajectory of the particle; the starting position is the starting point of the first movement trajectory.

[0006] Optionally, the obtaining the target position reached by the particle after moving along the shoreline based on the rebound position includes: Obtaining a first coordinate difference based on a second position and a first position on the shoreline; Obtaining a second coordinate difference based on the rebound position and the starting position; Obtaining the target position based on a first coefficient, a second coefficient, the first coordinate difference, and the second coordinate difference; the first coefficient is the coefficient of the first coordinate difference, and the second coefficient is the coefficient of the second coordinate difference.

[0007] Optionally, the first coordinate difference includes a first vertical coordinate difference, and the second coordinate difference includes a second vertical coordinate difference; the obtaining the target position based on a first coefficient, a second coefficient, the first coordinate difference, and the second coordinate difference includes: Obtaining a first target value based on the first coefficient and the first vertical coordinate difference; Obtaining a second target value based on the second coefficient and the second vertical coordinate difference; Based on the first target value and the second target value, obtain the abscissa in the target position.

[0008] Optionally, the first coordinate difference includes a first abscissa difference, and the second coordinate difference includes a second abscissa difference; the obtaining of the target position according to the first coefficient, the second coefficient, the first coordinate difference, and the second coordinate difference includes: Based on the first coefficient and the first abscissa difference, obtain a third target value; Based on the second coefficient and the second abscissa difference, obtain a fourth target value; Based on the third target value and the fourth target value, obtain the ordinate in the target position.

[0009] Optionally, the method further includes: Based on a first distance, a second distance, the starting position, and the bouncing position, obtain a first value; the first distance is the distance from the starting position to the bouncing position, the second distance is the distance from the bouncing position to the end position, and the end position is the end position on the first motion trajectory; Based on the first coordinate difference and the second coordinate difference, obtain a second value; Based on the first value and the second value, obtain the first coefficient.

[0010] Optionally, the method further includes: Based on the second distance, a third distance, the bouncing position, and the first coordinate difference, obtain a third value; the second distance is the distance from the bouncing position to the end position, the end position is the end position on the first motion trajectory, and the third distance is the distance between the second position and the first position; Based on the first coordinate difference and the second coordinate difference, obtain a second value; Based on the third value and the second value, obtain the second coefficient.

[0011] Optionally, the first position and the second position are the positions of two adjacent shore points on the shoreline.

[0012] Optionally, the method further includes: In the case where the target position is outside the shoreline range, use the bouncing position as the initial position of the particle, and re-determine the target position reached by the particle after moving along the shoreline from the updated initial position until the updated target position is within the shoreline range.

[0013] Optionally, when the first movement trajectory of the particle intersects with the shoreline, determining the rebound position of the particle after colliding with the shoreline includes: Dividing the space where the shoreline is located into a plurality of grids; Determining target grids in the plurality of grids whose distances from the initial position of the particle are less than a preset distance; When there is the shoreline in the target grid and the first movement trajectory of the particle intersects with the shoreline, determining the rebound position of the particle after colliding with the shoreline.

[0014] To achieve the above object, the present disclosure provides an electronic device, including: A memory storing a computer program thereon; A processor configured to execute the computer program in the memory to implement the steps of the method for simulating the movement trajectory of a particle proposed by the present disclosure.

[0015] Through the above technical solution, the rebound position of the particle after the first movement trajectory of the particle collides with the shoreline can be determined, and the target position reached by the particle after moving along the shoreline can be obtained based on the rebound position; finally, the corrected second movement trajectory of the particle can be obtained based on the starting position and the target position of the particle. In this process, the target position where the particle moves along the direction of the longshore current is obtained based on the rebound position of the particle, taking into account the influence of the longshore current of the shoreline on the movement trajectory of the particle, making the obtained second movement trajectory more accurate.

[0016] Other features and advantages of the present disclosure will be described in detail in the subsequent specific implementation section. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The drawings are used to provide a further understanding of the present disclosure and constitute a part of the specification, and are used to explain the present disclosure together with the following specific implementation manners, but do not constitute a limitation to the present disclosure. In the drawings: Figure 1 is a flowchart of the steps of a method for simulating the movement trajectory of a particle shown according to an exemplary embodiment.

[0018] Figure 2 is a schematic diagram of the first movement trajectory and the second movement trajectory of a particle shown according to an exemplary embodiment.

[0019] Figure 3 is a flowchart of the steps of a method for simulating the movement trajectory of a particle shown according to an exemplary embodiment.

[0020] Figure 4 is a schematic diagram of the first movement trajectory and the second movement trajectory of a particle shown according to an exemplary embodiment.

[0021] Figure 5 It is a flowchart of the steps of a method for simulating the movement trajectory of a particle shown according to an exemplary embodiment.

[0022] Figure 6 It is a flowchart of the steps of a method for simulating the movement trajectory of a particle shown according to an exemplary embodiment.

[0023] Figure 7 It is a flowchart of the steps of a method for simulating the movement trajectory of a particle shown according to an exemplary embodiment.

[0024] Figure 8 It is a schematic diagram of the first movement trajectory and the second movement trajectory of a particle shown according to an exemplary embodiment.

[0025] Figure 9 It is a block diagram of a device for simulating the movement trajectory of a particle shown according to an exemplary embodiment.

[0026] Figure 10 It is a block diagram of an electronic device shown according to an exemplary embodiment.

[0027] Figure 11 It is a block diagram of an electronic device shown according to an exemplary embodiment. Detailed implementation manners

[0028] The following will describe in detail the specific implementation manners of the present disclosure with reference to the accompanying drawings. It should be understood that the specific implementation manners described herein are only for explaining and interpreting the present disclosure, and are not used to limit the present disclosure.

[0029] Floating objects such as pollutants floating in water bodies can be regarded as a collection of a series of particles, and the movement trajectories of the particles in the pollutants can be predicted through the Lagrangian particle tracking algorithm.

[0030] As an important tool for simulating the diffusion and migration of pollutants in water bodies, the Lagrangian particle tracking algorithm can accurately track the movement trajectories of each particle, thereby comprehensively depicting the dynamic behavior of pollutants in water bodies. The Lagrangian particle tracking algorithm not only depends on the fine description of the water flow field, but also highly depends on the accurate processing of complex shoreline boundary conditions. Among them, the shoreline, as the boundary between the water body and the land, its shape, position and characteristics have a direct impact on the water flow pattern, and thus deeply affect the diffusion path, deposition rate and end position of the particles.

[0031] When simulating sudden water environment pollution accidents, the treatment of the shoreline is particularly crucial. The shoreline has diverse forms, including straight river shorelines and meandering coastlines. These different forms will lead to significant differences in water flow velocity, direction, and vortex structure. This requires that the Lagrangian particle tracking algorithm not only be able to accurately simulate the movement of particles on the water surface but also consider complex physical processes such as rebound, adsorption, and deposition when particles interact with the shoreline.

[0032] In related technologies, for non-transmissive shorelines, there are the following three methods for dealing with particles: The first method of treatment, reflection: After the particle collides with the shoreline, it rebounds back into the water body.

[0033] The second method of treatment, absorption: After the particle collides with the shoreline, it deposits or adheres to the shoreline.

[0034] The third method of treatment, partial reflection: Some particles rebound back into the water body after colliding with the shoreline, and some particles deposit or adhere to the shoreline after colliding with the shoreline.

[0035] For the method of dealing with the particle rebounding back into the water body after colliding with the shoreline, the reflection angle and displacement of the particle after colliding with the shoreline are calculated based on the incident angle of the particle towards the shoreline, and then the target position reached by the particle after colliding with the shoreline is simulated. However, in actual situations, affected by complex water flows such as the longshore current of the shoreline, after the particle collides with the shoreline, it will not simply return to the water body according to the reflection angle but will have a certain displacement along the shoreline direction. Then, there is an error in the accuracy of the target position of the particle simulated according to the reflection angle and displacement, and there is also a deviation in the movement trajectory of the particle obtained based on this target position.

[0036] Based on this, the present disclosure proposes a method for simulating the movement trajectory of particles, Figure 1 which is a flowchart of the steps of the method for simulating the movement trajectory of particles proposed by the present disclosure. The method for simulating the movement trajectory of particles includes the following steps: In step S10, when the first movement trajectory of the particle intersects with the shoreline, determine the rebound position of the particle after colliding with the shoreline.

[0037] Among them, the first movement trajectory of the particle is the initial movement trajectory of the particle, and can also be understood as the movement trajectory before the particle is corrected. There are a starting position and an ending position on this first movement trajectory. The starting position is the starting point of the movement on the first movement trajectory, and the ending position is the ending point of the movement on the first movement trajectory. This ending position can also be understood as the position that the particle naturally reaches under the action of the water flow and without colliding with the shoreline. If the first movement trajectory intersects with the shoreline, then the starting position and the ending position will be on both sides of the shoreline respectively.

[0038] See also Figure 2 As shown, taking the coastline being the CD segment and the first motion trajectory being the AB segment as an example, since the first motion trajectory AB intersects with the coastline CD, the starting position A and the end position B on the first motion trajectory are located on both sides of the coastline CD.

[0039] Among them, there are at least two particle positions on the first motion trajectory. If the first motion trajectory includes two particle positions, then the first motion trajectory is composed of a starting position and an end position, and no other particle positions are included between the starting position and the end position on the first motion trajectory; if the first motion trajectory includes three or more particle positions, then the starting position and the end position on the first motion trajectory are separated by at least one particle position. It can be understood that if the first motion trajectory of the particle in each time period is to be accurately obtained, then each first motion trajectory can be composed of two particle positions, the starting position and the end position. The starting position is the position of the particle at the previous time point, and the end position is the position of the particle at the next time point. The particle positions at two adjacent time points constitute a first motion trajectory, and the first motion trajectory of the particle at each adjacent time point can be corrected later to obtain the corrected second motion trajectory of the particle at each adjacent time point.

[0040] Among them, if the first motion trajectory is composed of two particle positions, a starting position and an end position, the end position of the particle can be obtained according to the starting position, the moving speed of the particle in the water body, and the interval between two adjacent time points. The end position is the position that the particle reaches at the next time point when the particle does not turn or rebound during the flow of the water body. For example, the moving distance of the particle is obtained based on the moving speed of the particle in the water body and the interval between two adjacent time points, and then the horizontal axis component of the moving distance on the horizontal axis is superimposed on the horizontal coordinate of the starting position of the particle to obtain the horizontal coordinate of the end position, and the vertical axis component of the moving distance on the vertical axis is superimposed on the vertical coordinate of the starting position of the particle to obtain the vertical coordinate of the end position, and the horizontal coordinate of the end position and the vertical coordinate of the end position constitute the end position.

[0041] Optionally, the intersection point position of the first motion trajectory and the shoreline may be determined, and then the rebound position may be obtained based on the intersection point position and the rebound coefficient. For example, the rebound position may be obtained by multiplying the intersection point position by the rebound coefficient.

[0042] Among them, the rebound position after the particle collides with the shoreline, the starting position of the particle and the end position of the particle are all located on the same straight line. The rebound position is the rebound position after the particle collides with the shoreline without considering the influence of complex water currents such as coastal currents.

[0043] Among them, the influencing factors of the bounce coefficient include the shoreline type and the material properties of the shoreline, etc. The shoreline types include waterweeds, stones, mud, etc. The bounce coefficients of waterweeds, mud, etc. are relatively small. After the particles collide with waterweeds and mud, they are easily adsorbed by waterweeds and mud and are not easy to bounce. Therefore, the distance between the bounce position of the particles after colliding with waterweeds, mud and other shorelines and the shoreline is relatively close; while the bounce coefficient of stones is relatively large, and the particles are easy to bounce after colliding with stones. Therefore, the distance between the bounce position of the particles after colliding with stones and the shoreline is relatively far. The material properties of the shoreline include properties such as oiliness, viscosity, rigidity, elasticity, etc. The bounce coefficients of shorelines with material properties such as oiliness and viscosity are relatively small. When the particles collide with shorelines with material properties such as oiliness and viscosity, the distance between the bounce position after the collision and the shoreline is relatively close; while the bounce coefficients of material properties such as rigidity and elasticity are relatively large. After the particles collide with shorelines with material properties such as rigidity and elasticity, the distance between the bounce position after the collision and the shoreline is far.

[0044] Optionally, when it is determined that the shoreline is closed and it is determined that the first motion trajectory of the particle intersects the closed shoreline, the bounce position of the particle after colliding with the shoreline is determined.

[0045] Among them, since the shoreline is composed of a series of discrete data points, if a series of received data points cannot form a continuous closed shoreline, these data points are considered abnormal data, and new data points will be obtained again, and a closed shoreline will be obtained again based on the new data points.

[0046] Among them, when judging whether the shoreline is closed, if the number of multiple obtained data points is greater than the preset number and there are no data points with the same coordinates among the multiple data points, then the line segment formed by the multiple obtained data points is considered to be the shoreline. The fact that there are no data points with the same coordinates among the multiple data points includes: the coordinates of the first data point and the tail data point among the multiple data points are different, and the coordinates of the other data points except the first data point and the tail data point among the multiple data points are different from each other.

[0047] Since the shoreline is a curve or a straight line, the number of data points on the shoreline will be greater than the preset number (for example, 3), and the first data point and the tail data point on the shoreline are not the same coordinate, and there are no data points with the same coordinates among the other data points except the first data point and the tail data point on the shoreline. When these conditions are met, it is considered that the multiple obtained data points can represent a shoreline.

[0048] In step S20, according to the bounce position, the target position reached by the particle after moving along the shoreline is obtained.

[0049] Among them, affected by complex water flows such as the longshore current along the shoreline, after the particle rebounds, it will move along the shoreline in the downstream direction of the shoreline. Then, the line connecting the rebound position and the target position is parallel to the shoreline. The target position is the position estimated to be reached after the particle collides with the shoreline and moves along the longshore current.

[0050] Optionally, the target position can be obtained based on the starting position of the particle on the first movement trajectory, the end position of the first movement trajectory, the first position on the shoreline, the second position on the shoreline, and the rebound position of the particle.

[0051] Among them, the first position and the second position on the shoreline can be two adjacent shore point positions on the shoreline. There are multiple shore point positions on the shoreline. Two shore point positions closest to the initial position of the particle can be selected from the multiple shore point positions, and these two shore point positions are used as the first position and the second position respectively. For example, please refer to Figure 2 As shown, point C, which is the closest to the initial position A of the particle among the multiple shore point positions, can be selected as the first position, and point D as the second position.

[0052] It can be understood that although the shoreline of a water body is a closed geometric shape, the shoreline in the obtained shoreline data of the water body is not a complete closed shoreline, but multiple shore line segments, and each shore line segment is composed of multiple shore point positions. When the particle collides with the shoreline, it collides with the shore line segment that is closer to the particle among the multiple shore line segments. Therefore, the line connecting the two shore point positions that are closest to the particle among the multiple shore point positions on the shoreline is used as the shore line segment intersecting with the particle. Of course, the line connecting the two shore point positions that are closest to the rebound position among the multiple shore point positions can also be used as the shore line segment intersecting with the particle.

[0053] Optionally, the target position can be obtained based on the rebound position of the particle, the second distance between the rebound position of the particle and the end position, and the first angle between the shoreline and the horizontal axis.

[0054] In step S30, based on the starting position and the target position of the particle, the corrected second movement trajectory of the particle is obtained.

[0055] Among them, the second movement trajectory is the movement trajectory of the particle corrected considering the influence of complex water flows such as the longshore current. The second movement trajectory includes two movement trajectories, one of which is from the initial position of the particle towards the shoreline, and the other is the same as the trend of the shoreline.

[0056] Optionally, the line connecting the starting position and the target position of the particle can be used as the corrected second movement trajectory of the particle. This second movement trajectory is the movement trajectory of the particle within the time period composed of two adjacent time points.

[0057] Optionally, the starting position, the bouncing position, and the target position of the particle can also be connected into a line as the second corrected movement trajectory of the particle.

[0058] Among them, the accurate second movement trajectory should be the movement trajectory from the starting position, passing through the intersection point between the first movement trajectory and the shoreline, the bouncing position, and then reaching the target position in sequence. Although connecting the starting position of the particle to the target position as the second movement trajectory has low accuracy from a microscopic perspective, when viewing the flow trajectory diagram of pollutants in the water body, it is not to infinitely magnify the flow trajectory diagram to microscopically observe the movement trajectory of each particle in the pollutants, but to view the flow trajectory of pollutants in the water body from a macroscopic perspective. Therefore, even if the starting position of the particle is connected to the target position as the second movement trajectory, the flow trajectory of pollutants in the water body viewed from a macroscopic perspective is accurate. Of course, if a more accurate flow trajectory diagram of the particle is to be shown, the starting position, the intersection point, the bouncing position, and the target position can also be connected into a line to construct an accurate flow trajectory diagram.

[0059] Through the above technical solution, the bouncing position after the first movement trajectory of the particle collides with the shoreline can be determined, and the target position reached by the particle after moving along the shoreline can be obtained based on the bouncing position; finally, the second corrected movement trajectory of the particle can be obtained based on the starting position and the target position of the particle. In this process, the target position where the particle moves along the direction of the coastal current is obtained based on the bouncing position of the particle, considering the influence of the coastal current of the shoreline on the movement trajectory of the particle, making the obtained second movement trajectory more accurate.

[0060] Figure 3 This is an exemplary embodiment involved in the above step S20, which is used to interpret an exemplary embodiment of obtaining the target position where the particle moves along the shoreline based on the bouncing position, and includes the following steps: In step S21, the target position can be obtained according to the bouncing position of the particle, the second distance between the bouncing position of the particle and the end position, and the first included angle between the shoreline and the horizontal axis.

[0061] Among them, the second distance between the bouncing position of the particle and the end position is equal to the distance between the bouncing position of the particle and the target position. Therefore, the second distance can be used as the distance between the bouncing position of the particle and the target position. For example, refer to Figure 4 As shown, taking the bouncing position as point E, the end position as point B, and the target position as point F as an example, the second distance between E and B is equal to the distance between E and F.

[0062] Among them, the line connecting the rebound position of the particle to the target position is parallel to the shoreline. Therefore, the first angle between the shoreline and the horizontal axis can be used as the angle between the line connecting the rebound position of the particle to the target position and the horizontal axis. For example, referring to Figure 4 As shown, taking the shoreline as CD and the line from the rebound position of the particle to the target position as EF, CD is parallel to EF.

[0063] Therefore, obtaining the target position based on the rebound position of the particle, the second distance between the rebound position of the particle and the end position, and the first angle between the shoreline and the horizontal axis can be regarded as obtaining the target position based on the rebound position of the particle, the second distance between the rebound position of the particle and the target position, and the angle between the line connecting the rebound position of the particle to the target position and the horizontal axis.

[0064] For example, the horizontal component of the second distance projected on the horizontal axis can be obtained according to the second distance between the rebound position of the particle and the end position and the cosine value of the first angle between the shoreline and the horizontal axis, and then the horizontal component is added to or subtracted from the abscissa of the rebound position to obtain the abscissa of the target position; the vertical component of the second distance projected on the vertical axis can be obtained according to the second distance between the rebound position of the particle and the end position and the sine value of the first angle between the shoreline and the horizontal axis, and then the vertical component is added to or subtracted from the ordinate of the rebound position to obtain the ordinate of the target position; finally, the abscissa and ordinate of the target position form the target position.

[0065] Taking the rebound position of the particle as Figure 4 the E point in Figure 4 and the target position as Figure 4 the F point in Figure 4 and the second distance as L2 in (1) In formula (1), is the abscissa of the target position, is the abscissa of the rebound position, is the ordinate of the target position, is the ordinate of the rebound position, L2 is the second distance between the rebound position and the target position F, and is also the distance between the rebound position and the end position B, is the first angle between the EF line segment and the horizontal axis of the coordinate system, and is also the angle between the shoreline and the horizontal axis of the coordinate system.

[0066] Through the above technical solution, the target position of the particle after moving along the shoreline can be obtained based on the rebound position of the particle, the second distance between the rebound position of the particle and the end position, and the first angle between the shoreline and the horizontal axis. Thus, the obtained target position of the particle is the target position reached after being affected by the longshore current along the shoreline of the particle. Then, the second movement trajectory obtained based on the initial position and the target position of the particle will be more accurate.

[0067] Figure 5 This is an exemplary embodiment related to the above step S20, which is used to interpret an exemplary solution for obtaining the target position according to the starting position of the particle on the first movement trajectory, the end position of the first movement trajectory, the first position on the shoreline, the second position on the shoreline, and the rebound position of the particle, and includes the following steps: In step S22, a first coordinate difference is obtained according to the second position and the first position on the shoreline.

[0068] The first coordinate difference includes a first abscissa difference and a first ordinate difference. The first abscissa difference can be obtained according to the abscissa of the second position and the abscissa of the first position on the shoreline; the first ordinate difference can be obtained according to the ordinate of the second position and the ordinate of the first position on the shoreline.

[0069] For example, please refer to Figure 2 as shown. Taking the second position as point D ( , ) and the first position as point C ( , ) as an example, the first coordinate difference includes the first abscissa difference and the first ordinate difference .

[0070] In step S23, a second coordinate difference is obtained according to the rebound position and the starting position.

[0071] The second coordinate difference can be obtained by subtracting the starting position from the rebound position. The second coordinate difference includes a second abscissa difference and a second ordinate difference. The second abscissa difference can be obtained according to the abscissa of the rebound position and the abscissa of the starting position; the second ordinate difference can be obtained according to the ordinate of the rebound position and the ordinate of the starting position.

[0072] For example, please refer to Figure 2 as shown. Taking the starting position as point A ( , ) and the rebound position as point E ( , ) as an example, the second coordinate difference includes the second abscissa difference and the second ordinate difference .

[0073] In step S24, the target position is obtained according to the first coefficient, the second coefficient, the first coordinate difference, and the second coordinate difference.

[0074] Optionally, the first coordinate difference can be adjusted by using the first coefficient, and the second coordinate difference can be adjusted by using the second coefficient to obtain the target position.

[0075] Exemplarily, the first abscissa difference and the first ordinate difference are adjusted by using the first coefficient, and the second abscissa difference and the second ordinate difference are adjusted by using the second coefficient to obtain the target position. Specifically, it includes: obtaining a first target value according to the first coefficient and the first ordinate difference; obtaining a second target value according to the second coefficient and the second ordinate difference; obtaining the abscissa in the target position according to the first target value and the second target value; obtaining a third target value according to the first coefficient and the first abscissa difference; obtaining a fourth target value according to the second coefficient and the second abscissa difference; obtaining the ordinate in the target position according to the third target value and the fourth target value.

[0076] For example, the calculation formulas for obtaining the abscissa and ordinate of the target position are as follows: (2) In formula (2), is the first coefficient, is the first ordinate difference, is the first abscissa difference; is the second coefficient, is the second ordinate difference, is the second abscissa difference.

[0077] It can be seen from the above formula (2) that the first coefficient can be used to adjust the first ordinate difference to obtain the first target value , and the second coefficient can be used to adjust the second ordinate difference to obtain the second target value ; then the abscissa of the target position is obtained according to the first target value and the second target value; the first coefficient can also be used to adjust the first abscissa difference to obtain the third target value , and the second coefficient can be used to adjust the second abscissa difference to obtain the fourth target value , and then the ordinate of the target position is obtained according to the third target value and the fourth target value.

[0078] Optionally, for the first coefficient, the following sub-steps A1 to A3 can be used to obtain the first coefficient: Sub-step A1: Obtain a first value based on the first distance, the second distance, the starting position, and the rebounding position.

[0079] Among them, the first distance is the distance from the starting position to the rebounding position, the second distance is the distance from the rebounding position to the end position, and the end position is the end position on the first movement trajectory.

[0080] Optionally, the first value can be obtained based on the first distance, the second distance, the cosine value of the second angle, the second abscissa difference between the abscissa of the rebounding position and the abscissa of the starting position, the second ordinate difference between the ordinate of the rebounding position and the ordinate of the starting position, the abscissa of the rebounding position, and the ordinate of the rebounding position. Among them, the second angle is the angle between the line connecting the starting position and the rebounding position and the line connecting the rebounding position and the target position. For example Figure 2 the angle between L1 and L2. Since the angle between L1 and L2 is an obtuse angle and is a negative value when substituted into the cosine formula, the supplementary angle of this angle is used as the angle between L1 and L2; the first value is the numerator in the first coefficient.

[0081] For example, the calculation formula for obtaining the first value is as follows: (3) In formula (3), is the first value, is the first distance, is the second distance, is the second angle, is the second abscissa difference, is the second ordinate difference, is the abscissa of the rebounding position, is the ordinate of the rebounding position.

[0082] It can be seen from the above formula (3) that the first value can be obtained by superimposing the product of the first distance, the second distance, and the cosine value of the second angle, the product of the abscissa of the rebounding position and the second abscissa difference, and the product of the ordinate of the rebounding position and the second ordinate.

[0083] Sub-step A2: Obtain a second value based on the first coordinate difference and the second coordinate difference.

[0084] Optionally, the second value can be obtained based on the first abscissa difference, the first ordinate difference, the second abscissa difference, and the second ordinate difference. Among them, the second value is the denominator of the first coefficient.

[0085] For example, the calculation formula for obtaining the second value is as follows: (4) In formula (4), is the second value, is the first abscissa difference, is the first ordinate difference, is the second abscissa difference, is the second ordinate difference.

[0086] It can be seen from the above formula (4) that the second value can be obtained by subtracting the product of the first abscissa difference and the second ordinate difference from the product of the second abscissa difference and the first ordinate difference.

[0087] Sub-step A3, obtaining the first coefficient according to the first value and the second value.

[0088] The ratio between the first value and the second value can be used as the first coefficient. For example, the ratio between the first value and the second value is used as the first coefficient .

[0089] Optionally, for the second coefficient, the following sub-steps B1~B3 can be adopted to obtain the second coefficient: Sub-step B1, obtaining a third value according to the second distance, the third distance, the bounce position and the first coordinate difference.

[0090] Among them, the second distance is the distance from the bounce position to the end position, the end position is the end position on the first motion trajectory, and the third distance is the distance between the second position and the first position.

[0091] Optionally, the third value can be obtained according to the second distance, the third distance, the first abscissa difference, the first ordinate difference, the abscissa of the bounce position and the ordinate of the bounce position. Among them, the third value is the numerator in the second coefficient.

[0092] For example, the calculation formula for obtaining the third value is as follows: (5) In formula (5), is the third value, is the second distance, is the third distance, is the first abscissa difference, is the first ordinate difference, is the abscissa of the bounce position, is the ordinate of the bounce position.

[0093] As can be seen from the above formula (5), the third value can be obtained by subtracting the product of the abscissa difference between the bounce position and the first abscissa and the product of the ordinate difference between the bounce position and the first ordinate from the product of the second distance and the third distance.

[0094] Sub-step B2: Obtain a second value according to the first coordinate difference and the second coordinate difference.

[0095] The method for obtaining the second value refers to the above sub-step A2 and will not be elaborated here. The second value can also be used as the denominator of the second coefficient.

[0096] Sub-step B3: Obtain the second coefficient according to the third value and the second value.

[0097] The ratio between the third value and the second value can be used as the second coefficient. For example, taking the ratio between the third value and the second value as the second coefficient .

[0098] In some scenarios, as shown in Figure 2 , taking the initial position as point A, the end position as point B, the second position on the shoreline as point D, the first position on the shoreline as point C, and the bounce position as point E as an example, the following vector formula can be obtained: (6) By arranging formula (6), the following formula can be obtained:

[0099] (7) By further transforming formula (7), the above formula (2) can be obtained.

[0100] Through the above technical solution, for any particle in the pollutant, the starting position of the particle at the previous time point, the end position at the next time point, the first position on the shoreline, the second position on the shoreline, and the bounce position of the particle after colliding with the shoreline can be obtained. By inputting these parameters into the above formula (2), the target position of the particle after moving along the shoreline can be obtained. Thus, the obtained target position of the particle is the target position reached after being affected by the longshore current along the shoreline of the particle. Then, the second movement trajectory obtained based on the initial position and the target position of the particle will be more accurate.

[0101] Figure 6 is an exemplary embodiment involved in the above step S10, which is used to interpret an exemplary solution for obtaining the bounce position of the particle after colliding with the shoreline, and includes the following steps: In step S11, the space where the shoreline is located is divided into a plurality of grids.

[0102] For example, the space where the shoreline is located can be divided into a plurality of grids such as grid A to grid F, and the position of each grid represents the position of a certain sub-region in the space where the shoreline is located.

[0103] In step S12, target grids in the plurality of grids whose distance from the initial position of the particle is less than a preset distance are determined.

[0104] Target grids whose distance from the initial position of the particle is less than a preset distance can be screened out from the plurality of grids, and then it is determined whether there is a shoreline in this part of the target grids.

[0105] In step S13, when there is a shoreline in the target grid and the first movement trajectory of the particle intersects the shoreline, the rebound position of the particle after colliding with the shoreline is determined.

[0106] When there is a shoreline in the target grid, it is determined whether the shoreline in the target grid intersects the first movement trajectory of the particle. If they intersect, the rebound position of the particle after colliding with the shoreline is determined.

[0107] It can be understood that the shoreline proposed in any embodiment of the present disclosure can refer to a shoreline segment in the shoreline of a water body.

[0108] Since there are a large number of shoreline segments in the space where the shoreline is located, if it is calculated whether all the shoreline segments intersect the first movement trajectory of the particle, the calculation amount is large.

[0109] Through the above technical solution, the space where the shoreline is located can be divided into a plurality of grids, then target grids whose distance from the initial position of the particle is less than a preset distance are screened out from the plurality of grids, and then it is determined whether the shoreline segment in the target grid intersects the first movement trajectory of the particle, that is, it is determined whether the shoreline segment in the target grid intersects the line connecting the particle positions at two adjacent time points, so as to eliminate the shoreline segments in the grids whose distance from the initial position of the particle is above the preset distance, without involving all the shoreline segments in the intersection calculation, reducing the data calculation amount.

[0110] Figure 7 This is an exemplary embodiment involved in the present disclosure, which is used to interpret the processing solution when the obtained target position is outside the shoreline range, and includes the following steps: In step S40, when the target position is outside the shoreline range, the bounce position is taken as the initial position of the particle, and the target position that the particle reaches after moving along the shoreline starting from the updated initial position is re-determined until the updated target position is within the shoreline range.

[0111] The particles in the pollutant drift in the water body, which can also be understood as the particles in the pollutant drift within the shoreline range. If the calculated target position is outside the shoreline range, it means that there is still an error in the simulated target position F. At this time, the bounce position can be taken as the updated initial position of the particle, the target position can be taken as the updated end position of the particle, and the updated first motion trajectory of the particle can be obtained based on the updated initial position and end position of the particle; then, when the updated first motion trajectory of the particle intersects with the updated shoreline, a new bounce position after the particle collides with the updated shoreline is determined; then, based on the new bounce position, the target position that the particle reaches after moving along the new shoreline is obtained, and this is repeatedly executed until the obtained target position is within the shoreline range, and then the starting position of the particle is connected to the target position within the shoreline range to obtain the corrected second motion trajectory of the particle.

[0112] Please refer to Figure 8 As shown, if the obtained target position F is outside the shoreline range, the bounce position E point of the particle can be taken as the updated initial position of the particle, and the target position F point can be taken as the updated end position, so as to obtain the updated first motion trajectory from point E to point F; then, determine the target position that the particle reaches after moving along the Figure 8 new shoreline shown ( Figure 8 the dotted line in), and judge whether the target position is within the shoreline range. If not, repeat the above steps. If so, the target position within the shoreline range is taken as the particle position of the particle at the next time point. Among them, the initial position is the particle position of the particle at the previous time point.

[0113] Through the above technical solution, when the simulated target position is outside the shoreline range, the new target position of the particle can be re-simulated based on the bounce position and the target position until the target position within the shoreline range is obtained. Only in this way will the second motion trajectory obtained based on the starting position of the particle and the target position within the shoreline range be more accurate.

[0114] Figure 9 FIG. is a block diagram of a particle motion trajectory simulation device according to an exemplary embodiment. The particle motion trajectory simulation device 900 includes: a bounce position calculation module 910, a target position calculation module 920, and a trajectory simulation module 930.

[0115] A bounce position calculation module 910, configured to determine a bounce position after the particle collides with the shoreline in the case where a first motion trajectory of the particle intersects the shoreline; A target position calculation module 920, configured to obtain a target position reached by the particle after moving along the shoreline according to the bounce position; A trajectory simulation module 930, configured to obtain a corrected second motion trajectory of the particle according to the starting position and the target position of the particle; the starting position is the starting point of the first motion trajectory.

[0116] Optionally, the target position calculation module 920 is further configured to obtain a first coordinate difference according to a second position and a first position on the shoreline; obtain a second coordinate difference according to the bounce position and the starting position; obtain the target position according to a first coefficient, a second coefficient, the first coordinate difference, and the second coordinate difference; the first coefficient is a coefficient of the first coordinate difference, and the second coefficient is a coefficient of the second coordinate difference.

[0117] Optionally, the first coordinate difference includes a first vertical coordinate difference, and the second coordinate difference includes a second vertical coordinate difference; the target position calculation module 920 is further configured to obtain a first target value according to the first coefficient and the first vertical coordinate difference; obtain a second target value according to the second coefficient and the second vertical coordinate difference; obtain the abscissa in the target position according to the first target value and the second target value.

[0118] Optionally, the first coordinate difference includes a first horizontal coordinate difference, and the second coordinate difference includes a second horizontal coordinate difference; the target position calculation module 920 is further configured to obtain a third target value according to the first coefficient and the first horizontal coordinate difference; obtain a fourth target value according to the second coefficient and the second horizontal coordinate difference; obtain the ordinate in the target position according to the third target value and the fourth target value.

[0119] Optionally, the motion trajectory simulation device 900 of the particle further includes: A first value module, configured to obtain a first value according to a first distance, a second distance, the starting position, and the bounce position; the first distance is the distance from the starting position to the bounce position, the second distance is the distance from the bounce position to an end position, and the end position is the end position on the first motion trajectory; A second value module, configured to obtain a second value according to the first coordinate difference and the second coordinate difference; A first coefficient module, configured to obtain the first coefficient according to the first value and the second value.

[0120] Optionally, the particle motion trajectory simulation device 900 further includes: A third value module, configured to obtain a third value according to the second distance, the third distance, the bounce position, and the first coordinate difference; the second distance is the distance from the bounce position to the end position, the end position is the end position on the first motion trajectory, and the third distance is the distance between the second position and the first position; A second value module, configured to obtain a second value according to the first coordinate difference and the second coordinate difference; A second coefficient module, configured to obtain the second coefficient according to the third value and the second value.

[0121] Optionally, the first position and the second position are the positions of two adjacent shore points on the shoreline.

[0122] Optionally, the particle motion trajectory simulation device 900 further includes: An iteration module, configured to, when the target position is outside the shoreline range, use the bounce position as the initial position of the particle, and re-determine the target position reached by the particle after moving along the shoreline from the updated initial position until the updated target position is within the shoreline range.

[0123] Optionally, the bounce position calculation module 910 is further configured to divide the space where the shoreline is located into a plurality of grids; determine the target grid whose distance from the initial position of the particle is less than a preset distance among the plurality of grids; and when there is a shoreline in the target grid and the first motion trajectory of the particle intersects the shoreline, determine the bounce position after the particle collides with the shoreline.

[0124] Regarding the device in the above embodiments, the specific manners in which each module performs operations have been described in detail in the embodiments related to the method, and will not be elaborated here.

[0125] Figure 10 is a block diagram of an electronic device 1000 shown according to an exemplary embodiment. As Figure 10 shown, the electronic device 1000 may include: a first processor 1001, a first memory 1002. The electronic device 1000 may further include one or more of a multimedia component 1003, a first input / output (I / O) interface 1004, and a first communication component 1005.

[0126] Among them, the first processor 1001 is used to control the overall operation of the electronic device 1000 to complete all or part of the steps in the above-mentioned particle motion trajectory simulation method. The first memory 1002 is used to store various types of data to support the operation of the electronic device 1000. These data may include, for example, instructions for any application or method operating on the electronic device 1000, as well as application-related data, such as contact data, sent and received messages, pictures, audio, video, and so on. The first memory 1002 can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as Static Random Access Memory (SRAM), Electrically Erasable Programmable Read-Only Memory (EEPROM), Erasable Programmable Read-Only Memory (EPROM), Programmable Read-Only Memory (PROM), Read-Only Memory (ROM), magnetic memory, flash memory, a magnetic disk, or an optical disc. The multimedia component 1003 may include a screen and an audio component. The screen may be, for example, a touch screen, and the audio component is used to output and / or input audio signals. For example, the audio component may include a microphone for receiving external audio signals. The received audio signal may be further stored in the first memory 1002 or sent through the first communication component 1005. The audio component further includes at least one speaker for outputting audio signals. The first input / output interface 1004 provides an interface between the first processor 1001 and other interface modules, and the above-mentioned other interface modules may be a keyboard, a mouse, buttons, etc. These buttons may be virtual buttons or physical buttons. The first communication component 1005 is used for wired or wireless communication between the electronic device 1000 and other devices. Wireless communication, such as Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, or 4G, or a combination of one or more of them. Accordingly, the first communication component 1005 may include: a Wi-Fi module, a Bluetooth module, and an NFC module.

[0127] In one exemplary embodiment, the electronic device 1000 may be implemented by one or more application specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components, and is used to execute the above-mentioned method for simulating the motion trajectory of particles.

[0128] In another exemplary embodiment, a computer-readable storage medium including program instructions is further provided. When the program instructions are executed by a processor, the steps of the above-mentioned method for simulating the motion trajectory of particles are implemented. For example, the computer-readable storage medium may be the first memory 1002 including the program instructions, and the above-mentioned program instructions may be executed by the first processor 1001 of the electronic device 1000 to complete the above-mentioned method for simulating the motion trajectory of particles.

[0129] In another exemplary embodiment, a computer program product is further provided. The computer program product includes a computer program that can be executed by a processor. When the computer program is executed by the processor, the steps of the above-mentioned method for simulating the motion trajectory of particles are implemented.

[0130] Figure 11 is a block diagram of an electronic device 1100 shown according to an exemplary embodiment. For example, the electronic device 1100 may be provided as a server. Referring to Figure 11 , the electronic device 1100 includes a second processor 1122, the number of which may be one or more, and a second memory 1132 for storing computer programs that can be executed by the second processor 1122. The computer programs stored in the second memory 1132 may include one or more modules each corresponding to a set of instructions. In addition, the second processor 1122 may be configured to execute the computer program to execute the above-mentioned method for simulating the motion trajectory of particles.

[0131] In addition, the electronic device 1100 may further include a power supply component 1126 and a second communication component 1150. The power supply component 1126 may be configured to perform power management of the electronic device 1100, and the second communication component 1150 may be configured to implement communication of the electronic device 1100, for example, wired or wireless communication. In addition, the electronic device 1100 may further include a second input / output (I / O) interface 1158. The electronic device 1100 may operate based on an operating system stored in the second memory 1132.

[0132] In another exemplary embodiment, there is also provided a computer-readable storage medium including program instructions. When the program instructions are executed by a processor, the steps of the above-described method for simulating the movement trajectory of particles are implemented. For example, the computer-readable storage medium may be the above-described second memory 1132 including program instructions, and the above program instructions may be executed by the second processor 1122 of the electronic device 1100 to complete the above-described method for simulating the movement trajectory of particles.

[0133] In another exemplary embodiment, there is also provided a computer program product. The computer program product includes a computer program capable of being executed by a processor. When the computer program is executed by the processor, the steps of the above-described method for simulating the movement trajectory of particles are implemented.

[0134] The preferred embodiments of the present disclosure have been described in detail above with reference to the accompanying drawings. However, the present disclosure is not limited to the specific details in the above embodiments. Within the scope of the technical concept of the present disclosure, various simple modifications can be made to the technical solutions of the present disclosure, and these simple modifications all fall within the protection scope of the present disclosure.

[0135] In addition, it should be noted that, in the above specific embodiments, the various specific technical features described can be combined in any suitable manner without conflict. To avoid unnecessary repetition, the present disclosure does not separately describe various possible combination manners.

[0136] In addition, any combination can be made between various different embodiments of the present disclosure as long as it does not violate the idea of the present disclosure, and it should also be regarded as the content disclosed by the present disclosure.

Claims

1. A method for simulating the motion trajectory of particles, characterized in that, Including: In the case where the first motion trajectory of the particle intersects the shoreline, determining the rebound position of the particle after colliding with the shoreline; Based on the rebound position, obtaining the target position reached by the particle after moving along the shoreline; Based on the starting position and the target position of the particle, obtaining the corrected second motion trajectory of the particle; the starting position is the starting point of the first motion trajectory.

2. The method according to claim 1, wherein The obtaining the target position reached by the particle after moving along the shoreline based on the rebound position includes: Obtaining a first coordinate difference based on a second position and a first position on the shoreline; Obtaining a second coordinate difference based on the rebound position and the starting position; Obtaining the target position based on a first coefficient, a second coefficient, the first coordinate difference, and the second coordinate difference; the first coefficient is the coefficient of the first coordinate difference, and the second coefficient is the coefficient of the second coordinate difference.

3. The method according to claim 2, characterized in that The first coordinate difference includes a first vertical coordinate difference, and the second coordinate difference includes a second vertical coordinate difference; the obtaining the target position based on a first coefficient, a second coefficient, the first coordinate difference, and the second coordinate difference includes: Obtaining a first target value based on the first coefficient and the first vertical coordinate difference; Obtaining a second target value based on the second coefficient and the second vertical coordinate difference; Obtaining the abscissa in the target position based on the first target value and the second target value.

4. The method according to claim 2, wherein The first coordinate difference includes a first horizontal coordinate difference, and the second coordinate difference includes a second horizontal coordinate difference; the obtaining the target position based on a first coefficient, a second coefficient, the first coordinate difference, and the second coordinate difference includes: Obtaining a third target value based on the first coefficient and the first horizontal coordinate difference; Obtaining a fourth target value based on the second coefficient and the second horizontal coordinate difference; Obtaining the ordinate in the target position based on the third target value and the fourth target value.

5. The method according to any one of claims 2 to 4, characterized in that, The method further includes: Obtaining a first value based on a first distance, a second distance, the starting position, and the rebound position; the first distance is the distance from the starting position to the rebound position, the second distance is the distance from the rebound position to the end position, and the end position is the end position on the first motion trajectory; Obtaining a second value based on the first coordinate difference and the second coordinate difference; Obtaining the first coefficient based on the first value and the second value.

6. The method according to any one of claims 2 to 4, characterized in that The method further includes: Obtaining a third value based on the second distance, a third distance, the rebound position, and the first coordinate difference; the second distance is the distance from the rebound position to the end position, the end position is the end position on the first motion trajectory, and the third distance is the distance between the second position and the first position; Obtaining a second value based on the first coordinate difference and the second coordinate difference; Obtaining the second coefficient based on the third value and the second value.

7. The method according to claim 2, wherein The first position and the second position are the positions of two adjacent shore points on the shoreline.

8. The method according to claim 1, characterized in that, The method further includes: When the target position is outside the shoreline range, use the rebound position as the initial position of the particle, and re-determine the target position reached by the particle after moving along the shoreline starting from the updated initial position until the updated target position is within the shoreline range.

9. The method according to claim 1, wherein When the first movement trajectory of the particle intersects the shoreline, determining the rebound position after the particle collides with the shoreline includes: Dividing the space where the shoreline is located into multiple grids; Determining the target grids in the multiple grids whose distance from the initial position of the particle is less than a preset distance; When there is a shoreline in the target grid and the first movement trajectory of the particle intersects the shoreline, determining the rebound position after the particle collides with the shoreline.

10. An electronic device, characterized in that, including: A memory storing a computer program thereon; A processor for executing the computer program in the memory to implement the steps of the method according to any one of claims 1 to 9.

Images