CN111597679B - Dynamic calculation method for external characteristic parameters of absorption heat pump for comprehensive energy network - Google Patents
Dynamic calculation method for external characteristic parameters of absorption heat pump for comprehensive energy network Download PDFInfo
- Publication number
- CN111597679B CN111597679B CN202010261434.2A CN202010261434A CN111597679B CN 111597679 B CN111597679 B CN 111597679B CN 202010261434 A CN202010261434 A CN 202010261434A CN 111597679 B CN111597679 B CN 111597679B
- Authority
- CN
- China
- Prior art keywords
- temperature
- refrigerant
- pressure
- calculating
- evaporator
- 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
Links
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 37
- 238000004364 calculation method Methods 0.000 title claims abstract description 37
- 239000003507 refrigerant Substances 0.000 claims abstract description 55
- 238000000034 method Methods 0.000 claims abstract description 32
- 230000000704 physical effect Effects 0.000 claims abstract description 16
- AMXOYNBUYSYVKV-UHFFFAOYSA-M lithium bromide Chemical compound [Li+].[Br-] AMXOYNBUYSYVKV-UHFFFAOYSA-M 0.000 claims description 72
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 57
- 239000006096 absorbing agent Substances 0.000 claims description 41
- 239000000498 cooling water Substances 0.000 claims description 25
- 238000001704 evaporation Methods 0.000 claims description 23
- 230000008020 evaporation Effects 0.000 claims description 20
- 238000009834 vaporization Methods 0.000 claims description 12
- 230000008016 vaporization Effects 0.000 claims description 12
- 238000009833 condensation Methods 0.000 claims description 4
- 230000005494 condensation Effects 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 239000002826 coolant Substances 0.000 claims description 3
- 238000005057 refrigeration Methods 0.000 claims description 3
- 230000008929 regeneration Effects 0.000 claims description 3
- 238000011069 regeneration method Methods 0.000 claims description 3
- 238000005094 computer simulation Methods 0.000 abstract description 10
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Sorption Type Refrigeration Machines (AREA)
Abstract
The invention relates to a dynamic calculation method for external characteristic parameters of an absorption heat pump for an integrated energy network, and belongs to the technical field of digital simulation of the integrated energy network. The method analyzes the dynamic characteristic of the absorption heat pump, ignores the dynamic process of the fast process through the comparison of the response time and the response time, simplifies the dynamic partial differential equation of the fast process into an algebraic equation, flexibly uses the physical property table look-up function of the refrigerant, reduces the complexity of the partial differential equation of the system dynamic model, avoids the complicated calculation of the physical property parameters of the refrigerant, ensures that the model not only keeps the external dynamic characteristic of the absorption heat pump, but also reduces the number of the partial differential equations for the dynamic modeling of the system, greatly reduces the calculation time and meets the requirement of calculating the real-time property in the dynamic simulation process of the comprehensive energy system. Therefore, the method is a good operation parameter calculation method of the absorption heat pump suitable for the system dynamic simulation of the comprehensive energy system.
Description
Technical Field
The invention relates to a dynamic calculation method for external characteristic parameters of an absorption heat pump for an integrated energy network, and belongs to the technical field of digital simulation of the integrated energy network.
Background
The comprehensive energy system covers different types of energy production and conversion equipment such as cold, heat and electricity, the absorption heat pump is low in energy consumption and taste, the utilization method is flexible, the proportion of the absorption heat pump in the comprehensive energy system utilizing waste heat is larger and larger, and dynamic simulation modeling of the external characteristics of the absorption heat pump plays an important role in simulation (control) of the comprehensive energy system.
The existing dynamic simulation modeling of the absorption heat pump is mainly based on the structural principle, utilizes complex differential equations, partial differential equations and algebraic equations to solve simultaneously, has high latitude and high nonlinear characteristics, is high in solving difficulty and high in calculation cost, and cannot meet the real-time requirement of dynamic simulation calculation of an integrated energy system. And the method is lack of a uniform processing means for factors such as iteration of working substances on physical properties, high latitude of partial differential equations caused by complex systems, and involved in complex phase change heat exchange.
Disclosure of Invention
The invention aims to provide a method for dynamic simulation modeling of external characteristics of an absorption heat pump, which is characterized in that the dynamic characteristics of the operation of a system are analyzed, the time of a dynamic response link is compared, the quick response process is ignored, a main dynamic link is seized, the number of partial differential equations for dynamic modeling of the system is reduced, iterative calculation is reduced, and the real-time requirement of dynamic simulation calculation of a comprehensive energy system is met.
The invention provides a dynamic calculation method for external characteristic parameters of an absorption heat pump of a comprehensive energy network, which comprises the following steps:
setting the calculation step length of the external characteristic parameter of the absorption heat pump as delta t, and initializing the absorptionTemperature t of dilute lithium bromide solution in receiverxsInitializing the number of calculation cycles i to 0;
(1) according to the lithium bromide solution temperature t in the time step of t-delta t calculation of the generator in the absorption heat pumpfsqAnd driving the heat source steam outlet temperature tqdcCalculating the heat exchange quantity Q of the driving heat source steam to the lithium bromide concentrated solution in the generator by using the following formulafsq：
In the formula, wqdTo drive the mass flow of heat source steam, alphafsTo drive the equivalent heat transfer coefficient, alpha, of the heat source steam and the lithium bromide concentrated solutionfsHas a value in the range of 2 to 15, AfsIs the heat exchange area of the generator, tqdjDriving the heat source steam inlet temperature;
(2) q calculated according to step (1)fsqCalculating and updating t + delta t by using the following formula, and calculating the temperature t of the steam outlet of the driving heat source in the time stepqdj：
In the formula, ρwTo drive the density of the heat source vapor, cwSpecific heat, V, for driving heat source steamgIs the inner volume of a heat exchange tube bundle in the generator;
(3) calculating the pressure P in the generator in time step according to t-delta tfsqCalculating a function f using the physical property parameter of the water vapor1(P) determining the pressure PfsqLower corresponding water vapor saturation temperature tsatThe temperature t of the lithium bromide concentrated solution in the generator is obtained by the following formulafsq：
tfsq＝124.937-7.71649*100*ξ+0.152286*104*ξ2-7.9509*102*ξ3+(-2.00775+16.976*ξ-31.33362*ξ2-795.09*ξ3)
In the formula, xi is the lithium bromide concentrated solution in the generatorLiquid concentration, f1(p) is a function of the saturation temperature of the water vapour as a function of the pressure;
(4) according to the temperature t of the lithium bromide concentrated solution in the step (3)fsqAnd t- Δ t calculating the generator pressure P in time stepfsqLooking up a function f by using the physical property parameters of the water vapor2(p, t), calculating the enthalpy h of the water vaporfsThe evaporation amount w of the refrigerant vapor in the generator is obtained by the following equationzzzf：
wzzzf＝Qfsq/hfs
Wherein f is2(p, t) is a table look-up function for solving the enthalpy of the water vapor according to the pressure and the temperature;
(5) calculating the inlet temperature and the outlet temperature t of the circulating cooling water of the condenser in the absorption heat pump in the time step according to the t-delta tlcol1And tlcol2The condensing temperature t of the condenser is determined by the following equationln：
In the formula,. DELTA.tlndIs the heat exchange end difference of the condenser, delta tlndThe value of (a) is 2 to 3;
(6) the condensing temperature t of the condenser obtained in the step (5)lnLooking up the function f from the physical property parameters of the water vapor3(t) obtaining the refrigerant water vapor pressure in the condenser, and updating the refrigerant water vapor pressure p in the generator at the current calculation time step by using the following formulafsq：
In the formula,. DELTA.pzlFor the design of the pressure loss from the generator to the condenser under operating conditions, wzzzf0For the evaporation capacity of refrigerant steam under design conditions, f3(t) solving a table look-up function of the saturation pressure of the water vapor according to the temperature;
(7) the condensation temperature t obtained according to the step (5)lnCalculating internal condenser in time step by combining t-delta tInlet temperature t of intermediate cooling waterlcol1The cooling water outlet temperature t in the condenser is updated using the following formulalcol2：
In the formula, wcolFor cooling water flow, VlnIs the internal volume of a heat exchange tube bundle in a condenser, f4A physical property lookup function for looking up the latent heat of vaporization of water vapor based on the pressure, hrlnThe latent heat of vaporization of the water vapor under the corresponding current pressure;
(8) according to the self-balancing dynamic process of the throttle valve in the process that the refrigerant flows from the condenser to the evaporator, the evaporation amount of the refrigerant in the evaporator in the time step is calculated by using the following formula to obtain t + delta t
In the formula (I), the compound is shown in the specification,
and
calculating the amount of refrigerant evaporated in the evaporator in time steps for t + Δ t and t- Δ t, respectively, tczfIs the dynamic response time constant, t, of the refrigerant passing from the condenser to the evaporatorczfIs 10-50;
(9) combined with refrigerant water inlet temperature t in the evaporatorlm1And t- Δ t calculating refrigerant water outlet temperature in evaporator in time step
The latent heat of vaporization h of the refrigerant steam of the evaporator under the current pressure is obtained by solving the following simultaneous equationrzfRefrigeration workRate QclAnd t + delta t calculating the outlet temperature of the coolant water in the time step
In the formula,. DELTA.tlmThe phase change heat exchange end difference in the evaporator is 2-3, tzfIs the evaporation temperature, p, of the refrigerant in the evaporatorzfIs the evaporation pressure in the evaporator, hrzfFor pressure p corresponding to water vapour in the evaporatorzfLatent heat of vaporization, VlmIs the inner volume of a refrigerant water heat exchange tube bundle, wlmMass flow rate of refrigerant water;
(10) evaporating pressure p in the evaporator according to step (9)zfThe refrigerant water vapor pressure p in the absorber is calculated using the following formulaxs：
In the formula, wzf0To the evaporation capacity, Δ p, of refrigerant in the evaporator under design conditionszfxsThe pressure loss from the evaporator to the absorber under the design working condition is realized;
(11) according to the temperature t of the lithium bromide dilute solution in the absorberxsThe specific heat c of the dilute lithium bromide solution in the absorber was determined by the following equationxs：
cxs＝4.1868*(0.9928285+txs*(-3.18742*10-5-3.0105*10-6*txs)+(-1.3169179+txs*(2.9856*10-3-1.7172*10-6*txs))*ξ1+(0.6481006+txs*(-4.0198*10-3+8.3641*10-6*txs))*ξ1 2)
In the formula, xi1Is the concentration of the dilute lithium bromide solution in the absorber;
(12) calculating the temperature t of the lithium bromide dilute solution in the absorber in the time step according to the t-delta t in the step (13)xsCirculating cooling water outlet temperature tlcol1Boundary-given recirculated cooling water inlet temperature tcoljThe heat Q obtained by absorbing refrigerant vapor with dilute lithium bromide solution in the absorber is determined by the following formulaxs：
In the formula, alphaxsIs the heat transfer coefficient between the absorber solution and the cooling water, AxsIs the heat exchange area;
(13) the heat exchange power Q obtained according to the step (12)xsAnd the specific heat c of the solution in the absorber calculated in the step (11)xsThe solution temperature t in the regeneration absorber is calculated using the following formulaxsAnd the outlet temperature t of the circulating cooling water in the absorberlcol1：
In the formula, MxsThe mass of the solution in the absorber;
(14) and (4) returning to the step (1), and entering t +2 delta t time step cycle calculation to realize dynamic calculation of external characteristic parameters of the absorption heat pump for the comprehensive energy network.
The invention provides a dynamic calculation method for external characteristic parameters of an absorption heat pump for a comprehensive energy network, which has the characteristics and advantages that:
the method is used for modeling the external characteristic of the absorption heat pump meeting the real-time requirement of the dynamic simulation of the comprehensive energy network, and the dynamic characteristic of the absorption heat pump mainly comprises the dynamic characteristic of a generator, the dynamic process of phase change evaporation of the generator, the condensation process of a condenser, the phase change evaporation process of an evaporator, the dynamic process of heat absorption of an absorber and the like from the external characteristic, so that the dynamic process of phase change heat exchange is neglected according to the dynamic response time, the dynamic process of heat exchange of the generator and the absorber and the dynamic process of heat exchange of cooling water in the condenser and the evaporator are mainly considered, a partial differential equation of the phase change heat exchange is changed into an algebraic equation, and the method is characterized in that the method comprises the following steps: the generator, the condenser, the evaporator and the absorber are sequentially used for establishing a dynamic model.
The dynamic calculation method for the external characteristic parameters of the absorption heat pump for the comprehensive energy network analyzes the dynamic characteristics of the absorption heat pump, ignores the dynamic process of the fast process through comparison of response time speed, simplifies the dynamic partial differential equation of the fast process into an algebraic equation, flexibly uses a physical property table look-up function of a refrigerant, reduces the complexity of the partial differential equation of a system dynamic model, avoids complicated calculation of the physical property parameters of the refrigerant, ensures that the model not only keeps the external dynamic characteristics of the absorption heat pump, but also reduces the number of the partial differential equations of system dynamic modeling, greatly reduces the calculation time, and meets the requirement of calculating the real-time performance in the dynamic simulation process of the comprehensive energy system. Therefore, the method is a good operation parameter calculation method of the absorption heat pump suitable for the system dynamic simulation of the comprehensive energy system.
Drawings
Figure 1 is a schematic diagram of the construction of an absorption heat pump involved in the process of the present invention.
Detailed Description
The invention provides a dynamic calculation method for external characteristic parameters of an absorption heat pump of an integrated energy network, wherein the related absorption heat pump has a structural schematic diagram shown in figure 1, and the method comprises the following steps:
setting the calculation step length of the external characteristic parameter of the absorption heat pump as delta t, and initializing the temperature t of the lithium bromide dilute solution in the absorberxsInitializing the number of calculation cycles i to 0;
(1) according to the lithium bromide solution temperature t in the time step of t-delta t calculation of the generator in the absorption heat pumpfsqAnd driving the heat source steam outlet temperature tqdcCalculating the heat exchange quantity Q of the driving heat source steam to the lithium bromide concentrated solution in the generator by using the following formulafsq：
In the formula, wqdTo drive the mass flow of heat source steam, alphafsTo drive the equivalent heat transfer coefficient, alpha, of the heat source steam and the lithium bromide concentrated solutionfsHas a value in the range of 2 to 15, AfsIs the heat exchange area of the generator, tqdjDriving the heat source steam inlet temperature;
(2) q calculated according to step (1)fsqCalculating and updating t + delta t by using the following formula, and calculating the temperature t of the steam outlet of the driving heat source in the time stepqdj：
In the formula, ρwTo drive the density of the heat source vapor, cwSpecific heat, V, for driving heat source steamgIs the inner volume of a heat exchange tube bundle in the generator;
(3) calculating the pressure P in the generator in time step according to t-delta tfsqCalculating a function f using the physical property parameter of the water vapor1(P) determining the pressure PfsqLower corresponding water vapor saturation temperature tsatThe temperature t of the lithium bromide concentrated solution in the generator is obtained by the following formulafsq：
tfsq＝124.937-7.71649*100*ξ+0.152286*104*ξ2-7.9509*102*ξ3+(-2.00775+16.976*ξ-31.33362*ξ2-795.09*ξ3)
Where xi is the concentration of the lithium bromide solution in the generator, f1(p) is a function of the saturation temperature of the water vapour as a function of the pressure;
(4) according to the temperature t of the lithium bromide concentrated solution in the step (3)fsqAnd t- Δ t calculating the generator pressure P in time stepfsqLooking up a function f by using the physical property parameters of the water vapor2(p, t), calculating the enthalpy h of the water vaporfsThe evaporation amount w of the refrigerant vapor in the generator is obtained by the following equationzzzf：
wzzzf＝Qfsq/hfs
Wherein f is2(p, t) is a table look-up function for solving the enthalpy of the water vapor according to the pressure and the temperature;
(5) calculating the inlet temperature and the outlet temperature t of the circulating cooling water of the condenser in the absorption heat pump in the time step according to the t-delta tlcol1And tlcol2The condensing temperature t of the condenser is determined by the following equationln：
In the formula,. DELTA.tlndIs the heat exchange end difference of the condenser, delta tlndThe value of (a) is 2 to 3;
(6) the condensing temperature t of the condenser obtained in the step (5)lnLooking up the function f from the physical property parameters of the water vapor3(t) obtaining the refrigerant water vapor pressure in the condenser, and updating the refrigerant water vapor pressure p in the generator at the current calculation time step by using the following formulafsq：
In the formula,. DELTA.pzlFor the design of the pressure loss from the generator to the condenser under operating conditions, wzzzf0For the evaporation capacity of refrigerant steam under design conditions, f3(t) solving a table look-up function of the saturation pressure of the water vapor according to the temperature;
(7) the condensation temperature t obtained according to the step (5)lnAnd calculating the inlet temperature t of cooling water in the condenser in the time step by combining t-delta tlcol1The cooling water outlet temperature t in the condenser is updated using the following formulalcol2：
In the formula, wcolFor cooling water flow, VlnIs the internal volume of a heat exchange tube bundle in a condenser, f4For investigating water vapour by pressurePhysical property lookup function of latent heat of vaporization, hrlnThe latent heat of vaporization of the water vapor under the corresponding current pressure;
(8) according to the self-balancing dynamic process of the throttle valve in the process that the refrigerant flows from the condenser to the evaporator, the evaporation amount of the refrigerant in the evaporator in the time step is calculated by using the following formula to obtain t + delta t
In the formula (I), the compound is shown in the specification,
and
calculating the amount of refrigerant evaporated in the evaporator in time steps for t + Δ t and t- Δ t, respectively, tczfIs the dynamic response time constant, t, of the refrigerant passing from the condenser to the evaporatorczfIs 10-50;
(9) combined with refrigerant water inlet temperature t in the evaporatorlm1And t- Δ t calculating refrigerant water outlet temperature in evaporator in time step
The latent heat of vaporization h of the refrigerant steam of the evaporator under the current pressure is obtained by solving the following simultaneous equationrzfRefrigeration power QclAnd t + delta t calculating the outlet temperature of the coolant water in the time step
In the formula,. DELTA.tlmIs the phase change heat exchange end difference and phase in the evaporatorThe value of the heat transfer end difference is 2-3, tzfIs the evaporation temperature, p, of the refrigerant in the evaporatorzfIs the evaporation pressure in the evaporator, hrzfFor pressure p corresponding to water vapour in the evaporatorzfLatent heat of vaporization, VlmIs the inner volume of a refrigerant water heat exchange tube bundle, wlmMass flow rate of refrigerant water;
(10) evaporating pressure p in the evaporator according to step (9)zfThe refrigerant water vapor pressure p in the absorber is calculated using the following formulaxs：
In the formula, wzf0To the evaporation capacity, Δ p, of refrigerant in the evaporator under design conditionszfxsThe pressure loss from the evaporator to the absorber under the design working condition is realized;
(11) according to the temperature t of the lithium bromide dilute solution in the absorberxsThe specific heat c of the dilute lithium bromide solution in the absorber was determined by the following equationxs：
cxs＝4.1868*(0.9928285+txs*(-3.18742*10-5-3.0105*10-6*txs)+(-1.3169179+txs*(2.9856*10-3-1.7172*10-6*txs))*ξ1+(0.6481006+txs*(-4.0198*10-3+8.3641*10-6*txs))*ξ1 2)
In the formula, xi1Is the concentration of the dilute lithium bromide solution in the absorber;
(12) calculating the temperature t of the lithium bromide dilute solution in the absorber in the time step according to the t-delta t in the step (13)xsCirculating cooling water outlet temperature tlcol1Boundary-given recirculated cooling water inlet temperature tcoljThe heat Q obtained by absorbing refrigerant vapor with dilute lithium bromide solution in the absorber is determined by the following formulaxs：
In the formula, alphaxsIs the heat transfer coefficient between the absorber solution and the cooling water, AxsIs the heat exchange area;
(13) the heat exchange power Q obtained according to the step (12)xsAnd the specific heat c of the solution in the absorber calculated in the step (11)xsThe solution temperature t in the regeneration absorber is calculated using the following formulaxsAnd the outlet temperature t of the circulating cooling water in the absorberlcol1：
In the formula, MxsThe mass of the solution in the absorber;
(14) and (4) returning to the step (1), and entering t +2 delta t time step cycle calculation to realize dynamic calculation of external characteristic parameters of the absorption heat pump for the comprehensive energy network.
Claims (1)
1. A dynamic calculation method for external characteristic parameters of an absorption heat pump of an integrated energy network is characterized by comprising the following steps:
setting the calculation step length of the external characteristic parameter of the absorption heat pump as delta t, and initializing the temperature t of the lithium bromide dilute solution in the absorberxsInitializing the number of calculation cycles i to 0;
(1) according to the lithium bromide solution temperature t in the time step of t-delta t calculation of the generator in the absorption heat pumpfsqAnd driving the heat source steam outlet temperature tqdcCalculating the heat exchange quantity Q of the driving heat source steam to the lithium bromide concentrated solution in the generator by using the following formulafsq：
In the formula, wqdTo drive the mass flow of heat source steam, alphafsTo drive the equivalent heat transfer coefficient, alpha, of the heat source steam and the lithium bromide concentrated solutionfsHas a value in the range of 2 to 15, AfsFor heat exchange area of generator，tqdjDriving the heat source steam inlet temperature;
(2) q calculated according to step (1)fsqCalculating and updating t + delta t by using the following formula, and calculating the temperature t of the steam outlet of the driving heat source in the time stepqdj：
In the formula, ρwTo drive the density of the heat source vapor, cwSpecific heat, V, for driving heat source steamgIs the inner volume of a heat exchange tube bundle in the generator;
(3) calculating the pressure P in the generator in time step according to t-delta tfsqCalculating a function f using the physical property parameter of the water vapor1(P) determining the pressure PfsqLower corresponding water vapor saturation temperature tsatThe temperature t of the lithium bromide concentrated solution in the generator is obtained by the following formulafsq：
tfsq＝124.937-7.71649*100*ξ+0.152286*104*ξ2-7.9509*102*ξ3+(-2.00775+16.976*ξ-31.33362*ξ2-795.09*ξ3)
Where xi is the concentration of the lithium bromide solution in the generator, f1(p) is a function of the saturation temperature of the water vapour as a function of the pressure;
(4) according to the temperature t of the lithium bromide concentrated solution in the step (3)fsqAnd t- Δ t calculating the generator pressure P in time stepfsqLooking up a function f by using the physical property parameters of the water vapor2(p, t), calculating the enthalpy h of the water vaporfsThe evaporation amount w of the refrigerant vapor in the generator is obtained by the following equationzzzf：
wzzzf＝Qfsq/hfs
Wherein f is2(p, t) is a table look-up function for solving the enthalpy of the water vapor according to the pressure and the temperature;
(5) calculating the inlet temperature and the outlet temperature t of the circulating cooling water of the condenser in the absorption heat pump in the time step according to the t-delta tlcol1And tlcol2The cold is determined by the following formulaCondensation temperature t of the condenserln：
In the formula,. DELTA.tlndIs the heat exchange end difference of the condenser, delta tlndThe value of (a) is 2 to 3;
(6) the condensing temperature t of the condenser obtained in the step (5)lnLooking up the function f from the physical property parameters of the water vapor3(t) obtaining the refrigerant water vapor pressure in the condenser, and updating the refrigerant water vapor pressure p in the generator at the current calculation time step by using the following formulafsq：
In the formula,. DELTA.pzlFor the design of the pressure loss from the generator to the condenser under operating conditions, wzzzf0For the evaporation capacity of refrigerant steam under design conditions, f3(t) solving a table look-up function of the saturation pressure of the water vapor according to the temperature;
(7) the condensation temperature t obtained according to the step (5)lnAnd calculating the inlet temperature t of cooling water in the condenser in the time step by combining t-delta tlcol1The cooling water outlet temperature t in the condenser is updated using the following formulalcol2：
In the formula, wcolFor cooling water flow, VlnIs the internal volume of a heat exchange tube bundle in a condenser, f4A physical property lookup function for looking up the latent heat of vaporization of water vapor based on the pressure, hrlnThe latent heat of vaporization of the water vapor under the corresponding current pressure;
(8) the evaporation of refrigerant in the evaporator in a time step is calculated using the following formula to obtain t + Δ t from the self-balancing dynamic process of the throttle valve during the refrigerant flow from the condenser to the evaporatorAmount of hair
In the formula (I), the compound is shown in the specification,
and
calculating the amount of refrigerant evaporated in the evaporator in time steps for t + Δ t and t- Δ t, respectively, tczfIs the dynamic response time constant, t, of the refrigerant passing from the condenser to the evaporatorczfIs 10-50;
(9) combined with refrigerant water inlet temperature t in the evaporatorlm1And t- Δ t calculating refrigerant water outlet temperature in evaporator in time step
The latent heat of vaporization h of the refrigerant steam of the evaporator under the current pressure is obtained by solving the following simultaneous equationrzfRefrigeration power QclAnd t + delta t calculating the outlet temperature of the coolant water in the time step
In the formula,. DELTA.tlmThe phase change heat exchange end difference in the evaporator is 2-3, tzfIs the evaporation temperature, p, of the refrigerant in the evaporatorzfIs the evaporation pressure in the evaporator, hrzfFor pressure p corresponding to water vapour in the evaporatorzfLatent heat of vaporization, VlmIs the inner volume of a refrigerant water heat exchange tube bundle, wlmMass flow rate of refrigerant water;
(10) evaporating pressure p in the evaporator according to step (9)zfThe refrigerant water vapor pressure p in the absorber is calculated using the following formulaxs：
In the formula, wzf0To the evaporation capacity, Δ p, of refrigerant in the evaporator under design conditionszfxsThe pressure loss from the evaporator to the absorber under the design working condition is realized;
(11) according to the temperature t of the lithium bromide dilute solution in the absorberxsThe specific heat c of the dilute lithium bromide solution in the absorber was determined by the following equationxs：
cxs＝4.1868*(0.9928285+txs*(-3.18742*10-5-3.0105*10-6*txs)+(-1.3169179+txs*(2.9856*10-3-1.7172*10-6*txs))*ξ1+(0.6481006+txs*(-4.0198*10-3+8.3641*10-6*txs))*ξ1 2)
In the formula, xi1Is the concentration of the dilute lithium bromide solution in the absorber;
(12) calculating the temperature t of the lithium bromide dilute solution in the absorber in the time step according to the t-delta t in the step (13)xsCirculating cooling water outlet temperature tlcol1Boundary-given recirculated cooling water inlet temperature tcoljThe heat Q obtained by absorbing refrigerant vapor with dilute lithium bromide solution in the absorber is determined by the following formulaxs：
In the formula, alphaxsIs the heat transfer coefficient between the absorber solution and the cooling water, AxsIs the heat exchange area;
(13) the heat exchange power Q obtained according to the step (12)xsAnd the specific heat c of the solution in the absorber calculated in the step (11)xsThe solution temperature t in the regeneration absorber is calculated using the following formulaxsAnd the outlet temperature t of the circulating cooling water in the absorberlcol1：
In the formula, MxsThe mass of the solution in the absorber;
(14) and (4) returning to the step (1), and entering t +2 delta t time step cycle calculation to realize dynamic calculation of external characteristic parameters of the absorption heat pump for the comprehensive energy network.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010261434.2A CN111597679B (en) | 2020-04-03 | 2020-04-03 | Dynamic calculation method for external characteristic parameters of absorption heat pump for comprehensive energy network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010261434.2A CN111597679B (en) | 2020-04-03 | 2020-04-03 | Dynamic calculation method for external characteristic parameters of absorption heat pump for comprehensive energy network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111597679A CN111597679A (en) | 2020-08-28 |
| CN111597679B true CN111597679B (en) | 2021-06-22 |
Family
ID=72181770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010261434.2A Expired - Fee Related CN111597679B (en) | 2020-04-03 | 2020-04-03 | Dynamic calculation method for external characteristic parameters of absorption heat pump for comprehensive energy network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111597679B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003090655A (en) * | 2001-09-14 | 2003-03-28 | Nippon Steel Corp | System and method for simulating suction heat pump and recording medium |
| CN106529021A (en) * | 2016-11-09 | 2017-03-22 | 东南大学 | Air conditioning system simulation method based on feature recognition |
| CN107016187A (en) * | 2017-03-31 | 2017-08-04 | 天津大学 | A kind of optimization method for determining fume afterheat organic rankine cycle system parameter |
| CN107025519A (en) * | 2017-03-30 | 2017-08-08 | 厦门大学 | Area distribution formula cool and thermal power polygenerations systeme MIXED INTEGER nonlinear model optimization method |
| CN108021729A (en) * | 2017-03-07 | 2018-05-11 | 林振娴 | A kind of computational methods of definite hot water two-period form single-effect lithiumbromide chiller heat force characteristic |
| CN108019978A (en) * | 2017-03-07 | 2018-05-11 | 林振娴 | A kind of computational methods of definite steam single-effect lithiumbromide refrigeration mechanism cold |
| CN109711614A (en) * | 2018-12-24 | 2019-05-03 | 新奥数能科技有限公司 | A kind of the dynamic optimization progress control method and system of distributed busbar protection |
| CN109886523A (en) * | 2018-12-25 | 2019-06-14 | 清华大学 | A kind of comprehensive energy net dynamic model multi tate calculation method |
| CN110458353A (en) * | 2019-08-08 | 2019-11-15 | 上海交通大学 | Electric heating combines microgrid energy step optimization method and system |
| CN110610037A (en) * | 2019-08-30 | 2019-12-24 | 中国科学院工程热物理研究所 | Thermodynamic simulation method for counter-flow heat regenerator |
-
2020
- 2020-04-03 CN CN202010261434.2A patent/CN111597679B/en not_active Expired - Fee Related
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003090655A (en) * | 2001-09-14 | 2003-03-28 | Nippon Steel Corp | System and method for simulating suction heat pump and recording medium |
| CN106529021A (en) * | 2016-11-09 | 2017-03-22 | 东南大学 | Air conditioning system simulation method based on feature recognition |
| CN108021729A (en) * | 2017-03-07 | 2018-05-11 | 林振娴 | A kind of computational methods of definite hot water two-period form single-effect lithiumbromide chiller heat force characteristic |
| CN108019978A (en) * | 2017-03-07 | 2018-05-11 | 林振娴 | A kind of computational methods of definite steam single-effect lithiumbromide refrigeration mechanism cold |
| CN107025519A (en) * | 2017-03-30 | 2017-08-08 | 厦门大学 | Area distribution formula cool and thermal power polygenerations systeme MIXED INTEGER nonlinear model optimization method |
| CN107016187A (en) * | 2017-03-31 | 2017-08-04 | 天津大学 | A kind of optimization method for determining fume afterheat organic rankine cycle system parameter |
| CN109711614A (en) * | 2018-12-24 | 2019-05-03 | 新奥数能科技有限公司 | A kind of the dynamic optimization progress control method and system of distributed busbar protection |
| CN109886523A (en) * | 2018-12-25 | 2019-06-14 | 清华大学 | A kind of comprehensive energy net dynamic model multi tate calculation method |
| CN110458353A (en) * | 2019-08-08 | 2019-11-15 | 上海交通大学 | Electric heating combines microgrid energy step optimization method and system |
| CN110610037A (en) * | 2019-08-30 | 2019-12-24 | 中国科学院工程热物理研究所 | Thermodynamic simulation method for counter-flow heat regenerator |
Non-Patent Citations (2)
| Title |
|---|
| 吸收式热泵中吸收器传热传质与匹配特性研究;李静原;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20171215;第C038-46页 * |
| 多能微网阶段化多模式混合仿真关键技术综述;张树卿等;《计算机研究与发展》;20170415;第54卷(第4期);第683-694页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111597679A (en) | 2020-08-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101650095B (en) | Multistage absorption refrigerating/heat pump unit | |
| CN105042931B (en) | A kind of trans critical cycle and the combined heat-pump system of absorption heat pump coproduction | |
| WO2017020767A1 (en) | Multi-stage plate-type evaporation absorption cooling device and method | |
| Xu et al. | Theoretical analysis and optimization of a double-effect parallel-flow-type absorption chiller | |
| CN107091542A (en) | A kind of coupling circulation system and control method for solar energy thermal-power-generating | |
| CN109858163B (en) | Active disturbance rejection feedforward control method of lithium bromide absorption refrigerator | |
| Wang | A new style solar-driven diffusion absorption refrigerator and its operating characteristics | |
| CN207230989U (en) | Spray the dual temperature sweat cooling system of synergy | |
| CN102491440B (en) | Solar energy jet-refrigeration seawater desalting plant | |
| CN102226447B (en) | Medium-low temperature terrestrial heat power generating set system device | |
| CN111597679B (en) | Dynamic calculation method for external characteristic parameters of absorption heat pump for comprehensive energy network | |
| Mushtaq et al. | Exergy destruction rate minimization in the absorber of a double effect vapor absorption system | |
| CN109269143B (en) | Novel absorption heat pump and application method thereof | |
| CN108895708B (en) | Waste heat gradient recycling device and working method | |
| Singh et al. | A study on performance comparison of air cooled and water cooled condenser in vapour absorption and compression refrigeration systems | |
| Anand et al. | Exergy analysis of a LiBr–H 2 O vapor absorption refrigeration plant: a case study | |
| CN212299535U (en) | Absorption type refrigerating unit for preparing low-temperature cold water | |
| CN211233437U (en) | Single-effect heating and double-effect refrigerating steam type lithium bromide absorption heat pump unit | |
| CN211233438U (en) | Steam type lithium bromide absorption heat pump unit capable of achieving double-effect refrigeration | |
| CN211233436U (en) | Two-section lithium bromide absorption heat pump unit capable of achieving double-effect refrigeration | |
| Rameshkumar | Study On The Optimization Of Solar Assisted Vapour Absorption Refrigeration System | |
| CN108387024B (en) | Low-temperature waste heat recovery type lithium bromide absorption type hot and cold water system | |
| Adenane et al. | Effect of Operating Temperatures and Working Pairs on Performance of Solar Adsorption Cooling System. | |
| Qian et al. | Thermal analysis of a novel single-effect absorption refrigeration system using water/ionic liquid as working fluids | |
| CN104315747B (en) | A kind of utilize low-temperature waste heat for refrigeration equipment |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210622 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |