Particle motion trajectory simulation method and electronic equipment

Abstract

The present disclosure relates to a particle trajectory simulation method and electronic device, and relates to the technical field of pollutant trajectory simulation. The method includes: determining the rebound position of the particle after colliding with the shoreline when the particle's first trajectory intersects the shoreline; obtaining the target position reached by the particle after moving along the shoreline based on the rebound position; and obtaining a corrected second trajectory of the particle based on the particle's starting position and the target position; the starting position is the starting point of the first trajectory. Using the particle trajectory simulation method proposed in the present disclosure, a more accurate particle trajectory can be obtained.

Description

Technical Field

[0001] The present disclosure relates to the technical field of pollutant trajectory simulation, and in particular, to a particle motion trajectory simulation method and electronic equipment. Background Art

[0002] Currently, sudden water pollution accidents occur due to industrial emissions, transportation accidents, and other reasons. Simulating the diffusion and migration of pollutants in water bodies can track the movement of pollutants in water bodies, predict the distribution and concentration changes of pollutants in water bodies, and provide a basis for pollutant control.

[0003] In the related art, when simulating the motion trajectories of particles in pollutants, errors exist in the simulated motion trajectories of the particles. Summary of the Invention

[0004] The purpose of the present disclosure is to provide a particle motion trajectory simulation method and electronic device to solve the above technical problems.

[0005] In order to achieve the above objectives, the present disclosure provides a method for simulating the motion trajectory of particles, comprising:

[0006] When the first motion trajectory of the particle intersects the shoreline, determining a rebound position of the particle after colliding with the shoreline;

[0007] According to the rebound position, a target position reached by the particle after moving along the shoreline is obtained;

[0008] A corrected second motion trajectory of the particle is obtained according to the starting position and the target position of the particle; the starting position is the starting point of the first motion trajectory.

[0009] Optionally, obtaining the target position reached by the particle after moving along the shoreline according to the rebound position includes:

[0010] Obtaining a first coordinate difference based on a second position and a first position on the shoreline;

[0011] Obtaining a second coordinate difference according to the rebound position and the starting position;

[0012] The target position is obtained according to the first coefficient, the 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.

[0013] Optionally, the first coordinate difference includes a first ordinate difference, and the second coordinate difference includes a second ordinate difference; and obtaining the target position according to the first coefficient, the second coefficient, the first coordinate difference, and the second coordinate difference includes:

[0014] Obtaining a first target value according to a difference between the first coefficient and the first ordinate;

[0015] Obtaining a second target value according to a difference between the second coefficient and the second ordinate;

[0016] The horizontal coordinate of the target position is obtained according to the first target value and the second target value.

[0017] Optionally, the first coordinate difference includes a first horizontal coordinate difference, and the second coordinate difference includes a second horizontal coordinate difference; and obtaining the target position according to the first coefficient, the second coefficient, the first coordinate difference, and the second coordinate difference includes:

[0018] Obtaining a third target value according to a difference between the first coefficient and the first horizontal coordinate;

[0019] obtaining a fourth target value according to a difference between the second coefficient and the second abscissa;

[0020] The vertical coordinate of the target position is obtained according to the third target value and the fourth target value.

[0021] Optionally, the method further includes:

[0022] 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;

[0023] Obtaining a second value according to the first coordinate difference and the second coordinate difference;

[0024] The first coefficient is obtained according to the first value and the second value.

[0025] Optionally, the method further includes:

[0026] Obtaining a third value based on a second distance, a third distance, and a difference between the rebound position and the first coordinate; the second distance is the distance from the rebound position to an end position, the end position being the end position on the first motion trajectory; and the third distance is the distance between the second position and the first position;

[0027] Obtaining a second value according to the first coordinate difference and the second coordinate difference;

[0028] The second coefficient is obtained according to the third value and the second value.

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

[0030] Optionally, the method further includes:

[0031] When the target position is outside the coastline, the rebound position is used as the initial position of the particle, and the target position reached by the particle after moving along the coastline from the updated initial position is re-determined until the updated target position is within the coastline.

[0032] Optionally, when the first motion trajectory of the particle intersects a shoreline, determining a rebound position of the particle after colliding with the shoreline includes:

[0033] Dividing the space where the shoreline is located into a plurality of grids;

[0034] determining a target grid among the plurality of grids, the distance between which and the initial position of the particle is less than a preset distance;

[0035] When the coastline exists in the target grid and the first motion trajectory of the particle intersects the coastline, a rebound position of the particle after colliding with the coastline is determined.

[0036] In order to achieve the above objectives, the present disclosure provides an electronic device, comprising:

[0037] a memory having a computer program stored thereon;

[0038] A processor is used to execute the computer program in the memory to implement the steps of the particle motion trajectory simulation method proposed in the present disclosure.

[0039] The above technical solution can determine the rebound position of the particle's first trajectory after colliding with the shoreline, and based on the rebound position, the target position reached by the particle after moving along the shoreline is obtained. Finally, the particle's corrected second trajectory is obtained based on the particle's starting position and target position. In this process, the target position of the particle moving along the longshore current is obtained based on the particle's rebound position, taking into account the influence of the shoreline's longshore current on the particle's trajectory, making the obtained second trajectory more accurate.

[0040] Other features and advantages of the present disclosure will be described in detail in the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The accompanying drawings are used to provide a further understanding of the present disclosure and constitute a part of the specification. Together with the following detailed description, they are used to explain the present disclosure but do not constitute a limitation of the present disclosure. In the accompanying drawings:

[0042] Figure 1 The figure is a flowchart of a method for simulating the motion trajectory of a particle according to an exemplary embodiment.

[0043] Figure 2 The figure is a schematic diagram showing a first motion trajectory and a second motion trajectory of a particle according to an exemplary embodiment.

[0044] Figure 3 The figure is a flowchart of a method for simulating the motion trajectory of a particle according to an exemplary embodiment.

[0045] Figure 4 The figure is a schematic diagram showing a first motion trajectory and a second motion trajectory of a particle according to an exemplary embodiment.

[0046] Figure 5 The figure is a flowchart of a method for simulating the motion trajectory of a particle according to an exemplary embodiment.

[0047] Figure 6 The figure is a flowchart of a method for simulating the motion trajectory of a particle according to an exemplary embodiment.

[0048] Figure 7 The figure is a flowchart of a method for simulating the motion trajectory of a particle according to an exemplary embodiment.

[0049] Figure 8 The figure is a schematic diagram showing a first motion trajectory and a second motion trajectory of a particle according to an exemplary embodiment.

[0050] Figure 9 The figure is a block diagram of a device for simulating the motion trajectory of a particle according to an exemplary embodiment.

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

[0052] Figure 11 It is a block diagram of an electronic device according to an exemplary embodiment. DETAILED DESCRIPTION

[0053] The following describes the specific embodiments of the present disclosure in detail with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and explain the present disclosure and are not intended to limit the present disclosure.

[0054] Floating objects such as pollutants floating in water can be regarded as a collection of particles, and the Lagrangian particle tracking algorithm can be used to predict the motion trajectories of particles in the pollutants.

[0055] The Lagrangian particle tracking algorithm, a crucial tool for simulating the diffusion and migration of pollutants in water, precisely tracks the trajectory of each particle, thereby comprehensively depicting the dynamic behavior of pollutants in water. This algorithm relies not only on a detailed description of the flow field but also on the precise handling of complex shoreline boundary conditions. The shoreline, as the boundary between water and land, has a direct impact on flow patterns, profoundly influencing the diffusion path, deposition rate, and final location of particles.

[0056] When simulating sudden aquatic pollution incidents, shoreline management is crucial. Shorelines vary in form, from straight riverbanks to winding coastlines. These diverse shapes lead to significant differences in water velocity, direction, and vortex structure. This requires the Lagrangian particle tracking algorithm to not only accurately simulate particle motion in the water but also account for complex physical processes such as rebound, adsorption, and deposition when particles interact with the shoreline.

[0057] In related technologies, there are three methods for processing particles on non-transmissive coastlines:

[0058] The first processing method is reflection: particles hit the shoreline and bounce back into the water.

[0059] The second treatment method is absorption: particles collide with the shoreline and then settle or adhere to the shoreline.

[0060] The third processing method is partial reflection: some particles hit the shoreline and bounce back into the water, while other particles hit the shoreline and deposit or adhere to the shoreline.

[0061] The method for processing particles that bounce back into the water after colliding with the shoreline calculates the reflection angle and displacement of the particles after colliding with the shoreline based on the particle's incident angle toward the shoreline, thereby simulating the target position reached by the particles after colliding with the shoreline. However, in reality, due to the influence of complex water currents such as coastal currents, particles do not simply return to the water body according to the reflection angle after colliding with the shoreline. Instead, they undergo a certain displacement along the shoreline. As a result, the target position simulated based on the reflection angle and displacement may be inaccurate, and the particle trajectory derived based on this target position may also be biased.

[0062] Based on this, the present disclosure proposes a method for simulating the motion trajectory of particles. Figure 1 The following is a flow chart of the method for simulating the motion trajectory of particles proposed in the present disclosure, which includes the following steps:

[0063] In step S10 , when the first motion trajectory of the particle intersects with the coastline, a rebound position of the particle after colliding with the coastline is determined.

[0064] The first trajectory of a particle is its initial trajectory, which can also be understood as the trajectory of the particle before correction. This first trajectory has a starting position and an ending position. The starting position is the starting point of the movement on the first trajectory, and the ending position is the ending point of the movement on the first trajectory. The ending position can also be understood as the position that the particle naturally reaches if it is affected by the water flow and does not collide with the shoreline. If the first trajectory intersects the shoreline, the starting position and the ending position will be located on either side of the shoreline, respectively.

[0065] 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 ending position B on the first motion trajectory are located on both sides of the coastline CD.

[0066] Wherein, 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 consists 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, a starting position and an 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. Only then can the first motion trajectory of the particle at each two adjacent time points be corrected to obtain the corrected second motion trajectory of the particle at each two adjacent time points.

[0067] Among them, if the first motion trajectory consists of two particle positions, a starting position and an end position, the end position of the particle can be obtained based on the starting position, the particle's movement speed 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 as the water body flows. For example, the particle's movement distance is obtained based on the particle's movement speed in the water body and the interval between two adjacent time points. The horizontal axis component of the movement distance on the horizontal axis is then superimposed on the horizontal coordinate of the particle's starting position to obtain the horizontal coordinate of the end position. The vertical axis component of the movement distance on the vertical axis is superimposed on the vertical coordinate of the particle's starting position to obtain the vertical coordinate of the end position. The horizontal coordinate of the end position and the vertical coordinate of the end position constitute the end position.

[0068] Alternatively, the intersection point 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.

[0069] The rebound position after the particle collides with the shoreline, the starting position of the particle, and the ending position of the particle are all located on the same straight line. This rebound position is the position where the particle rebounds after colliding with the shoreline without considering the influence of complex water flows such as coastal currents.

[0070] Factors influencing the rebound coefficient include shoreline type and shoreline material properties. Shoreline types include aquatic plants, rocks, and soil. These types have relatively low rebound coefficients, making particles easily attracted to them and less likely to rebound. Therefore, the rebound point after colliding with these types of shorelines is relatively close to the shoreline. Stone, on the other hand, has a relatively high rebound coefficient, making particles easily rebound after colliding with them. Therefore, the rebound point after colliding with these types of shorelines is relatively far from the shoreline. Shoreline material properties include oiliness, viscosity, rigidity, and elasticity. Shorelines with oily or viscous properties have relatively low rebound coefficients, resulting in particles being closer to the shoreline after colliding with them. However, those with rigid or elastic properties have relatively high rebound coefficients, resulting in particles being farther away from the shoreline after colliding with these types of shorelines.

[0071] Optionally, when it is determined that the coastline is closed and the first motion trajectory of the particle intersects the closed coastline, the rebound position of the particle after the collision with the coastline is determined.

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

[0073] When determining whether a shoreline is closed, if the number of acquired data points exceeds a preset number and no data points have the same coordinates among the data points, the line segment formed by the acquired data points is considered to be the shoreline. The absence of data points with the same coordinates among the data points includes: the coordinates of the first and last data points among the data points are different, or the coordinates of the data points other than the first and last data points among the data points are different.

[0074] Since the coastline is a curve or a straight line, the number of data points on the coastline will be greater than the preset number (for example, 3), and the first data point and the last data point on the coastline do not have the same coordinates, and other data points on the coastline except the first data point and the last data point do not have the same coordinates. When these conditions are met, it is considered that the multiple data points obtained can represent a coastline.

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

[0076] Influenced by complex currents such as coastal currents, particles bounce off the shoreline and move downstream. The line connecting the rebound location and the target location is parallel to the shoreline. The target location is the estimated location the particle will reach after colliding with the shoreline and moving along the coastal current.

[0077] Optionally, the target position may be obtained according to the starting position of the particle on the first motion trajectory, the ending position of the first motion trajectory, the first position on the shoreline, the second position on the shoreline, and the rebound position of the particle.

[0078] The first position and the second position on the shoreline can be two adjacent shore points on the shoreline. There are multiple shore points on the shoreline, and the two shore points closest to the initial position of the particle can be selected from the multiple shore points and used as the first position and the second position respectively. For example, see Figure 2 As shown, point C closest to the initial position A of the particle can be screened out from multiple shore point positions as the first position, and point D as the second position.

[0079] It is understandable that although the shoreline of a water body is a closed geometric shape, the shoreline in the acquired water body shoreline data is not a completely closed shoreline, but rather multiple shoreline segments, each of which is composed of multiple shoreline points. When a particle collides with a shoreline, it collides with the shoreline segment that is closest to the particle among the multiple shoreline segments. Therefore, the line connecting the two shoreline points closest to the particle is selected from the multiple shoreline points as the shoreline segment that intersects with the particle. Of course, the line connecting the two shoreline points closest to the rebound position among the multiple shoreline points can also be used as the shoreline segment that intersects with the particle.

[0080] Optionally, the target position may be obtained according to 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.

[0081] In step S30 , a corrected second motion trajectory of the particle is obtained according to the starting position and the target position of the particle.

[0082] Among them, the second motion trajectory is the motion trajectory of the particle after being corrected to take into account the influence of complex water currents such as coastal currents. The second motion trajectory includes two motion trajectories, one of which is from the initial position of the particle toward the shoreline, and the other motion trajectory is the same as the shoreline.

[0083] Alternatively, the line between the starting position and the target position of the particle can be used as the corrected second motion trajectory of the particle. The second motion trajectory is the motion trajectory of the particle in a time period consisting of two adjacent time points.

[0084] Optionally, the starting position, rebound position, and target position of the particle may be connected into a line as the corrected second motion trajectory of the particle.

[0085] Among them, the accurate second motion trajectory should be the motion trajectory from the starting position, passing through the intersection between the first motion trajectory and the shoreline, the rebound position, and then reaching the target position. Although connecting the starting position of the particle to the target position into the second motion trajectory is not very accurate from a microscopic perspective, when viewing the flow trajectory diagram of the pollutant in the water body, it is not necessary to infinitely magnify the flow trajectory diagram to observe the motion trajectory of each particle in the pollutant from a microscopic perspective, but to view the flow trajectory of the pollutant in the water body from a macroscopic perspective. Therefore, even if the starting position of the particle to the target position is connected into the second motion trajectory, the flow trajectory of the pollutant in the water body viewed from a macroscopic perspective is accurate. Of course, if you want to show a more accurate particle flow trajectory diagram, you can also connect the starting position, intersection point, rebound position, and target position into a line to construct an accurate flow trajectory diagram.

[0086] The above technical solution can determine the rebound position of the particle's first trajectory after colliding with the shoreline, and based on the rebound position, the target position reached by the particle after moving along the shoreline is obtained. Finally, the particle's corrected second trajectory is obtained based on the particle's starting position and target position. In this process, the target position of the particle moving along the longshore current is obtained based on the particle's rebound position, taking into account the influence of the shoreline's longshore current on the particle's trajectory, making the obtained second trajectory more accurate.

[0087] Figure 3 This is an exemplary embodiment involved in the above step S20, which is used to explain an exemplary embodiment of obtaining the target position reached by the particle after moving along the coastline according to the rebound position, including the following steps:

[0088] In step S21 , the target position may be obtained according to 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.

[0089] The second distance between the particle's rebound position and the end position is equal to the distance between the particle's rebound position and the target position, so the second distance can be used as the distance between the particle's rebound position and the target position. Figure 4 As shown, taking the rebound position as point E, the end position as point B, and the target position as point F as an example, the second distance between EB is equal to the distance between EF.

[0090] The line connecting the particle's rebound position to the target position is parallel to the shoreline, so the first angle between the shoreline and the horizontal axis can be used as the angle between the line connecting the particle's rebound position to the target position and the horizontal axis. Figure 4 As shown, taking the shoreline as CD and the particle's rebound position to the target position as EF as an example, CD is parallel to EF.

[0091] Therefore, the target position is obtained according to 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. This can be regarded as obtaining the target position according to 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.

[0092] For example, the horizontal component of the second distance projected on the horizontal axis can be obtained based on 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 the horizontal coordinate of the target position can be obtained by superimposing or subtracting the horizontal component on the basis of the horizontal coordinate of the rebound position; the vertical component of the second distance projected on the vertical axis can be obtained based on the second distance between the rebound position of the particle and the end position and the sine value between the shoreline and the horizontal axis, and the vertical coordinate of the target position can be obtained by superimposing or subtracting the vertical component on the basis of the vertical coordinate of the rebound position; finally, the horizontal coordinate and the vertical coordinate of the target position constitute the target position.

[0093] The rebound position of the particle is Figure 4 Point E in the target position is Figure 4 Point F in the middle, the second distance is Figure 4 The first angle between the shoreline and the horizontal axis is L2. Figure 4 in For example, the calculation formulas for the horizontal and vertical coordinates of the target position are as follows:

[0094] (1)

[0095] In formula (1), is the horizontal coordinate of the target position, is the horizontal coordinate of the rebound position, is the ordinate of the target position, is the vertical coordinate of the rebound position, L2 is the second distance between the rebound position and the target position F, and also the distance between the rebound position and the end position B. It is the first angle between the EF segment and the horizontal axis of the coordinate system, and it is also the angle between the coastline and the horizontal axis of the coordinate system.

[0096] 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 said end position, and the first angle between the shoreline and the horizontal axis. The target position of the particle obtained in this way is the target position reached by the particle after passing through the influence of the coastal current of the shoreline. Then, the second motion trajectory obtained based on the initial position of the particle and the target position of the particle will be more accurate.

[0097] Figure 5 This is an exemplary embodiment involved in the above step S20, which is used to explain an exemplary solution for obtaining a target position based on the starting position of the particle on the first motion trajectory, the end position of the first motion trajectory, the first position on the shoreline, the second position on the shoreline, and the rebound position of the particle, including the following steps:

[0098] In step S22, a first coordinate difference is obtained according to the second position and the first position on the shoreline.

[0099] The first coordinate difference includes a first horizontal coordinate difference and a first vertical coordinate difference. The first horizontal coordinate difference can be obtained based on the horizontal coordinate of the second position on the coastline and the horizontal coordinate of the first position; the first vertical coordinate difference can be obtained based on the vertical coordinate of the second position on the coastline and the vertical coordinate of the first position.

[0100] For example, see Figure 2 As shown, the second position is point D ( , ), the first position is point C ( , ), for example, the first coordinate difference includes the first horizontal coordinate difference Difference from the first ordinate .

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

[0102] 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 based on the abscissa of the rebound position and the abscissa of the starting position; the second ordinate difference can be obtained based on the ordinate of the rebound position and the ordinate of the starting position.

[0103] For example, see Figure 2 As shown, the starting position is point A ( , ), the rebound position is point E ( , ), the second coordinate difference includes the second horizontal coordinate difference Difference from the second ordinate .

[0104] 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.

[0105] Optionally, the first coordinate difference may be adjusted using a first coefficient, and the second coordinate difference may be adjusted using a second coefficient to obtain the target position.

[0106] For example, the target position is obtained by adjusting the first horizontal coordinate difference and the first vertical coordinate difference using the first coefficient, and adjusting the second horizontal coordinate difference and the second vertical coordinate difference using the second coefficient. Specifically, the method 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 horizontal coordinate of the target position based on the first target value and the second target value; 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; and obtaining the vertical coordinate of the target position based on the third target value and the fourth target value.

[0107] For example, the calculation formula for obtaining the horizontal and vertical coordinates of the target position is as follows:

[0108] (2)

[0109] 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.

[0110] From the above formula (2), we can see that the first coefficient To adjust the first ordinate difference Get the first target value , using the second coefficient To adjust the second ordinate difference Get the second target value ; Then get the horizontal coordinate of the target position based on the first target value and the second target value ; You can also use the first coefficient To adjust the first horizontal coordinate difference Get the third target value , using the second coefficient To adjust the second horizontal coordinate difference Get the fourth target value , and then obtain the vertical coordinate of the target position according to the third target value and the fourth target value .

[0111] Optionally, for the first coefficient, the following sub-steps A1 to A3 may be used to obtain the first coefficient:

[0112] Sub-step A1: obtaining a first value according to the first distance, the second distance, the starting position and the rebound position.

[0113] 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.

[0114] Alternatively, the first value may 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 rebound position and the abscissa of the starting position, the second ordinate difference between the ordinate of the rebound position and the ordinate of the starting position, the abscissa of the rebound position, and the ordinate of the rebound position. The second angle is the angle between the line connecting the starting position to the rebound position and the line connecting the rebound position to the target position, for example Figure 2 The angle between L1 and L2 is an obtuse angle, so the angle between L1 and L2 is a negative value when it is substituted into the cosine formula. Therefore, the complementary angle of the angle is As the angle between L1 and L2; the first value is the numerator in the first coefficient.

[0115] For example, the calculation formula for obtaining the first value is as follows:

[0116] (3)

[0117] 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 horizontal coordinate of the rebound position, is the vertical coordinate of the rebound position.

[0118] It can be seen from the above formula (3) that the first value can be obtained by superimposing the product between the horizontal coordinate of the rebound position and the difference between the second horizontal coordinates on the basis of the product between the first distance, the second distance and the cosine value of the second angle, and then superimposing the product between the vertical coordinate of the rebound position and the second vertical coordinate.

[0119] Sub-step A2: obtaining a second value according to the first coordinate difference and the second coordinate difference.

[0120] Optionally, a second value may be obtained according to the first abscissa difference, the first ordinate difference, the second abscissa difference, and the second ordinate difference, wherein the second value is the denominator of the first coefficient.

[0121] For example, the calculation formula for obtaining the second value is as follows:

[0122] (4)

[0123] 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.

[0124] 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.

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

[0126] The ratio between the first value and the second value can be used as the first coefficient. With the second value The ratio between them is the first coefficient .

[0127] Optionally, for the second coefficient, the following sub-steps B1 to B3 may be used to obtain the second coefficient:

[0128] Sub-step B1, obtaining a third value according to the second distance, the third distance, the rebound position and the first coordinate difference.

[0129] 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 from the second position to the first position.

[0130] Alternatively, a third value may be obtained according to the second distance, the third distance, the first horizontal coordinate difference, the first vertical coordinate difference, the horizontal coordinate of the rebound position, and the vertical coordinate of the rebound position, wherein the third value is the numerator of the second coefficient.

[0131] For example, the calculation formula for obtaining the third value is as follows:

[0132] (5)

[0133] 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 horizontal coordinate of the rebound position, is the vertical coordinate of the rebound position.

[0134] It can be seen from the above formula (5) that the third value can be obtained by subtracting the product of the difference between the horizontal coordinate of the rebound position and the first horizontal coordinate from the product of the second distance and the third distance, and then subtracting the product of the difference between the vertical coordinate of the rebound position and the first vertical coordinate.

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

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

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

[0138] The ratio between the third value and the second value can be used as the second coefficient. With the second value The ratio between them is the second coefficient .

[0139] In some cases, see Figure 2 As shown, 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 rebound position as point E as an example, the following vector formula can be obtained:

[0140] (6)

[0141] By reorganizing formula (6), we can get the following formula:

[0142]

[0143] (7)

[0144] By further transforming formula (7), we can obtain the above formula (2).

[0145] 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 rebound position after the particle collides 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. The target position of the particle obtained in this way is the target position reached by the particle after passing through the influence of the coastal current of the shoreline. Then, the second motion trajectory obtained based on the initial position of the particle and the target position of the particle will be more accurate.

[0146] Figure 6This is an exemplary embodiment involved in the above step S10, which is used to explain an exemplary solution for obtaining the rebound position of particles after colliding with the shoreline, including the following steps:

[0147] In step S11, the space where the coastline is located is divided into a plurality of grids.

[0148] For example, the space where the coastline is located can be divided into multiple grids such as grid A to grid F, and each grid position represents the position of a sub-area in the space where the coastline is located.

[0149] In step S12, a target grid in the plurality of grids is determined, the distance between the target grid and the initial position of the particle being less than a preset distance.

[0150] Target grids whose distances from the initial positions of the particles are less than a preset distance may be screened out from the multiple grids, and then it is determined whether a shoreline exists in these target grids.

[0151] In step S13 , when the coastline exists in the target grid and the first motion trajectory of the particle intersects the coastline, a rebound position of the particle after colliding with the coastline is determined.

[0152] If there is a coastline in the target grid, it is determined whether the coastline in the target grid intersects with the first motion trajectory of the particle. If so, the rebound position of the particle after colliding with the coastline is determined.

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

[0154] Since the space where the coastline is located includes a large number of coastline segments, if it is calculated whether the multiple coastline segments intersect with the first motion trajectory of the particle, the calculation amount is large.

[0155] Through the above technical solution, the space where the shoreline is located can be divided into multiple grids, and then the target grid whose distance from the initial position of the particle is less than a preset distance is screened out from the multiple grids, and then it is judged whether the shoreline segment in the target grid will intersect with the first motion trajectory of the particle, that is, it is judged whether the shoreline segment in the target grid will intersect with the line connecting the particle positions at two adjacent time points, thereby eliminating the shoreline segments in the grid whose distance from the initial position of the particle is more than the preset distance, without having to involve all shoreline segments in the intersection calculation, thereby reducing the amount of data calculation.

[0156] Figure 7 This is an exemplary embodiment of the present disclosure, which is used to explain a processing solution when the obtained target position is outside the shoreline range, including the following steps:

[0157] In step S40, when the target position is outside the shoreline range, the rebound position is used as the initial position of the particle, and the target position reached by the particle after moving along the shoreline from the updated initial position is re-determined until the updated target position is within the shoreline range.

[0158] The particles in the pollutants drift in the water body, which can also be understood as the particles in the pollutants drifting within the coastline. If the calculated target position is outside the coastline, it means that the simulated target position F still has errors. At this time, the rebound position can be used as the initial position of the particle after update, and the target position can be used as the end position of the particle after update. The first motion trajectory of the particle after update is 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 coastline, the new rebound position of the particle after collision with the updated coastline is determined; then, based on the new rebound position, the target position reached by the particle after moving along the new coastline is obtained, and this process is repeated until the target position obtained is within the coastline range. Then, the starting position of the particle is connected to the target position within the coastline range to obtain the corrected second motion trajectory of the particle.

[0159] See also Figure 8 As shown in the figure, if the target position F is outside the shoreline range, the rebound position E of the particle can be used as the updated initial position of the particle, and the target position F can be used as the updated end position to obtain the updated first motion trajectory from point E to point F. Then determine the particle along Figure 8 The new shoreline shown ( Figure 8 The particle moves to the target position (the dotted line in the figure). The particle determines whether the target position is within the shoreline. If not, the above steps are repeated. If so, the target position within the shoreline is used as the particle's position at the next time point. The initial position is the particle's position at the previous time point.

[0160] 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 rebound position and the target position until the target position within the shoreline range is obtained. Then the second motion trajectory obtained based on the starting position of the particle and the target position within the shoreline range will be more accurate.

[0161] Figure 9 1 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 rebound position calculation module 910 , a target position calculation module 920 and a trajectory simulation module 930 .

[0162] The rebound position calculation module 910 is configured to determine the rebound position of the particle after the particle collides with the shoreline when the first motion trajectory of the particle intersects the shoreline;

[0163] A target position calculation module 920 is configured to obtain a target position reached by the particle after moving along the shoreline according to the rebound position;

[0164] The trajectory simulation module 930 is 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.

[0165] Optionally, the target position calculation module 920 is also configured to obtain a first coordinate difference based on the second position and the first position on the shoreline; obtain a second coordinate difference based on the rebound position and the starting position; obtain the target position based on the first coefficient, the 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.

[0166] Optionally, the first coordinate difference includes a first ordinate difference, and the second coordinate difference includes a second ordinate difference; the target position calculation module 920 is also configured to obtain a first target value based on the first coefficient and the first ordinate difference; obtain a second target value based on the second coefficient and the second ordinate difference; and obtain the horizontal coordinate in the target position based on the first target value and the second target value.

[0167] 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 also configured to obtain a third target value based on the first coefficient and the first horizontal coordinate difference; obtain a fourth target value based on the second coefficient and the second horizontal coordinate difference; and obtain the vertical coordinate in the target position based on the third target value and the fourth target value.

[0168] Optionally, the particle motion trajectory simulation device 900 further includes:

[0169] a first value module configured to obtain 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;

[0170] a second value module configured to obtain a second value according to the first coordinate difference and the second coordinate difference;

[0171] The first coefficient module is configured to obtain the first coefficient according to the first value and the second value.

[0172] Optionally, the particle motion trajectory simulation device 900 further includes:

[0173] a third value module configured to obtain a third value based on a second distance, a third distance, and a difference between the rebound position and the first coordinate; the second distance is the distance from the rebound position to an end position, the end position is an end position on the first motion trajectory, and the third distance is the distance between the second position and the first position;

[0174] a second value module configured to obtain a second value according to the first coordinate difference and the second coordinate difference;

[0175] The second coefficient module is configured to obtain the second coefficient according to the third value and the second value.

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

[0177] Optionally, the particle motion trajectory simulation device 900 further includes:

[0178] The iterative module is configured to use the rebound position as the initial position of the particle when the target position is outside the coastline range, and to re-determine the target position reached by the particle after moving along the coastline starting from the updated initial position until the updated target position is within the coastline range.

[0179] Optionally, the rebound position calculation module 910 is also configured to divide the space where the coastline is located into multiple grids; determine a target grid in the multiple grids whose distance from the initial position of the particle is less than a preset distance; and determine the rebound position of the particle after colliding with the coastline when the coastline exists in the target grid and the first motion trajectory of the particle intersects with the coastline.

[0180] Regarding the apparatus in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

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

[0182] 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. This 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 information, sent and received messages, pictures, audio, video, etc. 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 storage, flash memory, magnetic disk, or optical disk. 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 signals may be further stored in the first memory 1002 or transmitted via the first communication component 1005. The audio component also 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. These 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 may include Wi-Fi, Bluetooth, Near Field Communication (NFC), 2G, 3G, or 4G, or a combination of one or more thereof. Accordingly, the first communication component 1005 may include: a Wi-Fi module, a Bluetooth module, and an NFC module.

[0183] In an 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 to execute the above-mentioned particle motion trajectory simulation method.

[0184] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided. When executed by a processor, the program instructions implement the steps of the above-described particle motion trajectory simulation method. For example, the computer-readable storage medium may be the above-described first memory 1002 including the program instructions. The above-described program instructions may be executed by the first processor 1001 of the electronic device 1000 to implement the above-described particle motion trajectory simulation method.

[0185] 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 particle motion trajectory simulation method are implemented.

[0186] Figure 11 1 is a block diagram of an electronic device 1100 according to an exemplary embodiment. For example, the electronic device 1100 may be provided as a server. Figure 11 The electronic device 1100 includes one or more second processors 1122, and a second memory 1132 for storing a computer program executable by the second processor 1122. The computer program stored in the second memory 1132 may include one or more modules, each corresponding to a set of instructions. Furthermore, the second processor 1122 may be configured to execute the computer program to perform the above-described particle motion trajectory simulation method.

[0187] 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 for the electronic device 1100, and the second communication component 1150 may be configured to enable communication, such as wired or wireless communication, for the electronic device 1100. Furthermore, 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.

[0188] In another exemplary embodiment, a computer-readable storage medium including program instructions is also provided. When executed by a processor, the program instructions implement the steps of the above-described particle motion trajectory simulation method. For example, the computer-readable storage medium may be the above-described second memory 1132 including the program instructions. The above-described program instructions may be executed by the second processor 1122 of the electronic device 1100 to implement the above-described particle motion trajectory simulation method.

[0189] 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 particle motion trajectory simulation method are implemented.

[0190] The preferred embodiments of the present disclosure are described in detail above in conjunction with the accompanying drawings. However, the present disclosure is not limited to the specific details of the above embodiments. Within 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 scope of protection of the present disclosure.

[0191] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any appropriate manner without contradiction. In order to avoid unnecessary repetition, the present disclosure will not further describe various possible combinations.

[0192] In addition, the various embodiments of the present disclosure may be arbitrarily combined, and as long as they do not violate the concept of the present disclosure, they should also be regarded as the contents disclosed by the present disclosure.

Claims

1. A method for simulating particle motion trajectories, characterized in that: include: When the first motion trajectory of the particle intersects the shoreline, determining a rebound position of the particle after colliding with the shoreline; According to the rebound position, a target position reached by the particle after moving along the shoreline is obtained; Obtaining 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; Wherein, obtaining the target position reached by the particle after moving along the shoreline according to the rebound position includes: Obtaining a first coordinate difference based on a second position and a first position on the shoreline, wherein the first position and the second position are positions of two adjacent shore points on the shoreline; Obtaining a second coordinate difference according to the rebound position and the starting position; The target position is obtained according to the first coefficient, the 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.

2. The method according to claim 1, characterized in that The first coordinate difference includes a first ordinate difference, and the second coordinate difference includes a second ordinate difference; and obtaining the target position according to the first coefficient, the second coefficient, the first coordinate difference, and the second coordinate difference includes: Obtaining a first target value according to a difference between the first coefficient and the first ordinate; Obtaining a second target value according to a difference between the second coefficient and the second ordinate; The horizontal coordinate of the target position is obtained according to the first target value and the second target value.

3. The method according to claim 1, characterized in that The first coordinate difference includes a first horizontal coordinate difference, and the second coordinate difference includes a second horizontal coordinate difference; and obtaining the target position according to the first coefficient, the second coefficient, the first coordinate difference, and the second coordinate difference includes: Obtaining a third target value according to a difference between the first coefficient and the first horizontal coordinate; obtaining a fourth target value according to a difference between the second coefficient and the second abscissa; The vertical coordinate of the target position is obtained according to the third target value and the fourth target value.

4. The method according to any one of claims 1 to 3, characterized in that The method further comprises: 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 according to the first coordinate difference and the second coordinate difference; The first coefficient is obtained according to the first value and the second value.

5. The method according to any one of claims 1 to 3, characterized in that The method further comprises: Obtaining a third value based on a second distance, a third distance, and a difference between the rebound position and the first coordinate; the second distance is the distance from the rebound position to an end position, the end position being 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 according to the first coordinate difference and the second coordinate difference; The second coefficient is obtained according to the third value and the second value.

6. The method according to claim 1, characterized in that The method further comprises: When the target position is outside the coastline, the rebound position is used as the initial position of the particle, and the target position reached by the particle after moving along the coastline from the updated initial position is re-determined until the updated target position is within the coastline.

7. The method according to claim 1, characterized in that When the first motion trajectory of the particle intersects the shoreline, determining a 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 a target grid among the plurality of grids, the distance between which and the initial position of the particle is less than a preset distance; When the coastline exists in the target grid and the first motion trajectory of the particle intersects the coastline, a rebound position of the particle after colliding with the coastline is determined.

8. An electronic device, characterized in that: include: a memory having a computer program stored thereon; A processor, configured to execute the computer program in the memory to implement the steps of the method according to any one of claims 1 to 7.

Images