CN116680838A

CN116680838A - Heat transfer calculation method of plate-fin heat exchanger

Info

Publication number: CN116680838A
Application number: CN202310933467.0A
Authority: CN
Inventors: 谭军华; 林学星; 李林华; 张铭
Original assignee: Dongguan Pengjin Machinery Technology Co ltd
Current assignee: Dongguan Pengjin Machinery Technology Co ltd
Priority date: 2023-07-27
Filing date: 2023-07-27
Publication date: 2023-09-01
Anticipated expiration: 2043-07-27
Also published as: CN116680838B

Abstract

The invention belongs to the technical field of heat exchangers, and particularly relates to a heat transfer calculation method of a plate-fin heat exchanger, which comprises the following steps of: s1: acquiring fin structure parameters of different fluid working conditions or different plate-fin heat exchanger types; s2: respectively constructing calculation expressions of fin channel parameters, fluid physical property parameters and convection heat transfer parameters of the plate-fin heat exchanger; s3: based on different fluid working conditions or different plate-fin heat exchanger types, respectively establishing a corresponding parallel flow model, a corresponding countercurrent model and a corresponding cross flow model, and solving the outlet temperature of a fin channel of the fluid under different fluid working conditions or different plate-fin heat exchanger types by adopting programming software; s4: and respectively solving the heat exchange efficiency and the total heat transfer coefficient of different fluid working conditions or different plate-fin heat exchanger types of lower plate-fin heat exchangers based on the outlet temperature of the fluid in the fin channels.

Description

Heat transfer calculation method of plate-fin heat exchanger

Technical Field

The invention belongs to the technical field of heat exchangers, and particularly relates to a heat transfer calculation method of a plate-fin heat exchanger.

Background

In the prior art, the heat exchanger is a general device for chemical industry, petroleum, steel, automobile, video and many other industrial departments, plays an important role in production, and particularly in chemical industry production, and can be widely applied as a heater, a cooler, a condenser, an evaporator and the like. The heat transfer efficiency and the resistance characteristic of the heat exchanger can directly influence the investment cost, the operation cost and the occupied area of the whole system, so that the heat exchanger has great significance for researching the heat exchange performance of the heat exchanger.

At present, the plate-fin heat exchanger can form strong turbulence for fluid due to the special structure of the fins, so that the thermal resistance is effectively reduced, the heat transfer efficiency is improved, and meanwhile, the fins are thin and compact in structure and can be made of aluminum alloy, so that the plate-fin heat exchanger is small in size, light in weight and low in cost, and further, the plate-fin heat exchanger is widely applied to industries such as aerospace, chemical industry and refrigeration, and particularly in the field of waste heat recovery. However, in the research and development process of actual industrial products, the determination of sample parameters needs to be determined through continuous experimental analysis, and the process is complicated; meanwhile, as shown in fig. 1, in the existing NMP waste gas waste heat recovery device, a plate-fin heat exchanger with higher heat transfer efficiency mostly adopts a design mode of 'inlet and outlet section cross flow-middle section countercurrent' (i.e. an 'i+z' type or an 'l+l' type) so as to realize efficient heat exchange of two fluids without mixing. The simulation research of the plate-fin heat exchanger mostly adopts the more commonly used model software fluid dynamics (CFD) simulation software Ansys Fluent to perform model establishment, grid division and model solving, and the software has a certain degree of expertise and is difficult to directly operate by non-technical personnel; meanwhile, the systematic parameter optimization solution of the fluid dynamics (CFD) simulation software Ansys Fluent on the heat exchange mode of 'inlet and outlet section cross flow-middle section countercurrent' (namely 'I+Z' type or 'L+L' type) is difficult to realize, so that the design and optimization of the plate-fin heat exchanger capable of carrying out efficient heat exchange consume very much time cost and experiment cost.

Wherein, countercurrent section all sets up to the sawtooth structure in the plate-fin heat exchanger to strengthen heat transfer, triangle-shaped cross flow area is flat structure, in order to reduce windage pressure drop. As shown in fig. 2, in the "i+z" type plate-fin heat exchanger, the inlet and outlet directions of the hot fluid are "I", and the inlet and outlet directions of the cold fluid are "Z"; as shown in fig. 3, in the "l+l" type plate-fin heat exchanger, the inlet and outlet directions of the hot fluid and the cold fluid are both "L"; in addition, the inlets for the cold fluid and the hot fluid are respectively located at the opposite ends of the plate-fin heat exchanger.

Therefore, it is necessary to design a heat transfer algorithm of the plate-fin heat exchanger for the "i+z" type and the "l+l" type plate-fin heat exchangers capable of performing efficient heat exchange, so that the heat transfer algorithm can optimally solve the fin structure of the heat exchanger by adopting programming software, and further improve the heat transfer performance of the plate-fin heat exchanger.

Disclosure of Invention

The invention provides a heat transfer calculation method of a plate-fin heat exchanger, which is characterized in that corresponding models are respectively established through different fin types, fin structure parameters and fluid working conditions of the fins of the plate-fin heat exchanger, different fluid working conditions and different types of lower fin channel outlet temperatures of the plate-fin heat exchanger are solved, corresponding heat exchange efficiency and a heat transfer system are obtained, the heat exchange effect of the current plate-fin heat exchanger structure can be obtained, further, the experiment and labor cost are reduced, and further, a reference is provided for the design of a high-efficiency low-resistance heat exchanger.

The heat transfer calculation method of the plate-fin heat exchanger comprises the following steps:

s1: acquiring fin structure parameters of different fluid working conditions or different plate-fin heat exchanger types;

s2: respectively constructing calculation expressions of fin channel parameters, fluid physical property parameters and convection heat transfer parameters of the plate-fin heat exchanger;

s3: based on different fluid working conditions or different plate-fin heat exchanger types, respectively establishing a corresponding parallel flow model, a corresponding countercurrent model and a corresponding cross flow model, and solving the outlet temperature of a fin channel of the fluid under different fluid working conditions or different plate-fin heat exchanger types by adopting programming software;

s4: and respectively solving the heat exchange efficiency and the total heat transfer coefficient of different fluid working conditions or different plate-fin heat exchanger types of lower plate-fin heat exchangers based on the outlet temperature of the fluid in the fin channels.

Corresponding models are respectively established through different fin types, fin structure parameters and fluid working conditions of the plate-fin heat exchanger fins, different fluid working conditions and different plate-fin heat exchanger lower fin channel outlet temperatures are solved, corresponding heat exchange efficiency and a corresponding heat transfer system are obtained, the heat exchange effect of the current plate-fin heat exchanger structure can be obtained, further the experiment and labor cost are reduced, and further reference is provided for the design of the efficient low-resistance heat exchanger.

Further, in the step S1,

the fluid conditions include: parallel flow, counter flow and cross flow;

the plate fin heat exchanger types include: an 'I+Z' -shaped plate-fin heat exchanger and an 'L+L' -shaped plate-fin heat exchanger;

the fin types include: straight fins, corrugated fins, saw-tooth fins;

the fin structure parameters include: fin spacing, fin width, number of fin channel layers, fin core width, fin height, fin thickness, separator thickness, seal thickness.

Further, in the step S2, the fin channels of the plate-fin heat exchanger are formed by continuously stamping a whole aluminum plate, and the adjacent fins have a height dislocation of the fin thickness t, and the current fin and the adjacent baffle are regarded as a whole, namely, the upper side and the lower side of each fin channel can be divided into two thicknesses of the baffle, the fin and the whole baffle; the calculation of fin channel parameters of the plate-fin heat exchanger comprises the following steps:

(1) Calculation of equivalent diameter of individual Fin channels；

Equivalent diameter of straight finThe calculated expression of (2) is: />；

Equivalent diameter of corrugated finThe calculated expression of (2) is: />；

Equivalent diameter of sawtooth finThe calculated expression of (2) is: />；

in the formula ,is fin pitch, which includes cross-flow section hot side fin pitch +. >Cross-flow section cold side fin spacing +.>Fin pitch ∈>；/>Is fin height, which includes cross-flow section hot side fin height +>Cross-flow section Cold side Fin height +.>；/>Is fin thickness; />Is the ratio of the corrugated fin flow length to the radial length; />Is the pitch of the sawtooth fins;

(2) Counting the number of channels of a monolayer fin；

Number of single-layer fin channels on hot side of cross-flow sectionThe calculated expression of (2) is: />；

Number of single-layer fin channels on cold side of cross-flow sectionThe calculated expression of (2) is: />；

Number of single layer fin channels of the counterflow sectionThe calculated expression of (2) is: />；

In [ of ]]Is a rounding symbol;fin width at hot side of cross flow section; />Fin width is the cold side of the cross flow section;the thickness of the seal is the thickness of the seal;

(3) Calculating fin core height；

Height of fin coreThe calculated expression of (2) is: />；

in the formula ,the number of layers of fin channels; />Is the thickness of the separator;

(4) Calculating the windward area of the fin channel；

Hot side fin channel windward areaThe calculated expression of (2) is: />；

Windward area of cold side fin channelThe calculated expression of (2) is: />；

(5) Calculating the flow area of each fin channel；

Hot side fin channel flow areaThe calculated expression of (2) is: />；

Flow area of cold side fin channel The calculated expression of (2) is: />；

(6) Calculating the ratio of the flow area to the windward area；

Hot side fin channel flow areaIs>Ratio of->The calculated expression of (2) is:；

flow area of cold side fin channelIs>Ratio of->The calculated expression of (2) is:；

(7) Calculating the thickness of each fin channel separator；

Thickness of separatorThe calculated expression of (2) is: />；

The integral thickness of the fin and the baffleThe calculation expression of the relative position is: />；

Wherein use is made ofRepresents the thickness of the heat transfer partition in each fin channel, < > of>。

Further, in the step S2, the fluid is dry air, that is, the calculation of the corresponding fluid physical property parameter includes:

(1) Calculating the density of air in the fin channels；

Density of air in hot side fin channelsThe calculated expression of (2) is: />；

Density of air in cold side fin channelsThe calculated expression of (2) is: />；

in the formula ,to thermally measure the fin channel temperature; />The temperature of the fin channel is measured cold;

(2) Calculating the air in the wingSpecific heat capacity of sheet channel；

Specific heat capacity of air in hot side fin channelsThe calculated expression of (2) is: />；

Specific heat capacity of air in cold side fin channelsThe calculated expression of (2) is: />；

(3) Calculating the thermal conductivity of air in the Cold side Fin channel ；

Thermal conductivity of air in hot side fin channelsThe calculated expression of (2) is:；

thermal conductivity of air in cold side fin channelsThe calculated expression of (2) is:；

(4) Calculating the viscosity coefficient of air in fin channels；

Viscosity coefficient of air in hot side fin channelThe calculated expression of (2) is:；

viscosity coefficient of air in cold side fin channelThe calculated expression of (2) is:；

(5) Calculating the inlet density of air in the fin channels；

Air inlet density at hot side fin channelThe calculated expression of (2) is: />；

Air inlet density in cold side fin channelThe calculated expression of (2) is: />；

in the formula ,inlet temperature for air at the hot side fin channel; />Inlet temperature for air at the hot side fin channel;

(6) Calculating the inlet flow velocity of air in the fin channels；

Air inlet flow rate at hot side fin channelThe calculated expression of (2) is: />；

Air inlet flow rate at cold side fin channelThe calculated expression of (2) is: />；

in the formula ,the inlet air quantity of the hot side fin channel; />The air return ratio is the cold side air return ratio;

(7) Calculating the mass flow rate of air in fin channels；

Air on hot side fin channel mass flow rateThe calculated expression of (2) is: />；

Air on the cold side fin channel mass flow rateThe calculated expression of (2) is: / >；

(8) Calculating the flow velocity of air in the fin channels；

Air flow rate in hot side fin channelThe calculated expression of (2) is: />；

Air flow rate in cold side fin channelThe calculated expression of (2) is: />；

(9) Calculating Reynolds number of air in fin channel；

Reynolds number of air on hot side fin channelThe calculated expression of (2) is: />；

Reynolds number of air on cold side fin channelThe calculated expression of (2) is: />；

in the formula ,equivalent diameter of fin channels for hot side fluid flow; />Equivalent diameter of fin channels for cold side fluid flow;

(10) Calculating the prandtl number of air in fin channels；

Air-conditionerPrandtl number of fin channels on hot sideThe calculated expression of (2) is: />；

Prandtl number of air on cold side fin channelThe calculated expression of (2) is: />。

Further, in S2, the calculating the convective heat transfer parameter includes:

(1) Calculating the convective heat transfer coefficient of air in fin channels；

Convection heat transfer coefficient of air in hot side fin channelThe calculated expression of (2) is: />；

Convective heat transfer coefficient of air in cold side fin channelThe calculated expression of (2) is: />；

wherein ,；

。

(2) Calculating wind resistance of air in fin channels；

Windage of air on hot side fin channelThe calculated expression of (2) is: / >；

Windage of air on cold side fin channelThe calculated expression of (2) is: />；

in the formula ,，/>，，/>is the coefficient of foreshortening of the hot-side fin channel, +.>The sudden expansion coefficient of the hot side fin channel; />，/>，，/>Is the coefficient of foreshortening of the cold-side fin channel, +.>Is the sudden expansion coefficient of the cold side fin channel;

wherein ,；

；

further, in the step S3, the establishing of the parallel flow model and the counter flow model includes:

(1) Establishing a parallel flow model;

dividing air into infinite microelements along the flowing direction of the fin channel, and establishing a parallel flow model equation within the temperature change range of the hot fluid and the cold fluid at the corresponding temperature, namely:

；

in the formula ,for primary heat transfer of air through the partition plate in the fin passages, and (2)>Secondary heat transfer of air in the fin channels through the fins is performed; namely:

；

in the formula ,for hot side fin efficiency at the current fin infinitesimal length, its computational expression is: />, wherein ,/>；/>For the cold side fin efficiency at the current fin infinitesimal length, the calculation expression is: />, wherein ,/>；

Wherein, the boundary conditions are:

；

establishing the temperature of air in the fin passages、/>Regarding fin length->The differential equation set of (2) is:

；

(2) Establishing a countercurrent model;

；

Wherein, the boundary conditions are:

；

。

further, in the step S3, based on different fluid working conditions and different plate-fin heat exchanger types, the process of establishing the cross flow model specifically includes:

single layer fin channel number based on hot side and cold side of plate fin heat exchanger、/>Dividing a single-layer fin channel into +.>Is defined, i.e. each grid cell +.>In the region of the hot-side temperature matrix element +.>Cold side temperature matrix element->And the corresponding thermal resistance->, wherein />,/>；

In grid cellsEstablishing heat balance in an areaEquation, namely:

；

in the formula ,for grid cells->The heat transfer of the region is calculated as:

；

wherein, the boundary conditions of the 'I+Z' -shaped plate-fin heat exchanger are as follows:

；

Wherein, the boundary conditions of the L+L-shaped plate-fin heat exchanger are as follows:

。

further, in the step S3, the process of solving the outlet temperature of the fin channel of the fluid under different fluid working conditions by using programming software specifically includes:

s311: when the fluid working condition is parallel flow, based on a parallel flow model, adopting programming software to obtain a numerical solution, and respectively obtaining the outlet temperature of the hot side fin channel and the outlet temperature of the cold side fin channel under the parallel flow working condition;

s312: when the fluid working condition is countercurrent, based on a countercurrent model, adopting programming software to obtain a numerical solution, and respectively obtaining the outlet temperature of the hot side fin channel and the outlet temperature of the cold side fin channel under the countercurrent working condition;

s313: when the fluid working condition is cross flow, based on a cross flow model, adopting programming software to iteratively acquire the outlet temperature of the hot side fin channel and the outlet temperature of the cold side fin channel under the cross flow working condition;

S3131：divided based on construction of cross-flow modelsIs made->，/>Acquiring grid cell->Hot-side temperature matrix element in the region +.>Cold side temperature matrix element->And the corresponding thermal resistance->Establishing a heat balance equation to realize iteration of the whole heat exchange grid matrix;

S3132: based on the heat balance equation, solving the grid cell areaHot side temperature matrix element in (a)Grid cell area->Cold-side temperature matrix element->；

S3133: order theIterative calculation in column order until +.>Stopping iteration;

s3134: order theIterative computation is sequentially carried out according to the rows until +.>Stopping iteration;

s3135: acquiring a hot side temperature matrixLast line->Obtaining the outlet temperature of the hot-side fin channel in the case of cross-flow +.>；

Acquiring a cold side temperature matrixIs->The column average temperature was obtained to obtain the outlet temperature +.>。

Further, in the step S3, the process of solving the outlet temperature of the fin channels of the fluid under different plate-fin heat exchanger types by using programming software specifically includes:

s321: obtaining inlet temperature of hot side fin channelPresetting inlet temperature of cold side fin channel>Solving to obtain the approximate outlet temperature of the first fin channel of the cold side close to the hot side as +.>；

S322: inlet temperature based on hot side fin channelsAnd the approximate outlet temperature of the first fin channel of the cold side close to the hot side is +.>Establishing a heat exchange grid matrix and a heat balance equation, and solving the outlet temperature of the first fin channel at the cold side >；

S3221: based on the first fin channel with the cold side close to the hot side, gradually calculating and solving the inlet temperature of the first fin channel with the cold side close to the hot side according to the sequence of the first cross flow section, the counter flow section and the second cross flow section；

S3222: calculation ofWhen the value is less than +.>When the accuracy condition is satisfied, i.e., the inlet temperature of the first fin channel of the current cold side close to the hot side is obtained +.>；

S323: sequentially solving the outlet temperature of the next fin channel of the cold side close to the hot side until the outlet temperature of the nth fin channel of the cold side close to the hot side is completedIs solved;

s324: acquiring a temperature matrix of each fluid working condition stage, sequentially assigning values, and solving the outlet temperature of the hot side fin channelOutlet temperature of Cold side Fin channel +.>；

S3241: inlet temperature based on acquired hot side fin channelAnd a preset inlet temperature of the cold-side fin channel +.>；

S3242: solving a first cross-flow section hot side temperature matrix by combining a cross-flow modelFirst cross-flow section cold side temperature matrix +.>And assign the diagonal elements to the hot-side temperature matrix of the countercurrent section>Is a first row and counter-current section cold side temperature matrix +.>Respectively serving as an inlet temperature of a hot side fin channel of the countercurrent section and an outlet temperature of a cold side fin channel of the countercurrent section;

S3243: solving the outlet temperature of the hot side fin channel of the countercurrent section and the inlet temperature of the cold side fin channel of the countercurrent section by combining the countercurrent model, and respectively assigning the outlet temperature and the inlet temperature to a hot side temperature matrix of the countercurrent sectionIs a second row and countercurrent section cold side temperature matrixIn turn, the hot side temperature matrix of the countercurrent section +.>Is a second row and countercurrent section cold side temperature matrix +.>The elements of the second column of (a) are assigned to the hot-side temperature matrix of the second cross-flow section, respectively +.>Second cross-flow section cold side temperature matrix +.>Respectively serving as the inlet temperature of the hot side fin channel of the second cross-flow section and the outlet temperature of the cold side fin channel of the second cross-flow section;

s3244: solving the outlet temperature of the second cross-flow section hot side fin channel and the inlet temperature of the second cross-flow section cold side fin channel by combining the cross-flow model, and respectively assigning the outlet temperature and the inlet temperature to a second cross-flow section hot side temperature matrixThe last row, second cross-flow section cold-side temperature matrix>Is the last column of (2);

s3245: solving a second cross-flow section hot side temperature matrixThe row average temperature of the last row or the column average temperature of the last column of the heat exchanger to obtain the outlet temperature of the hot side fin channels of the I+Z-type or L+L-type plate-fin heat exchanger >；

Solving a second cross-flow section cold side temperature matrixThe row average temperature of the first row, or the column average temperature of the first column, to obtain "I+ZOr outlet temperature of hot side fin channel of L+L type plate-fin heat exchanger ≡>。

Further, in the step S4, the heat exchange efficiency, the hot side total heat transfer coefficient and the cold side total heat transfer coefficient are calculated based on the hot side outlet temperature and the cold side outlet temperature; namely:

efficiency of heat exchangeThe calculated expression of (2) is: />；

Total heat transfer coefficient of hot sideThe calculated expression of (2) is: />；

Total heat transfer coefficient of cold sideThe calculated expression of (2) is: />；

in the formula ,，/>；/>。

a computer device, comprising: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, the processor and the memory in communication over the bus when the computer device is running, the machine-readable instructions when executed by the processor performing the method of any of the above.

A computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the method of any of the preceding claims.

The beneficial effects of the invention are as follows:

according to the invention, through different fin types, fin structural parameters and fluid working conditions of the fin of the plate-fin heat exchanger, corresponding models are respectively established, different fluid working conditions and the outlet temperatures of fin channels under different plate-fin heat exchanger types are solved, corresponding heat exchange efficiency and a corresponding heat transfer system are obtained, the heat exchange effect of the current plate-fin heat exchanger structure can be obtained, further, the experiment and labor cost are reduced, and further, reference is provided for the design of the high-efficiency low-resistance heat exchanger.

Drawings

FIG. 1 is a schematic diagram of an "I+Z" or "L+L" plate-fin heat exchanger;

FIG. 2 is a schematic diagram of an "I+Z" plate-fin heat exchanger;

FIG. 3 is a schematic diagram of an "L+L" plate-fin heat exchanger;

FIG. 4 is a flow chart of the present invention;

FIG. 5 is a schematic illustration of a heat exchange grid matrix in a cross-flow model;

FIG. 6 is a flow chart for calculating the fin channel outlet temperature of a plate-fin heat exchanger under cross-flow conditions;

FIG. 7 is a flow chart for calculating the fin channel outlet temperature of a plate fin heat exchanger for different plate fin heat exchanger types;

FIG. 8 is a flow chart of a method of inlet and outlet temperature correction for cross-flow and counter-flow sections;

FIG. 9 (a) is a schematic diagram showing the distribution of the hot side temperature field of the plate-fin heat exchanger under the cross-flow condition in example 2;

FIG. 9 (b) is a schematic diagram showing the distribution of the cold side temperature field of the plate-fin heat exchanger under the cross-flow condition in example 2;

FIG. 10 (a) is a thermal side temperature field distribution diagram of an "I+Z" plate-fin heat exchanger of example 2;

FIG. 10 (b) is a graph showing the cold side temperature field profile of the "I+Z" plate-fin heat exchanger of example 2;

FIG. 11 is a system of a heat transfer calculation method for a plate-fin heat exchanger of example 3;

fig. 12 is a schematic structural diagram of a computer device.

Detailed Description

The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.

It is noted that various aspects of the embodiments are described below within the scope of the following claims. It should be apparent that the aspects described herein may be embodied in a wide variety of forms and that any specific structure and/or function described herein is merely illustrative. Based on the present disclosure, one skilled in the art will appreciate that one aspect described herein may be implemented independently of any other aspect, and that two or more of these aspects may be combined in various ways. For example, an apparatus may be implemented and/or a method practiced using any number of the aspects set forth herein. In addition, such apparatus may be implemented and/or such methods practiced using other structure and/or functionality in addition to one or more of the aspects set forth herein.

Furthermore, in the following description, specific details are provided for the purpose of providing a thorough understanding of the examples, and it will be apparent to one skilled in the art that the specific meaning of the terms described above in the present application will be practiced with specificity.

Example 1

Fig. 4 shows a heat transfer calculation method of a plate-fin heat exchanger, wherein corresponding models are respectively built through different fin types, fin structure parameters and fluid working conditions of the plate-fin heat exchanger fins, solution of different fluid working conditions and fin channel outlet temperatures under different plate-fin heat exchanger types is carried out, corresponding heat exchange efficiency and a heat transfer system are obtained, the heat exchange effect of the current plate-fin heat exchanger structure can be obtained, further experiment and labor cost are reduced, and reference is provided for design of a high-efficiency low-resistance heat exchanger. The method specifically comprises the following steps:

wherein, fluid operating mode includes: parallel flow, counter flow and cross flow;

wherein the plate fin heat exchanger types include: an 'I+Z' -shaped plate-fin heat exchanger and an 'L+L' -shaped plate-fin heat exchanger;

wherein, the fin includes: straight fins, corrugated fins, saw-tooth fins;

Wherein, fin structure parameters include: fin spacing, fin width, number of fin channel layers, fin core width, fin height, fin thickness, separator thickness, seal thickness.

the fin channels of the plate-fin heat exchanger are formed by continuously stamping a whole aluminum plate, the adjacent fins are staggered in height by the fin thickness t, and the current fins and the adjacent baffle plates are regarded as a whole, namely, the upper side and the lower side of each fin channel can be divided into two thicknesses of the baffle plates, the fins and the whole baffle plates; the calculation of fin channel parameters of the plate-fin heat exchanger comprises the following steps:

(1) Calculation of equivalent diameter of individual Fin channels；

Equivalent diameter of straight finThe calculated expression of (2) is: />；

Equivalent diameter of corrugated finThe calculated expression of (2) is: />；

Equivalent diameter of sawtooth finThe calculated expression of (2) is: />；

in the formula ,is fin pitch, which includes cross-flow section hot side fin pitch +.>Cross-flow section cold side fin spacing +.>Fin pitch ∈>；/>Is fin height, which includes cross-flow section hot side fin height +>Cross-flow section Cold side Fin height +. >；/>Is fin thickness; />Is the ratio of the corrugated fin flow length to the radial length; />Is the pitch of the sawtooth fins;

the cold side fin spacing of the countercurrent section is the same as the hot side fin spacing of the countercurrent section, and is all。

(2) Counting the number of channels of a monolayer fin；

Number of single-layer fin channels on hot side of cross-flow sectionThe calculated expression of (2) is: />；/>

Wherein, in the formula []Is a rounding symbol;fin width at hot side of cross flow section; />Fin width is the cold side of the cross flow section; />The thickness of the seal is the thickness of the seal;

(3) Calculating fin core heightCalculating the fin core height +.>；

Height of fin coreThe calculated expression of (2) is: />；

(4) Calculating the windward area of the fin channel；

Hot side fin channel windward areaThe calculated expression of (2) is: />；

(5) Calculating the flow area of each fin channel；

Hot side fin channel flow areaThe calculated expression of (2) is: />；

Flow area of cold side fin channelThe calculated expression of (2) is: / >；

(6) Calculating the ratio of the flow area to the windward area；

(7) Calculating the thickness of each fin channel separator；

Thickness of separatorThe calculated expression of (2) is: />；

The integral thickness of the fin and the baffleThe calculated expression of (2) is: />；

The fluid is dry air, physical parameters with the minimum relative error are obtained by adopting least square fitting, and the calculation of the corresponding physical parameters of the fluid at the temperature of 0-120 ℃ comprises the following steps:

(1) Calculating the density of air in the fin channels；

in the formula ,to thermally measure the fin channel temperature; />The temperature of the fin channel is measured cold; />

(2) Calculating specific heat capacity of air in fin channels；

Specific heat capacity of air in cold side fin channelsThe calculated expression of (2) is: / >；

(3) Calculating the thermal conductivity of air in the Cold side Fin channel；

(4) Calculating the viscosity coefficient of air in fin channels；

(5) Calculating the inlet density of air in the fin channels；

(6) Calculating the inlet flow velocity of air in the fin channels；

(7) Calculating the mass flow rate of air in fin channels；

Air on hot side fin channel mass flow rateThe calculated expression of (2) is: / >；

Air on the cold side fin channel mass flow rateThe calculated expression of (2) is: />；

(8) Calculating the flow velocity of air in the fin channels；

Air flow rate in hot side fin channelThe calculated expression of (2) is: />；/>

(9) Calculating Reynolds number of air in fin channel；

(10) Calculating the prandtl number of air in fin channels；

Prandtl number of air on hot side fin channelThe calculated expression of (2) is: />；

Wherein, the calculation of the convection heat transfer parameter comprises:

wherein ,；

；

(2) Calculating wind resistance of air in fin channels；

Windage of air on hot side fin channel The calculated expression of (2) is: />；

in the formula ,，/>，，/>is the coefficient of foreshortening of the hot-side fin channel, +.>The sudden expansion coefficient of the hot side fin channel; />，/>，/>，/>Is the coefficient of foreshortening of the cold-side fin channel, +.>Is the sudden expansion coefficient of the cold side fin channel;

wherein ,；

。

(1) Establishing a parallel flow model;

；

in the formula ,for hot side fin efficiency at the current fin infinitesimal length, its computational expression is: />, wherein ,/>；/>For the cold side fin efficiency at the current fin infinitesimal length, the calculation expression is: / >, wherein ,/>；

Wherein, the boundary conditions are:

；

；/>

(2) Establishing a countercurrent model;

；

Wherein, the boundary conditions are:

；

。

(3) Establishing a cross-flow model;

as shown in FIG. 5, the number of single layer fin channels based on the hot side and the cold side of the plate fin heat exchanger、/>Dividing a single-layer fin channel into +.>Is defined, i.e. each grid cell +.>Representing hot side temperature matrix elements in a regionCold side temperature matrix element- >And correspondingHeat resistance->, wherein />,；

In grid cellsEstablishing a heat balance equation in the region, namely:

；

in the formula ,for grid cells->The heat transfer of the region is calculated as: />

；

Wherein the grid cellsRegional thermal resistance->The calculated expression of (2) is:

；

。

the process of iteratively solving the outlet temperature of the fin channel of the fluid under different fluid working conditions by adopting programming software specifically comprises the following steps:

S3131: divided based on construction of cross-flow modelsIs made->，/>Acquiring grid cell->Hot-side temperature matrix element in the region +.>Cold side temperature matrix element->And the corresponding thermal resistance->Establishing a heat balance equation to realize iteration of the whole heat exchange grid matrix;

s3132: based on heat balance equation, solveHot side temperature matrix element in (a)Grid cell area->Cold-side temperature matrix element->；

s3135: acquiring a hot side temperature matrixLast line->Obtaining the row average temperature to obtain the outlet temperature of the hot side fin channel under the cross flow condition +.>；/>

Acquiring a cold side temperature matrixLast column->The column average temperature was obtained to obtain the outlet temperature +.>。

FIG. 6 is a flow chart for calculating the outlet temperature of fin channels of the plate-fin heat exchanger under the cross-flow working condition.

The process of iteratively solving the outlet temperature of the fin channels of the fluid under different plate-fin heat exchanger types by adopting programming software specifically comprises the following steps:

It should be noted that the "I+Z" type plate-fin heat exchanger can fix the inlet temperature of the hot-side fin channels onlyPresetting inlet temperature of cold side fin channel>The method comprises the steps of carrying out a first treatment on the surface of the The L+L-type fin heat exchanger can fix the inlet temperature of the hot side fin channel>Presetting inlet temperature of cold side fin channel>It is also possible to fix the inlet temperature of the cold-side fin passage +.>Presetting inlet temperature of hot side fin channel>。

S322: inlet temperature based on hot side fin channelsAnd the approximate outlet temperature of the first fin channel of the cold side close to the hot side is +.>Establishing a heat exchange grid matrix and a heat balance equation, and solving the outlet temperature of the first fin channel at the cold side>；

S3222: calculation ofWhen the value is less than +. >When the accuracy condition is satisfied, i.e., the inlet temperature of the first fin channel of the current cold side close to the hot side is obtained +.>；

Wherein after the first iterationIn the time-course of which the first and second contact surfaces,first pair->Gradually decreasing by 1 ℃ until->Then +.>Gradually increasing the temperature by 0.1 ℃ until +.>Then +.>Gradually decreasing by 0.01 ℃ until +.>Sequentially until the precision +.>Until that, the accuracy condition is satisfied.

Solving a second cross-flow section cold side temperature matrixThe row average temperature of the first row or the column average temperature of the first column to obtain the outlet temperature of the hot side fin channel of the I+Z-type or L+L-type plate-fin heat exchanger +.>。

FIG. 7 is a flow chart for calculating the fin channel outlet temperature of a plate-fin heat exchanger for different plate-fin types; wherein,for the number of iterations->The value is 0.01 as the precision condition.

In this embodiment, when the fluid flows between the fin channels with different fin pitches in the same plate-fin heat exchanger, the process of solving the temperature assignment of the fin channels between the first cross-flow section and the counter-flow section and between the counter-flow section and the second cross-flow section is further deepened, which is: as the fin spacing between the cross-flow and counter-flow sections is different, and the fin spacing between the hot and cold sides of the cross-flow sections may also be different, this results in a different number of their fin channels. To facilitate calculation of temperature change of fluid flowing through different numbers of fin channels to cross-flow hot sideThe fin channels are used as reference, the cold side of the cross-flow section is +.>The individual fin channels are subdivided into +_ according to flow area>Equal parts, but equivalent diameter is regarded as unchanged, whereby +. >And (5) a cross-flow model of the heat exchange grid matrix. After the temperature of all fin channels has been determined, the determined cold side of the cross-flow section is then subjected to +.>The individual fin channel temperature is actually +.>Correcting the inlet and outlet temperatures of each fin channel, and correcting the obtained +.>The individual fin channel temperatures are practical for the hot and cold sides of the counterflow sectionThe inlet and outlet temperatures of the individual fin passages are temperature corrected. The specific process comprises the following steps:

t1: assume that the hot side fin spacing of the cross-flow section isThe fin spacing of the cold side of the cross flow section is +.>The fin spacing of the countercurrent section isCross-flow section hot side to cold side core width ratio +.>The method comprises the steps of carrying out a first treatment on the surface of the Solving a temperature matrix of the hot side fin channels and the cold side fin channels; />

T2: the first to cold side of cross flow section, hot side and cold side of countercurrent sectionThe temperature of each fin channel is corrected to obtain the corrected +.>Individual fin channel temperature->。

Assume the first of the cold side of the cross-flow section, the hot side and the cold side of the counter-flow sectionnTemperature of individual fin channelsBy corresponding +.>Temperature->To->Temperature->The heat exchanger is obtained by numerical average mixing according to the proportion of fin spacing, namely:

cold side of cross-flow sectionnTemperature of individual fin channels The calculated expression of (2) is:

when (when)When (I)>；

Otherwise, and whenWhen (I)>；

In [ of ]]In order to round the symbol in its entirety,；/>is->Temperature of individual fin channels->The value range is +.>To->。

The hot side and the cold side of the countercurrent sectionnTemperature of individual fin channelsThe calculated expression of (2) is:

when (when)When (I)>；

Otherwise, and whenWhen (I)>；

FIG. 8 is a flow chart of a method for correcting inlet and outlet temperatures of a cross flow section and a counter flow section; wherein,the iteration times; />Is->At->，/>) Within (I) of>Is calculated.

S4: and respectively solving the heat exchange efficiency and the total heat transfer coefficient of different fluid working conditions or different plate-fin heat exchanger types of lower plate-fin heat exchangers based on the outlet temperature of the fluid in the fin channels. Namely:

efficiency of heat exchangeThe calculated expression of (2) is: />；

in the formula ,，/>；/>。

in this embodiment, MATLAB software is used as the programming software.

Example 2

In this embodiment, the difference from embodiment 1 is that: the fluid is two drying air streams, and the mass flow rate of the air entering the hot side fin channels Inlet temperature of air at hot side fin channel +.>Inlet temperature of air at hot side fin channel +.>Cold side return air ratio->。

In this embodiment, under different fluid working conditions, the fin structure parameters adopted are:

；

in this embodiment, the "i+z" type plate-fin heat exchanger, or the "l+l" type plate-fin heat exchanger, adopts the fin structure parameters:

；

after being calculated by MATLAB programming software, the obtained different fluid working conditions, the outlet temperatures of fin channels of different types of plate-fin heat exchangers, the heat exchange efficiency of the plate-fin heat exchangers, the total heat transfer coefficient and the wind resistance are specifically as follows:

；

in this embodiment, based on the above results, the OriginPro software is used to map the temperature field profiles for different fluid conditions, or different plate-fin heat exchangers.

The distribution shown in fig. 9 (a) and 9 (b) is a hot side and cold side temperature field distribution diagram of the plate-fin heat exchanger under the cross-flow working condition; FIGS. 10 (a) and 10 (b) are graphs showing hot side and cold side temperature field profiles of an "I+Z" type plate-fin heat exchanger, respectively, wherein each 25 scale on the abscissa represents a length of 1.25 m.

Example 3

As shown in fig. 11, the embodiment provides a system of a heat transfer calculation method of a plate-fin heat exchanger, which comprises a data acquisition module, a data calculation module, a model building module and a data solving module, wherein the data acquisition module is used for acquiring fin structure data of the plate-fin heat exchanger, and performing related parameter calculation to build models under different fluid working conditions so as to solve heat exchange efficiency and total heat transfer coefficient of the plate-fin heat exchanger under different fluid working conditions or different plate-fin heat exchanger types, and further acquire heat exchange conditions of the plate-fin heat exchanger, so that in practical operation, time cost and sample own cost are saved, and optimization of a fin structure is facilitated.

Specifically, the data acquisition module is used for acquiring fin structure parameters of different fluid working conditions or different plate-fin heat exchanger types;

specifically, the data calculation module is used for calculating fin channel parameters, fluid physical property parameters and convection heat transfer parameters of the plate-fin heat exchanger;

specifically, the model building module is used for respectively building a corresponding parallel flow model, a corresponding countercurrent model and a corresponding cross flow model based on different fluid working conditions or different plate-fin heat exchanger types;

specifically, the data solving module is used for solving the outlet temperature of the fin channels of the fluid under different fluid working conditions or different plate-fin heat exchanger types, and the heat exchange efficiency and the total heat transfer coefficient of the plate-fin heat exchanger under different fluid working conditions or different plate-fin heat exchanger types.

Example 4

Based on the same technical concept, as shown in fig. 12, the present embodiment further provides a computer device corresponding to the method provided in the foregoing embodiment, including a processor 2, a memory 1, and a bus, where the memory stores machine-readable instructions executable by the processor, and when the computer device is running, the processor and the memory communicate through the bus, and the machine-readable instructions are executed by the processor to perform any one of the methods described in the foregoing.

The memory 1 includes at least one type of readable storage medium including flash memory, a hard disk, a multimedia card, a card type memory (e.g., SD or DX memory, etc.), a magnetic memory, a magnetic disk, an optical disk, etc.

The memory 1 may in some embodiments be an internal storage unit of a heat transfer computing system of a plate-fin heat exchanger, such as a hard disk. The memory 1 may in other embodiments also be an external storage device of a heat transfer computing system of a plate-fin heat exchanger, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) or the like. Further, the memory 1 may also include both internal and external storage devices of the heat transfer computing system of the substrate fin heat exchanger. The memory 1 may be used not only for storing application software of a heat transfer computing system installed in a plate-fin heat exchanger and various types of data, such as codes of a heat transfer computing system program of a plate-fin heat exchanger, etc., but also for temporarily storing data that has been output or is to be output.

The processor 2 may in some embodiments be a central processing unit (Central Processing Unit, CPU), controller, microcontroller, microprocessor or other data processing chip for running program code or processing data stored in the memory 1, for example performing a heat transfer calculation program for a plate fin heat exchanger or the like.

The disclosed embodiments also provide a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs the steps of the method described in the method embodiments above. Wherein the storage medium may be a volatile or nonvolatile computer readable storage medium.

The computer program product for applying the page content refreshing method provided by the embodiment of the present invention includes a computer readable storage medium storing program codes, and the instructions included in the program codes may be used to execute the steps of the method described in the method embodiment, specifically, refer to the method embodiment and are not repeated herein.

The disclosed embodiments also provide a computer program which, when executed by a processor, implements any of the methods of the previous embodiments. The computer program product may be realized in particular by means of hardware, software or a combination thereof. In an alternative embodiment, the computer program product is embodied as a computer storage medium, and in another alternative embodiment, the computer program product is embodied as a software product, such as a software development kit (Software Development Kit, SDK), or the like.

It is to be understood that the same or similar parts in the above embodiments may be referred to each other, and that in some embodiments, the same or similar parts in other embodiments may be referred to.

It should be noted that in the description of the present invention, the terms first, second, and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Furthermore, in the description of the present invention, unless otherwise indicated, a plurality of meanings means at least two.

Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and further implementations are included within the scope of the preferred embodiment of the present invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the present invention.

It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.

Those of ordinary skill in the art will appreciate that all or a portion of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and where the program, when executed, includes one or a combination of the steps of the method embodiments.

In addition, each functional unit in the embodiments of the present invention may be integrated in one processing module, or each unit may exist alone physically, or two or more units may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented in the form of software functional modules and sold or used as a stand-alone product.

The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like.

In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.

Claims

1. The heat transfer calculation method of the plate-fin heat exchanger is characterized by comprising the following steps of:

2. A heat transfer calculation method of a plate-fin heat exchanger according to claim 1, wherein, in S1,

The fluid conditions include: parallel flow, counter flow and cross flow;

the fin types include: straight fins, corrugated fins, saw-tooth fins;

3. The method for calculating heat transfer of a plate-fin heat exchanger according to claim 1, wherein in S2, fin channels of the plate-fin heat exchanger are formed by continuously stamping a whole aluminum plate, and the adjacent fins have a height dislocation of fin thickness t, and the current fins and the adjacent baffles are regarded as a whole, namely, the upper side and the lower side of each fin channel can be divided into two thicknesses of the baffles, the fins and the whole baffles; the calculation of fin channel parameters of the plate-fin heat exchanger comprises the following steps:

(1) Calculation of equivalent diameter of individual Fin channels；

Equivalent diameter of straight finThe calculated expression of (2) is: />；

Equivalent diameter of corrugated finThe calculated expression of (2) is: />；

Equivalent diameter of sawtooth finThe calculated expression of (2) is: />；

in the formula ,Is fin pitch, which includes cross-flow section hot side fin pitch +.>Cross-flow section cold side fin spacing +.>Fin pitch ∈>；/>Is fin height, which includes cross-flow section hot side fin height +>Cross-flow section Cold side Fin height +.>；/>Is fin thickness; />Is the ratio of the corrugated fin flow length to the radial length; />Is the pitch of the sawtooth fins;

(2) Counting the number of channels of a monolayer fin；

In [ of ]]Is a rounding symbol;fin width at hot side of cross flow section; />Fin width is the cold side of the cross flow section; />The thickness of the seal is the thickness of the seal;

(3) Calculating fin core height；

Height of fin coreThe calculated expression of (2) is: />；

(4) Calculating the windward area of the fin channel；

Hot side fin channel windward areaThe calculated expression of (2) is: />；

(5) Calculating the flow area of each fin channel；

Hot side fin channel flow area The calculated expression of (2) is: />；

Flow area of cold side fin channelThe calculated expression of (2) is: />；

(6) Calculating the ratio of the flow area to the windward area；

Hot side fin channel flow areaIs>Ratio of->The calculated expression of (2) is: />；

Flow area of cold side fin channelIs>Ratio of->The calculated expression of (2) is: />；

(7) Calculating the thickness of each fin channel separator；

Thickness of separatorThe calculated expression of (2) is: />；

4. A method for calculating heat transfer of a plate-fin heat exchanger according to claim 3, wherein in S2, the fluid is dry air, i.e. the calculation of the corresponding fluid physical property parameter comprises:

(1) Calculating the density of air in the fin channels；

(2) Calculating specific heat capacity of air in fin channels；

Specific heat capacity of air in hot side fin channelsThe calculated expression of (2) is: / >；

(3) Calculating the thermal conductivity of air in the Cold side Fin channel；

(4) Calculating the viscosity coefficient of air in fin channels；

(5) Calculating the inlet density of air in the fin channels；

(6) Calculating the inlet flow velocity of air in the fin channels；

(7) Calculating the mass flow rate of air in fin channels ；

(8) Calculating the flow velocity of air in the fin channels；

Air flow rate in hot side fin channelThe calculated expression of (2) is: />；

(9) Calculating Reynolds number of air in fin channel；

(10) Calculating the prandtl number of air in fin channels；

5. The method for calculating heat transfer of a plate-fin heat exchanger according to claim 4, wherein in S2, calculating the convective heat transfer parameter comprises:

Convection heat transfer coefficient of air in hot side fin channelThe calculated expression of (2) is: / >；

wherein ,；

；

(2) Calculating wind resistance of air in fin channels；

Windage of air on hot side fin channelThe calculated expression of (2) is: />；

wherein ,；

。

6. the method for calculating heat transfer of a plate-fin heat exchanger according to claim 5, wherein in S3, the establishing of the parallel flow model and the counter flow model includes:

(1) Establishing a parallel flow model;

；

Wherein, the boundary conditions are:

；

(2) Establishing a countercurrent model;

；

Wherein, the boundary conditions are:

；

。

7. the heat transfer calculation method of a plate-fin heat exchanger according to claim 6, wherein in S3, based on different fluid working conditions and different plate-fin heat exchanger types, the cross flow model building process specifically includes:

Single layer fin channel number based on hot side and cold side of plate fin heat exchanger、/>Dividing a single-layer fin channel intoIs defined, i.e. each grid cell +.>In the region of the hot-side temperature matrix element +.>Cold side temperature matrix element->And the corresponding thermal resistance->, wherein />,/>；

In grid cellsEstablishing a heat balance equation in the region, namely:

；

。

8. the method for calculating heat transfer of a plate-fin heat exchanger according to claim 7, wherein in S3, the process of solving the outlet temperature of the fin channels of the fluid under different fluid working conditions by using programming software specifically comprises:

s3131: divided based on construction of cross-flow modelsIs made->，/>Acquiring grid cell->Hot-side temperature matrix element in the region +.>Cold side temperature matrix element->And corresponding thermal resistanceEstablishing a heat balance equation to realize iteration of the whole heat exchange grid matrix;

s3135: acquiring a hot side temperature matrixLast line->Obtaining the row average temperature to obtain the outlet temperature of the hot side fin channel under the cross flow condition +.>；

9. The method for calculating heat transfer of a plate-fin heat exchanger according to claim 8, wherein in S3, the process of solving the outlet temperature of the fin channels of the fluid under different plate-fin heat exchanger types by using programming software specifically comprises:

S322: inlet temperature based on hot side fin channelsAnd the approximate outlet temperature of the first fin channel of the cold side close to the hot side is +.>Establishing a heat exchange grid matrix and a heat balance equation, and solving the outlet temperature of a first fin channel at the cold side；

s324: acquiring a temperature matrix of each fluid working condition stage, sequentially assigning values, and solving the outlet temperature of the hot side fin channel Outlet temperature of Cold side Fin channel +.>；

S3241: inlet temperature based on acquired hot side fin channelInlet temperature of preset cold side fin channel；

s3243: solving the outlet temperature of the hot side fin channel of the countercurrent section and the inlet temperature of the cold side fin channel of the countercurrent section by combining the countercurrent model, and respectively assigning the outlet temperature and the inlet temperature to a hot side temperature matrix of the countercurrent sectionIs a second row and countercurrent section cold side temperature matrix +.>In turn, the hot side temperature matrix of the countercurrent section +.>Is a second row and countercurrent section cold side temperature matrix +.>The elements of the second column of (a) are assigned to the hot-side temperature matrix of the second cross-flow section, respectively +.>Second cross-flow section cold side temperature matrix +.>Respectively serving as the inlet temperature of the hot side fin channel of the second cross-flow section and the outlet temperature of the cold side fin channel of the second cross-flow section;

s3245: solving a second cross-flow section hot side temperature matrixThe row average temperature of the last row or the column average temperature of the last column of the heat exchanger to obtain the outlet temperature of the hot side fin channels of the I+Z-type or L+L-type plate-fin heat exchanger>；

10. The method for calculating heat transfer of a plate-fin heat exchanger according to claim 9, wherein in S4, the outlet temperature of the fin channels on the hot side is based on the fluidOutlet temperature of Cold side Fin channel +.>Calculating heat exchange efficiency, a hot side total heat transfer coefficient and a cold side total heat transfer coefficient; namely:

efficiency of heat exchangeThe calculated expression of (2) is: />；

Total heat transfer coefficient of hot sideThe calculated expression of (2) is: / >；

in the formula ,，/>；/>。

11. a computer device, comprising: a processor, a memory and a bus, the memory storing machine-readable instructions executable by the processor, the processor and the memory in communication over the bus when the computer device is running, the machine-readable instructions when executed by the processor performing the method of any one of claims 1 to 10.

12. A computer-readable storage medium, characterized in that it has stored thereon a computer program which, when executed by a processor, performs the method according to any of claims 1 to 10.