CN105956329B - The modelling by mechanism calculation method of each channel gain of heat exchanger - Google Patents

The modelling by mechanism calculation method of each channel gain of heat exchanger Download PDF

Info

Publication number
CN105956329B
CN105956329B CN201610369509.2A CN201610369509A CN105956329B CN 105956329 B CN105956329 B CN 105956329B CN 201610369509 A CN201610369509 A CN 201610369509A CN 105956329 B CN105956329 B CN 105956329B
Authority
CN
China
Prior art keywords
formula
heat
heat exchanger
channel gain
gain
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.)
Expired - Fee Related
Application number
CN201610369509.2A
Other languages
Chinese (zh)
Other versions
CN105956329A (en
Inventor
康英伟
归数
归一数
杨平
陈欢乐
李芹
王松
于会群
王念龙
徐春梅
吴周晶
余洁
张经纬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai University of Electric Power
Shanghai Minghua Electric Power Technology and Engineering Co Ltd
Original Assignee
Shanghai University of Electric Power
Shanghai Minghua Electric Power Technology and Engineering Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai University of Electric Power, Shanghai Minghua Electric Power Technology and Engineering Co Ltd filed Critical Shanghai University of Electric Power
Priority to CN201610369509.2A priority Critical patent/CN105956329B/en
Publication of CN105956329A publication Critical patent/CN105956329A/en
Application granted granted Critical
Publication of CN105956329B publication Critical patent/CN105956329B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Feedback Control In General (AREA)

Abstract

The present invention relates to a kind of modelling by mechanism calculation methods of each channel gain of heat exchanger, modelling by mechanism for heat exchanger, or in System Discrimination modeling for determining the value range of each channel gain of heat exchanger, technical process mechanism of the calculation method based on heat exchanger, according to the structural parameters and thermal parameter of heat exchanger, the gain in each channel is calculated.Compared with prior art, the present invention has many advantages, such as simple, efficient.

Description

The modelling by mechanism calculation method of each channel gain of heat exchanger
Technical field
The present invention relates to a kind of heat exchanger modeling techniques, more particularly, to a kind of modelling by mechanism of each channel gain of heat exchanger Calculation method.
Background technique
Heat exchanger is a kind of important industrial equipment, in chemical industry, petroleum, power, food and other many industrial productions It is all widely used, its role is to heat cold medium or cooling thermal medium to suitable technological parameter.In order to design heat exchanger Control system, needs to establish the mathematical model of heat exchanger, and the gain of each process channel is the important parameter of heat exchanger mathematical model.
The method for establishing heat exchanger dynamic mathematical models can be generally divided into two classes: modelling by mechanism method and system identification side Method.Modelling by mechanism method is the physics followed according to heat exchanger technical process, chemistry rule, establishes its mathematical model；System Heat exchanger is then considered as a "black box" by discrimination method, according to the inputoutput data of measurement, is existed according to pre-determined criterion One is selected in Models Sets to coincide best model with data.
It is particularly significant to resolve following two problem for the mathematical model that heat exchanger is established with system identifying method: 1) really Determine the structure of heat-exchanger model；2) range of model parameter to be estimated is determined.Selecting reasonable model structure is to establish accurately The important prerequisite of heat-exchanger model；And properly determine model parameter range, can significantly improve parameter Estimation precision and Speed.Currently, both of these problems are still lacked with perfect solution.
Summary of the invention
It simply, is efficiently changed it is an object of the present invention to overcome the above-mentioned drawbacks of the prior art and provide a kind of The modelling by mechanism calculation method of hot each channel gain of device.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of modelling by mechanism calculation method of each channel gain of heat exchanger, for the modelling by mechanism of heat exchanger, or in system For determining the value range of each channel gain of heat exchanger, technique mistake of the calculation method based on heat exchanger when identification modeling Journey mechanism calculates the gain in each channel according to the structural parameters and thermal parameter of heat exchanger.
The structural parameters and thermal parameter of the heat exchanger include: the physical parameter, heat exchange area, heat exchange system of heat-transfer surface Number；The physical parameter of hot and cold medium；The operating point parameter of heat exchanger, hot and cold rate-of flow and enthalpy.
The physical parameter of the hot and cold medium includes specific heat at constant pressure.
The characteristics of calculation method heat exchanging device, is made as defined below:
(1) the hot and cold medium of heat exchanger is single-phase, is occurred without phase-change；
(2) heat transfer type between hot and cold medium and heat-transfer surface is heat convection.
The technical process of the heat exchanger is as shown in Fig. 1, process channel gain K involved in the calculation methodi(i= 1,2,3,4,5,6,7,8)、KjThe concrete meaning of (j=a, b, c, d) is shown in Table 1 and attached drawing 2.
The calculation method specifically includes the following steps:
Cold medium in heat exchanging device, which distinguishes application quality law of conservation and law of conservation of energy, to be obtained:
Dcs,in-Dcs,out=0 (1)
Dcs,inHcs,in-Dcs,inHcs,out+Qcs=0 (2)
In formula (1) and formula (2), D is rate-of flow, and H is medium specific enthalpy, and Q is heat exchange amount；Lower footnote cs represents cold medium, In represents entrance, and out represents outlet；
Formula (1) is substituted into formula (2), and will be obtained after formula (2) linearisation:
Cp,cs(Tcs,in-Tcs,out)|0ΔDcs,in+Cp,csDcs|0ΔTcs,in+ΔQcs=Cp,csDcs|0ΔTcs,out (3)
In formula, CpFor specific heat at constant pressure, T is temperature；Lower footnote 0 represents operating point；
And then obtain KcCalculation formula:
The quantity of heat convection in formula (2) is calculated according to the cooling formula of Newton:
Qcs=αcsA(Tw-Tcs,out) (5)
In formula, α is the coefficient of heat transfer, and A is heat exchange area；Lower footnote w represents heat-transfer surface；
Formula (5) and formula (1) are substituted into formula (2), and take hot-face temperature TwIt is constant, it will be obtained after formula (2) linearisation:
Cp,cs(Tcs,in-Tcs,out)|0ΔDcs,in+Cp,csDcs|0ΔTcs,in=(Cp,csDcs+αcsA)|0ΔTcs,out (6)
And then obtain K1And K2Calculation formula:
Take hot-face temperature TwIt is constant, it will be obtained after formula (5) linearisation:
ΔQcs=-αcsAΔTcs,out (9)
And then obtain KaCalculation formula:
Thermal medium side is deduced accordingly, the calculation formula of following gain is obtained:
For formula (11) into formula (14), lower footnote hs represents thermal medium；
According to fig. 2, K is obtained5To K8Calculation formula:
K5=K1KaKb (15)
K6=K2KaKb (16)
K7=K3KdKc (17)
K8=K4KdKc (18)
Compared with prior art, present invention is primarily based on the technical process mechanism of heat exchanger, are joined according to the structure of heat exchanger Several and thermal parameter can calculate the gain in each channel, and calculated result is for determining that the range of each channel gain of heat exchanger has There are stronger directive function and practical value.
Detailed description of the invention
Fig. 1 is heat exchanger technical process schematic diagram；
Fig. 2 is each channel gain calculated relationship figure of heat exchanger.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment
The technical process of heat exchanger passes through simultaneously as shown in Figure 1, hot and cold medium respectively flows through the different channels of heat exchanger Heat-transfer surface carries out heat exchange.The purpose of heat exchanger technique, or cold medium is heated using thermal medium, or utilize cold Jie Matter cools down thermal medium.The characteristics of to the heat exchanger considered, is made as defined below:
(1) the hot and cold medium of heat exchanger is single-phase, is occurred without phase-change.
(2) heat transfer type between hot and cold medium and heat-transfer surface is heat convection.
Table 1 lists the main process channel of heat exchanger.From the technical process mechanism of heat exchanger, according to heat exchanger Structural parameters and thermal parameter can calculate the gain in these channels.Fig. 2 gives the calculated relationship figure of each channel gain.
Table 1
The invention proposes the calculation method of each process channel gain of a heat exchanger, the specific implementation steps of this method Are as follows:
Obtain the structural parameters and thermal parameter of heat exchanger.Required parameter mainly has: physical parameter, the heat-transfer surface of heat-transfer surface Product, the coefficient of heat transfer；The physical parameter of hot and cold medium, such as specific heat at constant pressure；The operating point parameter of heat exchanger, hot and cold rate-of flow And enthalpy, etc..
According to formula (4), formula (7), formula (8), formula (10) and formula (11) to formula (18), each process channel of heat exchanger is calculated Gain Ki(i=1,2,3,4,5,6,7,8).

Claims (1)

1. a kind of modelling by mechanism calculation method of each channel gain of heat exchanger, is distinguished for the modelling by mechanism of heat exchanger, or in system For determining the value range of each channel gain of heat exchanger when knowing modeling, which is characterized in that the calculation method is based on heat exchange The technical process mechanism of device calculates the gain in each channel according to the structural parameters and thermal parameter of heat exchanger；
The structural parameters and thermal parameter of the heat exchanger include: the physical parameter, heat exchange area, the coefficient of heat transfer of heat-transfer surface； The physical parameter of hot and cold medium；The operating point parameter of heat exchanger, hot and cold rate-of flow and enthalpy；
The physical parameter of the hot and cold medium includes specific heat at constant pressure；
The characteristics of calculation method heat exchanging device, is made as defined below:
(1) the hot and cold medium of heat exchanger is single-phase, is occurred without phase-change；
(2) heat transfer type between hot and cold medium and heat-transfer surface is heat convection；
The calculation method specifically includes the following steps:
Cold medium in heat exchanging device, which distinguishes application quality law of conservation and law of conservation of energy, to be obtained:
Dcs,in-Dcs,out=0 (1)
Dcs,inHcs,in-Dcs,inHcs,out+Qcs=0 (2)
In formula (1) and formula (2), D is rate-of flow, and H is medium specific enthalpy, and Q is heat exchange amount；Lower footnote cs represents cold medium, in generation Table entry, out represent outlet；
Formula (1) is substituted into formula (2), and will be obtained after formula (2) linearisation:
Cp,cs(Tcs,in-Tcs,out)|0ΔDcs,in+Cp,csDcs|0ΔTcs,in+ΔQcs=Cp,csDcs|0ΔTcs,out (3)
In formula, CpFor specific heat at constant pressure, T is temperature；Lower footnote 0 represents operating point；
And then obtain gain KcCalculation formula:
In formula, KcIt is with QcsFor input variable Tcs,outFor the channel gain of output variable；
The quantity of heat convection in formula (2) is calculated according to the cooling formula of Newton:
Qcs=αcsA(Tw-Tcs,out) (5)
In formula, α is the coefficient of heat transfer, and A is heat exchange area；Lower footnote w represents heat-transfer surface；
Formula (5) and formula (1) are substituted into formula (2), and take hot-face temperature TwIt is constant, it will be obtained after formula (2) linearisation:
Cp,cs(Tcs,in-Tcs,out)|0ΔDcs,in+Cp,csDcs|0ΔTcs,in=(Cp,csDcs+αcsA)|0ΔTcs,out (6)
And then obtain gain K1And K2Calculation formula:
In formula, K1It is with Dcs,inFor input variable Tcs,outFor the channel gain of output variable, K2It is with Tcs,inFor input variable Tcs,outFor the channel gain of output variable；
Take hot-face temperature TwIt is constant, it will be obtained after formula (5) linearisation:
ΔQcs=-αcsAΔTcs,out (9)
And then obtain gain KaCalculation formula:
In formula, KaIt is with Tcs,outFor input variable QcsFor the channel gain of output variable；
Thermal medium side is deduced accordingly, the calculation formula of following gain is obtained:
Formula (11) is into formula (14), KbIt is with QhsFor input variable Ths,outFor the channel gain of output variable, K3It is with Dhs,inFor Input variable Ths,outFor the channel gain of output variable, K4It is with Ths,inFor input variable Ths,outIncrease for the channel of output variable Benefit, KdIt is with Ths,outFor input variable QhsFor the channel gain of output variable；Lower footnote hs represents thermal medium；
Obtain gain K5To K8Calculation formula:
K5=K1KaKb (15)
K6=K2KaKb (16)
K7=K3KdKc (17)
K8=K4KdKc (18)
In formula, K5It is with Dcs,inFor input variable Ths,outFor the channel gain of output variable, K6It is with Tcs,inFor input variable Ths,outFor the channel gain of output variable, K7It is with Dhs,inFor input variable Tcs,outFor the channel gain of output variable, K8It is With Ths,inFor input variable Tcs,outFor the channel gain of output variable.
CN201610369509.2A 2016-05-30 2016-05-30 The modelling by mechanism calculation method of each channel gain of heat exchanger Expired - Fee Related CN105956329B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610369509.2A CN105956329B (en) 2016-05-30 2016-05-30 The modelling by mechanism calculation method of each channel gain of heat exchanger

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610369509.2A CN105956329B (en) 2016-05-30 2016-05-30 The modelling by mechanism calculation method of each channel gain of heat exchanger

Publications (2)

Publication Number Publication Date
CN105956329A CN105956329A (en) 2016-09-21
CN105956329B true CN105956329B (en) 2019-04-09

Family

ID=56909926

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610369509.2A Expired - Fee Related CN105956329B (en) 2016-05-30 2016-05-30 The modelling by mechanism calculation method of each channel gain of heat exchanger

Country Status (1)

Country Link
CN (1) CN105956329B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065515B (en) * 2017-04-17 2020-04-24 河南工程学院 Plate heat exchanger model construction method based on fuzzy PID control
CN107085384B (en) * 2017-04-25 2020-07-31 杭州百子尖科技股份有限公司 Heat exchanger modeling method and device for dynamic simulation system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5177489A (en) * 1989-09-26 1993-01-05 Magnavox Electronic Systems Company Pseudolite-aided method for precision kinematic positioning
US5317514A (en) * 1992-05-29 1994-05-31 Alliedsignal Inc. Integrity monitoring of navigation systems using Baye's rule
CN103822758A (en) * 2014-03-06 2014-05-28 中国石油大学(北京) Online diagnosis and selective control method and device for leakage current unusual service conditions of heat exchanger
CN104536292A (en) * 2014-12-05 2015-04-22 北京航空航天大学 Method for conducting fault diagnosis on heat exchanger of aircraft environmental control system based on STF (Strong Tracking Filter) and MB

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5177489A (en) * 1989-09-26 1993-01-05 Magnavox Electronic Systems Company Pseudolite-aided method for precision kinematic positioning
US5317514A (en) * 1992-05-29 1994-05-31 Alliedsignal Inc. Integrity monitoring of navigation systems using Baye's rule
CN103822758A (en) * 2014-03-06 2014-05-28 中国石油大学(北京) Online diagnosis and selective control method and device for leakage current unusual service conditions of heat exchanger
CN104536292A (en) * 2014-12-05 2015-04-22 北京航空航天大学 Method for conducting fault diagnosis on heat exchanger of aircraft environmental control system based on STF (Strong Tracking Filter) and MB

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《A High Precision Simulation Model for Single Phase Heated Tubes of Power Plant Boilers》;YW Kang;《Springer Berlin Heidelberg》;20141231;第81-86页
《A reduced 1D dynamic model of a planar direct internal reforming solid oxide fuel cell for system research》;YW Kang;《Journal of Power Sources》;20091231;第170-176页
《电站锅炉过热器的集总参数动态建模与仿真》;康英伟 等;《计算机仿真》;20120930;第29卷(第9期);第332-334、342页

Also Published As

Publication number Publication date
CN105956329A (en) 2016-09-21

Similar Documents

Publication Publication Date Title
Jamshidi et al. Experimental analysis of heat transfer enhancement in shell and helical tube heat exchangers
Zhang et al. Modeling and performance analysis of a two-phase thermosyphon loop with partially/fully liquid-filled downcomer
US20170130982A1 (en) Inferential sensor for internal heat exchanger parameters
Gao et al. Dynamic analysis of cross flow heat exchangers in data centers using transient effectiveness method
Carl et al. The theoretical and experimental investigation of the heat transfer process of an automobile radiator
CN105956329B (en) The modelling by mechanism calculation method of each channel gain of heat exchanger
Zhou et al. Prediction of temperature distribution in shell-and-tube heat exchangers
Du et al. Experimental study and numerical simulation of flow and heat transfer performance on an offset plate-fin heat exchanger
CN105159138B (en) A kind of radiometer temperature control system heat transfer mathematical model establishing method
Kishan et al. CFD Analysis of heat exchanger models design using ansys fluent
Torregrosa-Jaime et al. Thermal characterisation of compact heat exchangers for air heating and cooling in electric vehicles
Kammuang-Lue et al. Correlation to predict the maximum heat flux of a vertical closed-loop pulsating heat pipe
Ranjbar et al. The effect of geometrical parameters on heat transfer coefficient in helical double tube exchangers
Baymetova et al. Optimization of the developed outer surface of an industrial oil cooler
Atmaca et al. Investigation of transient behaviour of combi boiler type appliances for domestic hot water
Alam et al. Experimental and numerical investigation of an air to water heat exchanger
Chen et al. Thermal sizing design and experimental evaluation of molten salt-to-air heat exchanger
Subhani et al. The Systematic Comparision on Analysis of Parallel Flow and Counter Flow Heat Exchanger by using CFD and Practicle Methods
Lubis et al. Design of a heat exchanger of three concentric tube layer on contrary flow
Chopra et al. Thermal performance analysis of cross-flow unmixed-unmixed heat exchanger by the variation of inlet condition of hot fluid
Ahamad et al. Experimental investigation of three fluid heat exchangers using roughness on the outer surface of a helical coil
Sarmiento et al. A numerical method for shell and thermosyphon heat exchanger analysis
Patil et al. CFD analysis of heat transfer enhancement in double pipe heat exchanger using single and double triangular tape as inserts
Turki et al. Experimentally validated bond graph model of a brazed-plate heat exchanger (BPHE)
Tanish et al. Experimental study of heat transfer in concentric triple pipe heat exchanger

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190409

Termination date: 20200530

CF01 Termination of patent right due to non-payment of annual fee