Combustion System Design Method Based on Target Heat Release Rate

Abstract

Disclosed is a combustion system design method based on a target heat release rate, which belongs to the technical field of diesel engine combustion chamber design. The method includes: obtaining an ideal heat release rate based on Sabathe-Miller cycle; simulating the ideal heat release rate based on a double-Wiebe function and obtaining the target heat release rate; constructing a mapping relation among the heat release rate, piston geometric parameters and fuel injection parameters, which includes target start of combustion being an function of fuel injection timing and ignition delay, premixed combustion parameters being functions of throat radius, injection pressure and nozzle diameter, and diffusion combustion being a function of piston pit depth; solving target piston geometric parameters and target fuel injection parameters based on the mapping relation; and then designing a combustion system. The method does not depend on experience and multi-scheme design, greatly shortens the combustion system design.

Claims

1. A combustion system design method based on a target heat release rate, comprising following steps: step S1: obtaining an optimal and ideal heat release rate based on a Miller and Sabathe cycle coupled model; step S2: simulating the optimal and ideal heat release rate based on a double-Wiebe function and obtaining the target heat release rate; step S3: constructing a mapping relation among the target heat release rate, piston geometric parameters and oil injection parameters, wherein, the mapping relation comprises a target combustion starting point being an oil injection timing and ignition delay period function, a premixed combustion duration being a combustion chamber throat radius, oil injection pressure and injection orifice diameter function, a premixed combustion mass being an oil injection pressure and injection orifice diameter function, and diffusion combustion being a piston pit depth function; step S4: solving target piston geometric parameters and target oil injection parameters based on the mapping relation and the target heat release rate; and step S5: designing a combustion system according to the target piston geometric parameters and the target oil injection parameters.

2. The combustion system design method based on the target heat release rate according to claim 1, wherein the step S1 specifically comprises: step S101: constructing an effective thermal efficiency calculating zero-dimensional model based on a Miller and Sabathe coupled cycle; and step S102: analyzing an influence rule of an isochoric degree and a Miller degree on effective thermal efficiency based on the effective thermal efficiency calculating zero-dimensional model, and obtaining the optimal and ideal heat release rate.

3. The combustion system design method based on the target heat release rate according to claim 2, wherein the isochoric degree is a ratio of an isochoric combustion fuel mass to a total fuel mass, and the Miller degree is an IVC time volume ratio of a conventional cycle to a Miller cycle.

4. The combustion system design method based on the target heat release rate according to claim 1, wherein the step S2 specifically comprises: step S201: simulating the ideal heat release rate by the double-Wiebe function based on a one-dimensional GT simulation model and obtaining Wiebe function characteristic parameters; step S202: optimizing the Wiebe function characteristic parameters in combination with preset diesel engine limiting conditions; and step S203: obtaining the target heat release rate according to the optimized Wiebe function characteristic parameters.

5. The combustion system design method based on the target heat release rate according to claim 4, wherein the Wiebe function is: $\frac{{\mathrm{dx}}_1}{d\phi }=\left[\left(m_p+1\right)*6.908*{\left(\frac{1}{2\tau }\right)}^{m_p+1}{\left(\phi -{\phi }_B\right)}^{m_p}{e}^{-6.908{\left(\frac{1}{2\tau }\right)}^{m_p+1}{\left(\phi -{\phi }_B\right)}^{m_p+1}}\right]\left(1-Q_d\right),$ x.sub.1 is a heat release volume, m.sub.p is a constant, τ is ignition delay, φ is a crank angle, φ.sub.B is a crank angle corresponding to a combustion starting point, and Q.sub.d is a combusted fuel mass fraction.

6. The combustion system design method based on the target heat release rate according to claim 1, wherein the step S3 specifically comprises: step S301: dividing a heat release rate into a premixed stage and a diffuse stage with a spray wall-impact time as a dividing point, and respectively studying relations among characteristic parameters in the premixed stage, oil injection parameters and combustion chamber shape parameters and among characteristic parameters in the diffuse stage, the oil injection parameters and the combustion chamber shape parameters; and step S302: respectively determining a target combustion starting point being an oil injection timing and ignition delay period function, a premixed combustion duration being a combustion chamber throat radius, oil injection pressure and injection orifice diameter function, and a premixed combustion mass being an oil injection pressure and injection orifice diameter function in the premixed stage, and determining diffusion combustion being a piston pit depth function in the diffuse stage.

7. The combustion system design method based on the target heat release rate according to claim 6, wherein the target combustion starting point is the oil injection timing and ignition delay period function:
θ.sub.SOC=θ.sub.SOI+θ.sub.id
θ.sub.id=C.sub.1 exp(C.sub.2/T.sub.a)d.sub.0.sup.C.sup.3ρ.sub.A.sup.C.sup.4P.sub.inj.sup.C.sup.5O.sub.2.sup.C.sup.6, θ.sub.SOC is a combustion starting point, θ.sub.SOI is an oil injection starting point, θ.sub.id is an ignition delay period, C.sub.1-C.sub.6 are constants, T.sub.a is an environmental temperature, ρ.sub.A is an environmental density, d.sub.0 is an injection orifice diameter, P.sub.inj is an injection pressure, and O.sub.2 is an oxygen concentration.

8. The combustion system design method based on the target heat release rate according to claim 6, wherein when a combustion chamber throat is a wall-impact position, the premixed combustion duration is the combustion chamber throat radius, oil injection pressure and injection orifice diameter function:
√{square root over ((r.sub.th).sup.2+(x(φ)+I.sub.cl).sup.2)}=K.sub.p((P.sub.i−P.sub.a)/ρ.sub.A).sup.0.25d.sub.0.sup.0.5t.sup.0.5, r.sub.th is a combustion chamber throat radius, x(φ) is a piston displacement at a crank corner, I.sub.cl is a clearance height, K.sub.p is a constant, P.sub.i is an injection pressure, P.sub.a is a backpressure, ρ.sub.A is an environmental density, d.sub.0 is an injection orifice diameter, and t is a duration from an oil injection starting point to a wall-impact time, namely the sum of the premixed combustion duration and an ignition delay period.

9. The combustion system design method based on the target heat release rate according to claim 6, wherein the premixed combustion mass is the oil injection pressure and injection orifice diameter function: $m_{\mathrm{net}}=k_{\rho }{d_0^2\left(P_i-P_a\right)}^{0.5}.Math.\left(0.5a.Math.\ln \left(\frac{P_i-P_a}{K_p{\rho }_A{d}_0}\right)-b\right),$ m.sub.net is a premixed combustion mass, k.sub.ρ is a model parameter, d.sub.0 is an injection orifice diameter, P.sub.i is an injection pressure, P.sub.a is a backpressure, a, b and K.sub.p are constants, and ρ.sub.A is an environmental density.

10. The combustion system design method based on the target heat release rate according to claim 6, wherein the diffusion combustion is the piston pit depth function: $\mathrm{HRR}\propto m\left(\frac{\sqrt{0.5k_{f_{w2}}}}{\sqrt[3]{V_2}}+\frac{\sqrt{0.5k_{f_{w1}}}}{\sqrt[3]{V_1}}\right),$ HRR is a diffusion combustion heat release rate, m is a diffusion combustion fuel mass, k.sub.f.sub.w1 and k.sub.f.sub.w2 are turbulence intensities corresponding to discrete volumes respectively, and V.sub.1 and V.sub.2 are discrete volumes respectively.

Description

BRIEF DESCRIPTION OF FIGURES

[0023] The above contents and/or additional aspects and advantages of the disclosure become obvious and easy to understand in descriptions of embodiments by combining drawings:

[0024] FIG. 1 is a flowchart of a combustion system design method based on a target heat release rate according to one embodiment of the disclosure;

[0025] FIG. 2 is a schematic diagram based on a Miller and Sabathe coupled cycle according to one embodiment of the disclosure;

[0026] FIG. 3 is a schematic diagram of a heat release rate according to one embodiment of the disclosure;

[0027] FIG. 4 is a schematic diagram of jet flow and a combustion chamber according to one embodiment of the disclosure; and

[0028] FIG. 5 is a schematic diagram of combustion chamber discretization according to one embodiment of the disclosure.

DETAILED DESCRIPTION

[0029] Embodiments of the disclosure are described in detail below, examples of the embodiments are shown in drawings, and mark numbers being the same or similar from beginning to end show same or similar members or members with same or similar functions. The following embodiments described in reference to the drawings are exemplary, and aim to explain the disclosure instead of limiting the disclosure.

[0030] A combustion system design method based on a target heat release rate according to one embodiment of the disclosure is described by referring drawings.

[0031] FIG. 1 is a flowchart of a combustion system design method based on a target heat release rate according to one embodiment of the disclosure.

[0032] As shown in FIG. 1, the combustion system design method based on the target heat release rate includes following steps:

[0033] Step S1: An optimal and ideal heat release rate is obtained based on a Miller and Sabathe cycle coupled model.

[0034] Furthermore, in one embodiment of the disclosure, the step S1 specifically includes:

[0035] Step S101: An effective thermal efficiency calculating zero-dimensional model is constructed based on a Miller and Sabathe coupled cycle.

[0036] Step S102: An influence rule of an isochoric degree and a Miller degree on effective thermal efficiency is analyzed based on the effective thermal efficiency calculating zero-dimensional model, and the optimal and ideal heat release rate is obtained.

[0037] Specifically, as shown in FIG. 2, optimized independent variables of the effective thermal efficiency calculating zero-dimensional model constructed on the basis of the Miller and Sabathe coupled cycle according to the embodiment of the disclosure include an isochoric degree x.sub.iv being a ratio of an isochoric combustion fuel mass to a total fuel mass and a Miller degree r.sub.m being an IVC time volume ratio of a conventional cycle to a Miller cycle. Matching between the isochoric degree and the Miller degree corresponding to the optimal thermal efficiency can be obtained by analyzing the influence rule of the isochoric degree and the Miller degree on effective thermal efficiency, thereby obtaining an optimal heat release rate for ideal cycle.

[00004]$x_{\mathrm{iv}}=\frac{m_{34}}{m_{\mathrm{tot}}}$$r_m=\frac{v_o}{v_m}$

[0038] In the step S2, the optimal and ideal heat release rate based on a double-Wiebe function is simulated to obtain a target heat release rate.

[0039] Furthermore, in one embodiment of the disclosure, the step S2 specifically includes:

[0040] Step S201: An ideal heat release rate by a double-Wiebe function based on a one-dimensional GT simulation model is simulated to obtain Wiebe function characteristic parameters.

[0041] Step S202: The Wiebe function characteristic parameters are optimized in combination with preset diesel engine limiting conditions.

[0042] Step S203: A target heat release rate is obtained according to the optimized Wiebe function characteristic parameters.

[0043] In other words, the final target heat release rate, namely the target heat release rate is obtained by simulating the ideal heat release rate by the double-Wiebe function based on the one-dimensional GT simulation model and then optimizing the Wiebe function characteristic parameters in combination with the limiting conditions. The Wiebe function is:

[00005]$\frac{{\mathrm{dx}}_1}{d\phi }=\left[\left(m_p+1\right)*6.908*{\left(\frac{1}{2\tau }\right)}^{m_p+1}{\left(\phi -{\phi }_B\right)}^{m_p}{e}^{-6.908{\left(\frac{1}{2\tau }\right)}^{m_p+1}{\left(\phi -{\phi }_B\right)}^{m_p+1}}\right]\left(1-Q_d\right)$

[0044] x.sub.l is a heat release volume, m.sub.p is a constant, τ is ignition delay, φ is a crank angle, φ.sub.B is a crank angle corresponding to a combustion starting point, and Q.sub.d is a combusted fuel mass fraction.

[0045] In a step S3, a mapping relation among the heat release rate, piston geometric parameters and oil injection parameters is constructed, where, the mapping relation includes a target combustion starting point being an oil injection timing and ignition delay period function, a premixed combustion duration being a combustion chamber throat radius, oil injection pressure and injection orifice diameter function, a premixed combustion mass being an oil injection pressure and injection orifice diameter function, and diffusion combustion being a piston pit depth function.

[0046] Furthermore, in one embodiment of the disclosure, the step S3 specifically includes:

[0047] Step S301: A heat release rate is divided into a premixed stage and a diffuse stage with a spray wall-impact time as a dividing point, and relations among characteristic parameters in the premixed stage, oil injection parameters and combustion chamber shape parameters and among characteristic parameters in the diffuse stage, the oil injection parameters and the combustion chamber shape parameters are respectively studied.

[0048] Step S302: A target combustion starting point being an oil injection timing and ignition delay period function, a premixed combustion duration being a combustion chamber throat radius, oil injection pressure and injection orifice diameter function, and a premixed combustion mass being an oil injection pressure and injection orifice diameter function in the premixed stage is determined, and diffusion combustion being a piston pit depth function in the diffuse stage is determined.

[0049] In other words, the heat release rate is divided into the premixed stage and the diffuse stage with the spray wall-impact time as the dividing point, and the relations among the characteristic parameters in the premixed stage, the oil injection parameters and the combustion chamber shape parameters and among the characteristic parameters in the diffuse stage, the oil injection parameters and the combustion chamber shape parameters are respectively studied, where, the target combustion starting point is the oil injection timing and ignition delay period function, the premixed combustion duration and the combustion chamber throat radius are in the oil injection pressure and injection orifice diameter function, the premixed combustion mass is the oil injection pressure and injection orifice diameter function, and the diffusion combustion is the piston pit depth function.

[0050] Furthermore, the combustion starting point, oil injection timing and ignition delay period function is:

θ.sub.SOC=θ.sub.SOI+θ.sub.id

θ.sub.id=C.sub.1 exp(C.sub.2/T.sub.a)d.sub.0.sup.C.sup.3ρ.sub.A.sup.C.sup.4P.sub.inj.sup.C.sup.5O.sub.2.sup.C.sup.6

[0051] θ.sub.SOC is a combustion starting point, θ.sub.SOI is an oil injection starting point, θ.sub.id is an ignition delay period, C.sub.1-C.sub.6 are constants, T.sub.a is an environmental temperature, ρ.sub.A is an environmental density, d.sub.0 is an injection orifice diameter, P.sub.inj is an injection pressure, and O.sub.2 is an oxygen concentration.

[0052] As shown in FIG. 3, it is theoretically analyzed that a premixed combustion end point is a spray wall-impact point, and when a wall-impact position is a throat position, following relations are established:

l.sub.imp=√{square root over ((r.sub.th).sup.2+(x(φ)+l.sub.cl).sup.2)}

S=K.sub.p((P.sub.i−P.sub.a)ρ.sub.A).sup.0.25d.sub.0.sup.0.5t.sup.0.5

[0053] l.sub.imp is a distance from an injection orifice to a wall-impact position, S is a spray penetration distance, r.sub.th is a combustion chamber throat radius, x(φ) is a piston displacement at a crank corner, l.sub.cl is a clearance height, K.sub.p is a constant, P.sub.i is an injection pressure, P.sub.a is a backpressure, ρ.sub.A is an environmental density, d.sub.0 is an injection orifice diameter, and t is a duration from an oil injection starting point to a spray wall-impact time.

[0054] The two equations are made to be equal, and a premixed combustion duration from the combustion starting point to the spray wall-impact time is established, and is the combustion chamber throat radius, oil injection pressure and injection orifice diameter function. Namely, the spray penetration distance is equal to the distance from the injection orifice to the wall-impact position at the premixed combustion ending time. When a combustion chamber throat is a spray wall-impact position, the spray wall-impact position and the combustion chamber throat radius are in a function corresponding to the oil injection pressure and the injection orifice diameter, which is specifically shown as below:

√{square root over ((r.sub.th).sup.2+(x(φ)+l.sub.cl).sup.2)}=K.sub.p((P.sub.i−P.sub.a)/ρ.sub.A).sup.0.25d.sub.0.sup.0.5t.sup.0.5

[0055] r.sub.th is a combustion chamber throat radius, x(φ) is a piston displacement at a crank corner, l.sub.cl is a clearance height, K.sub.p is a constant, P.sub.i is an injection pressure, P.sub.a is a backpressure, ρ.sub.A is an environmental density, d.sub.0 is an injection orifice diameter, and t is a duration from an oil injection starting point to a wall-impact time, namely the sum of the premixed combustion duration and an ignition delay period.

[0056] As shown in FIG. 3, the premixed combustion mass, oil injection pressure and injection orifice diameter function is established according to a jet flow theoretical analysis, where, the premixed combustion mass is defined as a fuel mass for premixed combustion within an equivalent ratio range on a spray periphery at the spray wall-impact time, and is solved in a manner of establishing a conical simplified one-dimensional model on the spray and performing integration on an axial direction and a radial direction, which is specifically shown as below:

[00006]$m_{\mathrm{net}}=k_{\rho }{d_0^2\left(P_i-P_a\right)}^{0.5}.Math.\left(0.5a.Math.\ln \left(\frac{P_i-P_a}{K_p{\rho }_A{d}_0}\right)-b\right)$

[0057] m.sub.net is a premixed combustion mass, k.sub.ρ is a model parameter, d.sub.0 is an injection orifice diameter, P.sub.i is an injection pressure, P.sub.a is a backpressure, a, b and K.sub.p are constants, and ρ.sub.A is an environmental density.

[0058] Furthermore, as shown in FIG. 4, it can be obtained, according to a jet flow theory, that a characteristic hybrid parameter representing a diffusion combustion mixing process is a ratio

[00007]$\frac{\sqrt{k_f}}{\sqrt[3]{V_t}}$

of a turbulence intensity to a turbulence length. As shown in FIG. 5, a piston combustion chamber is discretized, then, a characteristic hybrid parameter representing a diffusion combustion mixing process is a weighted sum of ratios of discrete turbulence intensities to a turbulence length, and thus, diffusion combustion is a piston pit depth function:

[00008]$\mathrm{HRR}\propto m\left(\frac{\sqrt{0.5k_{f_{w2}}}}{\sqrt[3]{V_2}}+\frac{\sqrt{0.5k_{f_{w1}}}}{\sqrt[3]{V_1}}\right)$

[0059] HRR is a diffusion combustion heat release rate, m is a diffusion combustion fuel mass, k.sub.f.sub.w1 and k.sub.f.sub.w2 are turbulence intensities corresponding to discrete volumes respectively, V.sub.1 and V.sub.2 are discrete volumes respectively, and the function is about a piston pit depth.

[0060] Computational formulas corresponding to V.sub.1 and V.sub.2 are as below:

[00009]$V_1={0.135}^2\pi .Math.\left(l_{\mathrm{cl}}+x\left(\phi \right)\right)-\frac{1}{3}\pi l_2^2\left(l_{\mathrm{cl}}+x\left(\phi \right)\right)$$V_2=\pi l_2^2{h}_2-\frac{1}{2}\left(\pi l_2^2{h}_2-\pi l_3^2{h}_3\right)+\frac{1}{3}\pi l_2^2\left(l_{\mathrm{cl}}+x\left(\phi \right)\right)-\frac{1}{3}\pi l_3^2{h}_3$$h_1=\cos \alpha .Math.l_{\mathrm{imp}}$$l_2=\sin \alpha .Math.l_{\mathrm{imp}}$

[0061] In a step S4, target piston geometric parameters and target oil injection parameters are solved based on the mapping relation and the target heat release rate.

[0062] In a step S5, a combustion system is designed according to the target piston geometric parameters and the target oil injection parameters.

[0063] In other words, parameter values are reversely deduced according to the mapping relation among the heat release rate, the oil injection characteristic parameters and the piston geometric characteristic parameters on the basis of the target heat release rate, and a combustion system design scheme capable of realizing the target heat release rate can be obtained according to the parameter values.

[0064] In conclusion, the combustion system design method based on the target heat release rate according to the embodiment of the disclosure is a combustion system forward design method, focuses on the mechanism relation between the combustion system characteristic parameters and combustion, does not depend on experience, does not need multi-scheme design, greatly shortens calculation time, shortens a combustion system research and development period, and overcomes defects of a conventional combustion system reverse design method, and a scientific and accurate combustion system design method is provided for combustion system designers.

[0065] In addition, terms “first” and “second” are merely used for describing the purpose but not understood as indicating or implying relative importance or implying to indicate the number of indicated technical features. Thus, features limited with “first” and “second” may explicitly or invisibly include at least one feature. In the description of the disclosure, “a plurality of” means at least two, such as two or three, except additional clear and specific limitations.

[0066] In the description of the specification, descriptions of reference terms “one embodiment”, “some embodiments”, “example”, “specific example” or “some examples” or the like imply that specific features, structures, materials or characteristics combined with the embodiment or the example are included in at least one embodiment or example of the disclosure. In the specification, schematic statements for the above terms are not necessarily specific to the same embodiment or example. Moreover, the described specific features, structures, materials or characteristics can be combined in a proper manner in any one or more embodiments or examples. In addition, those skilled in the art can combine different embodiments or examples described in the specification and features of different embodiments or examples under the situation that mutual contradiction does not exist.

[0067] Although the embodiments of the disclosure are shown and described above, it can be understood that the above embodiments are illustrative and cannot be understood as limitations on the disclosure, and those ordinarily skilled in the art can change, modify, replace and transform the above embodiments in the scope of the disclosure.