A pre-cooling combined engine closed helium cycle system and digital simulation modeling method

Abstract

The application discloses a pre-cooling combined engine closed helium circulation system and a digital simulation modeling method, determines the matching working relationship of each component according to the working characteristics of the air circulation, helium circulation and hydrogen circulation of the pre-cooling combined engine, establishes an engine component-level model by using a component method, builds an engine common working equation (a balance equation) on the basis, solves the equation by using a chord division method, and obtains a total digital simulation model of the closed helium circulation system, thereby providing reliable data support for subsequent high-performance wide-envelop pre-cooling combined engine closed helium circulation system tests. The application establishes an engine component-level model by using a component method, wherein the component characteristics are ground test data, the model accuracy is ensured, seven balance equations are built on the basis, the chord division method is used for iterative solving, and the convergence performance of the iterative algorithm is improved.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of simulation modeling, and relates to a pre-cooling combined engine closed helium cycle system and a digital simulation modeling method. BACKGROUND

[0002] To meet the research and development needs of future wide flight envelope hypersonic aerospace vehicles, combined engines (turbine-ramjet-rocket) have become a research hotspot in the field of hypersonic power. The pre-cooling combined engine greatly improves the working efficiency and specific impulse of the combined engine by adding a pre-cooler to cool the high-speed incoming air, and simplifies the mode conversion problem in the combined engine. Therefore, the pre-cooling combined engine is the most promising combined engine for achieving Ma0-25 horizontal take-off and landing and air-space round-trip missions.

[0003] The current research process of the pre-cooling combined engine is stuck at the development of a closed helium cycle system in the engine for realizing heat exchange and pre-cooling. The closed helium cycle system essentially uses helium as an intermediate medium to extract the cold energy from low-temperature liquid hydrogen fuel to efficiently cool the incoming air. The helium gas, which absorbs heat and warms up, reuses the heat obtained from the high-temperature air to drive the turbine to work through expansion, Figure 1 The pre-cooling combined engine closed helium cycle system is shown in the figure. To meet the above characteristics, the closed helium cycle system couples the air, helium, and hydrogen three flow paths to work together, involving multiple components such as an air pre-cooler, a hydrogen-helium recuperator, a helium heater, an air compressor, a helium turbine, a helium compressor, a hydrogen turbine, a combustion chamber, and a nozzle, and the working process is complex.

[0004] Product testing is a necessary step for a mature engine, but the development of a combined engine closed helium cycle system test needs to be based on clear system characteristics and complete design point and off-design point matching characteristics. Therefore, a complete numerical simulation model of the combined engine closed helium cycle system needs to be established to analyze the characteristics of each component and the system matching characteristics. However, there is still no complete numerical simulation model for performance analysis of the pre-cooling combined engine closed helium cycle system. SUMMARY

[0005] The application aims to solve the problems in the prior art and provide a pre-cooling combined engine closed helium cycle system and a digital simulation modeling method.

[0006] To achieve the above purpose, the application adopts the following technical solutions:

[0007] In a first aspect, the application provides a pre-cooling combined engine closed helium cycle system and a digital simulation modeling method, including the following steps:

[0008] Step 1, according to the air compressor flow ratio function Z, the air compressor inlet total temperature T0_a4 and the air compressor inlet total pressure P0_a4, the air compressor outlet total temperature T_a5, the air compressor outlet total pressure P_a5, the air compressor power N_C and the air compressor outlet flow qm_a5 are calculated;

[0009] Step 2, according to the helium turbine inlet flow qm0_h1 and the helium turbine inlet total pressure P0_h1, the calculated value qm_h1 of the helium turbine inlet flow is calculated; the value of the helium turbine inlet total pressure P0_h1 is adjusted to balance the helium turbine inlet flow qm0_h1 and the calculated value qm_h1 of the helium turbine inlet flow; the helium turbine power N_T1 is obtained according to the helium turbine inlet total pressure P0_h1 and the calculated value qm_h1 of the helium turbine inlet flow;

[0010] Step 3, the value of the helium turbine inlet flow qm0_h1 is adjusted to balance the helium turbine power N_T1 and the air compressor power N_C;

[0011] Step 4, according to the helium compressor pressure ratio Pr_S1 and the hydrogen storage tank outlet hydrogen total pressure P_q0, the hydrogen storage tank outlet hydrogen flow qm_q1 is calculated, and according to the hydrogen storage tank outlet hydrogen flow qm_q1, the hydrogen-helium recuperator outlet helium side flow qm_h3, the hydrogen-helium recuperator outlet helium side total pressure P_h3, the hydrogen-helium recuperator outlet helium side total temperature T_h3, the hydrogen-helium recuperator outlet hydrogen side flow qm_q2, the hydrogen-helium recuperator outlet hydrogen side total pressure P_q2 and the hydrogen-helium recuperator outlet hydrogen side total temperature T_q2 are calculated; according to the helium compressor pressure ratio Pr_S1, the hydrogen-helium recuperator outlet helium side flow qm_h3, the hydrogen-helium recuperator outlet helium side total pressure P_h3 and the hydrogen-helium recuperator outlet helium side total temperature T_h3, the helium compressor outlet total temperature T_h4, the helium compressor outlet total pressure P_h4, the helium compressor outlet flow qm_h4 and the helium compressor power N_S1 are calculated;

[0012] Step 5, according to the hydrogen turbine drop pressure ratio Pr_TH1, the hydrogen-helium recuperator outlet hydrogen side flow qm_q2, the hydrogen-helium recuperator outlet hydrogen side total pressure P_q2 and the hydrogen-helium recuperator outlet hydrogen side total temperature T_q2, the hydrogen turbine outlet flow qm_q4, the hydrogen turbine outlet total temperature T_q4, the hydrogen turbine outlet total pressure P_q4 and the hydrogen turbine power N_TH1 are calculated; the value of the hydrogen turbine drop pressure ratio Pr_TH1 is adjusted to balance the hydrogen turbine power N_TH1 and the helium compressor power N_S1; the given hydrogen storage tank outlet hydrogen total pressure P_q0 is adjusted to balance the helium turbine outlet total pressure P_q4 and the air compressor outlet total pressure P_a5;

[0013] Step 6, according to the pre-cooler inlet air side flow qm_a2, the pre-cooler inlet air side total pressure P_a2, the pre-cooler inlet air side total temperature T_a2, the helium pressure compressor outlet flow qm_h4, the helium pressure compressor outlet total pressure P_h4 and the helium pressure compressor outlet total temperature T_h4, the pre-cooler outlet helium side total pressure P_h5, the pre-cooler outlet helium side total temperature T_h5, the pre-cooler outlet air side total pressure P_a3 and the pre-cooler outlet air side total temperature T_a3 are calculated; according to the pre-cooler outlet helium side total pressure P_h5, the pre-cooler outlet helium side total temperature T_h5, the pre-cooler outlet air side total pressure P_a3 and the pre-cooler outlet air side total temperature T_a3, the helium heater outlet helium total pressure P_h6 is calculated; the value of the helium pressure compressor pressure ratio Pr_S1 is adjusted to balance the helium heater outlet helium total pressure P_h6 and the helium turbine inlet total pressure P0_h1; the values of the air pressure compressor inlet total temperature T0_a4 and the air pressure compressor inlet total pressure P0_a4 are adjusted to balance the air pressure compressor inlet total temperature T0_a4 and the pre-cooler outlet air side total temperature T_a3, and balance the air pressure compressor inlet total pressure P0_a4 and the pre-cooler outlet air side total pressure P_a3;

[0014] Step 7, according to the air pressure compressor outlet flow qm_a5, the air pressure compressor outlet total pressure P_a5, the air pressure compressor outlet total temperature T_a5, the hydrogen turbine outlet flow qm_q4, the hydrogen turbine outlet total pressure P_q4 and the hydrogen turbine outlet total temperature T_q4, the main combustion chamber outlet total temperature T_b, the main combustion chamber outlet total pressure P_b and the main combustion chamber outlet flow qm_b are calculated; according to the main combustion chamber outlet total temperature T_b, the main combustion chamber outlet total pressure P_b and the main combustion chamber outlet flow qm_b, the calculated value of the tail nozzle outlet flow qm_e, the tail nozzle outlet flow velocity V_e and the tail nozzle outlet static pressure P_e are calculated; according to the flow conservation, the theoretical value of the tail nozzle outlet flow qm_t is obtained; the air pressure compressor flow ratio function Z is adjusted to balance the calculated value of the tail nozzle outlet flow qm_e and the theoretical value of the tail nozzle outlet flow qm_t;

[0015] Step 8, according to the calculated value of the tail nozzle outlet flow qm_e, the tail nozzle outlet flow velocity V_e, the tail nozzle outlet static pressure P_e, the atmospheric environment static pressure P_a0 and the tail nozzle throat area A_e, the engine performance parameters are calculated.

[0016] In a second aspect, the application provides a closed helium cycle system of a pre-cooled combined engine, comprising:

[0017] An initial value input module is configured to input an air compressor flow ratio function Z, an air compressor inlet total temperature T0_a4, an air compressor inlet total pressure P0_a4, a helium turbine inlet flow qm0_h1, a helium turbine inlet total pressure P0_h1, a helium compressor pressure ratio Pr_S1, a hydrogen storage tank outlet hydrogen total pressure P_q0, and a hydrogen turbine drop pressure ratio Pr_TH1;

[0018] A component calculation module includes an air compressor calculation module, a helium turbine calculation module, a hydrogen-helium recuperator calculation module, a helium compressor calculation module, a hydrogen turbine calculation module, a precooler calculation module, a gas generator calculation module, a helium heater calculation module, a main combustion chamber calculation module, a tail nozzle calculation module, and an engine performance calculation module;

[0019] A residual error judgment module is configured to judge whether a tail nozzle outlet flow reaches equilibrium, whether an air compressor inlet total temperature / total pressure reaches equilibrium, whether helium turbine and air compressor power reaches equilibrium, whether a helium turbine inlet flow reaches equilibrium, whether a helium heater outlet total pressure is closed loop, whether a combustion chamber inlet medium total pressure reaches equilibrium, and whether hydrogen turbine and helium compressor power reaches equilibrium.

[0020] Compared with the prior art, the present application has the following beneficial effects:

[0021] The present application determines the matching working relationship of each component according to the working characteristics of the air cycle, the helium cycle and the hydrogen cycle of the precooled combined engine, establishes an engine component-level model by using the component method, and on this basis, constructs an engine common working equation (equilibrium equation), iteratively solves by using the secant method, and obtains a closed helium cycle system overall digital simulation model, which provides reliable data support for subsequent high-performance wide envelope precooled combined engine closed helium cycle system tests. The digital simulation model of the closed helium cycle system of the precooled combined engine provided in the present application can effectively predict the performance of the engine under actual working conditions, provide reliable data support for engine ground and flight tests, provide a component performance selection idea for the construction of a test system, guide the control law of parameters in the actual test process, speed up the test planning process, and reduce the test risk. The present application establishes an engine component-level model by using the component method, wherein the component characteristics are ground test data, which ensures the accuracy of the model, and on this basis, seven equilibrium equations are established, and the secant method is used for iterative solution, which improves the convergence performance of the iterative algorithm. BRIEF DESCRIPTION OF DRAWINGS

[0022] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0023] Figure 1 This is a schematic diagram of the logic operation principle of the pre-cooled combined engine closed helium cycle system model in this application. Detailed Implementation

[0024] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. The components of the embodiments of this application described and shown in the accompanying drawings can generally be arranged and designed in various different configurations.

[0025] Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely to illustrate selected embodiments of the application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application.

[0026] It should be noted that similar labels and letters in the following figures indicate similar items. Therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures.

[0027] In the description of the embodiments of this application, it should be noted that if terms such as "upper," "lower," "horizontal," or "inner" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings, or the orientation or positional relationship commonly used when the product of the invention is in use, they are only for the convenience of describing this application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation on this application. In addition, terms such as "first" and "second" are only used to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0028] Furthermore, the use of the term "horizontal" does not imply that the component must be absolutely horizontal, but rather that it can be slightly tilted. For example, "horizontal" simply means that its direction is more horizontal than "vertical," and does not mean that the structure must be completely horizontal, but can be slightly tilted.

[0029] In the description of the embodiments of the present application, it also needs to be explained that, unless explicitly defined and limited, if the terms "set", "install", "connect", "join" appear, they should be understood in a broad sense, for example, can be fixedly connected, can be detachably connected, or integrally connected; can be mechanically connected, can be electrically connected; can be directly connected, can be indirectly connected through an intermediate medium, can be the communication inside two elements. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0030] The present application will be further described in detail below with reference to the accompanying drawings:

[0031] Referring to Figure 1 The present application discloses a pre-cooling combined engine closed helium cycle system and a digital simulation modeling method, comprising the following steps:

[0032] Step 1: Taking the air compressor inlet as the operation starting point, each component is calculated in turn.

[0033] Given the air compressor flow ratio function Z, the air compressor inlet total temperature T0_a4 and the air compressor inlet total pressure P0_a4, the air compressor outlet total temperature T_a5, the air compressor outlet total pressure P_a5, the air compressor power N_C and the air compressor outlet flow qm_a5 are calculated;

[0034] Step 2:

[0035] Given the helium turbine inlet flow qm0_h1 and the helium turbine inlet total pressure P0_h1, the calculated value qm_h1 of the helium turbine inlet flow is calculated, the value of the helium turbine inlet total pressure P0_h1 is adjusted to balance the helium turbine inlet flow qm0_h1 and the calculated value qm_h1 of the helium turbine inlet flow, and the helium turbine power N_T1 is obtained according to the helium turbine inlet total pressure P0_h1 and the calculated value qm_h1 of the helium turbine inlet flow;

[0036] Step 3:

[0037] The value of the helium turbine inlet flow qm0_h1 is adjusted to balance the helium turbine power N_T1 and the air compressor power N_C;

[0038] Step 4:

[0039] Given the helium pressure ratio Pr_S1, given the hydrogen tank outlet hydrogen total pressure P_q0, the hydrogen tank outlet hydrogen flow rate qm_q1 is calculated, and the hydrogen helium regenerator outlet helium side flow rate qm_h3, the hydrogen helium regenerator outlet helium side total pressure P_h3, the hydrogen helium regenerator outlet helium side total temperature T_h3, and the hydrogen helium regenerator outlet hydrogen side flow rate qm_q2, the hydrogen helium regenerator outlet hydrogen side total pressure P_q2, and the hydrogen helium regenerator outlet hydrogen side total temperature T_q2 are calculated according to the hydrogen helium regenerator outlet hydrogen flow rate qm_q1; according to the helium pressure ratio Pr_S1 and the hydrogen helium regenerator outlet helium side flow rate qm_h3, the hydrogen helium regenerator outlet helium side total pressure P_h3, and the hydrogen helium regenerator outlet helium side total temperature T_h3, the helium compressor outlet total temperature T_h4, the helium compressor outlet total pressure P_h4, the helium compressor outlet flow rate qm_h4, and the helium compressor power N_S1 are calculated;

[0040] Step5:

[0041] Given the hydrogen turbine drop pressure ratio Pr_TH1, and according to the hydrogen helium regenerator outlet hydrogen side flow rate qm_q2, the hydrogen helium regenerator outlet hydrogen side total pressure P_q2, and the hydrogen helium regenerator outlet hydrogen side total temperature T_q2, the hydrogen turbine outlet flow rate qm_q4, the hydrogen turbine outlet total temperature T_q4, the hydrogen turbine outlet total pressure P_q4, and the hydrogen turbine power N_TH1 are calculated; the value of the hydrogen turbine drop pressure ratio Pr_TH1 is adjusted so that the hydrogen turbine power N_TH1 and the helium compressor power N_S1 are balanced, and the given hydrogen tank outlet hydrogen total pressure P_q0 is adjusted so that the helium turbine outlet total pressure P_q4 and the air compressor outlet total pressure P_a5 are balanced;

[0042] Step6:

[0043] According to the precooler inlet air side flow rate qm_a2, the precooler inlet air side total pressure P_a2, the precooler inlet air side total temperature T_a2, the helium compressor outlet flow rate qm_h4, the helium compressor outlet total pressure P_h4, and the helium compressor outlet total temperature T_h4, the precooler outlet helium side total pressure P_h5, the precooler outlet helium side total temperature T_h5, the precooler outlet air side total pressure P_a3, and the precooler outlet air side total temperature T_a3 are calculated; according to the precooler outlet helium side total pressure P_h5, the precooler outlet helium side total temperature T_h5, the precooler outlet air side total pressure P_a3, and the precooler outlet air side total temperature T_a3, the helium heater outlet helium total pressure P_h6 is calculated, the value of the helium compressor pressure ratio Pr_S1 is adjusted so that the helium heater outlet helium total pressure P_h6 and the helium turbine inlet total pressure P0_h1 are balanced, the values of the air compressor inlet total temperature T0_a4 and the air compressor inlet total pressure P0_a4 are adjusted so that the air compressor inlet total temperature T0_a4 and the precooler outlet air side total temperature T_a3 are balanced, and the air compressor inlet total pressure T0_a4 and the precooler outlet air side total pressure T_a3 are balanced;

[0044] Step7:

[0045] According to the air compressor outlet flow qm_a5, air compressor outlet total pressure P_a5, air compressor outlet total temperature T_a5, and hydrogen turbine outlet flow qm_q4, hydrogen turbine outlet total pressure P_q4, hydrogen turbine outlet total temperature T_q4, the main combustion chamber outlet total temperature T_b, the main combustion chamber outlet total pressure P_b, and the main combustion chamber outlet flow qm_b are calculated; according to the main combustion chamber outlet total temperature T_b, the main combustion chamber outlet total pressure P_b, and the main combustion chamber outlet flow qm_b, the calculated value qm_e of the exhaust nozzle outlet flow, the exhaust nozzle outlet flow velocity V_e, and the exhaust nozzle outlet static pressure P_e are calculated, and according to the flow conservation, the theoretical value qm_t of the exhaust nozzle outlet flow is obtained, and the air compressor flow ratio function Z is adjusted to balance the calculated value qm_e of the exhaust nozzle outlet flow and the theoretical value qm_t of the exhaust nozzle outlet flow;

[0046] Step8:

[0047] According to the exhaust nozzle outlet flow calculated value qm_e, the exhaust nozzle outlet flow velocity V_e, the exhaust nozzle outlet static pressure P_e, and the atmospheric environment static pressure P_a0, the exhaust nozzle throat area A_e, the engine performance parameters are calculated, and the operation is ended.

[0048] Principles of the present application:

[0049] Figure 1 The pre-cooled combined engine closed helium cycle system model logical operation principle diagram includes 7 initial value input modules, 11 component calculation modules and 7 residual error judgment modules (balance equations), specifically including:

[0050] (1) The initial value input module is:

[0051] The air compressor flow ratio function Z (Step1 input, Step7 to determine whether to revalue);

[0052] The air compressor inlet total temperature T0_a4 and the air compressor inlet total pressure P0_a4 (Step1 input, Step6 to determine whether to revalue);

[0053] The helium turbine inlet flow qm0_h1 (Step2 input, Step3 to determine whether to revalue);

[0054] The helium turbine inlet total pressure P0_h1 (Step2 input, Step2 to determine whether to revalue);

[0055] The helium compressor pressure ratio Pr_S1 (Step4 input, Step6 to determine whether to revalue);

[0056] Hydrogen tank outlet total pressure P_q0 (Step 4 input, Step 5 judge whether to revalue)

[0057] Hydrogen turbine pressure ratio Pr_TH1 (Step 5 input, Step 5 judge whether to revalue).

[0058] (2) Component calculation module is:

[0059] Air compressor calculation module, helium turbine calculation module, hydrogen helium regenerator calculation module, helium compressor calculation module, hydrogen turbine calculation module, pre-cooler calculation module, gas generator calculation module, helium heater calculation module, main combustion chamber calculation module, nozzle calculation module and engine performance calculation module.

[0060] Residual error judgment module (balance equation) is:

[0061] 1. Nozzle outlet flow balance, that is

[0062] Delta_Z=(qm_t-qm_e) / qm_t

[0063] Located in Step7, if Delta_Z≮1×10 -3 , return to Step 1, revalue the air compressor flow ratio function Z;

[0064] 2. Air compressor inlet total temperature / total pressure balance, that is

[0065] Delta_T_a4=(T0_a4-T_a3) / T0_a4

[0066] Delta_P_a4=(P0_a4-P_a3) / P0_a4

[0067] Located in Step6, if Delta_T_a4≮1×10 -3 Or Delta_P_a4≮1×10 -3 , return to Step 1, revalue the air compressor inlet total temperature T0_a4, air compressor inlet total pressure P0_a4;

[0068] 3. Helium turbine and air compressor power balance, that is

[0069] Delta_N_T1=(N_T1-N_C) / N_C

[0070] Located in Step3, if Delta_N_T1≮1×10 -3 , return to Step 2, revalue the helium turbine inlet flow qm0_h1;

[0071] 4. Helium turbine inlet flow balance, that is

[0072] Delta_qm_h1=(qm_h1-qm0_h1) / qm0_h1

[0073] At Step2, if Delta_qm_h1>1x10 -3 , return to Step2, and revalue the helium turbine inlet total pressure P0_h1;

[0074] 5. Helium heater outlet total pressure closed loop, i.e.

[0075] Delta_P_h6=(P0_h1-P_h6) / P0_h1

[0076] At Step6, if Delta_P_h6>1x10 -3 , return to Step4, and revalue the helium compressor pressure ratio Pr_S1;

[0077] 6. Combustion chamber inlet medium total pressure balance, i.e.

[0078] Delta_P_pb=(P_a5-P_q4) / P_a5

[0079] At Step5, if Delta_P_pb>1x10 -3 , return to Step4, and revalue the hydrogen tank outlet hydrogen total pressure P_q0;

[0080] 7. Hydrogen turbine and helium compressor power balance, i.e.

[0081] Delta_N_TH1=(N_TH1-N_S1) / N_TH1

[0082] At Step5, if Delta_N_TH1>1x10 -3 , return to Step5, and revalue the hydrogen turbine drop pressure ratio Pr_TH1.

[0083] The application also has the following advantages:

[0084] When the turbine performance is calculated, the ground test data is used as the component characteristics of the components, wherein the component characteristics of the air compressor and the helium compressor adopt three-dimensional interpolation tables of flow rate, pressure ratio and efficiency; the component characteristics of the helium turbine and the hydrogen turbine are three-dimensional fitting functions of flow rate, pressure ratio and efficiency; the application solves the nonlinear equation set through the secant iteration method. The residual error judgment standard is:

[0085] 1) Tail nozzle outlet flow balance, i.e.

[0086] Delta_Z=(qm_t-qm_e) / qm_t

[0087] 2) Air compressor inlet total temperature / total pressure balance, i.e.

[0088] Delta_T_a4 = (T0_a4 - T_a3) / T0_a4

[0089] Delta_P_a4 = (P0_a4 - P_a3) / P0_a4

[0090] 3) Helium turbine and air compressor power balance, i.e.

[0091] Delta_N_T1 = (N_T1 - N_C) / N_C

[0092] 4) Helium turbine inlet flow balance, i.e.

[0093] Delta_qm_h1 = (qm_h1 - qm0_h1) / qm0_h1

[0094] 5) Helium heater outlet total pressure closed loop, i.e.

[0095] Delta_P_h6 = (P_h1 - P_h6) / P_h1

[0096] 6) Combustion chamber inlet medium total pressure balance, i.e.

[0097] Delta_P_pb = (P_a5 - P_q4) / P_a5

[0098] 7) Hydrogen turbine and helium compressor power balance, i.e.

[0099] Delta_N_TH1 = (N_TH1 - N_S1) / N_TH1

[0100] Error criteria are:

[0101] 1) Delta_Z > 1 x 10 -3

[0102] 2) Delta_T_a4 > 1 x 10 -3 or Delta_P_a4 > 1 x 10 -3

[0103] 3) Delta_N_T1 > 1 x 10 -3

[0104] 4) Delta_qm_h1 > 1 x 10 -3

[0105] 5) Delta_P_h6 > 1 x 10 -3

[0106] 6) Delta_P_pb > 1 x 10-3

[0107] 7) Delta_N_TH1 > 1 x 10 -3

[0108] Finally, the present application incorporates the helium cycle into the engine core by solving the balance equations of the precooler, compressor, and helium heater components as a whole.

[0109] The above only is the preferred embodiment of the present application, and is not used to limit the present application, for the person skilled in the art, the present application can have various changes and changes. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A pre-cooling combined engine closed helium cycle system and digital simulation modeling method, characterized in that, The method comprises the following steps: Step 1, according to the air compressor flow ratio function Z, the air compressor inlet total temperature T0_a4 and the air compressor inlet total pressure P0_a4, the air compressor outlet total temperature T_a5, the air compressor outlet total pressure P_a5, the air compressor power N_C and the air compressor outlet flow qm_a5 are calculated; Step 2, according to the helium turbine inlet flow qm0_h1 and the helium turbine inlet total pressure P0_h1, the calculated value qm_h1 of the helium turbine inlet flow is calculated; the value of the helium turbine inlet total pressure P0_h1 is adjusted to balance the helium turbine inlet flow qm0_h1 and the calculated value qm_h1 of the helium turbine inlet flow; the helium turbine power N_T1 is obtained according to the helium turbine inlet total pressure P0_h1 and the calculated value qm_h1 of the helium turbine inlet flow; Step 3, the value of the helium turbine inlet flow qm0_h1 is adjusted to balance the helium turbine power N_T1 and the air compressor power N_C; Step 4, according to the helium compressor pressure ratio Pr_S1 and the hydrogen storage tank outlet hydrogen total pressure P_q0, the hydrogen storage tank outlet hydrogen flow qm_q1 is calculated, and according to the hydrogen storage tank outlet hydrogen flow qm_q1, the hydrogen-helium recuperator outlet helium side flow qm_h3, the hydrogen-helium recuperator outlet helium side total pressure P_h3, the hydrogen-helium recuperator outlet helium side total temperature T_h3, the hydrogen-helium recuperator outlet hydrogen side flow qm_q2, the hydrogen-helium recuperator outlet hydrogen side total pressure P_q2 and the hydrogen-helium recuperator outlet hydrogen side total temperature T_q2 are calculated; according to the helium compressor pressure ratio Pr_S1, the hydrogen-helium recuperator outlet helium side flow qm_h3, the hydrogen-helium recuperator outlet helium side total pressure P_h3 and the hydrogen-helium recuperator outlet helium side total temperature T_h3, the helium compressor outlet total temperature T_h4, the helium compressor outlet total pressure P_h4, the helium compressor outlet flow qm_h4 and the helium compressor power N_S1 are calculated; Step 5, according to the hydrogen turbine drop pressure ratio Pr_TH1, the hydrogen-helium recuperator outlet hydrogen side flow qm_q2, the hydrogen-helium recuperator outlet hydrogen side total pressure P_q2 and the hydrogen-helium recuperator outlet hydrogen side total temperature T_q2, the hydrogen turbine outlet flow qm_q4, the hydrogen turbine outlet total temperature T_q4, the hydrogen turbine outlet total pressure P_q4 and the hydrogen turbine power N_TH1 are calculated; the value of the hydrogen turbine drop pressure ratio Pr_TH1 is adjusted to balance the hydrogen turbine power N_TH1 and the helium compressor power N_S1; the given hydrogen storage tank outlet hydrogen total pressure P_q0 is adjusted to balance the helium turbine outlet total pressure P_q4 and the air compressor outlet total pressure P_a5; Step 6, according to the pre-cooler inlet air side flow rate qm_a2, the pre-cooler inlet air side total pressure P_a2, the pre-cooler inlet air side total temperature T_a2, the helium pressure compressor outlet flow rate qm_h4, the helium pressure compressor outlet total pressure P_h4 and the helium pressure compressor outlet total temperature T_h4, the pre-cooler outlet helium side total pressure P_h5, the pre-cooler outlet helium side total temperature T_h5, the pre-cooler outlet air side total pressure P_a3 and the pre-cooler outlet air side total temperature T_a3 are calculated; according to the pre-cooler outlet helium side total pressure P_h5, the pre-cooler outlet helium side total temperature T_h5, the pre-cooler outlet air side total pressure P_a3 and the pre-cooler outlet air side total temperature T_a3, the helium heater outlet helium total pressure P_h6 is calculated; the value of the helium pressure compressor pressure ratio Pr_S1 is adjusted so that the helium heater outlet helium total pressure P_h6 and the helium turbine inlet total pressure P0_h1 reach equilibrium; the values of the air pressure compressor inlet total temperature T0_a4 and the air pressure compressor inlet total pressure P0_a4 are adjusted so that the air pressure compressor inlet total temperature T0_a4 and the pre-cooler outlet air side total temperature T_a3 reach equilibrium, and the air pressure compressor inlet total pressure P0_a4 and the pre-cooler outlet air side total pressure P_a3 reach equilibrium; Step 7, according to the air pressure compressor outlet flow rate qm_a5, the air pressure compressor outlet total pressure P_a5, the air pressure compressor outlet total temperature T_a5, the hydrogen turbine outlet flow rate qm_q4, the hydrogen turbine outlet total pressure P_q4 and the hydrogen turbine outlet total temperature T_q4, the main combustion chamber outlet total temperature T_b, the main combustion chamber outlet total pressure P_b and the main combustion chamber outlet flow rate qm_b are calculated; According to the main combustion chamber outlet total temperature T_b, the main combustion chamber outlet total pressure P_b and the main combustion chamber outlet flow rate qm_b, the calculated value of the tail nozzle outlet flow rate qm_e, the tail nozzle outlet flow velocity V_e and the tail nozzle outlet static pressure P_e are calculated; According to the flow conservation, the theoretical value of the tail nozzle outlet flow rate qm_t is obtained; the air pressure compressor flow ratio function Z is adjusted so that the calculated value of the tail nozzle outlet flow rate qm_e and the theoretical value of the tail nozzle outlet flow rate qm_t reach equilibrium; Step 8, according to the calculated value of the tail nozzle outlet flow rate qm_e, the tail nozzle outlet flow velocity V_e, the tail nozzle outlet static pressure P_e, the atmospheric environment static pressure P_a0 and the tail nozzle throat area A_e, the engine performance parameters are calculated.

2. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method according to claim 1, characterized in that, The value of the helium turbine inlet total pressure P0_h1 is adjusted so that the helium turbine inlet flow rate qm0_h1 and the calculated value of the helium turbine inlet flow rate qm_h1 reach equilibrium, and the specific method is as follows: Delta_qm_h1=(qm_h1-qm0_h1) / qm0_h1 If Delta_qm_h1 > 1 x 10 -3 Then return to revalue the helium turbine inlet total pressure P0_h1.

3. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method of claim 1, wherein, The value of the helium turbine inlet flow rate qm0_h1 is adjusted so that the helium turbine power N_T1 and the air pressure compressor power N_C reach equilibrium, and the specific method is as follows: Delta_N_T1=(N_T1-N_C) / N_C If Delta_N_T1 > 1 x 10 -3 Then return to revalue the helium turbine inlet flow qm0_h1.

4. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method of claim 1, wherein, The value of the hydrogen turbine drop pressure ratio Pr_TH1 is adjusted so that the hydrogen turbine power N_TH1 and the helium pressure compressor power N_S1 reach equilibrium, and the specific method is as follows: Delta_N_TH1 = (N_TH1 - N_S1) / N_TH1 If Delta_N_TH1 > 1 x 10 -3 Then return to revalue the hydrogen turbine pressure ratio Pr_TH1.

5. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method of claim 1, wherein, The given hydrogen storage tank outlet hydrogen total pressure P_q0 is adjusted to balance the helium turbine outlet total pressure P_q4 and the air compressor outlet total pressure P_a5, specifically as follows: Delta_P_pb = (P_a5 - P_q4) / P_a5 If Delta_P_pb> 1 x 10 -3 Then return to revalue the total hydrogen pressure P_q0 at the outlet of the hydrogen storage tank.

6. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method of claim 1, wherein, The value of the helium compressor pressure ratio Pr_S1 is adjusted to balance the helium heater outlet helium total pressure P_h6 and the helium turbine inlet total pressure P0_h1, specifically as follows: Delta_P_h6 = (P0_h1 - P_h6) / P0_h1 If Delta_P_h6> 1 x 10 -3 Then return to revalue the helium compressor pressure ratio Pr_S1.

7. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method of claim 1, wherein, The values of the air compressor inlet total temperature T0_a4 and the air compressor inlet total pressure P0_a4 are adjusted to balance the air compressor inlet total temperature T0_a4 and the precooler outlet air side total temperature T_a3, and to balance the air compressor inlet total pressure P0_a4 and the precooler outlet air side total pressure P_a3, specifically as follows: Delta_T_a4 = (T0_a4 - T_a3) / T0_a4 Delta_P_a4 = (P0_a4 - P_a3) / P0_a4 If Delta_T_a4 > 1 x 10 -3 or Delta_P_a4 > 1 x 10 -3 then return to revalue the air compressor inlet total temperature T0_a4 or the air compressor inlet total pressure P0_a4.

8. The pre-cooling combined engine closed helium cycle system and digital simulation modeling method of claim 1, wherein, The air compressor flow ratio function Z is adjusted to balance the calculated value qm_e of the nozzle outlet flow and the theoretical value qm_t of the nozzle outlet flow, specifically as follows: Delta_Z = (qm_t - qm_e) / qm_t If Delta_Z > 1 x 10 -3 Then re-evaluate the air compressor flow ratio function Z.

9. A pre-cooling combined engine closed helium cycle system based on the method of any one of claims 1-8, characterized in that, It comprises: An initial value input module for inputting the air compressor flow ratio function Z, the air compressor inlet total temperature T0_a4, the air compressor inlet total pressure P0_a4, the helium turbine inlet flow qm0_h1, the helium turbine inlet total pressure P0_h1, the helium compressor pressure ratio Pr_S1, the hydrogen storage tank outlet hydrogen total pressure P_q0, and the hydrogen turbine drop pressure ratio Pr_TH1; A component calculation module comprising an air compressor calculation module, a helium turbine calculation module, a hydrogen-helium recuperator calculation module, a helium compressor calculation module, a hydrogen turbine calculation module, a precooler calculation module, a gas generator calculation module, a helium heater calculation module, a main combustion chamber calculation module, a nozzle calculation module, and an engine performance calculation module; A residual error judgment module for judging whether the nozzle outlet flow is balanced, judging whether the air compressor inlet total temperature / total pressure is balanced, judging whether the helium turbine and air compressor power is balanced, judging whether the helium turbine inlet flow is balanced, judging whether the helium heater outlet total pressure is closed-loop, judging whether the combustion chamber inlet medium total pressure is balanced, and judging whether the hydrogen turbine and helium compressor power is balanced.

10. The pre-cooling combined engine closed helium cycle system of claim 9, wherein, The air compressor calculation module is configured to calculate the air compressor outlet total temperature T_a5, the air compressor outlet total pressure P_a5, the air compressor power N_C, and the air compressor outlet flow qm_a5. The helium turbine calculation module is configured to calculate the calculated value qm_h1 of the helium turbine inlet flow. The hydrogen-helium recuperator calculation module is configured to calculate the hydrogen turbine inlet flow qm_TH1, the hydrogen turbine inlet total pressure P0_TH1, the hydrogen turbine inlet temperature T0_TH1, the hydrogen turbine outlet flow qm_TH2, the hydrogen turbine outlet total pressure P_TH2, the hydrogen turbine outlet temperature T_TH2, the hydrogen turbine power N_TH2, the hydrogen turbine drop pressure ratio Pr_TH1, the hydrogen turbine drop temperature ratio T_TH1, and the hydrogen turbine drop pressure ratio Pr_TH1. The helium compressor calculation module is configured to calculate the helium compressor inlet flow qm_S1, the helium compressor inlet total pressure P0_S1, the helium compressor inlet temperature T0_S1, the helium compressor outlet flow qm_S2, the helium compressor outlet total pressure P_S2, the helium compressor outlet temperature T_S2, the helium compressor power N_S2, the helium compressor pressure ratio Pr_S1, the helium compressor drop pressure ratio Pr_S1, and the helium compressor drop temperature ratio T_S1. The helium heater calculation module is configured to calculate the helium heater inlet flow qm_H1, the helium heater inlet total pressure P0_H1, the helium heater inlet temperature T0_H1, the helium heater outlet flow qm_H2, the helium heater outlet total pressure P_H2, the helium heater outlet temperature T_H2, the helium heater power N_H2, the helium heater drop pressure ratio Pr_H1, the helium heater drop temperature ratio T_H1, and the helium heater drop pressure ratio Pr_H1. The main combustion chamber calculation module is configured to calculate the main combustion chamber inlet flow qm_C1, the main combustion chamber inlet total pressure P0_C1, the main combustion chamber inlet temperature T0_C1, the main combustion chamber outlet flow qm_C2, the main combustion chamber outlet total pressure P_C2, the main combustion chamber outlet temperature T_C2, the main combustion chamber power N_C2, the main combustion chamber drop pressure ratio Pr_C1, the main combustion chamber drop temperature ratio T_C1, and the main combustion chamber drop pressure ratio Pr_C1. The nozzle calculation module is configured to calculate the nozzle inlet flow qm_N1, the nozzle inlet total pressure P0_N1, the nozzle inlet temperature T0_N1, the nozzle outlet flow qm_N2, the nozzle outlet total pressure P_N2, the nozzle outlet temperature T_N2, the nozzle power N_N2, the nozzle drop pressure ratio Pr_N1, the nozzle drop temperature ratio T_N1, and the nozzle drop pressure ratio Pr_N1. The engine performance calculation module is configured to calculate the engine inlet flow qm_E1, the engine inlet total pressure P0_E1, the engine inlet temperature T0_E1, the engine outlet flow qm_E2, the engine outlet total pressure P_E2, the engine outlet temperature T_E2, the engine power N_E2, the engine drop pressure ratio Pr_E1, the engine drop temperature ratio T_E1, and the engine drop pressure ratio Pr_E1. A hydrogen-helium recuperator calculation module is configured to calculate a hydrogen-helium recuperator outlet helium side flow rate qm_h3, a hydrogen-helium recuperator outlet helium side total pressure P_h3, a hydrogen-helium recuperator outlet helium side total temperature T_h3, a hydrogen-helium recuperator outlet hydrogen side flow rate qm_q2, a hydrogen-helium recuperator outlet hydrogen side total pressure P_q2, and a hydrogen-helium recuperator outlet hydrogen side total temperature T_q2; A helium compressor calculation module is configured to calculate a helium compressor outlet total temperature T_h4, a helium compressor outlet total pressure P_h4, a helium compressor outlet flow rate qm_h4, and a helium compressor power N_S1; A hydrogen turbine calculation module is configured to calculate a hydrogen turbine outlet flow rate qm_q4, a hydrogen turbine outlet total temperature T_q4, a hydrogen turbine outlet total pressure P_q4, and a hydrogen turbine power N_TH1; A precooler calculation module is configured to calculate a precooler outlet helium side total pressure P_h5, a precooler outlet helium side total temperature T_h5, a precooler outlet air side total pressure P_a3, and a precooler outlet air side total temperature T_a3; A helium heater calculation module is configured to calculate a helium heater outlet helium total pressure P_h6; A main combustion chamber calculation module is configured to calculate a main combustion chamber outlet total temperature T_b, a main combustion chamber outlet total pressure P_b, and a main combustion chamber outlet flow rate qm_b; A nozzle calculation module and an engine performance calculation module are configured to calculate engine performance parameters.

Images