Oil-gas mixture ratio calculation method, device, equipment and medium

Abstract

The application discloses a kind of oil-gas mixture ratio calculation method, device, equipment and medium, it is related to two-phase flow numerical algorithm field, comprising: using KD tree algorithm, geometric positioning algorithm and the current coordinate of droplet in combustion chamber, determine the target control body of droplet in combustion chamber at current time from all control bodies;The next coordinate of droplet in combustion chamber is calculated, so as to determine the target control body at next time, the target control body at next time is updated as the target control body at current time, if first preset condition is not satisfied, then jump to the step of calculating the next coordinate of droplet in combustion chamber again, until first preset condition is satisfied, to obtain the motion trajectory of droplet in combustion chamber;The distribution result of droplet in combustion chamber is calculated using motion trajectory, and the oil-gas mixture ratio in fuel chamber is calculated based on distribution result. It can reduce the amount of calculation, and improve the speed and efficiency of calculating oil-gas mixture ratio.

Description

Technical Field

[0001] This invention relates to the field of numerical algorithms for two-phase flow, and particularly to methods, apparatus, equipment and media for calculating oil-gas mixture ratio. Background Technology

[0002] The air-fuel mixture ratio in the combustion chamber is one of the key factors in improving the performance of aero-engines. Studying the fuel spray atomization mechanism helps to improve fuel atomization by controlling injection conditions. Numerical simulation has significant advantages in in-depth exploration of the fuel atomization mechanism; therefore, efficient atomization numerical algorithms are effective tools for revealing the atomization mechanism. In the study of fuel atomization mechanism, the fuel droplets in the secondary breakup region are small in size and numerous. Using Eulerian meshes to capture droplets would result in a large computational load, making it impractical. To improve the numerical calculation speed, parallel computing is currently widely used for simulation. After droplets are generated in the computational domain, droplet localization must first be performed to determine the correspondence between discrete phase droplets and the mesh. The first drawback is that the most direct method is to judge all meshes one by one based on the droplet coordinates and the spatial range of the mesh, which leads to a large computational load and a significant reduction in calculation speed. The second drawback is that the efficiency of repositioning droplets after position updates is very low.

[0003] In summary, how to reduce the computational load and improve the speed and efficiency of calculating the oil-gas mixture ratio is a problem that needs to be solved in this field. Summary of the Invention

[0004] In view of this, the purpose of this invention is to provide a method, apparatus, device, and medium for calculating the oil-gas mixture ratio, which can reduce the computational workload and improve the speed and efficiency of calculating the oil-gas mixture ratio. The specific solution is as follows:

[0005] In a first aspect, this application discloses a method for calculating the oil-gas mixture ratio, including:

[0006] Using the KD tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber, the target control volume where the droplet in the combustion chamber is located at the current moment is determined from all control volumes;

[0007] The next coordinate of the droplet in the combustion chamber is calculated based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber.

[0008] Using the current coordinates and the next coordinates, the current target control surface is determined from the target control volume at the current moment, and the target control volume at the next moment of the liquid droplet in the combustion chamber is determined based on the current target control surface;

[0009] The target control body at the next moment is updated to the target control body at the current moment, and it is determined whether the first preset condition is met. If not, the process jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber.

[0010] The distribution of droplets in the combustion chamber is calculated using the motion trajectory, and the air-fuel mixture ratio in the fuel chamber is calculated based on the distribution results.

[0011] Optionally, the step of determining the target control volume where the droplet in the combustion chamber is located at the current moment from all control volumes using the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber includes:

[0012] The KD-tree algorithm is used to determine the control volume to be screened from all control volumes, and the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber are used to screen the target control volume where the droplet in the combustion chamber is located at the current moment from the control volume to be screened.

[0013] Optionally, the step of using a geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber to filter out the target control body where the droplet in the combustion chamber is located at the current moment from the control bodies to be screened includes:

[0014] Determine whether the number of the control bodies to be screened is greater than a preset threshold;

[0015] If the value is greater than the target control body, the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber are used to sequentially traverse each control surface of each control body to be screened, so as to screen out the target control body where the droplet in the combustion chamber is located at the current moment.

[0016] Optionally, the step of sequentially traversing each control surface of each of the control bodies to be screened in order to screen out the target control body where the droplet in the combustion chamber is located at the current moment includes:

[0017] Calculate the vector between the droplet in the combustion chamber and the current control surface of the control body to be screened, and obtain the dot product between the normal vector of the droplet in the combustion chamber and the current control surface;

[0018] Determine whether the vector and the dot product are consistent. If they are consistent, determine that the current control body to be screened is the target control body where the droplet in the combustion chamber is located at the current moment.

[0019] Optionally, calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber includes:

[0020] The next coordinate of the droplet in the combustion chamber is calculated using the first motion parameters stored in the target control body at the current time and the second motion parameters stored in the current container corresponding to the target control body at the current time.

[0021] Optionally, before calculating the next coordinate of the droplet in the combustion chamber using the current first motion parameters stored in the target control body at the current time and the current second motion parameters of the droplet in the combustion chamber stored in the current container corresponding to the target control body at the current time, the method further includes:

[0022] Store the current second motion parameter into the current container corresponding to the target control body at the current time.

[0023] Optionally, after determining the target control body where the droplet in the combustion chamber will be located at the next moment based on the current target control surface, the method further includes:

[0024] Determine whether the target control body at the next moment is the same as the target control body at the current moment under the second preset condition. If not, delete the current second motion parameter stored in the current container, obtain the next second motion parameter of the droplet in the combustion chamber, and then save the next second motion parameter to the next container corresponding to the target control body at the next moment.

[0025] Based on preset parallel boundary nodes, determine whether the process at the current moment and the process at the next moment meet preset cross-process communication conditions. If they meet, update the process at the next moment to the process at the current moment.

[0026] Accordingly, updating the target control entity at the next moment to the target control entity at the current moment includes:

[0027] The target control body at the next moment is updated to the target control body at the current moment, and the next second motion parameter and the next container are updated to the current second motion parameter and the current container, respectively.

[0028] Secondly, this application discloses an oil-gas mixture ratio calculation device, comprising:

[0029] The first determining module is used to determine the target control body where the droplet in the combustion chamber is located at the current moment from all control bodies by using the KD tree algorithm, the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber;

[0030] The next coordinate calculation module is used to calculate the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber;

[0031] The second determining module is used to determine the current target control surface from the target control body at the current moment using the current coordinates and the next coordinates, and to determine the target control body at the next moment where the liquid droplet in the combustion chamber is located based on the current target control surface;

[0032] The trajectory acquisition module is used to update the target control body at the next moment to the target control body at the current moment, and determine whether the first preset condition is met. If not, it jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber.

[0033] The mixture ratio calculation module is used to calculate the distribution of liquid droplets in the combustion chamber using the motion trajectory, and to calculate the air-fuel mixture ratio in the fuel chamber based on the distribution results.

[0034] Thirdly, this application discloses an electronic device, comprising:

[0035] Memory, used to store computer programs;

[0036] A processor is used to execute the computer program to implement the steps of the aforementioned disclosed method for calculating the oil-gas mixture ratio.

[0037] Fourthly, this application discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the steps of the aforementioned disclosed method for calculating the oil-gas mixture ratio.

[0038] As can be seen, this application utilizes the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber to determine the target control body where the droplet is located at the current moment from all control bodies; calculates the next coordinate of the droplet based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber; uses the current coordinate and the next coordinate to determine the current target control surface from the target control body at the current moment, and determines the target control body where the droplet is located at the next moment based on the current target control surface; updates the target control body at the next moment to the target control body at the current moment, and determines whether the current condition is met. If not, it jumps back to the step of calculating the next coordinate of the droplet based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber; calculates the distribution result of the droplet in the combustion chamber using the motion trajectory, and calculates the air-fuel mixture ratio in the fuel chamber based on the distribution result. Therefore, this application eliminates the need to judge each control body individually. Instead, it utilizes the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet within the combustion chamber to determine the target control body where the droplet is located at the current moment from among all control bodies. This reduces the computational load and significantly improves the computational speed. By calculating the next coordinate and using it to determine the current control surface (i.e., which side of the target control body the droplet exited from at the current moment), the target control body where the droplet will be located at the next moment can be directly determined without having to judge each control body individually. This greatly improves the efficiency of obtaining the trajectory of the droplet within the combustion chamber and allows for a faster calculation of the air-fuel mixture ratio within the fuel chamber. Attached Figure Description

[0039] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0040] Figure 1 This is a flowchart of a method for calculating the oil-gas mixture ratio disclosed in this application;

[0041] Figure 2 This is a schematic diagram of a specific control body disclosed in this application;

[0042] Figure 3 This is a specific mesh decomposition diagram disclosed in this application;

[0043] Figure 4 This is a specific grid partitioning diagram disclosed in this application;

[0044] Figure 5 This is a flowchart of a specific oil-gas mixture ratio calculation method disclosed in this application;

[0045] Figure 6 This is a schematic diagram of a specific control body and droplet disclosed in this application;

[0046] Figure 7 This is a schematic diagram of a specific droplet motion disclosed in this application;

[0047] Figure 8 This is a flowchart of a specific oil-gas mixture ratio calculation method disclosed in this application;

[0048] Figure 9 This is a schematic diagram of a specific droplet transprocess disclosed in this application;

[0049] Figure 10 This is a schematic diagram of the structure of an oil-gas mixture ratio calculation device disclosed in this application;

[0050] Figure 11 This is a structural diagram of an electronic device disclosed in this application. Detailed Implementation

[0051] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0052] The air-fuel mixture ratio in the combustion chamber is one of the key factors in improving the performance of aero-engines. Studying the fuel spray atomization mechanism helps to improve fuel atomization by controlling injection conditions. Numerical simulation has significant advantages in in-depth exploration of the fuel atomization mechanism; therefore, efficient atomization numerical algorithms are effective tools for revealing the atomization mechanism. In the study of fuel atomization mechanism, the fuel droplets in the secondary breakup region are small in size and numerous. Using Eulerian meshes to capture droplets would result in a large computational load, making it impractical. To improve the numerical calculation speed, parallel computing is currently widely used for simulation. After droplets are generated in the computational domain, droplet localization must first be performed to determine the correspondence between discrete phase droplets and the mesh. The first drawback is that the most direct method is to judge all meshes one by one based on the droplet coordinates and the spatial range of the mesh, which leads to a large computational load and a significant reduction in calculation speed. The second drawback is that the efficiency of repositioning droplets after position updates is very low.

[0053] Therefore, this application provides a corresponding scheme for calculating the oil-gas mixture ratio, which can reduce the amount of calculation and improve the speed and efficiency of calculating the oil-gas mixture ratio.

[0054] See Figure 1 As shown in the figure, this application discloses a method for calculating the oil-gas mixture ratio, including:

[0055] Step S11: Using the KD tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber, determine the target control volume where the droplet in the combustion chamber is located at the current moment from all control volumes.

[0056] In this embodiment, the KD-tree algorithm is first used to determine the n nearest neighboring control bodies to the droplet in the combustion chamber from all control bodies. Due to the irregularity of the grid shape, the nearest neighboring control body may not be the control body where the droplet is located, so n is generally greater than 1. Existing technologies check all control bodies one by one. For example, if there are 100 control bodies (numbers 1 to 100), then when determining the target control body where the droplet is currently located in the combustion chamber, it is necessary to start checking from number 1 until the target control body where the droplet is currently located is determined. If the droplet is currently located in the combustion chamber... If the target control body is number 50, then 50 judgments are required. However, this application uses the KD-tree algorithm to first determine the nearest control bodies to be screened from 100 control bodies, such as control bodies 46 to 51. Then, the geometric positioning algorithm is used to screen the target control body where the droplet in the combustion chamber is located at the current moment from control bodies 46 to 51. Thus, the computational load is reduced from 50 to 6. In reality, the number of all control bodies may be very large, but the number of nearest control bodies to be screened is generally small. Therefore, the computational load is greatly reduced and the computational speed is significantly improved.

[0057] In this embodiment, before determining the target control volume where the droplet in the combustion chamber is currently located from all control volumes using the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber, the target mesh needs to be decomposed accordingly. For example... Figure 2 The diagram shows a specific control volume. The dashed area of ​​grid node 11 represents control volume 11, and the black dot represents a droplet within the combustion chamber. An Euler-Lagrange solver is implemented using the MPI (Mean Interface) programming environment, and the target mesh is decomposed using the Euler-Lagrange solver. For example... Figure 3The diagram illustrates a specific mesh decomposition method, dividing 42 mesh nodes into four regions. Each of these four regions is processed by a different CPU (central processing unit) process. Since determining the current location of the droplet within the combustion chamber requires data from the next adjacent mesh node, the adjacent mesh nodes at the boundary points must also be included in the partitioning process. Figure 4 The diagram illustrates a specific grid partitioning process. For example, calculating the value of grid node 4 requires data from node 5, and vice versa. The target grid is decomposed into regions. Adjacent grid sub-regions communicate via parallel boundary nodes. These parallel boundary nodes are both external nodes of the current partition and internal nodes of their corresponding adjacent sub-regions. Therefore, when a droplet is located at a parallel boundary node in this partition, its second velocity parameter is transferred to the partition where its control volume is an internal node to complete droplet communication, enabling subsequent calculation of the droplet's next coordinate within the combustion chamber. It is important to note that, because calculations need to be performed using the process corresponding to the target control body where the droplet is currently located in the combustion chamber, the corresponding process must first be determined. First, the process where the nozzle is located is determined based on the nozzle coordinates and defined as the injection process. The nozzle injection range is calculated based on the nozzle model to determine the process that the initially injected droplet may enter. Second, the injection process executes droplet injection based on the nozzle model, distributing the injected droplets to the possible processes. The droplets are located based on the initial droplet coordinates and a rapid droplet positioning method. If the positioning control body has parallel boundary nodes and the final internal nodes are not the control body where the droplet is located, then the droplet does not belong to the current process. If the droplet control body is an internal node of the current process, then the droplet belongs to the current process. The droplet is then stored in the corresponding container of the corresponding control body, completing the distribution of the injected droplets, i.e., the droplet positioning.

[0058] Step S12: Calculate the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber.

[0059] In this embodiment, when tracking the droplet trajectory, the droplet's coordinates are first obtained, and then the control volume where the droplet is located is determined based on the coordinates. It is important to note that, in order to accurately calculate the droplet's next coordinate, the first motion parameters within the control volume also affect the droplet's motion. For example, if the air within the target control volume at the current moment moves to the left, it will exert a force to the left on the droplet. The force may also vary depending on the temperature. Therefore, the first motion parameters may include the control volume's temperature and velocity parameters, while the second motion parameters may include the droplet's velocity, temperature, volume, and mass parameters.

[0060] Step S13: Using the current coordinates and the next coordinates, determine the current target control surface from the target control body at the current moment, and determine the target control body at the next moment where the droplet in the combustion chamber is located based on the current target control surface.

[0061] In this embodiment, it is understood that the coordinates of the droplet can be obtained first, and then the control body where the droplet is located can be determined. In the prior art, when determining the target control body where the droplet is located in the combustion chamber at the next moment, it is still necessary to re-determine whether each control body is the target control body where the droplet is located in the combustion chamber at the next moment. For example, if there are 100 control bodies in total, and the target control body at the next moment is number 51, then it is necessary to make 51 determinations to determine it, which is computationally intensive and inefficient. This application determines the current target control surface from the target control body at the current moment based on the current coordinates and the next coordinates, that is, it determines which side of the target control body the droplet will pass through to reach the target control body at the next moment. The amount of computation is greatly reduced and the computational efficiency is significantly improved. The time interval between the previous and next moments is the computation time step. To accurately track the droplet's trajectory, the computation time step is controlled so that the fastest-moving droplet cannot cross the smallest control volume within the grid. For example, if the target control volume where the droplet was located at moment 10 at the previous moment, then the target control volume where the droplet will be located at moment 2, 3, 4, 9, 10, 11, 16, 17, or 18 at the next moment. When the droplet's position changes due to motion, the droplet repositioning first determines whether the droplet is still in the current control volume using a geometric positioning method. If it is, the repositioning is complete. If not, the droplet's position at the two moments determines which face of the control volume the droplet exits from. The droplet's latest control volume is updated to the control volume sharing the exit control face, and the repositioning is complete.

[0062] Step S14: Update the target control body at the next moment to the target control body at the current moment, and determine whether the first preset condition is met. If not, jump back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber.

[0063] In this embodiment, the first preset condition can be a pre-set time node, such as whether the current time has reached a pre-set time node. If it has, the process stops; if it has not, the process jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber. The first preset condition can also be a pre-set distance point, such as determining whether the current coordinate is not less than a pre-set distance point. If it is not less than a pre-set distance point, the process stops; if it is less than a pre-set distance point, the process jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber.

[0064] Step S15: Calculate the distribution of droplets in the combustion chamber using the motion trajectory, and calculate the air-fuel mixture ratio in the fuel chamber based on the distribution results.

[0065] In this embodiment, since the movement trajectory of each droplet in the combustion chamber is obtained, the distribution result of the droplets in the combustion chamber can be calculated using the movement trajectory, and then the fuel-air mixture ratio of the droplets in the fuel chamber and the air in the fuel chamber can be obtained.

[0066] As can be seen, this application utilizes the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber to determine the target control body where the droplet is located at the current moment from all control bodies; calculates the next coordinate of the droplet based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber; uses the current coordinate and the next coordinate to determine the current target control surface from the target control body at the current moment, and determines the target control body where the droplet is located at the next moment based on the current target control surface; updates the target control body at the next moment to the target control body at the current moment, and determines whether the current condition is met. If not, it jumps back to the step of calculating the next coordinate of the droplet based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber; calculates the distribution result of the droplet in the combustion chamber using the motion trajectory, and calculates the air-fuel mixture ratio in the fuel chamber based on the distribution result. Therefore, this application eliminates the need to judge each control body individually. Instead, it utilizes the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet within the combustion chamber to determine the target control body where the droplet is located at the current moment from among all control bodies. This reduces the computational load and significantly improves the computational speed. By calculating the next coordinate and using it to determine the current control surface (i.e., which side of the target control body the droplet exited from at the current moment), the target control body where the droplet will be located at the next moment can be directly determined without having to judge each control body individually. This greatly improves the efficiency of obtaining the trajectory of the droplet within the combustion chamber and allows for a faster calculation of the air-fuel mixture ratio within the fuel chamber.

[0067] See Figure 5 As shown in the figure, this application discloses a specific method for calculating the oil-gas mixture ratio, including:

[0068] Step S21: Use the KD tree algorithm to determine the control body to be screened from all control bodies, and use the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber to screen the target control body where the droplet in the combustion chamber is located at the current moment.

[0069] In this embodiment, the step of using a geometric positioning algorithm and the current coordinates of the droplets in the combustion chamber to filter out the target control body where the droplets in the combustion chamber are located at the current moment specifically includes: determining whether the number of control bodies to be filtered is greater than a preset threshold; if it is greater, then using a geometric positioning algorithm and the current coordinates of the droplets in the combustion chamber, sequentially traversing each control surface of each of the control bodies to be filtered out in order to filter out the target control body where the droplets in the combustion chamber are located at the current moment.

[0070] In this embodiment, the step of sequentially traversing each control surface of each of the control bodies to be screened in order to screen out the target control body where the droplet in the combustion chamber is located at the current moment specifically includes: calculating the vector between the droplet in the combustion chamber and the current control surface of the current control body to be screened, and obtaining the dot product between the normal vector of the droplet in the combustion chamber and the normal vector of the current control surface; determining whether the vector and the dot product are consistent, and if they are consistent, determining that the current control body to be screened is the target control body where the droplet in the combustion chamber is located at the current moment. For example Figure 6 The diagram shows a specific control body and droplet. Vp is the vector of the droplet in the combustion chamber and the current control surface of the control body to be screened, and n is the normal vector of the droplet in the combustion chamber and the current control surface.

[0071] Step S22: Calculate the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber.

[0072] In this embodiment, Figure 7 The diagram illustrates a specific droplet motion, where P2 represents the current droplet position, P1 represents the next droplet position, and the intersection point r of the droplet's trajectory and the control surface is determined using the principle of similar triangles. c .

[0073] Step S23: Using the current coordinates and the next coordinates, determine the current target control surface from the target control body at the current moment, and determine the target control body at the next moment where the droplet in the combustion chamber is located based on the current target control surface.

[0074] Step S24: Update the target control body at the next moment to the target control body at the current moment, and determine whether the first preset condition is met. If not, jump back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber.

[0075] Step S25: Calculate the distribution of droplets in the combustion chamber using the motion trajectory, and calculate the air-fuel mixture ratio in the fuel chamber based on the distribution results.

[0076] Therefore, this application uses the KD tree algorithm to determine the control body to be screened from all control bodies, and uses the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber to screen the target control body where the droplet in the combustion chamber is located at the current moment. It does not require judging all control bodies one by one, which greatly reduces the workload and significantly improves the calculation speed.

[0077] See Figure 8 As shown in the figure, this application discloses a specific method for calculating the oil-gas mixture ratio, including:

[0078] Step S31: Using the KD tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber, determine the target control volume where the droplet in the combustion chamber is located at the current moment from all control volumes.

[0079] Step S32: Calculate the next coordinate of the droplet in the combustion chamber using the current first motion parameters stored in the target control body at the current time and the current second motion parameters of the droplet in the current container corresponding to the target control body at the current time.

[0080] In this embodiment, before calculating the next coordinate of the droplet in the combustion chamber using the current first motion parameters stored in the target control body at the current time and the current second motion parameters of the droplet in the combustion chamber stored in the current container corresponding to the target control body at the current time, the method further includes: storing the current second motion parameters in the current container corresponding to the target control body at the current time. The storage method for the droplet's second motion parameters is to establish a linked list data structure container. The length of the container can be freely increased or decreased. In order to facilitate the exchange of mass, momentum, and energy between the droplet and the gas phase, the container corresponds to the control body, that is, the droplet is stored in the container of its respective control body.

[0081] Step S33: Using the current coordinates and the next coordinates, determine the current target control surface from the target control body at the current moment, and determine the target control body at the next moment where the droplet in the combustion chamber is located based on the current target control surface.

[0082] In this embodiment, after determining the target control body where the droplet in the combustion chamber is located at the next moment based on the current target control surface, the method further includes: determining whether the target control body at the next moment and the target control body at the current moment meet a second preset condition; if not, deleting the current second motion parameters stored in the current container, obtaining the next second motion parameters of the droplet in the combustion chamber, and then saving the next second motion parameters to the next container corresponding to the target control body at the next moment; based on a preset parallel boundary node, determining whether the process at the current moment and the process at the next moment meet a preset cross-process communication condition; if so, updating the process at the next moment to the process at the current moment. The second preset condition is whether the target control body at the current moment and the target control body at the next moment are the same control body; if so, the second preset condition is met; otherwise, it is not met. It should be noted that it is also necessary to determine whether the current process and the next process are the same process, for example... Figure 9 The diagram illustrates a specific droplet process. Control bodies 1 to 4, 8 to 11, and 15 to 18 constitute the first process, while control bodies 5 to 7, 12 to 14, and 19 to 21 constitute the second process. The target control body at the current moment is 11, meaning the current process is the first process. The target control body at the next moment is 12, meaning the current process is the second process. Therefore, the current second motion parameters stored in the current container corresponding to target control body 11 need to be forwarded to the current container corresponding to target control body 12 so that the second process can perform relevant calculations. Furthermore, the current second motion parameters in the current container corresponding to target control body 11 need to be deleted, i.e., redundant data needs to be removed.

[0083] Step S34: Update the target control body at the next moment to the target control body at the current moment, and determine whether the first preset condition is met. If not, jump back to the step of calculating the next coordinate of the droplet in the combustion chamber using the current first motion parameters stored in the current container corresponding to the target control body at the current moment and the current second motion parameters of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber.

[0084] In this embodiment, updating the target control body at the next moment to the target control body at the current moment specifically includes: updating the target control body at the next moment to the target control body at the current moment, and updating the next second motion parameter and the next container to the current second motion parameter and the current container, respectively.

[0085] Step S35: Calculate the distribution of droplets in the combustion chamber using the motion trajectory, and calculate the air-fuel mixture ratio in the fuel chamber based on the distribution results.

[0086] Therefore, when this application uses a container to store the second motion parameters of a droplet, it needs to be based on the target control body where the droplet is located at the current moment. If the target control body at the current moment is different from the target control body at the next moment, the current second motion parameters stored in the container are deleted, that is, redundant data needs to be deleted to save resource space. Moreover, communication only needs to be carried out between the target control body at the current moment and the target control body at the next moment, without global communication, thereby reducing communication time.

[0087] See Figure 10 As shown in the figure, this application discloses an oil-gas mixture ratio calculation device, including:

[0088] The first determining module 11 is used to determine the target control body where the droplet in the combustion chamber is located at the current moment from all control bodies by using the KD tree algorithm, the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber;

[0089] The next coordinate calculation module 12 is used to calculate the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber;

[0090] The second determining module 13 is used to determine the current target control surface from the target control body at the current moment using the current coordinates and the next coordinates, and to determine the target control body at the next moment where the liquid droplet in the combustion chamber is located based on the current target control surface;

[0091] The trajectory acquisition module 14 is used to update the target control body at the next moment to the target control body at the current moment, and determine whether the first preset condition is met. If not, it jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber.

[0092] The mixture ratio calculation module 15 is used to calculate the distribution result of liquid droplets in the combustion chamber using the motion trajectory, and to calculate the air-fuel mixture ratio in the fuel chamber based on the distribution result.

[0093] As can be seen, this application utilizes the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber to determine the target control body where the droplet is located at the current moment from all control bodies; calculates the next coordinate of the droplet based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber; uses the current coordinate and the next coordinate to determine the current target control surface from the target control body at the current moment, and determines the target control body where the droplet is located at the next moment based on the current target control surface; updates the target control body at the next moment to the target control body at the current moment, and determines whether the current condition is met. If not, it jumps back to the step of calculating the next coordinate of the droplet based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber; calculates the distribution result of the droplet in the combustion chamber using the motion trajectory, and calculates the air-fuel mixture ratio in the fuel chamber based on the distribution result. Therefore, this application eliminates the need to judge each control body individually. Instead, it utilizes the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet within the combustion chamber to determine the target control body where the droplet is located at the current moment from among all control bodies. This reduces the computational load and significantly improves the computational speed. By calculating the next coordinate and using it to determine the current control surface (i.e., which side of the target control body the droplet exited from at the current moment), the target control body where the droplet will be located at the next moment can be directly determined without having to judge each control body individually. This greatly improves the efficiency of obtaining the trajectory of the droplet within the combustion chamber and allows for a faster calculation of the air-fuel mixture ratio within the fuel chamber.

[0094] In some specific embodiments, the first determining module 11 includes:

[0095] The target control volume determination unit is used to determine the control volume to be screened from all control volumes using the KD tree algorithm, and to screen the target control volume where the droplet in the combustion chamber is located at the current moment from the control volumes to be screened using the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber.

[0096] In some specific embodiments, the target control body determining unit includes:

[0097] The first judgment unit is used to determine whether the number of the control bodies to be screened is greater than a preset threshold.

[0098] The filtering unit is used to sequentially traverse each control surface of each of the control bodies to be filtered by using a geometric positioning algorithm and the current coordinates of the droplets in the combustion chamber if the target control body is larger than the target control body at the current moment.

[0099] In some specific embodiments, the filtering unit includes:

[0100] The vector calculation unit is used to calculate the vector between the liquid droplet in the combustion chamber and the current control surface of the current control body to be screened;

[0101] The dot product acquisition unit is used to acquire the dot product between the liquid droplet in the combustion chamber and the normal vector of the current control surface;

[0102] The second judgment unit is used to determine whether the vector is consistent with the dot product. If they are consistent, the current control body to be screened is determined to be the target control body where the droplet in the combustion chamber is located at the current moment.

[0103] In some specific embodiments, the next coordinate calculation module 12 includes:

[0104] The droplet next coordinate calculation unit is used to calculate the next coordinate of the droplet in the combustion chamber by using the current first motion parameters stored in the target control body at the current time and the current second motion parameters of the droplet in the combustion chamber stored in the current container corresponding to the target control body at the current time.

[0105] In some specific embodiments, the droplet lower coordinate calculation unit includes:

[0106] The first storage unit is used to store the current second motion parameter into the current container corresponding to the target control body at the current time.

[0107] In some specific embodiments, the second determining module 13 includes:

[0108] The second storage unit is used to determine whether the target control body at the next moment and the target control body at the current moment meet the second preset condition; if not, the current second motion parameter stored in the current container is deleted, and the next second motion parameter of the droplet in the combustion chamber is obtained, and then the next second motion parameter is saved to the next container corresponding to the target control body at the next moment.

[0109] The process update unit is used to determine, based on a preset parallel boundary node, whether the process at the current time and the process at the next time time meet the preset cross-process communication conditions. If they are met, the process at the next time time is updated to the process at the current time time.

[0110] In some specific embodiments, the trajectory acquisition module 14 includes:

[0111] The update unit is used to update the target control body at the next moment to the target control body at the current moment, and to update the next second motion parameter and the next container to the current second motion parameter and the current container, respectively.

[0112] Figure 11 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Specifically, it may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the oil-gas mixture ratio calculation method performed by the electronic device disclosed in any of the foregoing embodiments.

[0113] In this embodiment, the power supply 23 is used to provide operating voltage for various hardware devices on the electronic device; the communication interface 24 can create a data transmission channel between the electronic device and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0114] The processor 21 may include one or more processing cores, such as a quad-core processor or an octa-core processor. The processor 21 may be implemented using at least one hardware form selected from DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). The processor 21 may also include a main processor and a coprocessor. The main processor, also known as a CPU (Central Processing Unit), is used to process data in the wake-up state; the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, the processor 21 may integrate a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the screen. In some embodiments, the processor 21 may also include an AI (Artificial Intelligence) processor, which is used to handle computational operations related to machine learning.

[0115] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored on it include operating system 221, computer program 222 and data 223, etc., and the storage method can be temporary storage or permanent storage.

[0116] The operating system 221 manages and controls the various hardware devices and computer programs 222 on the electronic device to enable the processor 21 to perform calculations and processing on the massive amounts of data 223 in the memory 22. The operating system can be Windows, Unix, Linux, etc. The computer program 222, in addition to including a computer program capable of performing the oil-gas mixture ratio calculation method executed by the electronic device as disclosed in any of the foregoing embodiments, may further include computer programs capable of performing other specific tasks. The data 223 may include data received by the electronic device from external devices, as well as data collected by its own input / output interface 25.

[0117] Furthermore, embodiments of this application also disclose a computer-readable storage medium storing a computer program. When the computer program is loaded and executed by a processor, it implements the method steps performed during the calculation of the oil-gas mixture ratio disclosed in any of the foregoing embodiments.

[0118] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0119] The above provides a detailed description of the oil-gas mixture ratio calculation method, apparatus, equipment, and medium provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for calculating the oil-gas mixture ratio, characterized in that, include: Using the KD tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber, the target control volume where the droplet in the combustion chamber is located at the current moment is determined from all control volumes; The next coordinate of the droplet in the combustion chamber is calculated based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber. Using the current coordinates and the next coordinates, the current target control surface is determined from the target control volume at the current moment, and the target control volume at the next moment of the liquid droplet in the combustion chamber is determined based on the current target control surface; The target control body at the next moment is updated to the target control body at the current moment, and it is determined whether the first preset condition is met. If not, the process jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber. The distribution of liquid droplets in the combustion chamber is calculated using the motion trajectory, and the air-fuel mixture ratio in the fuel chamber is calculated based on the distribution results. The step of determining the target control volume where the droplet in the combustion chamber is located at the current moment from all control volumes using the KD-tree algorithm, geometric positioning algorithm, and the current coordinates of the droplet in the combustion chamber includes: The KD tree algorithm is used to determine the n control bodies that are closest to the droplet in the combustion chamber from all control bodies. Then, the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber are used to select the target control body where the droplet in the combustion chamber is located at the current moment from the control bodies to be screened. The step of using a geometric positioning algorithm and the current coordinates of the droplets inside the combustion chamber to filter out the target control volume where the droplets are currently located from the control volumes to be screened includes: Determine whether the number of the control bodies to be screened is greater than a preset threshold; if it is, use a geometric positioning algorithm and the current coordinates of the droplets in the combustion chamber to traverse each control surface of each control body to be screened in order to screen out the target control body where the droplets in the combustion chamber are located at the current moment. The step of sequentially traversing each control surface of each of the control bodies to be screened in order to screen out the target control body where the droplet in the combustion chamber is located at the current moment includes: Calculate the vector between the droplet in the combustion chamber and the current control surface of the control body to be screened, and obtain the dot product between the normal vector of the droplet in the combustion chamber and the normal vector of the current control surface; determine whether the vector and the dot product are consistent, and if they are consistent, determine that the control body to be screened is the target control body where the droplet in the combustion chamber is located at the current moment.

2. The method for calculating the oil-gas mixture ratio according to claim 1, characterized in that, The calculation of the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber includes: The next coordinate of the droplet in the combustion chamber is calculated using the first motion parameters stored in the target control body at the current time and the second motion parameters stored in the current container corresponding to the target control body at the current time.

3. The method for calculating the oil-gas mixture ratio according to claim 2, characterized in that, Before calculating the next coordinate of the droplet in the combustion chamber using the current first motion parameters stored in the target control body at the current time and the current second motion parameters of the droplet in the combustion chamber stored in the current container corresponding to the target control body at the current time, the method further includes: Store the current second motion parameter into the current container corresponding to the target control body at the current time.

4. The method for calculating the oil-gas mixture ratio according to claim 3, characterized in that, After determining the target control body where the droplet in the combustion chamber will be located at the next moment based on the current target control surface, the method further includes: Determine whether the target control body at the next moment and the target control body at the current moment meet the second preset condition. If not, delete the current second motion parameter stored in the current container, obtain the next second motion parameter of the droplet in the combustion chamber, and then save the next second motion parameter to the next container corresponding to the target control body at the next moment. Based on preset parallel boundary nodes, determine whether the process at the current moment and the process at the next moment meet preset cross-process communication conditions. If they meet, update the process at the next moment to the process at the current moment. Accordingly, updating the target control entity at the next moment to the target control entity at the current moment includes: The target control body at the next moment is updated to the target control body at the current moment, and the next second motion parameter and the next container are updated to the current second motion parameter and the current container, respectively.

5. A device for calculating the oil-gas mixture ratio, characterized in that, The steps for implementing the oil-gas mixture ratio calculation method as described in any one of claims 1 to 4 include: The first determining module is used to determine the target control body where the droplet in the combustion chamber is located at the current moment from all control bodies by using the KD tree algorithm, the geometric positioning algorithm and the current coordinates of the droplet in the combustion chamber; The next coordinate calculation module is used to calculate the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current time and the current second motion parameter of the droplet in the combustion chamber; The second determining module is used to determine the current target control surface from the target control body at the current moment using the current coordinates and the next coordinates, and to determine the target control body at the next moment where the liquid droplet in the combustion chamber is located based on the current target control surface; The trajectory acquisition module is used to update the target control body at the next moment to the target control body at the current moment, and determine whether the first preset condition is met. If not, it jumps back to the step of calculating the next coordinate of the droplet in the combustion chamber based on the current first motion parameter corresponding to the target control body at the current moment and the current second motion parameter of the droplet in the combustion chamber, until the first preset condition is met, so as to obtain the motion trajectory of the droplet in the combustion chamber. The mixture ratio calculation module is used to calculate the distribution of liquid droplets in the combustion chamber using the motion trajectory, and to calculate the air-fuel mixture ratio in the fuel chamber based on the distribution results.

6. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the steps of the oil-gas mixture ratio calculation method as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, Used to store a computer program; wherein, when the computer program is executed by a processor, it implements the steps of the oil-gas mixture ratio calculation method as described in any one of claims 1 to 4.

Images