CN108090307A - Plate-fin heat exchanger channel layout design method under a kind of multi-state based on integral mean temperature differential method - Google Patents
Plate-fin heat exchanger channel layout design method under a kind of multi-state based on integral mean temperature differential method Download PDFInfo
- Publication number
- CN108090307A CN108090307A CN201810041244.2A CN201810041244A CN108090307A CN 108090307 A CN108090307 A CN 108090307A CN 201810041244 A CN201810041244 A CN 201810041244A CN 108090307 A CN108090307 A CN 108090307A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- channel layout
- mfrac
- design
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Abstract
The invention discloses plate-fin heat exchanger channel layout design methods under a kind of multi-state based on integral mean temperature differential method.It is primarily based on designing a model between integral mean temperature differential method structure channel layout, fin structure parameter, cold fluid and hot fluid entrance condition parameter and heat exchanger heat exchanger efficiency;Then multi-state channel layout coefficient is defined, constructs multi-state channel layout coordinating intervals;Finally using channel layout coordinating intervals as design domain, it is up to design object with multi-state weighting heat exchange amount, realizes the plate-fin heat exchanger channel layout optimal design for considering multi-state design requirement.The characteristics of this method is the channel layout design that can be suitably used for plate-fin heat exchanger under multi-state, strengthens the heat transfer process that ultra-large type surpasses plate-fin heat exchanger, excavates its energy-saving potential.
Description
Technical field
The present invention relates to a kind of plate-fin heat exchanger channel layout design method, more particularly, to one kind towards under multi-state
Ultra-large type plate-fin heat exchanger channel layout design method.
Background technology
Cold fluid and hot fluid species is more in plate-fin heat exchanger, number of active lanes is more, and the species of channel layout mode is with port number
Purpose increases to be increased into exponential type, and the heat exchange amount of different channel layouts is different, and therefore, channel layout design is that plate-fin changes
One key issue of hot device design.
Existing channel layout mode mostly according to the number of active lanes of cold fluid and hot fluid, the entrance condition parameter of cold fluid and hot fluid,
It balances 2 principles according to cold fluid and hot fluid isolation arrangement, local heat load and arranges optimal design into row of channels.These methods can obtain
Optimal channel arrangement under given duty parameter, and heat exchange property is improved, promote enhanced heat exchange, the property of plate-fin heat exchanger
It can improve.But these methods all can be only applied to the ultra-large type plate-fin heat exchanger design under either simplex condition, it is impossible to suitable for multiplexing
The channel layout design of heat exchanger under condition.
Often operating condition is complicated for existing ultra-large type plate-fin heat exchanger, the heat exchanger used in being equipped such as Large Air Separation Devices
The demand of adaptation 80%-110% oxygen load variations is needed, the heat exchange amount of heat exchanger is needed according to yield in the production of ethylene enterprise
Variation adaptively adjusted, therefore, how to account for complex working condition design requirement channel layout design be ultra-large type
The key that plate-fin heat exchanger enhanced heat exchange, energy-saving potential excavate.
The content of the invention
For the above-mentioned problems in the prior art, in order to strengthen ultra-large type plate-fin heat exchanger under complex working condition
Heat exchange property excavates its energy-saving potential, and the object of the invention is with providing a kind of ultra-large type plate-fin heat exchanger towards multi-state
Channel layout design method.
Plate-fin heat exchanger channel layout design method under a kind of multi-state based on integral mean temperature differential method,
It is characterized in that comprising the following steps:
1) based on integral mean temperature differential method structure channel layout, fin structure parameter, cold fluid and hot fluid entrance condition parameter with
Designing a model between heat exchanger heat exchanger efficiency;
2) channel layout coefficient and channel layout design space are defined;
3) multiple either simplex condition channel layout design results are carried out based on Hybrid Particle Swarm, determines multi-state channel layout
Coordinating intervals;
4) using channel layout coordinating intervals as design domain, it is up to design object with multi-state weighting heat exchange amount, is examined
Consider the plate-fin heat exchanger channel layout optimal design of multi-state design requirement.
Plate-fin heat exchanger channel layout design method under a kind of multi-state based on integral mean temperature differential method,
It is characterized in that the model construction of step 1) uses following steps:
1.1) plate-fin heat exchanger is seen as 3 stream heat exchangers by interchannel heat exchange series connection to form, for by N number of
The plate-fin heat exchanger of passage fluid composition, wherein 3 strands be made of i-th strand, (i+1) stock and (i+2) stock passage fluid
It flows sub- heat exchanger and is defined as i-th (i ∈ [1, N-2]) a sub- heat exchanger;
1.2) by taking i-th is sub- heat exchanger as an example, thermodynamic analysis is carried out, heat exchange amount Q is established and is flowed with heat exchanger fluid
Direction j, heat exchange area A, the specific heat at constant pressure C of each fluidp, quantitative relationship between overall heat-transfer coefficient U, the i-th fluid streams and the
(i+1) the temperature difference T between fluid streamsi+1,i, the temperature difference T of (i+1) between fluid streams and (i+2) fluid streamsi+2,i+1Respectively such as
Shown in formula (1) and (2):
Wherein,
Wherein, r1And r2It is equation r2+(-a1-a5)r+(a1a5-a2a4Two real roots of)=0, C1~C4Be according to x=0 or
4 constants of the differential thermal calculation of stock i~(i+2) are flowed at person x=L;For the mass flow of fluid, CpIt is specific heat at constant pressure, T is
Temperature, L are heat exchanger length, and j is heat exchanger fluid flow direction, A is heat exchange area, U is overall heat-transfer coefficient, and e is nature pair
Number;Subscript i represents the i-th fluid streams;QiIt is heat output of i-th fluid streams to (i+1) fluid streams, QiDuring > 0, the i-th fluid streams to
(i+1) fluid streams are conducted heat, QiDuring < 0, (i+1) fluid streams are conducted heat to the i-th fluid streams;
1.3) from x=0 to x=L, formula (1) is integrated using integral mean temperature differential method, then divided by heat exchanger length L,
Determine the integration heat transfer temperature difference of stream stock i and (i+1)As shown in formula (4);Stream stock (i+1) is determined using same method
The integration heat transfer temperature difference of (i+2)As shown in formula (5):
1.4) according to the stream stock i and the integration heat transfer temperature difference of (i+1) being calculatedFlow the biography between stock i and (i+1)
Heat QiShown in computational methods such as formula (6);
Wherein εiIt is channel factor, shown in the outlet temperature such as formula (7) for the passage i being calculated:
1.5) total heat exchange amount Q of heat exchanger is the sum of total amount of heat of all hot fluid outflows, as shown in formula (8):
Wherein subscript hot represents all hot fluids in all heat exchangers.
Ultra-large type plate-fin heat exchanger channel layout is set under a kind of complex working condition based on integral mean temperature differential method
Meter method, it is characterised in that the channel layout coefficient of step 2) is defined as follows:
The demand of heat exchange amount under the given operating mode of channel layout coefficient reflection, for any design conditions, channel layout system
Number η is demand heat exchange amount minimum value Q under the operating modereqWith attainable heat exchange amount maximum Q under the operating modemaxRatio, according to
The structural behaviour principle of continuity, channel layout coefficient η reflect the size in the operating mode lower channel layout designs space, when demand most
Small heat exchange amount QreqWith heat exchange amount maximum QmaxUnanimously, i.e. Qreq=Qmax, channel layout coefficient η=1, and channel layout design is empty
Between be only a design point;As the minimum heat exchange amount Q of demandreq>=0, channel layout coefficient η=0, channel layout design space is
The global design space of channel layout.
Ultra-large type plate-fin heat exchanger channel layout is set under a kind of complex working condition based on integral mean temperature differential method
Meter method, it is characterised in that step 3) based on multiple either simplex condition channel layout design results, determine that multi-state channel layout is assisted
Section is adjusted to comprise the steps of:
3.1) the maximum heat exchange amount Q under either simplex condition is calculated using Hybrid Particle Swarmmax,iWith optimal channel layout type,
In Hybrid Particle Swarm, design object is the total heat exchange amount maximized under either simplex condition;Heredity is determined according to total number of channels amount
Algebraical sum is often for population quantity;Cross method is intersected at random for particle preferably individual and population preferably individual;Variation method
For the random variation of individual;
If 3.2) the optimal channel layout result under all operating modes is the same, the channel layout coefficient of each operating mode is
1, channel layout design space is 1 design point, which corresponds to the optimal channel layout of either simplex condition;
3.3) if the optimal channel layout result of each operating mode is different, by coordinating the channel layout coefficient of each operating mode, structure
The channel layout design space of multi-state is built, is as follows:
3.3.1 the channel layout coefficient for) assuming each operating mode is the same, and by reducing channel layout coefficient, structure is initial
Channel layout design space, if there is m key Design operating mode, it is assumed that the coefficient of each operating mode is η*, then i-th crucial operating mode
Heat exchange amount design requirement is Qi≥η*Qmax,i, η is made first*=1, according to the heat exchange amount demand under each operating mode, channel layout design area
Between be empty set Φ;Then η is reduced with the step delta η of a very little*, i.e. η*=η*- Δ η extends the channel layout design of each operating mode
Section;If the coordinate design section of multi-state remains as empty set Φ, continue to reduce η with step delta η*, until multi-state is coordinated
Channel layout design space is not empty set, remembers that channel layout coefficient at this time is η ';
3.3.2 the channel layout coefficient of each operating mode) is improved, further shrinks multi-state lower channel layout designs space;From
Crucial operating mode 1 arrives crucial operating mode m, increases channel layout coefficient η using small step delta ηi;Under crucial operating mode i, work as channel layout
Coefficient ηiIncrease to η 'i, multi-state coordinate path layout designs space is not empty set Φ;But ηiIncrease to η 'i+ Δ η, multiplexing
Condition coordinate path layout designs space is empty set Φ, then η 'iFor the optimal channel location coefficient under crucial operating mode i;
3.4) the optimal channel location coefficient η ' under each crucial operating mode is determinediAfterwards, the heat exchange amount range of needs under each operating mode
For Qhot,i≥η′iQmax,i, channel layout design space is all channel layout modes for meeting heat exchange amount demand under each operating mode
Set.
Ultra-large type plate-fin heat exchanger channel layout is set under a kind of complex working condition based on integral mean temperature differential method
Meter method, the channel layout optimal design of the multi-state design requirement of the step 4) use following steps:
4.1) build multi-state channel layout majorized function and consider m crucial operating condition design requirement, it is excellent to establish multi-state
Change function G (x) as shown in formula (9), the coefficient w of each operating modeiAs shown in formula (10):
In formula, x is passage arrangement mode, yiIt is the design parameter of i-th of crucial operating mode.The weight coefficient w roots of each operating mode
It is determined according to the Estimated Time Of Operation t of each operating mode;Consider two special circumstances, when heat exchanger always works at i-th of crucial operating mode,
Then wi=1 and wj=0 (1≤j≤m, j ≠ i);When heat exchanger needs frequently to carry out Off-design operation, then heat exchanger is m pass
Run time under key operating mode is essentially identical, then
4.2) the overall heat exchange amount that multi-state channel layout mathematical optimization models majorized function is maximization formula (9) is built,
The heat exchange amount constraint for each operating mode that design constraint determines for channel layout coefficient, as shown in formula (11):
Qhot,i(x,yi)≥η′iQi,maxi∈[1,m] (11);
4.3) multi-state lower channel layout optimization design design parameter is carried out using hybrid particle swarm method to arrange for passage
Mode, maximum genetic algebra and often determines for population invariable number according to number of active lanes, particle individual and the optimal particle individual of history into
Row intersects, and realizes the optimization design of channel layout.
By using above-mentioned technology, compared with prior art, the invention has the advantages that:
(1) present invention establishes the Heat Exchangers thermodynamical model based on integral mean temperature differential method, realizes basis
Total heat exchange amount of heat exchanger structure parameter, duty parameter and channel layout mode calculates;
(2) present invention proposes multi-state channel layout cooperation index, it is proposed that " is first reduced, rear to increase " under complex working condition
Multi-state channel layout coordinating intervals construction method, it is determined that consider complex working condition design requirement channel layout design domain,
Channel layout optimization has been carried out using particle cluster algorithm, has realized channel layout optimization design;
(3) invention applies the complex working condition channel layout optimization design of 24 stream heat exchangers, and total heat exchange amount is used
Result verification analysis has been carried out with two methods of accumulative heat load distribution, it is shown that the channel layout that the method for the present invention obtains has
Effect property.
Description of the drawings
Fig. 1 is the flow chart of the present invention;
Fig. 2 is multi-state lower channel layout optimization design flow chart;
The optimal channel that Fig. 3 is to determine is laid out the accumulative heat load distribution figure under 80% operating mode;
The optimal channel that Fig. 4 is to determine is laid out the accumulative heat load distribution figure under 90% operating mode;
The optimal channel that Fig. 5 is to determine is laid out the accumulative heat load distribution figure under standard condition;
The optimal channel that Fig. 6 is to determine is laid out the accumulative heat load distribution figure under 110% operating mode.
Specific embodiment
With reference to plate-fin heat exchanger channel layout under complex working condition, the invention will be further described.
As shown in Figs. 1-2, plate-fin heat exchanger channel layout under the multi-state of the invention based on integral mean temperature differential method
Design method comprises the following steps:
1) based on integral mean temperature differential method structure channel layout, fin structure parameter, cold fluid and hot fluid entrance condition parameter with
Designing a model between heat exchanger heat exchanger efficiency, model construction process use following steps:
1.1) plate-fin heat exchanger is seen as 3 stream heat exchangers by interchannel heat exchange series connection to form, for by N number of
The plate-fin heat exchanger of passage fluid composition, wherein 3 strands be made of i-th strand, (i+1) stock and (i+2) stock passage fluid
It flows sub- heat exchanger and is defined as i-th (i ∈ [1, N-2]) a sub- heat exchanger;
1.2) by taking i-th is sub- heat exchanger as an example, thermodynamic analysis is carried out, heat exchange amount Q is established and is flowed with heat exchanger fluid
Direction j, heat exchange area A, the specific heat at constant pressure C of each fluidp, quantitative relationship between overall heat-transfer coefficient U, the i-th fluid streams and the
(i+1) the temperature difference T between fluid streamsi+1,i, the temperature difference T of (i+1) between fluid streams and (i+2) fluid streamsi+2,i+1Respectively such as
Shown in formula (1) and (2):
Wherein,
Wherein, r1And r2It is equation r2+(-a1-a5)r+(a1a5-a2a4Two real roots of)=0, C1~C4Be according to x=0 or
4 constants of the differential thermal calculation of stock i~(i+2) are flowed at person x=L;For the mass flow of fluid, CpIt is specific heat at constant pressure, T is
Temperature, L are heat exchanger length, and j is heat exchanger fluid flow direction, A is heat exchange area, U is overall heat-transfer coefficient, and e is nature pair
Number;Subscript i represents the i-th fluid streams;QiIt is heat output of i-th fluid streams to (i+1) fluid streams, QiDuring > 0, the i-th fluid streams to
(i+1) fluid streams are conducted heat, QiDuring < 0, (i+1) fluid streams are conducted heat to the i-th fluid streams;
1.3) from x=0 to x=L, formula (1) is integrated using integral mean temperature differential method, then divided by heat exchanger length L,
Determine the integration heat transfer temperature difference of stream stock i and (i+1)As shown in formula (4);Stream stock (i+1) is determined using same method
The integration heat transfer temperature difference of (i+2)As shown in formula (5):
1.4) according to the stream stock i and the integration heat transfer temperature difference of (i+1) being calculatedFlow the biography between stock i and (i+1)
Heat QiShown in computational methods such as formula (6);
Wherein εiIt is channel factor, shown in the outlet temperature such as formula (7) for the passage i being calculated:
1.5) total heat exchange amount Q of heat exchanger is the sum of total amount of heat of all hot fluid outflows, as shown in formula (8):
Wherein subscript hot represents all hot fluids in all heat exchangers;
2) define channel layout coefficient and channel layout design space, the channel layout coefficient are defined as follows:Passage
The demand of heat exchange amount under the given operating mode of location coefficient reflection, for any design conditions, channel layout coefficient η is under the operating mode
Demand heat exchange amount minimum value QreqWith attainable heat exchange amount maximum Q under the operating modemaxRatio, according to structural behaviour continuity
Principle, channel layout coefficient η reflect the size in the operating mode lower channel layout designs space, as the minimum heat exchange amount Q of demandreqWith
Heat exchange amount maximum QmaxUnanimously, i.e. Qreq=Qmax, channel layout coefficient η=1, and channel layout design space is only one and sets
Enumeration;As the minimum heat exchange amount Q of demandreq>=0, channel layout coefficient η=0, channel layout design space is the whole of channel layout
Body design space;
3) multiple either simplex condition channel layout design results are carried out based on Hybrid Particle Swarm, determines multi-state channel layout
Coordinating intervals specifically comprise the steps of:
3.1) the maximum heat exchange amount Q under either simplex condition is calculated using Hybrid Particle Swarmmax,iWith optimal channel layout type,
In Hybrid Particle Swarm, design object is the total heat exchange amount maximized under either simplex condition;Heredity is determined according to total number of channels amount
Algebraical sum is often for population quantity;Cross method is intersected at random for particle preferably individual and population preferably individual;Variation method
For the random variation of individual;
If 3.2) the optimal channel layout result under all operating modes is the same, the channel layout coefficient of each operating mode is
1, channel layout design space is 1 design point, which corresponds to the optimal channel layout of either simplex condition;
3.3) if the optimal channel layout result of each operating mode is different, by coordinating the channel layout coefficient of each operating mode, structure
The channel layout design space of multi-state is built, is as follows:
3.3.1 the channel layout coefficient for) assuming each operating mode is the same, and by reducing channel layout coefficient, structure is initial
Channel layout design space, if there is m key Design operating mode, it is assumed that the coefficient of each operating mode is η*, then i-th crucial operating mode
Heat exchange amount design requirement is Qi≥η*Qmax,i, η is made first*=1, according to the heat exchange amount demand under each operating mode, channel layout design area
Between be empty set Φ;Then η is reduced with the step delta η of a very little*, i.e. η*=η*- Δ η extends the channel layout design of each operating mode
Section;If the coordinate design section of multi-state remains as empty set Φ, continue to reduce η with step delta η*, until multi-state is coordinated
Channel layout design space is not empty set, remembers that channel layout coefficient at this time is η ';
3.3.2 the channel layout coefficient of each operating mode) is improved, further shrinks multi-state lower channel layout designs space;From
Crucial operating mode 1 arrives crucial operating mode m, increases channel layout coefficient η using small step delta ηi;Under crucial operating mode i, work as channel layout
Coefficient ηiIncrease to η 'i, multi-state coordinate path layout designs space is not empty set Φ;But ηiIncrease to η 'i+ Δ η, multiplexing
Condition coordinate path layout designs space is empty set Φ, then η 'iFor the optimal channel location coefficient under crucial operating mode i;
3.4) the optimal channel location coefficient η ' under each crucial operating mode is determinediAfterwards, the heat exchange amount range of needs under each operating mode
For Qhot,i≥η′iQmax,i, channel layout design space is all channel layout modes for meeting heat exchange amount demand under each operating mode
Set;
4) using channel layout coordinating intervals as design domain, it is up to design object with multi-state weighting heat exchange amount, is examined
The plate-fin heat exchanger channel layout optimal design of multi-state design requirement is considered, specifically using following steps:
4.1) build multi-state channel layout majorized function and consider m crucial operating condition design requirement, it is excellent to establish multi-state
Change function G (x) as shown in formula (9), the coefficient w of each operating modeiAs shown in formula (10):
In formula, x is passage arrangement mode, yiIt is the design parameter of i-th of crucial operating mode.The weight coefficient w roots of each operating mode
It is determined according to the Estimated Time Of Operation t of each operating mode;Consider two special circumstances, when heat exchanger always works at i-th of crucial operating mode,
Then wi=1 and wj=0 (1≤j≤m, j ≠ i);When heat exchanger needs frequently to carry out Off-design operation, then heat exchanger is m pass
Run time under key operating mode is essentially identical, then
4.2) the overall heat exchange amount that multi-state channel layout mathematical optimization models majorized function is maximization formula (9) is built,
The heat exchange amount constraint for each operating mode that design constraint determines for channel layout coefficient, as shown in formula (11):
Qhot,i(x,yi)≥η′iQi,maxi∈[1,m] (11);
4.3) multi-state lower channel layout optimization design design parameter is carried out using hybrid particle swarm method to arrange for passage
Mode, maximum genetic algebra and often determines for population invariable number according to number of active lanes, particle individual and the optimal particle individual of history into
Row intersects, and realizes the optimization design of channel layout.
Embodiment:
The complex working condition lower channel layout for carrying out 24 stream heat exchangers comprising each 3 strands of A, B, C, I, J, K, L and M is set
Meter, the fluid that 24 stream heat exchangers use and fin design parameter are as shown in table 1:
The fluid that 1 24 stream heat exchanger of table uses and fin design parameter
Used fin material be aluminium, heat transfer coefficient 0.19158kW/mK.Hot fluid is flowed along the increased directions of x
Dynamic, cold fluid is flowed along the direction of x reductions, and therefore, the flow direction of cold fluid and hot fluid is represented respectively with 1 and -1.Heat exchanger
Total length and width are respectively 1m and 0.5m.
Account for 24 stream heat exchangers under 80%-110% variable working condition requirements (including flowing stock A, B, C, F, G, H and N)
Channel layout design.Either simplex condition channel layout result under standard condition, 80% operating mode, 90% operating mode and 110% operating mode is such as
Shown in table 2, layout 1, layout 2, layout 3 and layout 4 are respectively designated as.Coordinate criterion using multi-state, according to specific embodiment party
3.3 in formula) the step of, show that the channel layout cooperation index of each operating mode is respectively 0.945,0.875,0.915 and 0.890, structure
Multi-state passage coordinate design space is built, optimal channel arrangement mode is the layout 5 in table 2 under complex working condition.
The channel layout design result of 2 24 stream heat exchanger of table
So far, there are no the Optimality Criteria of complex working condition lower channel arrangement design, therefore, 2 kinds of indirect methods are used
Detect the correctness and validity of the complex working condition channel layout design result of this chapter.
First method is channel layout that the traditional either simplex condition method of comparison and the method for the present invention obtain in various operative employees
The heat exchange amount of condition, as shown in table 3.
Heat exchange amount of the different passage arrangement modes of table 3 under each operation operating mode
From table 3 it can be seen that although heat exchange amount of the layout 5 that optimizes of complex working condition under each operating mode is not most
Excellent, it is optimal that the layout 1- of tradition list operating condition method, which is laid out 4 heat exchange amount under corresponding operation operating mode, but is laid out 5
Good heat exchange amount can be obtained under each operating mode, but it is bad to be laid out 1- 4 heat exchange amounts under other operating modes of layout.Therefore, originally
The layout 5 that inventive method obtains meets heat exchanger channel layout designs requirement under multi-state.
Second method is the accumulative heat load distribution situation under each operating mode of analysis.Accumulative thermic load method calculating passage i,
From passage 1 to the incidence relation the accumulative thermic load of passage i.It is designed for a good channel layout, adds up thermic load one
Directly oscillated about 0;It is designed for a bad channel layout, accumulative thermic load does not stop to accumulate in one direction, Ran Hou
Negative direction does not stop to accumulate.If cold fluid and hot fluid arranged crosswise, the accumulative thermic load between cold fluid and hot fluid is ceaselessly compensated, and adds up heat
Load is always near 0, then this kind of channel layout mode can obtain a good heat exchanger efficiency.According to the criterion, analyze
Under standard condition, 80% operating mode, 90% operating mode and 110% operating mode, 5 accumulative thermic load such as Fig. 3, Fig. 4, Fig. 5, Fig. 6 institute is laid out
Show, from figure and table the results show that the obtained layout 5 of the method for the present invention uniformly obtains well under each operation operating mode
Accumulative heat load distribution also shows the validity for the channel layout that the method for the present invention obtains.
Claims (5)
1. a kind of plate-fin heat exchanger channel layout design method under multi-state based on integral mean temperature differential method, it is characterised in that
Comprise the following steps:
1) based on integral mean temperature differential method structure channel layout, fin structure parameter, cold fluid and hot fluid entrance condition parameter and heat exchange
Designing a model between device heat exchanger efficiency;
2) channel layout coefficient and channel layout design space are defined;
3) multiple either simplex condition channel layout design results are carried out based on Hybrid Particle Swarm, determines that multi-state channel layout is coordinated
Section;
4) using channel layout coordinating intervals as design domain, it is up to design object with multi-state weighting heat exchange amount, accounts for more
The plate-fin heat exchanger channel layout optimal design of operating condition design requirement.
2. plate-fin heat exchanger channel layout under a kind of multi-state based on integral mean temperature differential method according to claim 1
Design method, it is characterised in that the model construction of step 1) uses following steps:
1.1) plate-fin heat exchanger is seen as 3 stream heat exchangers by interchannel heat exchange series connection to form, for by N number of passage
The plate-fin heat exchanger of fluid composition, wherein 3 plumes being made of i-th strand, (i+1) stock and (i+2) stock passage fluid
Heat exchanger is defined as i-th (i ∈ [1, N-2]) a sub- heat exchanger;
1.2) by taking i-th is sub- heat exchanger as an example, thermodynamic analysis is carried out, establishes heat exchange amount Q and heat exchanger fluid flow direction
J, the specific heat at constant pressure C of heat exchange area A, each fluidp, quantitative relationship between overall heat-transfer coefficient U, the i-th fluid streams and (i+1)
Temperature difference T between fluid streamsi+1,i, the temperature difference T of (i+1) between fluid streams and (i+2) fluid streamsi+2,i+1Respectively such as formula (1)
(2) shown in:
<mrow>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mi>x</mi>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mi>x</mi>
</mrow>
</msup>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mi>x</mi>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>4</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mi>x</mi>
</mrow>
</msup>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>U</mi>
<mi>A</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mi>i</mi>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>U</mi>
<mi>A</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>U</mi>
<mi>A</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mi>i</mi>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>U</mi>
<mi>A</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>U</mi>
<mi>A</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>U</mi>
<mi>A</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<msub>
<mi>j</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</msub>
<mi>L</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, r1And r2It is equation r2+(-a1-a5)r+(a1a5-a2a4Two real roots of)=0, C1~C4It is according to x=0 or x
4 constants of the differential thermal calculation of stock i~(i+2) are flowed at=L;For the mass flow of fluid, CpIt is specific heat at constant pressure, T is temperature
Degree, L are heat exchanger length, and j is heat exchanger fluid flow direction, A is heat exchange area, U is overall heat-transfer coefficient, and e is natural logrithm;
Subscript i represents the i-th fluid streams;QiIt is heat output of i-th fluid streams to (i+1) fluid streams, QiDuring > 0, the i-th fluid streams are to
(i+1) fluid streams are conducted heat, QiDuring < 0, (i+1) fluid streams are conducted heat to the i-th fluid streams;
Then divided by heat exchanger length L 1.3) from x=0 to x=L, formula (1) is integrated using integral mean temperature differential method, is determined
Flow the integration heat transfer temperature difference of stock i and (i+1)As shown in formula (4);Stream stock (i+1) and (i+ are determined using same method
2) integration heat transfer temperature differenceAs shown in formula (5):
<mrow>
<msubsup>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>I</mi>
<mi>M</mi>
<mi>T</mi>
<mi>D</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>L</mi>
</msubsup>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mi>L</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>L</mi>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>I</mi>
<mi>M</mi>
<mi>T</mi>
<mi>D</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>L</mi>
</msubsup>
<msub>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>L</mi>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>r</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
<mi>L</mi>
</mrow>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>r</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>L</mi>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
</mrow>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
1.4) according to the stream stock i and the integration heat transfer temperature difference of (i+1) being calculatedFlow the heat output between stock i and (i+1)
QiShown in computational methods such as formula (6);
<mrow>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mi>i</mi>
</msub>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>L</mi>
</msubsup>
<msub>
<mrow>
<mo>(</mo>
<mi>U</mi>
<mi>A</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>x</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>U</mi>
<mi>A</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<msubsup>
<mi>&Delta;T</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>I</mi>
<mi>M</mi>
<mi>T</mi>
<mi>D</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein εiIt is channel factor, shown in the outlet temperature such as formula (7) for the passage i being calculated:
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mrow>
<mo>(</mo>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
1.5) total heat exchange amount Q of heat exchanger is the sum of total amount of heat of all hot fluid outflows, as shown in formula (8):
<mrow>
<mi>Q</mi>
<mo>=</mo>
<mi>&Sigma;</mi>
<msub>
<mrow>
<mo>(</mo>
<mover>
<mi>m</mi>
<mo>&CenterDot;</mo>
</mover>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mi>h</mi>
<mi>o</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mrow>
<mo>{</mo>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
<mrow>
<mi>h</mi>
<mi>o</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein subscript hot represents all hot fluids in all heat exchangers.
3. ultra-large type plate-fin heat exchanger under a kind of complex working condition based on integral mean temperature differential method according to claim 1
Channel layout design method, it is characterised in that the channel layout coefficient of step 2) is defined as follows:
The demand of heat exchange amount under the given operating mode of channel layout coefficient reflection, for any design conditions, channel layout coefficient η is
Demand heat exchange amount minimum value Q under the operating modereqWith attainable heat exchange amount maximum Q under the operating modemaxRatio, according to structural
The energy principle of continuity, channel layout coefficient η reflect the size in the operating mode lower channel layout designs space, when the minimum heat exchange of demand
Measure QreqWith heat exchange amount maximum QmaxUnanimously, i.e. Qreq=Qmax, channel layout coefficient η=1, and channel layout design space is only
One design point;As the minimum heat exchange amount Q of demandreq>=0, channel layout coefficient η=0, channel layout design space is passage cloth
The global design space of office.
4. ultra-large type plate-fin heat exchanger under a kind of complex working condition based on integral mean temperature differential method according to claim 1
Channel layout design method, it is characterised in that step 3) based on multiple either simplex condition channel layout design results, determine multi-state
Channel layout coordinating intervals comprise the steps of:
3.1) the maximum heat exchange amount Q under either simplex condition is calculated using Hybrid Particle Swarmmax,iWith optimal channel layout type, mixing
It closes in particle cluster algorithm, design object is the total heat exchange amount maximized under either simplex condition;Genetic algebra is determined according to total number of channels amount
With per generation population quantity;Cross method is intersected at random for particle preferably individual and population preferably individual;Variation method is a
The random variation of body;
If 3.2) the optimal channel layout result under all operating modes is the same, the channel layout coefficient of each operating mode is 1, is led to
Road layout designs space is 1 design point, which corresponds to the optimal channel layout of either simplex condition;
If 3.3) the optimal channel layout result of each operating mode is different, by coordinating the channel layout coefficient of each operating mode, structure is more
The channel layout design space of operating mode, is as follows:
3.3.1 the channel layout coefficient for) assuming each operating mode is the same, by reducing channel layout coefficient, builds initial channel
Layout designs space, if there is m key Design operating mode, it is assumed that the coefficient of each operating mode is η*, then the heat exchange of i-th of crucial operating mode
Amount design requirement is Qi≥η*Qmax,i, η is made first*=1, according to the heat exchange amount demand under each operating mode, channel layout design section is
Empty set Φ;Then η is reduced with the step delta η of a very little*, i.e. η*=η*- Δ η extends the channel layout design area of each operating mode
Between;If the coordinate design section of multi-state remains as empty set Φ, continue to reduce η with step delta η*, until multi-state is coordinated to lead to
Road layout designs space is not empty set, remembers that channel layout coefficient at this time is η ';
3.3.2 the channel layout coefficient of each operating mode) is improved, further shrinks multi-state lower channel layout designs space;From key
Operating mode 1 arrives crucial operating mode m, increases channel layout coefficient η using small step delta ηi;Under crucial operating mode i, when channel layout coefficient
ηiIncrease to ηi', multi-state coordinate path layout designs space is not empty set Φ;But ηiIncrease to ηi'+Δ η, multi-state association
It is empty set Φ to adjust channel layout design space, then ηi' for the optimal channel location coefficient under crucial operating mode i;
3.4) the optimal channel location coefficient η under each crucial operating mode is determinedi' after, the heat exchange amount range of needs under each operating mode is
Qhot,i≥ηi′Qmax,i, channel layout design space is the collection for meeting all channel layout modes of heat exchange amount demand under each operating mode
It closes.
5. ultra-large type plate-fin heat exchanger under a kind of complex working condition based on integral mean temperature differential method according to claim 1
Channel layout design method, the channel layout optimal design of the multi-state design requirement of the step 4) use following steps:
4.1) build multi-state channel layout majorized function and consider m crucial operating condition design requirement, establish multi-point optimization letter
Number G (x) is as shown in formula (9), the coefficient w of each operating modeiAs shown in formula (10):
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<msub>
<mi>Q</mi>
<mrow>
<mi>h</mi>
<mi>o</mi>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>m</mi>
</munderover>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, x is passage arrangement mode, yiIt is the design parameter of i-th of crucial operating mode.The weight coefficient w of each operating mode is according to each work
The Estimated Time Of Operation t of condition is determined;Consider two special circumstances, when heat exchanger always works at i-th of crucial operating mode, then wi=
1 and wj=0 (1≤j≤m, j ≠ i);When heat exchanger needs frequently to carry out Off-design operation, then heat exchanger is in m crucial operating mode
Under run time it is essentially identical, then
4.2) the overall heat exchange amount that multi-state channel layout mathematical optimization models majorized function is maximization formula (9), design are built
The heat exchange amount constraint for each operating mode that channel layout coefficient determines is constrained to, as shown in formula (11):
Qhot,i(x,yi)≥ηi′Qi,maxi∈[1,m] (11);
4.3) it is passage arrangement mode to carry out multi-state lower channel layout optimization design design parameter using hybrid particle swarm method,
It maximum genetic algebra and is often determined for population invariable number according to number of active lanes, particle individual and the optimal particle individual of history are handed over
Fork realizes the optimization design of channel layout.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810041244.2A CN108090307B (en) | 2018-01-16 | 2018-01-16 | Multi-working-condition plate-fin heat exchanger channel layout design method based on integral average temperature difference method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810041244.2A CN108090307B (en) | 2018-01-16 | 2018-01-16 | Multi-working-condition plate-fin heat exchanger channel layout design method based on integral average temperature difference method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108090307A true CN108090307A (en) | 2018-05-29 |
CN108090307B CN108090307B (en) | 2021-02-19 |
Family
ID=62183038
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810041244.2A Active CN108090307B (en) | 2018-01-16 | 2018-01-16 | Multi-working-condition plate-fin heat exchanger channel layout design method based on integral average temperature difference method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108090307B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109190141A (en) * | 2018-07-11 | 2019-01-11 | 杭州电子科技大学 | A kind of heat exchange efficiency calculation method of Heat Exchangers |
CN110008579A (en) * | 2019-03-29 | 2019-07-12 | 中国原子能科学研究院 | The design method of vertical fins tubing heat exchanger |
CN110160380A (en) * | 2019-06-03 | 2019-08-23 | 中国矿业大学 | A kind of broad passage plate heat exchanger and heat exchanger particle group optimizing construction design method |
CN110779378A (en) * | 2018-07-31 | 2020-02-11 | 中国科学院工程热物理研究所 | Method for intensifying heat exchange |
CN111159903A (en) * | 2019-12-31 | 2020-05-15 | 重庆邮电大学 | Design and manufacturing method of compact type multi-channel multi-fluid heat exchange device |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106539A (en) * | 2011-11-14 | 2013-05-15 | 谢月 | Optimization method of gathering and transportation pipe network system |
CN104657551A (en) * | 2015-02-09 | 2015-05-27 | 浙江大学 | Plate-fin heat exchanger core structure optimization method based on dynamic pixel granularity |
-
2018
- 2018-01-16 CN CN201810041244.2A patent/CN108090307B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106539A (en) * | 2011-11-14 | 2013-05-15 | 谢月 | Optimization method of gathering and transportation pipe network system |
CN104657551A (en) * | 2015-02-09 | 2015-05-27 | 浙江大学 | Plate-fin heat exchanger core structure optimization method based on dynamic pixel granularity |
Non-Patent Citations (7)
Title |
---|
GEORG FIEG 等: "A monogenetic algorithm for optimal design of large-scale heat exchanger", 《CHEMICAL ENGINEERING AND PROCESSING》 * |
M. YOUSEFI 等: "Optimal design of plate-fin heat exchangers by a hybrid evolutionary algorithm", 《INTERNATIONAL COMMUNICATIONS IN HEAT AND MASS TRANSFER》 * |
S. GHOSH 等: "Optimum stacking pattern for multi-stream plate-fin heat exchanger", 《INTERNATIONAL JOURNAL OF THERMAL SCIENCES》 * |
XIANG PENG 等: "Effect of inlet flow maldistribution on the passage arrangement design of multi-stream plate-fin heat exchanger", 《APPLIED THERMAL ENGINEERING》 * |
XING LUO 等: "A hybrid genetic algorithm for synthesis of heat exchanger networks", 《COMPUTERS AND CHEMICAL ENGINEERING》 * |
宋长华: "换热器优化设计中的平均温度差计算", 《重庆市电机工程学会2006年学术会议论文》 * |
彭翔: "复杂产品设计中参数关联和等效简化方法及其应用", 《万方数据知识服务平台》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109190141A (en) * | 2018-07-11 | 2019-01-11 | 杭州电子科技大学 | A kind of heat exchange efficiency calculation method of Heat Exchangers |
CN109190141B (en) * | 2018-07-11 | 2023-04-07 | 杭州电子科技大学 | Heat exchange efficiency calculation method of multi-strand heat exchanger |
CN110779378A (en) * | 2018-07-31 | 2020-02-11 | 中国科学院工程热物理研究所 | Method for intensifying heat exchange |
CN110779378B (en) * | 2018-07-31 | 2021-02-19 | 中国科学院工程热物理研究所 | Method for intensifying heat exchange |
CN110008579A (en) * | 2019-03-29 | 2019-07-12 | 中国原子能科学研究院 | The design method of vertical fins tubing heat exchanger |
CN110008579B (en) * | 2019-03-29 | 2021-06-25 | 中国原子能科学研究院 | Design method of vertical fin tube type heat exchanger |
CN110160380A (en) * | 2019-06-03 | 2019-08-23 | 中国矿业大学 | A kind of broad passage plate heat exchanger and heat exchanger particle group optimizing construction design method |
CN111159903A (en) * | 2019-12-31 | 2020-05-15 | 重庆邮电大学 | Design and manufacturing method of compact type multi-channel multi-fluid heat exchange device |
CN111159903B (en) * | 2019-12-31 | 2023-07-21 | 重庆邮电大学 | Design and manufacturing method of compact multi-channel multi-fluid heat exchange device |
Also Published As
Publication number | Publication date |
---|---|
CN108090307B (en) | 2021-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108090307A (en) | Plate-fin heat exchanger channel layout design method under a kind of multi-state based on integral mean temperature differential method | |
Chu et al. | Study on hydraulic and thermal performance of printed circuit heat transfer surface with distributed airfoil fins | |
Siddiqui et al. | Efficient energy utilization through proper design of microchannel heat exchanger manifolds: A comprehensive review | |
Maharudrayya et al. | Flow distribution and pressure drop in parallel-channel configurations of planar fuel cells | |
CN111222213B (en) | Dynamic simulation method and device for thermodynamic network | |
CN104131983B (en) | Petrochemical industry recirculating cooling water system pump valve optimum combination operating scheme defining method | |
Wang et al. | An investigation of influence factor including different tube bundles on inclined elliptical fin-tube heat exchanger | |
Du et al. | Correlations and optimization of a heat exchanger with offset fins by genetic algorithm combining orthogonal design | |
US20040083012A1 (en) | Method of modeling and sizing a heat exchanger | |
CN102155860B (en) | Method for constructing heat exchange network based on exergy consumption cost | |
CN108842713A (en) | Heavily silt-carrying river water quality and quantity unified collocation and dispatching method and its system | |
CN105864016B (en) | A kind of more water pump distributing system variable water volume runing adjustment methods of open type | |
CN107679294B (en) | A kind of board-like heat exchanger inlet and outlet design method of multichannel | |
CN102305553B (en) | Determination method of total heat transfer coefficient of condenser of thermal generator set | |
CN105302984A (en) | Modeling and simulation method for ground source heat pump set | |
CN106650125A (en) | Method and system for optimizing centrifugal compressor impeller | |
CN109242267A (en) | A kind of crude pipeline transportation * consumption appraisement system analysis method based on analytic hierarchy process (AHP) | |
CN110059372A (en) | A kind of objective design method of the shell-and-tube heat exchanger based on differential evolution algorithm | |
Zhao et al. | Cost-effective optimal design of groundwater source heat pumps | |
CN109522644A (en) | A kind of enhanced heat exchange surface comprehensive performance evaluation method | |
CN116150993A (en) | Heat exchanger design method based on frequency domain characteristics and heat exchanger | |
CN105160062B (en) | A kind of same journey hydraulic pipeline check method | |
CN114358604A (en) | Urban centralized heating system time-sharing and zone-dividing regulation and control method considering building heat storage characteristics | |
CN111928911A (en) | Multi-path wide-range flow measuring device and measuring method thereof | |
CN109883741A (en) | A kind of existing heat exchanger and composed heat-exchange network efficiency evaluation method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |