Thermal environment calculation method of power device and related equipment

Abstract

The invention discloses a thermal environment calculation method of a power device and related equipment, and the method comprises the steps: determining the operation condition of a target power device through the attribute information of the target power device, and determining the incoming flow condition of the target power device according to the operation condition; calculating a main flow parameter in a combustion chamber in the target power device based on a size parameter and an incoming flow condition in the attribute information; calculating the wall surface heat flux density of the target power device according to the mainstream parameters and the preset temperature by using an Ehmite reference enthalpy method; based on the wall surface heat flux density, the outlet temperature of a cooling working medium in the target power device is calculated through an energy conservation equation; and through a preset convective heat transfer equation, the wall surface temperature of the target power device is calculated according to the outlet temperature, temperature distribution data and heat flow distribution data of the target power device are determined according to the wall surface temperature and the preset temperature, and the temperature distribution data and the heat flow distribution data serve as the thermal environment of the target power device.

Description

Technical Field

[0001] The present invention relates to the technical field of thermal environment calculation, and in particular to a thermal environment calculation method for a power device and related equipment. Background Art

[0002] During operation, power units often need to withstand heat loads. As the power provided by the power unit increases, the heat load on the power unit increases. As the demand for faster speeds increases, the heat load on the power unit far exceeds the allowable temperature limit of the material structure. Therefore, thermal protection measures need to be taken for the power unit. Currently, the most effective thermal protection measure is active cooling technology, which uses a cooling medium to remove heat from the wall of the power unit, thereby maintaining the temperature within a safe range and ensuring long-term stable operation of the power unit.

[0003] However, the heat transfer process of the coolant in the wall cooling channel is strongly coupled with the heat transfer process of the high-temperature combustion gas inside the power unit. The distribution of the wall temperature and heat flux is subject to the dual constraints of the cooling condition and combustion. On the one hand, the heat flux of the power unit wall is required as the thermal boundary condition of the cooling channel to calculate the temperature distribution of the cooling channel wall. On the other hand, the wall temperature of the power unit needs to be known to calculate the wall heat flux.

[0004] In the prior art, the research on the thermal environment of the wall of the actively cooled power plant can be divided into two categories. One is to measure the wall temperature of the power plant under a given cooling condition through a test system, and the other is to establish an analytical model of the combustion chamber and cooling channel of the power plant through CFD (computational fluid dynamics) numerical simulation software for solution. Since the flow field in the combustion chamber of the power plant is very complex, there are complex physical phenomena such as shock waves and boundary layer separation, and the working fluid in the cooling channel may also have heat transfer phenomena such as phase change and cracking, the use of two-dimensional or three-dimensional simulation has a large amount of calculation, a long time consumption, and low efficiency, which is not conducive to the preliminary analysis and design of the active cooling structure of the power plant. The use of experimental means is costly, has a longer cycle, and has many uncontrollable factors. Therefore, the existing research methods for the thermal environment of the wall of the actively cooled power plant have limitations, which is not conducive to the research on the thermal environment of the wall of the actively cooled power plant. Summary of the invention

[0005] Based on this, it is necessary to propose a thermal environment calculation method and related equipment for a power plant in order to address the technical problems of existing research methods on the thermal environment of the wall of an actively cooled power plant, such as large amount of calculation, long time consumption and low efficiency.

[0006] In a first aspect, a method for calculating a thermal environment of a power device is provided, the method comprising:

[0007] Determining an operating condition of the target power device through attribute information of the target power device, and determining an incoming flow condition of the target power device according to the operating condition;

[0008] Establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation;

[0009] Calculating the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature by using the Eckert reference enthalpy method;

[0010] Based on the wall heat flux density, calculating the outlet temperature of the cooling medium in the target power device by using an energy conservation equation;

[0011] By presetting the convection heat transfer equation, the wall temperature of the target power device is calculated according to the outlet temperature, and the temperature distribution data and the heat flux distribution data of the target power device are determined according to the wall temperature and the preset temperature, and the temperature distribution data and the heat flux distribution data are used as the thermal environment of the target power device.

[0012] Optionally, the step of establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation includes:

[0013] Establishing a first control equation of the one-dimensional control equation in the combustion chamber of the target power device according to the size parameters in the attribute information by means of a one-dimensional Euler equation;

[0014] Obtaining the gas flow rate in the combustion chamber, and if the gas flow rate is greater than or equal to a threshold, calculating the wall friction coefficient in the first control equation by a preset formula, and substituting the wall friction coefficient into the first control equation to obtain a second control equation of the one-dimensional control equation;

[0015] Solving the second control equation and a preset mixing model for characterizing the turbulent mixing process in the combustion chamber according to the incoming flow condition to obtain mainstream parameters in the combustion chamber of the target power device;

[0016] If the gas flow rate is less than the threshold, the first control equation and the preset hybrid model are solved according to the incoming flow condition to obtain the mainstream parameters in the combustion chamber of the target power device.

[0017] Optionally, the step of establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation, further includes:

[0018] Establishing a third control equation of the one-dimensional control equation in the combustion chamber of the target power device according to the size parameters in the attribute information by an influence coefficient method;

[0019] Obtaining the gas flow rate in the combustion chamber, and if the gas flow rate is greater than or equal to a threshold, calculating the wall friction coefficient in the third control equation by a preset formula, and substituting the wall friction coefficient into the third control equation to obtain a fourth control equation of the one-dimensional control equation;

[0020] Solving the fourth control equation and a preset mixing model for characterizing the turbulent mixing process in the combustion chamber according to the incoming flow condition to obtain mainstream parameters in the combustion chamber of the target power device;

[0021] If the gas flow rate is less than the threshold, the third control equation and the preset hybrid model are solved according to the incoming flow condition to obtain the mainstream parameters in the combustion chamber of the target power device.

[0022] Optionally, before the step of calculating the wall temperature of the target power device according to the outlet temperature by presetting the convection heat transfer equation, the step further includes:

[0023] Acquire the type information of the cooling medium and the working condition data of the cooling medium, and select a convective heat transfer coefficient calculation formula corresponding to the type information and the working condition data from preset data according to the type information and the working condition data;

[0024] Calculating the convective heat transfer coefficient of the cooling medium based on the convective heat transfer coefficient calculation formula;

[0025] The initial convection heat transfer equation is updated according to the convection heat transfer coefficient to obtain the preset convection heat transfer equation.

[0026] Optionally, the step of determining the temperature distribution data and the heat flux distribution data of the target power device according to the wall temperature and the preset temperature includes:

[0027] Calculating a first difference between the wall temperature and the preset temperature, and if an absolute value of the first difference is greater than or equal to a preset error value, adjusting the preset temperature based on the absolute value of the first difference to obtain an updated temperature;

[0028] Calculating the updated wall heat flux density of the target power device according to the mainstream parameters and the updated temperature by using the Eckert reference enthalpy method;

[0029] Based on the updated wall heat flux density, calculating the updated outlet temperature of the cooling medium in the target power device through an energy conservation equation;

[0030] Calculating the updated wall temperature of the target power unit according to the updated outlet temperature by presetting a convection heat transfer equation;

[0031] Calculating a second difference between the updated wall temperature and the updated temperature, and if an absolute value of the second difference is greater than or equal to the preset error value, adjusting the updated temperature based on the absolute value of the second difference until the absolute value of the difference between the updated wall temperature and the updated temperature is less than the preset error value;

[0032] If the absolute value of the first difference is smaller than a preset error value, the temperature distribution data and the heat flux distribution data corresponding to the wall temperature are used as the temperature distribution data and the heat flux distribution data of the target power device.

[0033] Optionally, the method further includes:

[0034] Determine the maximum limit temperature of the wall of the target power device according to the material type information in the attribute information;

[0035] Acquiring operating parameters of the cooling medium in the target power device, and using the operating parameters as target inflow conditions;

[0036] Calculating target mainstream parameters in the combustion chamber of the target power device based on the size parameters in the attribute information and the target incoming flow conditions;

[0037] Calculating the target wall heat flux density of the target power device according to the target mainstream parameter and the preset temperature by using the Eckert reference enthalpy method;

[0038] Based on the target wall heat flux density, calculating the target outlet temperature of the cooling medium in the target power device by using an energy conservation equation;

[0039] By presetting a convection heat transfer equation, a target wall temperature of the target power device is calculated according to the target outlet temperature, and temperature distribution data of the target power device is determined according to the target wall temperature and the preset temperature;

[0040] Selecting a maximum temperature value from the temperature distribution data, and matching the maximum temperature value with the maximum limit temperature of the wall surface to obtain a matching result;

[0041] The cooling medium is evaluated according to the matching result.

[0042] Optionally, the method further includes:

[0043] Determine the maximum limit temperature of the wall of the target power device according to the material type information in the attribute information;

[0044] Taking the maximum limit temperature of the wall as the wall temperature of the target power device, and calculating a first outlet temperature of the first cooling medium according to the wall temperature and a first convection heat transfer equation corresponding to the first cooling medium;

[0045] Calculating the first wall heat flux density according to the first outlet temperature by using the energy conservation equation;

[0046] Calculating the first main flow parameter according to the first wall heat flux and the preset temperature by using the Eckert reference enthalpy method;

[0047] Calculate a first incoming flow condition based on the size parameter in the attribute information and the first mainstream parameter;

[0048] Taking the maximum limit temperature of the wall as the wall temperature of the target power device, and calculating the second outlet temperature of the second cooling medium according to the wall temperature and a second convection heat transfer equation corresponding to the second cooling medium;

[0049] Calculating the second wall heat flux density according to the second outlet temperature by using an energy conservation equation;

[0050] Calculating the second main flow parameter according to the second wall heat flux and the preset temperature by using the Eckert reference enthalpy method;

[0051] Calculate a second incoming flow condition based on the size parameter in the attribute information and the second main flow parameter;

[0052] Taking the medium flow rate of the cooling medium as a single variable, comparing the first medium flow rate in the first flow condition and the second medium flow rate in the second flow condition to obtain a comparison result;

[0053] The first cooling medium and the second cooling medium are evaluated according to the comparison result.

[0054] In a second aspect, the present application provides a thermal environment calculation device for a power device, the device comprising:

[0055] A data acquisition module, used to determine the operating condition of the target power device through the attribute information of the target power device, and determine the incoming flow condition of the target power device according to the operating condition;

[0056] A first calculation module, for establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation;

[0057] A second calculation module is used to calculate the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature by using the Eckert reference enthalpy method;

[0058] A third calculation module, configured to calculate the outlet temperature of the cooling medium in the target power device based on the wall heat flux density by using an energy conservation equation;

[0059] The fourth calculation module is used to calculate the wall temperature of the target power device according to the outlet temperature by using a preset convection heat transfer equation, and to determine the temperature distribution data and heat flux distribution data of the target power device according to the wall temperature and the preset temperature, and to use the temperature distribution data and heat flux distribution data as the thermal environment of the target power device.

[0060] In a third aspect, a computer device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the thermal environment calculation method of the power device as described above when executing the computer program.

[0061] In a fourth aspect, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the thermal environment calculation method of the power device are implemented.

[0062] The present application determines the operating condition of the target power device through the attribute information of the target power device, and determines the incoming flow condition of the target power device according to the operating condition; establishes a one-dimensional control equation in the combustion chamber of the target power device based on the dimensional parameters in the attribute information, and calculates the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation; uses the Eckert reference enthalpy method to calculate the wall heat flux density of the target power device according to the mainstream parameters and a preset temperature; based on the wall heat flux density, calculates the outlet temperature of the cooling medium in the target power device through the energy conservation equation; calculates the wall temperature of the target power device according to the outlet temperature through a preset convection heat transfer equation, and determines the temperature distribution data and heat flux distribution data of the target power device according to the wall temperature and the preset temperature, and uses the temperature distribution data and heat flux distribution data as the thermal environment of the target power device. By establishing a one-dimensional control equation to calculate the mainstream parameters in the combustion chamber of the target power device, the amount of calculation is reduced and the calculation efficiency is improved. The wall temperature and heat flux density of the power device under active cooling conditions can be quickly calculated, thereby improving the efficiency of the study of the thermal environment of the wall of the active cooling power device. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0064] in:

[0065] Figure 1 is a flow chart of a method for calculating the thermal environment of a power device in one embodiment;

[0066] Figure 2 is a flow chart for evaluating the cooling performance of a single cooling medium in one embodiment;

[0067] Figure 3 is a flow chart for evaluating the cooling performance of different cooling media in one embodiment;

[0068] Figure 4 is a structural block diagram of a thermal environment calculation device for a power device in one embodiment;

[0069] Figure 5 is a structural block diagram of a computer device in one embodiment;

[0070] Figure 6 It is a structural block diagram of a computer device in another embodiment. DETAILED DESCRIPTION

[0071] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0072] The present invention is described in detail below through specific embodiments.

[0073] See also Figure 1 As shown, Figure 1 A schematic flow chart of a method for calculating the thermal environment of a power device provided in an embodiment of the present invention includes the following steps:

[0074] S101, determining an operating condition of a target power device through attribute information of the target power device, and determining an incoming flow condition of the target power device according to the operating condition;

[0075] Exemplarily, the attribute information includes but is not limited to the type of the target power device and the dimensional parameters of the target power device. The type of the target power device can be understood as the type of the target power device, the dimensional parameters of the target power device can be understood as the geometric parameters of the target power device, and the incoming flow conditions can be understood as the airflow environment in which the target power device is located. The operating conditions of the target power device are interrelated with the incoming flow conditions, that is, when the target power device is in different operating conditions, it is in different airflow environments and corresponds to different incoming flow conditions.

[0076] Exemplarily, the process of determining the operating condition of the target power device through the attribute information of the target power device may be to determine the operating condition of the target power device through operating parameters in the attribute information of the target power device, such as speed, operating angle, etc.

[0077] S102, establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation;

[0078] In a possible implementation manner, the step of establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation includes:

[0079] Establishing a first control equation of the one-dimensional control equation in the combustion chamber of the target power device according to the size parameters in the attribute information by means of a one-dimensional Euler equation;

[0080] Obtaining the gas flow rate in the combustion chamber, and if the gas flow rate is greater than or equal to a threshold, calculating the wall friction coefficient in the first control equation by a preset formula, and substituting the wall friction coefficient into the first control equation to obtain a second control equation of the one-dimensional control equation;

[0081] Solving the second control equation and a preset mixing model for characterizing the turbulent mixing process in the combustion chamber according to the incoming flow condition to obtain mainstream parameters in the combustion chamber of the target power device;

[0082] If the gas flow rate is less than the threshold, the first control equation and the preset hybrid model are solved according to the incoming flow condition to obtain the mainstream parameters in the combustion chamber of the target power device.

[0083] Exemplarily, when establishing a one-dimensional control equation in the combustion chamber of the target power device, the one-dimensional control equation in the combustion chamber of the target power device can be established using the one-dimensional Euler equation. The Euler equation has high accuracy and comprehensiveness, and the efficiency of the mainstream parameter solution process is guaranteed by adopting the one-dimensional Euler equation.

[0084] Exemplarily, by considering the influence of the gas flow rate on the wall friction coefficient, the accuracy of mainstream parameter calculation is ensured.

[0085] Exemplarily, the turbulent mixing process in the combustion chamber is characterized by a preset mixing model, so that the present solution can be applied to different types of power devices.

[0086] S103, using the Eckert reference enthalpy method to calculate the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature;

[0087] Exemplarily, the Eckert reference enthalpy method uses the Reynolds analogy principle to perform qualitative temperature calculations, so a preset temperature needs to be set, and the enthalpy value at the wall temperature is calculated by the preset temperature, and then the wall heat flux density of the target power device is calculated based on the enthalpy value at the wall temperature.

[0088] Specifically,

[0089]

[0090] *Indicates reference temperature T * Variables under reference temperature T * The formula can be:

[0091]

[0092] Get, where h w is the enthalpy at the wall temperature, h aw is the adiabatic wall enthalpy, r is the airflow recovery coefficient, μ is the viscosity, and λ is the thermal conductivity. The wall heat flux density q at this temperature is obtained using the above formula w , Pr* is the Prandtl number, P e is the Peclet number, R is the gas constant, u e is the flow rate, r* is the reference temperature T * The airflow recovery coefficient under the condition of , μ is the viscosity of the airflow, λ* is the thermal conductivity, C p * is the specific heat at constant pressure, T is the real temperature, h is the enthalpy value at the real temperature, h * is the reference temperature T * enthalpy value under, Re is the Reynolds number, Re x is the Reynolds number at the location of the target power plant.

[0093] S104, calculating the outlet temperature of the cooling medium in the target power device based on the wall heat flux density by using an energy conservation equation;

[0094] Exemplarily, the high-temperature combustion gas generated in the combustion chamber of the target power device transfers heat to the wall of the target power device, and the cooling fluid in the wall channel absorbs the heat, thereby maintaining the wall temperature of the target power device within a safe range. That is, the cooling fluid in the wall channel exchanges energy with the wall of the target power device, and then the outlet temperature of the cooling fluid in the target power device is calculated by the energy conservation equation and the wall heat flux density.

[0095] Exemplarily, the wall heat flux density is: the amount of heat passing through the cross-sectional area of ​​the wall per unit time, that is, the heat passing through the cross-sectional area of ​​the wall can be understood as the heat generated by energy exchange between the cooling fluid in the wall channel and the wall of the target power device. It can be understood that the heat generated by energy exchange between the cooling fluid in the wall channel and the wall of the target power device is absorbed by the cooling fluid, and then the outlet temperature of the cooling fluid can be calculated.

[0096] S105. By presetting a convection heat transfer equation, the wall temperature of the target power device is calculated according to the outlet temperature, and the temperature distribution data and the heat flux distribution data of the target power device are determined according to the wall temperature and the preset temperature, and the temperature distribution data and the heat flux distribution data are used as the thermal environment of the target power device.

[0097] In a possible implementation manner, before the step of calculating the wall temperature of the target power device according to the outlet temperature by presetting the convection heat transfer equation, the step further includes:

[0098] Acquire the type information of the cooling medium and the working condition data of the cooling medium, and select a convective heat transfer coefficient calculation formula corresponding to the type information and the working condition data from preset data according to the type information and the working condition data;

[0099] Calculating the convective heat transfer coefficient of the cooling medium based on the convective heat transfer coefficient calculation formula;

[0100] The initial convection heat transfer equation is updated according to the convection heat transfer coefficient to obtain the preset convection heat transfer equation.

[0101] For example, for different cooling media, the calculation formulas for the convective heat transfer coefficient h are different, and the calculation formulas for the convective heat transfer coefficient of the same cooling media under different working conditions will also be different. Therefore, it is necessary to select the convective heat transfer coefficient calculation formula corresponding to the type information and the working condition data according to the type of cooling media and the working condition data of the cooling media to calculate the convective heat transfer coefficient of the cooling media.

[0102] Exemplarily, in a scenario where the type of cooling medium and / or the operating data of the cooling medium cannot be determined, the Gnielinski formula is selected:

[0103]

[0104] f=[1.82log(Re)-1.64] -2

[0105] Calculate the convective heat transfer coefficient of the cooling medium, where N u is the Nusselt number, Re is the Reynolds number, P r is the Prandtl number, d is the equivalent diameter, l is the length, f is the calculation parameter, and log() represents the logarithmic function. The Gnielinski formula has the advantage of high accuracy and can ensure the reliability of the calculation results.

[0106] The operating condition of the target power device is determined through the attribute information of the target power device, and the incoming flow condition of the target power device is determined according to the operating condition; a one-dimensional control equation in a combustion chamber in the target power device is established based on the size parameters in the attribute information, and the mainstream parameters in the combustion chamber in the target power device are calculated through the incoming flow condition and the one-dimensional control equation; the wall heat flux density of the target power device is calculated according to the mainstream parameters and a preset temperature by using the Eckert reference enthalpy method; based on the wall heat flux density, the outlet temperature of the cooling medium in the target power device is calculated by the energy conservation equation; the wall temperature of the target power device is calculated according to the outlet temperature by using a preset convection heat transfer equation, and the temperature distribution data and heat flux distribution data of the target power device are determined according to the wall temperature and the preset temperature, and the temperature distribution data and heat flux distribution data are used as the thermal environment of the target power device. By establishing a one-dimensional control equation to calculate the mainstream parameters in the combustion chamber of the target power device, the amount of calculation is reduced and the calculation efficiency is improved. The wall temperature and heat flux density of the power device under active cooling conditions can be quickly calculated, thereby improving the efficiency of the study of the thermal environment of the wall of the active cooling power device.

[0107] In a possible implementation manner, the step of establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation, further includes:

[0108] Establishing a third control equation of the one-dimensional control equation in the combustion chamber of the target power device according to the size parameters in the attribute information by an influence coefficient method;

[0109] Obtaining the gas flow rate in the combustion chamber, and if the gas flow rate is greater than or equal to a threshold, calculating the wall friction coefficient in the third control equation by a preset formula, and substituting the wall friction coefficient into the third control equation to obtain a fourth control equation of the one-dimensional control equation;

[0110] Solving the fourth control equation and a preset mixing model for characterizing the turbulent mixing process in the combustion chamber according to the incoming flow condition to obtain mainstream parameters in the combustion chamber of the target power device;

[0111] If the gas flow rate is less than the threshold, the third control equation and the preset hybrid model are solved according to the incoming flow condition to obtain the mainstream parameters in the combustion chamber of the target power device.

[0112] Exemplarily, the influence coefficient method is used to consider factors such as combustion, friction, and cross-sectional changes, and the basic control equation (the third control equation of the one-dimensional control equation) is constructed according to the size parameters in the attribute information:

[0113]

[0114] Where dA is the cross-sectional change, C f is the wall friction coefficient, T t is the total temperature; γ is the specific heat ratio; x is the length of the power unit component, Ma is the Mach number, A is the cross-sectional area of ​​the combustion chamber, D H is the equivalent diameter, T t is the total temperature of the incoming flow, is the fuel flow rate, γ is the specific heat ratio, y is the ratio of the streamwise velocity component of the injected fluid to the mainstream velocity, d represents the differential calculation, p is the pressure, and a, b, and c are preset parameters used to associate the above formulas.

[0115] When the gas velocity in the combustion chamber is high, such as supersonic combustion, the wall friction coefficient C f is large and cannot be ignored and can be solved by a preset formula.

[0116]

[0117] Calculate the wall friction coefficient C in the third control equation f , lg() is the log() function with base 10.

[0118] Considering that the actual turbulent mixing process cannot be simulated in one-dimensional calculation, a mixing model can be used to describe the preset mixing model. For example, for a slope nozzle, the mathematical expression of the preset mixing model is:

[0119]

[0120] Among them, η m is the mixing efficiency, exp() is the natural exponential function, x represents the length of different structures in the slope nozzle, A1 is a constant, usually ranging from 1 to 5, L m Indicates the mixing length, which can be selected according to the specific power unit configuration.

[0121] When the operating conditions and geometric dimensions of the power unit are known, the above differential equations are solved according to the incoming flow conditions to calculate the changes in the temperature, pressure and velocity of the fuel gas inside the power unit along the flow direction (mainstream parameters).

[0122] In a possible implementation manner, the step of determining the temperature distribution data and the heat flux distribution data of the target power device according to the wall temperature and the preset temperature includes:

[0123] Calculating a first difference between the wall temperature and the preset temperature, and if an absolute value of the first difference is greater than or equal to a preset error value, adjusting the preset temperature based on the absolute value of the first difference to obtain an updated temperature;

[0124] Calculating the updated wall heat flux density of the target power device according to the mainstream parameters and the updated temperature by using the Eckert reference enthalpy method;

[0125] Based on the updated wall heat flux density, calculating the updated outlet temperature of the cooling medium in the target power device through an energy conservation equation;

[0126] Calculating the updated wall temperature of the target power unit according to the updated outlet temperature by presetting a convection heat transfer equation;

[0127] Calculating a second difference between the updated wall temperature and the updated temperature, and if an absolute value of the second difference is greater than or equal to the preset error value, adjusting the updated temperature based on the absolute value of the second difference until the absolute value of the difference between the updated wall temperature and the updated temperature is less than the preset error value;

[0128] If the absolute value of the first difference is smaller than a preset error value, the temperature distribution data and the heat flux distribution data corresponding to the wall temperature are used as the temperature distribution data and the heat flux distribution data of the target power device.

[0129] For example, the wall heat flux q is obtained w After that, the outlet temperature T of the cooling medium can be obtained by using the energy conservation equation. c,2 .

[0130] Q w =m c C P,c (T c,1 -T c,2 )

[0131] Where m c is the cooling medium flow rate, T c,1 is the initial temperature of the cooling medium, C p,c Q is the specific heat capacity of the cooling medium, which is updated according to the temperature change. w is the total heat, which is equal to the product of heat flux and heat transfer area. For example, when the channel is rectangular, the total heat can be expressed as:

[0132] Q w =q w WL

[0133] Wherein, W is the length of the channel rectangle, and L is the width of the channel rectangle.

[0134] After obtaining the cooling medium temperature, the wall temperature T can be obtained according to the convection heat transfer equation between the wall and the medium. w2 :

[0135] Qw=A1h c (T w2 -T c )+η f A f h c (T w2 -T c )

[0136] Where T c is the average temperature of the cooling medium, A r , A f is the area of ​​the rib root and the rib surface when the rib effect is considered, η is the fin efficiency, and the fin efficiency is calculated by selecting the calculation formula corresponding to the channel shape according to the different channel shapes, and h c is the convective heat transfer coefficient of the cooling medium.

[0137] Exemplarily, the preset temperature is recorded as T W1 , and the cooling channel wall temperature T W2 After that, it should be compared with the assumed power plant wall temperature T W1 If the two are different, it is necessary to re-assume the power plant wall temperature T W1 , continue to repeat the steps of using the Eckert reference enthalpy method to calculate the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature; based on the wall heat flux density, calculate the outlet temperature of the cooling medium in the target power device through the energy conservation equation, and continuously iterate until the calculated wall temperature satisfies the error. At this time, the temperature distribution and heat flux distribution obtained are the thermal environment of the power device under active cooling conditions.

[0138] In a possible implementation, the method further includes:

[0139] Determine the maximum limit temperature of the wall of the target power device according to the material type information in the attribute information;

[0140] Acquiring operating parameters of the cooling medium in the target power device, and using the operating parameters as target inflow conditions;

[0141] Calculating target mainstream parameters in the combustion chamber of the target power device based on the size parameters in the attribute information and the target incoming flow conditions;

[0142] Calculating the target wall heat flux density of the target power device according to the target mainstream parameter and the preset temperature by using the Eckert reference enthalpy method;

[0143] Based on the target wall heat flux density, calculating the target outlet temperature of the cooling medium in the target power device by using an energy conservation equation;

[0144] By presetting a convection heat transfer equation, a target wall temperature of the target power device is calculated according to the target outlet temperature, and temperature distribution data of the target power device is determined according to the target wall temperature and the preset temperature;

[0145] Selecting a maximum temperature value from the temperature distribution data, and matching the maximum temperature value with the maximum limit temperature of the wall surface to obtain a matching result;

[0146] The cooling medium is evaluated according to the matching result.

[0147] For example, when evaluating the cooling effect of any cooling medium, considering that the primary goal of active cooling is to protect the wall of the power unit from being damaged by high temperature, the maximum wall temperature is selected as the main evaluation indicator of cooling performance. For the evaluation of cooling performance between different working fluids, considering the sensitivity of the aircraft to the overall mass, the mass flow rate of the cooling working fluid is selected as the main evaluation indicator. In addition, depending on the purpose, other evaluation indicators can also be selected, such as the maximum flow rate of the cooling working fluid, the average wall temperature of the power unit, or the maximum heat flux density of the wall.

[0148] Evaluation of the cooling performance of any working fluid, such as Figure 2As shown, the flow rate of the working fluid is first determined, and the initial parameters of the working fluid at the inlet of the cooling channel are selected, mainly including the inlet pressure and temperature, etc. (equivalent to taking the working parameters as the target inflow conditions, and the working parameters include the initial parameters), and the maximum wall temperature Tmax is calculated and solved (equivalent to calculating the target mainstream parameters in the combustion chamber of the target power device based on the size parameters in the attribute information and the target inflow conditions; using the Eckert reference enthalpy method, the target wall heat flux density of the target power device is calculated according to the target mainstream parameters and the preset temperature; based on the target wall heat flux density, the target outlet temperature of the cooling working fluid in the target power device is calculated by the energy conservation equation; by the preset convection heat transfer equation, the target wall temperature of the target power device is calculated according to the target outlet temperature, and the temperature distribution data of the target power device is determined according to the target wall temperature and the preset temperature; the maximum temperature value is selected in the temperature distribution data). Then, the allowable maximum limit temperature Tlimit of the wall is determined according to the type of the selected power unit wall material (equivalent to determining the maximum limit temperature of the wall of the target power unit according to the material type information in the attribute information), and the calculated maximum wall temperature Tmax is compared with the allowable limit temperature Tlimit of the wall. If the calculated maximum wall temperature Tmax is greater than the allowable limit temperature Tlimit of the wall, the cooling fluid cannot meet the cooling requirements of the power unit, and the fluid flow rate or inlet parameters need to be changed and redesigned. If the calculated maximum wall temperature Tmax is less than the allowable limit temperature Tlimi of the wall, the surface cooling fluid can meet the cooling requirements (equivalent to matching the maximum temperature value with the maximum limit temperature of the wall to obtain a matching result; and completing the evaluation of the cooling fluid based on the matching result).

[0149] In a possible implementation, the method further includes:

[0150] Determine the maximum limit temperature of the wall of the target power device according to the material type information in the attribute information;

[0151] Taking the maximum limit temperature of the wall as the wall temperature of the target power device, and calculating a first outlet temperature of the first cooling medium according to the wall temperature and a first convection heat transfer equation corresponding to the first cooling medium;

[0152] Calculating the first wall heat flux density according to the first outlet temperature by using the energy conservation equation;

[0153] Calculating the first main flow parameter according to the first wall heat flux and the preset temperature by using the Eckert reference enthalpy method;

[0154] Calculate a first incoming flow condition based on the size parameter in the attribute information and the first mainstream parameter;

[0155] Taking the maximum limit temperature of the wall as the wall temperature of the target power device, and calculating the second outlet temperature of the second cooling medium according to the wall temperature and a second convection heat transfer equation corresponding to the second cooling medium;

[0156] Calculating the second wall heat flux density according to the second outlet temperature by using an energy conservation equation;

[0157] Calculating the second main flow parameter according to the second wall heat flux and the preset temperature by using the Eckert reference enthalpy method;

[0158] Calculate a second incoming flow condition based on the size parameter in the attribute information and the second main flow parameter;

[0159] Taking the medium flow rate of the cooling medium as a single variable, comparing the first medium flow rate in the first flow condition and the second medium flow rate in the second flow condition to obtain a comparison result;

[0160] The first cooling medium and the second cooling medium are evaluated according to the comparison result.

[0161] For example, Figure 3 As shown, select working fluid A (first cooling medium), determine the working fluid inlet parameters of the cooling channel, and calculate the mass flow rate m of working fluid A required when the maximum wall temperature of the power unit is Tlimit. A Select working fluid B (second cooling medium), use the same cooling channel working fluid inlet parameters, and calculate the mass flow rate m of working fluid B required when the maximum wall temperature of the power unit is Tlimit B . Compare the calculated mass flow rate m of working fluid A A and the mass flow rate m of working fluid B B Size, if the mass flow rate of working fluid A is m A Less than the mass flow rate m of working fluid B B , it means that the cooling effect of working fluid A is better than that of working fluid B. If the mass flow rate of working fluid A is m A Greater than the mass flow rate m of working fluid B B , it means that the cooling effect of working fluid B is better than that of working fluid A.

[0162] Second, as Figure 4 As shown, the present application provides a thermal environment calculation device for a power device, the device comprising:

[0163] The data acquisition module 201 is used to determine the operating condition of the target power device through the property information of the target power device, and determine the incoming flow condition of the target power device according to the operating condition;

[0164] A first calculation module 202 is used to establish a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculate the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation;

[0165] A second calculation module 203 is used to calculate the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature by using the Eckert reference enthalpy method;

[0166] A third calculation module 204 is used to calculate the outlet temperature of the cooling medium in the target power device based on the wall heat flux density through the energy conservation equation;

[0167] The fourth calculation module 205 is used to calculate the wall temperature of the target power device according to the outlet temperature by using a preset convection heat transfer equation, and to determine the temperature distribution data and heat flux distribution data of the target power device according to the wall temperature and the preset temperature, and to use the temperature distribution data and heat flux distribution data as the thermal environment of the target power device.

[0168] The operating condition of the target power device is determined through the attribute information of the target power device, and the incoming flow condition of the target power device is determined according to the operating condition; a one-dimensional control equation in a combustion chamber in the target power device is established based on the size parameters in the attribute information, and the mainstream parameters in the combustion chamber in the target power device are calculated through the incoming flow condition and the one-dimensional control equation; the wall heat flux density of the target power device is calculated according to the mainstream parameters and a preset temperature by using the Eckert reference enthalpy method; based on the wall heat flux density, the outlet temperature of the cooling medium in the target power device is calculated by the energy conservation equation; the wall temperature of the target power device is calculated according to the outlet temperature by using a preset convection heat transfer equation, and the temperature distribution data and heat flux distribution data of the target power device are determined according to the wall temperature and the preset temperature, and the temperature distribution data and heat flux distribution data are used as the thermal environment of the target power device. By establishing a one-dimensional control equation to calculate the mainstream parameters in the combustion chamber of the target power device, the amount of calculation is reduced and the calculation efficiency is improved. The wall temperature and heat flux density of the power device under active cooling conditions can be quickly calculated, thereby improving the efficiency of the study of the thermal environment of the wall of the active cooling power device.

[0169] In one embodiment, a computer device is provided. The computer device may be a server, and its internal structure diagram may be as follows: Figure 5As shown. The computer device includes a processor, a memory, a network interface and a database connected via a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile and / or volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external client via a network connection. When the computer program is executed by the processor, it realizes the functions or steps on the server side of a thermal environment calculation method for a power device.

[0170] In one embodiment, a computer device is provided. The computer device may be a client, and its internal structure diagram may be as follows: Figure 5 As shown. The computer device includes a processor, a memory, a network interface, a display screen and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external server through a network connection. When the computer program is executed by the processor, the functions or steps on the client side of a thermal environment calculation method of a power device are implemented.

[0171] In one embodiment, an electronic device is proposed, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the following steps when executing the computer program: determining the operating condition of the target power device through attribute information of the target power device, and determining the incoming flow condition of the target power device according to the operating condition; establishing a one-dimensional control equation in a combustion chamber in the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber in the target power device through the incoming flow condition and the one-dimensional control equation; calculating the wall heat flux density of the target power device according to the mainstream parameters and a preset temperature by using the Eckert reference enthalpy method; calculating the outlet temperature of the cooling medium in the target power device by using the energy conservation equation based on the wall heat flux density; calculating the wall temperature of the target power device according to the outlet temperature by using a preset convection heat transfer equation, and determining the temperature distribution data and heat flux distribution data of the target power device according to the wall temperature and the preset temperature, and using the temperature distribution data and heat flux distribution data as the thermal environment of the target power device.

[0172] In one embodiment, a computer-readable storage medium is proposed, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the following steps are implemented: determining the operating condition of the target power device through attribute information of the target power device, and determining the incoming flow condition of the target power device according to the operating condition; establishing a one-dimensional control equation in a combustion chamber in the target power device based on the dimensional parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber in the target power device through the incoming flow condition and the one-dimensional control equation; calculating the wall heat flux density of the target power device according to the mainstream parameters and a preset temperature using the Eckert reference enthalpy method; calculating the outlet temperature of the cooling medium in the target power device through the energy conservation equation based on the wall heat flux density; calculating the wall temperature of the target power device according to the outlet temperature through a preset convection heat transfer equation, and determining the temperature distribution data and heat flux distribution data of the target power device according to the wall temperature and the preset temperature, and using the temperature distribution data and heat flux distribution data as the thermal environment of the target power device.

[0173] It should be noted that the above functions or steps that can be implemented by the computer-readable storage medium or computer device can refer to the relevant descriptions on the server side and the client side in the aforementioned method embodiment. To avoid repetition, they will not be described one by one here.

[0174] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0175] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example. In actual applications, the above-mentioned functions can be distributed and completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above.

[0176] The embodiments described above are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A method for calculating the thermal environment of a power plant, characterized in that: The method comprises: Determining an operating condition of the target power device through attribute information of the target power device, and determining an incoming flow condition of the target power device according to the operating condition; Establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation; Calculating the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature by using the Eckert reference enthalpy method; Based on the wall heat flux density, calculating the outlet temperature of the cooling medium in the target power device by using an energy conservation equation; By presetting the convection heat transfer equation, the wall temperature of the target power device is calculated according to the outlet temperature, and the temperature distribution data and the heat flux distribution data of the target power device are determined according to the wall temperature and the preset temperature, and the temperature distribution data and the heat flux distribution data are used as the thermal environment of the target power device.

2. The thermal environment calculation method of a power plant according to claim 1, characterized in that: The step of establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation includes: Establishing a first control equation of the one-dimensional control equation in the combustion chamber of the target power device according to the size parameters in the attribute information by means of a one-dimensional Euler equation; Obtaining the gas flow rate in the combustion chamber, and if the gas flow rate is greater than or equal to a threshold, calculating the wall friction coefficient in the first control equation by a preset formula, and substituting the wall friction coefficient into the first control equation to obtain a second control equation of the one-dimensional control equation; Solving the second control equation and a preset mixing model for characterizing the turbulent mixing process in the combustion chamber according to the incoming flow condition to obtain mainstream parameters in the combustion chamber of the target power device; If the gas flow rate is less than the threshold, the first control equation and the preset hybrid model are solved according to the incoming flow condition to obtain the mainstream parameters in the combustion chamber of the target power device.

3. The thermal environment calculation method of a power plant according to claim 1, characterized in that: The step of establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation, further includes: Establishing a third control equation of the one-dimensional control equation in the combustion chamber of the target power device according to the size parameters in the attribute information by an influence coefficient method; Obtaining the gas flow rate in the combustion chamber, and if the gas flow rate is greater than or equal to a threshold, calculating the wall friction coefficient in the third control equation by a preset formula, and substituting the wall friction coefficient into the third control equation to obtain a fourth control equation of the one-dimensional control equation; Solving the fourth control equation and a preset mixing model for characterizing the turbulent mixing process in the combustion chamber according to the incoming flow condition to obtain mainstream parameters in the combustion chamber of the target power device; If the gas flow rate is less than the threshold, the third control equation and the preset hybrid model are solved according to the incoming flow condition to obtain the mainstream parameters in the combustion chamber of the target power device.

4. The thermal environment calculation method of a power plant according to claim 1, characterized in that: Before the step of calculating the wall temperature of the target power device according to the outlet temperature by presetting the convection heat transfer equation, the method further includes: Acquire the type information of the cooling medium and the working condition data of the cooling medium, and select a convective heat transfer coefficient calculation formula corresponding to the type information and the working condition data from preset data according to the type information and the working condition data; Calculating the convective heat transfer coefficient of the cooling medium based on the convective heat transfer coefficient calculation formula; The initial convection heat transfer equation is updated according to the convection heat transfer coefficient to obtain the preset convection heat transfer equation.

5. The thermal environment calculation method of a power plant according to claim 1, characterized in that: The step of determining the temperature distribution data and the heat flux distribution data of the target power device according to the wall temperature and the preset temperature comprises: Calculating a first difference between the wall temperature and the preset temperature, and if an absolute value of the first difference is greater than or equal to a preset error value, adjusting the preset temperature based on the absolute value of the first difference to obtain an updated temperature; Calculating the updated wall heat flux density of the target power device according to the mainstream parameters and the updated temperature by using the Eckert reference enthalpy method; Based on the updated wall heat flux density, calculating the updated outlet temperature of the cooling medium in the target power device through an energy conservation equation; Calculating the updated wall temperature of the target power unit according to the updated outlet temperature by presetting a convection heat transfer equation; Calculating a second difference between the updated wall temperature and the updated temperature, and if an absolute value of the second difference is greater than or equal to the preset error value, adjusting the updated temperature based on the absolute value of the second difference until the absolute value of the difference between the updated wall temperature and the updated temperature is less than the preset error value; If the absolute value of the first difference is smaller than a preset error value, the temperature distribution data and the heat flux distribution data corresponding to the wall temperature are used as the temperature distribution data and the heat flux distribution data of the target power device.

6. The thermal environment calculation method of a power plant according to claim 1, characterized in that: The method further comprises: Determine the maximum limit temperature of the wall of the target power device according to the material type information in the attribute information; Acquiring operating parameters of a cooling medium in the target power device, and using the operating parameters as target inflow conditions; Calculating target mainstream parameters in the combustion chamber of the target power device based on the size parameters in the attribute information and the target incoming flow conditions; Calculating the target wall heat flux density of the target power device according to the target mainstream parameter and the preset temperature by using the Eckert reference enthalpy method; Based on the target wall heat flux density, calculating the target outlet temperature of the cooling medium in the target power device by using an energy conservation equation; By presetting a convection heat transfer equation, a target wall temperature of the target power device is calculated according to the target outlet temperature, and temperature distribution data of the target power device is determined according to the target wall temperature and the preset temperature; Selecting a maximum temperature value from the temperature distribution data, and matching the maximum temperature value with the maximum limit temperature of the wall surface to obtain a matching result; The cooling medium is evaluated according to the matching result.

7. The thermal environment calculation method of a power plant according to claim 1, characterized in that: The method further comprises: Determine the maximum limit temperature of the wall of the target power device according to the material type information in the attribute information; Taking the maximum limit temperature of the wall as the wall temperature of the target power device, and calculating a first outlet temperature of the first cooling medium according to the wall temperature and a first convection heat transfer equation corresponding to the first cooling medium; Calculating the first wall heat flux density according to the first outlet temperature by using the energy conservation equation; Calculating the first main flow parameter according to the first wall heat flux and the preset temperature by using the Eckert reference enthalpy method; Calculate a first incoming flow condition based on the size parameter in the attribute information and the first mainstream parameter; Taking the maximum limit temperature of the wall as the wall temperature of the target power device, and calculating the second outlet temperature of the second cooling medium according to the wall temperature and a second convection heat transfer equation corresponding to the second cooling medium; Calculating the second wall heat flux density according to the second outlet temperature by using an energy conservation equation; Calculating the second main flow parameter according to the second wall heat flux and the preset temperature by using the Eckert reference enthalpy method; Calculate a second incoming flow condition based on the size parameter in the attribute information and the second main flow parameter; Taking the medium flow rate of the cooling medium as a single variable, comparing the first medium flow rate in the first incoming flow condition and the second medium flow rate in the second incoming flow condition to obtain a comparison result; The first cooling medium and the second cooling medium are evaluated according to the comparison result.

8. A thermal environment calculation device for a power plant, characterized in that: The device comprises: A data acquisition module, used to determine the operating condition of the target power device through the attribute information of the target power device, and determine the incoming flow condition of the target power device according to the operating condition; A first calculation module, for establishing a one-dimensional control equation in the combustion chamber of the target power device based on the size parameters in the attribute information, and calculating the mainstream parameters in the combustion chamber of the target power device through the incoming flow condition and the one-dimensional control equation; A second calculation module is used to calculate the wall heat flux density of the target power device according to the mainstream parameters and the preset temperature by using the Eckert reference enthalpy method; A third calculation module, configured to calculate the outlet temperature of the cooling medium in the target power device based on the wall heat flux density by using an energy conservation equation; The fourth calculation module is used to calculate the wall temperature of the target power device according to the outlet temperature by using a preset convection heat transfer equation, and to determine the temperature distribution data and heat flux distribution data of the target power device according to the wall temperature and the preset temperature, and to use the temperature distribution data and heat flux distribution data as the thermal environment of the target power device.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the thermal environment calculation method of the power plant according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the thermal environment calculation method of a power plant according to any one of claims 1 to 7 are implemented.

Images