Computing resource allocation method, device and equipment for detonation tube calculation

Abstract

The invention provides a calculation resource allocation method, device and equipment for detonation tube calculation, and relates to the technical field of computers, and the method comprises the steps: obtaining calculation characteristic parameters and to-be-calculated parameters of a detonation tube, and to-be-allocated calculation resources; wherein the calculation characteristic parameters are used for representing the physical characteristics of the detonation process in the detonation tube calculation; determining an induction area and a non-induction area of the detonation tube according to the calculated characteristic parameters; matching the weights of the sensing area and the non-sensing area when the parameters to be calculated are calculated from the configured comparison tables of different parameters to be calculated, different areas and different weights; and according to the weights of the sensing area and the non-sensing area when the to-be-calculated parameters are calculated and the to-be-allocated computing resources, determining the computing resources allocated to the sensing area and the non-sensing area when the to-be-calculated parameters are calculated. Different computing resources are allocated for the sensing area and the non-sensing area according to the to-be-calculated parameters, and the utilization rate of the computing resources and the efficiency of detonation tube calculation are improved.

Description

Technical Field

[0001] This application relates to the field of computer technology, and more specifically, to a method, apparatus, and device for allocating computing resources for detonation tube calculations. Background Technology

[0002] Detonation tubes are typically millimeter-wide, and require a sufficiently long distance to fully develop their detonation waves after detonation. Therefore, detonation tubes exhibit spatial multi-scale phenomena. Furthermore, the timescale of chemical reactions within the detonation wave is extremely small, with a typical computational step size of approximately 10⁻⁸ seconds, while the time it takes for the detonation wave to traverse the entire computational domain is on the order of milliseconds. Thus, detonation tube computation also exhibits temporal multi-scale phenomena. These spatial and temporal multi-scale phenomena mean that detonation tube computation often consumes significant computational resources. Current detonation tube computations generally employ multi-node parallel methods, dividing the overall computational task into multiple regions, with each node responsible for one region's computation, without allocating computational resources based on the calculated parameters and the corresponding region. This results in insufficient computational resources in the inductive region when calculating some parameters, while non-inductive regions have idle computational resources when calculating the same parameters, leading to not only wasted computational resources but also reduced efficiency in detonation tube computation. Summary of the Invention

[0003] The purpose of this application is to provide a method, apparatus, and device for allocating computing resources for detonation tube calculation, in order to solve the above-mentioned problems existing in the prior art. Different computing resources can be allocated to the sensing area and non-sensing area according to the parameters to be calculated, thereby improving the utilization rate of computing resources and the efficiency of detonation tube calculation.

[0004] Firstly, a method for allocating computational resources for detonation tube calculations is provided, which may include: Obtain the computational characteristic parameters and parameters to be calculated of the detonation tube, as well as the computational resources to be allocated; wherein, the computational characteristic parameters are used to characterize the physical properties of the detonation process in the calculation of the detonation tube; Based on the calculated characteristic parameters, the induction zone and non-induction zone of the detonation tube are determined; From the configured lookup table of different parameters to be calculated, different areas, and different weights, match the weights of the sensing area and the non-sensing area when calculating the parameters to be calculated. Based on the weights of the sensing area and the non-sensing area in calculating the parameter to be calculated, and the computing resources to be allocated, the computing resources to be allocated when calculating the parameter to be calculated are determined.

[0005] In an optional implementation, the calculated characteristic parameters include: the pressure value and airflow velocity at each coordinate point of the detonation tube, and the first speed of sound of the environment in which the detonation tube is located in the ground reference frame. In an optional implementation, the induction zone and non-induction zone of the detonation tube are determined based on the calculated characteristic parameters, including: For any moment after the detonation tube generates a detonation wave, the first detonation wave surface coordinates at that moment are determined based on the comparison between the pressure values ​​at each coordinate point of the detonation tube and the configured pressure threshold. Obtain the wavefront coordinates of the second detonation wave at the previous time. The sound velocity surface coordinates at the given moment are determined based on the first sound velocity, the first detonation wave surface coordinates, and the second detonation wave surface coordinates. The region inside the detonation tube located between the coordinates of the first detonation wave surface and the coordinates of the sound velocity surface is taken as the sensing zone of the detonation tube at that time. The area outside the sensing zone inside the detonation tube is considered as the non-sensing zone of the detonation tube at that time.

[0006] In an optional implementation, the coordinates of the first detonation wave surface at the specified moment are determined based on a comparison between the pressure values ​​at each coordinate point of the detonation tube and the configured pressure threshold, including: For any virtual sampling line inside the detonation tube at that moment, extract the coordinates of the points on the virtual sampling line where the pressure value is greater than the configured pressure threshold to obtain the first target coordinates; wherein, any virtual sampling line is parallel to the axis of the detonation tube and the axis of the detonation tube is taken as the positive direction of the horizontal axis; The first target coordinate point with the largest x-coordinate is selected as the wavefront coordinate of the virtual sampling line; Calculate the arithmetic mean of the wavefront coordinates of each virtual sampling line to obtain the wavefront coordinates of the first detonation wave at that moment.

[0007] In an optional implementation, determining the sound velocity surface coordinates at the given moment based on the first sound velocity, the first detonation wave surface coordinates, and the second detonation wave surface coordinates includes: The ratio of the displacement difference between the first and second detonation wave surface coordinates to the time interval between the stated time and the time before the stated time is taken as the detonation wave velocity. Based on the detonation wave velocity, the first sound velocity is transformed into a coordinate system to obtain the second sound velocity in the detonation wave coordinate system; For any virtual sampling line inside the detonation tube at that moment, extract the coordinate point of the airflow velocity on the virtual sampling line at the second speed of sound to obtain the second target coordinate point; The second target coordinate point with the largest abscissa is selected as the sound velocity surface coordinate of the virtual sampling line; Calculate the arithmetic mean of the sound velocity surface coordinates of each virtual sampling line to obtain the sound velocity surface coordinates at that moment.

[0008] In an optional implementation, the parameters to be calculated include: convection flux, diffusion flux, chemical reaction source term, physical property variables, velocity components, pressure, density, and temperature.

[0009] In an optional implementation, the method further includes: When the parameters to be calculated are chemical reaction source terms, physical property variables, velocity components, pressure, density, and temperature, the weight of the non-inductive region is 0.

[0010] Secondly, a computational resource allocation device for detonation tube calculations is provided, the device comprising: The acquisition unit is used to acquire the calculation characteristic parameters and parameters to be calculated of the detonation tube, as well as the computing resources to be allocated; wherein, the calculation characteristic parameters are used to characterize the physical properties of the detonation process in the calculation of the detonation tube. The determining unit is used to determine the induction zone and non-induction zone of the detonation tube based on the calculated characteristic parameters; The matching unit is used to match the weights of the sensing area and the non-sensing area when calculating the parameters to be calculated from a configured lookup table of different parameters to be calculated, different areas and different weights. An allocation unit is configured to determine the computing resources to be allocated to the sensing area and the non-sensing area when calculating the parameter to be calculated, based on the weights of the sensing area and the non-sensing area in calculating the parameter to be calculated and the computing resources to be allocated.

[0011] Thirdly, an electronic device is provided, which includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; When a processor executes a program stored in memory, it implements any of the steps described in the first aspect above.

[0012] Fourthly, a computer-readable storage medium is provided, wherein a computer program is stored therein, and when executed by a processor, the computer program implements the steps of any of the methods described in the first aspect above.

[0013] This application determines the induction zone and non-induction zone of the detonation tube by calculating characteristic parameters, and determines the weight of the induction zone and non-induction zone when calculating parameters based on different parameters to be calculated, different zones, and different weight comparison tables. This determines the computing resources allocated to the induction zone and non-induction zone when calculating parameters. On the premise that the calculation of the non-induction zone will not affect the calculation of the induction zone and the non-induction zone occupies a large space, the calculation of the non-induction zone is ignored, and more computing resources are allocated to the induction zone. This makes the load of the computing resources to be allocated balanced, thereby accelerating the overall calculation, reducing the amount of calculation, and not affecting the calculation results. Attached Figure Description

[0014] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments of this application will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0015] Figure 1 An architecture diagram of a computing resource allocation system for detonation tube calculation provided in this application embodiment; Figure 2 A flowchart illustrating a computational resource allocation method for detonation tube calculation provided in an embodiment of this application; Figure 3 A schematic diagram of the distribution of sensing and non-sensing areas provided in an embodiment of this application; Figure 4 A schematic diagram of a computing resource allocation device for detonation tube calculation provided in an embodiment of this application; Figure 5 This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Detailed Implementation

[0016] 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 a part of the embodiments of this application, and not all of the embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application. Unless otherwise defined, the technical or scientific terms used in this application should have the ordinary meaning understood by those skilled in the art. The words "first," "second," and similar terms used in this application do not indicate any order, quantity, or importance, but are only used to distinguish different components. The words "comprising" or "including," etc., mean that the element or object preceding the word covers the element or object listed after the word and its equivalents, but do not exclude other elements or objects. The words "connected," "coupled," or "connected," etc., are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. "Up," "down," "left," "right," etc., are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0017] Detonation tubes are a fundamental platform for studying the initiation and propagation characteristics of detonation waves. They are pipes with rectangular or circular cross-sections, filled with a mixture of fuel and oxidizer. After being ignited at one end by a high-energy ignition device, the mixture detonates, generating a detonation wave. The detonation wave continues to propagate along the pipe, igniting the remaining mixture.

[0018] The computational resource allocation method for detonation tube calculation provided in this application embodiment can be applied to... Figure 1 In the system architecture shown, such as Figure 1 As shown, the system may include: a server, multiple terminals with computing resources, and acquisition devices; The server is used to execute the computing resource allocation method for detonation tube calculation provided in the embodiments of this application, and to generate computing tasks based on the parameters to be calculated and the computing resources allocated to the sensing area and non-sensing area, and to distribute them to terminals whose computing resources match the computing resources required by the computing tasks. A terminal is used to receive and execute computing tasks issued by the server. The data acquisition equipment includes: a pressure acquisition device, a sound velocity acquisition device, and an airflow velocity acquisition device; these are used to acquire the pressure and airflow velocity at different coordinate points inside the detonation tube, as well as the first sound velocity of the environment in which the detonation tube is located in the ground reference frame.

[0019] Specifically, the server can be a physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, content delivery networks (CDNs), and big data and artificial intelligence platforms. The terminal can be a user equipment (UE) such as a mobile phone, smartphone, laptop, digital radio receiver, personal digital assistant (PDA), tablet computer (PAD), handheld device, in-vehicle device, wearable device, computing device or other processing device connected to a wireless modem, mobile station (MS), mobile terminal, etc. The terminal and server can be directly or indirectly connected through wired or wireless communication methods, which is not limited herein.

[0020] The preferred embodiments of this application are described below with reference to the accompanying drawings. It should be understood that the preferred embodiments described herein are for illustration and explanation only and are not intended to limit this application. Furthermore, the embodiments and features in the embodiments of this application can be combined with each other without conflict.

[0021] Figure 2 This is a flowchart illustrating a computational resource allocation method for detonation tube calculation provided in an embodiment of this application. Figure 2 As shown, the method may include: Step S210: Obtain the calculation characteristic parameters and parameters to be calculated of the detonation tube, as well as the calculation resources to be allocated.

[0022] The computational characteristic parameters are used to characterize the physical properties of the detonation process in the detonation tube calculation. These parameters quantify the key states of the flow field and environmental references within the detonation tube, including: pressure values ​​and airflow velocities at each coordinate point of the detonation tube, and the first sound velocity of the environment in the ground reference frame. The ground reference frame is an inertial reference frame, a coordinate system established with the ground as the reference. The parameters to be calculated are the fundamental physical quantities to be solved in the numerical calculation of the detonation tube, including: convection flux, diffusion flux, chemical reaction source terms, property variables, velocity components, pressure, density, and temperature. Convection flux includes: mass convection flux, momentum convection flux, energy convection flux, and component convection flux. Diffusion flux includes: momentum diffusion flux, energy diffusion flux, and component diffusion flux. Chemical reaction source terms include: component source terms and energy source terms. Property variables include: thermodynamic properties and transport properties. The computational resources to be allocated are the idle computational resources of each terminal used to calculate the parameters to be calculated in the numerical calculation of the detonation tube. The detonation tube calculation is a CFD simulation of the detonation tube.

[0023] Step S220: Determine the induction zone and non-induction zone of the detonation tube based on the calculated characteristic parameters.

[0024] In this embodiment of the application, the induction zone and non-induction zone of the detonation tube are determined based on calculated characteristic parameters, including: (1) For any moment after the detonation tube generates a detonation wave, the first detonation wave surface coordinates at that moment are determined based on the comparison between the pressure values ​​at each coordinate point of the detonation tube and the configured pressure threshold. The detonation wave is a meteorological detonation wave. The pressure threshold can be 1 MPa. The determination of the first detonation wave surface coordinates at that moment includes: for any virtual sampling line in the detonation tube at that moment, extract the coordinate points on the virtual sampling line whose pressure values ​​are greater than the configured pressure threshold to obtain the first target coordinate point. The virtual sampling line is parallel to the axis of the detonation tube and the axis of the detonation tube is taken as the positive x-axis. With the axis of the detonation tube as the positive x-axis, several virtual sampling lines parallel to x are uniformly arranged in the cross section yz. The first target coordinate point with the largest x-axis is selected as the surface coordinate of the virtual sampling line. The arithmetic mean of the surface coordinates of each virtual sampling line is calculated to obtain the first detonation wave surface coordinates at that moment.

[0025] (2) Obtain the second detonation wave surface coordinates of the previous time of the current time; wherein, at each time after the detonation tube generates the detonation wave, the corresponding first detonation wave surface coordinates are determined and stored as a detonation wave surface time series dataset; obtain the first detonation wave surface coordinates of the previous time of the current time from the detonation wave surface time series dataset and use them as the second detonation wave surface coordinates.

[0026] (3) Determine the sound velocity surface coordinates at this moment based on the first sound velocity, the first detonation wave surface coordinates, and the second detonation wave surface coordinates; wherein, the sound velocity surface is the surface with the detonation wave as the reference frame and the gas velocity after the detonation wave in the detonation tube is the local sound velocity; the sound velocity surface coordinates are the coordinates of the sound velocity surface in the detonation tube. Specifically, the ratio of the displacement difference between the first and second detonation wave surface coordinates to the time interval between the current moment and the previous moment is taken as the detonation wave velocity. Since the propagation velocity of the detonation wave often exhibits unsteady fluctuations after initiation, and considering that instantaneous wave surface jitter may introduce errors, the wave velocity is statistically analyzed at relatively long intervals, such as 50mm. The average wave velocity is determined by the ratio of the displacement difference of the detonation wave surface to the corresponding time interval, thereby suppressing the influence of local disturbances on velocity assessment. Based on the detonation wave velocity, the first sound velocity is transformed to obtain the second sound velocity in the detonation wave coordinate system. The detonation wave coordinate system has its origin at the geometric center of the closed end of the detonation tube (e.g., the ignition end, detonation position), with the positive x-axis along the tube axis (flow direction) and the radial direction as the r-axis (or y-axis). (Axis); For any virtual sampling line inside the detonation tube at that moment, extract the coordinates of the airflow velocity on the virtual sampling line at the second sound speed to obtain the second target coordinate point; select the second target coordinate point with the largest abscissa as the sound speed surface coordinate of the virtual sampling line; calculate the arithmetic mean of the sound speed surface coordinates of each virtual sampling line to obtain the sound speed surface coordinate at that moment.

[0027] (4) The region within the detonation tube located between the coordinates of the first detonation wavefront and the coordinates of the sound velocity surface is taken as the induction zone of the detonation tube at that moment; the other regions within the detonation tube outside the induction zone are taken as the non-induction zone of the detonation tube at that moment; wherein, the induction zone within the detonation tube includes the induction zone and the combustion zone, which is the region between the detonation wavefront (leading shock wave front) and the sound velocity surface (CJ surface); the non-induction zone includes the region after the sound velocity surface or before the detonation wavefront, that is, the detonated or non-detonated region where the chemical reaction is basically completed and the flow field reaches equilibrium; in practical applications, the non-induction zone is taken as a redundant region or the region where the non-induction zone is more than a preset distance threshold from the sound velocity surface or the detonation wavefront is taken as a redundant region, thereby avoiding inaccurate calculations due to inaccurate wavefront or sound velocity surface statistics; the preset distance threshold can be 50 mm; the induction zone and non-induction zone within the detonation tube are as follows Figure 3 As shown.

[0028] Step S230: Match the weights of the sensing area and the non-sensing area when calculating the parameters from the configured tables of different parameters to be calculated, different areas, and different weights.

[0029] The table comparing different parameters to be calculated, different zones, and different weights includes the weights of the sensing zone and the non-sensing zone when calculating different parameters to be calculated.

[0030] In this embodiment of the application, the weights of different regions when calculating different parameters to be calculated are determined based on the following method: The historical computational complexity, total historical computational resources, and historical computational error of the detonation tube when calculating various historical parameters are obtained under different historical division methods. The volume of the sensing and non-sensing zones is different under different historical division methods. The obtained historical data is used to construct a training dataset. The training dataset is input into a pre-built weight prediction model to obtain the weights of different zones when calculating different parameters. The weight prediction model includes: The input layer is used to input the historical computational complexity, total historical computational resources, and historical computational error of the detonation tube when calculating each historical parameter to be calculated using different historical division methods, and to normalize the input data to obtain input features; The LSTM layer consists of two layers; it is used to collect the temporal dependencies of the input features to obtain the hidden state features. The Attention layer is used to calculate attention weights on the hidden state features output by the LSTM layer using an additive attention mechanism, resulting in an attention-weighted feature vector. The fully connected layer consists of two layers. It is used to determine the weights of different regions when calculating different parameters based on the attention-weighted feature vectors, and to use the Sigmoid activation function to ensure that the output weights are in the range of [0,1]. It is also constrained that the sum of the weights of different regions is 1 when calculating any parameter.

[0031] Step S240: Determine the computing resources to be allocated to the sensing area and the non-sensing area when calculating the parameters to be calculated, based on the weights of the sensing area and the non-sensing area in calculating the parameters to be calculated and the computing resources to be allocated.

[0032] When the parameters to be calculated are chemical reaction source terms, physical property variables, velocity components, pressure, density, and temperature, the weight of the non-inductive region is 0.

[0033] Specifically, in the calculation process of the detonation tube, parameters in the non-inductive or redundant regions can be omitted when calculating some parameters, thus greatly accelerating the overall calculation. However, when calculating other parameters, parameters in both the inductive and non-inductive regions must be calculated. Therefore, this application assigns different weights to different parameters and regions. When the non-inductive region does not need to calculate parameters, it is assigned a weight of 0, thus avoiding the allocation of limited computing resources to the non-inductive region, which would lead to idle computing resources, while the computing resource load in the inductive region is large, resulting in full computing efficiency. When both the non-inductive and inductive regions need to calculate parameters, since the computing load and amount of calculation in the inductive region are much greater than those in the non-inductive region, the weight of the non-inductive region is less than that of the inductive region, so that the inductive region can be allocated more computing resources, thereby balancing the load of different computing resources and effectively improving the calculation efficiency of the detonation tube.

[0034] In one embodiment of this application, the pressure threshold is determined by the following method: Obtain the range of fuel characteristic parameters within the detonation tube and the initial pressure range before detonation; wherein, fuel characteristic parameters , R is the activation energy of the reaction, and R is the gas constant. The initial temperature inside the detonation tube; Based on the range of fuel characteristic parameters and the range of initial pressure before detonation, multiple operating condition combinations are obtained; each operating condition combination includes a fuel characteristic parameter and an initial pressure. Obtain the simulated pressure of the detonation tube after detonation simulation under multiple operating conditions; Based on multiple operating condition combinations and corresponding simulated pressures, the dynamic pressure threshold function is fitted using the least squares method, as follows:

[0035] in, Indicates the pressure threshold; This indicates the initial pressure of the detonation tube before it detonates, in MPa. The working condition correction coefficient is obtained by fitting the solution using the least squares method. The pressure values ​​at different coordinate points are extracted from the calculated characteristic parameters of the detonation tube, and the real-time fuel characteristic parameters of the detonation tube are obtained. The pressure values ​​at different coordinate points of the detonation tube and real-time fuel characteristic parameters are input into the dynamic pressure threshold function to obtain the pressure threshold.

[0036] In another embodiment of this application, determining the computing resources allocated to the sensing area and non-sensing area when calculating the parameters to be calculated further includes: Obtain the resource type and resource performance parameters of the computing resources to be allocated; wherein, the resource performance parameters include: the floating-point operation capability, communication latency and load factor of each computing resource to be allocated; when the computing resources to be allocated are heterogeneous resources, their corresponding resource types are different; Based on the resource performance parameters of each computing resource to be allocated, calculate the performance coefficient of each computing resource to be allocated, and sort the computing resources to be allocated in descending order of performance coefficient. Obtain the historical time spent calculating different historical parameters for each computing resource to be allocated; For any historical parameter to be calculated and any resource to be allocated, the ratio of the historical time taken by the resource to be allocated when calculating the historical parameter to be calculated to the minimum historical time taken by different resources to be allocated when calculating the historical parameter to be calculated is used as the matching degree between the resource to be allocated and the historical parameter to be calculated. Based on the sum of the performance coefficients of all computing resources to be allocated, and based on the weights of the sensing area and the non-sensing area in calculating the parameters to be calculated, determine the performance coefficients required by the sensing area and the non-sensing area in calculating the parameters to be calculated. Resources with a matching degree greater than the configured first matching degree threshold are allocated to the sensing area first, until the computing resources allocated to the sensing area are greater than the required computing resources; resources with a matching degree greater than the configured second matching degree threshold are allocated to the non-sensing area, until the computing resources allocated to the non-sensing area are greater than the required computing resources.

[0037] Corresponding to the above method, this application also provides a computing resource allocation device for detonation tube calculation, such as... Figure 4 As shown, the device includes: The acquisition unit 410 is used to acquire the calculation characteristic parameters and parameters to be calculated of the detonation tube, as well as the computing resources to be allocated; wherein, the calculation characteristic parameters are used to characterize the physical properties of the detonation process in the calculation of the detonation tube. The determining unit 420 is used to determine the induction zone and non-induction zone of the detonation tube based on the calculated characteristic parameters; Matching unit 430 is used to match the weights of sensing area and non-sensing area in calculating the parameters from different configured lookup tables of different parameters to be calculated, different areas and different weights. The allocation unit 440 is used to determine the computing resources allocated to the sensing area and the non-sensing area when calculating the parameters to be calculated, based on the weights of the sensing area and the non-sensing area in calculating the parameters to be calculated and the computing resources to be allocated.

[0038] The functions of each functional unit in the computational resource allocation device for detonation tube calculation provided in the above embodiments of this application can be implemented through the above-described method steps. Therefore, the specific working process and beneficial effects of each unit in the computational resource allocation device for detonation tube calculation provided in the embodiments of this application will not be repeated here.

[0039] This application also provides an electronic device, such as... Figure 5 As shown, it includes a processor 510, a communication interface 520, a memory 530, and a communication bus 540, wherein the processor 510, the communication interface 520, and the memory 530 communicate with each other through the communication bus 540.

[0040] Memory 530 is used to store computer programs; When the processor 510 executes the program stored in the memory 530, it performs the following steps: Obtain the computational characteristic parameters and parameters to be calculated of the detonation tube, as well as the computational resources to be allocated; among them, the computational characteristic parameters are used to characterize the physical properties of the detonation process in the calculation of the detonation tube. Based on the calculated characteristic parameters, the induction zone and non-induction zone of the detonation tube are determined; From the different parameters to be calculated, different areas, and different weights in the configuration table, match the weights of the sensing area and the non-sensing area when calculating the parameters to be calculated. Based on the weights of the sensing area and the non-sensing area in calculating the parameters to be calculated, as well as the computing resources to be allocated, the computing resources allocated to the sensing area and the non-sensing area in calculating the parameters to be calculated are determined.

[0041] The communication bus mentioned above can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. This communication bus can be divided into address bus, data bus, control bus, etc. For ease of illustration, only one thick line is used to represent it in the diagram, but this does not mean that there is only one bus or one type of bus.

[0042] The communication interface is used for communication between the aforementioned electronic devices and other devices.

[0043] The memory may include random access memory (RAM) or non-volatile memory (NVM), such as at least one disk storage device. Optionally, the memory may also be at least one storage device located remotely from the aforementioned processor.

[0044] The processors mentioned above can be general-purpose processors, including central processing units (CPUs), network processors (NPs), etc.; they can also be digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0045] The implementation methods and beneficial effects of the various components of the electronic device in the above embodiments for solving the problem can be found in [reference needed]. Figure 2 The steps in the illustrated embodiments are used to implement the electronic device. Therefore, the specific working process and beneficial effects of the electronic device provided in this application will not be repeated here.

[0046] In another embodiment provided in this application, a computer-readable storage medium is also provided, which stores instructions that, when executed on a computer, cause the computer to perform the computing resource allocation method for detonation tube calculation as described in any of the above embodiments.

[0047] In another embodiment provided in this application, a computer program product containing instructions is also provided, which, when run on a computer, causes the computer to execute the computing resource allocation method for detonation tube calculation as described in any of the above embodiments.

[0048] Those skilled in the art will understand that the embodiments in this application can be provided as methods, systems, or computer program products. Therefore, the embodiments in this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the embodiments in this application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0049] This application describes embodiments of methods, apparatus (systems), and computer program products according to embodiments of this application with reference to flowchart illustrations and / or block diagrams. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0050] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0051] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0052] Although preferred embodiments have been described in this application, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of this application.

[0053] Obviously, those skilled in the art can make various modifications and variations to the embodiments of this application without departing from the spirit and scope of the embodiments of this application. Therefore, if these modifications and variations to the embodiments of this application fall within the scope of this application and its equivalents, then these modifications and variations are also intended to be included in the embodiments of this application.

Claims

1. A method for allocating computational resources for detonation tube calculations, characterized in that, The method includes: Obtain the computational characteristic parameters and parameters to be calculated of the detonation tube, as well as the computational resources to be allocated; wherein, the computational characteristic parameters are used to characterize the physical properties of the detonation process in the calculation of the detonation tube; Based on the calculated characteristic parameters, the induction zone and non-induction zone of the detonation tube are determined; From the configured lookup table of different parameters to be calculated, different areas, and different weights, match the weights of the sensing area and the non-sensing area when calculating the parameters to be calculated. Based on the weights of the sensing area and the non-sensing area in calculating the parameter to be calculated, and the computing resources to be allocated, the computing resources to be allocated when calculating the parameter to be calculated are determined.

2. The method as described in claim 1, characterized in that, The calculated characteristic parameters include: the pressure value and airflow velocity at each coordinate point of the detonation tube, and the first sound velocity of the environment in which the detonation tube is located in the ground reference frame.

3. The method as described in claim 2, characterized in that, Based on the calculated characteristic parameters, the induction zone and non-induction zone of the detonation tube are determined, including: For any moment after the detonation tube generates a detonation wave, the first detonation wave surface coordinates at that moment are determined based on the comparison between the pressure values ​​at each coordinate point of the detonation tube and the configured pressure threshold. Obtain the wavefront coordinates of the second detonation wave at the previous time. The sound velocity surface coordinates at the given moment are determined based on the first sound velocity, the first detonation wave surface coordinates, and the second detonation wave surface coordinates. The region inside the detonation tube located between the coordinates of the first detonation wave surface and the coordinates of the sound velocity surface is taken as the sensing zone of the detonation tube at that time. The area outside the sensing zone inside the detonation tube is considered as the non-sensing zone of the detonation tube at that time.

4. The method as described in claim 3, characterized in that, Based on the comparison between the pressure values ​​at each coordinate point of the detonation tube and the configured pressure threshold, the coordinates of the first detonation wave surface at the specified moment are determined, including: For any virtual sampling line inside the detonation tube at that moment, extract the coordinates of the points on the virtual sampling line where the pressure value is greater than the configured pressure threshold to obtain the first target coordinates; wherein, any virtual sampling line is parallel to the axis of the detonation tube and the axis of the detonation tube is taken as the positive direction of the horizontal axis; The first target coordinate point with the largest x-coordinate is selected as the wavefront coordinate of the virtual sampling line; Calculate the arithmetic mean of the wavefront coordinates of each virtual sampling line to obtain the wavefront coordinates of the first detonation wave at that moment.

5. The method as described in claim 3, characterized in that, Determining the sound velocity surface coordinates at the given moment based on the first sound velocity, the first detonation wave surface coordinates, and the second detonation wave surface coordinates includes: The ratio of the displacement difference between the first and second detonation wave surface coordinates to the time interval between the stated time and the time before the stated time is taken as the detonation wave velocity. Based on the detonation wave velocity, the first sound velocity is transformed into a coordinate system to obtain the second sound velocity in the detonation wave coordinate system; For any virtual sampling line inside the detonation tube at that moment, extract the coordinate point of the airflow velocity on the virtual sampling line at the second speed of sound to obtain the second target coordinate point; The second target coordinate point with the largest abscissa is selected as the sound velocity surface coordinate of the virtual sampling line; Calculate the arithmetic mean of the sound velocity surface coordinates of each virtual sampling line to obtain the sound velocity surface coordinates at that moment.

6. The method as described in claim 1, characterized in that, The parameters to be calculated include: convection flux, diffusion flux, chemical reaction source term, physical property variables, velocity component, pressure, density, and temperature.

7. The method as described in claim 6, characterized in that, The method further includes: When the parameters to be calculated are chemical reaction source terms, physical property variables, velocity components, pressure, density, and temperature, the weight of the non-inductive region is 0.

8. A computational resource allocation device for detonation tube calculation, characterized in that, The device includes: The acquisition unit is used to acquire the calculation characteristic parameters and parameters to be calculated of the detonation tube, as well as the computing resources to be allocated; wherein, the calculation characteristic parameters are used to characterize the physical properties of the detonation process in the calculation of the detonation tube. The determining unit is used to determine the induction zone and non-induction zone of the detonation tube based on the calculated characteristic parameters; The matching unit is used to match the weights of the sensing area and the non-sensing area when calculating the parameters to be calculated from a configured lookup table of different parameters to be calculated, different areas and different weights. An allocation unit is configured to determine the computing resources to be allocated to the sensing area and the non-sensing area when calculating the parameter to be calculated, based on the weights of the sensing area and the non-sensing area in calculating the parameter to be calculated and the computing resources to be allocated.

9. An electronic device, characterized in that, The electronic device includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, used to store computer programs; A processor, when executing a program stored in memory, implements the method of any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the method described in any one of claims 1-7.

Images