A flow field shock wave detection method, device, equipment and storage medium

Abstract

The present application discloses a flow field shock wave detection method, device, equipment and storage medium, which relates to the field of fluid mechanics technology, including: using a preset simulation tool to simulate the flow field distribution of a target device under a target operating environment to obtain the flow field variable value of each flow field grid area corresponding to the target device; determining the pressure gradient of the target flow field grid area based on the flow field variable value of the target flow field grid area, and determining the vector dot product value of the pressure gradient and the velocity vector; performing dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or two-dimensional flow field area of ​​the target flow field grid area to obtain processed data; determining whether there is a shock wave in the target flow field grid area based on the processed data. In this way, shock waves can be judged by the velocity and pressure distribution characteristics of the flow field, and shock waves in complex flow fields can be accurately identified.

Description

Technical Field

[0001] The present invention relates to the field of fluid mechanics technology, and in particular to a flow field shock wave detection method, device, equipment and storage medium. Background Art

[0002] In supersonic compressible airflow, there is a flow discontinuity phenomenon called shock wave, which is characterized by a sudden jump in parameters when the airflow passes over it, and there is also a loss of mechanical energy, which is an irreversible process. Figure 1 The figure shows a shock wave appearing in the middle of the upper surface of an airfoil in transonic flow. A shock wave is a strong compression wave, a typical example of nonlinear waves in fluid mechanics, and a strong discontinuous wave. Its abrupt and energy dissipation characteristics are crucial to the design of high-speed engineering projects, such as aircraft.

[0003] In CFD (Computational Fluid Dynamics) numerical simulations and other applications involving supersonic aerodynamics (such as aircraft design and engine flow field research), it is often necessary to determine the location of shock waves in the flow field, that is, shock wave detection. Therefore, how to detect shock waves is a problem that needs to be solved. Summary of the Invention

[0004] In view of this, the purpose of the present invention is to provide a flow field shock wave detection method, device, equipment and storage medium that can accurately identify shock waves in complex flow fields. The specific solution is as follows:

[0005] In a first aspect, the present application discloses a flow field shock wave detection method, comprising:

[0006] Using a preset simulation tool to simulate the flow field distribution of the target device in the target operating environment to obtain the flow field variable value of each flow field grid area corresponding to the target device;

[0007] Determining a pressure gradient of a target flow field grid area based on a flow field variable value of the target flow field grid area, and determining a vector dot product value of the pressure gradient and a velocity vector;

[0008] Performing dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or the two-dimensional flow field area of ​​the target flow field grid region to obtain processed data;

[0009] Determine whether a shock wave exists in the target flow field grid area based on the processed data.

[0010] Optionally, simulating the flow field distribution of the target device under the target operating environment using a preset simulation tool to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device includes:

[0011] Obtain the boundary conditions and incoming flow conditions of the target operating environment and target device, and input the boundary conditions and the incoming flow conditions into a preset simulation tool to simulate the flow field distribution to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device.

[0012] Optionally, the step of obtaining boundary conditions and incoming flow conditions of the target operating environment and target device, and inputting the boundary conditions and incoming flow conditions into a preset simulation tool to simulate the flow field distribution to obtain flow field variable values ​​of each flow field grid area corresponding to the target device, includes:

[0013] Obtain the boundary conditions and incoming flow conditions of the target operating environment and target device, and process the boundary conditions and the incoming flow conditions using computational fluid dynamics methods or wind tunnel test methods to obtain the velocity vector and flow field pressure of each flow field grid area corresponding to the target device.

[0014] Optionally, determining the pressure gradient of the target flow field grid area based on the flow field variable value of the target flow field grid area, and determining the vector dot product value of the pressure gradient and the velocity vector, includes:

[0015] Determining a first flow field pressure of a target flow field grid area and a second flow field pressure of other flow field grid areas adjacent to the target flow field grid area from the flow field variable values ​​of each flow field grid area;

[0016] Processing the first flow field pressure and the second flow field pressure based on a least squares method or a Gaussian method to obtain a pressure gradient in the target flow field grid area;

[0017] The velocity vector of the target flow field grid area is determined from the flow field variable values ​​of each flow field grid area, and the vector dot product value of the pressure gradient and the velocity vector is determined.

[0018] Optionally, determining whether a shock wave exists in the target flow field grid area based on the processed data includes:

[0019] Processing the processed data based on a preset switching function model to obtain a corresponding output value, and determining whether the output value satisfies a preset condition; the preset condition being a condition that the preset output value is approximately equal to zero;

[0020] If the output value satisfies the preset condition, it is determined that a shock wave exists in the target flow field grid area.

[0021] Optionally, processing the processed data based on a preset switching function model to obtain a corresponding output value includes:

[0022] Obtaining processed data of other flow field grid areas related to the target flow field grid area, and inputting the processed data of the target flow field grid area and the processed data of the other flow field grid areas into a preset maximum function model to obtain target analysis data;

[0023] The target analysis data is processed based on a preset switching function model to obtain a corresponding output value.

[0024] Optionally, determining whether a shock wave exists in the target flow field grid area based on the processed data includes:

[0025] Determine a target preset threshold based on a preset numerical judgment condition, and determine whether the processed data is greater than the target preset threshold;

[0026] If the processed data is greater than the target preset threshold, it is determined that a shock wave exists in the target flow field grid area.

[0027] In a second aspect, the present application discloses a flow field shock wave detection device, comprising:

[0028] A flow field variable value determination module is used to simulate the flow field distribution of the target device under the target operating environment using a preset simulation tool to obtain the flow field variable value of each flow field grid area corresponding to the target device;

[0029] a variable value processing module, configured to determine a pressure gradient of a target flow field grid region based on a flow field variable value of the target flow field grid region, and determine a vector dot product value of the pressure gradient and a velocity vector;

[0030] a data processing module, configured to perform dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or the two-dimensional flow field area of ​​the target flow field grid region to obtain processed data;

[0031] A shock wave detection module is used to determine whether a shock wave exists in the target flow field grid area based on the processed data.

[0032] In a third aspect, the present application discloses an electronic device, comprising:

[0033] Memory, used to store computer programs;

[0034] A processor is used to execute the computer program to implement the aforementioned flow field shock wave detection method.

[0035] In a fourth aspect, the present application discloses a computer-readable storage medium for storing a computer program, which implements the aforementioned flow field shock wave detection method when executed by a processor.

[0036] It can be seen that in this application, a preset simulation tool is used to simulate the flow field distribution of the target device under the target operating environment to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device; the pressure gradient of the target flow field grid area is determined based on the flow field variable value of the target flow field grid area, and the vector dot product value of the pressure gradient and the velocity vector is determined; the vector dot product value is dimensionlessly processed based on the three-dimensional flow field volume or two-dimensional flow field area of ​​the target flow field grid area to obtain processed data; based on the processed data, it is determined whether there is a shock wave in the target flow field grid area. In this way, by using the velocity and pressure distribution characteristics of the flow field to judge shock waves, shock waves in complex flow fields can be accurately identified and will not be misjudged by expansion waves. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0038] Figure 1 A schematic diagram of shock waves on an upper wing surface of an airfoil disclosed in this application;

[0039] Figure 2 This is a flow chart of a flow field shock wave detection method disclosed in this application;

[0040] Figure 3 A schematic diagram of a flow field grid disclosed in this application;

[0041] Figure 4 This is a schematic structural diagram of a flow field shock wave detection device disclosed in this application;

[0042] Figure 5 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0043] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0044] In supersonic aerodynamics, determining the location of shock waves in the flow field is often necessary, a process known as shock wave detection. For example, in aircraft design and optimization, the distribution of shock waves on the vehicle surface is detected to optimize airfoils or air inlets to reduce drag and heat load. During re-entry, shock wave detection is used to assess the intensity of aerodynamic heating and design high-temperature-resistant cooling systems. In medical applications, shock lithotripsy is commonly used to break up kidney or gallstones. This process requires precise control of shock wave intensity and positioning, which is also essential for shock wave detection. In CFD, shock wave detection is used to verify the accuracy of shock wave fitting algorithms, such as HLLC (HLL with Contact) and HLLE++ (a format used for numerical simulation). The grid resolution can also be dynamically adjusted based on the detected shock wave location to improve the simulation accuracy of complex flows, such as shock wave / boundary layer interaction. Therefore, this application specifically describes a flow field shock wave detection method that can identify shock waves in complex flow fields.

[0045] See also Figure 2 As shown, the embodiment of the present application discloses a flow field shock wave detection method, comprising:

[0046] Step S11: using a preset simulation tool to simulate the flow field distribution of the target device in the target operating environment to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device.

[0047] In this embodiment, the method of simulating the flow field distribution of the target device under the target operating environment using a preset simulation tool to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device includes: obtaining the boundary conditions and incoming flow conditions of the target operating environment and the target device, and inputting the boundary conditions and the incoming flow conditions into the preset simulation tool to simulate the flow field distribution to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device. In this embodiment, the method of obtaining the boundary conditions and incoming flow conditions of the target operating environment and the target device, and inputting the boundary conditions and the incoming flow conditions into the preset simulation tool to simulate the flow field distribution to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device includes: obtaining the boundary conditions and incoming flow conditions of the target operating environment and the target device, and processing the boundary conditions and the incoming flow conditions using a computational fluid dynamics method or a wind tunnel test method to obtain the velocity vector and flow field pressure of each flow field grid area corresponding to the target device. Specifically, from actual application scenarios, such as aircraft flow field, pipeline flow field, engine flow field and other scenarios, certain boundary conditions and incoming flow conditions are determined, wherein the boundary conditions include but are not limited to the aircraft's shape parameters, the temperature of the surface, and the flow conditions of the nozzle nozzle. The incoming flow conditions include but are not limited to the Reynolds number, the incoming flow velocity, the aircraft's flight altitude, the aircraft's attitude angle, etc. Then, the CFD method or the wind tunnel test method is used to simulate the determined boundary conditions and incoming flow conditions, and some flow variables in the flow field (including density , velocity vector The distribution of flow field grid area obtained by simulation is as follows: Figure 3 Taking CFD as an example, after each numerical iteration step of CFD is completed, the variable value of the flow field unit center is obtained, including density , velocity vector and pressure p.

[0048] Step S12: determining the pressure gradient of the target flow field grid area based on the flow field variable value of the target flow field grid area, and determining the vector dot product value of the pressure gradient and the velocity vector.

[0049] In this embodiment, the method of determining the pressure gradient of the target flow field grid area based on the flow field variable value of the target flow field grid area, and determining the vector dot product value of the pressure gradient and the velocity vector, includes: determining the first flow field pressure of the target flow field grid area and the second flow field pressure of other flow field grid areas adjacent to the target flow field grid area from the flow field variable values ​​of each flow field grid area; processing the first flow field pressure and the second flow field pressure based on the least squares method or the Gaussian method to obtain the pressure gradient of the target flow field grid area; determining the velocity vector of the target flow field grid area from the flow field variable values ​​of each flow field grid area, and determining the vector dot product value of the pressure gradient and the velocity vector. Figure 3 In the example, flow field grid area C2, flow field grid area C3 and flow field grid area C4 are adjacent flow field grid areas of flow field grid area C1. Using the least square method or Gaussian method, the gradient of the target flow field grid area can be calculated. , here we take Gaussian method as an example:

[0050] ;

[0051] in is the volume of the flow field grid area C1, is the number of faces in the flow field grid area C1, and also the number of adjacent grids of C1. is the unit surface normal vector of the i-th surface, pointing to the outside of the unit cell in the flow field grid area, is the area of ​​the ith face. Then, calculate the velocity vector With pressure gradient The dot product of:

[0052] ;

[0053] It should be noted here that in the shock wave region, The value of is positive; and for expansion waves, The value of is negative. For smooth areas, The value of is small; while for the area where there are compression waves or expansion waves, The value of is larger.

[0054] Step S13: performing dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or the two-dimensional flow field area of ​​the target flow field grid region to obtain processed data.

[0055] In this embodiment, the vector dot product value is dimensionless processed using the following formula:

[0056] ;

[0057] Where c is the speed of sound, It is the cube root of the volume. When the flow field is two-dimensional, it is the square root of the area. It is a dimension of length.

[0058] Step S14: determining whether shock waves exist in the target flow field grid area based on the processed data.

[0059] In this embodiment, the determination of whether shock waves exist in the target flow field grid area based on the processed data includes: determining a target preset threshold based on a preset numerical judgment condition, and judging whether the processed data is greater than the target preset threshold; if the processed data is greater than the target preset threshold, then determining that shock waves exist in the target flow field grid area. That is, after the dimensionless calculation is obtained, a threshold parameter can be given. ,when When , there is a shock wave here. According to numerical simulation experience, The possible value is 75.

[0060] However, considering the actual situation, there may be a situation where the target flow field grid area is a non-structured network or the shock wave is not aligned with the grid. Therefore, in this embodiment, the determination of whether there is a shock wave in the target flow field grid area based on the processed data includes: processing the processed data based on a preset switching function model to obtain a corresponding output value, and judging whether the output value meets a preset condition; the preset condition is a condition that the preset output value is approximately equal to zero; if the output value meets the preset condition, it is determined that there is a shock wave in the target flow field grid area. Among them, the processing of the processed data based on the preset switching function model to obtain a corresponding output value includes: obtaining the processed data of other flow field grid areas with the target flow field grid area, and inputting the processed data of the target flow field grid area and the processed data of the other flow field grid areas into a preset maximum function model to obtain target analysis data; processing the target analysis data based on the preset switching function model to obtain a corresponding output value. Specifically, the processed data of the target flow field grid area is processed in the manner of the maximum value of adjacent units. That is, based on the following formula:

[0061] ;

[0062] Where j is the adjacent flow field grid area of ​​the flow field grid area C1. Figure 3 These are cells C2, C3, and C4.

[0063] Then define a preset switch function model :

[0064] ;

[0065] in, The range of is [0,1]. Based on this, the above calculation can be performed on all flow field grid areas under the entire flow field. The position where the value is equal to 0 is the position of the shock wave.

[0066] As can be seen, in this embodiment, a preset simulation tool is used to simulate the flow field distribution of the target device under the target operating environment to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device; the pressure gradient of the target flow field grid area is determined based on the flow field variable values ​​of the target flow field grid area, and the vector dot product value of the pressure gradient and the velocity vector is determined; the vector dot product value is dimensionlessly processed based on the three-dimensional flow field volume or two-dimensional flow field area of ​​the target flow field grid area to obtain processed data; and based on the processed data, it is determined whether shock waves exist in the target flow field grid area. In this way, by using the velocity and pressure distribution characteristics of the flow field to judge shock waves, shock waves in complex flow fields can be accurately identified and will not be misjudged by expansion waves.

[0067] refer to Figure 4 The present application also discloses a flow field shock wave detection device, including:

[0068] The flow field variable value determination module 11 is used to simulate the flow field distribution of the target device under the target operating environment using a preset simulation tool to obtain the flow field variable value of each flow field grid area corresponding to the target device;

[0069] a variable value processing module 12, configured to determine a pressure gradient of a target flow field grid region based on a flow field variable value of the target flow field grid region, and determine a vector dot product value of the pressure gradient and a velocity vector;

[0070] A data processing module 13 is configured to perform dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or the two-dimensional flow field area of ​​the target flow field grid region to obtain processed data;

[0071] The shock wave detection module 14 is configured to determine whether a shock wave exists in the target flow field grid area based on the processed data.

[0072] It can be seen that in this embodiment, by utilizing the velocity and pressure distribution characteristics of the flow field to perform shock wave judgment, shock waves in complex flow fields can be accurately identified without being misjudged by expansion waves.

[0073] In some specific embodiments, the flow field variable value determination module 11 may specifically include:

[0074] The flow field simulation submodule is used to obtain the boundary conditions and incoming flow conditions of the target operating environment and target device, and input the boundary conditions and the incoming flow conditions into a preset simulation tool to simulate the flow field distribution to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device.

[0075] In some specific embodiments, the flow field simulation submodule may specifically include:

[0076] The flow field variable processing unit is used to obtain the boundary conditions and incoming flow conditions of the target operating environment and the target device, and process the boundary conditions and the incoming flow conditions using computational fluid dynamics methods or wind tunnel test methods to obtain the velocity vector and flow field pressure of each flow field grid area corresponding to the target device.

[0077] In some specific embodiments, the variable value processing module 12 may specifically include:

[0078] a pressure determining unit, configured to determine a first flow field pressure of a target flow field grid area and a second flow field pressure of other flow field grid areas adjacent to the target flow field grid area from the flow field variable values ​​of the flow field grid areas;

[0079] a pressure gradient acquiring unit, configured to process the first flow field pressure and the second flow field pressure based on a least squares method or a Gaussian method to obtain a pressure gradient in the target flow field grid area;

[0080] The vector dot product value determining unit is used to determine the velocity vector of the target flow field grid area from the flow field variable values ​​of each flow field grid area, and determine the vector dot product value of the pressure gradient and the velocity vector.

[0081] In some specific embodiments, the shock wave detection module 14 may specifically include:

[0082] An output value acquisition submodule is used to process the processed data based on a preset switching function model to obtain a corresponding output value, and to determine whether the output value meets a preset condition; the preset condition is that the preset output value is approximately equal to zero;

[0083] The first shock wave determination unit is configured to determine whether a shock wave exists in the target flow field grid area if the output value satisfies the preset condition.

[0084] In some specific embodiments, the output value acquisition submodule may specifically include:

[0085] a data processing unit, configured to obtain processed data of other flow field grid areas than the target flow field grid area, and input the processed data of the target flow field grid area and the processed data of the other flow field grid areas into a preset maximum function model to obtain target analysis data;

[0086] The data analysis unit is used to process the target analysis data based on a preset switching function model to obtain a corresponding output value.

[0087] In some specific embodiments, the shock wave detection module 14 may specifically include:

[0088] A threshold judgment unit, configured to determine a target preset threshold based on a preset numerical judgment condition, and to judge whether the processed data is greater than the target preset threshold;

[0089] The second shock wave determination unit is configured to determine that a shock wave exists in the target flow field grid area if the processed data is greater than the target preset threshold.

[0090] Furthermore, the embodiment of the present application also discloses an electronic device, Figure 5 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram should not be considered as any limitation to the scope of application of the present application.

[0091] Figure 5 This is a schematic diagram of the structure of an electronic device 20 provided in an embodiment of the present application. The electronic device 20 may specifically 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 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the flow field shock wave detection method disclosed in any of the aforementioned embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0092] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device. The communication protocol it follows is any communication protocol that can be applied to the technical solution of this application and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world. Its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0093] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or CD, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0094] The operating system 221 is used to manage and control the hardware devices and computer program 222 on the electronic device 20, and can be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program capable of implementing the flow field shock wave detection method performed by the electronic device 20 disclosed in any of the aforementioned embodiments, the computer program 222 can further include a computer program capable of implementing other specific tasks.

[0095] Furthermore, this application discloses a computer-readable storage medium for storing a computer program; wherein, when executed by a processor, the computer program implements the aforementioned flow field shock wave detection method. The specific steps of this method can be referred to the corresponding contents disclosed in the aforementioned embodiments and will not be repeated here.

[0096] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. Reference can be made to the descriptions of the identical or similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple, and the relevant parts can be referred to the descriptions of the methods.

[0097] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to their functions. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0098] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in random access memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

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

[0100] The above is a detailed introduction to the technical solution provided by the present application. Specific examples are used herein to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea. At the same time, for those skilled in the art, according to the ideas of the present application, there may be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A flow field shock wave detection method, characterized in that: include: Using a preset simulation tool to simulate the flow field distribution of the target device in the target operating environment to obtain the flow field variable value of each flow field grid area corresponding to the target device; Determining a pressure gradient of a target flow field grid area based on a flow field variable value of the target flow field grid area, and determining a vector dot product value of the pressure gradient and a velocity vector; Performing dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or the two-dimensional flow field area of ​​the target flow field grid region to obtain processed data; The vector dot product value is dimensionless processed using the following formula: ； Where c is the speed of sound, It is the cube root of volume. When the flow field is two-dimensional, it is the square root of area. It is a dimension of length. determining whether a shock wave exists in the target flow field grid area based on the processed data; Wherein, determining whether a shock wave exists in the target flow field grid area based on the processed data includes: processing the processed data based on a preset switching function model to obtain a corresponding output value, and judging whether the output value satisfies a preset condition; the preset condition is a condition that the preset output value is approximately equal to zero; if the output value satisfies the preset condition, then judging that a shock wave exists in the target flow field grid area; Among them, the preset switch function model : ； in, The range of is [0,1]. There is a shock wave.

2. The flow field shock wave detection method according to claim 1, characterized in that: The method of simulating the flow field distribution of the target device under the target operating environment by using a preset simulation tool to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device includes: Obtain the boundary conditions and incoming flow conditions of the target operating environment and target device, and input the boundary conditions and the incoming flow conditions into a preset simulation tool to simulate the flow field distribution to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device.

3. The flow field shock wave detection method according to claim 2, characterized in that: The step of obtaining the boundary conditions and incoming flow conditions of the target operating environment and the target device, and inputting the boundary conditions and the incoming flow conditions into a preset simulation tool to simulate the flow field distribution to obtain the flow field variable values ​​of each flow field grid area corresponding to the target device includes: Obtain the boundary conditions and incoming flow conditions of the target operating environment and target device, and process the boundary conditions and the incoming flow conditions using computational fluid dynamics methods or wind tunnel test methods to obtain the velocity vector and flow field pressure of each flow field grid area corresponding to the target device.

4. The flow field shock wave detection method according to claim 1, characterized in that: The determining of the pressure gradient of the target flow field grid area based on the flow field variable value of the target flow field grid area, and determining the vector dot product value of the pressure gradient and the velocity vector, includes: Determining a first flow field pressure of a target flow field grid area and a second flow field pressure of other flow field grid areas adjacent to the target flow field grid area from the flow field variable values ​​of each flow field grid area; Processing the first flow field pressure and the second flow field pressure based on a least squares method or a Gaussian method to obtain a pressure gradient in the target flow field grid area; The velocity vector of the target flow field grid area is determined from the flow field variable values ​​of each flow field grid area, and the vector dot product value of the pressure gradient and the velocity vector is determined.

5. The flow field shock wave detection method according to claim 1, characterized in that: The processing of the processed data based on a preset switching function model to obtain a corresponding output value includes: Obtaining processed data of other flow field grid areas related to the target flow field grid area, and inputting the processed data of the target flow field grid area and the processed data of the other flow field grid areas into a preset maximum function model to obtain target analysis data; The target analysis data is processed based on a preset switching function model to obtain a corresponding output value.

6. The flow field shock wave detection method according to any one of claims 1 to 4, characterized in that: Determining whether a shock wave exists in the target flow field grid area based on the processed data includes: Determine a target preset threshold based on a preset numerical judgment condition, and determine whether the processed data is greater than the target preset threshold; If the processed data is greater than the target preset threshold, it is determined that a shock wave exists in the target flow field grid area.

7. A flow field shock wave detection device, characterized in that: include: A flow field variable value determination module is used to simulate the flow field distribution of the target device under the target operating environment using a preset simulation tool to obtain the flow field variable value of each flow field grid area corresponding to the target device; a variable value processing module, configured to determine a pressure gradient of a target flow field grid region based on a flow field variable value of the target flow field grid region, and determine a vector dot product value of the pressure gradient and a velocity vector; a data processing module, configured to perform dimensionless processing on the vector dot product value based on the three-dimensional flow field volume or the two-dimensional flow field area of ​​the target flow field grid region to obtain processed data; The vector dot product value is dimensionless processed using the following formula: ； Where c is the speed of sound, It is the cube root of volume. When the flow field is two-dimensional, it is the square root of area. It is a dimension of length. a shock wave detection module, configured to determine whether a shock wave exists in the target flow field grid area based on the processed data; Wherein, determining whether a shock wave exists in the target flow field grid area based on the processed data includes: processing the processed data based on a preset switching function model to obtain a corresponding output value, and judging whether the output value satisfies a preset condition; the preset condition is a condition that the preset output value is approximately equal to zero; if the output value satisfies the preset condition, then judging that a shock wave exists in the target flow field grid area; Among them, the preset switch function model : ； in, The range of is [0,1]. There is a shock wave.

8. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor, configured to execute the computer program to implement the flow field shock wave detection method according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that Used to store a computer program, which, when executed by a processor, implements the flow field shock wave detection method according to any one of claims 1 to 6.

Images