CN115898697A - Open expansion cycle engine and parameter estimation method and device thereof - Google Patents

Open expansion cycle engine and parameter estimation method and device thereof Download PDF

Info

Publication number
CN115898697A
CN115898697A CN202211621701.8A CN202211621701A CN115898697A CN 115898697 A CN115898697 A CN 115898697A CN 202211621701 A CN202211621701 A CN 202211621701A CN 115898697 A CN115898697 A CN 115898697A
Authority
CN
China
Prior art keywords
fuel
oxidant
turbine
power
pump
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211621701.8A
Other languages
Chinese (zh)
Inventor
崔朋
周城宏
刘鹰
刘阳
谭云涛
朱雄峰
雍子豪
韩秋龙
谷建光
吴非霏
王铁兵
李晨阳
王一杉
谭胜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
63921 Troops of PLA
Original Assignee
63921 Troops of PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 63921 Troops of PLA filed Critical 63921 Troops of PLA
Priority to CN202211621701.8A priority Critical patent/CN115898697A/en
Publication of CN115898697A publication Critical patent/CN115898697A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Control Of Turbines (AREA)

Abstract

The invention relates to an open expansion cycle engine and a parameter estimation method and device thereof, belonging to the technical field of liquid rocket engine overall scheme design. The engine disclosed by the invention divides fuel, wherein one path of fuel is heated by the cooling channel while cooling the combustion chamber, and then works on the fuel pump and the oxidant pump through the fuel turbine and the oxidant turbine in sequence to improve the pressure of the fuel and the oxidant, and the other path of fuel directly enters the thrust chamber to provide thrust through combustion work; the coupling between the chamber pressure and the turbine pressure ratio is greatly reduced, and the engine thrust is improved. On the basis of the engine structure, the method further takes turbine pump power matching as a condition, carries out key node parameter estimation based on given engine design parameters, facilitates visual understanding of different node state parameter distribution rules of the engine by people, and facilitates overall scheme design and parameter optimization of the engine, thereby developing structural parameter design of expansion cycle engine components (such as a turbine, a pump, a thrust chamber and the like).

Description

Open expansion cycle engine and parameter estimation method and device thereof
Technical Field
The disclosure relates to the technical field of expansion cycle engines in liquid rocket engines, in particular to an open expansion cycle engine and a parameter estimation method and device thereof.
Background
The circulation mode is one of the most important technical characteristics of the liquid rocket engine, and directly determines the performance, the application range and the corresponding technical scheme of the liquid rocket engine. The expansion cycle is one of three typical cycles of the liquid rocket engine, and has certain advantages compared with the other two cycle modes (gas generator cycle and afterburning cycle). Compared with the circulation of a fuel gas generator, the closed circulation is adopted, on one hand, complex thermodynamic components such as the fuel gas generator and the like are not needed, on the other hand, the waste of fuel gas which does not fully work is avoided, and the specific impulse performance is higher. Compared with afterburning circulation, complex thermodynamic components such as a precombustion chamber and the like are not needed, and the structure is simple.
Although the expansion cycle method has certain advantages, the defects are also obvious. The literature (research on a parallel type electric heating synergistic pressure-increasing variable thrust rocket engine scheme. Manned space, 2020, 26 (6): 702-9) indicates that in the mode, a propellant absorbs heat from a cooling channel, then drives a turbine to do work, and finally enters a combustion chamber for combustion, the propellant has limited work capacity after heat absorption compared with high-temperature and high-pressure gas generated after combustion, the turbine ratio is not too large, so that the expansion cycle engine chamber pressure is low, the thrust is low, and the parallel type electric heating synergistic pressure-increasing variable thrust rocket engine scheme is mainly applied to the upper stage and cannot be used in the core stage or the boosting of a carrier rocket. Therefore, the full open expansion cycle engine with the patent number of 202211255422.4 is provided with the technical scheme that the fuel absorbing heat through a cooling channel is divided into two paths, one path drives a turbine to do work and then is discharged into the environment, and the other path enters a thrust chamber to be combusted to generate thrust; so as to improve the pressure of the combustion chamber, increase the specific impulse and improve the thrust. However, the engine still has the problem of limited work capacity of the propellant after heat absorption.
Disclosure of Invention
The purpose of this disclosure is to provide an open expansion cycle engine for the problem that the structural chamber pressure of the existing full open expansion cycle engine is not high enough, thrust is not big enough; on the basis, in order to enable people to intuitively know the engine and know the change rule of the state parameter, a state parameter estimation method and a state parameter estimation device based on the engine are provided.
The purpose of the present disclosure is achieved by the following technical solutions.
In a first aspect, the present disclosure provides an open expansion cycle engine comprising a fuel turbine, a fuel pump, a fuel valve, a cooling channel, an injector, a combustion chamber, a nozzle, an oxidant turbine, an oxidant pump, an oxidant valve; the fuel is pressurized by a fuel pump and then is divided into two paths through a fuel valve, one path of the fuel enters a cooling channel to absorb heat from a combustion chamber and then sequentially drives a fuel turbine and an oxidant turbine to do work, and the other path of the fuel and the oxidant pressurized by an oxidant pump are injected into the combustion chamber through an injector through an oxidant of an oxidant valve and then are ejected out through a spray pipe to push a rocket to fly; the fuel turbine and the oxidant turbine respectively drive the fuel pump and the oxidant pump to work.
In a second aspect, the present disclosure provides a parameter estimation method, based on the engine provided in the first aspect, including:
thrust chamber pressure based on enginep c Thrust chamber mixing ratioMRAnd nozzle outlet pressurep e Calculating thrust chamber specific impulseI spc
Based on engine thrustFMRFlow ratio of thrust chamber of engine
Figure 100002_DEST_PATH_IMAGE002
AndI spc the flow of the oxidant in the thrust chamber is calculated by>
Figure 100002_DEST_PATH_IMAGE004
And the fuel flow is greater or less>
Figure 100002_DEST_PATH_IMAGE006
And turbo functional flow>
Figure 100002_DEST_PATH_IMAGE008
:/>
Figure 100002_DEST_PATH_IMAGE010
Initializing a fuel turbo-pressure ratio and an oxidant turbo-pressure ratio;
obtaining cooling channel exit temperatureT rc
Calculating the power of a fuel pump, calculating the power of a fuel turbine based on the outlet temperature of the cooling channel, comparing the power of the fuel pump with the power of the fuel turbine, if the power of the fuel pump is not matched with the power of the fuel turbine, adjusting the pressure ratio of the fuel turbine, and repeating the steps until the power of the fuel pump is matched with the power of the fuel turbine;
and calculating and comparing the power of the oxidant pump and the power of the oxidant turbine, and if the power of the oxidant pump is not matched with the power of the oxidant turbine, adjusting the pressure ratio of the oxidant turbine until the power of the oxidant pump is matched with the power of the oxidant turbine.
In a third aspect, the present disclosure provides a parameter estimation apparatus, based on the engine provided in the first aspect, including:
input module for obtaining engine thrustFNozzle outlet pressurep e And the working condition parameters are as follows: room pressurep c Mixing ratio ofMRFlow ratio of
Figure 100002_DEST_PATH_IMAGE002A
Cooling channel outlet temperatureT rc
Thrust chamber specific impulse calculation module for calculating the thrust based onp cMRAndp e calculating thrust chamber specific impulseI spc
A flow calculation module for calculating a flow based onFMR
Figure 100002_DEST_PATH_IMAGE002AA
AndI spc the flow of the oxidant in the thrust chamber is calculated by>
Figure 100002_DEST_PATH_IMAGE004A
And the fuel flow is greater or less>
Figure 100002_DEST_PATH_IMAGE006A
And turbo functional flow->
Figure 100002_DEST_PATH_IMAGE008A
Figure 100002_DEST_PATH_IMAGE010A
The turbine pressure ratio initialization module: for initializing a fuel turbo-pressure ratio and an oxidant turbo-pressure ratio;
the fuel turbine pump matching module is used for calculating the power of the fuel pump, giving the outlet temperature of the cooling channel, calculating the power of the fuel turbine, comparing the power of the fuel pump with the power of the fuel turbine, if the power of the fuel pump is not matched with the power of the fuel turbine, adjusting the pressure ratio of the fuel turbine, and repeating the module until the power of the fuel pump is matched with the power of the fuel turbine;
and the oxidant turbine pump matching module is used for calculating and comparing the power of the oxidant pump and the power of the oxidant turbine, and if the power of the oxidant pump is not matched with the power of the oxidant turbine, adjusting the pressure ratio of the oxidant turbine until the power of the oxidant pump is matched with the power of the oxidant turbine.
In a fourth aspect, the present disclosure provides an electronic device comprising:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of any of the embodiments of the second aspect.
In a fifth aspect, the present disclosure provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, causes the processor to perform the method according to any of the embodiments of the second aspect.
In a sixth aspect, the present disclosure provides a computer program product comprising a computer program/instructions which, when executed by a processor, cause the processor to perform the method of any of the embodiments of the second aspect.
Advantageous effects
Compared with the prior art, the fuel is divided, one path of the fuel is heated by self heat absorption while cooling the combustion chamber through the cooling channel, the fuel pump and the oxidant pump are driven to rotate by the fuel turbine and the oxidant turbine in sequence, the pressure of the fuel and the oxidant is increased, and the other path of the fuel directly enters the thrust chamber to provide thrust through combustion expansion work; when the coupling between the chamber pressure and the turbine pressure ratio is reduced, the working medium acting capacity of the turbine driving the fuel pump and the oxidant pump to rotate is improved, so that the pressure of fuel and oxidant entering the thrust chamber is improved, the chamber pressure is further improved, the specific impulse is increased, and the engine thrust is improved. On the basis of the engine, in order to facilitate visual understanding of the change rule of different node state parameters of the engine, a parameter estimation method is provided, based on given engine design and working condition parameters, node state parameter estimation is carried out by taking turbine pump power matching as a condition and taking a turbine pressure ratio as an adjustment variable (different from the determination of overall structure parameters during thrust adjustment of the existing expansion cycle engine, and the estimation of the node state parameters is carried out by taking the flow of a driving turbine as the adjustment variable, such as Trypop, liqing, cheng, et al. Liquid oxygen methane expansion cycle variable thrust engine system scheme contrast research [ J ]. University of defense science and technology, 2020, 42 (3): 106-15), overall scheme design and parameter optimization are conveniently carried out on the engine, so that the design of the structure parameters of subsequent expansion cycle engine components (such as a turbine, a pump, a thrust chamber and the like) is facilitated.
Drawings
FIG. 1 is a schematic illustration of an open expansion cycle engine according to an embodiment of the present disclosure;
FIG. 2 is a schematic flow chart of a parameter estimation method based on the engine shown in FIG. 1 according to an embodiment of the present disclosure;
FIG. 3 is a schematic illustration of a distribution of state parameters for an open expansion cycle engine provided in accordance with an embodiment of the present disclosure;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present disclosure.
Detailed Description
The present disclosure will be described in detail below with reference to specific embodiments shown in the drawings. These embodiments are not limited to the disclosure, and structural, methodological, or functional changes made by those skilled in the art according to these embodiments are all included in the scope of the disclosure.
In the description of the present disclosure, it is to be understood that the terms "central," "longitudinal," "lateral," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," and the like are used in the orientations and positional relationships indicated in the drawings, which are based on the orientations and positional relationships indicated in the drawings, and are used merely to facilitate the description of the disclosure and to simplify the description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed in a particular orientation, and be operated, and therefore should not be construed as limiting the disclosure. Furthermore, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or to implicitly indicate a number of the indicated technical features. Thus, a feature defined as "first," "second," etc. may explicitly or implicitly include one or more of that feature. In the description of the present disclosure, the meaning of "a plurality" is two or more unless otherwise specified.
For the purpose of illustrating the objects, technical solutions and advantages of the embodiments of the present disclosure, the technical solutions in the embodiments of the present disclosure will be clearly and completely described below with reference to the drawings in the embodiments of the present disclosure.
Fig. 1 is an open expansion cycle engine provided by the present disclosure, which as shown comprises a fuel turbine 2, a fuel pump 1, a fuel valve 3, a cooling channel 5, an injector 7, a combustion chamber 6, a nozzle 4, an oxidant turbine 10, an oxidant pump 9, an oxidant valve 8. Wherein the cooling channel 5, the injector 7, the combustion chamber 6 and the nozzle 4 form a thrust chamber. The fuel pump 1 is respectively connected with the fuel turbine 2 and the fuel valve 3, the oxidant pump 9 is respectively connected with the oxidant turbine 10 and the oxidant valve 8, the cooling channel 5 is arranged at the periphery of the thrust chamber to cool the thrust chamber, the outlet of the cooling channel 5 is connected with the fuel turbine 2, the fuel turbine 2 is connected with the oxidant turbine 10, the fuel is divided into two paths at the outlet of the fuel valve 3, one path enters the cooling channel 5, the other path and the oxidant at the outlet of the oxidant valve 8 enter the injector 7 together, the injector 7 injects the fuel and the oxidant into the combustion chamber 6, and high-temperature and high-pressure gas generated by combustion is discharged from the nozzle 4.
As can be seen from the above, the fuel is pressurized by the fuel pump 1 and divided into two paths by the fuel valve 3, one path enters the cooling channel 5 to absorb heat from the combustion chamber 6 and then sequentially drives the fuel turbine 2 and the oxidant turbine 10 to do work, and the other path and the oxidant pressurized by the oxidant pump 9 and passing through the oxidant valve 8 are injected into the combustion chamber 6 by the injector 7 and then are injected by the nozzle 4 to push the rocket to fly; the fuel turbine 2 and the oxidizer turbine 10 drive the fuel pump 1 and the oxidizer pump 9 to operate, respectively.
The engine chamber pressure is not high enough to current full open type expansion cycle engine chamber, thrust big problem inadequately, this disclosure proposes the improvement scheme to current full open type expansion cycle, through the promotion acting ability that can be better after the coolant temperature of analysis drive turbine acting promotes, only introduce the fuel of driving turbine acting one way into cooling channel and absorb heat, another way directly gets into the combustion chamber burning, thereby when keeping reducing the coupling between room pressure and the turbine pressure ratio, the working medium acting ability of the turbine of drive fuel pump and oxidant pump promotes, and then promote the pressure of fuel and oxidant, reach the effect that promotes the combustion chamber pressure, increase specific impulse, promote engine thrust. In the circulation mode, most of propellant (such as liquid hydrogen) directly enters the thrust chamber after passing through the main valve, and a small part of propellant absorbs heat to become high-temperature and high-pressure gas after passing through the cooling channel, so that the working capacity is enhanced, and the high-temperature and high-pressure gas is directly discharged into the atmosphere after driving the turbine.
In order to have an intuitive overall understanding of the open expansion cycle engine, state parameter calculations are first performed by power balancing to obtain the state parameters (flow, pressure, temperature, power) of the key nodes of the engine. The scheme contrast research [ J ] of a liquid oxygen methane expansion cycle variable thrust engine system, the university of national defense science and technology, 2020, 42 (3): 106-15.) shows a state parameter calculation method of an expansion cycle engine, because the flow rate of a coolant in an expansion cycle is the flow rate of fuel in a thrust chamber, and therefore working condition parameters in the method only need to consider the chamber pressure and the mixing ratio and do not need to consider the flow rate ratio, and the power balance is realized by adjusting the turbine flow rate ratio (the ratio of the flow rate of methane driving a turbine to the total flow rate of methane). Because of their different structures, the above state parameter calculation methods cannot be applied to the expansion cycle engine proposed by the present disclosure, and a parameter estimation method suitable for an open expansion cycle engine is needed to develop a general scheme design and parameter optimization based on room pressure, mixture ratio, flow ratio (coolant flow rate/thrust chamber flow rate) for the open expansion engine, so as to further develop a structural parameter design of expansion cycle engine components (including turbine, pump, thrust chamber, etc.).
FIG. 2 illustrates a parameter estimation method based on the open expansion cycle engine of FIG. 1, including the following:
1. thrust chamber pressure based on enginep c Thrust chamber mixing ratioMRAnd nozzle outlet pressurep e Calculating thrust chamber specific impulseI spc
I spc Can be obtained by a thermodynamic calculation method and can be referred to as a reference (liquid rocket engine design published in Beijing aerospace university Press in 2011 of Caizian national Biao, li Jia Wen, tian Aimei, zhang lihui and the like); software may also be available, such as rocket engine performance analysis software. This example was obtained using Rocket engine performance Analysis software, socket performance Analysis (RPA, www.
2. Based on engine thrustFMRFlow ratio of thrust chamber of engine
Figure DEST_PATH_IMAGE002AAA
AndI spc the flow of the oxidant in the thrust chamber is calculated by>
Figure DEST_PATH_IMAGE004AA
And the fuel flow is greater or less>
Figure DEST_PATH_IMAGE006AA
And turbo functional flow->
Figure DEST_PATH_IMAGE008AA
:/>
Figure DEST_PATH_IMAGE010AA
3. A fuel turbo-pressure ratio and an oxidant turbo-pressure ratio are initialized.
It can be initialized according to a preset value, or can be an externally given value, i.e. an input value. This value is continually adjusted in subsequent steps to accommodate the relevant design parameters given in 1.
4. Calculating fuel pump power based on cooling gallery outlet temperatureT rc And calculating the power of the fuel turbine, comparing the power of the fuel pump with the power of the fuel turbine, if the power of the fuel pump is not matched with the power of the fuel turbine, adjusting the pressure ratio of the fuel turbine, and repeating the steps until the power of the fuel pump is matched with the power of the fuel turbine.
In particular, fuel pump powerP pf Can be calculated by the following formula:
Figure 100002_DEST_PATH_IMAGE012
wherein, the first and the second end of the pipe are connected with each other,
Figure 100002_DEST_PATH_IMAGE014
for fuel injection pressure drop, by formula>
Figure 100002_DEST_PATH_IMAGE016
Obtaining; />
Figure 100002_DEST_PATH_IMAGE018
For the fuel valve pressure drop, by means of the formula>
Figure 100002_DEST_PATH_IMAGE020
Obtaining; />
Figure 100002_DEST_PATH_IMAGE022
And &>
Figure 100002_DEST_PATH_IMAGE024
Is a constant coefficient; />
Figure 100002_DEST_PATH_IMAGE026
Is the fuel pump inlet pressure, at a preset constant; />
Figure 100002_DEST_PATH_IMAGE028
Is the fuel density; />
Figure 100002_DEST_PATH_IMAGE030
Is the fuel pump efficiency, is a predetermined constant.
In particular, fuel turbine powerP tf Can be calculated by the following formula:
Figure 100002_DEST_PATH_IMAGE032
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE034
based on pressure, for a fuel turbine gas specific heat ratio>
Figure 100002_DEST_PATH_IMAGE036
And cooling channel exit temperatureT rc Get->
Figure 100002_DEST_PATH_IMAGE038
For cooling the channel pressure drop, by means of the formula>
Figure 100002_DEST_PATH_IMAGE040
Get and->
Figure 100002_DEST_PATH_IMAGE042
Is a constant coefficient; />
Figure 100002_DEST_PATH_IMAGE044
Based on pressure->
Figure 100002_DEST_PATH_IMAGE036A
AndT rc obtaining; />
Figure 100002_DEST_PATH_IMAGE046
Is a predetermined constant for fuel turbine efficiency.
When in useP pf And withP tf When not equal, i.e. not matching, the fuel turbo pressure ratio needs to be adjustedπ f Post-recalculationP tf Until the two match. The adjustment mode can be external resetting or adjustment according to a preset step length. Furthermore, in order to improve the matching efficiency, the step length is set to be a variable step length. Step size according toP pf AndP tf the difference value of (a) is calculated according to a preset proportion.
The gas specific heat ratio and the gas constant can be obtained by using query software, in this example
Figure 100002_DEST_PATH_IMAGE048
And &>
Figure DEST_PATH_IMAGE044A
Utilizing thermophysical property query software REFPROP (w. Boulder. Last. Gov/div838/the same/refop/frequencyly) according to the input
Figure DEST_PATH_IMAGE036AA
AndT rc and (4) obtaining.
5. And calculating and comparing the power of the oxidant pump and the power of the oxidant turbine, if the power of the oxidant pump is not matched with the power of the oxidant turbine, adjusting the pressure ratio of the oxidant turbine, and repeating the steps until the power of the oxidant pump is matched with the power of the oxidant turbine.
In particular, the oxidant pump powerP pox Can be calculated by the following formula:
Figure 100002_DEST_PATH_IMAGE050
wherein, the first and the second end of the pipe are connected with each other,
Figure 100002_DEST_PATH_IMAGE052
for oxidant injection pressure drop, by formula>
Figure 100002_DEST_PATH_IMAGE054
Obtaining; />
Figure 100002_DEST_PATH_IMAGE056
For oxidant valve pressure drop, by formula>
Figure 100002_DEST_PATH_IMAGE058
Obtaining; />
Figure 100002_DEST_PATH_IMAGE060
And &>
Figure 100002_DEST_PATH_IMAGE062
Is a constant coefficient; />
Figure 100002_DEST_PATH_IMAGE064
Is the oxidizer pump inlet pressure, at a predetermined constant; />
Figure 100002_DEST_PATH_IMAGE066
Is the oxidant density; />
Figure 100002_DEST_PATH_IMAGE068
Is the oxidizer pump efficiency, is a predetermined constant.
In particular, oxidant turbine powerP tox Can be calculated by the following formula:
Figure 100002_DEST_PATH_IMAGE070
wherein the content of the first and second substances,
Figure 100002_DEST_PATH_IMAGE072
based on pressure, is the oxidant turbo gas specific heat ratio>
Figure 100002_DEST_PATH_IMAGE074
And oxidant turbine inlet temperature->
Figure 100002_DEST_PATH_IMAGE076
Get and->
Figure DEST_PATH_IMAGE076A
Pass and/or>
Figure 100002_DEST_PATH_IMAGE078
Calculating a formula; />
Figure 100002_DEST_PATH_IMAGE080
Is an oxidant turbine gas constant, based on pressure
Figure 100002_DEST_PATH_IMAGE074A
And &>
Figure DEST_PATH_IMAGE076AA
Obtaining; />
Figure DEST_PATH_IMAGE082
For oxidant turbine efficiency, a predetermined constant is used.
When in useP pox AndP tox when not equal, i.e. not matching, the oxidant turbo pressure ratio needs to be adjustedπ ox Post-recalculationP tox Until the two match. The adjustment mode can be external resetting or adjustment according to a preset step length. In the same way as above, to improve the matching efficiency, the step length is set to be a variable step length. Step size is according toP pox AndP tox the difference value of (a) is calculated according to a preset proportion.
The same as the above, this example
Figure 100002_DEST_PATH_IMAGE072A
And &>
Figure DEST_PATH_IMAGE080A
Using thermophysical property query software REFPROP to input according to
Figure 100002_DEST_PATH_IMAGE074AA
And &>
Figure DEST_PATH_IMAGE076AAA
And (4) obtaining.
Through the process, after the respective turbopump powers of the fuel and the oxidant are equal, the engine reaches the optimal state under the given parameter combination, and the state parameters of each key node can be obtained according to the optimal state parameters.
In a specific embodiment, the method further includes a step of outputting the key node state parameter. Specifically, the key node state parameters include flow, pressure, temperature and turbopump power at the following nodes: fuel pump front, fuel pump back, fuel valve back (injector front), cooling channel inlet, cooling channel outlet, fuel turbine back, oxidant valve back, oxidant pump front, combustion chamber.
The state parameters of each key node can be obtained by a person skilled in the art according to the relevant knowledge of the engine design. It can also be obtained from the following provided in this example:
for the flow rate:
the flow rates of the oxidant behind the oxidant valve, the oxidant pump and the oxidant pump are equal to the flow rate of the oxidant in the thrust chamber
Figure DEST_PATH_IMAGE004AAA
Equal;
flow rate of combustion chamber
Figure DEST_PATH_IMAGE084
The flow rates before and after the fuel pump are the same, is
Figure DEST_PATH_IMAGE086
The inlet of the cooling channel, the outlet of the cooling channel, the flow behind the fuel turbine and the flow behind the oxidant turbine are all equal to the flow of the turbine work mass
Figure DEST_PATH_IMAGE008AAA
Equal;
fuel valve back (injector front) flow and fuel flow
Figure DEST_PATH_IMAGE006AAA
Are equal.
For the pressure:
the pressure before the fuel pump is the above
Figure DEST_PATH_IMAGE026A
The pressure before the oxidant pump is equal to>
Figure DEST_PATH_IMAGE064A
The pressure of the combustion chamber is as abovep c
Post fuel pump pressure of
Figure DEST_PATH_IMAGE088
Fuel valveBoth door back (injector front) and cooling channel entrance pressures
Figure DEST_PATH_IMAGE090
Cooling channel outlet pressure of
Figure DEST_PATH_IMAGE036AAA
Oxidant post-pump pressure of
Figure DEST_PATH_IMAGE092
Oxidant valve back pressure of
Figure DEST_PATH_IMAGE094
;/>
Fuel turbine back pressure of
Figure DEST_PATH_IMAGE074AAA
Oxidant turbine back pressure of
Figure DEST_PATH_IMAGE096
For the temperature:
combustion chamber temperature based onp cMRAndp e calculated by a thermodynamic calculation method or software (such as RPA);
oxidant after-turbine temperature of
Figure DEST_PATH_IMAGE098
The outlet temperature of the cooling channel is as described aboveT rc
The inlet temperature of the cooling channel, the temperature behind the fuel valve (in front of the injector) and the temperature behind the fuel pump are the same as the temperature in front of the fuel pump and are set values;
the temperature behind the oxidant valve and the temperature behind the oxidant pump are the same as the temperature before the oxidant pump and are set values;
fuel post-turbine temperature as above
Figure DEST_PATH_IMAGE076AAAA
Specific calculation examples are given below according to the above-described method.
The input parameters are shown in the following table:
TABLE 1 input parameters
Serial number Parameter (Unit) Numerical value Serial number Parameter(s) Numerical value
1 F(kN) 1000 11
Figure DEST_PATH_IMAGE062A
0.1
2 p e (MPa) 0.06 12 η tf 0.6
3 p c (MPa) 10 13 η tox 0.6
4 MR 6 14 η pf 0.7
5 α 0.17 15 η pox 0.7
6
Figure DEST_PATH_IMAGE022A
0.1 16
Figure DEST_PATH_IMAGE028A
(kg·m -3
70.37
7
Figure DEST_PATH_IMAGE060A
0.2 17
Figure DEST_PATH_IMAGE066A
(kg·m -3
1142.6
8
Figure DEST_PATH_IMAGE024A
0.1 18
Figure DEST_PATH_IMAGE026AA
(MPa)
0.3
9
Figure DEST_PATH_IMAGE042A
0.2 19
Figure DEST_PATH_IMAGE064AA
(MPa)
0.3
10 T rc (K) 500
Using the above method, the open expansion cycle engine state parameter profiles, including flow, pressure, temperature at each node, and turbopump power, can be calculated according to the conditions in the table, as shown in FIG. 3. When the turbine and pump powers are matched separately, the oxygen turbine power reaches 3670kW and the fuel turbine power reaches 10706kW. The overall recipe design and parameter optimization can then be adjusted based on the parameter distribution.
The present disclosure also provides a parameter estimation device based on the engine shown in fig. 1, including:
input module for obtaining engine thrustFNozzle outlet pressurep e And cooling channel exit temperatureT rc And the working condition parameters are as follows: room pressurep c And mixing ratio ofMRFlow ratio of
Figure DEST_PATH_IMAGE002AAAA
A thrust chamber specific impulse calculating module for calculating a thrust based onp cMRAndp e calculating thrust chamber specific impulseI spc
A flow calculation module for calculating a flow based onFMR
Figure DEST_PATH_IMAGE002_5A
AndI spc the flow of the oxidant in the thrust chamber is calculated by>
Figure DEST_PATH_IMAGE004AAAA
And the fuel flow is greater or less>
Figure DEST_PATH_IMAGE006AAAA
And turbo functional flow->
Figure DEST_PATH_IMAGE008AAAA
:/>
Figure DEST_PATH_IMAGE010AAA
The turbine pressure ratio initialization module: for initializing a fuel turbo-pressure ratio and an oxidant turbo-pressure ratio;
the fuel turbine pump matching module is used for calculating the power of the fuel pump, calculating the power of the fuel turbine based on the outlet temperature of the cooling channel, comparing the power of the fuel pump with the power of the fuel turbine, if the power of the fuel pump is not matched with the power of the fuel turbine, adjusting the pressure ratio of the fuel turbine, and repeating the module until the power of the fuel pump is matched with the power of the fuel turbine;
and the oxidant turbine pump matching module is used for calculating and comparing the power of the oxidant pump and the power of the oxidant turbine, if the power of the oxidant pump and the power of the oxidant turbine are not matched, adjusting the pressure ratio of the oxidant turbine, and repeating the module until the power of the oxidant pump and the power of the oxidant turbine are matched.
Alternatively, the fuel and oxidant turbine pressure ratios are initialized according to preset values or inputs.
Optionally, the adjustment is according to a preset step length.
Optional, fuel pump powerP pf Calculated by the following formula:
Figure DEST_PATH_IMAGE012A
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE014A
for fuel injection pressure drop, by formula>
Figure DEST_PATH_IMAGE016A
Obtaining; />
Figure DEST_PATH_IMAGE018A
For the fuel valve pressure drop, by means of the formula>
Figure DEST_PATH_IMAGE020A
Obtaining; />
Figure DEST_PATH_IMAGE022AA
And &>
Figure DEST_PATH_IMAGE024AA
Is a constant coefficient; />
Figure DEST_PATH_IMAGE026AAA
Is the fuel pump inlet pressure, at a preset constant; />
Figure DEST_PATH_IMAGE028AA
Is the fuel density; />
Figure DEST_PATH_IMAGE030A
Is the fuel pump efficiency, is a preset constant.
Alternatively, said fuel turbine powerP tf Calculated by the following formula:
Figure DEST_PATH_IMAGE032A
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE034A
based on pressure, for a fuel turbine gas specific heat ratio>
Figure DEST_PATH_IMAGE036AAAA
And cooling channel exit temperatureT rc Get and->
Figure DEST_PATH_IMAGE038A
For cooling the channel pressure drop, by means of the formula>
Figure DEST_PATH_IMAGE040A
Get and->
Figure DEST_PATH_IMAGE042AA
Is a constant coefficient; />
Figure DEST_PATH_IMAGE044AA
Based on pressure->
Figure DEST_PATH_IMAGE036_5A
AndT rc obtaining; />
Figure DEST_PATH_IMAGE046A
Is a predetermined constant for fuel turbine efficiency.
Optionally, oxidant pump powerP pox Calculated by the following formula:
Figure DEST_PATH_IMAGE050A
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE052A
for oxidant injection pressure drop, by formula>
Figure DEST_PATH_IMAGE054A
Obtaining; />
Figure DEST_PATH_IMAGE056A
For oxidant valve pressure drop, by formula>
Figure DEST_PATH_IMAGE058A
Obtaining; />
Figure DEST_PATH_IMAGE060AA
And &>
Figure DEST_PATH_IMAGE062AA
Is a constant coefficient; />
Figure DEST_PATH_IMAGE064AAA
The oxidant pump inlet pressure is a preset constant; />
Figure DEST_PATH_IMAGE066AA
Is the oxidant density; />
Figure DEST_PATH_IMAGE068A
Is the oxidizer pump efficiency, is a predetermined constant.
Optionally, said oxidant turbine powerP tox Calculated by the following formula:
Figure DEST_PATH_IMAGE070A
/>
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE072AA
based on pressure, is the oxidant turbo gas specific heat ratio>
Figure DEST_PATH_IMAGE074AAAA
And oxidant turbine inlet temperature->
Figure DEST_PATH_IMAGE076_5A
Get and->
Figure DEST_PATH_IMAGE076_6A
Pass and/or>
Figure DEST_PATH_IMAGE078A
Calculating a formula; />
Figure DEST_PATH_IMAGE080AA
Is an oxidant turbine gas constant, based on pressure
Figure DEST_PATH_IMAGE074_5A
And &>
Figure DEST_PATH_IMAGE076_7A
Obtaining; />
Figure DEST_PATH_IMAGE082A
For oxidant turbine efficiency, a predetermined constant is set.
Optionally, the apparatus further includes a key node state parameter output module, where the key node state parameters include flow, pressure, temperature, and power of the following nodes: fuel pump front, fuel pump back, fuel valve back (injector front), cooling channel inlet, cooling channel outlet, fuel turbine back, oxidant valve back, oxidant pump front, combustion chamber.
For the apparatus embodiment, since it is substantially similar to the method embodiment, the description is relatively simple, and reference may be made to the partial description of the method embodiment for relevant points.
Fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present disclosure, where the electronic device may execute the processing flow provided by the foregoing method embodiment, and as shown in fig. 4, the electronic device 110 includes: memory 111, processor 112, computer programs, and communications interface 113; wherein the computer program is stored in the memory 111 and is configured to be executed by the processor 112 for performing the method as described above.
In addition, the embodiment of the present disclosure also provides a computer readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to implement the method described in the above embodiment.
Embodiments of the present disclosure also provide a computer program product comprising a computer program/instructions which, when executed by a processor, cause the processor to perform the method as described above.
Those of ordinary skill in the art will understand that: all or a portion of the steps of implementing the above-described method embodiments may be performed by hardware associated with program instructions. The foregoing program may be stored in a computer-readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
Specific examples are given in this specification for the purpose of illustrating the disclosure and the manner of practicing the disclosure. The details introduced in the examples are not intended to limit the scope of the claims but rather to aid in understanding the present disclosure. Those skilled in the art will understand that: although the description is given in terms of embodiments, not every embodiment includes only a single embodiment, and such description is for clarity only, and those skilled in the art will recognize that the embodiments described herein may be combined as a whole to form other embodiments as would be understood by those skilled in the art. And that various modifications, changes, or alterations to the steps of the preferred embodiments are possible without departing from the spirit and scope of this disclosure and the appended claims. Therefore, the disclosure should not be limited to the disclosure of the preferred embodiment and the drawings.

Claims (10)

1. An open expansion cycle engine characterized by: the method comprises the following steps:
the fuel injection device comprises a fuel turbine (2), a fuel pump (1), a fuel valve (3), a cooling channel (5), an injector (7), a combustion chamber (6), a spray pipe (4), an oxidant turbine (10), an oxidant pump (9) and an oxidant valve (8);
the fuel is pressurized by a fuel pump (1) and then is divided into two paths after passing through a fuel valve (3); one path enters a cooling channel (5) to absorb heat from a combustion chamber (6) and then sequentially drives a fuel turbine (2) and an oxidant turbine (10) to do work and discharge the work into the environment; the other path of the mixed gas is pressurized by an oxidant pump (9) and is injected into a combustion chamber (6) by an oxidant valve (8) through an injector (7) and then is sprayed out by a spray pipe (4); the fuel turbine (2) and the oxidant turbine (10) respectively drive the fuel pump (1) and the oxidant pump (9) to work.
2. A parameter estimation method based on the open expansion cycle engine of claim 1, characterized in that: the method comprises the following steps:
thrust chamber pressure based on enginep c Thrust chamber mixing ratioMRAnd nozzle outlet pressurep e Calculating thrust chamber specific impulseI spc
Based on engine thrustFMRFlow ratio of thrust chamber of engine
Figure DEST_PATH_IMAGE002
AndI spc calculating the oxidant flow in the thrust chamber by>
Figure DEST_PATH_IMAGE004
And the fuel flow is greater or less>
Figure DEST_PATH_IMAGE006
And turbo functional flow->
Figure DEST_PATH_IMAGE008
Figure DEST_PATH_IMAGE010
Initializing fuel turbo pressure ratioπ f And turbo pressure ratio of oxidantπ ox
Calculating fuel pump power based on cooling gallery outlet temperatureT rc Calculating the power of a fuel turbine, comparing the power of the fuel pump with the power of the fuel turbine, if the power of the fuel pump is not matched with the power of the fuel turbine, adjusting the pressure ratio of the fuel turbine, and repeating the contents of the sections until the power of the fuel pump is matched with the power of the fuel turbine;
and calculating and comparing the power of the oxidant pump and the power of the oxidant turbine, if the power of the oxidant pump and the power of the oxidant turbine are not matched, adjusting the pressure ratio of the oxidant turbine, and repeating the content of the segments until the power of the oxidant pump and the power of the oxidant turbine are matched.
3. The parameter estimation method according to claim 2, characterized in that: the initialization fuel and oxidant turbine pressure ratios are initialized according to preset values or input values.
4. The parameter estimation method according to claim 2, characterized in that: the adjustment is according to the preset step length adjustment.
5. The parameter estimation method according to claim 2, characterized in that: the fuel pump powerP pf Calculated by the following formula:
Figure DEST_PATH_IMAGE012
wherein the content of the first and second substances,
Figure DEST_PATH_IMAGE014
for fuel injection pressure drop, by formula>
Figure DEST_PATH_IMAGE016
Obtaining; />
Figure DEST_PATH_IMAGE018
For fuel valve pressure drop, by formula
Figure DEST_PATH_IMAGE020
Obtaining; />
Figure DEST_PATH_IMAGE022
And &>
Figure DEST_PATH_IMAGE024
Is a constant coefficient; />
Figure DEST_PATH_IMAGE026
Is the fuel pump inlet pressure, at a preset constant; />
Figure DEST_PATH_IMAGE028
Is the fuel density; />
Figure DEST_PATH_IMAGE030
Is the fuel pump efficiency, is a predetermined constant.
6. The parameter estimation method according to claim 5, characterized in that: the fuel turbine powerP tf Calculated by the following formula:
Figure DEST_PATH_IMAGE032
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE034
based on pressure, for a fuel turbine gas specific heat ratio>
Figure DEST_PATH_IMAGE036
AndT rc get and->
Figure DEST_PATH_IMAGE038
For cooling the channel pressure drop, by means of the formula>
Figure DEST_PATH_IMAGE040
Get and->
Figure DEST_PATH_IMAGE042
Is a constant coefficient; />
Figure DEST_PATH_IMAGE044
Is a fuel turbine gas constant based on pressure
Figure DEST_PATH_IMAGE036A
AndT rc obtaining; />
Figure DEST_PATH_IMAGE046
Is a predetermined constant for fuel turbine efficiency.
7. The parameter estimation method according to claim 6, characterized in that: the oxidant pump powerP pox Calculated by the following formula:
Figure DEST_PATH_IMAGE048
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE050
for oxidant injection pressure drop, by formula>
Figure DEST_PATH_IMAGE052
Obtaining; />
Figure DEST_PATH_IMAGE054
For oxidant valve pressure drop, by formula>
Figure DEST_PATH_IMAGE056
Obtaining; />
Figure DEST_PATH_IMAGE058
And &>
Figure DEST_PATH_IMAGE060
Is a constant coefficient; />
Figure DEST_PATH_IMAGE062
The oxidant pump inlet pressure is a preset constant; />
Figure DEST_PATH_IMAGE064
Is the oxidant density; />
Figure DEST_PATH_IMAGE066
Is the oxidizer pump efficiency, is a predetermined constant.
8. The parameter estimation method according to claim 7, characterized in that: the oxidant turbine powerP tox Calculated by the following formula:
Figure DEST_PATH_IMAGE068
wherein, the first and the second end of the pipe are connected with each other,
Figure DEST_PATH_IMAGE070
based on pressure on the specific heat ratio of the oxidant turbine gas>
Figure DEST_PATH_IMAGE072
And oxidant turbine inlet temperature>
Figure DEST_PATH_IMAGE074
Get->
Figure DEST_PATH_IMAGE074A
By>
Figure DEST_PATH_IMAGE076
Calculating a formula; />
Figure DEST_PATH_IMAGE078
Is oxidant turbine gas constant, based on pressure
Figure DEST_PATH_IMAGE072A
And &>
Figure DEST_PATH_IMAGE074AA
Obtaining; />
Figure DEST_PATH_IMAGE080
For oxidant turbine efficiency, a predetermined constant is set.
9. The parameter estimation method according to claim 8, characterized in that: the method also comprises a step of outputting the state parameters of the key nodes.
10. A parameter estimation device based on the open expansion cycle engine of claim 1, characterized in that: the method comprises the following steps:
input module for obtaining engine thrustFNozzle outlet pressurep e And cooling channel exit temperatureT rc And the working condition parameters are as follows: room pressurep c Mixing ratio ofMRFlow rate ratio
Figure DEST_PATH_IMAGE002A
Thrust chamber specific impulse calculation module for calculating the thrust based onp cMRAndp e calculating thrust chamber specific impulseI spc
A flow calculation module for calculating a flow based onFMR
Figure DEST_PATH_IMAGE002AA
AndI spc calculating the oxidant flow in the thrust chamber by>
Figure DEST_PATH_IMAGE004A
And the fuel flow is greater or less>
Figure DEST_PATH_IMAGE006A
And turbo functional flow->
Figure DEST_PATH_IMAGE008A
Figure DEST_PATH_IMAGE010A
/>
The turbine pressure ratio initialization module: for initializing a fuel turbo-pressure ratio and an oxidant turbo-pressure ratio;
the fuel turbine pump matching module is used for calculating the power of the fuel pump, calculating the power of the fuel turbine based on the outlet temperature of the cooling channel, comparing the power of the fuel pump with the power of the fuel turbine, if the power of the fuel pump is not matched with the power of the fuel turbine, adjusting the pressure ratio of the fuel turbine, and repeating the module until the power of the fuel pump is matched with the power of the fuel turbine;
and the oxidant turbine pump matching module is used for calculating and comparing the power of the oxidant pump and the power of the oxidant turbine, if the power of the oxidant pump is not matched with the power of the oxidant turbine, adjusting the pressure ratio of the oxidant turbine, and repeating the module until the power of the oxidant pump is matched with the power of the oxidant turbine.
CN202211621701.8A 2022-12-16 2022-12-16 Open expansion cycle engine and parameter estimation method and device thereof Pending CN115898697A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211621701.8A CN115898697A (en) 2022-12-16 2022-12-16 Open expansion cycle engine and parameter estimation method and device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211621701.8A CN115898697A (en) 2022-12-16 2022-12-16 Open expansion cycle engine and parameter estimation method and device thereof

Publications (1)

Publication Number Publication Date
CN115898697A true CN115898697A (en) 2023-04-04

Family

ID=86478087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211621701.8A Pending CN115898697A (en) 2022-12-16 2022-12-16 Open expansion cycle engine and parameter estimation method and device thereof

Country Status (1)

Country Link
CN (1) CN115898697A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117740391A (en) * 2024-02-21 2024-03-22 中国航发四川燃气涡轮研究院 Tail nozzle thrust test platform based on electric pump type gas generator

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100024386A1 (en) * 2007-07-20 2010-02-04 Greene William D Gas-Generator Augmented Expander Cycle Rocket Engine
US8250853B1 (en) * 2011-02-16 2012-08-28 Florida Turbine Technologies, Inc. Hybrid expander cycle rocket engine
US8407981B1 (en) * 2010-02-05 2013-04-02 Florida Turbine Technologies, Inc. Johnson Sexton cycle rocket engine
CN115324773A (en) * 2022-10-13 2022-11-11 中国人民解放军63921部队 Full open expansion cycle engine

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100024386A1 (en) * 2007-07-20 2010-02-04 Greene William D Gas-Generator Augmented Expander Cycle Rocket Engine
US8407981B1 (en) * 2010-02-05 2013-04-02 Florida Turbine Technologies, Inc. Johnson Sexton cycle rocket engine
US8250853B1 (en) * 2011-02-16 2012-08-28 Florida Turbine Technologies, Inc. Hybrid expander cycle rocket engine
CN115324773A (en) * 2022-10-13 2022-11-11 中国人民解放军63921部队 Full open expansion cycle engine

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
崔朋;李清廉;成鹏;张北辰;: "液氧甲烷膨胀循环变推力发动机系统方案对比研究", 国防科技大学学报, no. 03, 28 June 2020 (2020-06-28), pages 106 - 115 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117740391A (en) * 2024-02-21 2024-03-22 中国航发四川燃气涡轮研究院 Tail nozzle thrust test platform based on electric pump type gas generator

Similar Documents

Publication Publication Date Title
US11041439B2 (en) Hybrid expander cycle with turbo-generator and cooled power electronics
EP3623602B1 (en) Hybrid expander cycle with intercooling and turbo-generator
CA2507977C (en) Method and apparatus for operating an intercooler for a gas turbine engine
US20120039701A1 (en) Closed Cycle Brayton Cycle System and Method
CN106014637A (en) Air precooling compression aircraft engine and hypersonic velocity aircraft
CN115898697A (en) Open expansion cycle engine and parameter estimation method and device thereof
CN107035529B (en) Micro gas turbine system
Wang et al. Thermodynamic analysis of chemical precooled turbine combined engine cycle
Lv et al. Thermodynamic modeling and analysis of ammonia injection pre-compressor cooling cycle: A novel scheme for high Mach number turbine engines
RU2661427C1 (en) Bypass turbojet engine
WO2006085944A2 (en) Laser with brayton cycle outlet pump
Rähse et al. Gas dynamic simulation of shockless explosion combustion for gas turbine power cycles
CN112983681A (en) Method for rapidly calculating mass of high-thrust liquid rocket engine
CN115758627A (en) Variable working condition parameter estimation method and device for open expansion cycle engine
CN114320666B (en) Air turbine rocket engine system
CN115726903A (en) Method and device for designing overall scheme of open expansion cycle engine
CN116127815A (en) Modeling method of turbofan engine with injection nozzle
Lee et al. Modeling of the performance of a F100-PW229 equivalent engine under sea-level static conditions
Jin et al. Numerical and experimental study on shear coaxial injectors with hot hydrogen-rich gas/oxygen-rich gas and GH2/GO2
RU2197628C2 (en) Method of operation of liquid-propellant rocket engine with turbopump delivery of cryogenic propellant on basis of oxygen oxidizer and hydrocarbon fuel and liquid- propellant rocket engine implementing said method
Zhang et al. Thermodynamic analysis of a typical compressed air energy storage system coupled with a fully automatic ejector under slip pressure conditions
Delimont et al. Computational Modeling of a Direct Fired Oxy-Fuel Combustor for sCO2 Power Cycles
RU2594828C1 (en) Propulsion engine of supersonic aircraft
CN112228248A (en) Combined circulating liquid rocket engine
Lv et al. Parametric optimization and exergy analysis of a high mach number aeroengine with an ammonia mass injection pre-compressor cooling cycle

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination