Combustion simulation method, device and equipment for combustion chamber and medium

Abstract

The invention relates to the technical field of engine combustion chamber combustion simulation, and particularly provides a combustion simulation method, device, equipment and medium for a combustion chamber, and the method comprises the steps: obtaining a current chemical reaction parameter and a current flow characteristic parameter in the combustion chamber in a current simulation iteration process; based on the current chemical reaction parameters or the current flow characteristic parameters and a pre-constructed model parameter equation, current model parameter values in the current simulation iteration process are determined, flame factors in a thickening factor equation are used for distinguishing combustion modes, and turbulence diffusion coefficients in a diffusion coefficient equation are used for representing turbulence flow characteristics; substituting the current model parameter value into the dynamic thickened flame surface model to obtain a target dynamic thickened flame surface model to be solved; and solving the target dynamic thickened flame surface model to obtain a combustion simulation result in the current simulation iteration process. The accuracy and reliability of the obtained combustion chamber simulation calculation result can be improved.

Description

Technical Field

[0001] The present disclosure relates to the technical field of combustion simulation of engine combustion chambers, and in particular to a combustion simulation method, device, equipment and medium for combustion chambers. Background Art

[0002] The combustion chamber is one of the core components of aircraft engines and gas turbines for achieving energy conversion. It usually works under extreme conditions of high temperature and high pressure. When designing the combustion chamber, the combustion process of the combustion chamber can be simulated and calculated to facilitate the reasonable design of the combustion chamber structure based on the simulation results.

[0003] In the related art, the combustion simulation of the combustion chamber is usually implemented based on the combustion model, among which the dynamic thickened flame surface model based on the thickened flame surface idea is widely used due to its small amount of calculation.

[0004] However, the dynamic thickening flame surface model provided in the related art usually does not consider the influence of turbulent flow and different flame combustion modes on the combustion simulation, resulting in low accuracy of the obtained combustion chamber simulation calculation results. Summary of the invention

[0005] The present disclosure is proposed in view of the above problems, and provides a combustion simulation method, device, equipment and medium for a combustion chamber, which can improve the accuracy and reliability of the obtained combustion chamber simulation calculation results.

[0006] According to one aspect of the present disclosure, a combustion simulation method for a combustion chamber is provided, comprising:

[0007] Obtain the current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during the current simulation iteration;

[0008] Based on the current chemical reaction parameters or the current flow characteristic parameters, and the pre-constructed model parameter equation, determine the current model parameter value in the current simulation iteration process, wherein the model parameter equation includes a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation, the equation parameters in the thickening factor equation include a flame factor for distinguishing combustion modes, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics;

[0009] Substituting the current model parameter value into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved;

[0010] The target dynamic thickening flame surface model is solved to obtain the combustion simulation result in the current simulation iteration process.

[0011] According to another aspect of the present disclosure, there is provided a combustion simulation device for a combustion chamber, comprising:

[0012] An acquisition module is configured to acquire current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during a current simulation iteration;

[0013] A determination module is configured to determine the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and a pre-constructed model parameter equation, wherein the model parameter equation includes a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation, the equation parameters in the thickening factor equation include a flame factor for distinguishing a combustion mode, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics;

[0014] A processing module is configured to substitute the current model parameter value into the dynamic thickening flame surface model to obtain a target dynamic thickening flame surface model to be solved;

[0015] The solving module is configured to solve the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process.

[0016] According to another aspect of the present disclosure, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the combustion simulation method of the combustion chamber is implemented.

[0017] According to another aspect of the present disclosure, a computer program product is provided, including a computer program, which implements the above-mentioned combustion simulation method of the combustion chamber when executed by a processor.

[0018] The present invention provides a combustion simulation method, device, equipment and medium for a combustion chamber. During the combustion simulation process of the combustion chamber, the current model parameter values ​​of the dynamic thickening flame surface model in the current simulation iteration process can be determined through pre-constructed detection function equations, flame index equations, thickening factor equations and diffusion coefficient equations, and the current model parameter values ​​are substituted into the dynamic thickening flame surface model for solving to obtain the combustion simulation results in the current simulation iteration process. Since the equation parameters in the thickening factor equation include the flame factor for distinguishing the combustion mode, and the equation parameters in the diffusion coefficient equation include the turbulent diffusion coefficient for characterizing the turbulent flow characteristics, the combustion simulation of the combustion chamber under the influence of turbulent flow and different flame combustion modes is realized, and the accuracy and reliability of the obtained combustion chamber simulation calculation results are improved.

[0019] It is to be understood that both the foregoing general description and the following detailed description are exemplary, and are intended to provide further explanation of the technology as claimed. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The above and other purposes, features and advantages of the present disclosure will become more apparent by describing the embodiments of the present disclosure in more detail in conjunction with the accompanying drawings. The accompanying drawings are used to provide a further understanding of the embodiments of the present disclosure and constitute a part of the specification. Together with the embodiments of the present disclosure, they are used to explain the present disclosure and do not constitute a limitation of the present disclosure. In the accompanying drawings, the same reference numerals generally represent the same components or steps.

[0021] Figure 1 is a flow chart illustrating a combustion simulation method for applying a combustion chamber according to an embodiment of the present disclosure.

[0022] Figure 2 It is a schematic diagram of a simulation result illustrating an embodiment of the present disclosure.

[0023] Figure 3 It is another schematic diagram of simulation results illustrating an embodiment of the present disclosure.

[0024] Figure 4 It is another schematic diagram of simulation results illustrating an embodiment of the present disclosure.

[0025] Figure 5 It is a schematic diagram of another simulation result of an embodiment of the present disclosure.

[0026] Figure 6 is a thickening factor distribution cloud diagram illustrating an embodiment of the present disclosure.

[0027] Figure 7 It is a cloud diagram of efficiency function distribution of an embodiment of the present disclosure.

[0028] Figure 8 is a methane component distribution cloud diagram illustrating an embodiment of the present disclosure.

[0029] Fig. 9 Graph showing an oxygen component distribution cloud diagram according to an embodiment of the present disclosure.

[0030] Fig.10 It is a temperature distribution cloud diagram illustrating an embodiment of the present disclosure.

[0031] Fig.11 is a block diagram illustrating a combustion simulation apparatus to which a combustion chamber according to an embodiment of the present disclosure is applied.

[0032] Fig.12 is a schematic diagram illustrating a computer program product according to an embodiment of the present disclosure.

[0033] Fig.13 is a hardware block diagram illustrating an electronic device according to an embodiment of the present disclosure. DETAILED DESCRIPTION

[0034] In order to make the purpose, technical solution and advantages of the present disclosure more obvious, the exemplary embodiments according to the present disclosure will be described in detail with reference to the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present disclosure, rather than all the embodiments of the present disclosure, and it should be understood that the present disclosure is not limited to the exemplary embodiments described here.

[0035] In order to solve the above problems, the embodiment of the present disclosure provides a combustion chamber combustion simulation method, which can be applied to a terminal device, which can be a computer, a notebook or a tablet computer, etc. Figure 1 As shown, the combustion simulation method of the combustion chamber includes:

[0036] Step S101, obtaining current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during the current simulation iteration;

[0037] Step S102, determining the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation;

[0038] Among them, the model parameter equations include a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation. The equation parameters in the thickening factor equation include a flame factor for distinguishing combustion modes, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics.

[0039] Step S103, substituting the current model parameter values ​​into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved;

[0040] Step S104, solving the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process.

[0041] In summary, the combustion simulation method for the combustion chamber provided in the example of the present disclosure can, during the combustion simulation process of the combustion chamber, determine the current model parameter values ​​of the dynamic thickening flame surface model in the current simulation iteration process through the pre-constructed detection function equation, flame index equation, thickening factor equation and diffusion coefficient equation, and substitute the current model parameter values ​​into the dynamic thickening flame surface model for solving to obtain the combustion simulation result in the current simulation iteration process. Since the equation parameters in the thickening factor equation include the flame factor for distinguishing the combustion mode, and the equation parameters in the diffusion coefficient equation include the turbulent diffusion coefficient for characterizing the turbulent flow characteristics, the combustion simulation of the combustion chamber under the influence of turbulent flow and different flame combustion modes is realized, and the accuracy and reliability of the obtained combustion chamber simulation calculation results are improved.

[0042] The following Figure 1The specific implementation methods of each step in the embodiment shown are described in detail:

[0043] In step S101, the terminal device obtains current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during the current simulation iteration.

[0044] In the disclosed embodiment, chemical reaction parameters are used to characterize the chemical reaction characteristics of the combustion chamber. The chemical reaction parameters may include: chemical reaction activation energy, reaction index of fuel and reaction index of oxidant. The chemical reaction parameters are determined by a chemical reaction simulation module participating in the combustion simulation of the combustion chamber; flow characteristic parameters are used to characterize the flow characteristics of the combustion chamber. The flow characteristic parameters may include efficiency function values. The flow characteristic parameters are determined by a flow simulation module participating in the combustion simulation of the combustion chamber; wherein, the efficiency function value is determined by the flow simulation module based on the efficiency function proposed by Colin. The efficiency function value is used to characterize the intensity of the wrinkling effect of small-scale turbulence on the flame surface.

[0045] In an optional embodiment, the process by which the terminal device obtains the current chemical reaction parameters and the current flow characteristic parameters in the combustion chamber during the current simulation iteration may include: obtaining the current chemical reaction parameters during the current simulation iteration in response to receiving the current chemical reaction parameters during the current simulation iteration sent by the chemical reaction simulation module; and obtaining the current flow characteristic parameters during the current simulation iteration in response to receiving the current flow characteristic parameters during the current simulation iteration sent by the flow simulation module.

[0046] It should be noted that the chemical reaction simulation module and the flow simulation module can be deployed in a terminal device, and the chemical reaction simulation module and the flow simulation module can also be deployed in other terminal devices, which is not limited in the embodiments of the present disclosure.

[0047] In step S102, the terminal device determines the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation.

[0048] In the embodiment of the present disclosure, the model parameter equations include a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation. The equation parameters in the thickening factor equation include a flame factor for distinguishing combustion modes, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics. It can be understood that the current model parameters include the current detection function value, the current flame index value, the current thickening factor value and the current diffusion coefficient value.

[0049] In an optional embodiment, the process in which the terminal device determines the current model parameter value in the current simulation iteration process based on the current chemical reaction parameters or the current flow characteristic parameters and the pre-constructed model parameter equation includes: inputting the reaction index of the current fuel, the reaction index of the current oxidant, the current reaction activation energy in the current chemical reaction parameters, the mass fraction of the fuel component, the mass fraction of the oxidant component and the temperature in the combustion chamber attribute field information, and the pre-configured first parameter and gas state constant into the detection function equation to obtain the current detection function value in the current simulation iteration process. The reaction activation energy in each simulation iteration process can be considered to determine the detection function value in each simulation iteration process to ensure that in the combustion simulation process of the combustion chamber based on the dynamic thickening flame surface model, thickening is performed only in the combustion reaction area, thereby improving the accuracy of the combustion simulation of the combustion chamber.

[0050] Among them, the detection function equation is:

[0051]

[0052] In formula 1, ω is the detection function value, Y F is the mass fraction of the fuel component, Y O is the mass fraction of the oxidant component, v F is the reactivity index of the fuel, v O is the reaction index of the oxidant, E a is the reaction activation energy, T is the temperature, R is the gas state constant, and Γ is the first parameter.

[0053] It should be noted that, in the embodiment of the present disclosure, the combustion chamber property field information may include velocity field information, mixture fraction field information, temperature field information, component field information, etc.; among which, the component field information includes the mass fraction of the fuel component and the mass fraction of the oxidizer component, and the temperature field information includes temperature.

[0054] It can be understood that the gas state constant and the first parameter can be determined based on actual needs, and the embodiments of the present disclosure are not limited to this, wherein the first parameter can be less than 1.

[0055] In an optional embodiment, the process in which the terminal device determines the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation includes: obtaining a pre-configured second parameter; then, inputting the second parameter and the current detection function value into the flame index equation to obtain the current flame index value in the current simulation iteration process; wherein the second parameter is used to control the size of the transition zone between the thickened area and the non-thickened area, and the second parameter can be determined based on actual needs, which is not limited in the embodiment of the present disclosure. For example, the second parameter can be 10. The flame index value in each simulation iteration process can be determined based on the detection function value in each simulation iteration process to ensure the accuracy of the flame area determined in the combustion simulation process of the combustion chamber based on the dynamic thickening flame surface model, so as to improve the accuracy of the combustion simulation of the combustion chamber.

[0056] The flame index equation is:

[0057]

[0058] In formula 2, Ω is the flame index value, C F is the second parameter.

[0059] In an optional embodiment, the process of determining the current model parameter value in the current simulation iteration process by the terminal device based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation includes: determining the current maximum value of the thickening factor according to the current flame surface thickness in the current chemical reaction parameter and the pre-configured number of grids and grid size; further, if the flame factor is greater than or equal to the preset parameter, the current maximum value of the thickening factor and the current flame index value are input into the thickening factor equation to obtain the current thickening factor value, wherein the flame factor is based on the combustion chamber attribute field information, The product of the gradient value of the mean mass fraction of the fuel component and the gradient value of the mean mass fraction of the oxidant component, wherein the preset parameter is 0, when the flame factor is greater than or equal to 0, the combustion belongs to the premixed combustion mode, and when the flame factor is less than 0, the combustion belongs to the non-premixed combustion mode; the flame factor determined according to the gradient value of the mean mass fraction of the fuel component and the oxidant during the combustion simulation and the comparison result with the preset parameters can be used to determine the combustion mode of the combustion chamber, so as to ensure that the combustion mode of the combustion chamber is taken into account during the combustion simulation of the combustion chamber based on the dynamic thickening flame surface model, thereby improving the simulation reliability and accuracy of the combustion chamber.

[0060] The thickening factor equation is:

[0061]

[0062] In formula 3, F is the thickening factor value, F max is the maximum value of the thickening factor, n is the number of grids, Δ x is the grid size, δ l is the flame surface thickness, FI is the flame factor, is the gradient value of the mean mass fraction of the fuel component, is the gradient value of the mean mass fraction of the oxidant component.

[0063] It should be noted that if the flame factor is less than the preset parameter, the current thickening factor value can be determined to be 1; the component field information also includes the gradient value of the mass fraction mean of the fuel component and the gradient value of the mass fraction mean of the oxidizer component.

[0064] In an optional embodiment, the process in which the terminal device determines the current model parameter value in the current simulation iteration process based on the current chemical reaction parameters or the current flow characteristic parameters, and the pre-constructed model parameter equation includes: inputting the current flame index value, the current thickening factor value, the current efficiency function value in the current flow characteristic parameters, and the turbulent diffusion coefficient and laminar diffusion coefficient in the combustion chamber property field information into the diffusion coefficient equation to obtain the current diffusion coefficient value; the diffusion coefficient value in each simulation iteration process can be determined based on the turbulent diffusion coefficient to ensure that in the combustion simulation process of the combustion chamber based on the dynamic thickening flame surface model, the turbulence effect is taken into account to improve the accuracy of the combustion simulation of the combustion chamber.

[0065] The diffusion coefficient equation is:

[0066] D eff =EFD+(1-Ω)D t ; (Formula 4)

[0067] In Formula 4, D eff is the diffusion coefficient value, E is the efficiency function value, D is the laminar diffusion coefficient, D t is the turbulent diffusion coefficient.

[0068] In step S103, the terminal device substitutes the current model parameter value into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved.

[0069] In an embodiment of the present disclosure, the dynamic thickening flame surface model may include a transport equation for mass fraction and a transport equation for enthalpy.

[0070] The transport equation for mass fraction is:

[0071]

[0072] In formula 5, is the mean value of the mass fraction of the ith component, is the mean velocity of the flow field, is the mean flow field density, is the source phase of the ith component, is the mean value of the diffusion coefficient;

[0073] The transport equation for enthalpy is:

[0074]

[0075] In Formula 6, is the mean enthalpy.

[0076] It should be noted that in the embodiments of the present disclosure, the mean of the parameters in the transport equation is the average of the parameter values ​​determined in the current simulation iteration and the mean of the parameters determined in the previous simulation iteration; for example, the mean of the diffusion coefficient value is the current diffusion coefficient value and the average of the mean of the diffusion coefficient determined in the previous simulation iteration.

[0077] In an optional implementation, the process in which the terminal device substitutes the current model parameter value into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved may include: substituting the current thickening factor value and the current diffusion coefficient value in the current model parameter value, as well as the current efficiency function value in the current flow characteristic parameter, into the mass fraction transport equation to obtain the mass fraction transport equation to be solved; and substituting the current diffusion coefficient value into the enthalpy transport equation to obtain the enthalpy transport equation to be solved; and at the same time, combining the mass fraction transport equation to be solved and the enthalpy transport equation to be solved to obtain the target dynamic thickening flame surface model to be solved. The model parameter values ​​determined in each simulation iteration process may be substituted into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved in each simulation iteration process, so as to solve the target dynamic thickening flame surface model and obtain more accurate combustion simulation results in each simulation iteration process.

[0078] In step S104, the terminal device solves the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process.

[0079] In the disclosed embodiment, the combustion simulation results may include distribution information of each component, enthalpy value, temperature field, and combustion products.

[0080] In an optional implementation, the process in which the terminal device solves the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process includes: solving the target dynamic thickening flame surface model by a partial differential solver to obtain the combustion simulation result in the current simulation iteration process. The target dynamic thickening flame surface model can be solved by a partial differential solver to improve the efficiency of solving the combustion simulation result.

[0081] It should be noted that in the embodiment of the present disclosure, the combustion simulation of the combustion chamber needs to be implemented after multiple iterations of the coordinated chemical reaction simulation module and the flow simulation module, wherein the iteration stop condition can be determined based on actual needs, and the embodiment of the present disclosure does not limit this.

[0082] For example, Figures 2 to 5 The figure shows the comparison results of the temperature at different radial positions of the combustion chamber in the calculation process of the unrestricted methane air premixed swirl flame by the three methods of the disclosed embodiment, the experiment and the literature calculation, wherein: Figure 2 Schematic diagram of temperature distribution determined by three methods at 10mm in the combustion chamber. Figure 3 The temperature distribution diagram at 20 mm in the combustion chamber is determined by three methods. Figure 4 The temperature distribution diagram is determined by three methods at 30mm in the combustion chamber. Figure 5 Schematic diagram of temperature distribution determined by three methods at 60mm in the combustion chamber.

[0083] It should be noted that Figures 2 to 5 In the figure, the curve composed of circles is the experimental result, the dotted line is the result calculated in the literature, and the solid line is the result determined in the embodiment of the present disclosure. Figures 2 to 5 It can be seen that the temperature distribution calculated by the combustion simulation scheme of the combustion chamber provided by the embodiment of the present disclosure is highly consistent with the experimental results at 10mm and 20mm upstream of the combustion chamber, and is also very consistent at 30mm and 60mm downstream of the combustion chamber, and the distribution within the range from the center axis to half the radius is also very consistent. Compared with the calculation results of published literature, it has higher accuracy.

[0084] at the same time, Figures 6 to 10 The distribution cloud diagrams of the thickening factor, efficiency function, methane component, oxygen component and temperature distribution along the axial section in the middle of the combustion chamber in the calculation process of the unrestricted methane-air premixed swirl flame according to the embodiment of the present disclosure are shown in sequence. Figure 6 The thickening factor distribution cloud diagram and Figure 7 It can be seen from the distribution cloud diagram of the efficiency function that the effective position of the thickening factor is consistent with the effective position of the efficiency function, both of which are near the swirl outlet where the components react to form the flame surface, proving the accuracy of the effective position of the thickening factor determined by the embodiment of the present disclosure; Figure 8The methane component distribution cloud diagram and Fig. 9 The oxygen component distribution cloud diagram shown is Fig.10 Temperature distribution cloud shown Figure 1 This proves the accuracy of the flame area determined by the embodiment of the present disclosure.

[0085] An exemplary embodiment of the present disclosure provides a combustion simulation device for a combustion chamber, which may be a chip of a terminal device. Fig.11 FIG. 1 shows a schematic block diagram of functional modules of a combustion simulation device for a combustion chamber according to an exemplary embodiment of the present disclosure. Fig.11 As shown, the combustion simulation device 1100 of the combustion chamber includes:

[0086] An acquisition module 1101 is configured to acquire current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during a current simulation iteration;

[0087] The determination module 1102 is configured to determine the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation, wherein the model parameter equation includes a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation, the equation parameters in the thickening factor equation include a flame factor for distinguishing a combustion mode, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics;

[0088] The processing module 1103 is configured to substitute the current model parameter value into the dynamic thickening flame surface model to obtain a target dynamic thickening flame surface model to be solved;

[0089] The solving module 1104 is configured to solve the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process.

[0090] Optionally, the determining module 1102 is configured to:

[0091] The reaction index of the current fuel, the reaction index of the current oxidant, the current reaction activation energy in the current chemical reaction parameters, the mass fraction of the fuel component, the mass fraction and temperature of the oxidant component in the combustion chamber attribute field information, and the pre-configured first parameter and gas state constant are input into the detection function equation to obtain the current detection function value in the current simulation iteration process. The detection function equation is:

[0092]

[0093] Among them, ω is the detection function value, Y F is the mass fraction of the fuel component, YO is the mass fraction of the oxidant component, v F is the reactivity index of the fuel, v O is the reaction index of the oxidant, E a is the reaction activation energy, T is the temperature, R is the gas state constant, and Γ is the first parameter.

[0094] Optionally, the determining module 1102 is further configured to:

[0095] Acquire a preconfigured second parameter, wherein the second parameter is used to control the size of a transition zone between a thickened area and a non-thickened area;

[0096] The second parameter and the current detection function value are input into the flame index equation to obtain the current flame index value in the current simulation iteration process. The flame index equation is:

[0097]

[0098] Where Ω is the flame index value, C F is the second parameter.

[0099] Optionally, the determining module 1102 is configured to:

[0100] Determine the maximum value of the current thickening factor according to the current flame surface thickness in the current chemical reaction parameters, and the pre-configured number of grids and grid size;

[0101] If the flame factor is greater than or equal to the preset parameter, the current maximum value of the thickening factor and the current flame index value are input into the thickening factor equation to obtain the current thickening factor value, wherein the flame factor is the product of the gradient value of the mass fraction mean of the fuel component and the gradient value of the mass fraction mean of the oxidant component in the combustion chamber attribute field information, and the thickening factor equation is:

[0102]

[0103] Among them, F is the thickening factor value, F max is the maximum value of the thickening factor, n is the number of grids, Δ x is the grid size, δ l is the flame surface thickness, FI is the flame factor, is the gradient value of the mean mass fraction of the fuel component, is the gradient value of the mean mass fraction of the oxidant component.

[0104] Optionally, the determining module 1102 is configured to:

[0105] The current flame index value, the current thickening factor value, the current efficiency function value in the current flow characteristic parameters, and the turbulent diffusion coefficient and laminar diffusion coefficient in the combustion chamber attribute field information are input into the diffusion coefficient equation to obtain the current diffusion coefficient value. The diffusion coefficient equation is:

[0106] D eff =EFD+(1-Ω)D t ;

[0107] Among them, D eff is the diffusion coefficient value, E is the efficiency function value, D is the laminar diffusion coefficient, D t is the turbulent diffusion coefficient.

[0108] Optionally, the processing module 1103 is configured to:

[0109] Substituting the current thickening factor value and the current diffusion coefficient value in the current model parameter value, and the current efficiency function value in the current flow characteristic parameter, into the transport equation of the mass fraction to obtain the transport equation of the mass fraction to be solved;

[0110] Substituting the current diffusion coefficient value into the enthalpy transport equation to obtain the enthalpy transport equation to be solved;

[0111] The transport equation of the mass fraction to be solved and the transport equation of enthalpy to be solved are combined to obtain the target dynamic thickening flame surface model to be solved. The transport equation of the mass fraction is:

[0112]

[0113] in, is the mean value of the mass fraction of the ith component, is the mean velocity of the flow field, is the mean flow field density, is the source phase of the ith component, is the mean value of the diffusion coefficient;

[0114] The transport equation for enthalpy is:

[0115]

[0116] in, is the mean enthalpy.

[0117] Optionally, the solution module 1104 is configured to:

[0118] The target dynamic thickening flame surface model is solved by a partial differential solver to obtain a combustion simulation result in the current simulation iteration process.

[0119] The exemplary embodiment of the present disclosure also provides an electronic device, comprising: at least one processor; and a memory connected to the at least one processor in communication. The memory stores a computer program that can be executed by the at least one processor, and the computer program is used to cause the electronic device to perform the method according to the embodiment of the present disclosure when executed by the at least one processor.

[0120] The exemplary embodiments of the present disclosure also 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 the method according to the embodiments of the present disclosure.

[0121] like Fig.12 As shown, the exemplary embodiment of the present disclosure further provides a computer program product 1200, including a computer program 1201, wherein when the computer program is executed by a processor of a computer, it is used to enable the computer to perform the method according to the embodiment of the present disclosure.

[0122] refer to Fig.13 , a block diagram of an electronic device 1300 that can be used as a terminal device of 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, workbenches, 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.

[0123] like Fig.13 As shown, the electronic device 1300 includes a computing unit 1301, which can perform various appropriate actions and processes according to a computer program stored in a read-only memory (ROM) 1302 or a computer program loaded from a storage unit 1308 into a random access memory (RAM) 1303. In the RAM 1303, various programs and data required for the operation of the electronic device 1300 can also be stored. The computing unit 1301, the ROM 1302, and the RAM 1303 are connected to each other via a bus 1304. An input / output (I / O) interface 1305 is also connected to the bus 1304.

[0124] A plurality of components in the electronic device 1300 are connected to the I / O interface 1305, including: an input unit 1306, an output unit 1307, a storage unit 1308, and a communication unit 1309. The input unit 1306 may be any type of device capable of inputting information to the electronic device 1300, and the input unit 1306 may receive input digital or character information, and generate key signal inputs related to user settings and / or function control of the electronic device. The output unit 1307 may 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. The storage unit 1308 may include but is not limited to a disk, an optical disk. The communication unit 1309 allows the electronic device 1300 to exchange information / data with other devices through a computer network 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 a chipset, such as a Bluetooth™ device, a WiFi device, a WiMax device, a cellular communication device, and / or the like.

[0125] The computing unit 1301 may be a variety of general and / or special processing components with processing and computing capabilities. Some examples of the computing unit 1301 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 running machine learning model algorithms, digital signal processors (DSPs), and any appropriate processors, controllers, microcontrollers, etc. The computing unit 1301 performs the various methods and processes described above. For example, in some embodiments, the method of the exemplary 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 1308. In some embodiments, part or all of the computer program may be loaded and / or installed on the electronic device 1300 via the ROM 1302 and / or the communication unit 1309. In some embodiments, the computing unit 1301 may be configured to perform the method of the exemplary embodiment of the present disclosure by any other appropriate means (e.g., by means of firmware).

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

[0127] In the context of the present disclosure, a machine-readable medium may be a tangible medium that may contain or store a program for use by or in conjunction with an instruction execution system, device, or equipment. A machine-readable medium may be a machine-readable signal medium or a machine-readable storage medium. A machine-readable medium may 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 may 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.

[0128] 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., disk, optical disk, memory, programmable logic device (PLD)) for providing machine instructions and / or data to a programmable processor, including a machine-readable medium 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.

[0129] 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).

[0130] The systems and techniques described herein may 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 with 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 may 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.

[0131] A computer system may include clients and servers. Clients and servers are generally remote from each other and usually interact through a communication network. The relationship of client and server is generated by computer programs running on respective computers and having a client-server relationship to each other.

[0132] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented by software, it 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 instruction is loaded and executed on a computer, the process or function described in the embodiment of the present disclosure is executed in whole or in part. The computer may 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 instruction may 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 instruction may be transmitted from one website site, computer, server or data center to another website site, computer, server or data center by wired or wireless means. The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server, data center, etc. that integrates one or more available media. The available medium may be a magnetic medium, for example, a floppy disk, a hard disk, a tape; it may also be an optical medium, for example, a digital video disc (DVD); it may also be a semiconductor medium, for example, a solid state drive (SSD).

[0133] Although the present disclosure has been described in conjunction with specific features and embodiments thereof, it is 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 exemplary illustrations of the present disclosure as defined by the appended claims, and are deemed to have covered 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, if these modifications and variations of the present disclosure fall within the scope of the claims of the present disclosure and their equivalents, the present disclosure is also intended to include these modifications and variations.

Claims

1. A combustion simulation method for a combustion chamber, characterized in that: include: Obtain the current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during the current simulation iteration; Based on the current chemical reaction parameters or the current flow characteristic parameters, and the pre-constructed model parameter equation, determine the current model parameter value in the current simulation iteration process, wherein the model parameter equation includes a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation, the equation parameters in the thickening factor equation include a flame factor for distinguishing combustion modes, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics; Substituting the current model parameter values ​​into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved; The target dynamic thickening flame surface model is solved to obtain the combustion simulation result in the current simulation iteration process.

2. The combustion simulation method of a combustion chamber according to claim 1, characterized in that: The determining of the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation includes: The reaction index of the current fuel, the reaction index of the current oxidant, the current reaction activation energy in the current chemical reaction parameters, the mass fraction of the fuel component, the mass fraction and temperature of the oxidant component in the combustion chamber attribute field information, and the pre-configured first parameter and gas state constant are input into the detection function equation to obtain the current detection function value in the current simulation iteration process. The detection function equation is: Among them, ω is the detection function value, Y F is the mass fraction of the fuel component, Y O is the mass fraction of the oxidant component, v F is the reactivity index of the fuel, v O is the reaction index of the oxidant, E a is the reaction activation energy, T is the temperature, R is the gas state constant, and Γ is the first parameter.

3. The combustion simulation method of a combustion chamber according to claim 2, characterized in that: The method further comprises: Acquire a preconfigured second parameter, wherein the second parameter is used to control the size of a transition zone between a thickened area and a non-thickened area; The second parameter and the current detection function value are input into the flame index equation to obtain the current flame index value in the current simulation iteration process. The flame index equation is: Where Ω is the flame index value, C F is the second parameter.

4. The combustion simulation method of a combustion chamber according to claim 3, characterized in that: The determining of the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation includes: Determine the maximum value of the current thickening factor according to the current flame surface thickness in the current chemical reaction parameters, and the pre-configured number of grids and grid size; If the flame factor is greater than or equal to the preset parameter, the current maximum value of the thickening factor and the current flame index value are input into the thickening factor equation to obtain the current thickening factor value, wherein the flame factor is the product of the gradient value of the mass fraction mean of the fuel component and the gradient value of the mass fraction mean of the oxidant component in the combustion chamber attribute field information, and the thickening factor equation is: Among them, F is the thickening factor value, F max is the maximum value of the thickening factor, n is the number of grids, Δ x is the grid size, δ l is the flame surface thickness, FI is the flame factor, is the gradient value of the mean mass fraction of the fuel component, is the gradient value of the mean mass fraction of the oxidant component.

5. The combustion simulation method of a combustion chamber according to claim 4, characterized in that: The determining of the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and the pre-constructed model parameter equation includes: The current flame index value, the current thickening factor value, the current efficiency function value in the current flow characteristic parameters, and the turbulent diffusion coefficient and laminar diffusion coefficient in the combustion chamber attribute field information are input into the diffusion coefficient equation to obtain the current diffusion coefficient value. The diffusion coefficient equation is: D eff =EFD+(1-Ω)D t ; Among them, D eff is the diffusion coefficient value, E is the efficiency function value, D is the laminar diffusion coefficient, D t is the turbulent diffusion coefficient.

6. The combustion simulation method of a combustion chamber according to claim 1, characterized in that: Substituting the current model parameter value into the dynamic thickening flame surface model to obtain the target dynamic thickening flame surface model to be solved includes: Substituting the current thickening factor value and the current diffusion coefficient value in the current model parameter value, and the current efficiency function value in the current flow characteristic parameter, into the transport equation of the mass fraction to obtain the transport equation of the mass fraction to be solved; Substituting the current diffusion coefficient value into the enthalpy transport equation to obtain the enthalpy transport equation to be solved; The transport equation of the mass fraction to be solved and the transport equation of enthalpy to be solved are combined to obtain the target dynamic thickening flame surface model to be solved. The transport equation of the mass fraction is: in, is the mean value of the mass fraction of the ith component, is the mean velocity of the flow field, is the mean flow field density, is the source phase of the ith component, is the mean value of the diffusion coefficient; The transport equation for enthalpy is: in, is the mean enthalpy.

7. The combustion simulation method of a combustion chamber according to claim 1, characterized in that: The step of solving the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process includes: The target dynamic thickening flame surface model is solved by a partial differential solver to obtain a combustion simulation result in the current simulation iteration process.

8. A combustion simulation device for a combustion chamber, characterized in that: include: An acquisition module is configured to acquire current chemical reaction parameters and current flow characteristic parameters in the combustion chamber during a current simulation iteration; A determination module is configured to determine the current model parameter value in the current simulation iteration process based on the current chemical reaction parameter or the current flow characteristic parameter and a pre-constructed model parameter equation, wherein the model parameter equation includes a detection function equation, a flame index equation, a thickening factor equation and a diffusion coefficient equation, the equation parameters in the thickening factor equation include a flame factor for distinguishing a combustion mode, and the equation parameters in the diffusion coefficient equation include a turbulent diffusion coefficient for characterizing turbulent flow characteristics; A processing module is configured to substitute the current model parameter value into the dynamic thickening flame surface model to obtain a target dynamic thickening flame surface model to be solved; The solving module is configured to solve the target dynamic thickening flame surface model to obtain the combustion simulation result in the current simulation iteration process.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the combustion simulation method of the combustion chamber described in any one of claims 1-7 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the combustion simulation method of the combustion chamber described in any one of claims 1-7 is implemented.

Images