Method for determining parameters of thermoacoustic plate in annular combustor and annular combustor

Abstract

The present disclosure provides a method for determining parameters of a thermoacoustic panel in an annular combustion chamber and an annular combustion chamber, relating to the field of aeroengine technology. The method comprises: constructing a first scattered acoustic field function formed by the thermoacoustic panel in the combustion chamber, constructing a second scattered acoustic field function formed by the thermoacoustic panel in a back cavity, and constructing a third scattered acoustic field function formed by a heat source in the combustion chamber; processing the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function by Laplace transform to obtain a characteristic equation; if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instabilities, then determining the current setting parameters of the thermoacoustic panel as target setting parameters of the thermoacoustic panel in the annular combustion chamber, wherein the setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel. The present disclosure can effectively suppress low-frequency modal combustion instabilities in an annular combustion chamber.

Description

Technical Field

[0001] The present disclosure relates to the technical field of aero-engines, and in particular to a method for determining parameters of a thermoacoustic panel in an annular combustion chamber and the annular combustion chamber. Background Art

[0002] Combustion instability, caused by the full coupling of acoustic waves and unsteady heat release, is a common problem in combustion chambers of aerospace propulsion systems. It is accompanied by large-scale pressure oscillations, which can, in severe cases, cause wall ablation and even damage the system.

[0003] In related technologies, passive control methods, unaffected by heat sources, have become the primary means of suppressing combustion instabilities in multi-operating environments of aircraft engines. The idea behind these passive control methods is to increase the sound absorption capacity of the boundary to suppress combustion oscillations. Traditional passive control methods, such as the classic Helmholtz resonator or a non-localized acoustic liner structure consisting of a perforated plate and back cavity that also functions as a cooling wall, can effectively suppress most high-frequency modal combustion instabilities within annular combustion chambers.

[0004] Conventional passive control methods satisfy the acoustic impedance matching condition, which requires that the acoustic wave satisfy the acoustic pressure jump or acoustic velocity jump condition on both sides of the Helmholtz resonator or non-local acoustic liner, thereby generating acoustic energy dissipation and suppressing high-frequency modal combustion instabilities in annular combustors. However, for low-frequency modes, based on the principle of acoustic impedance matching, traditional passive control methods have difficulty effectively suppressing low-frequency modal combustion instabilities in geometrically constrained annular combustors. Summary of the Invention

[0005] The present disclosure provides a method for determining parameters of a thermoacoustic panel in an annular combustion chamber and an annular combustion chamber, thereby effectively suppressing low-frequency modal combustion instabilities in the annular combustion chamber. The technical solutions of the present disclosure are as follows:

[0006] According to a first aspect of the present disclosure, a method for determining parameters of a thermoacoustic panel in an annular combustion chamber is provided. The method is used to determine the setting parameters of the thermoacoustic panel arranged axially on the outer wall and / or inner wall of the combustion chamber. The method comprises:

[0007] Constructing a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the combustion chamber, constructing a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity, and constructing a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber in the combustion chamber;

[0008] Processing the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function by Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source;

[0009] If the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instabilities, the current setting parameters of the thermoacoustic panel are determined as target setting parameters of the thermoacoustic panel in the annular combustion chamber, wherein the setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel.

[0010] According to a second aspect of the present disclosure, a device for determining parameters of a thermoacoustic panel in an annular combustion chamber is provided, the device being used to determine setting parameters of the thermoacoustic panel arranged axially on the outer wall and / or inner wall of the combustion chamber, the device comprising:

[0011] a construction module configured to construct a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber, construct a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the back cavity, and construct a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber;

[0012] a data processing module configured to process the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function through Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source;

[0013] The determination module is configured to determine the current setting parameters of the thermoacoustic panel as target setting parameters of the thermoacoustic panel in the annular combustion chamber if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the low-frequency modal combustion instability suppression requirements, wherein the setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel.

[0014] According to a third aspect of the present disclosure, there is provided an electronic device, including:

[0015] processor; and

[0016] Memory for storing programs,

[0017] The program includes instructions, which, when executed by the processor, cause the processor to perform the method as described in the first aspect.

[0018] According to a fourth aspect of the present disclosure, a non-transitory computer-readable storage medium storing computer instructions is provided, wherein the computer instructions are used to cause the computer to execute the method according to the first aspect.

[0019] According to a fifth aspect of the present disclosure, there is provided an annular combustion chamber, comprising:

[0020] The housing comprises an outer shell wall and an inner shell wall spaced apart in the axial direction;

[0021] A combustion chamber is arranged between the outer shell wall and the inner shell wall, and the combustion chamber includes an outer wall and an inner wall arranged at intervals along the axial direction, wherein the outer wall and / or the inner wall are axially provided with a thermoacoustic plate, and the setting parameters of the thermoacoustic plate are determined based on the method described in the first aspect.

[0022] The parameter determination method of the thermoacoustic plate in the annular combustion chamber and the annular combustion chamber provided in the embodiments of the present application can simultaneously consider the coupling effect of the thermoacoustic plate and the heat source, establish a combustion instability model of the coupling effect of the thermoacoustic plate and the heat source, and update the setting parameters of the thermoacoustic plate by comparing the solution results of the characteristic equation of the combustion instability model of the coupling effect of the thermoacoustic plate and the heat source, as well as the low-frequency modal combustion instability suppression conditions, to determine the setting parameters of the thermoacoustic plate in the annular combustion chamber that can more effectively suppress the low-frequency modal combustion instability, so as to facilitate the effectiveness and reliability of the suppression of low-frequency modal combustion instability of the annular combustion chamber containing the thermoacoustic plate during actual operation. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Further details, features and advantages of the present disclosure are disclosed in the following description of exemplary embodiments in conjunction with the accompanying drawings, in which:

[0024] Figure 1 A schematic cross-sectional view of an annular combustion chamber according to an exemplary embodiment of the present disclosure is shown;

[0025] Figure 2 A schematic diagram showing temperature changes on both sides of a thermoacoustic panel according to an exemplary embodiment of the present disclosure is shown;

[0026] Figure 3 A schematic cross-sectional view of another annular combustion chamber according to an exemplary embodiment of the present disclosure is shown;

[0027] Figure 4 A schematic cross-sectional view of another annular combustion chamber according to an exemplary embodiment of the present disclosure is shown;

[0028] Figure 5 A flow chart showing a method for determining parameters of a thermoacoustic panel in an annular combustion chamber according to an exemplary embodiment of the present disclosure is provided;

[0029] Figure 6 A schematic block diagram showing a device for determining parameters of a thermoacoustic panel in an annular combustion chamber according to an exemplary embodiment of the present disclosure is shown;

[0030] Figure 7A schematic block diagram of a chip according to an exemplary embodiment of the present disclosure is shown;

[0031] Figure 8 A structural block diagram of an exemplary electronic device that can be used to implement the embodiments of the present disclosure is shown. DETAILED DESCRIPTION

[0032] The following describes embodiments of the present disclosure in more detail with reference to the accompanying drawings. Although certain embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be construed as limited to the embodiments described herein. Rather, these embodiments are provided to provide a more thorough and complete understanding of the present disclosure. It should be understood that the drawings and embodiments of the present disclosure are for illustrative purposes only and are not intended to limit the scope of protection of the present disclosure.

[0033] It should be understood that the various steps described in the method embodiments of the present disclosure may be performed in different orders and / or in parallel. In addition, the method embodiments may include additional steps and / or omit the steps shown. The scope of the present disclosure is not limited in this respect.

[0034] The term "including" and its variations used in this document are open inclusions, that is, "including but not limited to". The term "based on" means "based at least in part on". The term "one embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one other embodiment"; the term "some embodiments" means "at least some embodiments". The relevant definitions of other terms will be given in the description below. It should be noted that the concepts of "first", "second", etc. mentioned in this disclosure are only used to distinguish different devices, modules or units, and are not used to limit the order or interdependence of the functions performed by these devices, modules or units.

[0035] It should be noted that the modifications of "one" and "multiple" mentioned in the present disclosure are illustrative rather than restrictive, and those skilled in the art should understand that unless otherwise clearly indicated in the context, they should be understood as "one or more".

[0036] The names of the messages or information exchanged between multiple devices in the embodiments of the present disclosure are only used for illustrative purposes and are not used to limit the scope of these messages or information.

[0037] In the related art, for the problem of combustion instability, traditional passive control methods that rely on acoustic impedance matching conditions, such as Helmholtz resonators, or non-local acoustic lining structures consisting of perforated plates and back cavities that also serve as cooling walls, can usually have a good suppressive effect on high-frequency modal combustion instabilities in annular combustion chambers with restricted geometric structures.

[0038] However, within the geometrically constrained combustion chamber, traditional passive control methods struggle to effectively control low-frequency modal combustion instabilities within very small geometries. This is because the fundamental principle behind traditional passive control methods for acoustic energy dissipation is to generate a jump in acoustic pressure or velocity across a Helmholtz resonator or acoustic liner. This leads to an acoustic energy absorption limit for these methods. Suppressing low-frequency combustion instabilities within smaller geometries clearly requires a fundamental breakthrough in principle to overcome this absorption limit.

[0039] To address the above-mentioned issues, recent research has found that when sound waves pass through a thermoacoustic panel, while the viscosity of the panel produces a pressure drop, the significant temperature difference on both sides of the panel also results in thermoacoustic energy conversion, which in turn produces an additional sound velocity jump. This means that both a sound pressure jump and a sound velocity jump are achieved simultaneously on both sides of the panel, thereby breaking the sound absorption limit. Based on the benefits of sound absorption performance brought about by breakthroughs in the sound absorption principle of thermoacoustic panels, and the physical characteristics of the significant temperature difference naturally existing between the high-temperature gas in the combustion chamber and the cooling gas in the back cavity, the present disclosure provides an annular combustion chamber including a thermoacoustic panel, wherein, for example, Figure 1 As shown, Figure 1 A schematic cross-sectional view of an annular combustion chamber provided by an embodiment of the present disclosure is shown, comprising:

[0040] The housing comprises an outer shell wall 1011 and an inner shell wall 1012 spaced apart in the axial direction;

[0041] The combustion chamber is arranged between the 1011 and the inner shell wall 1012, and the combustion chamber includes an outer wall 1021 and an inner wall 1022 spaced apart in the axial direction, wherein the outer wall 1021 and / or the inner wall 1022 includes a thermoacoustic plate 103 arranged in the axial direction. Specifically, it can be determined based on actual needs. The embodiment of the present disclosure does not limit this. For example, Figure 1 Schematic diagram of an annular combustion chamber is shown in FIG. 1 , in which the outer wall surface 1021 and the inner wall surface 1022 of the combustion chamber are both axially provided with thermoacoustic plates 103 .

[0042] It should be noted that if Figure 1 As shown, the annular combustion chamber also includes a swirler heat source channel 104 radially arranged on the outer wall 1011 of the shell and the outer wall 1021 of the combustion chamber, which is used to introduce the heat source provided by the swirler into the combustion chamber through the heat source mounting surface m; wherein, there is a back cavity between the shell and the combustion chamber, which is used to provide constant temperature cooling gas for the annular combustion chamber; it can be understood that in the embodiment of the present disclosure, the heat source mounting surface is a partial wall surface radially arranged in the combustion chamber.

[0043] For example, Figure 2As shown, Figure 2 The temperature variation diagram of the two sides of the thermoacoustic plate is shown, where the temperature of the constant high temperature gas on the combustion chamber side c is T c The temperature of the constant low-temperature gas on the back cavity side b is T b , where the temperature of the thermoacoustic panel near the combustion chamber is higher than the temperature of the back cavity, where the horizontal axis represents the radial coordinate and the vertical axis represents the temperature.

[0044] In summary, the annular combustion chamber provided by the embodiment of the present disclosure can utilize the natural temperature difference formed by the cooling environment in the back cavity and the high temperature conditions in the combustion chamber. By arranging thermoacoustic plates on the outer wall and / or inner wall of the combustion chamber, the sound waves can simultaneously generate sound pressure jumps and sound velocity jumps on both sides of the thermoacoustic plates, thereby realizing the dissipation of sound energy caused by the sound pressure jumps and sound velocity jumps generated by the sound waves on both sides of the thermoacoustic plates, and more effective control of low-frequency combustion instabilities in the annular combustion chamber; at the same time, since the natural temperature difference formed by the cooling environment in the back cavity and the high temperature conditions in the combustion chamber can be utilized, the sound pressure jumps and sound velocity jumps can be generated on both sides of the thermoacoustic plates at the same time, without the need for additional temperature control equipment, so that the suppression of low-frequency combustion instabilities is not restricted by the combustion chamber structure, thereby reducing the difficulty of suppressing low-frequency combustion instabilities.

[0045] In an optional embodiment, as Figure 1 As shown, the outer wall surface 1021 and the inner wall surface 1022 of the combustion chamber may further include a perforated plate 105 arranged along the axial direction to further suppress noise.

[0046] Optional, such as Figure 3 As shown, Figure 3 A schematic cross-sectional view of another annular combustion chamber provided by an embodiment of the present disclosure is shown. Figure 3 In the annular combustion chamber shown, the outer wall surface 1021 of the combustion chamber includes a thermoacoustic plate 103 arranged in the axial direction; Figure 4 As shown, Figure 4 A schematic cross-sectional view of another annular combustion chamber provided by an embodiment of the present disclosure is shown. Figure 4 In the annular combustion chamber shown, the inner wall surface 1022 of the combustion chamber includes a thermoacoustic plate 103 arranged in the axial direction.

[0047] It should be noted that, in the embodiment of the present disclosure, the thermoacoustic plate may be a porous material provided with gaps, wherein the material of the thermoacoustic plate may be determined based on actual needs, and the embodiment of the present disclosure does not limit this.

[0048] The present disclosure provides a method for determining parameters of a thermoacoustic plate in an annular combustion chamber. The method is used to determine the setting parameters of the thermoacoustic plate arranged axially on the outer wall surface and / or the inner wall surface of the combustion chamber. The method can be applied to electronic devices with data processing capabilities. The electronic devices can be terminal devices or servers, wherein the terminal devices can be computers, notebooks, tablet computers, etc. Figure 5 As shown, the method includes:

[0049] Step S501, constructing a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber, constructing a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the back cavity, and constructing a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber;

[0050] Step S502 , processing the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function by Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source;

[0051] Step S503 , if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instability, then the current setting parameters of the thermoacoustic panel are determined as the target setting parameters of the thermoacoustic panel in the annular combustion chamber;

[0052] The setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel.

[0053] In summary, the method for determining the parameters of the thermoacoustic plate in the annular combustion chamber provided by the embodiment of the present disclosure can simultaneously consider the coupling effect of the thermoacoustic plate and the heat source, establish a combustion instability model of the coupling effect of the thermoacoustic plate and the heat source, and update the setting parameters of the thermoacoustic plate by comparing the solution results of the characteristic equation of the combustion instability model of the coupling effect of the thermoacoustic plate and the heat source, as well as the low-frequency modal combustion instability suppression conditions, so as to determine the setting parameters of the thermoacoustic plate that can more effectively suppress the low-frequency modal combustion instability in the annular combustion chamber, so as to facilitate the effectiveness and reliability of the suppression of low-frequency modal combustion instability of the annular combustion chamber containing the thermoacoustic plate during actual operation.

[0054] Step S501, constructing a first scattered sound field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber, constructing a second scattered sound field function describing the equivalent monopole source of the thermoacoustic plate formed in the back cavity, and constructing a third scattered sound field function describing the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber.

[0055] In the embodiment of the present disclosure, the first scattered sound field function formed by the equivalent monopole source of the thermoacoustic plate in the combustion chamber is related to the sound velocity of the thermoacoustic plate on one side of the combustion chamber, and the second scattered sound field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity is related to the sound velocity of the thermoacoustic plate on one side of the back cavity; therefore, before determining the first scattered sound field function and the second scattered sound field function, the impedance boundary condition between the sound velocity and the sound pressure of the thermoacoustic plate on one side of the combustion chamber and the impedance boundary condition between the sound velocity and the sound pressure of the thermoacoustic plate on the back cavity side can be determined.

[0056] In an optional embodiment, the impedance boundary conditions between the sound velocity and sound pressure of the thermoacoustic plate on one side of the combustion chamber, as well as the impedance boundary conditions between the sound velocity and sound pressure of the thermoacoustic plate on the side of the back cavity can be determined based on the engineering parameters of the annular combustion chamber and the current setting parameters of the thermoacoustic plate. The impedance boundary conditions between the sound velocity and sound pressure of the thermoacoustic plate on one side of the combustion chamber, as well as the impedance boundary conditions between the sound velocity and sound pressure of the thermoacoustic plate on the side of the back cavity can be determined in combination with the engineering parameters of the annular combustion chamber and the current setting parameters of the thermoacoustic plate, so as to construct a more accurate combustion instability model that describes the coupling effect of the thermoacoustic plate and the heat source.

[0057] The process of determining the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the combustion chamber and the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on the back cavity side based on the engineering parameters of the annular combustion chamber and the current setting parameters of the thermoacoustic plate may include:

[0058] Considering the constant temperature difference between the constant high-temperature gas in the combustion chamber and the constant low-temperature cooling gas in the back cavity on both sides of the thermoacoustic plate, when the sound wave passes through the thermoacoustic plate, it can interact with the medium inside the thermoacoustic plate to produce thermoacoustic energy exchange. According to the linear thermoacoustic theory, the thermoacoustic plate can be assumed to be a plate-like structure composed of parallel thin tubes arranged along the radial direction. The sound propagation in the thermoacoustic plate can be described as:

[0059]

[0060] In Formula 1 and Formula 2, i represents the imaginary unit, ω represents the frequency, and ρ m represents the average density of the medium in the thermoacoustic plate, S represents the cross-sectional area of the waveguide, φ represents the porosity of the pores in the thermoacoustic plate, represents the velocity disturbance, represents the pressure disturbance, r represents the radial coordinate of the combustion chamber, d represents the ordinary differential operator, γ represents the specific heat ratio, p m represents the average pressure of the medium in the thermoacoustic plate, σ represents the Prandtl number, T m Indicates the average temperature.

[0061] f υ and f κ is a frequency-dependent function that describes the viscous and thermal coupling effects between the gas and the thermoacoustic plate, respectively. For a medium with kinematic viscosity υ and thermal diffusivity κ, the viscous and thermal boundary layer thicknesses can be defined as and Then for radius r m The frequency-dependent function f υ and f κ It can be:

[0062]

[0063] In Equation 3, J0 represents the zero-order Bessel function of the first kind, and J1 represents the first-order Bessel function. Term A in Equation (1) describes the inertial and viscous effects of the gas in the thermoacoustic plate, while term B in Equation (2) describes the compression effect of the gas in the thermoacoustic plate and the loss caused by heat conduction. The additional term C represents the additional velocity jump caused by the thermoacoustic energy conversion within the thermoacoustic plate due to the temperature difference.

[0064] Furthermore, by simplifying Formula 2 through the quasi-isothermal assumption, we can obtain:

[0065]

[0066] In formula 4, and represent the sound pressure and sound velocity of the thermoacoustic plate on the back cavity side, and denote the sound pressure and sound velocity of the thermoacoustic plate on one side of the combustion chamber, R denotes the acoustic resistivity of the thermoacoustic plate material, and η denotes the temperature ratio of the high-temperature gas in the combustion chamber to the cooling gas in the back cavity. The functional relationship between the sound pressure and sound velocity on both sides of the thermoacoustic plate can be described as:

[0067]

[0068] In formula 5, It represents the relationship between the sound pressure and the sound velocity on the side of the thermoacoustic plate in the combustion chamber. In formula 6, represents the relationship between the sound pressure and sound velocity on the back cavity side of the thermoacoustic plate. Therefore, Formula 5 and Formula 6 can be determined as the acoustic response function of the thermoacoustic plate, where Formula 5 is the impedance boundary condition between the sound velocity and sound pressure of the thermoacoustic plate on the combustion chamber side, and Formula 6 is the impedance boundary condition between the sound velocity and sound pressure of the thermoacoustic plate on the back cavity side.

[0069] In an optional embodiment, since the presence of the thermoacoustic plate forms an acoustic impedance boundary condition at this spatial location, its effect on the acoustic field is consistent with the response characteristics of a pulsating spherical source, i.e., a monopole source. Based on generalized function theory, the process of constructing a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate within the combustion chamber may include:

[0070] Based on the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the combustion chamber, a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber is constructed. The first scattered acoustic field function is:

[0071]

[0072] In formula 7, represents the pressure disturbance generated in the combustion chamber by the equivalent monopole source describing the thermoacoustic plate, represents the position information of the equivalent monopole source describing the thermoacoustic plate at the observation point in three-dimensional space, r, θ, x represent the radial coordinate, circumferential coordinate, and axial coordinate in the cylindrical coordinate system respectively, t represents the time of the observation point, N represents the number of equivalent monopole sources describing the thermoacoustic plate, j represents the sequence number of the equivalent monopole source describing the thermoacoustic plate, S represents the area of the equivalent monopole source describing the thermoacoustic plate, represents the gas density in the combustion chamber, represents the normal acoustic particle velocity of the equivalent monopole source of the thermoacoustic plate on the combustion chamber side, is determined based on the impedance boundary condition between the acoustic velocity and the acoustic pressure of the thermoacoustic panel on one side of the combustion chamber, represents the source coordinates of the equivalent monopole source, represents the source time of the equivalent monopole source, and G represents the Green's function; it can be understood that in the embodiment of the present disclosure, represents the source point coordinates of the equivalent monopole source describing the thermoacoustic plate in three-dimensional space, The source time of the equivalent monopole source describing the thermoacoustic plate is shown.

[0073] Similarly, in an optional embodiment, the process of constructing a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic panel in the back cavity may include:

[0074] Based on the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the back cavity, a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity is constructed, wherein the second scattered acoustic field function is:

[0075]

[0076] In formula 8, represents the pressure disturbance generated by the equivalent monopole source of the thermoacoustic plate in the back cavity, represents the normal acoustic particle velocity of the equivalent monopole source of the thermoacoustic plate on the back cavity side, It is determined based on the impedance boundary condition between the acoustic velocity and the acoustic pressure of the thermal acoustic panel on one side of the back cavity.

[0077] In an optional embodiment, the influence of the heat source in the combustion chamber on the sound field in the combustion chamber can also be described by an equivalent monopole source. Then, the process of constructing a third scattered sound field function formed in the combustion chamber by the equivalent monopole source describing the heat source in the combustion chamber may include:

[0078] A flame combustion function describing the unsteady heat release characteristics of the flame of the heat source is determined, and based on the flame combustion function, a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber is constructed in the combustion chamber. The third scattered acoustic field function is:

[0079]

[0080] in, represents the scattered sound pressure generated in the combustion chamber by the equivalent monopole source describing the heat source in the combustion chamber, N f represents the number of equivalent monopole sources describing the heat source in the combustion chamber, q represents the sequence number of the equivalent monopole source describing the heat source, S′ represents the area of the equivalent monopole source describing the heat source, represents the acoustic particle velocity at the interface between the heat source outlet and the combustion chamber inlet, where It is obtained by describing the unsteady heat release characteristics of the flame of the heat source through the flame combustion function. It can be understood that in the embodiment of the present disclosure, represents the source point coordinates of the equivalent monopole source describing the heat source in the combustion chamber in three-dimensional space, represents the source time of the equivalent monopole source describing the heat source in the combustion chamber.

[0081] It should be noted that, in the embodiment of the present disclosure, the first scattered sound field function, the second scattered sound field function and the third scattered sound field function may be determined as a combustion instability model of the coupling effect between the thermoacoustic panel and the heat source.

[0082] In an optional embodiment, a fourth scattered acoustic field function formed by the interaction between the thermoacoustic panel and the unsteady heat release of the flame of the heat source in the combustion chamber can be determined based on a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic panel in the combustion chamber and a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber. The fourth scattered acoustic field function is:

[0083]

[0084] In formula 10, It represents the pressure disturbance formed in the combustion chamber by the coupling between the equivalent monopole source describing the heat source in the combustion chamber and the equivalent monopole source describing the thermoacoustic plate.

[0085] In step S502, the first scattered acoustic field function, the second scattered acoustic field function and the third scattered acoustic field function are processed by Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source.

[0086] In an optional embodiment, the process of processing the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function by Laplace transform to obtain the characteristic equation of the coupling effect between the thermoacoustic panel and the heat source may include: performing Laplace transform on a fourth scattered acoustic field function formed by the interaction of the unsteady heat release of the flame of the thermoacoustic panel and the heat source in the combustion chamber, and a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic panel in the back cavity, to obtain the characteristic equation of the coupling effect between the thermoacoustic panel and the heat source, wherein the characteristic equation is:

[0087] X[ζ,ω]P=0;(Formula 11)

[0088] In Equation 11, X represents the coefficient matrix of the characteristic equation, ζ represents the amplitude of the nonlinear limit cycle, ω represents the complex frequency of the thermoacoustic oscillation, and the vector P contains the complex pressure at each monopole source position in the combustion chamber.

[0089] In step S503, if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instabilities, the current setting parameters of the thermoacoustic panel are determined as target setting parameters of the thermoacoustic panel in the annular combustion chamber.

[0090] In an embodiment of the present disclosure, the setting parameters include setting position parameters and / or size parameters of the thermoacoustic plate, wherein the setting position parameters may include setting position parameters along the axial direction on the outer wall surface and / or the inner wall surface of the combustion chamber, and the size parameters may include the thickness, porosity and pore radius of the thermoacoustic plate; the characteristic solution of the characteristic equation may include the real part, the imaginary part and the amplitude of the nonlinear limit cycle, wherein the real part is the thermoacoustic oscillation frequency.

[0091] In an optional embodiment, determining whether the characteristic solution of the characteristic equation, which indicates the current setting parameters of the thermoacoustic panel, meets the requirement for suppressing low-frequency modal combustion instabilities may include: wherein the characteristic solution of the characteristic equation may include the real part of the complex frequency of the thermoacoustic oscillation, the imaginary part of the complex frequency of the thermoacoustic oscillation, and the amplitude of the nonlinear limit cycle; if the real part of the complex frequency of the thermoacoustic oscillation is within a preset low-frequency range and the imaginary part of the complex frequency of the thermoacoustic oscillation is greater than 0, determining that the characteristic solution of the characteristic equation indicates the current setting parameters of the thermoacoustic panel meets the requirement for suppressing low-frequency modal combustion instabilities; or, if the real part of the complex frequency of the thermoacoustic oscillation is within the preset low-frequency range and the amplitude of the nonlinear limit cycle is less than a preset parameter value, determining that the characteristic solution of the characteristic equation indicates the current setting parameters of the thermoacoustic panel meets the requirement for suppressing low-frequency modal combustion instabilities; wherein the preset low-frequency range may be determined based on actual needs and is not limited in the embodiment of the present disclosure. For example, the preset low-frequency range may be 20 Hz to 200 Hz. The preset parameter value is used to characterize the stress that the machine can withstand. The preset parameter value may be determined based on actual engineering practice and is not limited in the embodiment of the present disclosure.

[0092] In an optional embodiment, if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel do not meet the requirements for suppressing low-frequency modal combustion instabilities, then updated setting parameters of the thermoacoustic panel are determined within a parameter threshold interval of the thermoacoustic panel, and based on the updated setting parameters of the thermoacoustic panel, an updated characteristic equation for the coupling between the thermoacoustic panel and the heat source is obtained. Furthermore, the above process of obtaining the updated characteristic equation is repeated until the characteristic solution of the updated characteristic equation indicates that the updated setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instabilities, and the updated setting parameters of the thermoacoustic panel are determined as the target setting parameters of the thermoacoustic panel in the annular combustion chamber. By traversing the parameter threshold interval, the parameters of the thermoacoustic panel that meet the requirements for suppressing low-frequency combustion instabilities can be determined, thereby improving the suppression effect of the annular combustion chamber equipped with the thermoacoustic panel on low-frequency combustion instabilities.

[0093] It can be understood that in the embodiment of the present disclosure, the parameter threshold interval may include a setting position parameter threshold interval and / or a size parameter interval. The setting position parameter threshold interval can be determined based on the actual size of the annular combustion chamber, and the size parameter interval can be determined based on actual needs. The embodiment of the present disclosure does not limit this.

[0094] The above mainly introduces the solution provided by the embodiment of the present disclosure from the perspective of an electronic device. It is understandable that, in order to realize the above functions, the electronic device includes a hardware structure and / or software module corresponding to the execution of each function. Those skilled in the art should easily realize that, in combination with the units and algorithm steps of each example described in the embodiments disclosed herein, the present disclosure can be implemented in the form of hardware or a combination of hardware and computer software. Whether a function is executed in a hardware or computer software driven hardware manner depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present disclosure.

[0095] The embodiments of the present disclosure can divide the functional units of the electronic device according to the above method examples. For example, each functional module can be divided according to each function, or two or more functions can be integrated into one processing module. The above integrated modules can be implemented in the form of hardware or in the form of software functional modules. It should be noted that the division of modules in the embodiments of the present disclosure is schematic and is only a logical function division. In actual implementation, there may be other division methods.

[0096] In the case of dividing each functional module according to each function, an exemplary embodiment of the present disclosure provides a parameter determination device for a thermoacoustic plate in an annular combustion chamber. The parameter determination device for a thermoacoustic plate in an annular combustion chamber can be an electronic device or a chip applied to an electronic device. Figure 6 The following is a schematic block diagram of the functional modules of a device for determining parameters of a thermoacoustic panel in an annular combustion chamber according to an exemplary embodiment of the present disclosure, wherein the device is used to determine the setting parameters of the thermoacoustic panel arranged axially on the outer wall and / or inner wall of the combustion chamber. Figure 6 As shown, the parameter determination device 600 of the thermoacoustic panel in the annular combustion chamber includes:

[0097] A construction module 601 is configured to construct a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic panel formed in the combustion chamber, construct a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic panel formed in the back cavity, and construct a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber;

[0098] A data processing module 602 is configured to process the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function through Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source;

[0099] The determination module 603 is configured to determine the current setting parameters of the thermoacoustic panel as target setting parameters of the thermoacoustic panel in the annular combustion chamber if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instabilities, wherein the setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel.

[0100] Optional, such as Figure 6 As shown, the apparatus further includes a pre-processing module 604 configured to:

[0101] Based on the engineering parameters of the annular combustion chamber and the current setting parameters of the thermoacoustic plate, the impedance boundary conditions between the sound velocity and the sound pressure of the thermoacoustic plate on one side of the combustion chamber and the impedance boundary conditions between the sound velocity and the sound pressure of the thermoacoustic plate on the side of the back cavity are determined.

[0102] Optionally, the building module 601 is configured to:

[0103] Based on the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the combustion chamber, a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber is constructed. The first scattered acoustic field function is:

[0104]

[0105] in, represents the pressure disturbance generated in the combustion chamber by the equivalent monopole source describing the thermoacoustic plate, represents the position information of the equivalent monopole source describing the thermoacoustic plate at the observation point in three-dimensional space, r, θ, x represent the radial coordinate, circumferential coordinate, and axial coordinate in the cylindrical coordinate system respectively, t represents the time of the observation point, N represents the number of equivalent monopole sources describing the thermoacoustic plate, j represents the sequence number of the equivalent monopole source describing the thermoacoustic plate, S represents the area of the equivalent monopole source describing the thermoacoustic plate, represents the gas density in the combustion chamber, represents the normal acoustic particle velocity of the equivalent monopole source of the thermoacoustic plate on the combustion chamber side, is determined based on the impedance boundary condition between the acoustic velocity and the acoustic pressure of the thermoacoustic panel on one side of the combustion chamber, represents the source coordinates of the equivalent monopole source, represents the source time of the equivalent monopole source, and G represents the Green's function.

[0106] Optionally, the construction module 601 is configured to:

[0107] Based on the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the back cavity, a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity is constructed. The second scattered acoustic field function is:

[0108]

[0109] in, Describe the pressure disturbance in the back cavity caused by the equivalent monopole source of the thermoacoustic plate, represents the normal acoustic particle velocity of the equivalent monopole source of the thermoacoustic plate on the back cavity side, It is determined based on the impedance boundary condition between the acoustic velocity and the acoustic pressure of the thermal acoustic panel on one side of the back cavity.

[0110] Optionally, the construction module 601 is configured to:

[0111] Determine the flame combustion function that describes the unsteady heat release characteristics of the flame of the heat source;

[0112] Based on the flame combustion function, a third scattered acoustic field function is constructed to describe the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber. The third scattered acoustic field function is:

[0113]

[0114] in, represents the scattered sound pressure generated in the combustion chamber by the equivalent monopole source describing the heat source in the combustion chamber, N f represents the number of equivalent monopole sources describing the heat source in the combustion chamber, q represents the sequence number of the equivalent monopole source describing the heat source, S′ represents the area of the equivalent monopole source describing the heat source, Represents the acoustic particle velocity at the interface between the heat source outlet and the combustion chamber inlet.

[0115] Optional, such as Figure 6 As shown, the apparatus further includes an updating module 605 configured to:

[0116] If the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel do not meet the low-frequency modal combustion instability suppression requirement, then determining updated setting parameters of the thermoacoustic panel within a parameter threshold interval of the thermoacoustic panel, and obtaining an updated characteristic equation of the coupling between the thermoacoustic panel and the heat source based on the updated setting parameters of the thermoacoustic panel;

[0117] The above process is repeated until the characteristic solution of the updated characteristic equation indicates the updated setting parameters of the thermoacoustic panel, which meets the requirements for suppressing low-frequency modal combustion instabilities, and the updated setting parameters of the thermoacoustic panel are determined as the target setting parameters of the thermoacoustic panel in the annular combustion chamber.

[0118] Figure 7 FIG. 1 shows a schematic block diagram of a chip according to an exemplary embodiment of the present disclosure. Figure 7 As shown, the chip 700 includes one or more (including two) processors 701 and a communication interface 702. The communication interface 702 can support the electronic device to execute the data transmission and reception steps in the above-mentioned method for determining the parameters of the thermoacoustic panel in the annular combustion chamber, and the processor 701 can support the electronic device to execute the data processing steps in the above-mentioned method for determining the parameters of the thermoacoustic panel in the annular combustion chamber.

[0119] Optional, such as Figure 7 As shown, the chip 700 also includes a memory 703, which may include a read-only memory and a random access memory, and provides operation instructions and data to the processor. Part of the memory may also include a non-volatile random access memory (NVRAM).

[0120] In some embodiments, as Figure 7 As shown, the processor 701 performs corresponding operations by calling the operation instructions stored in the memory (the operation instructions may be stored in the operating system). The processor 701 controls the processing operations of any electronic device, and the processor may also be called a central processing unit (CPU). The memory 703 may include a read-only memory and a random access memory, and provides instructions and data to the processor 701. A portion of the memory 703 may also include NVRAM. For example, in an application, the memory, the communication interface, and the memory are coupled together through a bus system, wherein the bus system may include a power bus, a control bus, and a status signal bus in addition to a data bus. However, for the sake of clarity, in Figure 7 Various buses are labeled as bus system 704 .

[0121] The methods disclosed in the above embodiments of the present disclosure can be applied to or implemented by a processor. The processor may be an integrated circuit chip with signal processing capabilities. During implementation, each step of the above method can be completed by hardware integrated logic circuits in the processor or by software instructions. The above processor may be a general-purpose processor, a digital signal processor (DSP), an ASIC, a field-programmable gate array (FPGA), or other programmable logic device, a discrete gate or transistor logic device, or a discrete hardware component. The methods, steps, and logic block diagrams disclosed in the embodiments of the present disclosure can be implemented or executed. The general-purpose processor may be a microprocessor or any conventional processor. The steps of the methods disclosed in conjunction with the embodiments of the present disclosure can be directly implemented and executed by a hardware decoding processor, or by a combination of hardware and software modules in the decoding processor. The software module can be located in a storage medium well-known in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, etc. The storage medium is located in the memory, and the processor reads the information in the memory and, in conjunction with its hardware, completes the steps of the above method.

[0122] The exemplary embodiments of the present disclosure further provide an electronic device, comprising: at least one processor; and a memory communicatively connected to the at least one processor. The memory stores a computer program executable by the at least one processor, the computer program being configured to cause the electronic device to perform a method according to an exemplary embodiment of the present disclosure when executed by the at least one processor.

[0123] Exemplary embodiments of the present disclosure further provide a non-transitory computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor of a computer, is used to cause the computer to perform a method according to an embodiment of the present disclosure.

[0124] Exemplary embodiments of the present disclosure further provide a computer program product, including a computer program, wherein when the computer program is executed by a processor of a computer, it is used to cause the computer to perform the method according to the embodiment of the present disclosure.

[0125] refer to Figure 8, a structural block diagram of an electronic device 800 that can be used as the present disclosure will now be described, which is an example of a hardware device that can be applied to various aspects of the present disclosure. The electronic device is intended to represent various forms of digital electronic computer devices, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processing, cellular phones, smart phones, wearable devices and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely examples and are not intended to limit the implementation of the present disclosure described and / or required herein.

[0126] like Figure 8 As shown, the electronic device 800 includes a computing unit 801, which can perform various appropriate actions and processes according to a computer program stored in a read-only memory (ROM) 802 or a computer program loaded from a storage unit 808 into a random access memory (RAM) 803. In the RAM 803, various programs and data required for the operation of the electronic device 800 can also be stored. The computing unit 801, the ROM 802, and the RAM 803 are connected to each other via a bus 804. An input / output (I / O) interface 805 is also connected to the bus 804.

[0127] Multiple components within electronic device 800 are connected to I / O interface 805, including an input unit 806, an output unit 807, a storage unit 808, and a communication unit 809. Input unit 806 can be any type of device capable of inputting information into electronic device 800. Input unit 806 can receive input digital or character information and generate signal inputs related to user settings and / or function control of the electronic device. Output unit 807 can be any type of device capable of presenting information and may include, but is not limited to, a display, a speaker, a video / audio output terminal, a vibrator, and / or a printer. Storage unit 808 may include, but is not limited to, a magnetic disk or an optical disk. Communication unit 809 allows electronic device 800 to exchange information / data with other devices via computer networks such as the Internet and / or various telecommunication networks and may include, but is not limited to, a modem, a network card, an infrared communication device, a wireless communication transceiver and / or chipset, such as a Bluetooth™ device, a WiFi device, a WiMax device, a cellular communication device, and / or the like.

[0128] The computing unit 801 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the computing unit 801 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various dedicated artificial intelligence (AI) computing chips, various computing units for running machine learning model algorithms, digital signal processors (DSPs), and any appropriate processors, controllers, microcontrollers, etc. The computing unit 801 performs the various methods and processes described above. For example, in some embodiments, the method of the embodiment of the present disclosure may be implemented as a computer software program, which is tangibly included in a machine-readable medium, such as a storage unit 808. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 800 via the ROM 802 and / or the communication unit 809. In some embodiments, the computing unit 801 may be configured to perform the method of the embodiment of the present disclosure by any other appropriate means (e.g., by means of firmware).

[0129] The program code for implementing the method of the present disclosure can be written in any combination of one or more programming languages. These program codes can be provided to a processor or controller of a general-purpose computer, a special-purpose computer, or other programmable data processing device so that when the program code is executed by the processor or controller, the functions / operations specified in the flow chart and / or block diagram are implemented. The program code can be executed entirely on the machine, partially on the machine, as a stand-alone software package, partially on the machine and partially on a remote machine, or entirely on a remote machine or server.

[0130] In the context of the present disclosure, a machine-readable medium can be a tangible medium that can contain or store a program for use by or in conjunction with an instruction execution system, device or equipment. A machine-readable medium can be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium can include, but is not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or equipment, or any suitable combination of the foregoing. A more specific example of a machine-readable storage medium can include an electrical connection based on one or more lines, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing.

[0131] As used in this disclosure, the terms "machine-readable medium" and "computer-readable medium" refer to any computer program product, apparatus, and / or device (e.g., magnetic disk, optical disk, memory, programmable logic device (PLD)) for providing machine instructions and / or data to a programmable processor, including machine-readable media that receives machine instructions as a machine-readable signal. The term "machine-readable signal" refers to any signal for providing machine instructions and / or data to a programmable processor.

[0132] To provide interaction with a user, the systems and techniques described herein can be implemented on a computer having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user can provide input to the computer. Other types of devices can also be used to provide interaction with the user; for example, the feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including acoustic input, voice input, or tactile input).

[0133] The systems and techniques described herein can be implemented in a computing system that includes back-end components (e.g., as a data server), or a computing system that includes middleware components (e.g., an application server), or a computing system that includes front-end components (e.g., a user computer having a graphical user interface or a web browser through which a user can interact with implementations of the systems and techniques described herein), or a computing system that includes any combination of such back-end components, middleware components, or front-end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network (LAN), a wide area network (WAN), and the Internet.

[0134] Computer systems may include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The client and server relationship arises through computer programs running on the respective computers and having a client-server relationship to each other.

[0135] In the above embodiments, they can be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented using software, they can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer programs or instructions. When the computer program or instructions are loaded and executed on a computer, the processes or functions described in the embodiments of the present disclosure are performed in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, a terminal, a user device, or other programmable device. The computer program or instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer program or instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired or wireless means. The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device such as a server or data center that integrates one or more available media. The available medium can be a magnetic medium, such as a floppy disk, a hard disk, or a tape; it can also be an optical medium, such as a digital video disc (DVD); it can also be a semiconductor medium, such as a solid state drive (SSD).

[0136] Although the present disclosure has been described with reference to specific features and embodiments thereof, it will be apparent that various modifications and combinations may be made thereto without departing from the spirit and scope of the present disclosure. Accordingly, this specification and the drawings are merely illustrative of the present disclosure as defined by the appended claims and are deemed to cover any and all modifications, variations, combinations or equivalents within the scope of the present disclosure. Obviously, those skilled in the art may make various modifications and variations to the present disclosure without departing from the spirit and scope of the present disclosure. Thus, the present disclosure is intended to include such modifications and variations if they fall within the scope of the claims of the present disclosure and their equivalents.

Claims

1. A method for determining parameters of a thermoacoustic panel in an annular combustion chamber, characterized in that: The annular combustion chamber comprises: a shell, comprising an outer shell wall and an inner shell wall spaced apart in the axial direction; a combustion chamber disposed between the outer shell wall and the inner shell wall, the combustion chamber comprising an outer wall and an inner wall spaced apart in the axial direction, wherein a thermoacoustic plate is disposed in the axial direction on the outer wall and / or the inner wall. The method is used to determine the setting parameters of the thermoacoustic plate, and the method comprises: Constructing a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the combustion chamber, constructing a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity, and constructing a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber in the combustion chamber; Processing the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function by Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source; If the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the requirements for suppressing low-frequency modal combustion instability, then determining the current setting parameters of the thermoacoustic panel as target setting parameters of the thermoacoustic panel in the annular combustion chamber, wherein the setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel; Before determining the first scattered sound field function and the second scattered sound field function, the method further includes: Determining, based on engineering parameters of the annular combustion chamber and current setting parameters of the thermoacoustic panel, an impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic panel on one side of the combustion chamber, and an impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic panel on one side of the back cavity, includes: According to the linear thermoacoustic theory, the thermoacoustic plate is assumed to be a plate-like structure composed of parallel thin tubes arranged along the radial direction. The sound propagation in the thermoacoustic plate is obtained as follows: Where i represents the imaginary unit, ω represents the frequency, and ρ m represents the average density of the medium in the thermoacoustic plate, S represents the cross-sectional area of the waveguide, φ represents the porosity of the pores in the thermoacoustic plate, represents the velocity disturbance, represents the pressure disturbance, r represents the radial coordinate of the combustion chamber, d represents the ordinary differential operator, γ represents the specific heat ratio, p m represents the average pressure of the medium in the thermoacoustic plate, σ represents the Prandtl number, T m represents the average temperature; f v and f κ is a frequency-dependent function that describes the viscous and thermal coupling effects between the gas and the thermoacoustic plate, respectively. In the sound propagation in the thermoacoustic plate, term A describes the inertial and viscous effects of the gas in the thermoacoustic plate, term B describes the compression effect of the gas in the thermoacoustic plate and the loss caused by heat conduction, and term C represents the additional velocity jump generated by the thermoacoustic energy conversion in the thermoacoustic plate due to the presence of temperature difference; The sound propagation in the thermoacoustic plate is simplified by the quasi-isothermal assumption, thereby obtaining the impedance boundary conditions between the sound velocity and the sound pressure of the thermoacoustic plate on the combustion chamber side, and the impedance boundary conditions between the sound velocity and the sound pressure of the thermoacoustic plate on the back cavity side.

2. The method for determining parameters of a thermoacoustic panel in an annular combustion chamber according to claim 1, wherein: The constructing of a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic panel formed in the combustion chamber comprises: Based on the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the combustion chamber, a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber is constructed. The first scattered acoustic field function is: in, represents the pressure disturbance generated in the combustion chamber by the equivalent monopole source describing the thermoacoustic plate, represents the position information of the equivalent monopole source describing the thermoacoustic plate at the observation point in three-dimensional space, r, θ, x represent the radial coordinate, circumferential coordinate, and axial coordinate in the cylindrical coordinate system respectively, t represents the time of the observation point, N represents the number of equivalent monopole sources describing the thermoacoustic plate, j represents the sequence number of the equivalent monopole source describing the thermoacoustic plate, S represents the area of the equivalent monopole source describing the thermoacoustic plate, represents the gas density in the combustion chamber, represents the normal acoustic particle velocity of the equivalent monopole source of the thermoacoustic plate on the combustion chamber side, is determined based on the impedance boundary condition between the acoustic velocity and the acoustic pressure of the thermoacoustic panel on one side of the combustion chamber, represents the source coordinates of the equivalent monopole source, represents the source time of the equivalent monopole source, and G represents the Green's function.

3. The method for determining parameters of a thermoacoustic panel in an annular combustion chamber according to claim 1, wherein: The method of constructing a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity includes: Based on the impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic plate on one side of the back cavity, a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate in the back cavity is constructed. The second scattered acoustic field function is: in, Describe the pressure disturbance in the back cavity caused by the equivalent monopole source of the thermoacoustic plate, represents the normal acoustic particle velocity of the equivalent monopole source of the thermoacoustic plate on the back cavity side, It is determined based on the impedance boundary condition between the acoustic velocity and the acoustic pressure of the thermal acoustic panel on one side of the back cavity.

4. The method for determining parameters of a thermoacoustic panel in an annular combustion chamber according to claim 1, wherein: The constructing of a third scattered acoustic field function formed in the combustion chamber by an equivalent monopole source describing the heat source in the combustion chamber comprises: Determine the flame combustion function that describes the unsteady heat release characteristics of the flame of the heat source; Based on the flame combustion function, a third scattered acoustic field function is constructed to describe the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber. The third scattered acoustic field function is: in, represents the equivalent monopole source describing the heat source in the combustion chamber, and the scattered sound pressure formed in the combustion chamber, N f represents the number of equivalent monopole sources describing the heat source in the combustion chamber, q represents the sequence number of the equivalent monopole source describing the heat source, S ′ represents the area of the equivalent monopole source describing the heat source, Represents the acoustic particle velocity at the interface between the heat source outlet and the combustion chamber inlet.

5. The method for determining parameters of a thermoacoustic panel in an annular combustion chamber according to claim 1, wherein: The method further comprises: If the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel do not meet the low-frequency modal combustion instability suppression requirement, then determining updated setting parameters of the thermoacoustic panel within a parameter threshold interval of the thermoacoustic panel, and obtaining an updated characteristic equation of the coupling between the thermoacoustic panel and the heat source based on the updated setting parameters of the thermoacoustic panel; The above process of obtaining the updated characteristic equation is repeated until the characteristic solution of the updated characteristic equation indicates the updated setting parameters of the thermoacoustic panel, which meets the requirements for suppressing low-frequency modal combustion instabilities, and the updated setting parameters of the thermoacoustic panel are determined as the target setting parameters of the thermoacoustic panel in the annular combustion chamber.

6. A device for determining parameters of a thermoacoustic panel in an annular combustion chamber, characterized in that: The annular combustion chamber comprises: a shell, comprising an outer shell wall and an inner shell wall spaced apart in the axial direction; a combustion chamber disposed between the outer shell wall and the inner shell wall, the combustion chamber comprising an outer wall and an inner wall spaced apart in the axial direction, wherein the outer wall and / or the inner wall are provided with thermoacoustic plates in the axial direction, and the device is used to determine the setting parameters of the thermoacoustic plates provided in the axial direction on the outer wall and / or the inner wall of the combustion chamber, the device comprising: a construction module configured to construct a first scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the combustion chamber, construct a second scattered acoustic field function describing the equivalent monopole source of the thermoacoustic plate formed in the back cavity, and construct a third scattered acoustic field function describing the equivalent monopole source of the heat source in the combustion chamber formed in the combustion chamber; a data processing module configured to process the first scattered acoustic field function, the second scattered acoustic field function, and the third scattered acoustic field function through Laplace transform to obtain a characteristic equation of the coupling between the thermoacoustic panel and the heat source; a determination module configured to determine the current setting parameters of the thermoacoustic panel as target setting parameters of the thermoacoustic panel in the annular combustion chamber if the characteristic solution of the characteristic equation indicates that the current setting parameters of the thermoacoustic panel meet the low-frequency modal combustion instability suppression requirements, wherein the setting parameters include setting position parameters and / or size parameters of the thermoacoustic panel; Before determining the first scattered acoustic field function and the second scattered acoustic field function, the method further includes: determining, based on engineering parameters of the annular combustion chamber and current setting parameters of the thermoacoustic panel, an impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic panel on one side of the combustion chamber, and an impedance boundary condition between the acoustic velocity and acoustic pressure of the thermoacoustic panel on the back cavity side, including: According to the linear thermoacoustic theory, the thermoacoustic plate is assumed to be a plate-like structure composed of parallel thin tubes arranged along the radial direction. The sound propagation in the thermoacoustic plate is obtained as follows: Where i represents the imaginary unit, ω represents the frequency, and ρ m represents the average density of the medium in the thermoacoustic plate, S represents the cross-sectional area of the waveguide, φ represents the porosity of the pores in the thermoacoustic plate, represents the velocity disturbance, represents the pressure disturbance, r represents the radial coordinate of the combustion chamber, d represents the ordinary differential operator, γ represents the specific heat ratio, p m represents the average pressure of the medium in the thermoacoustic plate, σ represents the Prandtl number, T m represents the average temperature; f υ and f κ is a frequency-dependent function that describes the viscous and thermal coupling effects between the gas and the thermoacoustic plate, respectively. In the sound propagation in the thermoacoustic plate, term A describes the inertial and viscous effects of the gas in the thermoacoustic plate, term B describes the compression effect of the gas in the thermoacoustic plate and the loss caused by heat conduction, and term C represents the additional velocity jump generated by the thermoacoustic energy conversion in the thermoacoustic plate due to the presence of temperature difference; The sound propagation in the thermoacoustic plate is simplified by the quasi-isothermal assumption, thereby obtaining the impedance boundary conditions between the sound velocity and the sound pressure of the thermoacoustic plate on the combustion chamber side, and the impedance boundary conditions between the sound velocity and the sound pressure of the thermoacoustic plate on the back cavity side.

7. An electronic device, characterized in that: include: processor; as well as, Memory for storing programs, The program includes instructions, which, when executed by the processor, cause the processor to perform the method according to any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium storing computer instructions, characterized in that: The computer instructions are used to cause the computer to execute the method according to any one of claims 1 to 5.

9. An annular combustion chamber, characterized in that: The annular combustion chamber comprises: The housing comprises an outer shell wall and an inner shell wall spaced apart in the axial direction; A combustion chamber is arranged between the outer shell wall and the inner shell wall, and the combustion chamber includes an outer wall and an inner wall arranged at intervals along the axial direction, wherein the outer wall and / or the inner wall are axially provided with a thermoacoustic plate, and the setting parameters of the thermoacoustic plate are determined based on the method described in any one of claims 1 to 5.

Images