CN117076830B - Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium - Google Patents

Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium Download PDF

Info

Publication number
CN117076830B
CN117076830B CN202311323376.1A CN202311323376A CN117076830B CN 117076830 B CN117076830 B CN 117076830B CN 202311323376 A CN202311323376 A CN 202311323376A CN 117076830 B CN117076830 B CN 117076830B
Authority
CN
China
Prior art keywords
nitrogen
compressor
inlet
outlet
relation
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.)
Active
Application number
CN202311323376.1A
Other languages
Chinese (zh)
Other versions
CN117076830A (en
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.)
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Original Assignee
Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
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 Computational Aerodynamics Institute of China Aerodynamics Research and Development Center filed Critical Computational Aerodynamics Institute of China Aerodynamics Research and Development Center
Priority to CN202311323376.1A priority Critical patent/CN117076830B/en
Publication of CN117076830A publication Critical patent/CN117076830A/en
Application granted granted Critical
Publication of CN117076830B publication Critical patent/CN117076830B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Abstract

The invention discloses a method and a device for determining the power of a nitrogen axial flow compressor, electronic equipment and a medium, and is applied to the technical field of axial flow compressors. The method comprises the following steps: acquiring a compressor inlet parameter of an axial flow compressor; based on the compressor inlet parameters and a thermodynamic relation established in advance based on a nitrogen three-time state equation, calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor; based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor, calculating the relation by combining the power of the compressor, and calculating to obtain the power of the compressor; the method provided by the invention calculates the power of the compressor based on the real nitrogen condition, can improve the calculation accuracy, has simple and feasible calculation process, and can improve the calculation efficiency.

Description

Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium
Technical Field
The invention relates to the technical field of axial flow compressors, in particular to a method and a device for determining the power of a nitrogen axial flow compressor, electronic equipment and a medium.
Background
The low-temperature wind tunnel is a wind tunnel which increases the Reynolds number of an experiment by a method of reducing the temperature of a working medium, and uses nitrogen as an operation medium and increases the Reynolds number in a low-temperature and pressurizing mode. The total pressure range is 1.15-4.50 multiplied by 105Pa, the total temperature range is 110-313K, and the Mach number range is 0.15-1.30. In the design of a low-temperature wind tunnel, the magnitude of the driving power of the axial flow compressor is related to the selection of a driving motor and a frequency converter and the range and efficiency of the running envelope curve of the wind tunnel, and is an important design parameter. When nitrogen is operated in a low temperature pressurization state, the actual gas effect must be considered in the calculation of the compressor power because the actual gas effect is larger than the actual gas assumption. At present, two methods exist for calculating the power of an axial flow compressor, one method is to calculate according to ideal gas and finally multiply the power by a compression factor, and the method has larger error under the working condition of low temperature and high pressure; the other is to calculate according to the real gas completely, look-up table according to the Weili formula is needed for calculation, and the calculation process is complicated.
In view of this, how to provide a nitrogen axial compressor power determination method, device, system and computer readable storage medium with high computational efficiency and high accuracy is a problem that needs to be solved by those skilled in the art.
Disclosure of Invention
The embodiment of the invention aims to provide a method and a device for determining the power of a nitrogen axial flow compressor, electronic equipment and a computer readable storage medium, which can improve the calculation accuracy and the calculation efficiency in the use process.
In order to solve the technical problems, the embodiment of the invention provides a method for determining the power of a nitrogen axial flow compressor, which comprises the following steps:
acquiring a compressor inlet parameter of an axial flow compressor;
based on the compressor inlet parameters and a thermodynamic relation established in advance based on a nitrogen three-time state equation, calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor;
and calculating a relation by combining the compressor power based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor, and obtaining the compressor power by calculation.
Optionally, the thermodynamic relation established in advance based on the nitrogen cubic state equation comprises:
and obtaining an entropy relation and an enthalpy relation based on the nitrogen three-time state equation and the thermodynamic first law.
Optionally, the compressor inlet parameters include compressor inlet nitrogen pressure, compressor inlet nitrogen temperature, and compressor isentropic compression ratio.
Optionally, calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor based on the parameters at the inlet of the compressor and the thermodynamic relation established in advance based on the equation of state of the nitrogen three times comprises:
based on the nitrogen pressure at the inlet of the compressor and the nitrogen temperature at the inlet of the compressor, combining with a nitrogen three-time state equation, and calculating to obtain the nitrogen specific volume at the inlet of the compressor;
based on the nitrogen pressure at the inlet of the compressor, the nitrogen temperature at the inlet of the compressor and the nitrogen specific volume at the inlet of the compressor, the total enthalpy of the nitrogen at the inlet of the compressor is calculated by combining the enthalpy relation;
based on the specific volume of the nitrogen at the inlet of the compressor and the temperature of the nitrogen at the inlet of the compressor, calculating to obtain the nitrogen entropy at the inlet of the compressor by combining an entropy relation;
obtaining the nitrogen entropy of the outlet of the compressor based on the nitrogen entropy of the inlet of the compressor;
obtaining the nitrogen pressure at the outlet of the compressor based on the nitrogen pressure at the inlet of the compressor and the isentropic expansion pressure ratio;
based on the nitrogen pressure at the outlet of the compressor and the nitrogen entropy at the outlet of the compressor, combining a nitrogen cubic state equation and an entropy relation to obtain the nitrogen temperature at the outlet of the compressor and the nitrogen specific volume at the outlet of the compressor;
and calculating the total enthalpy of the nitrogen at the outlet of the compressor based on the pressure of the nitrogen at the outlet of the compressor, the temperature of the nitrogen at the outlet of the compressor and the specific volume of the nitrogen at the outlet of the compressor by combining the enthalpy relation.
Optionally, the nitrogen three-time state equation is:
wherein->Represents the nitrogen pressure, T represents the nitrogen temperature, +.>Represents the specific volume of nitrogen>Is a first nitrogen gas constant, +.>Is a second nitrogen gas constant, R is a third nitrogen gas constant, ">Represents the nitrogen temperature function @, @>，/>Represents the fourth nitrogen gas constant, +.>Represents the critical temperature of nitrogen, and n represents the index constant.
Optionally, the entropy relation is:
wherein->The reference temperature is indicated as such,a constant pressure specific heat coefficient representing the total gas, < ->Representing gas entropy;
the enthalpy relation is:
wherein, the method comprises the steps of, wherein,/>indicating the enthalpy of the gas.
Optionally, the compressor power calculation relationship is:
，/>represents the compressor power of the i-th compressor,/->Indicating mass flow,/->Indicating total enthalpy of the compressor inlet nitrogen, +.>Indicating the total enthalpy of the compressor outlet nitrogen.
The embodiment of the invention also provides a power determining device of the nitrogen axial flow compressor, which comprises the following components:
the acquisition module is used for acquiring the inlet parameters of the compressor of the axial flow compressor;
the first calculation module is used for calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor based on the parameters at the inlet of the compressor and the thermodynamic relation established in advance based on the three-time state equation of the nitrogen;
and the second calculation module is used for calculating the compressor power based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor and combining the compressor power calculation relational expression.
The embodiment of the invention also provides electronic equipment, which comprises:
a memory for storing a computer program;
and a processor for implementing the steps of the nitrogen axial compressor power determination method as described above when executing a computer program.
The embodiment of the invention also provides a computer readable storage medium, and a computer program is stored on the computer readable storage medium, and when the computer program is executed by a processor, the steps of the nitrogen axial flow compressor power determining method are realized.
The embodiment of the invention provides a method and a device for determining the power of a nitrogen axial flow compressor, electronic equipment and a computer readable storage medium, wherein the method comprises the following steps: acquiring a compressor inlet parameter of an axial flow compressor; based on the compressor inlet parameters and a thermodynamic relation established in advance based on a nitrogen three-time state equation, calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor; and calculating a relation by combining the compressor power based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor, and obtaining the compressor power by calculation.
Therefore, in the embodiment of the invention, the thermodynamic relation is established in advance according to the nitrogen three-time state equation, then when the compressor power is calculated, the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor are calculated according to the compressor inlet parameter of the axial flow compressor, the nitrogen three-time state equation and the thermodynamic relation, and then the compressor power is further calculated according to the compressor power calculation relation.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required in the prior art and the embodiments will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for determining the power of a nitrogen axial flow compressor according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of an integral path curve according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of the entropy and enthalpy curves according to the pressure and temperature according to the embodiment of the present invention;
FIG. 4 is a schematic diagram of another entropy and enthalpy versus pressure and temperature curve provided by an embodiment of the present invention;
FIG. 5 is a schematic diagram of a compressor power versus pressure curve according to an embodiment of the present invention;
FIG. 6 is a graph showing a variation of power versus pressure for a compressor according to an embodiment of the present invention;
FIG. 7 is a graph showing a variation of power of a compressor according to another embodiment of the present invention;
FIG. 8 is a schematic diagram of a power determining device for a nitrogen axial compressor according to an embodiment of the present invention;
fig. 9 is a schematic structural diagram of an electronic device according to an embodiment of the present invention;
fig. 10 is a schematic structural diagram of a computer readable storage medium according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a method and a device for determining the power of a nitrogen axial flow compressor, electronic equipment and a computer readable storage medium, which can improve the calculation accuracy and the calculation efficiency in the use process.
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Referring to fig. 1, fig. 1 is a flow chart of a method for determining power of a nitrogen axial flow compressor according to an embodiment of the invention. The method comprises the following steps:
s110: acquiring a compressor inlet parameter of an axial flow compressor;
in the embodiment of the invention, a thermodynamic relation is established in advance based on a nitrogen cubic state equation, a compressor power calculation relation is established, and the nitrogen cubic state equation and each established relation are stored in advance. When the power of the nitrogen axial flow compressor needs to be calculated, the compressor inlet parameter of the axial flow compressor can be obtained, wherein the compressor inlet parameter of the axial flow compressor can comprise the compressor inlet nitrogen pressure, the compressor inlet nitrogen temperature and the compressor isentropic compression ratio, and the compressor inlet nitrogen temperature and the compressor isentropic compression ratio are not limited to the parameters, but can also comprise other parameters in practical application.
S120: based on the compressor inlet parameters and a thermodynamic relation established in advance based on a nitrogen three-time state equation, calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor;
after the compressor inlet parameters of the axial flow compressor are obtained, the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor can be obtained by calculation according to the compressor inlet parameters and the thermodynamic relation established based on the nitrogen three-time state equation.
It should be noted that, in the embodiment of the invention, the thermodynamic relation is established according to the three-time state equation of nitrogen, so that the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor can be accurately calculated based on the related parameters of the real gas of the nitrogen.
S130: and calculating a relation by combining the compressor power based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor, and obtaining the compressor power by calculation.
Specifically, after the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor are obtained, a prestored compressor power calculation relation is called, and the compressor power of the nitrogen axial flow compressor is further calculated.
The axial flow compressor is mainly characterized in that the working medium is nitrogen, the working condition is that the temperature is about 90K under most conditions, and the pressure can be in the range of 1-10 atmospheres. And isentropic compression is carried out on the low-temperature nitrogen through an axial flow compressor, so that the output is high-pressure gas. In practical application, an Aungier-Redlich-Kwongz equation in three state equations can be selected, namely, the nitrogen three state equation in the embodiment of the invention can be an ARK state equation, and the ARK state equation is used for describing the state change of low-temperature nitrogen.
Specifically, the third state equation of nitrogen (i.e., ARK state equation) in the embodiment of the present invention isWherein->Represents the nitrogen pressure, T represents the nitrogen temperature, +.>The specific volume of nitrogen is indicated,is a first nitrogen gas constant, +.>Is the second nitrogen gasThe number and R are the third nitrogen gas constant, < ->Represents the nitrogen temperature function @, @>，/>Represents the fourth nitrogen gas constant, +.>Represents the critical temperature of nitrogen, and n represents the index constant.
In general, in the case of a conventional,，/>，/>，，/>，/>，/>，wherein->Represents the critical temperature of nitrogen, ">Represents the fourth nitrogen gas constant, +.>Represents the critical pressure of nitrogen, +.>Represents a centrifugation factor, n represents an exponential constant, +.>Represents the critical specific volume of nitrogen gas,represents the critical density of nitrogen.
And (3) applying isentropic and thermodynamic first law according to an ARK state equation to obtain a sound velocity expression:
wherein->Representing sound speed,/->Representing the partial derivative of enthalpy with respect to temperature, +.>Represents the partial derivative of the gas density with respect to the pressure, +.>Indicating gas density, +.>Represents the partial derivative of the gas density with respect to temperature, +.>Representing the partial derivative of enthalpy with respect to pressure.
Further, in the embodiment of the present invention, the thermodynamic relation established in advance based on the nitrogen cubic state equation may include:
and obtaining an entropy relation and an enthalpy relation based on the nitrogen three-time state equation and the thermodynamic first law.
Specifically, the entropy relation in the embodiment of the invention is as follows:
wherein->Indicating a reference temperature->A constant pressure specific heat coefficient representing the total gas, < ->Representing gas entropy;
the enthalpy relation is:
wherein->Indicating the enthalpy of the gas.
Specifically, the compressor power calculation relational expression in the embodiment of the invention is as follows:
，/>represents the compressor power of the i-th compressor,/->Indicating mass flow,/->Indicating total enthalpy of the compressor inlet nitrogen, +.>Indicating the total enthalpy of the compressor outlet nitrogen.
Specifically, according to the first law of thermodynamics of the reversible process, the internal energy is regarded as a function of specific volume and temperature, and the differential expression of entropy and internal energy is obtained through simple deduction after differentiation:
wherein s is entropy,>to fix the specific heat->For internal energy (I/O)>Representing the partial derivative of pressure with respect to temperature.
Specifically, the selectable integration route may refer to fig. 2, wherein Pref, tref, vref respectively represents pressure, temperature and specific volume of an arbitrary selected reference point, the integration path starts from point a to point D through B, C, the AB segment is an isothermal process, the specific volume of point B is infinity and can be regarded as complete gas, the BC segment is an isovolumetric process, the complete gas state equation can be applied in this segment, the CD segment is an isothermal process, the AB and CD segments integrate by using an ARK equation, that is, the partial derivatives are obtained based on the ARK state equation, and the internal energy relation and the entropy relation are obtained by integrating, and further the enthalpy relation is obtained by the specific formula as follows:
the entropy relation is:；
internal energy relation:；
enthalpy relation:。
the expression of entropy, enthalpy and internal energy of the low-temperature nitrogen at the specified temperature and specific volume can be obtained through the above process.
When the compressor power of the nitrogen axial flow compressor needs to be calculated, the compressor inlet parameters including the compressor inlet nitrogen pressure, the compressor inlet nitrogen temperature and the compressor isentropic compression ratio may be obtained.
The process of calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor based on the compressor inlet parameters and the thermodynamic relation established in advance based on the three-time state equation of the nitrogen specifically may include:
based on the nitrogen pressure at the inlet of the compressor and the nitrogen temperature at the inlet of the compressor, combining with a nitrogen three-time state equation, and calculating to obtain the nitrogen specific volume at the inlet of the compressor;
that is, knowing the compressor inlet nitrogen pressure and the compressor inlet nitrogen temperature, the ARK state direction is combinedThe specific volume of the nitrogen at the inlet of the compressor can be obtained.
Based on the nitrogen pressure at the inlet of the compressor, the nitrogen temperature at the inlet of the compressor and the nitrogen specific volume at the inlet of the compressor, the total enthalpy of the nitrogen at the inlet of the compressor is calculated by combining the enthalpy relation;
after the compressor inlet nitrogen pressure, the compressor inlet nitrogen temperature and the compressor inlet nitrogen specific volume are obtained, according to the enthalpy relation,
the total enthalpy of the nitrogen at the inlet of the compressor can be obtained.
Based on the specific volume of the nitrogen at the inlet of the compressor and the temperature of the nitrogen at the inlet of the compressor, calculating to obtain the nitrogen entropy at the inlet of the compressor by combining an entropy relation;
specifically, after the specific volume of the nitrogen at the inlet of the compressor is obtained, the nitrogen temperature at the inlet of the compressor is combined, and the entropy relation is basedThe compressor inlet nitrogen entropy can be further obtained.
Obtaining the nitrogen entropy of the outlet of the compressor based on the nitrogen entropy of the inlet of the compressor;
specifically, since the axial flow compressor is isentropic compression, the nitrogen entropy of the compressor outlet is equal to the nitrogen entropy of the compressor inlet, so that the nitrogen entropy of the compressor outlet can be obtained after the nitrogen entropy of the compressor inlet is calculated.
Obtaining the nitrogen pressure at the outlet of the compressor based on the nitrogen pressure at the inlet of the compressor and the isentropic expansion pressure ratio;
the isentropic expansion pressure ratio of the compressor can be obtained based on the parameters of the compressor, and the outlet total pressure of the compressor is determined by the isentropic expansion pressure ratio and the inlet total pressure, so that the outlet nitrogen pressure of the compressor can be obtained by calculation according to the inlet nitrogen pressure of the compressor and the isentropic expansion pressure ratio.
Based on the nitrogen pressure at the outlet of the compressor and the nitrogen entropy at the outlet of the compressor, combining a nitrogen cubic state equation and an entropy relation to obtain the nitrogen temperature at the outlet of the compressor and the nitrogen specific volume at the outlet of the compressor;
specifically, the nitrogen pressure at the outlet of the compressor and the nitrogen entropy at the outlet of the compressor are obtained, and the nitrogen three-time state equation is combinedThe degree of entropy of the liquid medium is determined,the nitrogen temperature at the outlet of the compressor can be obtained, and then the nitrogen temperature at the outlet of the compressor is obtained, and then the nitrogen three-time state equation is utilized by combining the nitrogen temperature at the outlet of the compressor and the nitrogen pressure at the outlet of the compressor>The specific volume of the nitrogen at the outlet of the compressor can be further calculated.
And calculating the total enthalpy of the nitrogen at the outlet of the compressor based on the pressure of the nitrogen at the outlet of the compressor, the temperature of the nitrogen at the outlet of the compressor and the specific volume of the nitrogen at the outlet of the compressor by combining the enthalpy relation.
Specifically, the compressor outlet nitrogen pressure, the compressor outlet nitrogen temperature and the compressor outlet nitrogen specific volume are obtained through the above processes, and by utilizing the enthalpy relation,the total enthalpy of the nitrogen at the outlet of the compressor can be obtained.
The relation can be calculated based on the compressor power, considering that the kinetic energy and potential energy change are small in the compression process of the working mediumCalculating to obtain the compressor power +.>Wherein->Is the test value.
Therefore, in the embodiment of the invention, the thermodynamic relation is established in advance according to the nitrogen three-time state equation, then when the compressor power is calculated, the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor are calculated according to the compressor inlet parameter of the axial flow compressor, the nitrogen three-time state equation and the thermodynamic relation, and then the compressor power is further calculated according to the compressor power calculation relation.
It should be further noted that, in order to further verify the accuracy of the above calculation method in the embodiment of the present invention, the following details are:
and selecting nitrogen with the temperature of 110K and the pressure of 1 bar-10 bar for calculation. Specifically, the nitrogen pressure can be divided into 10 equal parts at 1bar intervals, the pressure p and the temperature T are known, and an ARK state equation is utilizedCalculating the specific volume of nitrogen and solving a nonlinear equation by utilizing Newton tangent method>The specific volume v is obtained (specifically, the specific volume can be obtained after iteration for a plurality of times); and obtaining the entropy and the enthalpy of the gas under the specified pressure according to the entropy relation and the enthalpy relation. By comparison with standard librariesComparison of the data can verify the accuracy of the formulas and algorithms herein.
Specifically, the entropy, enthalpy change curves with pressure and temperature calculated by the method according to the embodiment of the present invention are compared with NIST (National Institute of Standards and Technology national standard and technology database) data curves, as shown in fig. 3 to 4. As can be seen from fig. 3 to 4, the entropy and enthalpy change curves are matched with the NIST data curves under different pressure and temperature conditions.
In addition, the operation parameters of the ETW (European Transonic Wind Tunnel, european transonic) wind tunnel can be referred to, 450kpa is selected, the total temperature 110K state is verified and calculated, and the compressor efficiency is 0.85. Knowing the compressor inlet nitrogen pressure, temperature, compressor pressure isentropic shrinkage, the enthalpy difference between the compressor outlet and inlet can be found by the following procedure. Specific:
the pressure and temperature of nitrogen at the inlet of the compressor are known, and ARK state equation is utilizedCalculating the specific volume of nitrogen at the inlet of the compressor; further according to the inlet temperature and specific volume of the compressor, the inlet entropy of the compressor is calculated by utilizing the entropy relational expression, then the isentropic expansion pressure ratio is multiplied by the inlet nitrogen pressure of the compressor according to the outlet nitrogen pressure of the compressor, the outlet nitrogen entropy of the compressor is equal to the inlet nitrogen entropy, an ARK state equation and an entropy relational expression are combined, the outlet nitrogen temperature of the compressor is calculated by applying an iterative Newton iteration method, after the outlet nitrogen temperature of the compressor is obtained, the outlet nitrogen temperature of the compressor and the outlet nitrogen pressure data of the compressor are substituted into the ARK state equation to be iterated to obtain the outlet nitrogen specific volume of the compressor, finally the temperature, the pressure, the specific volume and other data of the inlet/outlet nitrogen are substituted into the enthalpy relational expression to calculate the total enthalpy of the inlet and the outlet, and the known flow is substituted into the nitrogen flow according to the prior artThe compression power is obtained through formula calculation, so that the type selection of the driving motor and the frequency converter can be better carried out, and the range and the efficiency of the wind tunnel operation envelope curve can be determined.
Through verification, the calculation result obtained by adopting the method is shown in fig. 5 to 7, and the calculation result obtained by the method is compared with the calculation result obtained based on the approximation algorithm, so that the result is consistent with the trend of the approximation curve.
Therefore, the isentropic process of the nitrogen under low temperature and high pressure is described by introducing the ARK state equation, and the expression based on the entropy, enthalpy and internal energy of the nitrogen under the ARK state equation is applied to the calculation of the axial flow compressor under the nitrogen under the low temperature and high pressure to calculate the power of the compressor, so that the calculation is simple and convenient, and the calculation precision is high.
On the basis of the above embodiment, the embodiment of the present invention further provides a device for determining the power of a nitrogen axial flow compressor, with reference to fig. 8, where the device includes:
an acquisition module 11 for acquiring compressor inlet parameters of the axial flow compressor;
the first calculation module 12 is configured to calculate a total enthalpy of nitrogen at the inlet of the compressor and a total enthalpy of nitrogen at the outlet of the compressor based on the parameters at the inlet of the compressor and a thermodynamic relation established in advance based on a three-time equation of state of nitrogen;
the second calculation module 13 is configured to calculate the compressor power based on the total enthalpy of the compressor inlet nitrogen and the total enthalpy of the compressor outlet nitrogen, in combination with the compressor power calculation relation.
It should be noted that, the nitrogen axial compressor power determining device provided in the embodiment of the present invention has the same advantages as the nitrogen axial compressor power determining method provided in the above embodiment, and for the specific description of the nitrogen axial compressor power determining method related in the embodiment of the present invention, please refer to the above embodiment, and the invention is not repeated herein.
Fig. 9 is a block diagram of an electronic device according to an embodiment of the present invention, where, as shown in fig. 9, the electronic device includes: a memory 20 for storing a computer program;
a processor 21 for implementing the steps of the method for determining the power of the nitrogen axial compressor according to the above embodiment when executing a computer program.
The electronic device provided in this embodiment may include, but is not limited to, a smart phone, a tablet computer, a notebook computer, a desktop computer, or the like.
Processor 21 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 21 may be implemented in at least one hardware form of DSP (Digital Signal Processing ), FPGA (Field-Programmable Gate Array, field programmable gate array), PLA (Programmable Logic Array ). The processor 21 may also comprise a main processor, which is a processor for processing data in an awake state, also called CPU (Central Processing Unit ); a coprocessor is a low-power processor for processing data in a standby state. In some embodiments, the processor 21 may integrate a GPU (Graphics Processing Unit, image processor) for rendering and drawing of content required to be displayed by the display screen. In some embodiments, the processor 21 may also include an AI (Artificial Intelligence ) processor for processing computing operations related to machine learning.
Memory 20 may include one or more computer-readable storage media, which may be non-transitory. Memory 20 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In this embodiment, the memory 20 is at least used for storing a computer program 201, wherein the computer program, when loaded and executed by the processor 21, is capable of implementing the relevant steps of the nitrogen axial compressor power determination method disclosed in any of the foregoing embodiments. In addition, the resources stored in the memory 20 may further include an operating system 202, data 203, and the like, where the storage manner may be transient storage or permanent storage. The operating system 202 may include Windows, unix, linux, among others. The data 203 may include, but is not limited to, a set offset, etc.
In some embodiments, the electronic device may further include a display 22, an input-output interface 23, a communication interface 24, a power supply 25, and a communication bus 26.
Those skilled in the art will appreciate that the structure shown in fig. 9 is not limiting of the electronic device and may include more or fewer components than shown.
It will be appreciated that the nitrogen axial compressor power determination method of the above embodiments, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in part or in whole or in part in the form of a software product stored in a storage medium for performing all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random-access Memory (Random Access Memory, RAM), an electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, a magnetic disk, or an optical disk, etc. various media capable of storing program codes.
Based on this, as shown in fig. 10, the embodiment of the present invention further provides a computer readable storage medium, on which a computer program 31 is stored in the computer readable storage medium 30, and the computer program 31 implements the steps of the above-mentioned nitrogen axial compressor power determining method when executed by a processor.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
It should also be noted that in this specification, relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
Those of skill would further appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the various illustrative elements and steps are described above generally in terms of functionality in order to clearly illustrate the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. The software modules may be disposed in Random Access Memory (RAM), memory, read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (6)

1. A method for determining the power of a nitrogen axial compressor, comprising:
acquiring a compressor inlet parameter of an axial flow compressor;
based on the compressor inlet parameters and a thermodynamic relation established in advance based on a nitrogen three-time state equation, calculating total enthalpy of nitrogen at the inlet of the compressor and total enthalpy of nitrogen at the outlet of the compressor;
based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor, calculating a relational expression by combining the compressor power, and calculating to obtain the compressor power; wherein:
the thermodynamic relation established in advance based on the nitrogen cubic state equation comprises the following steps:
based on a nitrogen three-time state equation and a thermodynamic first law, obtaining an entropy relation and an enthalpy relation;
the compressor inlet parameters comprise compressor inlet nitrogen pressure, compressor inlet nitrogen temperature and compressor isentropic compression ratio;
the method for calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor based on the parameters at the inlet of the compressor and the thermodynamic relation established in advance based on the three-time state equation of the nitrogen comprises the following steps:
based on the nitrogen pressure at the inlet of the compressor and the nitrogen temperature at the inlet of the compressor, combining the nitrogen three-time state equation, and calculating to obtain the nitrogen specific volume at the inlet of the compressor;
based on the compressor inlet nitrogen pressure, the compressor inlet nitrogen temperature and the compressor inlet nitrogen specific volume, calculating to obtain total enthalpy of the compressor inlet nitrogen by combining the enthalpy relation;
based on the specific volume of the nitrogen at the inlet of the compressor and the temperature of the nitrogen at the inlet of the compressor, calculating to obtain the entropy of the nitrogen at the inlet of the compressor by combining the entropy relation;
obtaining the nitrogen entropy of the compressor outlet based on the nitrogen entropy of the compressor inlet;
obtaining the nitrogen pressure at the outlet of the compressor based on the nitrogen pressure at the inlet of the compressor and the isentropic expansion pressure ratio;
based on the nitrogen pressure at the outlet of the compressor and the nitrogen entropy at the outlet of the compressor, combining the nitrogen cubic state equation and the entropy relation to obtain the nitrogen temperature at the outlet of the compressor and the nitrogen specific volume at the outlet of the compressor;
based on the compressor outlet nitrogen pressure, the compressor outlet nitrogen temperature and the compressor outlet nitrogen specific volume, calculating by combining the enthalpy relation to obtain total enthalpy of the compressor outlet nitrogen;
the compressor power calculation relationship is:
，/>represents the compressor power of the i-th compressor,/->Indicating mass flow,/->Indicating total enthalpy of the compressor inlet nitrogen, +.>Indicating the total enthalpy of the compressor outlet nitrogen, +.>Represents the specific volume of nitrogen>Indicating the nitrogen pressure.
2. The nitrogen axial compressor power determination method of claim 1, wherein the nitrogen cubic equation of state is:
wherein->Represents the nitrogen pressure, T represents the nitrogen temperature, +.>Represents the specific volume of nitrogen>Is a first nitrogen gas constant, +.>Is a second nitrogen gas constant, R is a third nitrogen gas constant, ">Represents the nitrogen temperature function @, @>，/>Represents the fourth nitrogen gas constant, +.>Represents the critical temperature of nitrogen, and n represents the index constant.
3. The nitrogen axial compressor power determination method of claim 2, wherein the entropy relationship is:
wherein->Indicating a reference temperature->A constant pressure specific heat coefficient representing the total gas, < ->Representing gas entropy;
the enthalpy relation is:
wherein->Indicating the enthalpy of the gas.
4. A nitrogen axial compressor power determination apparatus, comprising:
the acquisition module is used for acquiring the inlet parameters of the compressor of the axial flow compressor;
the first calculation module is used for calculating the total enthalpy of nitrogen at the inlet of the compressor and the total enthalpy of nitrogen at the outlet of the compressor based on the inlet parameters of the compressor and a thermodynamic relation established in advance based on a nitrogen three-time state equation;
the second calculation module is used for calculating the compressor power based on the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor by combining a compressor power calculation relation; wherein:
the thermodynamic relation established in advance based on the nitrogen cubic state equation comprises the following steps:
based on a nitrogen three-time state equation and a thermodynamic first law, obtaining an entropy relation and an enthalpy relation;
the compressor inlet parameters comprise compressor inlet nitrogen pressure, compressor inlet nitrogen temperature and compressor isentropic compression ratio;
the method for calculating the total enthalpy of the nitrogen at the inlet of the compressor and the total enthalpy of the nitrogen at the outlet of the compressor based on the parameters at the inlet of the compressor and the thermodynamic relation established in advance based on the three-time state equation of the nitrogen comprises the following steps:
based on the nitrogen pressure at the inlet of the compressor and the nitrogen temperature at the inlet of the compressor, combining the nitrogen three-time state equation, and calculating to obtain the nitrogen specific volume at the inlet of the compressor;
based on the compressor inlet nitrogen pressure, the compressor inlet nitrogen temperature and the compressor inlet nitrogen specific volume, calculating to obtain total enthalpy of the compressor inlet nitrogen by combining the enthalpy relation;
based on the specific volume of the nitrogen at the inlet of the compressor and the temperature of the nitrogen at the inlet of the compressor, calculating to obtain the entropy of the nitrogen at the inlet of the compressor by combining the entropy relation;
obtaining the nitrogen entropy of the compressor outlet based on the nitrogen entropy of the compressor inlet;
obtaining the nitrogen pressure at the outlet of the compressor based on the nitrogen pressure at the inlet of the compressor and the isentropic expansion pressure ratio;
based on the nitrogen pressure at the outlet of the compressor and the nitrogen entropy at the outlet of the compressor, combining the nitrogen cubic state equation and the entropy relation to obtain the nitrogen temperature at the outlet of the compressor and the nitrogen specific volume at the outlet of the compressor;
based on the compressor outlet nitrogen pressure, the compressor outlet nitrogen temperature and the compressor outlet nitrogen specific volume, calculating by combining the enthalpy relation to obtain total enthalpy of the compressor outlet nitrogen;
the compressor power calculation relationship is:
，/>represents the compressor power of the i-th compressor,/->Indicating mass flow,/->Indicating total enthalpy of the compressor inlet nitrogen, +.>Indicating the total enthalpy of the compressor outlet nitrogen, +.>Represents the specific volume of nitrogen>Indicating the nitrogen pressure.
5. An electronic device, comprising:
a memory for storing a computer program;
a processor for implementing the steps of the nitrogen axial compressor power determination method according to any one of claims 1 to 3 when executing the computer program.
6. A computer readable storage medium, characterized in that it has stored thereon a computer program which, when executed by a processor, implements the steps of the nitrogen axial compressor power determination method according to any one of claims 1 to 3.
CN202311323376.1A 2023-10-13 2023-10-13 Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium Active CN117076830B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311323376.1A CN117076830B (en) 2023-10-13 2023-10-13 Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311323376.1A CN117076830B (en) 2023-10-13 2023-10-13 Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium

Publications (2)

Publication Number Publication Date
CN117076830A CN117076830A (en) 2023-11-17
CN117076830B true CN117076830B (en) 2024-03-12

Family

ID=88702848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311323376.1A Active CN117076830B (en) 2023-10-13 2023-10-13 Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium

Country Status (1)

Country Link
CN (1) CN117076830B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1190727A (en) * 1996-11-28 1998-08-19 气体产品与化学公司 Working utilizing airflow of high pressure nitrogen-enriched air
JP2000356158A (en) * 1999-06-15 2000-12-26 Nissan Motor Co Ltd Control device for engine
CN101029788A (en) * 2006-02-27 2007-09-05 三洋电机株式会社 Heat exchanger and refrigeration cycle device using the same
CN101196349A (en) * 2006-12-06 2008-06-11 何君 Air circulation machine driven by floating bearing type high-speed electric machine
WO2008115101A1 (en) * 2007-03-20 2008-09-25 Volvo Technology Corp. Method for operating a variable displacement compressor and system with a variable displacement compressor
CN102481824A (en) * 2009-08-25 2012-05-30 标致·雪铁龙汽车公司 Air-conditioning having an onboard compressor power estimator
WO2013057314A2 (en) * 2011-10-21 2013-04-25 Single Buoy Moorings Inc. Multi nitrogen expansion process for lng production
DE102013101028B3 (en) * 2013-02-01 2014-04-10 Hochschule Kempten Fast and accurate method for calculating the caloric quantities and the isentropic state of gases and mixtures of gases
CN111832189A (en) * 2020-07-24 2020-10-27 重庆通用工业(集团)有限责任公司 Centrifugal vapor compressor model selection method
CN112960978A (en) * 2021-02-02 2021-06-15 西北工业大学 A-site high-entropy perovskite oxide MeTiO3Thermoelectric ceramic and preparation method thereof
CN113137774A (en) * 2021-05-13 2021-07-20 上海海洋大学 Software-combined two-stage compression refrigeration cycle thermodynamic analysis method
CN115288819A (en) * 2022-08-11 2022-11-04 山东大学 Supercritical CO under oxygen-enriched combustion of coal 2 Recompression Brayton cycle coupling carbon capture novel combined system and simulation method
WO2023173624A1 (en) * 2022-03-14 2023-09-21 西安热工研究院有限公司 Soft measurement method for efficiency parameters of heavy-duty gas turbine key components

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1190727A (en) * 1996-11-28 1998-08-19 气体产品与化学公司 Working utilizing airflow of high pressure nitrogen-enriched air
JP2000356158A (en) * 1999-06-15 2000-12-26 Nissan Motor Co Ltd Control device for engine
CN101029788A (en) * 2006-02-27 2007-09-05 三洋电机株式会社 Heat exchanger and refrigeration cycle device using the same
CN101196349A (en) * 2006-12-06 2008-06-11 何君 Air circulation machine driven by floating bearing type high-speed electric machine
WO2008115101A1 (en) * 2007-03-20 2008-09-25 Volvo Technology Corp. Method for operating a variable displacement compressor and system with a variable displacement compressor
CN102481824A (en) * 2009-08-25 2012-05-30 标致·雪铁龙汽车公司 Air-conditioning having an onboard compressor power estimator
WO2013057314A2 (en) * 2011-10-21 2013-04-25 Single Buoy Moorings Inc. Multi nitrogen expansion process for lng production
DE102013101028B3 (en) * 2013-02-01 2014-04-10 Hochschule Kempten Fast and accurate method for calculating the caloric quantities and the isentropic state of gases and mixtures of gases
CN111832189A (en) * 2020-07-24 2020-10-27 重庆通用工业(集团)有限责任公司 Centrifugal vapor compressor model selection method
CN112960978A (en) * 2021-02-02 2021-06-15 西北工业大学 A-site high-entropy perovskite oxide MeTiO3Thermoelectric ceramic and preparation method thereof
CN113137774A (en) * 2021-05-13 2021-07-20 上海海洋大学 Software-combined two-stage compression refrigeration cycle thermodynamic analysis method
WO2023173624A1 (en) * 2022-03-14 2023-09-21 西安热工研究院有限公司 Soft measurement method for efficiency parameters of heavy-duty gas turbine key components
CN115288819A (en) * 2022-08-11 2022-11-04 山东大学 Supercritical CO under oxygen-enriched combustion of coal 2 Recompression Brayton cycle coupling carbon capture novel combined system and simulation method

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Influence of different equations of state on simulation results of supercritical CO2 centrifugal compressor;Yueming yang等;《Mechanical engineering science》;第3卷(第2期);25-33 *
低温增压真实气体效应下轴流压缩机功率计算;牟斌等;《第九届全国流体力学学术会议论文摘要集》;335 *
低速非定常流动的预处理计算方法;肖中云等;《计算流体力学研究进展——第十二届全国计算流体力学会议论文集》;104-109 *
牟斌等.低温增压真实气体效应下轴流压缩机功率计算.《第九届全国流体力学学术会议论文摘要集》.2016,335. *
真实气体效应对高压涡轮理论绝热功的影响;秦艳平等;《中国航天第三专业信息网第四十届技术交流会暨第四届空天动力联合会议论文集》;527-532 *
考虑低温增压真实气体效应的运输机气动特性数值模拟研究;江雄等;《空气动力学报》;第34卷(第5期);652-658 *

Also Published As

Publication number Publication date
CN117076830A (en) 2023-11-17

Similar Documents

Publication Publication Date Title
Sirakov et al. Evaluation of heat transfer effects on turbocharger performance
Hardcastle et al. Numerical modelling of the lobes of radio galaxies in cluster environments
Diez et al. Global models for moving contact lines
Dickmann et al. Unsteady flow in a turbocharger centrifugal compressor: three-dimensional computational fluid dynamics simulation and numerical and experimental analysis of impeller blade vibration
Baumgärtner et al. The effect of isentropic exponent on transonic turbine performance
Nishida et al. Performance improvement of a return channel in a multistage centrifugal compressor using multiobjective optimization
US7992388B2 (en) Method for estimating compressor output temperature for a two-stage turbocharger
Costabile et al. Δ h-Appell sequences and related interpolation problem
Cui Numerical study of unsteady flows in a scroll compressor
CN110781626A (en) Simulation method of finite difference multiple resolution trigonometric function WENO format
CN113868982A (en) Numerical simulation method and system for supercritical carbon dioxide radial flow type turbomachine
CN117076830B (en) Method and device for determining power of nitrogen axial flow compressor, electronic equipment and medium
Lou et al. The effect of gas models on compressor efficiency including uncertainty
Templalexis et al. Development of a two-dimensional streamline curvature code
Xu et al. A unified performance conversion method for similar compressors working with different gases based on polytropic analysis and deep-learning improvement
NOVOTNÝ et al. Singular limit of the equations of magnetohydrodynamics in the presence of strong stratification
Schein et al. Scroll compressor simulation model
Aymanns et al. Turbocharger efficiencies in pulsating exhaust gas flow
Lou et al. Reconstructing Compressor Nonuniform Circumferential Flow Field From Spatially Undersampled Data—Part 2: Practical Application for Experiments
Humpherys et al. Stability of viscous detonations for Majda’s model
Hamidreza et al. An investigation on turbocharger turbine performance parameters under inlet pulsating flow
Turkyilmazoglu Heat and mass transfer on the unsteady magnetohydrodynamic flow due to a porous rotating disk subject to a uniform outer radial flow
Lou et al. A new approach for centrifugal impeller preliminary design and aerothermal analysis
CN108167205A (en) LNG compressors are with pressure to start the method for determining
CN114139345A (en) Turboshaft engine output shaft power prediction method based on test data correlation

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
GR01 Patent grant
GR01 Patent grant