Method for calculating thermal environment of power plant and related equipment

Abstract

The application discloses a thermal environment calculation method of a power device and related equipment, and the method comprises the following steps: determining the operation condition of a target power device through attribute information of the target power device, and determining the inflow condition of the target power device according to the operation condition; calculating the main flow parameter in the combustion chamber of the target power device based on the size parameter in the attribute information and the inflow condition; calculating the wall heat flux density of the target power device according to the main flow parameter and a preset temperature by using the Eckert reference enthalpy method; calculating the outlet temperature of the cooling working medium in the target power device through an energy conservation equation based on the wall heat flux density; calculating the wall temperature of the target power device through a preset convection heat transfer equation according to the outlet temperature, and determining the temperature distribution data and the heat flow distribution data of the target power device according to the wall temperature and the preset temperature, so as to take the temperature distribution data and the heat flow distribution data as the thermal environment of the target power device.

Description

Technical Field

[0001] This invention relates to the field of thermal environment calculation technology, and in particular to a method and related equipment for calculating the thermal environment of a power unit. Background Technology

[0002] Power units often experience thermal loads during operation. The greater the power output, the greater the thermal load. With increasing demands for higher speeds, the thermal load on power units far exceeds the allowable temperature limits of the materials used in their construction. Therefore, thermal protection measures are necessary. Currently, the most effective thermal protection measure is active cooling technology, which uses a cooling medium to remove heat from the power unit's walls, thereby maintaining the temperature within a safe range and ensuring stable operation over extended periods.

[0003] However, the heat transfer process of the cooling medium within 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 wall temperature and heat flux is constrained by both cooling conditions and combustion. On the one hand, the heat flux of the power unit wall is needed as the thermal boundary condition for the cooling channel to calculate the wall temperature distribution. On the other hand, the wall temperature of the power unit is also needed to calculate the wall heat flow.

[0004] In existing technologies, research on the wall thermal environment of actively cooled power units can be divided into two categories: one is to measure the wall temperature of the power unit under given cooling conditions through experimental systems, and the other is to establish an analytical model of the combustion chamber and cooling channel of the power unit using CFD (Computational Fluid Dynamics) numerical simulation software. Because the flow field inside the combustion chamber of the power unit is very complex, with complex physical phenomena such as shock waves and boundary layer separation, and the working fluid in the cooling channel may also experience heat transfer phenomena such as phase change and decomposition, two-dimensional or three-dimensional simulations are computationally intensive, time-consuming, and inefficient, which is not conducive to the early analysis and design of the active cooling structure of the power unit. On the other hand, experimental methods are costly, have a longer cycle, and involve many uncontrollable factors. Therefore, existing methods for studying the wall thermal environment of actively cooled power units have limitations and are not conducive to the research of this issue. Summary of the Invention

[0005] Based on this, it is necessary to address the limitations of existing methods for studying the thermal environment of the walls of actively cooled power units, such as large computational load, long time consumption, and low efficiency. Therefore, a method for calculating the thermal environment of power units and related equipment are proposed.

[0006] Firstly, a method for calculating the thermal environment of a power plant is provided, the method comprising:

[0007] The operating conditions of the target power unit are determined by the attribute information of the target power unit, and the incoming flow conditions of the target power unit are determined based on the operating conditions.

[0008] Based on the size parameters in the attribute information, a one-dimensional control equation is established in the combustion chamber of the target power unit, and the mainstream parameters in the combustion chamber of the target power unit are calculated through the incoming flow conditions and the one-dimensional control equation.

[0009] The wall heat flux density of the target power unit is calculated using the Eckert reference enthalpy method based on the mainstream parameters and preset temperature.

[0010] Based on the wall heat flux density, the outlet temperature of the cooling medium inside the target power unit is calculated using the energy conservation equation.

[0011] By using a preset convective heat transfer equation, the wall temperature of the target power unit is calculated based on the outlet temperature. The temperature distribution data and heat flow distribution data of the target power unit are determined based on the wall temperature and the preset temperature. The temperature distribution data and heat flow distribution data are used as the thermal environment of the target power unit.

[0012] Optionally, the step of establishing a one-dimensional control equation for the combustion chamber within the target power unit based on the dimensional parameters in the attribute information, and calculating the main flow parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation, includes:

[0013] Using the one-dimensional Euler equation, the first control equation of the one-dimensional control equation in the combustion chamber of the target power unit is established based on the size parameters in the attribute information.

[0014] The gas flow rate in the combustion chamber is obtained. If the gas flow rate is greater than or equal to a threshold, the wall friction coefficient in the first control equation is calculated by a preset formula, and the wall friction coefficient is substituted into the first control equation to obtain the second control equation of the one-dimensional control equation.

[0015] The second governing equation and the preset mixing model used to characterize the turbulent mixing process in the combustion chamber are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit.

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

[0017] Optionally, the step of establishing a one-dimensional control equation for the combustion chamber within the target power unit based on the dimensional parameters in the attribute information, and calculating the main flow parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation, further includes:

[0018] Using the influence coefficient method, a third control equation is established based on the size parameters in the attribute information to establish the one-dimensional control equation within the combustion chamber of the target power unit;

[0019] The gas flow rate in the combustion chamber is obtained. If the gas flow rate is greater than or equal to a threshold, the wall friction coefficient in the third control equation is calculated by a preset formula, and the wall friction coefficient is substituted into the third control equation to obtain the fourth control equation of the one-dimensional control equation.

[0020] The fourth governing equation and the preset mixing model used to characterize the turbulent mixing process in the combustion chamber are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit.

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

[0022] Optionally, before the step of calculating the wall temperature of the target power unit based on the outlet temperature using a preset convective heat transfer equation, the method further includes:

[0023] Obtain the type information of the cooling medium and the operating condition data of the cooling medium, and select the convective heat transfer coefficient calculation formula corresponding to the type information and the operating condition data from the preset data according to the type information and the operating condition data.

[0024] The convective heat transfer coefficient of the cooling medium is calculated based on the aforementioned formula.

[0025] The initial convective heat transfer equation is updated based on the convective heat transfer coefficient to obtain the preset convective heat transfer equation.

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

[0027] Calculate the first difference between the wall temperature and the preset temperature. If the absolute value of the first difference is greater than or equal to a preset error value, adjust the preset temperature based on the absolute value of the first difference to obtain an updated temperature.

[0028] The heat flux density of the target power unit's updated wall surface is calculated using the Eckert reference enthalpy method, based on the mainstream parameters and the updated temperature.

[0029] Based on the updated wall heat flux density, the updated outlet temperature of the cooling medium in the target power unit is calculated using the energy conservation equation.

[0030] The renewal wall temperature of the target power unit is calculated based on the renewal outlet temperature using a preset convection heat transfer equation.

[0031] Calculate the second difference between the updated wall temperature and the updated temperature. If the absolute value of the second difference is greater than or equal to the preset error value, adjust 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 less than the preset error value, then the temperature distribution data and heat flow distribution data corresponding to the wall temperature will be used as the temperature distribution data and heat flow distribution data of the target power device.

[0033] Optionally, the method further includes:

[0034] The maximum limiting temperature of the target power unit's wall surface is determined based on the material type information in the attribute information.

[0035] Obtain the operating parameters of the cooling working fluid inside the target power unit, and use the operating parameters as the target inflow conditions;

[0036] Calculate the target mainstream parameters in the combustion chamber of the target power unit based on the size parameters in the attribute information and the target inflow conditions;

[0037] Using the Eckert reference enthalpy method, the target wall heat flux density of the target power unit is calculated based on the target mainstream parameters and the preset temperature.

[0038] Based on the target wall heat flux density, the target outlet temperature of the cooling medium in the target power unit is calculated using the energy conservation equation.

[0039] The target wall temperature of the target power unit is calculated based on the target outlet temperature using a preset convection heat transfer equation, and the temperature distribution data of the target power unit is determined based on the target wall temperature and the preset temperature.

[0040] The maximum temperature value is selected from the temperature distribution data, and the maximum temperature value is matched with the highest limit temperature of the wall surface to obtain the matching result;

[0041] The evaluation of the cooling medium is completed based on the matching results.

[0042] Optionally, the method further includes:

[0043] The maximum limiting temperature of the target power unit's wall surface is determined based on the material type information in the attribute information.

[0044] The highest limit temperature of the wall surface is taken as the wall surface temperature of the target power device, and the first outlet temperature of the first cooling medium is calculated based on the wall surface temperature and the first convective heat transfer equation corresponding to the first cooling medium.

[0045] The heat flux density of the first wall is calculated based on the first outlet temperature using the energy conservation equation.

[0046] Using the Eckert reference enthalpy method, the first mainstream parameters are calculated based on the heat flux density of the first wall and the preset temperature;

[0047] The first incoming flow condition is calculated based on the size parameters in the attribute information and the first main flow parameter.

[0048] The highest limit temperature of the wall surface is taken as the wall temperature of the target power unit, and the second outlet temperature of the second cooling medium is calculated based on the wall temperature and the second convective heat transfer equation corresponding to the second cooling medium.

[0049] The heat flux density of the second wall is calculated based on the second outlet temperature using the energy conservation equation.

[0050] The second mainstream parameters are calculated using the Eckert reference enthalpy method based on the second wall heat flux density and the preset temperature.

[0051] The second incoming flow condition is calculated based on the size parameters and the second main flow parameters in the attribute information;

[0052] By using the cooling medium flow rate as a single variable, the first medium flow rate under the first incoming flow condition and the second medium flow rate under the second incoming flow condition are compared to obtain the comparison result;

[0053] The first cooling medium and the second cooling medium are evaluated based on the comparison results.

[0054] Secondly, this application provides a thermal environment calculation device for a power unit, the device comprising:

[0055] The data acquisition module is used to determine the operating conditions of the target power unit through the attribute information of the target power unit, and to determine the incoming flow conditions of the target power unit based on the operating conditions.

[0056] The first calculation module is used to establish a one-dimensional control equation for the combustion chamber inside the target power unit based on the size parameters in the attribute information, and to calculate the main flow parameters inside the combustion chamber inside the target power unit through the incoming flow conditions and the one-dimensional control equation.

[0057] The second calculation module is used to calculate the wall heat flux density of the target power unit based on the main parameters and preset temperature using the Eckert reference enthalpy method.

[0058] The third calculation module is used to calculate the outlet temperature of the cooling medium in the target power unit based on the wall heat flux density and the energy conservation equation.

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

[0060] Thirdly, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the thermal environment calculation method for the power unit as described above.

[0061] Fourthly, a computer-readable storage medium is provided, the computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the above-described method for calculating the thermal environment of a power device.

[0062] This application determines the operating conditions of the target power unit by means of its attribute information, and determines the incoming flow conditions of the target power unit based on the operating conditions; establishes a one-dimensional control equation for the combustion chamber of the target power unit based on the dimensional parameters in the attribute information, and calculates the mainstream parameters of the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation; calculates the wall heat flux density of the target power unit using the Eckert reference enthalpy method based on the mainstream parameters and a preset temperature; calculates the outlet temperature of the cooling medium in the target power unit using the energy conservation equation based on the wall heat flux density; calculates the wall temperature of the target power unit based on the outlet temperature using a preset convective heat transfer equation, and determines the temperature distribution data and heat flux distribution data of the target power unit based on 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 unit. By establishing one-dimensional control equations to calculate the main parameters in the combustion chamber of the target power unit, the amount of calculation is reduced and the calculation efficiency is improved. The wall temperature and heat flux density of the power unit under active cooling conditions can be calculated quickly, thereby improving the efficiency of the study on the thermal environment of the wall of the active cooling power unit. Attached Figure Description

[0063] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0064] in:

[0065] Figure 1 A flowchart of a method for calculating the thermal environment of a power unit in one embodiment;

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

[0067] Figure 3 This is a flowchart illustrating the evaluation of the cooling performance of different cooling media in one embodiment.

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

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

[0070] Figure 6 This is a structural block diagram of a computer device in another embodiment. Detailed Implementation

[0071] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0072] The present invention will now be described in detail through specific embodiments.

[0073] Please see Figure 1 As shown, Figure 1 A flowchart illustrating the thermal environment calculation method for a power unit provided in an embodiment of the present invention includes the following steps:

[0074] S101. Determine the operating conditions of the target power unit through the attribute information of the target power unit, and determine the incoming flow conditions of the target power unit based on the operating conditions;

[0075] For example, the attribute information includes, but is not limited to, the type of the target power device and the size parameters of the target power device. The type of the target power device can be understood as the category of the target power device, and the size parameters of the target power device can be understood as the geometric parameters of the target power device. 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 related to 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] For example, the process of determining the operating condition of the target power device through the attribute information of the target power device can be achieved by using the operating parameters in the attribute information of the target power device, such as speed and operating angle.

[0077] S102. Based on the size parameters in the attribute information, establish a one-dimensional control equation for the combustion chamber inside the target power unit, and calculate the main flow parameters inside the combustion chamber of the target power unit through the incoming flow conditions and the one-dimensional control equation.

[0078] In one possible implementation, the step of establishing a one-dimensional control equation for the combustion chamber within the target power unit based on the dimensional parameters in the attribute information, and calculating the mainstream parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation, includes:

[0079] Using the one-dimensional Euler equation, the first control equation of the one-dimensional control equation in the combustion chamber of the target power unit is established based on the size parameters in the attribute information.

[0080] The gas flow rate in the combustion chamber is obtained. If the gas flow rate is greater than or equal to a threshold, the wall friction coefficient in the first control equation is calculated by a preset formula, and the wall friction coefficient is substituted into the first control equation to obtain the second control equation of the one-dimensional control equation.

[0081] The second governing equation and the preset mixing model used to characterize the turbulent mixing process in the combustion chamber are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit.

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

[0083] For example, when establishing the one-dimensional control equations in the combustion chamber of the target power unit, the one-dimensional control equations in the combustion chamber of the target power unit can be established through the one-dimensional Euler equation. The Euler equation has high accuracy and comprehensiveness, and the efficiency of the solution process of the main parameters is guaranteed by adopting the one-dimensional Euler equation.

[0084] For example, by taking into account the influence of the gas flow rate on the wall friction coefficient, the accuracy of the calculation of the main parameters is ensured.

[0085] For example, by characterizing the turbulent mixing process in the combustion chamber using a preset mixing model, this solution can be applied to different types of power units.

[0086] S103. Using the Eckert reference enthalpy method, calculate the wall heat flux density of the target power unit based on the mainstream parameters and preset temperature;

[0087] For example, the Eckert reference enthalpy method uses the Reynolds analogy principle to perform qualitative temperature calculations. Therefore, a preset temperature needs to be set, and the enthalpy value at the wall temperature is calculated through the preset temperature. 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 * The variable below, reference temperature T * It can be done through the formula:

[0091]

[0092]

[0093] We obtain, where h w h is the enthalpy at the wall temperature. aw Let λ be the adiabatic wall enthalpy, r be the gas restitution coefficient, μ be the viscosity, and λ be the thermal conductivity. The wall heat flux density q at this temperature can be obtained using the above formula. w Pr* is the Prandtl number, P e U is the Peckley number, R is the gas constant, and u is the gas constant. e r* represents the flow rate, and r* represents the reference temperature T. * The airflow recovery coefficient, μ is the viscosity of the airflow, λ* is the thermal conductivity, and C p * represents the specific heat at constant pressure, T represents the actual temperature, and h represents the enthalpy at the actual temperature. * Reference temperature T * The enthalpy value below, Re is the Reynolds number, Re x The Reynolds number is the location of the target power unit.

[0094] S104. Based on the wall heat flux density, calculate the outlet temperature of the cooling medium in the target power unit using the energy conservation equation;

[0095] For example, the high-temperature gas generated in the combustion chamber of the target power device transfers heat to the wall of the target power device. The cooling medium 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 medium in the wall channel exchanges energy with the wall of the target power device, and then the outlet temperature of the cooling medium in the target power device is calculated by the energy conservation equation and the heat flux density of the wall.

[0096] For example, 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 exchanged between the cooling medium in the wall channel and the target power device wall. It can be understood that the heat exchanged between the cooling medium in the wall channel and the target power device wall is absorbed by the cooling medium, and the outlet temperature of the cooling medium can be calculated.

[0097] S105. Calculate the wall temperature of the target power unit based on the outlet temperature using a preset convection heat transfer equation, and determine the temperature distribution data and heat flow distribution data of the target power unit based on the wall temperature and the preset temperature, using the temperature distribution data and heat flow distribution data as the thermal environment of the target power unit.

[0098] In one possible implementation, prior to the step of calculating the wall temperature of the target power unit based on the outlet temperature using a preset convection heat transfer equation, the method further includes:

[0099] Obtain the type information of the cooling medium and the operating condition data of the cooling medium, and select the convective heat transfer coefficient calculation formula corresponding to the type information and the operating condition data from the preset data according to the type information and the operating condition data.

[0100] The convective heat transfer coefficient of the cooling medium is calculated based on the aforementioned formula.

[0101] The initial convective heat transfer equation is updated based on the convective heat transfer coefficient to obtain the preset convective heat transfer equation.

[0102] For example, the calculation formula for the convective heat transfer coefficient h differs for different cooling media, and the calculation formula for the convective heat transfer coefficient of the same media under different operating conditions will also differ. Therefore, it is necessary to select the convective heat transfer coefficient calculation formula corresponding to the type of cooling media and the operating condition data of the cooling media to calculate the convective heat transfer coefficient of the cooling media.

[0103] For example, in scenarios where the type of coolant and / or its operating data cannot be determined, the Gnielinski formula is selected:

[0104]

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

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

[0107] In one possible implementation, the step of establishing a one-dimensional control equation for the combustion chamber within the target power unit based on the dimensional parameters in the attribute information, and calculating the mainstream parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation, further includes:

[0108] Using the influence coefficient method, a third control equation is established based on the size parameters in the attribute information to establish the one-dimensional control equation within the combustion chamber of the target power unit;

[0109] The gas flow rate in the combustion chamber is obtained. If the gas flow rate is greater than or equal to a threshold, the wall friction coefficient in the third control equation is calculated by a preset formula, and the wall friction coefficient is substituted into the third control equation to obtain the fourth control equation of the one-dimensional control equation.

[0110] The fourth governing equation and the preset mixing model used to characterize the turbulent mixing process in the combustion chamber are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit.

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

[0112] For example, using the influence coefficient method, considering factors such as combustion, friction, and cross-sectional changes, the basic control equation (the third control equation of the one-dimensional control equation) is constructed based on the dimensional parameters in the attribute information as follows:

[0113]

[0114]

[0115] Where dA represents the cross-sectional change, and C f T is the coefficient of wall friction. 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 T is the equivalent diameter. t For the total incoming temperature, y is the fuel flow rate, γ is the specific heat ratio, y is the ratio of the flow 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 correlate the above formulas.

[0116] When the gas flow velocity in the combustion chamber is high, such as in supersonic combustion, the wall friction coefficient C f The value is relatively large and cannot be ignored; it can be addressed using a preset formula.

[0117]

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

[0119] Considering that one-dimensional calculations cannot simulate the actual turbulent mixing process, a mixing model can be used to describe a pre-defined mixing model. For example, for a ramp-type nozzle, the mathematical expression of the pre-defined mixing model is:

[0120]

[0121]

[0122] Where, η m For mixing efficiency, exp() is the natural exponential function, x represents the length of different structures in the ramp nozzle, A1 is a constant, usually ranging from 1 to 5, and L m This indicates the mixing length, which can be selected based on the specific power unit configuration.

[0123] 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 internal combustion gas along the flow direction (mainstream parameters).

[0124] In one possible implementation, the step of determining the temperature distribution data and heat flux distribution data of the target power unit based on the wall temperature and the preset temperature includes:

[0125] Calculate the first difference between the wall temperature and the preset temperature. If the absolute value of the first difference is greater than or equal to a preset error value, adjust the preset temperature based on the absolute value of the first difference to obtain an updated temperature.

[0126] The heat flux density of the target power unit's updated wall surface is calculated using the Eckert reference enthalpy method, based on the mainstream parameters and the updated temperature.

[0127] Based on the updated wall heat flux density, the updated outlet temperature of the cooling medium in the target power unit is calculated using the energy conservation equation.

[0128] The renewal wall temperature of the target power unit is calculated based on the renewal outlet temperature using a preset convection heat transfer equation.

[0129] Calculate the second difference between the updated wall temperature and the updated temperature. If the absolute value of the second difference is greater than or equal to the preset error value, adjust 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.

[0130] If the absolute value of the first difference is less than the preset error value, then the temperature distribution data and heat flow distribution data corresponding to the wall temperature will be used as the temperature distribution data and heat flow distribution data of the target power device.

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

[0132]

[0133] Where m c For the cooling working fluid flow rate, T c,1 C is the initial temperature of the cooling working fluid. p,c This represents the specific heat capacity of the cooling medium, which is updated numerically based on temperature changes. Q w The total heat is equal to the product of the heat flux density and the heat transfer area. For example, when the channel is rectangular, the total heat can be expressed as:

[0134]

[0135] Where W is the length of the channel rectangle and L is the width of the channel rectangle.

[0136] After obtaining the temperature of the cooling medium, the wall temperature T can be calculated using the convective heat transfer equation between the wall and the medium. w2 :

[0137]

[0138] In the formula T c To cool the average temperature of the working fluid, A r A f To account for the rib effect, the area of ​​the rib root and the rib surface, eta The rib efficiency is calculated using a formula corresponding to the channel shape, h, depending on the channel shape. c The convective heat transfer coefficient of the cooling working fluid.

[0139] For example, the preset temperature is denoted as T. W1 The wall temperature T of the cooling channel is obtained. W2 Then, it should be compared with the assumed wall temperature T of the power unit. W1 In contrast, if the two are different, it is necessary to re-assume the wall temperature T of the power unit. W1 The process continues by repeating the steps described above: using the Eckert reference enthalpy method to calculate the wall heat flux density of the target power unit based on the mainstream parameters and preset temperature; and using the energy conservation equation to calculate the outlet temperature of the cooling medium inside the target power unit based on the wall heat flux density. This process is iterated until the calculated wall temperature meets the error. The temperature distribution and heat flux distribution obtained at this point represent the thermal environment of the power unit under active cooling conditions.

[0140] In one possible implementation, the method further includes:

[0141] The maximum limiting temperature of the target power unit's wall surface is determined based on the material type information in the attribute information.

[0142] Obtain the operating parameters of the cooling working fluid inside the target power unit, and use the operating parameters as the target inflow conditions;

[0143] Calculate the target mainstream parameters in the combustion chamber of the target power unit based on the size parameters in the attribute information and the target inflow conditions;

[0144] Using the Eckert reference enthalpy method, the target wall heat flux density of the target power unit is calculated based on the target mainstream parameters and the preset temperature.

[0145] Based on the target wall heat flux density, the target outlet temperature of the cooling medium in the target power unit is calculated using the energy conservation equation.

[0146] The target wall temperature of the target power unit is calculated based on the target outlet temperature using a preset convection heat transfer equation, and the temperature distribution data of the target power unit is determined based on the target wall temperature and the preset temperature.

[0147] The maximum temperature value is selected from the temperature distribution data, and the maximum temperature value is matched with the highest limit temperature of the wall surface to obtain the matching result;

[0148] The evaluation of the cooling medium is completed based on the matching results.

[0149] For example, when evaluating the cooling effect of any cooling medium, considering that the primary goal of active cooling is to protect the propulsion unit walls from high-temperature damage, the highest wall temperature is selected as the main evaluation index for cooling performance. For evaluating the cooling performance of different cooling media, considering the aircraft's sensitivity to overall mass, the coolant mass flow rate is selected as the main evaluation index. Furthermore, depending on the purpose, other evaluation indices may be used, such as the maximum coolant velocity, the average temperature of the propulsion unit walls, or the maximum heat flux density of the walls.

[0150] Evaluation of the cooling performance of any working fluid, such as Figure 2The process begins by determining the flow rate of the working fluid and selecting initial parameters for the working fluid at the inlet of the cooling channel, primarily including inlet pressure and temperature (equivalent to using these operating parameters as target inflow conditions, which include the initial parameters). The maximum wall temperature Tmax is then calculated (equivalent to calculating the target mainstream parameters within the combustion chamber of the target power unit 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 unit is calculated based on 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 within the target power unit is calculated using the energy conservation equation; the target wall temperature of the target power unit is calculated based on the target outlet temperature using the preset convective heat transfer equation, and the temperature distribution data of the target power unit is determined based on the target wall temperature and the preset temperature; the maximum temperature value is selected from the temperature distribution data). Then, based on the selected wall material type of the power unit, the allowable maximum wall temperature limit Tlimit is determined (equivalent to determining the maximum wall temperature limit of the target power unit based on the material type information in the attribute information). The calculated maximum wall temperature Tmax is compared with the allowable wall temperature limit Tlimit. If the calculated maximum wall temperature Tmax is greater than the allowable wall temperature limit Tlimit, the cooling medium cannot meet the cooling requirements of the power unit, and the medium flow rate or inlet parameters need to be changed and redesigned. If the calculated maximum wall temperature Tmax is less than the allowable wall temperature limit Tlimit, the surface cooling medium can meet the cooling requirements (equivalent to matching the maximum temperature value with the maximum wall temperature limit to obtain a matching result; the evaluation of the cooling medium is completed based on the matching result).

[0151] In one possible implementation, the method further includes:

[0152] The maximum limiting temperature of the target power unit's wall surface is determined based on the material type information in the attribute information.

[0153] The highest limit temperature of the wall surface is taken as the wall surface temperature of the target power device, and the first outlet temperature of the first cooling medium is calculated based on the wall surface temperature and the first convective heat transfer equation corresponding to the first cooling medium.

[0154] The heat flux density of the first wall is calculated based on the first outlet temperature using the energy conservation equation.

[0155] Using the Eckert reference enthalpy method, the first mainstream parameters are calculated based on the heat flux density of the first wall and the preset temperature;

[0156] The first incoming flow condition is calculated based on the size parameters in the attribute information and the first main flow parameter.

[0157] The highest limit temperature of the wall surface is taken as the wall temperature of the target power unit, and the second outlet temperature of the second cooling medium is calculated based on the wall temperature and the second convective heat transfer equation corresponding to the second cooling medium.

[0158] The heat flux density of the second wall is calculated based on the second outlet temperature using the energy conservation equation.

[0159] The second mainstream parameters are calculated using the Eckert reference enthalpy method based on the second wall heat flux density and the preset temperature.

[0160] The second incoming flow condition is calculated based on the size parameters and the second main flow parameters in the attribute information;

[0161] By using the cooling medium flow rate as a single variable, the first medium flow rate under the first incoming flow condition and the second medium flow rate under the second incoming flow condition are compared to obtain the comparison result;

[0162] The first cooling medium and the second cooling medium are evaluated based on the comparison results.

[0163] For example, such as Figure 3 As shown, working medium A (the first cooling medium) is selected, the inlet parameters of the cooling channel are determined, and the required mass flow rate m of working medium A when the highest wall temperature of the power unit is Tlimit is calculated. A Selecting working medium B (second cooling medium), and using the same working medium inlet parameters for the cooling channel, calculate the required mass flow rate m of working medium B when the highest wall temperature of the power unit is Tlimit. B Compare the calculated mass flow rate m of working fluid A. A Mass flow rate m of working fluid B B Size, if the mass flow rate m of working fluid A A The mass flow rate m of working fluid B is less than B This indicates that the cooling effect of working fluid A is better than that of working fluid B. If the mass flow rate m of working fluid A is... A The mass flow rate m greater than working fluid B B This indicates that the cooling effect of working fluid B is better than that of working fluid A.

[0164] Secondly, such as Figure 4 As shown, this application provides a thermal environment calculation device for a power unit, the device comprising:

[0165] The data acquisition module 201 is used to determine the operating conditions of the target power device through the attribute information of the target power device, and to determine the incoming flow conditions of the target power device based on the operating conditions.

[0166] The first calculation module 202 is used to establish a one-dimensional control equation in the combustion chamber of the target power unit based on the size parameters in the attribute information, and to calculate the main flow parameters in the combustion chamber of the target power unit through the incoming flow conditions and the one-dimensional control equation.

[0167] The second calculation module 203 is used to calculate the wall heat flux density of the target power unit based on the mainstream parameters and the preset temperature using the Eckert reference enthalpy method.

[0168] The third calculation module 204 is used to calculate the outlet temperature of the cooling medium in the target power unit based on the wall heat flux density and the energy conservation equation.

[0169] The fourth calculation module 205 is used to calculate the wall temperature of the target power unit based on the outlet temperature using a preset convection heat transfer equation, and to determine the temperature distribution data and heat flow distribution data of the target power unit based on the wall temperature and the preset temperature, and to use the temperature distribution data and heat flow distribution data as the thermal environment of the target power unit.

[0170] The operating conditions of the target power unit are determined by analyzing its attribute information, and the incoming flow conditions are determined accordingly. A one-dimensional control equation is established within the combustion chamber of the target power unit based on the dimensional parameters in the attribute information. The mainstream parameters within the combustion chamber are calculated using the incoming flow conditions and the one-dimensional control equation. The wall heat flux density of the target power unit is calculated using the Eckert reference enthalpy method, based on the mainstream parameters and a preset temperature. The outlet temperature of the cooling medium within the target power unit is calculated using the energy conservation equation based on the wall heat flux density. The wall temperature of the target power unit is calculated based on the outlet temperature using a preset convective heat transfer equation. The temperature distribution data and heat flux distribution data of the target power unit are determined based on the wall temperature and the preset temperature, and these data are used as the thermal environment of the target power unit. By establishing one-dimensional control equations to calculate the main parameters in the combustion chamber of the target power unit, the amount of calculation is reduced and the calculation efficiency is improved. The wall temperature and heat flux density of the power unit under active cooling conditions can be calculated quickly, thereby improving the efficiency of the study on the thermal environment of the wall of the active cooling power unit.

[0171] In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows: Figure 5As shown, the computer device includes a processor, memory, network interface, and database connected via a system bus. The processor provides computational and control capabilities. The memory includes non-volatile and / or volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and database. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The network interface is used to communicate with external clients via a network connection. When the computer program is executed by the processor, it implements the functions or steps of a server-side method for calculating the thermal environment of a power unit.

[0172] In one embodiment, a computer device is provided, which may be a client, and its internal structure diagram may be as follows: Figure 5 As shown, the computer device includes a processor, memory, network interface, display screen, and input devices connected via a system bus. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The network interface is used to communicate with an external server via a network connection. When the computer program is executed by the processor, it implements the client-side functions or steps of a thermal environment calculation method for a power unit.

[0173] In one embodiment, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it performs the following steps: determining the operating conditions of the target power unit based on attribute information, and determining the incoming flow conditions of the target power unit based on the operating conditions; establishing a one-dimensional control equation within the combustion chamber of the target power unit based on the dimensional parameters in the attribute information, and calculating the mainstream parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation; calculating the wall heat flux density of the target power unit based on the mainstream parameters and a preset temperature using the Eckert reference enthalpy method; calculating the outlet temperature of the cooling medium within the target power unit using an energy conservation equation based on the wall heat flux density; calculating the wall temperature of the target power unit based on the outlet temperature using a preset convective heat transfer equation, and determining the temperature distribution data and heat flux distribution data of the target power unit based on the wall temperature and the preset temperature, using the temperature distribution data and heat flux distribution data as the thermal environment of the target power unit.

[0174] In one embodiment, a computer-readable storage medium is provided, the computer-readable storage medium storing a computer program, which, when executed by a processor, performs the following steps: determining the operating conditions of the target power unit based on attribute information of the target power unit, and determining the incoming flow conditions of the target power unit based on the operating conditions; establishing a one-dimensional control equation within the combustion chamber of the target power unit based on the dimensional parameters in the attribute information, and calculating the mainstream parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation; calculating the wall heat flux density of the target power unit based on the mainstream parameters and a preset temperature using the Eckert reference enthalpy method; calculating the outlet temperature of the cooling medium within the target power unit using an energy conservation equation based on the wall heat flux density; calculating the wall temperature of the target power unit based on the outlet temperature using a preset convective heat transfer equation, and determining the temperature distribution data and heat flux distribution data of the target power unit based on 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 unit.

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

[0176] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, 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. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0177] Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the above-described division of functional units and modules is used as an example. In practical applications, the above functions can be assigned to 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.

[0178] The above-described embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included within the protection scope of the present invention.

Claims

1. A method for calculating the thermal environment of a power unit, characterized in that, The method includes: The operating conditions of the target power unit are determined by the attribute information of the target power unit, and the incoming flow conditions of the target power unit are determined based on the operating conditions. Based on the size parameters in the attribute information, a one-dimensional control equation is established in the combustion chamber of the target power unit, and the mainstream parameters in the combustion chamber of the target power unit are calculated through the incoming flow conditions and the one-dimensional control equation. The wall heat flux density of the target power unit is calculated using the Eckert reference enthalpy method based on the mainstream parameters and preset temperature. Based on the wall heat flux density, the outlet temperature of the cooling medium inside the target power unit is calculated using the energy conservation equation. By using a preset convective heat transfer equation, the wall temperature of the target power unit is calculated based on the outlet temperature. The temperature distribution data and heat flow distribution data of the target power unit are determined based on the wall temperature and the preset temperature. The temperature distribution data and heat flow distribution data are used as the thermal environment of the target power unit. The step of determining the temperature distribution data and heat flux distribution data of the target power unit based on the wall temperature and the preset temperature includes: Calculate the first difference between the wall temperature and the preset temperature. If the absolute value of the first difference is greater than or equal to a preset error value, adjust the preset temperature based on the absolute value of the first difference to obtain an updated temperature. The heat flux density of the target power unit's updated wall surface is calculated using the Eckert reference enthalpy method, based on the mainstream parameters and the updated temperature. Based on the updated wall heat flux density, the updated outlet temperature of the cooling medium in the target power unit is calculated using the energy conservation equation. The renewal wall temperature of the target power unit is calculated based on the renewal outlet temperature using a preset convection heat transfer equation. Calculate the second difference between the updated wall temperature and the updated temperature. If the absolute value of the second difference is greater than or equal to the preset error value, adjust 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 less than the preset error value, then the temperature distribution data and heat flow distribution data corresponding to the wall temperature will be used as the temperature distribution data and heat flow distribution data of the target power device.

2. The method for calculating the thermal environment of a power unit according to claim 1, characterized in that, The step of establishing a one-dimensional control equation for the combustion chamber within the target power unit based on the dimensional parameters in the attribute information, and calculating the main flow parameters within the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation, includes: Using the one-dimensional Euler equation, the first control equation of the one-dimensional control equation in the combustion chamber of the target power unit is established based on the size parameters in the attribute information. The gas flow rate in the combustion chamber is obtained. If the gas flow rate is greater than or equal to a threshold, the wall friction coefficient in the first control equation is calculated by a preset formula, and the wall friction coefficient is substituted into the first control equation to obtain the second control equation of the one-dimensional control equation. The second governing equation and the preset mixing model used to characterize the turbulent mixing process in the combustion chamber are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit. If the gas flow rate is less than the threshold, the first control equation and the preset mixing model are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit.

3. The method for calculating the thermal environment of a power unit according to claim 1, characterized in that, The step of establishing a one-dimensional control equation for the combustion chamber of the target power unit based on the size parameters in the attribute information, and calculating the mainstream parameters of the combustion chamber of the target power unit using the incoming flow conditions and the one-dimensional control equation, further includes: Using the influence coefficient method, a third control equation is established based on the size parameters in the attribute information to establish the one-dimensional control equation within the combustion chamber of the target power unit; The gas flow rate in the combustion chamber is obtained. If the gas flow rate is greater than or equal to a threshold, the wall friction coefficient in the third control equation is calculated by a preset formula, and the wall friction coefficient is substituted into the third control equation to obtain the fourth control equation of the one-dimensional control equation. The fourth governing equation and the preset mixing model used to characterize the turbulent mixing process in the combustion chamber are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit. If the gas flow velocity is less than the threshold, the third control equation and the preset mixing model are solved according to the incoming flow conditions to obtain the mainstream parameters in the combustion chamber of the target power unit.

4. The method for calculating the thermal environment of a power unit according to claim 1, characterized in that, Before the step of calculating the wall temperature of the target power unit based on the outlet temperature using a preset convection heat transfer equation, the method further includes: Obtain the type information of the cooling medium and the operating condition data of the cooling medium, and select the convective heat transfer coefficient calculation formula corresponding to the type information and the operating condition data from the preset data according to the type information and the operating condition data. The convective heat transfer coefficient of the cooling medium is calculated based on the aforementioned formula. The initial convective heat transfer equation is updated based on the convective heat transfer coefficient to obtain the preset convective heat transfer equation.

5. The method for calculating the thermal environment of a power unit according to claim 1, characterized in that, The method further includes: The maximum limiting temperature of the target power unit's wall surface is determined based on the material type information in the attribute information. Obtain the operating parameters of the cooling working fluid inside the target power unit, and use the operating parameters as the target inflow conditions; Calculate the target mainstream parameters in the combustion chamber of the target power unit 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 unit is calculated based on the target mainstream parameters and the preset temperature. Based on the target wall heat flux density, the target outlet temperature of the cooling medium in the target power unit is calculated using the energy conservation equation. The target wall temperature of the target power unit is calculated based on the target outlet temperature using a preset convection heat transfer equation, and the temperature distribution data of the target power unit is determined based on the target wall temperature and the preset temperature. The maximum temperature value is selected from the temperature distribution data, and the maximum temperature value is matched with the highest limit temperature of the wall surface to obtain the matching result; The evaluation of the cooling medium is completed based on the matching results.

6. The method for calculating the thermal environment of a power unit according to claim 1, characterized in that, The method further includes: The maximum limiting temperature of the target power unit's wall surface is determined based on the material type information in the attribute information. The highest limit temperature of the wall surface is taken as the wall surface temperature of the target power device, and the first outlet temperature of the first cooling medium is calculated based on the wall surface temperature and the first convective heat transfer equation corresponding to the first cooling medium. The heat flux density of the first wall is calculated based on the first outlet temperature using the energy conservation equation. Using the Eckert reference enthalpy method, the first mainstream parameters are calculated based on the heat flux density of the first wall and the preset temperature; The first incoming flow condition is calculated based on the size parameters in the attribute information and the first main flow parameter. The highest limit temperature of the wall surface is taken as the wall temperature of the target power unit, and the second outlet temperature of the second cooling medium is calculated based on the wall temperature and the second convective heat transfer equation corresponding to the second cooling medium. The heat flux density of the second wall is calculated based on the second outlet temperature using the energy conservation equation. The second mainstream parameters are calculated using the Eckert reference enthalpy method based on the second wall heat flux density and the preset temperature. The second incoming flow condition is calculated based on the size parameters and the second main flow parameters in the attribute information; By using the cooling medium flow rate as a single variable, the first medium flow rate under the first incoming flow condition and the second medium flow rate under the second incoming flow condition are compared to obtain the comparison result; The first cooling medium and the second cooling medium are evaluated based on the comparison results.

7. A thermal environment calculation device for a power unit, characterized in that, The device includes: The data acquisition module is used to determine the operating conditions of the target power unit through the attribute information of the target power unit, and to determine the incoming flow conditions of the target power unit based on the operating conditions. The first calculation module is used to establish a one-dimensional control equation for the combustion chamber inside the target power unit based on the size parameters in the attribute information, and to calculate the main flow parameters inside the combustion chamber inside the target power unit through the incoming flow conditions and the one-dimensional control equation. The second calculation module is used to calculate the wall heat flux density of the target power unit based on the main parameters and preset temperature using the Eckert reference enthalpy method. The third calculation module is used to calculate the outlet temperature of the cooling medium in the target power unit based on the wall heat flux density and the energy conservation equation. The fourth calculation module is used to calculate the wall temperature of the target power unit based on the outlet temperature using a preset convective heat transfer equation, and to determine the temperature distribution data and heat flux distribution data of the target power unit based on the wall temperature and the preset temperature, using the temperature distribution data and heat flux distribution data as the thermal environment of the target power unit. The step of determining the temperature distribution data and heat flux distribution data of the target power unit based on the wall temperature and the preset temperature includes: calculating a first difference between the wall temperature and the preset temperature; if the 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; and using the Eckert reference enthalpy method to calculate the target power unit's temperature distribution data based on the main flow parameters and the updated temperature. The heat flux density of the updated wall surface of the power device is calculated; based on the updated wall surface heat flux density, the updated outlet temperature of the cooling medium in the target power device is calculated using the energy conservation equation; the updated wall surface temperature of the target power device is calculated based on the updated outlet temperature using a preset convective heat transfer equation; a second difference between the updated wall surface temperature and the updated temperature is calculated; if the absolute value of the second difference is greater than or equal to the preset error value, the updated temperature is adjusted based on the absolute value of the second difference until the absolute value of the difference between the updated wall surface temperature and the updated temperature is less than the preset error value; if the absolute value of the first difference is less than the preset error value, the temperature distribution data and heat flux distribution data corresponding to the wall surface temperature are used as the temperature distribution data and heat flux distribution data of the target power device.

8. 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, it implements the steps of the thermal environment calculation method for the power unit as described in any one of claims 1 to 6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the thermal environment calculation method for the power unit as described in any one of claims 1 to 6.

Images