A new flame surface turbulent combustion model

Abstract

The present invention belongs to the fields of engine combustion and computational fluid dynamics, and discloses a novel flamelet turbulent combustion model. The process of establishing the model is as follows: Solve the governing equations of the diffusion flamelet turbulent combustion model to obtain the diffusion flamelet database; solve the governing equations of the premixed flamelet turbulent combustion model to obtain the premixed flamelet database; solve the transport equations of velocity, turbulence, mixture fraction, and reaction progress variable to obtain the flow field parameters; obtain the temperature and components at each grid point of the flow field by means of table lookup and interpolation in the diffusion and premixed flamelet databases respectively according to the mixture fraction and reaction progress variable corresponding to each grid point of the flow field; calculate the weighting function from the flow field parameters of each grid point; perform weighted averaging on the obtained temperature and components through the weighting function; update the original temperature and components of the flow field with the weighted average temperature and components; repeat the iteration until the flow field converges. The model gives full play to the advantages of the diffusion and premixed flamelet turbulent combustion models and has engineering applicability.

Description

Technical Field

[0001] The present invention belongs to the fields of engine combustion and computational fluid dynamics, and particularly relates to a new flamelet turbulent combustion model. Background Art

[0002] The turbulent combustion process in an engine is very complex. Turbulence affects the average chemical reaction rate through enhanced mixing, while the chemical reaction heat release process in turn affects turbulent combustion. It is very difficult to achieve accurate numerical simulation of turbulent combustion. Up to now, a variety of turbulent combustion models have been developed. Among them, the flamelet turbulent combustion model decouples turbulent flow from chemical reactions, greatly reducing the computational amount and has been widely applied in engineering practice.

[0003] There are two types of flamelet turbulent combustion models: one is the diffusion flamelet turbulent combustion model, which establishes a diffusion flamelet database by solving the one-dimensional laminar counterflow flame; the other is the premixed flamelet turbulent combustion model, which establishes a premixed flamelet database by solving the one-dimensional laminar premixed flame. During numerical simulation, by querying the flamelet database with local flow field parameters, parameters such as components can be obtained. The two types of models are respectively applicable to diffusion combustion and premixed combustion, while the combustion in an engine is neither purely diffusion combustion nor purely premixed combustion, but partially premixed combustion. To address this problem, it is necessary to improve the existing flamelet turbulent combustion models so that they can improve the simulation ability of partially premixed combustion and enhance their engineering applicability.

[0004] Currently, there is an urgent need to develop a new flamelet turbulent combustion model. Summary of the Invention

[0005] The technical problem to be solved by the present invention is to provide a new flamelet turbulent combustion model to overcome the defects of the prior art.

[0006] The new flamelet turbulent combustion model of the present invention is established as follows:

[0007] Step1. Solve the control equations of the diffusion flamelet turbulent combustion model to obtain a diffusion flamelet database;

[0008] Step2. Solve the control equations of the premixed flamelet turbulent combustion model to obtain a premixed flamelet database;

[0009] Step3. Solve the transport equations of velocity, turbulence, mixture fraction, and reaction progress variable to obtain flow field parameters;

[0010] Step4. Obtain temperature and components by table lookup and interpolation in the diffusion and premixed flamelet databases respectively according to the mixture fraction and reaction progress variable corresponding to each grid point of the flow field;

[0011] Step 5. Calculate the weighting function from the flow field parameters of each grid point;

[0012] Step 6. Perform weighted averaging on the temperature and components obtained in Step 4 through the weighting function;

[0013] Step 7. Update the original temperature and components of the flow field with the weighted average temperature and components;

[0014] Step 8. Repeat Steps 3 to 7 until the flow field converges.

[0015] Furthermore, the governing equation of the diffusion flamelet turbulent combustion model in Step 1 is:

[0016] ;

[0017] where, ρ is the density, T is the temperature, c p is the specific heat at constant pressure, f is the mixture fraction, χ is the scalar dissipation rate, q R is the heat loss; ω i , Y i and h i are the chemical production rate, mass fraction and enthalpy of component i respectively;

[0018] Under different scalar dissipation rate χ conditions, solve Equation (1) to obtain:

[0019] ;

[0020] where, is the temperature or component of the diffusion flamelet turbulent combustion model;

[0021] Define the reaction progress variable C as a linear combination of components to obtain:

[0022] ;

[0023] Perform probability density integration on Equation (3) to obtain:

[0024] ;

[0025] where, is the probability density integration of , and are the mean and variance of the mixture fraction, respectively, and Equation (4) is the diffusion flamelet database.

[0026] Further, the governing equation of the premixed flamelet turbulent combustion model in Step 2 is:

[0027] ;

[0028] where λ is the thermal conductivity, is the density of the unburned gas, is the laminar flame speed, y is the spatial coordinate; and are the specific heat at constant pressure and the diffusion coefficient of component i , respectively;

[0029] Solve Equation (5) under different equivalence ratios Ф to obtain:

[0030] ;

[0031] where is the temperature or component of the premixed flamelet turbulent combustion model;

[0032] The equivalence ratio Ф and the mixture fraction f have the following relationship:

[0033] ;

[0034] where is the mixture fraction under stoichiometric conditions;

[0035] Combine Equation (6) and (7) to obtain:

[0036] ;

[0037] Perform transformation and probability density integration to obtain:

[0038] ;

[0039] where is the probability density integration of , and Equation (9) is the premixed flamelet database.

[0040] Further, the weighted average in Step 6 is based on the Damköhler number Da Define the weighting function F , and combine the diffusion flamelet database and the premixed flamelet database;

[0041] ;

[0042] wherein, Da is the Damköhler number, is the flow time scale, is the chemical reaction time scale, is the laminar flame thickness, ε is the turbulent dissipation rate, k is the turbulent kinetic energy, is the integrated value of temperature or components of the turbulent combustion model of the diffusion flame surface, is the integrated value of temperature or components of the turbulent combustion model of the premixed flame surface. Equation (14) is the flame surface database of the new turbulent combustion model of the flame surface.

[0043] The new turbulent combustion model of the flame surface of the present invention gives full play to the advantages of the turbulent combustion model of the diffusion flame surface and the turbulent combustion model of the premixed flame surface. By combining the diffusion flame surface database and the premixed flame surface database through a weighting function, it effectively solves the problem that the original turbulent combustion model of the flame surface is only applicable to diffusion combustion or premixed combustion alone, and has engineering applicability. Description of the Drawings

[0044] Figure 1 is a schematic diagram of the establishment process of the new turbulent combustion model of the flame surface of the present invention;

[0045] Figure 2 is a schematic diagram of the Cabra partially premixed flame of the embodiment;

[0046] Figure 2 in which, H is the height of flame lift; D is the inner diameter of the fuel nozzle;

[0047] Figure 3a is a comparison curve of the calculated and experimental temperature values of the typical cross-section of the Cabra partially premixed flame obtained by the new turbulent combustion model of the flame surface of the present invention ( x = D);

[0048] Figure 3b is a comparison curve of the calculated and experimental temperature values of the typical cross-section of the Cabra partially premixed flame obtained by the new turbulent combustion model of the flame surface of the present invention ( x = 15D);

[0049] Figure 3c is a comparison curve of the calculated and experimental temperature values of the typical cross-section of the Cabra partially premixed flame obtained by the new turbulent combustion model of the flame surface of the present invention ( x = 30D);

[0050] Figure 3dComparison curve of calculated and experimental temperature values of a typical cross-section of the Cabra partially premixed flame obtained by the new flame surface turbulent combustion model of the present invention ( x = 70D);

[0051] Figures 3a - 3d in x The Specific implementation mode

[0052] The present invention will be described in detail below with reference to the accompanying drawings and embodiments.

[0053] As Figure 1 shown, the establishment process of the new flame surface turbulent combustion model of the present invention is as follows:

[0054] Step1. Solve the control equations of the diffusion flame surface turbulent combustion model to obtain the diffusion flame surface database;

[0055] Step2. Solve the control equations of the premixed flame surface turbulent combustion model to obtain the premixed flame surface database;

[0056] Step3. Solve the transport equations of velocity, turbulence, mixture fraction and reaction progress variable to obtain the flow field parameters;

[0057] Step4. Obtain the temperature and components by looking up tables and interpolation in the diffusion and premixed flame surface databases respectively according to the mixture fraction and reaction progress variable corresponding to each grid point of the flow field;

[0058] Step5. Calculate the weighting function from the flow field parameters of each grid point;

[0059] Step6. Perform weighted averaging on the temperature and components obtained in Step4 through the weighting function;

[0060] Step7. Update the original temperature and components of the flow field with the weighted average temperature and components;

[0061] Step8. Repeat Steps Step3~Step7 until the flow field converges.

[0062] Furthermore, the control equations of the diffusion flame surface turbulent combustion model in Step1 are:

[0063] ;

[0064] Wherein, ρ is the density, T is the temperature, c p is the specific heat at constant pressure, f is the mixture fraction, χ is the scalar dissipation rate,q R is the heat loss; ω i , Y i and h i are the chemical production rate, mass fraction and enthalpy of component i respectively;

[0065] Under different scalar dissipation rate χ conditions, solving equation (1) gives:

[0066] ;

[0067] wherein, is the temperature or component of the diffusion flamelet turbulent combustion model;

[0068] Define the reaction progress variable C as a linear combination of components, and we get:

[0069] ;

[0070] Performing probability density integration on equation (3) gives:

[0071] ;

[0072] wherein, is the probability density integration of , and are the mean and variance of the mixture fraction respectively, and equation (4) is the diffusion flamelet database.

[0073] Furthermore, the control equation of the premixed flamelet turbulent combustion model in Step2 is:

[0074] ;

[0075] wherein, λ is the thermal conductivity, is the density of the unburned gas, is the laminar flame propagation speed, y is the spatial coordinate; and are the isobaric specific heat and diffusion coefficient of component i respectively;

[0076] Under different equivalence ratio Ф conditions, solving equation (5) gives:

[0077] ;

[0078] Among them, is the temperature or component of the premixed flamelet turbulent combustion model;

[0079] Equivalence ratio Ф and mixture fraction f have the following relational formula:

[0080] ;

[0081] Among them, is the mixture fraction under stoichiometric conditions;

[0082] Combining equations (6) and (7), we get:

[0083] ;

[0084] Performing transformation and probability density integration, we get:

[0085] ;

[0086] Among them, is the probability density integration of, and equation (9) is the premixed flamelet database.

[0087] Furthermore, the weighted average in the above-mentioned Step6 is based on the Damköhler number Da Define the weighting function F , and combine the diffusion flamelet database and the premixed flamelet database;

[0088] ;

[0089] Among them, Da is the Damköhler number, is the flow time scale, is the chemical reaction time scale, is the laminar flame thickness, ε is the turbulent dissipation rate, k is the turbulent kinetic energy, is the integral value of the temperature or component of the diffusion flamelet turbulent combustion model, is the integral value of the temperature or component of the premixed flamelet turbulent combustion model, and equation (14) is the flamelet database of the new flamelet turbulent combustion model.

[0090] Example: In this example, the new flamelet turbulent combustion model of the present invention is tested by the Cabra partially premixed flame.

[0091] Figure 2Figure 0 shows a schematic diagram of the Cabra partially premixed flame, which consists of a central fuel jet and a surrounding coflow. The inner diameter D of the nozzle of the fuel jet is 4.57 mm, and it is composed of a mixture of methane and air, with a temperature of 310 K and a jet velocity of 100 m / s; the diameter of the coflow is 210 mm, and it is composed of a mixture of water and air, with a temperature of 1355 K and a coflow velocity of 5.4 m / s.

[0092] Figures 3a - 3d Figure 4 shows the comparison curve of the calculated and experimental temperature values of the typical cross-section of the Cabra partially premixed flame obtained by the new flame surface turbulent combustion model of this embodiment. It can be seen that the calculated temperature values are in good agreement with the experimental values, indicating that the new flame surface turbulent combustion model of the present invention has a high prediction accuracy for partial premixed combustion and has engineering practical value.

[0093] Although the embodiments of the present invention have been disclosed as above, they are not limited to the applications listed in the specification and embodiments. For those skilled in the art, without departing from the principle of the present invention, all the features disclosed in the present invention, or all the steps in the disclosed methods or processes, except for mutually exclusive features and / or steps, can be combined in any way. The present invention is not limited to the specific details and the illustrated examples here.

Claims

1. A novel flame surface turbulent combustion model, characterized in that The establishment process of the described new flame surface turbulent combustion model is as follows: Step1. Solve the governing equations of the diffusion flame surface turbulent combustion model to obtain the diffusion flame surface database; Step2. Solve the governing equations of the premixed flame surface turbulent combustion model to obtain the premixed flame surface database; Step3. Solve the transport equations of velocity, turbulence, mixture fraction, and reaction progress variable to obtain the flow field parameters; Step4. Based on the mixture fraction and reaction progress variable corresponding to each grid point in the flow field, obtain the temperature and components by looking up tables and interpolation in the diffusion and premixed flame surface databases respectively; Step5. Calculate the weighting function from the flow field parameters of each grid point; Step6. Perform weighted averaging on the temperature and components obtained in Step4 through the weighting function; Step7. Update the original temperature and components in the flow field with the weighted average temperature and components; Step8. Repeat Steps Step3~Step7 until the flow field converges.

2. The novel turbulent combustion model of the flame surface according to claim 1, characterized in that The governing equations of the diffusion flame surface turbulent combustion model in Step1 described above are: ； Among them, ρ is the density, T is the temperature, c p is the specific heat at constant pressure, f is the mixture fraction, χ is the scalar dissipation rate, q R is the heat loss; ω i 、 Y i and h i are respectively the chemical production rate, mass fraction and enthalpy of component i ; Under different scalar dissipation rates χ Solve Equation (1) to obtain: ； Among them, is the temperature or component of the turbulent combustion model of the diffusion flame surface; Define the reaction progress variable C as a linear combination of the components, resulting in: ； Perform probability density integration on Equation (3) to obtain: ； wherein, is the probability density integral of and are the mean and variance of the mixture fraction respectively, and Equation (4) is the diffusion flamelet database.

3. The novel flame surface turbulent combustion model according to claim 2, wherein, The governing equations of the premixed flame surface turbulent combustion model in Step2 described above are: ； Among them, λ is the thermal conductivity coefficient, is the density of the unburned gas, is the laminar flame propagation speed, y is the space coordinate; and are respectively the isobaric specific heat and diffusion coefficient of component i ; Under different equivalence ratios Ф Solving Equation (5) under the conditions, we obtain: ； Among them, is the temperature or component of the premixed flamelet turbulent combustion model; Equivalence ratio Ф and mixture fraction f have the following relationship: ； wherein, is the mixture fraction under stoichiometric conditions; Combine Equation (6) and (7) to obtain: ； Perform transformation and probability density integration to obtain: ； Among them, is the probability density integral. Equation (9) is the premixed flamelet database.

4. The novel flame surface turbulent combustion model according to claim 3, wherein The weighted average in Step 6 is based on the Damköhler number Da Define a weighting function F , and combine the diffusion flamelet database and the premixed flamelet database; ； Among them, Da is the Damköhler number, is the flow time scale, is the chemical reaction time scale, is the laminar flame thickness, ε is the turbulent dissipation rate, k is the turbulent kinetic energy, is the integral value of temperature or component of the turbulent combustion model of the diffusion flame surface, is the integral value of temperature or component of the turbulent combustion model of the premixed flame surface. Equation (14) is the flame surface database of the new turbulent combustion model of the flame surface.

Images