Combustion oscillation characteristic determination method and device

Abstract

The invention discloses a combustion oscillation characteristic determination method and device, and the method comprises the steps: constructing flow field data of a combustor for simulating a flow field of the combustor in a combustion state, and collecting the pressure evolution data of each sampling point and the pressure snapshot data of a whole flow field; local combustion oscillation characteristics corresponding to the sampling points and global combustion oscillation characteristics corresponding to the whole flow field are obtained through Fourier transform and a dynamic mode decomposition mode, matching is carried out based on frequency data of the two characteristics, and according to the matching result, the local combustion oscillation characteristics corresponding to the sampling points and the global combustion oscillation characteristics corresponding to the whole flow field are obtained. And associating the local combustion oscillation characteristics with the global combustion oscillation characteristics, which are successfully matched in frequency, so as to obtain a characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics. Through the scheme provided by the invention, more perfect combustion oscillation characteristic data can be obtained to be applied to combustion oscillation analysis, so that the accuracy of combustion oscillation analysis work can be further improved.

Description

Technical Field

[0001] The present application relates to the field of power and energy technology, and in particular to a method and device for determining combustion oscillation characteristics. Background Art

[0002] In many applications in the field of energy and power, the combustion of fossil fuels is still the main source of power and energy for mankind. Combustion oscillation is an unstable combustion phenomenon in which a large pressure pulsation occurs in the combustion chamber or afterburner during the operation of the burner. It is more likely to occur when the flow velocity of the flow field in the combustion chamber of the burner is low, and it is relatively rare in high-flow rate burners where the flow velocity exceeds the speed of sound. However, according to the actual application, combustion oscillation may occur even when the flow velocity in the burner exceeds the speed of sound. Due to the complexity of combustion oscillation in high-speed flow, technicians lack understanding of combustion oscillation in high-flow environment. Therefore, in order to ensure the stable operation of the burner in a high-flow environment, it is a technical problem that technicians in this field must overcome at this stage to study combustion oscillation in high-speed flow and master the characteristics of combustion oscillation in high-speed flow.

[0003] At present, the characteristic analysis of combustion oscillation is mainly carried out by Fourier transform or dynamic analysis to determine the characteristic data for analysis and research. However, the characteristic data determined by the existing combustion oscillation analysis methods have obvious limitations, which restricts the further improvement of the accuracy of combustion oscillation analysis. Summary of the invention

[0004] The present application provides a method and device for determining combustion oscillation characteristics, which are used to solve the technical problem of low working accuracy of existing combustion oscillation analysis methods.

[0005] In order to solve the above technical problems, the first aspect of the present application provides a method for determining combustion oscillation characteristics, comprising:

[0006] According to the equipment structure information of the burner, the flow field data inside the burner in the unburned state is constructed by fluid mechanics simulation;

[0007] Based on the flow field data and in combination with preset fuel injection configuration information, the flow field of the burner in the combustion state is simulated, and pressure evolution data of each sampling point and pressure snapshot data of the entire flow field are collected;

[0008] Performing Fourier transformation on the pressure evolution data to obtain first frequency data of the pressure at each sampling point and a pressure evolution process image corresponding to the first frequency data, so as to obtain a local combustion oscillation feature corresponding to the sampling point;

[0009] Performing dynamic mode decomposition on the pressure snapshot data to obtain second frequency data of the full flow field pressure and a pressure dynamic mode image corresponding to the second frequency data, so as to obtain global combustion oscillation characteristics corresponding to the full flow field;

[0010] Matching is performed based on the first frequency data and the second frequency data, and based on the matching result, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain a characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0011] Preferably, the first frequency data and the second frequency data both include: a peak frequency and a dominant frequency, wherein the dominant frequency is a peak frequency whose amplitude corresponding to the frequency exceeds a preset amplitude threshold.

[0012] Preferably, determining the local combustion oscillation characteristics corresponding to the sampling point according to the frequency data specifically includes:

[0013] Determine, according to the amplitude corresponding to the peak frequency of the sampling point, a target peak frequency corresponding to the maximum amplitude, and determine a dominant oscillation period corresponding to the target peak frequency;

[0014] Extract the pressure evolution image information within the dominant oscillation period.

[0015] Preferably, matching the first frequency data with the second frequency data and associating the local combustion oscillation characteristics with successful frequency matching with the global combustion oscillation characteristics according to the matching result specifically includes:

[0016] According to a preset dominant frequency screening condition, performing data processing on the first frequency data to determine a dominant frequency in the first frequency data;

[0017] According to matching the dominant frequency in the first frequency data with the peak frequency of the second frequency data, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0018] Preferably, according to the equipment structure information of the burner, constructing the flow field data of the burner in the unburned state by means of fluid mechanics simulation specifically includes:

[0019] According to the equipment structure information of the burner, the steady flow field data of the inside of the burner in the unburned state is constructed by means of steady numerical simulation of fluid mechanics;

[0020] Using the steady flow field data as the initial field model, the flow field simulation data of a preset time length is calculated by an unsteady numerical simulation method of fluid mechanics;

[0021] Unsteady flow field data are constructed based on the flow field simulation data, and the unsteady flow field data are used as flow field data inside the burner in an unburned state.

[0022] At the same time, the second aspect of the present application provides a combustion oscillation characteristic determination device, comprising:

[0023] A flow field data construction unit, used to construct flow field data of the burner in an unburned state by fluid mechanics simulation according to the equipment structure information of the burner;

[0024] A combustion simulation unit, for simulating the flow field of the burner in a combustion state based on the flow field data and in combination with preset fuel injection configuration information, and collecting pressure evolution data of each sampling point and pressure snapshot data of the entire flow field;

[0025] A local feature acquisition unit, configured to perform Fourier transform on the pressure evolution data to obtain first frequency data of the pressure at each sampling point and a pressure evolution process image corresponding to the first frequency data, so as to obtain a local combustion oscillation feature corresponding to the sampling point;

[0026] A global feature acquisition unit, used for performing dynamic mode decomposition on the pressure snapshot data to obtain second frequency data of the full flow field pressure and a pressure dynamic mode image corresponding to the second frequency data, so as to obtain a global combustion oscillation feature corresponding to the full flow field;

[0027] A feature association unit is used to match the first frequency data with the second frequency data, and according to the matching result, associate the local combustion oscillation feature with the global combustion oscillation feature that successfully matches the frequency to obtain the feature relationship between the local combustion oscillation feature and the global combustion oscillation feature.

[0028] Preferably, the first frequency data and the second frequency data both include: a peak frequency and a dominant frequency, wherein the dominant frequency is a peak frequency whose amplitude corresponding to the frequency exceeds a preset amplitude threshold.

[0029] Preferably, the local feature acquisition unit is specifically used for:

[0030] Performing Fourier transformation on the pressure evolution data to obtain the peak frequency of each sampling point;

[0031] Determine, according to the amplitude corresponding to the peak frequency of the sampling point, a target peak frequency corresponding to the maximum amplitude, and determine a dominant oscillation period corresponding to the target peak frequency;

[0032] Extract the pressure evolution image information within the dominant oscillation period.

[0033] Preferably, the feature association unit is specifically used for:

[0034] According to a preset dominant frequency screening condition, performing data processing on the first frequency data to determine a dominant frequency in the first frequency data;

[0035] According to matching the dominant frequency in the first frequency data with the peak frequency of the second frequency data, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0036] Preferably, the flow field data construction unit is specifically used for:

[0037] According to the equipment structure information of the burner, the steady flow field data of the inside of the burner in the unburned state is constructed by means of steady numerical simulation of fluid mechanics;

[0038] Using the steady flow field data as the initial field model, the flow field simulation data of a preset time length is calculated by an unsteady numerical simulation method of fluid mechanics;

[0039] Unsteady flow field data are constructed based on the flow field simulation data, and the unsteady flow field data are used as flow field data inside the burner in an unburned state.

[0040] It can be seen from the above technical solutions that this application has the following advantages:

[0041] This application constructs the flow field data of the burner to simulate the flow field of the burner under the combustion state, and collects the pressure evolution data of each sampling point and the pressure snapshot data of the whole flow field at the same time, and obtains the local combustion oscillation characteristics corresponding to the sampling point and the global combustion oscillation characteristics corresponding to the whole flow field through Fourier transform and dynamic mode decomposition respectively, and matches the frequency data of the two characteristics. According to the matching result, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics, so as to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics. The scheme provided by this application can obtain more complete combustion oscillation characteristic data for application to combustion oscillation analysis, so that the accuracy of combustion oscillation analysis can be further improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0043] Figure 1 A schematic flow chart of a basic embodiment of a method for determining combustion oscillation characteristics provided in the present application.

[0044] Figure 2 A schematic flow chart of a specific embodiment of a method for determining combustion oscillation characteristics provided in the present application.

[0045] Figure 3 A schematic structural diagram of an embodiment of a combustion oscillation characteristic determination device provided in the present application. DETAILED DESCRIPTION

[0046] The embodiments of the present application provide a method and device for determining combustion oscillation characteristics, which are used to solve the technical problem of low working accuracy of existing combustion oscillation analysis methods.

[0047] In order to make the purpose, features, and advantages of the invention of this application more obvious and easy to understand, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the embodiments described below are only part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0048] First, a detailed description of an embodiment of a method for determining combustion oscillation characteristics provided by the present application is as follows:

[0049] See also Figure 1 , this embodiment provides a method for determining combustion oscillation characteristics, including:

[0050] Step 101: According to the equipment structure information of the burner, the flow field data inside the burner in the unburned state is constructed by fluid mechanics simulation.

[0051] Step 102: Based on the flow field data and in combination with the preset fuel injection configuration information, the flow field of the burner in the combustion state is simulated, and the pressure evolution data of each sampling point and the pressure snapshot data of the entire flow field are collected.

[0052] It should be noted that, first, based on the equipment structure information of the burner, the flow field data of the burner when no fuel is injected, that is, when no combustion occurs, is constructed. Then, based on this flow field data and in combination with the preset fuel injection configuration information, the flow field of the burner in the combustion state is simulated. At the same time, the pressure evolution data of each sampling point in the simulation process and the pressure snapshot data of the entire flow field are collected for use in subsequent steps. The sampling points mentioned in this embodiment are preferably arranged on the axis of the combustion chamber within the range from the fuel injection port to the combustion chamber outlet, with a number of 6 to 8 points. The pressure evolution data of the sampling points refers to the point data of the pressure taken once in each non-steady numerical simulation time step, and the pressure snapshot data of the flow field refers to the field data of the pressure taken once at regular intervals.

[0053] It should be noted that the collected data include: pressure evolution data of sampling points for 8 incoming flow cycles and full flow field pressure snapshot data for 2 incoming flow cycles. The incoming flow cycle refers to the time required for the fluid to flow from the inlet to the outlet at the incoming flow velocity. More specifically, the above-mentioned 2 incoming flow cycles are preferably the last 2 cycles among the 8 incoming flow cycles.

[0054] Step 103: Perform Fourier transform on the pressure evolution data to obtain first frequency data of the pressure at each sampling point and a pressure evolution process image corresponding to the first frequency data, so as to obtain a local combustion oscillation feature corresponding to the sampling point.

[0055] Step 104 , performing dynamic modal decomposition on the pressure snapshot data to obtain second frequency data of the full flow field pressure and a pressure dynamic modal image corresponding to the second frequency data, so as to obtain a global combustion oscillation feature corresponding to the full flow field.

[0056] It should be noted that based on the pressure evolution data of the sampling points obtained in step 102 and the pressure snapshot data of the entire flow field, the pressure evolution data is processed by Fourier transform to obtain the local combustion oscillation characteristics corresponding to each sampling point; the pressure snapshot data is processed by dynamic mode decomposition to obtain the global combustion oscillation characteristics corresponding to the entire flow field.

[0057] Among them, the dynamic modal image of the frequency of the full flow field pressure refers to an image drawn by plotting the mode corresponding to the frequency among the multiple modes obtained by performing dynamic modal decomposition on the flow field pressure snapshot data. The implementation steps of "obtaining the second frequency data of the full flow field pressure and the pressure dynamic modal image corresponding to the second frequency data" are: first perform dynamic modal decomposition on the flow field pressure snapshot data, and find the abscissa value when the peak value appears on the ordinate in the obtained frequency-amplitude diagram, that is, obtain the second frequency data of the full flow field pressure; then plot the mode corresponding to the peak frequency into a cloud diagram, that is, obtain the pressure dynamic modal image corresponding to the second frequency data.

[0058] Step 105 , matching the first frequency data with the second frequency data, and according to the matching result, associating the local combustion oscillation characteristics with successful frequency matching with the global combustion oscillation characteristics to obtain a characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0059] In summary, this embodiment constructs the flow field data of the burner to simulate the flow field of the burner under the combustion state, and collects the pressure evolution data of each sampling point and the pressure snapshot data of the whole flow field at the same time, and obtains the local combustion oscillation characteristics corresponding to the sampling point and the global combustion oscillation characteristics corresponding to the whole flow field through Fourier transform and dynamic mode decomposition respectively, and matches the frequency data of the two characteristics. According to the matching result, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics, so as to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics. The scheme provided by this application can obtain more complete combustion oscillation characteristic data for application to combustion oscillation analysis, so that the accuracy of combustion oscillation analysis can be further improved.

[0060] The above is a detailed description of a basic embodiment of a method for determining a combustion oscillation characteristic provided by the present application. The following is a detailed description of a specific embodiment of a method for determining a combustion oscillation characteristic provided by the present application based on the above basic embodiment, which is as follows:

[0061] Furthermore, if Figure 2 As shown, the step 101 mentioned in the above embodiment specifically includes:

[0062] Step 1011: According to the equipment structure information of the burner, the steady flow field data of the burner in the unburned state is constructed by means of steady numerical simulation of fluid mechanics.

[0063] Step 1012: Using the steady flow field data as the initial field model, the flow field simulation data of a preset time length is calculated by an unsteady numerical simulation method of fluid mechanics.

[0064] Step 1013: construct unsteady flow field data based on the flow field simulation data, and use the unsteady flow field data as the flow field data inside the burner in an unburned state.

[0065] It should be noted that the present embodiment provides a method for obtaining a steady flow field without fuel injection by using a conventional computational fluid dynamics steady numerical simulation method, and then using this steady flow field as the initial field, and continuing to calculate the flow field simulation data of four incoming flow cycles by using a conventional computational fluid dynamics unsteady numerical simulation method, and constructing unsteady flow field data based on the flow field simulation data, thereby obtaining an unsteady flow field without fuel injection, and using this unsteady flow field data as the flow field data used in subsequent simulation steps, wherein the fluid mechanics steady numerical simulation method mentioned in the present embodiment can be The steady solution method of the Reynolds-averaged Navier-Stokes equations and the turbulence model equations for the flow is performed using computational fluid dynamics software or program. The unsteady numerical simulation method of fluid mechanics can be an unsteady solution method of the Reynolds-averaged Navier-Stokes equations and the turbulence model equations for the flow using computational fluid dynamics software or program. Similarly, the unsteady simulation method of fluid mechanics under combustion conditions can be an unsteady solution method of the Reynolds-averaged Navier-Stokes equations and the component mass transport equations, the turbulence model equations and the turbulent combustion model equations using computational fluid dynamics software or program.

[0066] Furthermore, both the first frequency data and the second frequency data include: a peak frequency and a dominant frequency, wherein the dominant frequency is a peak frequency whose amplitude corresponding to the frequency exceeds a preset amplitude threshold, wherein the peak frequency of this embodiment specifically refers to the horizontal coordinate value when a peak appears on the vertical coordinate in the frequency (horizontal coordinate)-amplitude (vertical coordinate) diagram.

[0067] Furthermore, the step 103 mentioned in the above embodiment specifically includes:

[0068] According to the amplitude corresponding to the peak frequency of the sampling point, the target peak frequency corresponding to the maximum amplitude is determined, the dominant oscillation period corresponding to the target peak frequency is determined, and the pressure evolution image information within the dominant oscillation period is extracted as the local combustion oscillation feature corresponding to the sampling point.

[0069] It should be noted that the pressure evolution data of the sampling points are first subjected to a fast Fourier transform to obtain the peak frequency of each sampling point, and then the peak frequency with the largest amplitude is taken according to the amplitude corresponding to the peak frequency of each sampling point to determine its corresponding dominant oscillation period and obtain the pressure evolution process image near the sampling point within the duration of the dominant oscillation period. Among them, the reciprocal of the peak frequency is equal to the corresponding oscillation period. If the amplitude corresponding to a peak frequency is greater than a preset threshold, the peak frequency is the dominant frequency, and the reciprocal of the peak frequency is equal to the corresponding dominant oscillation period.

[0070] Furthermore, the step 105 mentioned in the above embodiment specifically includes:

[0071] Step 1051: Process the first frequency data according to a preset dominant frequency screening condition to determine the dominant frequency in the first frequency data.

[0072] Step 1052: Match the dominant frequency in the first frequency data with the peak frequency of the second frequency data, and associate the local combustion oscillation characteristics with the global combustion oscillation characteristics to obtain a characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0073] It should be noted that, based on the local combustion oscillation characteristics obtained in step 103, according to the frequency data corresponding to the combustion oscillation characteristics and the preset dominant frequency screening conditions, the sampling points with larger pressure amplitudes, that is, the sampling points with peak frequencies whose amplitudes exceed the preset threshold, are found, based on the dominant frequencies corresponding to these sampling points, the peak frequencies of the full flow field pressure obtained in step 104 are matched based on the dominant frequency information, and the dynamic modal images of these peak frequencies of the full flow field pressure are obtained, and the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0074] The relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics can be obtained through the above method steps. The obtained relationship characteristics can reflect two aspects of characteristics: (1) Whether the peak frequency with the largest amplitude in the local combustion oscillation characteristics is still the dominant frequency in the global combustion oscillation characteristics. If it is still the dominant frequency, it indicates that the local part is important in the global flow field. The larger the amplitude of the dominant frequency, the greater the importance of the local part corresponding to the dominant frequency; (2) For the dominant frequency in the local combustion oscillation characteristics of a certain part, whether there are other parts in the global flow field that also have combustion oscillations with the same frequency. If there are other parts that also have combustion oscillations with the same frequency, then the dominant frequency (important frequency) is important in multiple positions in the flow field. After obtaining this information, in the study of combustion oscillation, one is to determine which parts in the global flow field are important and how important they are; the other is to determine which parts of the global flow field are important for important frequencies.

[0075] The above content is a detailed description of a specific embodiment of a method for determining a combustion oscillation characteristic provided by the present application. The following is a detailed description of an embodiment of a device for determining a combustion oscillation characteristic provided by the present application, which is as follows:

[0076] See also Figure 3 , this embodiment provides a combustion oscillation characteristic determination device, comprising:

[0077] The flow field data construction unit 201 is used to construct the flow field data of the burner in the unburned state by means of fluid mechanics simulation according to the equipment structure information of the burner;

[0078] The combustion simulation unit 202 is used to simulate the flow field of the burner in the combustion state based on the flow field data and the preset fuel injection configuration information, and collect the pressure evolution data of each sampling point and the pressure snapshot data of the whole flow field;

[0079] The local feature acquisition unit 203 is used to perform Fourier transform on the pressure evolution data to obtain the first frequency data of the pressure at each sampling point and the pressure evolution process image corresponding to the first frequency data, so as to obtain the local combustion oscillation characteristics corresponding to the sampling point;

[0080] The global feature acquisition unit 204 is used to perform dynamic mode decomposition on the pressure snapshot data to obtain the second frequency data of the full flow field pressure and the pressure dynamic mode image corresponding to the second frequency data, so as to obtain the global combustion oscillation characteristics corresponding to the full flow field;

[0081] The feature association unit 205 is used to match the first frequency data with the second frequency data, and according to the matching result, associate the local combustion oscillation feature with the global combustion oscillation feature to obtain the feature relationship between the local combustion oscillation feature and the global combustion oscillation feature.

[0082] Furthermore, both the first frequency data and the second frequency data include: a peak frequency and a dominant frequency, wherein the dominant frequency is a peak frequency whose amplitude corresponding to the frequency exceeds a preset amplitude threshold.

[0083] Furthermore, the local feature acquisition unit is specifically used for:

[0084] Perform Fourier transform on the pressure evolution data to obtain the peak frequency of each sampling point;

[0085] According to the amplitude corresponding to the peak frequency of the sampling point, the target peak frequency corresponding to the maximum amplitude is determined, and the dominant oscillation period corresponding to the target peak frequency is determined;

[0086] Extract the image information of pressure evolution within the range of the dominant oscillation period.

[0087] Furthermore, the feature association unit is specifically used for:

[0088] According to a preset dominant frequency screening condition, data processing is performed on the first frequency data to determine the dominant frequency in the first frequency data;

[0089] According to matching the dominant frequency in the first frequency data with the peak frequency of the second frequency data, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain a characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

[0090] Furthermore, the flow field data construction unit is specifically used for:

[0091] According to the equipment structure information of the burner, the steady flow field data inside the burner in the unburned state is constructed through the steady numerical simulation method of fluid mechanics;

[0092] Using the steady flow field data as the initial field model, the flow field simulation data of the preset time length is calculated through the unsteady numerical simulation method of fluid mechanics;

[0093] Unsteady flow field data are constructed based on the flow field simulation data, and the unsteady flow field data are used as the flow field data inside the burner in an unburned state.

[0094] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described devices and units can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0095] In the several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0096] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can, for example, be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0097] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0098] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0099] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0100] As described above, the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for determining combustion oscillation characteristics, characterized in that: include: According to the equipment structure information of the burner, the flow field data inside the burner in the unburned state is constructed by fluid mechanics simulation; Based on the flow field data and in combination with preset fuel injection configuration information, the flow field of the burner in the combustion state is simulated, and pressure evolution data of each sampling point and pressure snapshot data of the entire flow field are collected; Performing Fourier transformation on the pressure evolution data to obtain first frequency data of the pressure at each sampling point and a pressure evolution process image corresponding to the first frequency data, so as to obtain a local combustion oscillation feature corresponding to the sampling point; Performing dynamic mode decomposition on the pressure snapshot data to obtain second frequency data of the full flow field pressure and a pressure dynamic mode image corresponding to the second frequency data, so as to obtain global combustion oscillation characteristics corresponding to the full flow field; Matching is performed based on the first frequency data and the second frequency data, and based on the matching result, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain a characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

2. A method for determining combustion oscillation characteristics according to claim 1, characterized in that: The first frequency data and the second frequency data both include: a peak frequency and a dominant frequency, wherein the dominant frequency is a peak frequency whose amplitude corresponding to the frequency exceeds a preset amplitude threshold.

3. A method for determining combustion oscillation characteristics according to claim 2, characterized in that: Determining the local combustion oscillation characteristics corresponding to the sampling point according to the frequency data specifically includes: Determine, according to the amplitude corresponding to the peak frequency of the sampling point, a target peak frequency corresponding to the maximum amplitude, and determine a dominant oscillation period corresponding to the target peak frequency; Extract the pressure evolution image information within the dominant oscillation period.

4. A method for determining combustion oscillation characteristics according to claim 2, characterized in that: Matching the first frequency data with the second frequency data and associating the local combustion oscillation characteristics with the global combustion oscillation characteristics after successful frequency matching according to the matching result specifically includes: According to a preset dominant frequency screening condition, performing data processing on the first frequency data to determine a dominant frequency in the first frequency data; According to matching the dominant frequency in the first frequency data with the peak frequency of the second frequency data, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

5. The method for determining combustion oscillation characteristics according to claim 1, characterized in that: According to the equipment structure information of the burner, the flow field data of the burner in the unburned state is constructed by fluid mechanics simulation, which specifically includes: According to the equipment structure information of the burner, the steady flow field data of the inside of the burner in the unburned state is constructed by means of steady numerical simulation of fluid mechanics; Using the steady flow field data as the initial field model, the flow field simulation data of a preset time length is calculated by an unsteady numerical simulation method of fluid mechanics; Unsteady flow field data are constructed based on the flow field simulation data, and the unsteady flow field data are used as flow field data inside the burner in an unburned state.

6. A combustion oscillation characteristic determination device, characterized in that: include: A flow field data construction unit, used to construct flow field data of the burner in an unburned state by fluid mechanics simulation according to the equipment structure information of the burner; A combustion simulation unit, for simulating the flow field of the burner in a combustion state based on the flow field data and in combination with preset fuel injection configuration information, and collecting pressure evolution data of each sampling point and pressure snapshot data of the entire flow field; A local feature acquisition unit, used to perform Fourier transform on the pressure evolution data to obtain first frequency data of the pressure at each sampling point and a pressure evolution process image corresponding to the first frequency data, so as to obtain a local combustion oscillation feature corresponding to the sampling point; A global feature acquisition unit, used for performing dynamic mode decomposition on the pressure snapshot data to obtain second frequency data of the full flow field pressure and a pressure dynamic mode image corresponding to the second frequency data, so as to obtain a global combustion oscillation feature corresponding to the full flow field; A feature association unit is used to match the first frequency data with the second frequency data, and according to the matching result, associate the local combustion oscillation feature with the global combustion oscillation feature that successfully matches the frequency to obtain the feature relationship between the local combustion oscillation feature and the global combustion oscillation feature.

7. The combustion oscillation characteristic determination device according to claim 6, characterized in that: The first frequency data and the second frequency data both include: a peak frequency and a dominant frequency, wherein the dominant frequency is a peak frequency whose amplitude corresponding to the frequency exceeds a preset amplitude threshold.

8. The combustion oscillation characteristic determination device according to claim 7, characterized in that: The local feature acquisition unit is specifically used for: Performing Fourier transformation on the pressure evolution data to obtain the peak frequency of each sampling point; Determine, according to the amplitude corresponding to the peak frequency of the sampling point, a target peak frequency corresponding to the maximum amplitude, and determine a dominant oscillation period corresponding to the target peak frequency; Extract the pressure evolution image information within the dominant oscillation period.

9. The combustion oscillation characteristic determination device according to claim 7, characterized in that: The feature association unit is specifically used for: According to a preset dominant frequency screening condition, performing data processing on the first frequency data to determine a dominant frequency in the first frequency data; According to matching the dominant frequency in the first frequency data with the peak frequency of the second frequency data, the local combustion oscillation characteristics with successful frequency matching are associated with the global combustion oscillation characteristics to obtain the characteristic relationship between the local combustion oscillation characteristics and the global combustion oscillation characteristics.

10. The combustion oscillation characteristic determination device according to claim 6, characterized in that: The flow field data construction unit is specifically used for: According to the equipment structure information of the burner, the steady flow field data of the inside of the burner in the unburned state is constructed by means of steady numerical simulation of fluid mechanics; Using the steady flow field data as the initial field model, the flow field simulation data of a preset time length is calculated by an unsteady numerical simulation method of fluid mechanics; Unsteady flow field data are constructed based on the flow field simulation data, and the unsteady flow field data are used as flow field data inside the burner in an unburned state.

Images